<div dir="ltr">Hi Jukka,<div><br></div><div>Thank you for your reply. I will play around with WMS + GetFeatureInfo and see if I can make that work. </div><div><br></div><div>I'm still pretty new to this stuff and WCS seemed the more natural service for this. </div><div><br></div><div>Thanks,</div><div><br></div><div>carl</div><div><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Mon, Nov 1, 2021 at 12:58 AM 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">
<div lang="FI" style="overflow-wrap: break-word;">
<div class="gmail-m_1315811357735716798WordSection1">
<p class="MsoNormal"><span>Hi,<u></u><u></u></span></p>
<p class="MsoNormal"><span><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">The WMS GetFeatureInfo should work and the answer in gis.stackexhange seems to contain all the details that you need for making it to work.<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">With WCS what you did is probably the best you can do with Mapserver. With some other servers you may get a single pixel output by using WCS slicing (subset=Long(-100)&subset=Lat(22))
or trimming with lowerlimit=upperlimit (subset=Long(-100,-100)&subset=Lat(22,22)) but Mapserver does not support either. See
<a href="https://osgeo-org.atlassian.net/browse/GEOS-9553" target="_blank">https://osgeo-org.atlassian.net/browse/GEOS-9553</a>.<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">I could not find a link to an answer by Stefan Meissl about why Mapserver does not support slicing but it was about dimensions. A raster image can be considered to be two dimensional
even when it has only one pixel, but by the WCS standard each trim is removing one dimension. Slicing DEM with two axes must return a zero dimensional result and raster formats do not support that. I still think that subsets with lower=upper limits should
works and return a single pixel output but that does not work with either Mapserver or Geoserver. See also
<a href="https://www.mail-archive.com/search?l=mapserver-users%40lists.osgeo.org&q=subject:%22%5C%5Bmapserver%5C-users%5C%5D+WCS+GetCoverage+request+one++by+one+pixel%22&o=newest&f=1" target="_blank">
https://www.mail-archive.com/search?l=mapserver-users%40lists.osgeo.org&q=subject:%22%5C%5Bmapserver%5C-users%5C%5D+WCS+GetCoverage+request+one++by+one+pixel%22&o=newest&f=1</a><u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">-Jukka Rahkonen-<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"><u></u> <u></u></span></p>
<div style="border-right:none;border-bottom:none;border-left:none;border-top:1pt solid rgb(225,225,225);padding:3pt 0cm 0cm">
<p class="MsoNormal"><b>Lähettäjä:</b> MapServer-users <<a href="mailto:mapserver-users-bounces@lists.osgeo.org" target="_blank">mapserver-users-bounces@lists.osgeo.org</a>>
<b>Puolesta </b>Carl Godkin<br>
<b>Lähetetty:</b> lauantai 30. lokakuuta 2021 2.10<br>
<b>Vastaanottaja:</b> <a href="mailto:mapserver-users@lists.osgeo.org" target="_blank">mapserver-users@lists.osgeo.org</a><br>
<b>Aihe:</b> [mapserver-users] Point query from Mapserver raster<u></u><u></u></p>
</div>
<p class="MsoNormal"><u></u> <u></u></p>
<div>
<p class="MsoNormal">Hi,<u></u><u></u></p>
<div>
<p class="MsoNormal">Now that I have my WMS and WCS servers working well (thanks, list!!), I have one more thing I would like to figure out how to do.<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal">I want to do point queries on [Lat,Lon] points for their elevation from the terrain layer I have.<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal">I have done some searching and come across various ideas such as Raster Query (<a href="https://mapserver.org/input/raster.html#raster-query" target="_blank">https://mapserver.org/input/raster.html#raster-query</a>) and a query template (<a href="https://gis.stackexchange.com/questions/167239/mapserver-raster-picking-elevation-data" target="_blank">https://gis.stackexchange.com/questions/167239/mapserver-raster-picking-elevation-data</a>)
but haven't found enough information for either to implement it.<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal">I came up with a solution that works fine, but might be considered as going "All Around the Barn" What works is to add an ASCII output format to my map file that serves terrain like this:<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<div>
<div>
<p class="MsoNormal" style="line-height:14.25pt"><span style="font-size:10.5pt;font-family:Consolas;color:black">
</span><span style="font-size:10.5pt;font-family:Consolas;color:rgb(121,94,38)">OUTPUTFORMAT</span><span style="font-size:10.5pt;font-family:Consolas;color:black"><u></u><u></u></span></p>
</div>
<div>
<p class="MsoNormal" style="line-height:14.25pt"><span style="font-size:10.5pt;font-family:Consolas;color:black">
</span><span style="font-size:10.5pt;font-family:Consolas;color:blue">NAME</span><span style="font-size:10.5pt;font-family:Consolas;color:black"> XYZ<u></u><u></u></span></p>
</div>
<div>
<p class="MsoNormal" style="line-height:14.25pt"><span style="font-size:10.5pt;font-family:Consolas;color:black">
</span><span style="font-size:10.5pt;font-family:Consolas;color:blue">DRIVER</span><span style="font-size:10.5pt;font-family:Consolas;color:black"> "GDAL/XYZ"<u></u><u></u></span></p>
</div>
<div>
<p class="MsoNormal" style="line-height:14.25pt"><span style="font-size:10.5pt;font-family:Consolas;color:black">
</span><span style="font-size:10.5pt;font-family:Consolas;color:blue">MIMETYPE</span><span style="font-size:10.5pt;font-family:Consolas;color:black"> "text/plain"<u></u><u></u></span></p>
</div>
<div>
<p class="MsoNormal" style="line-height:14.25pt"><span style="font-size:10.5pt;font-family:Consolas;color:black">
</span><span style="font-size:10.5pt;font-family:Consolas;color:blue">IMAGEMODE</span><span style="font-size:10.5pt;font-family:Consolas;color:black">
</span><span style="font-size:10.5pt;font-family:Consolas;color:rgb(0,16,128)">FLOAT32</span><span style="font-size:10.5pt;font-family:Consolas;color:black"><u></u><u></u></span></p>
</div>
<div>
<p class="MsoNormal" style="line-height:14.25pt"><span style="font-size:10.5pt;font-family:Consolas;color:black">
</span><span style="font-size:10.5pt;font-family:Consolas;color:blue">EXTENSION</span><span style="font-size:10.5pt;font-family:Consolas;color:black"> "txt"<u></u><u></u></span></p>
</div>
<div>
<p class="MsoNormal" style="line-height:14.25pt"><span style="font-size:10.5pt;font-family:Consolas;color:black">
</span><span style="font-size:10.5pt;font-family:Consolas;color:rgb(121,94,38)">END</span><span style="font-size:10.5pt;font-family:Consolas;color:black"><u></u><u></u></span></p>
</div>
<div>
<p class="MsoNormal" style="line-height:14.25pt"><span style="font-size:10.5pt;font-family:Consolas;color:black"><u></u> <u></u></span></p>
</div>
<div>
<p class="MsoNormal" style="line-height:14.25pt"><span style="font-size:10.5pt;font-family:Arial,sans-serif;color:rgb(121,94,38)">and then use GetCoverage on a small rectangle whose lower-left corner
</span><span style="font-size:10.5pt;color:black"><u></u><u></u></span></p>
</div>
<div>
<p class="MsoNormal" style="line-height:14.25pt"><span style="font-size:10.5pt;font-family:Arial,sans-serif;color:rgb(121,94,38)">is the point whose elevation I want:</span><span style="font-size:10.5pt;color:black"><u></u><u></u></span></p>
</div>
<div>
<p class="MsoNormal" style="line-height:14.25pt"><span style="font-size:10.5pt;color:black"><u></u> <u></u></span></p>
</div>
<div>
<p class="MsoNormal" style="line-height:14.25pt"><span style="font-size:9pt;font-family:"Courier New";color:rgb(33,33,33)"><a href="http://localhost:8080/wcs?SERVICE=WCS&VERSION=2.0.1&REQUEST=GetCoverage&coverageid=SRTM_3_arc-second_grid&format=XYZ&subset=Lat(42.948,42.958)&subset=Lon(-122.086,-122.076)" target="_blank">http://localhost:8080/wcs?SERVICE=WCS&VERSION=2.0.1&REQUEST=GetCoverage&coverageid=SRTM_3_arc-second_grid&format=XYZ&subset=Lat(42.948,42.958)&subset=Lon(-122.086,-122.076)</a></span><span style="font-size:10.5pt;color:black"><u></u><u></u></span></p>
</div>
<div>
<p class="MsoNormal" style="line-height:14.25pt"><span style="font-size:10.5pt;color:black"><u></u> <u></u></span></p>
</div>
<div>
<p class="MsoNormal" style="line-height:14.25pt"><span style="font-family:Arial,sans-serif;color:rgb(33,33,33)">and then take the first line of the returned output file.</span><span style="color:black"><u></u><u></u></span></p>
</div>
<div>
<p class="MsoNormal" style="line-height:14.25pt"><span style="color:black"><u></u> <u></u></span></p>
</div>
<div>
<p class="MsoNormal" style="line-height:14.25pt"><span style="font-family:Arial,sans-serif;color:rgb(33,33,33)">What would Mapserver best practice be for a problem like this?</span><u></u><u></u></p>
</div>
<div>
<p class="MsoNormal" style="line-height:14.25pt"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal" style="line-height:14.25pt"><span style="font-family:Arial,sans-serif;color:rgb(33,33,33)">Thanks very much,</span><u></u><u></u></p>
</div>
<div>
<p class="MsoNormal" style="line-height:14.25pt"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal" style="line-height:14.25pt"><span style="font-family:Arial,sans-serif;color:rgb(33,33,33)">carl</span><u></u><u></u></p>
</div>
<div>
<p class="MsoNormal" style="line-height:14.25pt"><span style="color:black"><u></u> <u></u></span></p>
</div>
<div>
<p class="MsoNormal" style="line-height:14.25pt"><span style="color:black"><u></u> <u></u></span></p>
</div>
</div>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
</div>
</div>
</div>
</blockquote></div>