<div dir="ltr"><div>Thank for the fast answers ! </div><div><br></div>@EvenRouault  : thanks for the snippet, I wasn't aware of this mechanism to transmit option trough python dictionary.<br>It should most likely be sufficient for my needs, even if I don't understand it fully yet.<div>I had read the tutorial but failed to translate it from C++ to python.</div><div>I'm still _very_ reticent to create XML for VRT in python, I don't see how it could be safe and age well. <br><br>@Anton :</div><div>thanks for the lib, I'll look into <a href="https://github.com/nansencenter/nansat/blob/master/nansat/vrt.py">vrt.py</a> to try to understand how to add an existing tiff file as a band to vrt.</div><div>I don't think I need "<span style="font-size:12.8000001907349px">SourceTransferType</span>", because the type of the vrt band will be the same as the type of the existing geotiff file.</div><div><br></div><div>Cheers,</div><div>Rémi-C<br><div><br></div></div></div><div class="gmail_extra"><br><div class="gmail_quote">2015-02-16 15:11 GMT+01:00 Anton Korosov <span dir="ltr"><<a href="mailto:anton.korosov@nersc.no" target="_blank">anton.korosov@nersc.no</a>></span>:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Dear Rémi,<br>
<br>
you may want to look at our Python package NANSAT (<a href="https://github.com/nansencenter/nansat" target="_blank">https://github.com/<u></u>nansencenter/nansat</a>) which we build around the GDAL and especially VRT. We even have a separate class VRT() for creating, changing and using VRT files.<br>
More specific regarding your question, you can try to use the SourceTransferType parameter of the VRT. Search for this keyword at<br>
<a href="http://gdal.org/gdal_vrttut.html" target="_blank">http://gdal.org/gdal_vrttut.<u></u>html</a><br>
<br>
Regards!<br>
Anton<div><div class="h5"><br>
<br>
On 16/02/15 14:13, Rémi Cura wrote:<br>
</div></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div class="h5">
Hey dear list,<br>
this is my first post here,<br>
so sorry if I'm asking something obvious.<br>
<br>
(I noticed the 2 previous mail thread regarding python and VRT).<br>
<br>
I wrote python functions to convert a numpy structured array to several<br>
gdal rasters.<br>
<br>
A numpy structured array is a 2D pixel array, but each pixel has in fact<br>
several values of possibly different type (like band).<br>
For instance a pixel dtype can be [float32,uint16,Float64].<br>
<br>
I want to convert this array to raster(s).<br>
First I tried to create one raster with several bands,<br>
but GeoTiff cannot handle bands with different pixel type (or I didn't<br>
find how).<br>
<br>
I concluded that I must create several rasters, and put them all as<br>
bands in a VRT.<br>
<br>
So far I have been able to create several geoTiff raster of various type<br>
corresponding to the numpy structured array.<br>
<br>
Now, I don't know how to create from python the VRT that unites all this<br>
rasters.<br>
<br>
When I try to write in a VRT raster band I get errors.<br>
<br>
I am well aware that I could create the XML directly from python, but it<br>
feels like a pretty bad architecture idea.<br>
<br>
Anybody knows how to do this /and/or has comments about this task?<br>
<br>
Many Thanks,<br>
Rémi-C<br>
<br>
<br>
<br>
<br></div></div>
______________________________<u></u>_________________<br>
gdal-dev mailing list<br>
<a href="mailto:gdal-dev@lists.osgeo.org" target="_blank">gdal-dev@lists.osgeo.org</a><br>
<a href="http://lists.osgeo.org/mailman/listinfo/gdal-dev" target="_blank">http://lists.osgeo.org/<u></u>mailman/listinfo/gdal-dev</a><br>
<br>
</blockquote>
<br>
</blockquote></div><br></div>