Hi I am calling the python method osgeo.gdal.VectorTranslate to convert a GeoJSON (with CRS) to shapefile. To keep consistency among EPSG versions, I have the pipeline of the conversion, using what in pyproj was "always_xy=True"
gdal.VectorTranslate( dst_file, src_file, coorindateOperation=coordinateOperation, dstSRS="EPSG:4326", ) with something like this: coordinateOperation="+proj=pipeline +step +inv +proj=tmerc +lat_0=44 +lon_0=144.25 +k=0.9999 +x_0=0 +y_0=0 +ellps=GRS80 +step +proj=unitconvert +xy_in=rad +xy_out=deg" the srcSRS is taken automatically from the GeoJSON. The first thing I noticed is that EPSG:4326 is left handed, and I have to swap the axes. Ok, it is easy to add " +step +proj=axisswap +order=2,1" In the case of EPSG:2455 (left handed CRS in Japan) I should do the swap also for the input values, that is a bit uglier. (I have to do some parsing of the pipeline that was generated time ago, and know that EPSG:2455 is left handed). Is there a way to tell gdal.VectorTranslate that transformation is "always_xy", and I do not have to do the swaps manually? Or maybe there is a better way to do it? What about the command line gdal_translate? Thank you. Javier .___ ._ ..._ .. . ._. .___ .. __ . _. . __.. ... .... ._ .__
_______________________________________________ gdal-dev mailing list gdal-dev@lists.osgeo.org https://lists.osgeo.org/mailman/listinfo/gdal-dev