<div dir="ltr"><div>Hi Even,</div><div>thank you for your response,</div><div>gdal version is: GDAL 3.6.4, released 2023/04/17</div><div>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.</div><div><a href="https://u.pcloud.link/publink/show?code=kZj5wh0ZtFWXtI9UGjQiD6KXxn5A9hiWDpEX">https://u.pcloud.link/publink/show?code=kZj5wh0ZtFWXtI9UGjQiD6KXxn5A9hiWDpEX</a></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Mon, Apr 29, 2024 at 4:23 PM Even Rouault <<a href="mailto:even.rouault@spatialys.com">even.rouault@spatialys.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><u></u>
<div>
<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>Le 29/04/2024 à 14:57, Liyuneh Gebre
via gdal-dev a écrit :<br>
</div>
<blockquote type="cite">
<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;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;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">
<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">https://www.linkedin.com/in/liyuneh-gebre-b842b440/</a><br>
</i></div>
<div><i><i><a>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></fieldset>
<pre>_______________________________________________
gdal-dev mailing list
<a href="mailto:gdal-dev@lists.osgeo.org" target="_blank">gdal-dev@lists.osgeo.org</a>
<a href="https://lists.osgeo.org/mailman/listinfo/gdal-dev" target="_blank">https://lists.osgeo.org/mailman/listinfo/gdal-dev</a>
</pre>
</blockquote>
<pre cols="72">--
<a href="http://www.spatialys.com" target="_blank">http://www.spatialys.com</a>
My software is free, but my time generally not.</pre>
</div>
</blockquote></div><br clear="all"><br><span class="gmail_signature_prefix">-- </span><br><div dir="ltr" class="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">https://www.linkedin.com/in/liyuneh-gebre-b842b440/</a><br></i></div><div><i><i>skype:liyenew_1</i></i></div><div><i></i></div><div><i>Addis Ababa, Ethioipa</i></div></div></div></div></div></div></div></div></div></div></div>