<!DOCTYPE html>
<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=UTF-8">
</head>
<body>
<p>Hi,</p>
<p>it is difficult to help you with just the elements you've
provided.</p>
<p>What could be helpful is:</p>
<p>- the output of gdalinfo on the input file</p>
<p>- the output of gdalinfo on the output file generate by ArcPy</p>
<p>- how you define the target projection in ArcPy<br>
</p>
<p>- the output of gdalinfo on the output file generate by GDAL</p>
<p>- links to all the above files<br>
</p>
<p>- the GDAL version you use</p>
<p><br>
</p>
<p>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.<br>
</p>
<p>Even<br>
</p>
<div class="moz-cite-prefix">Le 29/04/2024 à 14:57, Liyuneh Gebre
via gdal-dev a écrit :<br>
</div>
<blockquote type="cite"
cite="mid:CAE4zf46tV7mfiEAYKrMjhDQVdfRPnSpzdcje3a47FE0HvMgY3A@mail.gmail.com">
<meta http-equiv="content-type" content="text/html; charset=UTF-8">
<div>
<div>
<div
style="color:rgb(49,49,49);font-family:'-apple-system','helvetica neue';font-size:1rem;font-style:normal;font-weight:400;letter-spacing:normal;text-indent:0px;text-transform:none;white-space:normal;word-spacing:1px;text-decoration:none">Dear
all,</div>
<div
style="color:rgb(49,49,49);font-family:'-apple-system','helvetica neue';font-size:1rem;font-style:normal;font-weight:400;letter-spacing:normal;text-indent:0px;text-transform:none;white-space:normal;word-spacing:1px;text-decoration:none"
dir="auto">I am requesting again if you kind have some
responds or direction, thank you</div>
<div
style="color:rgb(49,49,49);font-family:'-apple-system','helvetica neue';font-size:1rem;font-style:normal;font-weight:400;letter-spacing:normal;text-indent:0px;text-transform:none;white-space:normal;word-spacing:1px;text-decoration:none"
dir="auto">I have difficulties in accurately reprojecting a
raster file to a user-defined projection.</div>
<div
style="color:rgb(49,49,49);font-family:'-apple-system','helvetica neue';font-size:1rem;font-style:normal;font-weight:400;letter-spacing:normal;text-indent:0px;text-transform:none;white-space:normal;word-spacing:1px;text-decoration:none">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. <br>
</div>
<span
style="color:rgb(49,49,49);font-family:'-apple-system','helvetica neue';font-size:16px;font-style:normal;font-weight:400;letter-spacing:normal;text-indent:0px;text-transform:none;white-space:normal;word-spacing:1px;text-decoration:none;display:inline!important;float:none">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 </span><span
style="color:rgb(49,49,49);font-family:'-apple-system','helvetica neue';font-size:1rem;font-style:normal;font-weight:400;letter-spacing:normal;text-indent:0px;text-transform:none;white-space:normal;word-spacing:1px;text-decoration:none">-wo
SOURCE_EXTRA=1000 </span><span
style="color:rgb(49,49,49);font-family:'-apple-system','helvetica neue';font-size:16px;font-style:normal;font-weight:400;letter-spacing:normal;text-indent:0px;text-transform:none;white-space:normal;word-spacing:1px;text-decoration:none;display:inline!important;float:none">-r
nearest rain_in.bil rain_projected.bil</span></div>
<br clear="all">
<br clear="all">
<div>
<div dir="ltr" class="gmail_signature"
data-smartmail="gmail_signature">
<div dir="ltr">
<div>
<div dir="ltr">
<div>
<div dir="ltr">
<div dir="ltr">
<div dir="ltr">
<div dir="ltr">
<div dir="ltr"><i><b>Liyuneh Gebre</b></i><i>,
</i><i><i><span><i>Associate Researcher</i></span></i></i></div>
<div><i>Ethiopian Institute of Agricultural
Research,EIAR<br>
</i></div>
<div dir="ltr">
<div><i>Cell Phone:+251 911858155<br>
</i></div>
<div><i>linkedIn: <a
href="https://www.linkedin.com/in/liyuneh-gebre-b842b440/"
target="_blank"
moz-do-not-send="true"
class="moz-txt-link-freetext">https://www.linkedin.com/in/liyuneh-gebre-b842b440/</a><br>
</i></div>
<div><i><i><a class="moz-txt-link-freetext" href="skype:liyenew_1">skype:liyenew_1</a></i></i></div>
<div><i>Addis Ababa, Ethioipa</i></div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
<br>
<fieldset class="moz-mime-attachment-header"></fieldset>
<pre class="moz-quote-pre" wrap="">_______________________________________________
gdal-dev mailing list
<a class="moz-txt-link-abbreviated" href="mailto:gdal-dev@lists.osgeo.org">gdal-dev@lists.osgeo.org</a>
<a class="moz-txt-link-freetext" href="https://lists.osgeo.org/mailman/listinfo/gdal-dev">https://lists.osgeo.org/mailman/listinfo/gdal-dev</a>
</pre>
</blockquote>
<pre class="moz-signature" cols="72">--
<a class="moz-txt-link-freetext" href="http://www.spatialys.com">http://www.spatialys.com</a>
My software is free, but my time generally not.</pre>
</body>
</html>