Hi,<br>I want to put together a time series of MODIS datasets in VRT, the originals being in HDF EOS format. I do that by defining the SDS for each date and sticking it into a file:<br>HDF4_EOS:EOS_GRID:&quot;../blablah.hdf&quot;:blah:blah<br>
HDF4_EOS:EOS_GRID:&quot;../blablah.hdf&quot;:blah:blah<br>HDF4_EOS:EOS_GRID:&quot;../blablah.hdf&quot;:blah:blah<br>HDF4_EOS:EOS_GRID:&quot;../blablah.hdf&quot;:blah:blah<br>[...]<br><br>I then use <br>    gdalbuildvrt -separate -input_files_list myfile.txt output_vrt.vrt<br>
<br>which creates the VRT file in output_vrt.vrt. The XML appears sensible, and all things are as they should. <br><br>&lt;VRTRasterBand dataType=&quot;Int16&quot; band=&quot;32&quot;&gt;<br>    &lt;Metadata&gt;<br>      &lt;MDI key=&quot;STATISTICS_MINIMUM&quot;&gt;-310&lt;/MDI&gt;<br>
      &lt;MDI key=&quot;STATISTICS_MAXIMUM&quot;&gt;32766&lt;/MDI&gt;<br>      &lt;MDI key=&quot;STATISTICS_MEAN&quot;&gt;321.45436944444&lt;/MDI&gt;<br>      &lt;MDI key=&quot;STATISTICS_STDDEV&quot;&gt;2990.0525001527&lt;/MDI&gt;<br>
    &lt;/Metadata&gt;<br>    &lt;NoDataValue&gt;3.27670000000000E+04&lt;/NoDataValue&gt;<br>    &lt;ComplexSource&gt;<br>      &lt;SourceFilename relativeToVRT=&quot;1&quot;&gt;HDF4_EOS:EOS_GRID:&amp;quot;../XXX&amp;quot;:XXXX:XXX&lt;/SourceFilename&gt;<br>
      &lt;SourceBand&gt;1&lt;/SourceBand&gt;<br>      &lt;SourceProperties RasterXSize=&quot;1200&quot; RasterYSize=&quot;1200&quot; DataType=&quot;Int16&quot; BlockXSize=&quot;1200&quot; BlockYSize=&quot;1&quot;/&gt;<br>
      &lt;SrcRect xOff=&quot;0&quot; yOff=&quot;0&quot; xSize=&quot;1200&quot; ySize=&quot;1200&quot;/&gt;<br>      &lt;DstRect xOff=&quot;0&quot; yOff=&quot;0&quot; xSize=&quot;1200&quot; ySize=&quot;1200&quot;/&gt;<br>      &lt;NODATA&gt;32767&lt;/NODATA&gt;<br>
    &lt;/ComplexSource&gt;<br>  &lt;/VRTRasterBand&gt;<br>  &lt;VRTRasterBand dataType=&quot;Int16&quot; band=&quot;33&quot;&gt;<br>    &lt;NoDataValue&gt;3.27670000000000E+04&lt;/NoDataValue&gt;<br>    &lt;ComplexSource&gt;<br>
      &lt;SourceFilename relativeToVRT=&quot;1&quot;&gt;HDF4_EOS:EOS_GRID:&amp;quot;../XXX&amp;quot;:XXXX:XXX&lt;/SourceFilename&gt;<br>      &lt;SourceBand&gt;1&lt;/SourceBand&gt;<br>      &lt;SourceProperties RasterXSize=&quot;1200&quot; RasterYSize=&quot;1200&quot; DataType=&quot;Int16&quot; BlockXSize=&quot;1200&quot; BlockYSize=&quot;1&quot;/&gt;<br>
      &lt;SrcRect xOff=&quot;0&quot; yOff=&quot;0&quot; xSize=&quot;1200&quot; ySize=&quot;1200&quot;/&gt;<br>      &lt;DstRect xOff=&quot;0&quot; yOff=&quot;0&quot; xSize=&quot;1200&quot; ySize=&quot;1200&quot;/&gt;<br>      &lt;NODATA&gt;32767&lt;/NODATA&gt;<br>
    &lt;/ComplexSource&gt;<br>  &lt;/VRTRasterBand&gt;<br><br><br>If I try to read this into python, using the GDAL bindings, it can deal with individual bands, but using <br><br>g = gdal.Open ( VRT_file )<br>data = g.ReadAsArray()<br>
<br>it fails claiming that it can&#39;t find the original HDF datafile. this error pops up with the band that corresponds to band number 33.<br><br>Running gdalinfo -stats on the new filename proceeds reporting the expected information, band by band... slowly but surely, until it reaches band 33 (it&#39;s band 33 everytime, even with different datasets):<br>
<br>Band 33 Block=128x128 Type=Int16, ColorInterp=Undefined               <br>ERROR 4: `HDF4_EOS:EOS_GRID:&quot;../BlahBlah.hdf&quot;:blah:blah&#39; does not exist in the file system,<br>and is not recognised as a supported dataset name.                                                                                            <br>
<br>ERROR 1: IReadBlock failed at X offset 0, Y offset 0<br><br>I can read individual bands OK (although I haven&#39;t tried looping through them), including band 33. I have stacks of memory, and data size should not be an issue. <br>
<br>Any clues?<br><br>Thanks!<br>Jose<br><br><br>