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 file generate by GDAL
- links to all the above files
- the GDAL version you use
Totally blind guess: try to add "+towgs84=0,0,0 " in your -t_srs PROJ
string, to ask PROJ to consider that your target CRS datum has the same
center and axis orientation as WGS84 (and thus go to geographic <-->
geocentric conversion when changing datums). Otherwise without datum
shifts, PROJ will consider that longitude,latitude coordinates of the
geographic CRS underlying the source and target CRS is the same.
Even
Le 29/04/2024 à 14:57, Liyuneh Gebre via gdal-dev a écrit :
Dear all,
I am requesting again if you kind have some responds or direction,
thank you
I have difficulties in accurately reprojecting a raster file to a
user-defined projection.
I used the following commands in gdal while the output is not similar
to arcpy results. kindly offer your guidance to achieve the same
result. The arcpy uses custom transformation and then reprojection.
gdalwarp -t_srs '+proj=laea +a=6378137.0 +f=0.0 +pm=0 +x_0=0.0
+y_0=0.0 +lon_0=20.0 +lat_0=5.0 +units=m +axis=enu +no_defs'
-dstalpha -wo SOURCE_EXTRA=1000 -r nearest rain_in.bil rain_projected.bil
/*Liyuneh Gebre*//, ////Associate Researcher///
/Ethiopian Institute of Agricultural Research,EIAR
/
/Cell Phone:+251 911858155
/
/linkedIn: https://www.linkedin.com/in/liyuneh-gebre-b842b440/
/
//skype:liyenew_1//
/Addis Ababa, Ethioipa/
_______________________________________________
gdal-dev mailing list
gdal-dev@lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/gdal-dev
--
http://www.spatialys.com
My software is free, but my time generally not.
_______________________________________________
gdal-dev mailing list
gdal-dev@lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/gdal-dev
