<html xmlns:v="urn:schemas-microsoft-com:vml" xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40">
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
<meta name="Generator" content="Microsoft Word 15 (filtered medium)">
<!--[if !mso]><style>v\:* {behavior:url(#default#VML);}
o\:* {behavior:url(#default#VML);}
w\:* {behavior:url(#default#VML);}
.shape {behavior:url(#default#VML);}
</style><![endif]--><style><!--
/* Font Definitions */
@font-face
        {font-family:Helvetica;
        panose-1:2 11 6 4 2 2 2 2 2 4;}
@font-face
        {font-family:"Cambria Math";
        panose-1:2 4 5 3 5 4 6 3 2 4;}
@font-face
        {font-family:Calibri;
        panose-1:2 15 5 2 2 2 4 3 2 4;}
@font-face
        {font-family:Consolas;
        panose-1:2 11 6 9 2 2 4 3 2 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
        {margin:0cm;
        margin-bottom:.0001pt;
        font-size:11.0pt;
        font-family:"Calibri",sans-serif;}
a:link, span.MsoHyperlink
        {mso-style-priority:99;
        color:blue;
        text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
        {mso-style-priority:99;
        color:purple;
        text-decoration:underline;}
p.msonormal0, li.msonormal0, div.msonormal0
        {mso-style-name:msonormal;
        mso-margin-top-alt:auto;
        margin-right:0cm;
        mso-margin-bottom-alt:auto;
        margin-left:0cm;
        font-size:11.0pt;
        font-family:"Calibri",sans-serif;}
p.xmsonormal, li.xmsonormal, div.xmsonormal
        {mso-style-name:x_msonormal;
        margin:0cm;
        margin-bottom:.0001pt;
        font-size:11.0pt;
        font-family:"Calibri",sans-serif;}
span.Shkpostityyli19
        {mso-style-type:personal-reply;
        font-family:"Calibri",sans-serif;
        color:windowtext;}
.MsoChpDefault
        {mso-style-type:export-only;
        font-size:10.0pt;}
@page WordSection1
        {size:612.0pt 792.0pt;
        margin:70.85pt 70.85pt 70.85pt 70.85pt;}
div.WordSection1
        {page:WordSection1;}
--></style><!--[if gte mso 9]><xml>
<o:shapedefaults v:ext="edit" spidmax="1026" />
</xml><![endif]--><!--[if gte mso 9]><xml>
<o:shapelayout v:ext="edit">
<o:idmap v:ext="edit" data="1" />
</o:shapelayout></xml><![endif]-->
</head>
<body lang="FI" link="blue" vlink="purple">
<div class="WordSection1">
<p class="MsoNormal"><span style="mso-fareast-language:EN-US">Hi Anton,<o:p></o:p></span></p>
<p class="MsoNormal"><span style="mso-fareast-language:EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US" style="mso-fareast-language:EN-US">What I meant was to have a try by commenting out this line from your mapfile<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="NL" style="font-size:10.0pt;font-family:Consolas">    PROCESSING "RESAMPLE=BILINEAR"<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="NL" style="font-size:10.0pt;font-family:Consolas"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="NL" style="font-size:10.0pt;font-family:Consolas">Actually it would be interesting to know what really happens if the interpolation/resampling method is configured in the mapfile and GetCoverage comes also with Interpolation
 parameter.<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="NL" style="font-size:10.0pt;font-family:Consolas"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US" style="mso-fareast-language:EN-US">-Jukka-<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US" style="mso-fareast-language:EN-US"><o:p> </o:p></span></p>
<div>
<div style="border:none;border-top:solid #E1E1E1 1.0pt;padding:3.0pt 0cm 0cm 0cm">
<p class="MsoNormal"><b>Lähettäjä:</b> Bakker, Anton <Anton.Bakker@kadaster.nl> <br>
<b>Lähetetty:</b> tiistai 7. huhtikuuta 2020 18.40<br>
<b>Vastaanottaja:</b> Rahkonen Jukka (MML) <jukka.rahkonen@maanmittauslaitos.fi>; mapserver-users@lists.osgeo.org<br>
<b>Aihe:</b> RE: WCS GetCoverage request one by one pixel<o:p></o:p></p>
</div>
</div>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"><span lang="NL" style="font-size:10.0pt;font-family:"Arial",sans-serif;color:black">Hi Jukka,<o:p></o:p></span></p>
<div>
<p class="MsoNormal"><span lang="NL" style="font-size:10.0pt;font-family:"Arial",sans-serif;color:black"><o:p> </o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span lang="NL" style="font-size:10.0pt;font-family:"Arial",sans-serif;color:black">Thanks for your response; do you mean to set the default resampling method in the mapfile? Setting the resampling method in the WCS request does not work
 (for both GET and POST requests):<o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span lang="NL" style="font-size:10.0pt;font-family:"Arial",sans-serif;color:black"><o:p> </o:p></span></p>
</div>
<div>
<p class="MsoNormal" style="margin-bottom:12.0pt"><span lang="NL" style="font-size:10.0pt;font-family:"Arial",sans-serif;color:black">POST:<o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span lang="NL" style="font-size:10.0pt;font-family:"Arial",sans-serif;color:black"><?xml version="1.0" encoding="UTF-8"?><o:p></o:p></span></p>
<div>
<p class="MsoNormal"><span lang="NL" style="font-size:10.0pt;font-family:"Arial",sans-serif;color:black"><wcs:GetCoverage<o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span lang="NL" style="font-size:10.0pt;font-family:"Arial",sans-serif;color:black">  xmlns:xsi='http://www.w3.org/2001/XMLSchema-instance'<o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span lang="NL" style="font-size:10.0pt;font-family:"Arial",sans-serif;color:black">  xsi:schemaLocation="http://www.opengis.net/wcs/2.0<o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span lang="NL" style="font-size:10.0pt;font-family:"Arial",sans-serif;color:black">   
<a href="http://schemas.opengis.net/wcs/2.0/wcsAll.xsd">http://schemas.opengis.net/wcs/2.0/wcsAll.xsd</a>"<o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span lang="NL" style="font-size:10.0pt;font-family:"Arial",sans-serif;color:black">  xmlns="<a href="http://www.opengis.net/wcs/2.0">http://www.opengis.net/wcs/2.0</a>"<o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span lang="NL" style="font-size:10.0pt;font-family:"Arial",sans-serif;color:black">  xmlns:wcs="<a href="http://www.opengis.net/wcs/2.0">http://www.opengis.net/wcs/2.0</a>"<o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span lang="NL" style="font-size:10.0pt;font-family:"Arial",sans-serif;color:black">  xmlns:wcscrs="<a href="http://www.opengis.net/wcs/crs/1.0">http://www.opengis.net/wcs/crs/1.0</a>"<o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span lang="NL" style="font-size:10.0pt;font-family:"Arial",sans-serif;color:black">  xmlns:scal="<a href="http://www.opengis.net/wcs/scaling/1.0">http://www.opengis.net/wcs/scaling/1.0</a>"<o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span lang="NL" style="font-size:10.0pt;font-family:"Arial",sans-serif;color:black">  xmlns:int="<a href="http://www.opengis.net/wcs/interpolation/1.0">http://www.opengis.net/wcs/interpolation/1.0</a>"<o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span lang="NL" style="font-size:10.0pt;font-family:"Arial",sans-serif;color:black">  service="WCS"<o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span lang="NL" style="font-size:10.0pt;font-family:"Arial",sans-serif;color:black">  version="2.0.1"><o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span lang="NL" style="font-size:10.0pt;font-family:"Arial",sans-serif;color:black">    <wcs:CoverageId>ahn3_05m_dsm</wcs:CoverageId><o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span lang="NL" style="font-size:10.0pt;font-family:"Arial",sans-serif;color:black">    <wcs:DimensionTrim><o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span lang="NL" style="font-size:10.0pt;font-family:"Arial",sans-serif;color:black">        <wcs:Dimension>x</wcs:Dimension><o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span lang="NL" style="font-size:10.0pt;font-family:"Arial",sans-serif;color:black">        <wcs:TrimLow>196396.75</wcs:TrimLow><o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span lang="NL" style="font-size:10.0pt;font-family:"Arial",sans-serif;color:black">        <wcs:TrimHigh>196396.75</wcs:TrimHigh><o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span lang="NL" style="font-size:10.0pt;font-family:"Arial",sans-serif;color:black">    </wcs:DimensionTrim><o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span lang="NL" style="font-size:10.0pt;font-family:"Arial",sans-serif;color:black">    <wcs:DimensionTrim><o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span lang="NL" style="font-size:10.0pt;font-family:"Arial",sans-serif;color:black">        <wcs:Dimension>y</wcs:Dimension><o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span lang="NL" style="font-size:10.0pt;font-family:"Arial",sans-serif;color:black">        <wcs:TrimLow>450031.75</wcs:TrimLow><o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span lang="NL" style="font-size:10.0pt;font-family:"Arial",sans-serif;color:black">        <wcs:TrimHigh>450031.75</wcs:TrimHigh><o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span lang="NL" style="font-size:10.0pt;font-family:"Arial",sans-serif;color:black">    </wcs:DimensionTrim><o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span lang="NL" style="font-size:10.0pt;font-family:"Arial",sans-serif;color:black">    <wcs:format>image/tiff</wcs:format><o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span lang="NL" style="font-size:10.0pt;font-family:"Arial",sans-serif;color:black">    <wcs:Extension><o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span lang="NL" style="font-size:10.0pt;font-family:"Arial",sans-serif;color:black">        <int:Interpolation><o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span lang="NL" style="font-size:10.0pt;font-family:"Arial",sans-serif;color:black">            <int:globalInterpolation>NEAREST</int:globalInterpolation><o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span lang="NL" style="font-size:10.0pt;font-family:"Arial",sans-serif;color:black">        </int:Interpolation><o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span lang="NL" style="font-size:10.0pt;font-family:"Arial",sans-serif;color:black">    </wcs:Extension><o:p></o:p></span></p>
</div>
<p class="MsoNormal"><span lang="NL" style="font-size:10.0pt;font-family:"Arial",sans-serif;color:black"></wcs:GetCoverage><o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span lang="NL" style="font-size:10.0pt;font-family:"Arial",sans-serif;color:black"><o:p> </o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span lang="NL" style="font-size:10.0pt;font-family:"Arial",sans-serif;color:black">GET: <br>
<br>
</span><span lang="NL" style="font-size:9.0pt;font-family:"Helvetica",sans-serif;color:#505050;background:white"><a href="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,196396.75)&subset=y(450031.75,450031.75)&INTERPOLATION=NEAREST">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,196396.75)&subset=y(450031.75,450031.75)&INTERPOLATION=NEAREST</a><br>
</span><span lang="NL" style="font-size:10.0pt;font-family:"Arial",sans-serif;color:black"><br>
Both POST and GET request fail with the same error: "msImageCreate(): Image handling error. Attempt to allocate raw image failed, out of memory." <o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span lang="NL" style="font-size:10.0pt;font-family:"Arial",sans-serif;color:black"><o:p> </o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span lang="NL" style="font-size:10.0pt;font-family:"Arial",sans-serif;color:black">Cheers,<o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span lang="NL" style="font-size:10.0pt;font-family:"Arial",sans-serif;color:black">Anton<o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span lang="NL" style="font-size:10.0pt;font-family:"Arial",sans-serif;color:black"><o:p> </o:p></span></p>
</div>
<div class="MsoNormal" align="center" style="text-align:center"><span lang="NL">
<hr size="2" width="98%" align="center">
</span></div>
<div id="divRplyFwdMsg">
<p class="MsoNormal"><b><span lang="NL" style="color:black">From:</span></b><span lang="NL" style="color:black"> Rahkonen Jukka (MML) <<a href="mailto:jukka.rahkonen@maanmittauslaitos.fi">jukka.rahkonen@maanmittauslaitos.fi</a>><br>
<b>Sent:</b> Monday, April 6, 2020 5:26 PM<br>
<b>To:</b> Bakker, Anton <<a href="mailto:Anton.Bakker@kadaster.nl">Anton.Bakker@kadaster.nl</a>>;
<a href="mailto:mapserver-users@lists.osgeo.org">mapserver-users@lists.osgeo.org</a> <<a href="mailto:mapserver-users@lists.osgeo.org">mapserver-users@lists.osgeo.org</a>><br>
<b>Subject:</b> Re: WCS GetCoverage request one by one pixel</span><span lang="NL">
<o:p></o:p></span></p>
<div>
<p class="MsoNormal"><span lang="NL"> <o:p></o:p></span></p>
</div>
</div>
<div>
<div>
<p class="xmsonormal">Hi,<o:p></o:p></p>
<p class="xmsonormal"> <o:p></o:p></p>
<p class="xmsonormal"><span lang="EN-US">Could you have a try with the default resampling “nearest”? And then also with subset having lower limit equal to upper. Now it fails, I tried already.</span><o:p></o:p></p>
<p class="xmsonormal"><span lang="EN-US"> </span><o:p></o:p></p>
<p class="xmsonormal"><a href="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,196396.75)&subset=y(450031.75,450031.75)"><span lang="EN-US">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,196396.75)&subset=y(450031.75,450031.75)</span></a><o:p></o:p></p>
<p class="xmsonormal"><span lang="EN-US"> </span><o:p></o:p></p>
<p class="xmsonormal"><span lang="EN-US">GetCoverage trim with equal low and high has at least worked in the past, see
</span><a href="http://osgeo-org.1560.x6.nabble.com/Why-WCS-2-0-1-slice-is-not-supported-td5343277.html"><span lang="EN-US">http://osgeo-org.1560.x6.nabble.com/Why-WCS-2-0-1-slice-is-not-supported-td5343277.html</span></a><o:p></o:p></p>
<p class="xmsonormal"><span lang="EN-US"> </span><o:p></o:p></p>
<p class="xmsonormal"><span lang="EN-US">It should be OK by the WCS standard because all intersected pixels should be included and in this case the intersection is a single pixel. But see also
<a href="https://osgeo-org.atlassian.net/browse/GEOS-9553">https://osgeo-org.atlassian.net/browse/GEOS-9553</a> and learn that Geoserver supports slicing but not trimming to one pixel.</span><o:p></o:p></p>
<p class="xmsonormal"><span lang="EN-US"> </span><o:p></o:p></p>
<p class="xmsonormal"><span lang="EN-US">-Jukka Rahkonen- </span><o:p></o:p></p>
<p class="xmsonormal"><o:p> </o:p></p>
</div>
</div>
</div>
</body>
</html>