Re: [gdal-dev] Discrepancy when tiling on MAC and Linux

2019-11-14 Thread Grégory Bataille
Hey Even, Actually a "simpler" solution for now for us is to remove the grid from the docker image we use to run dataset tiling. This way, the projection follows TOWGS84. It's suboptimal in term of absolute accuracy of the tiles, but it gives us consistency in our data which is more important to u

Re: [gdal-dev] Discrepancy when tiling on MAC and Linux

2019-11-14 Thread Grégory Bataille
Oh, I could gdalwarp the dataset myself to WGS84 and then tile it, so as to consistently "ignore" grids. That's an idea indeed. Thanks for the pointer, I'll take it into consideration --- Gregory Bataille On Thu, Nov 14, 2019 at 12:46 PM Even Rouault wrote: > On jeudi 14 novembre 2019 05:54:0

Re: [gdal-dev] Discrepancy when tiling on MAC and Linux

2019-11-14 Thread Even Rouault
On jeudi 14 novembre 2019 05:54:09 CET Grégory Bataille wrote: > wow, thanks. > > I'm not quite sure I'm following it all though. I'm reaching a bit my > knowledge limit here with those transforms. I'll need to spend a bit more > time on it to get more understanding and see what I can do, how I mi

Re: [gdal-dev] Discrepancy when tiling on MAC and Linux

2019-11-13 Thread Grégory Bataille
wow, thanks. I'm not quite sure I'm following it all though. I'm reaching a bit my knowledge limit here with those transforms. I'll need to spend a bit more time on it to get more understanding and see what I can do, how I might use your update. But thanks for the quick turnaround Cheers --- Gre

Re: [gdal-dev] Discrepancy when tiling on MAC and Linux

2019-11-13 Thread Even Rouault
> This is not a PROJ6 issue, but a GDAL one related to the logic starting at > https://github.com/OSGeo/gdal/blob/master/gdal/ogr/ogrct.cpp#L843 > > You could try to add a test that if OGRSpatialReference::GetTOWGS84() > returns OGRERR_NONE then the resolution to the WKT of the EPSG code is > skip

Re: [gdal-dev] Discrepancy when tiling on MAC and Linux

2019-11-13 Thread Even Rouault
On mercredi 13 novembre 2019 13:18:52 CET Grégory Bataille wrote: > Hey Even, > > So I can confirm that this is the issue > On Mac > OGRCT: Selecting transformation +proj=pipeline +step +proj=axisswap > +order=2,1 +step +inv +proj=tmerc +lat_0=0 +lon_0=12 +k=1 +x_0=450 > +y_0=0 +ellps=bessel +

Re: [gdal-dev] Discrepancy when tiling on MAC and Linux

2019-11-13 Thread Grégory Bataille
Hey Even, So I can confirm that this is the issue On Mac OGRCT: Selecting transformation +proj=pipeline +step +proj=axisswap +order=2,1 +step +inv +proj=tmerc +lat_0=0 +lon_0=12 +k=1 +x_0=450 +y_0=0 +ellps=bessel +step +proj=unitconvert +xy_in=rad +xy_out=deg +step +proj=axisswap +order=2,1 (I

Re: [gdal-dev] Discrepancy when tiling on MAC and Linux

2019-11-08 Thread Grégory Bataille
That’s a problem for us. With gdal 3 we started to force values for the towgs84 so that proj would not decide. This is wanted because our front end will place markers and it does coordinates transform. Without access to proj 6. Meaning it needs to have the transform Value. I’ll check what you ment

Re: [gdal-dev] Discrepancy when tiling on MAC and Linux

2019-11-08 Thread Even Rouault
On vendredi 8 novembre 2019 08:18:11 CET Grégory Bataille wrote: > Hey all, > > So I have a strange issue that I somewhat isolated, but I don't know how to > move further. > I have dataset in EPSG 31468-15949 I would check if that wouldn't be related to availability of the BETA2007.gsb grid. osg

[gdal-dev] Discrepancy when tiling on MAC and Linux

2019-11-07 Thread Grégory Bataille
Hey all, So I have a strange issue that I somewhat isolated, but I don't know how to move further. I have dataset in EPSG 31468-15949 My goal is to tile it, render it with openlayers, and add a marker on the GCP location. Trouble is that the GCP location was falling a bit away from the marker ima