<div dir="ltr"><div>Hi</div><div><br></div><div>I am calling the python method osgeo.gdal.VectorTranslate to convert a GeoJSON (with CRS) to shapefile.</div><div>To keep consistency among EPSG versions, I have the pipeline of the conversion, using what in pyproj was "always_xy=True"</div><div><br></div><div>            gdal.VectorTranslate(<br>                dst_file,<br>                src_file,</div><div>                coorindateOperation=coordinateOperation,</div><div>                dstSRS="EPSG:4326",<br>            )</div><div><br></div><div>with something like this:<br></div><div><span>coordinateOperation="</span>+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"</div><div><div><br></div><div>the srcSRS is taken automatically from the GeoJSON.<br></div><div><br></div></div><div>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"</div><div><br></div><div></div><div><div>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).</div><div><br></div><div>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?</div><div><br></div><div>What about the command line gdal_translate?<br></div><br><div>Thank you.</div></div><div>Javier</div><br><div><div><div><div dir="ltr" class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div></div><div>.___ ._ ..._ .. . ._.  .___ .. __ . _. . __..  ... .... ._ .__</div></div></div></div></div></div></div>