[gdal-dev] Trouble reprojecting feature that crosses antimeridian

[Even Rouault]

Was a bug. Fix queued in https://github.com/OSGeo/gdal/pull/5543

Le 29/03/2022 à 22:30, Robin Wilson a écrit :
> Hi,
>
> I’ve been running into problems creating STAC items for images that 
> cross the antimeridian. The geometries for the images are coming out 
> as spanning the width of the whole world (from -180 to +180 degrees). 
> With the help of the author of the tool I was using to create the STAC 
> items (rio-stac), I’ve narrowed it down to a simple test case of 
> reprojecting geometries (which is what the tool is doing internally).
>
> If you put the following GeoJSON for an antimeridian-crossing polygon 
> in WGS 84 / UTM zone 1N in a file called test_32601.geojson:
>
> ```
> {"type":"Polygon","coordinates":[[[377120,7577600],[418080,7577600],[418080,7618560],[377120,7618560],[377120,7577600]]]}
> ```
>
> and then run the following:
>
> ```
> ogr2ogr -f GeoJSON -s_srs epsg:32601 -t_srs epsg:4326 out_4326.geojson 
> test_32601.geojson
> ```
>
> The contents of the output file will be:
>
> ```
> { "type": "Polygon", "coordinates": [ [ [ -179.018099960087909, 
> 68.666857319910392 ], [ -178.985579869598382, 68.299719485120988 ], [ 
> -179.976982644969041, 68.284942315413858 ], [ 179.974309855567753, 
> 68.651801088224246 ], [ 180.0, 68.652259945574556 ], [ 180.0, 
> 68.652259945541758 ], [ 180.0, 68.652259944541754 ], [ -180.0, 
> 68.652259944541754 ], [ -180.0, 68.652259945541758 ], [ -180.0, 
> 68.652259945610126 ], [ -179.018099960087909, 68.666857319910392 ] ] ] }
> ```
>
> In that output geometry there are longitudes of 180 and -180, so when 
> plotted on a map it spans the whole width of the world.
>
> The author of the rio-stac library did a PR to try and fix this 
> problem (https://github.com/developmentseed/rio-stac/pull/30), but it 
> doesn’t seem to have fixed it for this example, as the underlying GDAL 
> code used by his library (through rasterio, I believe) is giving this 
> result.
>
> I find this result unexpected, but I have no knowledge of how GDAL 
> deals with the antimeridian. Is this to be expected, or is this 
> potentially a bug? If it isn’t a bug, is there any way around this?
>
> Best regards,
>
> Robin
>
> Dr Robin Wilson
> www.rtwilson.com <http://www.rtwilson.com>
>
>
>
> _______________________________________________
> gdal-dev mailing list
> gdal-dev at lists.osgeo.org
> https://lists.osgeo.org/mailman/listinfo/gdal-dev

-- 
http://www.spatialys.com
My software is free, but my time generally not.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/gdal-dev/attachments/20220329/58df4442/attachment-0001.html>