<html xmlns:v="urn:schemas-microsoft-com:vml" xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40">
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
<meta name="Generator" content="Microsoft Word 15 (filtered medium)">
<!--[if !mso]><style>v\:* {behavior:url(#default#VML);}
o\:* {behavior:url(#default#VML);}
w\:* {behavior:url(#default#VML);}
.shape {behavior:url(#default#VML);}
</style><![endif]--><style><!--
/* Font Definitions */
@font-face
        {font-family:"Cambria Math";
        panose-1:2 4 5 3 5 4 6 3 2 4;}
@font-face
        {font-family:Calibri;
        panose-1:2 15 5 2 2 2 4 3 2 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
        {margin:0cm;
        margin-bottom:.0001pt;
        font-size:12.0pt;
        font-family:"Times New Roman",serif;}
a:link, span.MsoHyperlink
        {mso-style-priority:99;
        color:blue;
        text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
        {mso-style-priority:99;
        color:purple;
        text-decoration:underline;}
p
        {mso-style-priority:99;
        margin:0cm;
        margin-bottom:.0001pt;
        font-size:12.0pt;
        font-family:"Times New Roman",serif;}
pre
        {mso-style-priority:99;
        mso-style-link:"HTML con formato previo Car";
        margin:0cm;
        margin-bottom:.0001pt;
        font-size:10.0pt;
        font-family:"Courier New";}
p.msonormal0, li.msonormal0, div.msonormal0
        {mso-style-name:msonormal;
        margin:0cm;
        margin-bottom:.0001pt;
        font-size:12.0pt;
        font-family:"Times New Roman",serif;}
span.EstiloCorreo19
        {mso-style-type:personal-reply;
        font-family:"Calibri",sans-serif;
        color:windowtext;}
span.HTMLconformatoprevioCar
        {mso-style-name:"HTML con formato previo Car";
        mso-style-priority:99;
        mso-style-link:"HTML con formato previo";
        font-family:"Courier New";}
.MsoChpDefault
        {mso-style-type:export-only;
        font-size:10.0pt;}
@page WordSection1
        {size:612.0pt 792.0pt;
        margin:72.0pt 72.0pt 72.0pt 72.0pt;}
div.WordSection1
        {page:WordSection1;}
--></style><!--[if gte mso 9]><xml>
<o:shapedefaults v:ext="edit" spidmax="1026" />
</xml><![endif]--><!--[if gte mso 9]><xml>
<o:shapelayout v:ext="edit">
<o:idmap v:ext="edit" data="1" />
</o:shapelayout></xml><![endif]-->
</head>
<body lang="ES" link="blue" vlink="purple">
<div class="WordSection1">
<p class="MsoNormal"><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:black">Anybody knows if this is possible?<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";color:black"><o:p> </o:p></span></p>
<p class="MsoNormal"><span style="font-size:10.0pt;font-family:"Courier New";color:black">gdal_rasterize -a value -b $COUNT -tr $RESOLUTION -te $BBOX -a_nodata $NODATA -a_srs $EPSG PG:'host=localhost dbname=databio user=postgres password=postgres' -sql "SELECT
 ST_Centroid(g.geom), AVG($COLUMN) as value FROM $GRID g, $OBSERVACION o WHERE utc_date='$DATE' AND vessel_cfr='$VESSEL' AND ST_Intersects(o.geom,g.geom) IS TRUE GROUP BY g.gid" ./bands/$DATE2"_"$COLUMN".tif"<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;mso-fareast-language:EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri",sans-serif;mso-fareast-language:EN-US">I mean using -b parameter to directly create a geotiff with different bands from a postgis database query.<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri",sans-serif;mso-fareast-language:EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri",sans-serif;mso-fareast-language:EN-US">I would like to avoid doing single rasterizations and then a merge.<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US" style="mso-fareast-language:EN-US">Thank you in advance.<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US" style="mso-fareast-language:EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US" style="mso-fareast-language:EN-US"><o:p> </o:p></span></p>
<div>
<div style="border:none;border-top:solid #E1E1E1 1.0pt;padding:3.0pt 0cm 0cm 0cm">
<p class="MsoNormal"><b><span style="font-size:11.0pt;font-family:"Calibri",sans-serif">From:</span></b><span style="font-size:11.0pt;font-family:"Calibri",sans-serif"> Nunez, Mario
<br>
<b>Sent:</b> Thursday, June 07, 2018 1:32 PM<br>
<b>To:</b> Even Rouault <even.rouault@spatialys.com>; gdal-dev@lists.osgeo.org<br>
<b>Subject:</b> Re: [gdal-dev] gdal_rasterize failing when trying multiband<o:p></o:p></span></p>
</div>
</div>
<p class="MsoNormal"><o:p> </o:p></p>
<div id="divtagdefaultwrapper">
<p><span style="font-family:"Calibri",sans-serif;color:black">Thank you for your help,<o:p></o:p></span></p>
<p><span style="font-family:"Calibri",sans-serif;color:black">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...<o:p></o:p></span></p>
<p><span style="font-family:"Calibri",sans-serif;color:black"><o:p> </o:p></span></p>
<p><span style="font-family:"Calibri",sans-serif;color:black">I would like to know if it is possible...<o:p></o:p></span></p>
<p><span style="font-family:"Calibri",sans-serif;color:black"><o:p> </o:p></span></p>
<p><span style="font-family:"Calibri",sans-serif;color:black">cheers,<o:p></o:p></span></p>
<p><span style="font-family:"Calibri",sans-serif;color:black">Mario<o:p></o:p></span></p>
</div>
<div class="MsoNormal" align="center" style="text-align:center">
<hr size="2" width="98%" align="center">
</div>
<div id="divRplyFwdMsg">
<p class="MsoNormal"><b><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:black">De:</span></b><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:black"> Even Rouault <<a href="mailto:even.rouault@spatialys.com">even.rouault@spatialys.com</a>><br>
<b>Enviado:</b> jueves, 7 de junio de 2018 13:05:03<br>
<b>Para:</b> <a href="mailto:gdal-dev@lists.osgeo.org">gdal-dev@lists.osgeo.org</a><br>
<b>Cc:</b> Nunez, Mario<br>
<b>Asunto:</b> Re: [gdal-dev] gdal_rasterize failing when trying multiband</span>
<o:p></o:p></p>
<div>
<p class="MsoNormal"> <o:p></o:p></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:11.0pt">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><o:p></o:p></span></p>
</div>
</div>
</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>