Using GeoDataFrames.jl I’m trying to reproject to Google Pseudo-Mercator (EPSG 3857) projection, but get wrong coordinates. In the original GeoJSON, coords are long/lat (see screenshot below): geometry 0 POINT (69.19970 41.36596) 1 POINT (69.62239 42.41753) 2 POINT (74.57231…

[image] culebron: Tried this, but the result is still the same :compliant is the default, can you pass order=:trad?

Ok. So GFT just provides wrappers for CRS so it can be shared between packages without re-parsing. It doesn’t really do anything with it. GeoDataFrames.jl uses gdal via ArchGDAL to both read and reproject data. It seems to me your data is being read as x/y rather than y/x at some point in the chai…

Indeed this looks like an axis order issue. The reproject function that you are using has a keyword argument order=:trad that you can use. The issue is that EPSG defines the argument order to be lat/lon, which is not the case here. And the axis order is taken based on the authority (in this case EPS…

It did the trick! Thanks! Hmmm… happy it worked, but sad I have to do such tricks. [изображение]

I’m still perplexed how cryptic the docs are. No info on return types. Also this example crashed with MethodError, couldn’t find how to call it on a vector. cities2 = GDF.read("/tmp/cities.geojson") sc = AG.importEPSG(4326) tc = AG.importEPSG(3857) AG.createcoordtrans(sc, tc) do transform AG.…

It’s best to make github issues for any problems you find. The docs can always be improved, but there are only a few people using small amounts of our spare time working on these packages. We all have other jobs. If you feel strongly about this you may like to get involved and for a start improve …

I fully agree that the order should always be lon,lat and not lat,lon but this is a PROJ thing that can be avoided in GDAL if OAMS_TRADITIONAL_GIS_ORDER is used in OSRSetAxisMappingStrategy() We have this in GMT (C lib) and AFAIK that is also the default in GMT.jl when it uses PROJ directly. #if G…

ArchGDAL transforming CRS incorrectly (probably wrong long/lat order)

New to Julia

Fliks November 11, 2022, 8:33am 3

Have you tried the order keyword argument to ArchGDAL.importEPSG? As far as I understand the docstring setting the order to trad could swap the axis to lon/lat.

Topic		Replies	Views
ArchGDAL.jl (or GDAL.jl not sure) - Undesired behaviour when reprojecting coordinates to EPSG 4326 Geo question , bug , geodesy	3	818	April 16, 2021
[ANN] CoordRefSystems.jl Package Announcements package , announcement , geo , geodesy	49	3722	December 9, 2025
How to convert a point from a given 2D crs to another one Geo coordinates	6	282	May 22, 2024
Status of geodesy packages in Julia (DGGS experiments) Geo geo , geodesy	12	371	October 25, 2025
How to transform CRS of entire shapefile (not a single point)? Geo geo , gdal , geodesy	6	998	June 18, 2024

ArchGDAL transforming CRS incorrectly (probably wrong long/lat order)

Related topics