<div dir="ltr">Hi again Jukka,<div><br></div><div>Thanks for the speedy response. I now think that I have the resolution and extent correct with these. (I was getting various errors trying other things but this </div><div>set seems okay for GetCoverage, GetCapabilities and DescribeCoverage):</div><div><br></div><div><div style="color:rgb(0,0,0);font-size:14px;line-height:19px;white-space:pre"><div style="font-family:Consolas,"Courier New",monospace"> "ows_extent" "-125 24 -66 50"</div><div style="font-family:Consolas,"Courier New",monospace"> "wcs_resolution" "0.0008333333 0.0008333333"</div><div style="font-family:Consolas,"Courier New",monospace"><br></div><div style=""><font face="arial, sans-serif">I actually retrieve a correctly located GeoTIFF (yay!!) out with this query for instance:</font></div><div style="font-family:Consolas,"Courier New",monospace"><br></div><div style="font-family:Consolas,"Courier New",monospace"><a href="http://localhost:8080/wcs?service=WCS&version=2.0.1&">http://localhost:8080/wcs?service=WCS&version=2.0.1&</a>
request=GetCoverage&
coverageid=SRTM_3_arc-second_grid&
format=geotiff_16&
subset=Lat(43.256274,43.721318)&
subset=Long(-112.357549,-111.719135)<br></div><div style="font-family:Consolas,"Courier New",monospace"><br></div><div style=""><font face="arial, sans-serif">but it's all zeros. </font></div><div style=""><font face="arial, sans-serif"><br></font></div><div style=""><font face="arial, sans-serif">My data files are all native SRTM ".hgt" files. (And yes, I did "gdalinfo -stats" on one tile for the above query and there are good values.)</font></div><div style=""><font face="arial, sans-serif"><br></font></div><div style=""><font face="arial, sans-serif">I could gdal_translate them all to GeoTIFF, but since GDAL supports HGT files I hoped I could use them as is.</font></div><div style=""><font face="arial, sans-serif"><br></font></div><div style=""><font face="arial, sans-serif">I have tried everything I can think of for the "wcs_native_format" metadata:</font></div><div style=""><font face="arial, sans-serif"><br></font></div><div style=""><div style="line-height:19px"><div style="font-family:Consolas,"Courier New",monospace"> "<span style="color:rgb(163,21,21)">wcs_native_format</span>" "application/srtmhgt"</div><div style="font-family:Consolas,"Courier New",monospace"><br></div><div style=""><font face="arial, sans-serif">and application/hgt, image/hgt, srtmhgt, etc. I don't see errors in mapserver log though no matter what I try. Is this my problem?</font></div><div style=""><font face="arial, sans-serif"><br></font></div><div style=""><font face="arial, sans-serif">What do you suggest?</font></div><div style=""><font face="arial, sans-serif"><br></font></div><div style=""><font face="arial, sans-serif">Thanks again,</font></div><div style=""><font face="arial, sans-serif">carl</font></div></div>
</div><div style=""><font face="arial, sans-serif"><br></font></div><div style=""><font face="arial, sans-serif"><br></font></div><div style=""><font face="arial, sans-serif"><br></font></div><div style=""><font face="arial, sans-serif"><br></font></div><div style="font-family:Consolas,"Courier New",monospace"><br></div><div style="font-family:Consolas,"Courier New",monospace"><br></div><div style="font-family:Consolas,"Courier New",monospace"><br></div><div style="font-family:Consolas,"Courier New",monospace"><br></div></div></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Wed, Oct 20, 2021 at 3:23 PM Rahkonen Jukka (MML) <<a href="mailto:jukka.rahkonen@maanmittauslaitos.fi">jukka.rahkonen@maanmittauslaitos.fi</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><br>
<br>
Lähettäjä: MapServer-users <<a href="mailto:mapserver-users-bounces@lists.osgeo.org" target="_blank">mapserver-users-bounces@lists.osgeo.org</a>> Puolesta Carl Godkin<br>
<br>
<br>
<br>
>> On Tue, Oct 19, 2021 at 5:37 PM Rahkonen Jukka (MML) <mailto:<a href="mailto:jukka.rahkonen@maanmittauslaitos.fi" target="_blank">jukka.rahkonen@maanmittauslaitos.fi</a>> wrote:<br>
>> Hi,<br>
<br>
>> Here is more accurate documentation but only in the text <br>
>> <a href="https://www.mapserver.org/ogc/wcs_server.html" rel="noreferrer" target="_blank">https://www.mapserver.org/ogc/wcs_server.html</a><br>
<br>
>> “The convention is that once (wcs|ows)_extent and one of <br>
>> (wcs|ows)_size and (wcs|ows)_resolution is set in the layer metadata, all the coverage specific metadata will be retrieved from there. Otherwise the source image is queried via GDAL, if possible.”<br>
<br>
>> It seems that only wcs_extent is documented in the list of metadata elements but this is what they do:<br>
>> • wcs_extent defines the bounding box or your coverage -> that will <br>
>> go into DescribeCoverage • wcs_size is the size of the coverage as <br>
>> pixels -> pixel size can be computed by extent and size • <br>
>> wcs_resolution tells the pixel size explicitly -> size of the <br>
>> coverage in pixels can be computed<br>
<br>
> Thanks a lot, Jukka. That all makes some sense as far as it goes, but <br>
> I'm still not clear on a few details. I have this range of data: 24N to 50N and 125W to 66W which is 26 degrees by 59 degrees.<br>
<br>
> Since the SRTM data has 3 arc-second spacing, that's 1200 pixels per <br>
> degree plus one for the extra edge so I have tried a number of things without success.<br>
<br>
What fails? Don't you get anything into DescribeCoverage or is all that you get wrong?<br>
<br>
> 1. Is wcs_size the size FULL range of the tiles that I have? And what is the order?<br>
<br>
You do not need to bother about the size if you set the resolution, but yes, it is the full range of the coverage. See this <a href="https://beta-karttakuva.maanmittauslaitos.fi/wcs/service/ows?service=WCS&REQUEST=DescribeCoverage&VERSION=2.0.1&COVERAGEID=ortokuva__ortokuva" rel="noreferrer" target="_blank">https://beta-karttakuva.maanmittauslaitos.fi/wcs/service/ows?service=WCS&REQUEST=DescribeCoverage&VERSION=2.0.1&COVERAGEID=ortokuva__ortokuva</a><br>
<br>
The size is converted into the pixel space of the coverage gml:GridEnvelope><br>
<gml:low>0 0</gml:low><br>
<gml:high>1348356 2316499</gml:high><br>
</gml:GridEnvelope<br>
<br>
I would say that the order is left-to-right - down-to-up but I fear that by the GML standard it can be something else. You can read about rectified grid from the OpenGIS Geography Markup Language (GML) Encoding Standard 3.2.1. I have never understood it but there are some images which may clarify this: "rectified grid: grid for which there is an affine transformation between the grid coordinates and the coordinates of an external coordinate reference system"<br>
<br>
> 2. Again, what is the order for wcs_resolution? Are the units degrees since the coordinate system of the original data is EPSG:4326?<br>
With Mapserver the order is easting/longitude - northing/latitude<br>
<br>
> I've tried all of these, one at a time, and get various errors. (1 /<br>
> 1200 = 0.00083333...)<br>
What errors?<br>
<br>
> #"wcs_resolution" "0.00083333 0.00083333"<br>
Looks right<br>
> #"wcs_size" "70801 31201"<br>
Looks right if you prefer to give size instead of resolution<br>
> "wcs_size" "31201 70801"<br>
This is wrong.<br>
<br>
> In case it's helpful, here's my LAYER definition (with probably extra things I've tried based on my searching and trial & error):<br>
<br>
LAYER<br>
NAME SRTM_3_arc-second_grid<br>
METADATA<br>
"wcs_label" "SRTM - 3 arc-second grid" ### required<br>
"wcs_rangeset_name" "Range 1" ### required to support DescribeCoverage request<br>
"wcs_rangeset_label" "Lower 48" ### required to support DescribeCoverage request<br>
"ows_extent" "-125 24 -66 50"<br>
#"wcs_resolution" "0.00083333 0.00083333" # 1/1200 for 3" spacing?<br>
#"wcs_size" "70801 31201" # Not sure of order<br>
"wcs_size" "31201 70801"<br>
"wcs_bandcount" "1"<br>
"wcs_native_format" "SRTMHGT" # What should this be? I can't find examples...<br>
<br>
It is supposed to be the native format of the layer data, format that theoretically requires no processing, just selecting the right pixels. In practice the format that you users get by default with GetCoverage. It goes into DescribeCoverage about this way:<br>
<wcs:nativeFormat>image/tiff</wcs:nativeFormat> <br>
and it is documented in the WCS 2.0.1 Core standard:<br>
"4.7 Native Format<br>
encoding format where, in a GetCoverage request, the range set values can be obtained unaltered"<br>
<br>
"wcs_srs" "EPSG:4326"<br>
END<br>
<br>
TYPE RASTER ### required<br>
STATUS ON<br>
TILEINDEX "../srtm_3_hgt/srtm_3_hgt-index.shp" # Path is relative to SHAPEPATH<br>
TILEITEM "Location"<br>
PROJECTION<br>
"init=epsg:4326"<br>
END<br>
END<br>
<br>
> I feel like I must be close. Thanks a lot for any further pointers,<br>
<br>
-Jukka-<br>
<br>
> carl<br>
<br>
</blockquote></div>