[mapserver-users] wcs 1.0.0 bbox problem

Wed Jan 29 00:32:37 PST 2014

Hi,

I wonder where BBOX (-2.5,-2.5,2.5,2.5) comes. If Mapserver is generating it automatically then it feels like a bug.

Your request does not contain BBOX at all and it is OK because it has TIME, and WCS 1.0.0 standard says  "One of BBOX or TIME is required".  CRS is required but because there is no BBOX it should not have any effect on the result. Standard says here
"GetCoverage requests must use this parameter to specify the coordinate reference system in which the request domain constraints are expressed (BBOX). The values of this re-quest parameter must be one of those defined in a requestResponseCRSs or requestCRSs element under the requested coverage." Thus you are allowed to use any CRS that your service supports and the result should be correct and identical. And your CRS list seems to be this:
"wcs_srs" "EPSG:4326 EPSG:900913 EPSG:3857 EPSG:32615 EPSG:32616 EPSG:32617"

Do you get the same error also if you make the request with for example &CRC=900913?
I guess that your system is sometimes generating requests with EPSG:4326 BBOX and therefore you just that as a compulsory dummy CRS parameter with your bare TIME requests and therefore using the native EPSG:32626 is not a good workaround.

-Jukka Rahkonen-

Lähettäjä: mapserver-users-bounces at lists.osgeo.org [mailto:mapserver-users-bounces at lists.osgeo.org] Puolesta Peter Freimuth
Lähetetty: 29. tammikuuta 2014 9:17
Vastaanottaja: mapserver-users at lists.osgeo.org
Aihe: [mapserver-users] wcs 1.0.0 bbox problem

Hi,
executing the following request
https://devcloud.blackbridge.com/wcs/PlenumsoftDemo?SERVICE=WCS&VERSION=1.0.0&REQUEST=GetCoverage&COVERAGE=testdata_raw&FORMAT=GEOTIFF_INT16&RESX=5&RESY=5&RESPONSE_CRS=EPSG:32616&CRS=EPSG:4326&TIME=2012-04-05

results in

<?xml version='1.0' encoding="ISO-8859-1" ?>
xmlns="http://www.opengis.net/ogc" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" xsi:schemaLocation="http://www.opengis.net/ogc http://schemas.opengeospatial.net/wcs/1.0.0/OGC-exception.xsd">
msWCSGetCoverage(): WCS server error. Requested BBOX (-2.5,-2.5,2.5,2.5) is outside requested coverage BBOX (210500,2230500,261500,2281500)
msProjectRect(): Projection library error. All points failed to reproject.
msProjectPoint(): Projection library error. proj says: latitude or longitude exceeded limits

The corresponding layer d3efinition is:

MAP
...
LAYER
# this is the tileindex that aggregates the single image products together and makes them accessible by date and extend
NAME "data_raw_idx"
STATUS ON

DEBUG 5
TYPE POLYGON
DATA "/mnt/mapserver/mapserver_wms/Raster/Plenumosft_RE_Sample_forWCS/plenumsoft_demo_data"
PROJECTION
"init=epsg:32616"
END
VALIDATION
# 'productid' '^([0-9])$'
END
#FILTER ( '[productid]'='%productid%')
METADATA
#'productid_validation_pattern' '^([0-9])$'
#'default_productid' '1234567'
#disable explicitly everything that shall not be supported for this Layer
"wms_enable_request" "!GetCapabilities !GetMap !GetLegendgraphic"
"wcs_enable_request" "!GetCapabilities !GetCoverage !DescribeCoverage"
"wfs_enable_request" "!GetCapabilities !GetFeature !DescribeFeatureType"
END
END
LAYER
# this layer serves the 16bit raw data
NAME "testdata_raw"
STATUS ON

DUMP TRUE
DEBUG 5
TYPE RASTER ### required
PROJECTION
"init=epsg:32616"
END
OFFSITE 0 0 0
UNITS meters
TILEINDEX "data_raw_idx"
TILEITEM "location"
#MASK "mask2"
METADATA
"wcs_title" "5 band 16bit RE Testdata with time support"
"wcs_extent" "210500.0 2230500.0 261500.0 2281500.0" # shall be set dynamically by the wrapper script on mapObj base
"wcs_srs" "EPSG:4326 EPSG:900913 EPSG:3857 EPSG:32615 EPSG:32616 EPSG:32617"
"wcs_label" "5 band 16bit RE Testdata" ### required
"wcs_rangeset_axes" "bands"
"wcs_rangeset_name" "Bands" ### required to support DescribeCoverage request
"wcs_rangeset_label" "RE Radiometric Bands" ### required to support DescribeCoverage request
"wcs_rangeset_nullvalue" "0"

#"wcs_size" "15000 15000" #if we have a global ows_extent set in web-metadata and a wcs_resolution this gets calculated by mapserver
"wcs_resolution" "5.0 5.0"
"wcs_bandcount" "5"
"wcs_band_names" "Blue Green Red Rededge Nir"
"wcs_imagemode" "INT16"
"wcs_formats" "GEOTIFF_INT16 AAIGRID"
"wcs_native_format" "image/tiff"
"wcs_nativeformat" "GTiff"
"wcs_enable_request" "*"
"wcs_timeposition" "2011-01-04,2011-04-20,2012-04-05" #this might get set by the mapscript wrapper accodingly to the dates avaialable for the service and is mandatory to allow temporal subsettings
"wcs_timeextent" "2011-01-04/2012-04-05" # shall be set dynamically by the wrapper script on mapObj base
"wcs_timeitem" "acquidate"
#disable explicitly everything that shall not be supported for this Layer
"wms_enable_request" "!GetCapabilities !GetMap !GetLegendgraphic"
"wfs_enable_request" "!GetCapabilities !GetFeature !DescribeFeatureType"
END
CLASS
NAME "5 Band (16bit) RE Imagedata"
KEYIMAGE "drgs_keyimage.png"
END
END
...
END

If i use CRS=EPSG:32616 which is the epsg code of the layer itself, everything works fine. I would expect that mapserver would reproject the wcs_extent or layer_extent or map_extent to the given CRS and use that BBOX by default.
related to:
https://github.com/mapserver/mapserver/issues/4853

--
Peter Freimuth

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/mapserver-users/attachments/20140129/d1fd896a/attachment.htm>