Dear Even,
Apologies for taking your time. I couldn’t locate the mentioned error in
arcpy. It runs smoothly . I’ve shared the arcpy projection python script
and associated data with the link
https://u.pcloud.link/publink/show?code=kZj5wh0ZtFWXtI9UGjQiD6KXxn5A9hiWDpEX.
If you find a way using gdal
Well the arcpy made transformation before re-projection. Appended the
sample script in the shared folder as txt. Anyway how can I achieve this
ArcPy output re-projection from gdal with the given gdalinfo on the input
file and by editing the PROJ.4 string (target projection )?
PARAMETER["Latitude o
At first sight, I'd say that the ArcPy output is incorrect. The SRS
definition of its output doesn't match what you want since it has
PARAMETER["Latitude of natural origin",20,
ANGLEUNIT["degree",0.0174532925199433],
ID["EPSG",8801]],
PARAMETER["Longitude of natural origin",6378137,
ANGLEUNIT["d
Hi Even,
thank you for your response,
gdal version is: GDAL 3.6.4, released 2023/04/17
I got the defined target projection included in the link. yes, it's
geocentric and for further exploration the custom project transformation
file .gtf extension is also included.
https://u.pcloud.link/publink/sho
Hi,
it is difficult to help you with just the elements you've provided.
What could be helpful is:
- the output of gdalinfo on the input file
- the output of gdalinfo on the output file generate by ArcPy
- how you define the target projection in ArcPy
- the output of gdalinfo on the output fi