[mapserver-users] Continuing problem with WCS and ENVI Multi-band files

[Peter Willis]

Hello,

This has actually been resolved. I just haven't got back here to
give an update until now.

The proper mapserver use case for ENVI type rasters and WCS
will be posted at mapserver.org by Frank Warmerdam once I
get my 'stuff' together.

There are some small sample ENVI test files available at:

http://filebin.ca/cpaux/ENVI_test_files.ZIP

NOTE: the files listed as ...int32...
are actually  ...int16... and the file name
is a typo.

These files (content) are not scientific data files
but 6 band binary files with headers.
The image data is a text burn in of the
band number in each band on a zero background.

The problem I was having with my WCS dataset
LAYER definition seemed to clear up with some
minor editing.

There is a sample map file for UMN mapserver
and ENVI format WCS located at:

http://filebin.ca/upkmxs/mapfile.txt

The geographic pixel size in the map file is incorrect.
It should be 2.8125 degrees not 0.08333. The file is only 128x64 spatially.

the band selection for GetCoverage URL is &bands=<number>

ie:(see end of URL line)

http://somewcsserver.net/cgi-bin/wcs?REQUEST=GetCoverage&SERVICE=WCS&VERSION=1.0.0&COVERAGE=BOGUStestData&CRS=EPSG:4326&BBOX=-180,-90,180,90&WIDTH=128&HEIGHT=64&FORMAT=GEOTIFF_FLOAT&bands=6

Hope that's useful.

If you need NDVI that's a whole other story. :)

Peter


Lime, Steve D (DNR) wrote:
> Any chance of getting a test dataset?
> 
> -----Original Message-----
> From: mapserver-users-bounces at lists.osgeo.org [mailto:mapserver-users-bounces at lists.osgeo.org] On Behalf Of Peter Willis
> Sent: Monday, February 08, 2010 12:37 PM
> To: mapserver-users at lists.osgeo.org
> Subject: [mapserver-users] Continuing problem with WCS and ENVI Multi-band files
> 
> Hello,
> 
> I continue to have a problem selecting
> a single band from a multi-band (563 bands)
> ENVI format BSQ file.
> 
> Using the following URL mapserver gives me a FLOAT32 GTiff
> file but it is always filled with zeros.
> 
> http://wcs.foo.com/cgi-bin/wcs?REQUEST=GetCoverage&SERVICE=WCS&VERSION=1.0.0&COVERAGE=NDVI&CRS=EPSG:4326&BBOX=-135,55,-121,46&WIDTH=256&HEIGHT=256&FORMAT=GEOTIFF_FLOAT&Band=122
> 
> Is there any documentation that clearly
> shows how WCS can be set up to do this?
> 
> My current LAYER definition is as follows:
> 
> LAYER
>    NAME VegetationIndex
>    STATUS OFF
>    DEBUG ON
>    TYPE RASTER
>    METADATA
>      "wcs_label"           "Data/NDVI"
>      "wcs_rangeset_name"   'bands'
>      "wcs_rangeset_label"  "NDVI"
>      "ows_extent" '-135 55 -121 46'
>      "wcs_resolution" '0.08333333333333332 -0.083333333333333332'
>      "ows_srs" 'EPSG:4326'
>      "wcs_srs" 'EPSG:4326'
>      "wcs_formats" 'GEOTIFF_FLOAT,GEOTIFF_INT16'
>      "wcs_nativeformat" 'ENVI'
> 
>    END
>    DATA /NDVI/ndvi.img
>    PROJECTION
>      "init=epsg:4623"
>    END
>    DUMP TRUE
> END