Greetings<div><br></div><div>A couple of days ago I have exchanged a few messages with Paul Kelly and Micha Silver regarding this issue. But now I have a new example and I want to know if I'm doing anything wrong and better understand what is GRASS behaviour.</div>
<div><br></div><div>I'm trying to have an aumatic procedure to reproject images using gdalwarp. First it checks image CRS (g.proj input=IMAGE -jf) and compares with Location's CRS (g.proj -jf).</div><div><br></div>
<div>As target_CRS I'm using the output from g.proj -j -f . </div><div>In this case, g.proj -jf give me the Location«s CRS:</div><div>+proj=utm +south +no_defs +zone=24 +a=6378137 +rf=298.257223563 +towgs84=0.000,0.000,0.000 +to_meter=1</div>
<div>(Sidenote. This Location was created based on an EPSG)</div><div><br></div><div>I use this expression to warp the image:</div><div><div>gdalwarp -t_srs "+proj=utm +south +no_defs +zone=24 +a=6378137 +rf=298.257223563 +towgs84=0.000,0.000,0.000 +to_meter=1" input.tif proj4.tif</div>
</div><div><br></div><div>But if I check proj4.tif CRS I get this:</div><div><div>+proj=utm +south +no_defs +zone=24 +a=6378137 +rf=298.257223563 +to_meter=1</div><div>Trying to open with OGR...</div><div>Trying to open with GDAL...</div>
<div>...succeeded.</div><div>Datum <unknown> not recognised by GRASS and no parameters found</div></div><div><br></div><div>Why towgs84 is not kept in the image?</div><div><br></div><div>THanks</div><div><br></div><div>
Kim</div><div><br></div>