<div dir="ltr">Hey dear list,<div>this is my first post here, </div><div>so sorry if I'm asking something obvious.</div><div><br></div><div>(I noticed the 2 previous mail thread regarding python and VRT).</div><div><br></div><div>I wrote python functions to convert a numpy structured array to several gdal rasters.</div><div><br></div><div>A numpy structured array is a 2D pixel array, but each pixel has in fact several values of possibly different type (like band).</div><div>For instance a pixel dtype can be [float32,uint16,Float64].</div><div><br></div><div>I want to convert this array to raster(s).</div><div>First I tried to create one raster with several bands,</div><div>but GeoTiff cannot handle bands with different pixel type (or I didn't find how).</div><div><br></div><div>I concluded that I must create several rasters, and put them all as bands in a VRT.</div><div><br></div><div>So far I have been able to create several geoTiff raster of various type corresponding to the numpy structured array.</div><div><br></div><div>Now, I don't know how to create from python the VRT that unites all this rasters.</div><div><br></div><div>When I try to write in a VRT raster band I get errors.</div><div><br></div><div>I am well aware that I could create the XML directly from python, but it feels like a pretty bad architecture idea.</div><div><br></div><div>Anybody knows how to do this /and/or has comments about this task?</div><div><br></div><div>Many Thanks,</div><div>Rémi-C</div><div><br></div><div><br></div></div>