<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
<style type="text/css" style="display:none;"><!-- P {margin-top:0;margin-bottom:0;} --></style>
</head>
<body dir="ltr">
<div id="divtagdefaultwrapper" style="font-size:12pt;color:#000000;font-family:Calibri,Helvetica,sans-serif;" dir="ltr">
<p style="margin-top:0;margin-bottom:0">Thank you for your help,</p>
<p style="margin-top:0;margin-bottom:0">I know that is a possibility but this is to process big amounts of data and I would like to avoid the merge processing time if possible by rasterizing over the same tif...</p>
<p style="margin-top:0;margin-bottom:0"><br>
</p>
<p style="margin-top:0;margin-bottom:0">I would like to know if it is possible...</p>
<p style="margin-top:0;margin-bottom:0"><br>
</p>
<p style="margin-top:0;margin-bottom:0">cheers,</p>
<p style="margin-top:0;margin-bottom:0">Mario</p>
</div>
<hr style="display:inline-block;width:98%" tabindex="-1">
<div id="divRplyFwdMsg" dir="ltr"><font face="Calibri, sans-serif" style="font-size:11pt" color="#000000"><b>De:</b> Even Rouault <even.rouault@spatialys.com><br>
<b>Enviado:</b> jueves, 7 de junio de 2018 13:05:03<br>
<b>Para:</b> gdal-dev@lists.osgeo.org<br>
<b>Cc:</b> Nunez, Mario<br>
<b>Asunto:</b> Re: [gdal-dev] gdal_rasterize failing when trying multiband</font>
<div> </div>
</div>
<div class="BodyFragment"><font size="2"><span style="font-size:11pt;">
<div class="PlainText">Mario,<br>
<br>
I suggest you use gdal_rasterize to produce as many GeoTIFF files as you have <br>
'burning sessions'<br>
And then use gdal_merge.py -separate to make a single GeoTIFF from all the <br>
above ones.<br>
<br>
Even<br>
<br>
> Hello,<br>
> <br>
> I was doing some point rasterization of vector data from postgis database by<br>
> using gdal_rasterize from a bash script, the command looks like this:<br>
> <br>
> <br>
> gdal_rasterize -a value -b $COUNT -tr $RESOLUTION -te $BBOX -a_nodata<br>
> $NODATA -a_srs $EPSG PG:'host=localhost dbname=databio user=postgres<br>
> password=postgres' -sql "SELECT ST_Centroid(g.geom), AVG($COLUMN) as value<br>
> FROM $GRID g, $OBSERVACION o WHERE utc_date='$DATE' AND<br>
> vessel_cfr='$VESSEL' AND ST_Intersects(o.geom,g.geom) IS TRUE GROUP BY<br>
> g.gid" ./bands/$DATE2"_"$COLUMN".tif"<br>
> <br>
> <br>
> The problem is that I want to generate a geotif file with multiple bands, so<br>
> if I remove the -b $COUNT parameter and value it works and I get a .tif<br>
> file with the -a value printed on it for band1, but the idea is to reuse<br>
> that command in a for loop to add further bands and I am gettin this error:<br>
> <br>
> <br>
> -b option cannot be used when creating a GDAL dataset.<br>
> <br>
> I tried to run the comand once to create the tif file without the -b option<br>
> and then run it additional times with the -b option when the tif is already<br>
> created but it makes no difference.<br>
> <br>
> <br>
> I am stucked with this and I am starting to think that -b option can only be<br>
> used with -burn parmeter and as far as I know -burn is only for constant<br>
> values...<br>
> <br>
> <br>
> thank you in advance for any help.<br>
> <br>
> cheers,<br>
> <br>
> Mario<br>
> <br>
> <br>
> This e-mail and the documents attached are confidential and intended solely<br>
> for the addressee; it may also be privileged. If you receive this e-mail in<br>
> error, please notify the sender immediately and destroy it. As its<br>
> integrity cannot be secured on the Internet, the Atos group liability<br>
> cannot be triggered for the message content. Although the sender endeavors<br>
> to maintain a computer virus-free network, the sender does not warrant that<br>
> this transmission is virus-free and will not be liable for any damages<br>
> resulting from any virus transmitted.<br>
> <br>
> Este mensaje y los ficheros adjuntos pueden contener información<br>
> confidencial destinada solamente a la(s) persona(s) mencionadas<br>
> anteriormente y pueden estar protegidos por secreto profesional. Si usted<br>
> recibe este correo electrónico por error, gracias por informar<br>
> inmediatamente al remitente y destruir el mensaje. Al no estar asegurada la<br>
> integridad de este mensaje sobre la red, Atos no se hace responsable por su<br>
> contenido. Su contenido no constituye ningún compromiso para el grupo Atos,<br>
> salvo ratificación escrita por ambas partes. Aunque se esfuerza al máximo<br>
> por mantener su red libre de virus, el emisor no puede garantizar nada al<br>
> respecto y no será responsable de cualesquiera daños que puedan resultar de<br>
> una transmisión de virus.<br>
<br>
<br>
-- <br>
Spatialys - Geospatial professional services<br>
<a href="http://www.spatialys.com">http://www.spatialys.com</a><br>
</div>
</span></font></div>
This e-mail and the documents attached are confidential and intended solely for the addressee; it may also be privileged. If you receive this e-mail in error, please notify the sender immediately and destroy it.
<br>
As its integrity cannot be secured on the Internet, the Atos group liability cannot be triggered for the message content. Although the sender endeavors to maintain a computer virus-free network, the sender does not warrant that this transmission is virus-free
 and will not be liable for any damages resulting from any virus transmitted. <br>
<br>
Este mensaje y los ficheros adjuntos pueden contener información confidencial destinada solamente a la(s) persona(s) mencionadas anteriormente y pueden estar protegidos por secreto profesional.
<br>
Si usted recibe este correo electrónico por error, gracias por informar inmediatamente al remitente y destruir el mensaje.
<br>
Al no estar asegurada la integridad de este mensaje sobre la red, Atos no se hace responsable por su contenido. Su contenido no constituye ningún compromiso para el grupo Atos, salvo ratificación escrita por ambas partes.
<br>
Aunque se esfuerza al máximo por mantener su red libre de virus, el emisor no puede garantizar nada al respecto y no será responsable de cualesquiera daños que puedan resultar de una transmisión de virus.
</body>
</html>