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