[mapserver-users] WCS GetCoverage request one by one pixel

[Bakker, Anton]

Hi all,

I am running into a problem when sending HTTP requests to a MapServer backed WCS service.  I am trying to request a GeoTIFF of exactly 1 by 1 pixel, however this request fails by responding with a huge pixel (here is the describeCoverage request of the coverage requested):

https://geodata.nationaalgeoregister.nl/ahn3/wcs?service=WCS&Request=GetCoverage&version=2.0.1&CoverageId=ahn3_05m_dsm&format=image/tiff&subset=x(196396.75,196397.25)&subset=y(450031.75,450032.25)

When I modify the request to cover an area of 2X2 pixels, the response turns out fine:

https://geodata.nationaalgeoregister.nl/ahn3/wcs?service=WCS&Request=GetCoverage&version=2.0.1&CoverageId=ahn3_05m_dsm&format=image/tiff&subset=x(196396.75,196397.75)&subset=y(450031.75,450032.75)

See the screenshots hosted on imgur; (https://imgur.com/a/8IBGs3G (first screenshot is the erronous response of the request with the 1X1 extent, second screenshot is the result of the 2X2 request).

The MapServer version we are running is 7.4.3. Here is an excerpt of the mapfile that is used to configure the WCS service:

LAYER
    NAME ahn3_05m_dsm
    METADATA
      "wcs_srs"              "EPSG:28992"
      "wcs_label"           "ahn3_05m_dsm"
      "wcs_rangeset_name"   "ahn3_05m_dsm"
      "wcs_rangeset_label"  "ahn3_05m_dsm"
      "wcs_extent"          "10000 250000 356250 618750"
      "wcs_resolution"      "0.5 0.5"
      "wcs_formats"         "GEOTIFF_FLOAT32 GEOTIFF_BYTE GEOTIFF_INT16 GEOTIFF_RGB"
      "wcs_imagemode"       "FLOAT32"
      "wcs_bandcount"       "1"
    END
    TYPE RASTER
    DUMP TRUE
    STATUS ON
    DATA /srv/data/ahn3_05m_dsm/ahn3_05m_dsm.vrt
    PROCESSING "RESAMPLE=BILINEAR"
    PROJECTION
      "init=epsg:28992"
    END
END

I quickly skimmed the issue tracker on Github but could not find any issue causing this problem. Any thoughts on the cause of the problem are highly appreciated.

Kind regards,

Anton Bakker, Software Engineer
Kadaster / Geo- en Vastgoedinformatie en Advies / PDOK
tel. 06 - 29 737 351
anton.bakker at kadaster.nl
www.kadaster.nl<http://www.kadaster.nl>



Disclaimer:
De inhoud van dit bericht is uitsluitend bestemd voor geadresseerde.
Gebruik van de inhoud van dit bericht door anderen zonder toestemming van het Kadaster 
is onrechtmatig. Mocht dit bericht ten onrechte bij u terecht komen, dan verzoeken wij u 
dit direct te melden aan de verzender en het bericht te vernietigen. 
Aan de inhoud van dit bericht kunnen geen rechten worden ontleend.

Disclaimer:
The content of this message is meant to be received by the addressee only.
Use of the content of this message by anyone other than the addressee without the consent 
of the Kadaster is unlawful. If you have received this message, but are not the addressee, 
please contact the sender immediately and destroy the message.
No rights can be derived from the content of this message
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/mapserver-users/attachments/20200406/2264424f/attachment-0001.html>