<br><tt><font size=2>gdalwarp -srcnodata -3.4028234663852886e+38 \<br>
 &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;EstRN.img EstSh.img icu.tif<br>
</font></tt>
<br><font size=2 face="sans-serif">I thought I must launch gdalwrap once
for each file. Didn't know I could specify several input files</font>
<br>
<br><font size=2 face="sans-serif">I think I need to specify the extent
and the resolution.</font>
<br><font size=2 face="sans-serif">- Is the resolution is the pixel size
get by gdalinfo?</font>
<br>
<br><font size=2 face="sans-serif">- The Pixel Size is not the same for
all images. Per example I got </font>
<br><font size=2 face="sans-serif">&nbsp;(30.000000000000099,-30.000000000000103)</font>
<br><font size=2 face="sans-serif">and </font>
<br><font size=2 face="sans-serif">this for another image</font>
<br><font size=2 face="sans-serif">Pixel Size = (30.000000000000000,-30.000000000000000)</font>
<br>
<br><font size=2 face="sans-serif">What will append if I specify -tr 30
-30 as resolution in this case</font>
<br>
<br><font size=2 face="sans-serif">thanks</font>
<br><font size=2 face="sans-serif">Steve<br>
</font><font size=3 color=#0066ff><i>Steve Toutant, M. Sc.</i></font><font size=3><br>
Analyste en géomatique<br>
Secteur environnement<br>
Direction des risques biologiques, environnementaux et occupationnels<br>
Institut national de santé publique du Québec<br>
945, avenue Wolfe<br>
Québec, Qc G1V 5B3 </font>
<p><font size=3>Tél.: (418) 650-5115 #5281<br>
Fax.: (418) 654-3144</font><font size=3 color=blue><u><br>
</u></font><a href=mailto:steve.toutant@inspq.qc.ca><font size=3 color=blue><u>steve.toutant@inspq.qc.ca</u></font></a><font size=3 color=blue><u><br>
</u></font><a href=http://www.inspq.qc.ca/><font size=3 color=blue><u>http://www.inspq.qc.ca</u></font></a>
<p><font size=3>&nbsp; </font>
<p>
<br>
<br>
<br>
<table width=100%>
<tr valign=top>
<td width=40%><font size=1 face="sans-serif"><b>Frank Warmerdam &lt;warmerdam@pobox.com&gt;</b>
</font>
<p><font size=1 face="sans-serif">30/07/2009 02:19 PM</font>
<td width=59%>
<table width=100%>
<tr valign=top>
<td>
<div align=right><font size=1 face="sans-serif">A</font></div>
<td><font size=1 face="sans-serif">Steve.Toutant@inspq.qc.ca</font>
<tr valign=top>
<td>
<div align=right><font size=1 face="sans-serif">cc</font></div>
<td><font size=1 face="sans-serif">gdal-dev@lists.osgeo.org</font>
<tr valign=top>
<td>
<div align=right><font size=1 face="sans-serif">Objet</font></div>
<td><font size=1 face="sans-serif">Re: [gdal-dev] merge with gdalwrap</font></table>
<br>
<table>
<tr valign=top>
<td>
<td></table>
<br></table>
<br>
<br>
<br>
<br>
<br><tt><font size=2>Steve.Toutant@inspq.qc.ca wrote:<br>
&gt; <br>
&gt; Hi,<br>
&gt; I tried gdal_merge.py but I got a memory Error because I use the -n
<br>
&gt; nodata value option. Several threads sugest that I must use gdalwrap
<br>
&gt; instead.<br>
&gt; This behavior of gdal_merge.py is very important to me in this <br>
&gt; particular case<br>
&gt; &quot;In areas of overlap, the last image will be copied over earlier
ones&quot;.<br>
&gt; <br>
&gt; Does gdalwrap has the same behavior? (I ask because it is difficult
to <br>
&gt; test and validate with my data)?<br>
<br>
Steve,<br>
<br>
Yes, gdalwarp will overlay valid new data over existing data with the<br>
last gdalwarp overlaying previous data.<br>
<br>
&gt; Also, the doc says &quot;Mosaicing into an existing output file is
supported <br>
&gt; if the output file already exists.&quot;<br>
&gt; But how? What is the option to append data? When launching a command
for <br>
&gt; another image it says that I must delete the output file.<br>
<br>
If you pass an output file to gdalwarp and it already exists then it<br>
will be mosaiced into. &nbsp;If the output file does not already exist<br>
then one is created. &nbsp;This is different than gdal_translate which<br>
always creates a new file.<br>
<br>
&gt; I would appreciate if someone could validate if this<br>
&gt; gdalwarp -srcnodata -3.4028234663852886e+38 -dstnodata <br>
&gt; -3.4028234663852886e+38 -of &quot;GTiff&quot; EstRN.img icu.tif<br>
&gt; <br>
&gt; is the equivalent of this<br>
&gt; gdal_merge.py -n -3.4028234663852886e+38 -o icu.tif -of &quot;GTiff&quot;
<br>
&gt; EstRN.img EstSh.img<br>
<br>
Well, it doesn't mention EstSH.img - so no it does not seem<br>
to be the same thing. &nbsp;If you did:<br>
<br>
gdalwarp -srcnodata -3.4028234663852886e+38 \<br>
 &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;EstRN.img EstSh.img icu.tif<br>
<br>
then it should be pretty similar. &nbsp;There are subtle differences<br>
of course.<br>
<br>
Best regards,<br>
-- <br>
---------------------------------------+--------------------------------------<br>
I set the clouds in motion - turn up &nbsp; | Frank Warmerdam, warmerdam@pobox.com<br>
light and sound - activate the windows | http://pobox.com/~warmerdam<br>
and watch the world go round - Rush &nbsp; &nbsp;| Geospatial Programmer
for Rent<br>
<br>
</font></tt>
<br>
<br>
<br>