using Mapserver for on-the-fly raster data?

Jeff de La Beaujardiere jeff2005 at SVS.GSFC.NASA.GOV
Thu Jul 7 07:28:42 PDT 2005


Frank Warmerdam wrote:
> I think the "simpliest" would be to have two VRT files for each target
> image, using the above method, and then another VRT file that
> mosaics them (using the SimpleSource method instead of raw bands). 

It works!! Thanks for the suggestion.  For the record, my 3 VRT files 
are reproduced below.

Now, is there any way to combine this arrangement with Mapserver's 
partial(?) support for WMS Time, Elevation and other sample Dimensions?

-Jeff DLB

-------- precpE.vrt - "East" image for longitudes [0,360]:

<VRTDataset rasterXSize="1000" rasterYSize="721">
   <SRS>EPSG:4326</SRS>
   <GeoTransform>0, 0.36, 0, 90, 0, -0.25</GeoTransform>
   <VRTRasterBand dataType="Float32" band="1" subClass="VRTRawRasterBand">
     <SourceFilename 
relativeToVRT="0">/svs/projects/map05/mapdata/private/\data/GEOS4NCEP/e32_G4ncep/Y2005/M07/D05/CNTL_NWP.e32_G4ncep_rout_prs.20050\705.000000-20050710.000000</SourceFilename>
     <ImageOffset>602756000</ImageOffset>
     <PixelOffset>4</PixelOffset>
     <LineOffset>4000</LineOffset>
     <ByteOrder>MSB</ByteOrder>
     <ColorInterp>Palette</ColorInterp>
   </VRTRasterBand>
</VRTDataset>

-------- precpW.vrt - "West" image for longitudes [-360,0]:

<VRTDataset rasterXSize="1000" rasterYSize="721">
   <SRS>EPSG:4326</SRS>
   <GeoTransform>-360, 0.36, 0, 90, 0, -0.25</GeoTransform>
   <VRTRasterBand dataType="Float32" band="1" subClass="VRTRawRasterBand">
     <SourceFilename>/path/rawbinaryfile</SourceFilename>
     <ImageOffset>602756000</ImageOffset>
     <PixelOffset>4</PixelOffset>
     <LineOffset>4000</LineOffset>
     <ByteOrder>MSB</ByteOrder>
     <ColorInterp>Palette</ColorInterp>
   </VRTRasterBand>
</VRTDataset>

-------- precp.vrt - Mosaic for longitudes [-360, 360]

<VRTDataset rasterXSize="2000" rasterYSize="721">
   <SRS>EPSG:4326</SRS>
   <GeoTransform>-360, 0.36, 0, 90, 0, -0.25</GeoTransform>
   <Simplesource>
     <SourceFilename relativeToVRT="1">precpW.vrt</SourceFilename>
     <SourceBand>1</SourceBand>
     <SrcRect xOff="0" yOff="0" xSize="1000" ySize="721"/>
     <DstRect xOff="0" yOff="0" xSize="1000" ySize="721"/>
   </Simplesource>
   <Simplesource>
     <SourceFilename relativeToVRT="1">precpE.vrt</SourceFilename>
     <SourceBand>1</SourceBand>
     <SrcRect xOff="0" yOff="0" xSize="1000" ySize="721"/>
     <DstRect xOff="1000" yOff="0" xSize="1000" ySize="721"/>
   </Simplesource>
</VRTDataset>

--------



More information about the MapServer-users mailing list