[gdal-dev] Re-post, user-defined

[Even Rouault]

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 at lists.osgeo.org
> https://lists.osgeo.org/mailman/listinfo/gdal-dev

-- 
http://www.spatialys.com
My software is free, but my time generally not.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/gdal-dev/attachments/20240429/be5ab434/attachment.htm>