<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=utf-8">
<meta name="Generator" content="Microsoft Word 15 (filtered medium)">
<style><!--
/* Font Definitions */
@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;}
@font-face
{font-family:"Lucida Console";
panose-1:2 11 6 9 4 5 4 2 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;}
pre
{mso-style-priority:99;
mso-style-link:"HTML-esimuotoiltu Char";
margin:0cm;
margin-bottom:.0001pt;
font-size:10.0pt;
font-family:"Courier New";}
span.gmail-n
{mso-style-name:gmail-n;}
span.gmail-s2
{mso-style-name:gmail-s2;}
span.gmail-o
{mso-style-name:gmail-o;}
span.HTML-esimuotoiltuChar
{mso-style-name:"HTML-esimuotoiltu Char";
mso-style-priority:99;
mso-style-link:HTML-esimuotoiltu;
font-family:Consolas;
mso-fareast-language:FI;}
span.Shkpostityyli24
{mso-style-type:personal-compose;
font-family:"Calibri",sans-serif;
color:windowtext;}
.MsoChpDefault
{mso-style-type:export-only;
font-family:"Calibri",sans-serif;
mso-fareast-language:EN-US;}
@page WordSection1
{size:612.0pt 792.0pt;
margin:70.85pt 2.0cm 70.85pt 2.0cm;}
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,<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">Sorry, this is something that I am not familiar with. The developer of the WCS plugin probably could tell if it is at all supported to use the XML file as input like you tried. I wonder
if Ari Jolma happens to follow the list. And what QGIS does is all mystery for me.<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>
<p class="MsoNormal"><span lang="EN-US" style="mso-fareast-language:EN-US">-Jukka Rahkonen-<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>
<p class="MsoNormal"><b>Lähettäjä:</b> 1520 gis <juliermeopensourcedeveloper@gmail.com>
<br>
<b>Lähetetty:</b> keskiviikko 11. marraskuuta 2020 16.55<br>
<b>Vastaanottaja:</b> Rahkonen Jukka (MML) <jukka.rahkonen@maanmittauslaitos.fi><br>
<b>Aihe:</b> Re: gdal-dev Digest, Vol 198, Issue 14<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<div>
<p class="MsoNormal">Dear Jukka Rahkonen,<o:p></o:p></p>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">Thank you very much for your patience, time, and detailed information given to my questions. I appreciated it all.<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">In your discussion in [1], I could get the a xml sample [2].<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">In [3], you will get the whole image available in geoserver. And in [4], You will get the Image [3] subset.<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">My goal is to load [3] e [4], mainly [4] subset in QGIS, by adding a XML raster. It is important to mention that I am not interested in accessing the raster data through a WCS connection, I am actually adding a xml as a raster (Ctrl + Shift +
R).<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">If you save the xml file in [2] or [5] and add to QGIS as raster data (Ctrl + Shift + R), you will get the image. <o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">Now, the xml you find in [5] was created by just changing the tag values according to your discussion in [1]. I could not find through gdal documentation how to generate this xml using gdal_translate. GDAL_WMTS driver stunningly does the
job. See <i><span style="font-size:11.5pt;font-family:"Arial",sans-serif;color:#888888"> </span></i><span class="gmail-n"><span style="font-size:9.0pt;font-family:"Lucida Console";color:#333333">gdal_translate</span></span><span style="font-size:9.0pt;font-family:"Lucida Console";color:#333333;background:#F8F8F8">
</span><span class="gmail-s2"><span style="font-size:9.0pt;font-family:"Lucida Console";color:#BB8844">"WMTS:<a href="http://maps.wien.gv.at/wmts/1.0.0/WMTSCapabilities.xml,layer=lb">http://maps.wien.gv.at/wmts/1.0.0/WMTSCapabilities.xml,layer=lb</a>"</span></span><span style="font-size:9.0pt;font-family:"Lucida Console";color:#333333;background:#F8F8F8">
</span><span class="gmail-n"><span style="font-size:9.0pt;font-family:"Lucida Console";color:#333333">wmts</span></span><span class="gmail-o"><b><span style="font-size:9.0pt;font-family:"Lucida Console";color:#333333">.</span></b></span><span class="gmail-n"><span style="font-size:9.0pt;font-family:"Lucida Console";color:#333333">xml</span></span><span style="font-size:9.0pt;font-family:"Lucida Console";color:#333333;background:#F8F8F8">
</span><span class="gmail-o"><b><span style="font-size:9.0pt;font-family:"Lucida Console";color:#333333">-</span></b></span><span class="gmail-n"><span style="font-size:9.0pt;font-family:"Lucida Console";color:#333333">of</span></span><span style="font-size:9.0pt;font-family:"Lucida Console";color:#333333;background:#F8F8F8">
</span><span class="gmail-n"><span style="font-size:9.0pt;font-family:"Lucida Console";color:#333333">WMTS.</span></span><o:p></o:p></p>
</div>
<p class="MsoNormal"><o:p> </o:p></p>
<div>
<p class="MsoNormal">My need is to create a xml for [4] just like the one available in [5]. I need to access the image in [4] subset using GDAL_WCS. <o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">The command line in [6] creates a tiff image, but I am looking for the xml that when loaded to QGIS ingest the WCS service and loads the image.<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">I wonder If there is a xml tag to the subset available in [4]. Or if there is a command line for gdal_translate that can generate XML in 3, and most importantly, generate a xml for the subset in [4], similar to gda GDAL_WMTS driver shown
above.<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">Sorry for the long e-mail. I have no words to thank you for your help and consideration of my messages.<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">Kinds Regards<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">Julierme<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal"> <o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">[1] <a href="https://geoserver-users.narkive.com/2FVEaNlm/define-geoserver-bounding-box-using-gdal-s-wcs-driver" target="_blank">https://geoserver-users.narkive.com/2FVEaNlm/define-geoserver-bounding-box-using-gdal-s-wcs-driver</a><o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">[2] XML Sample<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal"><WCS_GDAL><br>
<ServiceURL><a href="https://demo.geo-solutions.it/geoserver/wcs?version=2.0.1&">https://demo.geo-solutions.it/geoserver/wcs?version=2.0.1&</a>;</ServiceURL><br>
<CoverageName>nurc:mosaic</CoverageName><br>
<GetCoverageExtra>&BoundingBox=6.346,36.492,20.83,46.591</GetCoverageExtra><br>
</WCS_GDAL><o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:11.5pt;font-family:"Arial",sans-serif">[3]</span><a href="https://demo.geo-solutions.it/geoserver/wcs?Service=WCS">https://demo.geo-solutions.it/geoserver/wcs?Service=WCS</a><o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">&version=2.0.1<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">&bbox=-130.85168,20.7052,-62.0054,54.1141<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">&request=GetCoverage<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">&coverageId=nurc:Img_Sample<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">&format=image/png<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">[4] <a href="https://demo.geo-solutions.it/geoserver/wcs?service=WCS">
https://demo.geo-solutions.it/geoserver/wcs?service=WCS</a><o:p></o:p></p>
</div>
<p class="MsoNormal">&version=2.0.1<br>
&request=GetCoverage<br>
&SubsettingCRS=EPSG:4326<br>
&coverageId=nurc:Img_Sample<br>
&SUBSET=Long(-100,-80)&SUBSET=Lat(30,40)<br>
&format=image/png<o:p></o:p></p>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">[5] <WCS_GDAL><o:p></o:p></p>
</div>
<p class="MsoNormal"> <ServiceURL><a href="https://demo.geo-solutions.it/geoserver/wcs?version=2.0.1&">https://demo.geo-solutions.it/geoserver/wcs?version=2.0.1&</a>;</ServiceURL><br>
<CoverageName>nurc:Img_Sample</CoverageName><br>
<GetCoverageExtra>&BoundingBox=-130.85168,20.7052,-62.0054,54.1141</GetCoverageExtra><br>
</WCS_GDAL><o:p></o:p></p>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">[6] gdal_translate -oo CACHE=wcs_cache -oo CLEAR_CACHE -oo INTERLEAVE=PIXEL -projwin 377418 6683393.87938218 377717.879386966 6683094 -oo Subset="time(1985-01-01T00:00:00.000Z)" -outsize 60 0 "WCS:<a href="https://beta-karttakuva.maanmittauslaitos.fi/wcs/service/ows?version=2.0.1&coverage=ortokuva__ortokuva">https://beta-karttakuva.maanmittauslaitos.fi/wcs/service/ows?version=2.0.1&coverage=ortokuva__ortokuva</a>"
scaled.tiff<o:p></o:p></p>
</div>
</div>
<p class="MsoNormal"><o:p> </o:p></p>
<div>
<div>
<p class="MsoNormal">On Tue, Nov 10, 2020 at 1:23 PM Rahkonen Jukka (MML) <<a href="mailto:jukka.rahkonen@maanmittauslaitos.fi" target="_blank">jukka.rahkonen@maanmittauslaitos.fi</a>> wrote:<o:p></o:p></p>
</div>
<blockquote style="border:none;border-left:solid #CCCCCC 1.0pt;padding:0cm 0cm 0cm 6.0pt;margin-left:4.8pt;margin-right:0cm">
<div>
<div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto">Hi,<o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> <o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span lang="EN-US">Sending mail just to you because this is not a proper answer but rather some first-aid. Normally we keep the discussion on the list to other users can help, and
get help by the same.</span><o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span lang="EN-US"> </span><o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span lang="EN-US">Sorry, -spat is for ogr2ogr, in gdal_translate it is -projwin. Give the bbox coordinates as in the help:</span><o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span lang="EN-US"> [-projwin ulx uly lrx lry]</span><o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span lang="EN-US"> </span><o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span lang="EN-US">upper left x upper left y lower right x lower right y</span><o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span lang="EN-US"> </span><o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span lang="EN-US">Hope that you get it to work but I may look at this deeper when I have better time for that.</span><o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span lang="EN-US"> </span><o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span lang="EN-US">First step is to create the XML file. Then test is with gdalinfo, when it reports something that makes sense you can continue with gdal_translate and all the
standard options that it has.</span><o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span lang="EN-US"> </span><o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span lang="EN-US">WCS 2.0 is pretty easy to use also with plain hand written http requests, read the standard for learning. I would forget older versions.
</span><o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span lang="EN-US"> </span><o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span lang="EN-US">Test URL for you<br>
<a href="https://demo.geo-solutions.it/geoserver/wcs?SERVICE=WCS&REQUEST=GetCoverage&VERSION=2.0.1&CoverageId=nurc__mosaic&subset=Long(8.0,8.1)&subset=Lat(38,38.1)&format=tiff" target="_blank">https://demo.geo-solutions.it/geoserver/wcs?SERVICE=WCS&REQUEST=GetCoverage&VERSION=2.0.1&CoverageId=nurc__mosaic&subset=Long(8.0,8.1)&subset=Lat(38,38.1)&format=tiff</a></span><o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span lang="EN-US"> </span><o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span lang="EN-US">Unfortunately the server was down earlier today, I hope it works now.</span><o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span lang="EN-US"> </span><o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span lang="EN-US">-Jukka-</span><o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span lang="EN-US"> </span><o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><b>Lähettäjä:</b> 1520 gis <<a href="mailto:juliermeopensourcedeveloper@gmail.com" target="_blank">juliermeopensourcedeveloper@gmail.com</a>>
<br>
<b>Lähetetty:</b> tiistai 10. marraskuuta 2020 18.10<br>
<b>Vastaanottaja:</b> <a href="mailto:gdal-dev@lists.osgeo.org" target="_blank">gdal-dev@lists.osgeo.org</a>; Rahkonen Jukka (MML) <<a href="mailto:jukka.rahkonen@maanmittauslaitos.fi" target="_blank">jukka.rahkonen@maanmittauslaitos.fi</a>><br>
<b>Aihe:</b> Re: gdal-dev Digest, Vol 198, Issue 14<o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> <o:p></o:p></p>
<div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto">Dear Jukka Rahkonen,<o:p></o:p></p>
<div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> <o:p></o:p></p>
</div>
<div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto">Thank you very much for your reply.<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> <o:p></o:p></p>
</div>
<div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto">First of all, my apologies for my lack of knowledge about using GDAL_WCS driver. Let me give more details about the issues I am facing.<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> <o:p></o:p></p>
</div>
<div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto">Extract of <a href="http://172.21.14.45:8181/geoserver/ows/wcs?request=GetCapabilities" target="_blank">http://172.21.14.45:8181/geoserver/ows/wcs?request=GetCapabilities</a> is
shown in [1].<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> <o:p></o:p></p>
</div>
<div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto">Regarding gdal_translate GDAL_WCS driver Image subset extract, looking at <a href="https://gdal.org/drivers/raster/wcs.html#examples" target="_blank">https://gdal.org/drivers/raster/wcs.html#examples</a>,
I tried the following:<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> <o:p></o:p></p>
</div>
<div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto">Similar to GDAL_WMTS <<span style="font-size:9.0pt;font-family:"Lucida Console";color:#333333">gdal_translate<span style="background:#F8F8F8">
</span></span><span style="font-size:9.0pt;font-family:"Lucida Console";color:#BB8844">"WMTS:<a href="http://maps.wien.gv.at/wmts/1.0.0/WMTSCapabilities.xml,layer=lb" target="_blank">http://maps.wien.gv.at/wmts/1.0.0/WMTSCapabilities.xml,layer=lb</a>"</span><span style="font-size:9.0pt;font-family:"Lucida Console";color:#333333;background:#F8F8F8">
</span><span style="font-size:9.0pt;font-family:"Lucida Console";color:#333333">wmts<b>.</b>xml<span style="background:#F8F8F8">
</span><b>-</b>of<span style="background:#F8F8F8"> </span>WMTS>,</span><o:p></o:p></p>
</div>
<div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto">I tried running:<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto">(1)<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> gdal_translate "WCS:<a href="http://172.21.14.45:8181/geoserver/wcs?version=2.0.1&coverage=deter-amazonia__CBERS-4_AWFI_170_111_21052020" target="_blank">http://172.21.14.45:8181/geoserver/wcs?version=2.0.1&coverage=deter-amazonia__CBERS-4_AWFI_170_111_21052020</a>"
CBERS-4_AWFI_170_111_21052020.xml <o:p></o:p></p>
</div>
<div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> <o:p></o:p></p>
</div>
<div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><b>The response:</b> I did not get the xml as wmts.xml showed above<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> <o:p></o:p></p>
</div>
<div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto">Note: Including the -of option (WCS), the outcome fails<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> <o:p></o:p></p>
</div>
<div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto">(2)<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto">gdal_translate "WCS:<a href="http://172.21.14.45:8181/geoserver/ows/wcs?request=GetCapabilities&coverage=deter-amazonia__CBERS-4_AWFI_170_111_21052020" target="_blank">http://172.21.14.45:8181/geoserver/ows/wcs?request=GetCapabilities&coverage=deter-amazonia__CBERS-4_AWFI_170_111_21052020</a>"
\<br>
-spat -58 -11 -57.75 -10.7 \<br>
CBERS-4_AWFI_170_111_21052020.xml<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> <o:p></o:p></p>
</div>
<div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><b>The response:</b> There is no -spat option for gdal_translate<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> <o:p></o:p></p>
</div>
<div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><b>gdal_translate --help</b><o:p></o:p></p>
</div>
<div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> <o:p></o:p></p>
</div>
<div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto">Usage: gdal_translate [--help-general] [--long-usage]<br>
[-ot {Byte/Int16/UInt16/UInt32/Int32/Float32/Float64/<br>
CInt16/CInt32/CFloat32/CFloat64}] [-strict]<br>
[-of format] [-b band] [-mask band] [-expand {gray|rgb|rgba}]<br>
[-outsize xsize[%]|0 ysize[%]|0] [-tr xres yres]<br>
[-r {nearest,bilinear,cubic,cubicspline,lanczos,average,mode}]<br>
[-unscale] [-scale[_bn] [src_min src_max [dst_min dst_max]]]* [-exponent[_bn] exp_val]*<br>
[-srcwin xoff yoff xsize ysize] [-epo] [-eco]<br>
[-projwin ulx uly lrx lry] [-projwin_srs srs_def]<br>
[-a_srs srs_def] [-a_ullr ulx uly lrx lry] [-a_nodata value]<br>
[-a_scale value] [-a_offset value]<br>
[-nogcp] [-gcp pixel line easting northing [elevation]]*<br>
|-colorinterp{_bn} {red|green|blue|alpha|gray|undefined}]<br>
|-colorinterp {red|green|blue|alpha|gray|undefined},...]<br>
[-mo "META-TAG=VALUE"]* [-q] [-sds]<br>
[-co "NAME=VALUE"]* [-stats] [-norat]<br>
[-oo NAME=VALUE]*<br>
src_dataset dst_dataset<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> <o:p></o:p></p>
</div>
<div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> <o:p></o:p></p>
</div>
<div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> <o:p></o:p></p>
</div>
<div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto">[1] WCS GetCapabilities request <o:p></o:p></p>
</div>
<div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> <o:p></o:p></p>
</div>
<div>
<div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:10.0pt;font-family:"Courier New";color:black"><wcs:CoverageSummary></span><o:p></o:p></p>
</div>
<div style="margin-left:12.0pt">
<div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:10.0pt;font-family:"Courier New";color:black"><ows:Title>CBERS-4_AWFI_170_111_21052020</ows:Title></span><o:p></o:p></p>
</div>
<div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:10.0pt;font-family:"Courier New";color:black"><ows:Abstract>Generated from GeoTIFF</ows:Abstract></span><o:p></o:p></p>
</div>
<div id="gmail-m_7107943894366576226gmail-m_8371735024734129097gmail-folder28">
<div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:10.0pt;font-family:"Courier New";color:black"><ows:Keywords></span><o:p></o:p></p>
</div>
<div style="margin-left:12.0pt">
<div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:10.0pt;font-family:"Courier New";color:black"><ows:Keyword>deter-amazonia-CBERS-4_AWFI_170_111_21052020</ows:Keyword></span><o:p></o:p></p>
</div>
<div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:10.0pt;font-family:"Courier New";color:black"><ows:Keyword>WCS</ows:Keyword></span><o:p></o:p></p>
</div>
<div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:10.0pt;font-family:"Courier New";color:black"><ows:Keyword>GeoTIFF</ows:Keyword></span><o:p></o:p></p>
</div>
</div>
<div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:10.0pt;font-family:"Courier New";color:black"></ows:Keywords></span><o:p></o:p></p>
</div>
</div>
<div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:10.0pt;font-family:"Courier New";color:black"><wcs:CoverageId>deter-amazonia__CBERS-4_AWFI_170_111_21052020</wcs:CoverageId></span><o:p></o:p></p>
</div>
<div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:10.0pt;font-family:"Courier New";color:black"><wcs:CoverageSubtype>RectifiedGridCoverage</wcs:CoverageSubtype></span><o:p></o:p></p>
</div>
<div id="gmail-m_7107943894366576226gmail-m_8371735024734129097gmail-folder29">
<div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:10.0pt;font-family:"Courier New";color:black"><ows:WGS84BoundingBox></span><o:p></o:p></p>
</div>
<div style="margin-left:12.0pt">
<div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:10.0pt;font-family:"Courier New";color:black"><ows:LowerCorner>-63.280290728152664 -13.794296142085976</ows:LowerCorner></span><o:p></o:p></p>
</div>
<div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:10.0pt;font-family:"Courier New";color:black"><ows:UpperCorner>-53.624004749255874 -5.636644661059637</ows:UpperCorner></span><o:p></o:p></p>
</div>
</div>
<div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:10.0pt;font-family:"Courier New";color:black"></ows:WGS84BoundingBox></span><o:p></o:p></p>
</div>
</div>
<div id="gmail-m_7107943894366576226gmail-m_8371735024734129097gmail-folder30">
<div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:10.0pt;font-family:"Courier New";color:black"><ows:BoundingBox crs="<a href="http://www.opengis.net/def/crs/EPSG/0/EPSG:4326" target="_blank">http://www.opengis.net/def/crs/EPSG/0/EPSG:4326</a>"></span><o:p></o:p></p>
</div>
<div style="margin-left:12.0pt">
<div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:10.0pt;font-family:"Courier New";color:black"><ows:LowerCorner>-63.280290728152664 -13.794296142085976</ows:LowerCorner></span><o:p></o:p></p>
</div>
<div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:10.0pt;font-family:"Courier New";color:black"><ows:UpperCorner>-53.624004749255874 -5.636644661059637</ows:UpperCorner></span><o:p></o:p></p>
</div>
</div>
<div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:10.0pt;font-family:"Courier New";color:black"></ows:BoundingBox></span><o:p></o:p></p>
</div>
</div>
</div>
<div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:10.0pt;font-family:"Courier New";color:black"></wcs:CoverageSummary></span><o:p></o:p></p>
</div>
</div>
<div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> <o:p></o:p></p>
</div>
<div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> <o:p></o:p></p>
</div>
<div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto">Any comment will be very appreciated.<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> <o:p></o:p></p>
</div>
<div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto">GDAL
<a href="https://gdal.org/drivers/raster/wmts.html" target="_blank">website</a> provides examples on how to run GDAL_WMTS driver in order to get a wmts xml file. I have looked in <a href="https://gdal.org/drivers/raster/wcs.html" target="_blank">https://gdal.org/drivers/raster/wcs.html</a>
and could not find an example pointing to -spat option. I got an -spat option <o:p></o:p></p>
<div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto">to<a href="https://gdal.org/programs/gdal_grid.html?highlight=spat#cmdoption-gdal_grid-spat" target="_blank"> gdal_grid</a> though. <o:p></o:p></p>
</div>
<div>
<pre style="background:#F8F8F8;box-sizing:border-box;overflow:auto"><span style="font-size:9.0pt;font-family:"Lucida Console";color:#333333">gdal_translate </span><span style="font-size:9.0pt;font-family:"Lucida Console";color:#BB8844">"WMTS:<a href="http://maps.wien.gv.at/wmts/1.0.0/WMTSCapabilities.xml,layer=lb" target="_blank">http://maps.wien.gv.at/wmts/1.0.0/WMTSCapabilities.xml,layer=lb</a>"</span><span style="font-size:9.0pt;font-family:"Lucida Console";color:#333333"> wmts<b>.</b>xml <b>-</b>of WMTS</span><o:p></o:p></pre>
<pre style="background:#F8F8F8;box-sizing:border-box;overflow:auto"><span style="font-size:9.0pt;font-family:"Lucida Console";color:#333333">gdal_translate "WCS:<a href="http://www.dpi.inpe.br/fipcerrado-geoserver/ows/wcs?request=GetCapabilities" target="_blank">http://www.dpi.inpe.br/fipcerrado-geoserver/ows/wcs?request=GetCapabilities</a>" \</span><o:p></o:p></pre>
<pre style="margin-bottom:12.0pt;background:#F8F8F8;box-sizing:border-box;overflow:auto"><span style="font-size:9.0pt;font-family:"Lucida Console";color:#333333">-oo CoverageName=deter-cerrado:CBERS-4_AWFI_161_099_13072018 CBERS-4_AWFI_161_099_13072018.xml -of WCS</span><o:p></o:p></pre>
<pre style="background:#F8F8F8;box-sizing:border-box;overflow:auto"><span style="font-size:9.0pt;font-family:"Lucida Console";color:#333333"> </span><o:p></o:p></pre>
</div>
</div>
</div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> <o:p></o:p></p>
<div>
<div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto">On Mon, Nov 9, 2020 at 5:00 PM <<a href="mailto:gdal-dev-request@lists.osgeo.org" target="_blank">gdal-dev-request@lists.osgeo.org</a>> wrote:<o:p></o:p></p>
</div>
<blockquote style="border:none;border-left:solid #CCCCCC 1.0pt;padding:0cm 0cm 0cm 6.0pt;margin-left:4.8pt;margin-top:5.0pt;margin-right:0cm;margin-bottom:5.0pt">
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto">Send gdal-dev mailing list submissions to<br>
<a href="mailto:gdal-dev@lists.osgeo.org" target="_blank">gdal-dev@lists.osgeo.org</a><br>
<br>
To subscribe or unsubscribe via the World Wide Web, visit<br>
<a href="https://lists.osgeo.org/mailman/listinfo/gdal-dev" target="_blank">
https://lists.osgeo.org/mailman/listinfo/gdal-dev</a><br>
or, via email, send a message with subject or body 'help' to<br>
<a href="mailto:gdal-dev-request@lists.osgeo.org" target="_blank">gdal-dev-request@lists.osgeo.org</a><br>
<br>
You can reach the person managing the list at<br>
<a href="mailto:gdal-dev-owner@lists.osgeo.org" target="_blank">gdal-dev-owner@lists.osgeo.org</a><br>
<br>
When replying, please edit your Subject line so it is more specific<br>
than "Re: Contents of gdal-dev digest..."<br>
<br>
<br>
Today's Topics:<br>
<br>
1. Re: Selecting just the geometry column with OGR SQL (Even Rouault)<br>
2. Re: Selecting just the geometry column with OGR SQL<br>
(Rahkonen Jukka (MML))<br>
3. Re: WCS_GDAL Driver Image SubSet Request (jratike80)<br>
<br>
<br>
----------------------------------------------------------------------<br>
<br>
Message: 1<br>
Date: Mon, 09 Nov 2020 16:09:45 +0100<br>
From: Even Rouault <<a href="mailto:even.rouault@spatialys.com" target="_blank">even.rouault@spatialys.com</a>><br>
To: <a href="mailto:gdal-dev@lists.osgeo.org" target="_blank">gdal-dev@lists.osgeo.org</a><br>
Cc: "Rahkonen Jukka (MML)" <<a href="mailto:jukka.rahkonen@maanmittauslaitos.fi" target="_blank">jukka.rahkonen@maanmittauslaitos.fi</a>><br>
Subject: Re: [gdal-dev] Selecting just the geometry column with OGR<br>
SQL<br>
Message-ID: <<a href="mailto:1659214.2tRo4lr4Oo@even-i700" target="_blank">1659214.2tRo4lr4Oo@even-i700</a>><br>
Content-Type: text/plain; charset="us-ascii"<br>
<br>
Jukka,<br>
<br>
> I believe that OGR SQL dialect adds the geometry column into SQL selection<br>
> by default. Am I right with this? I volunteer to edit the documentation<br>
> <a href="https://gdal.org/user/ogr_sql_dialect.html" target="_blank">https://gdal.org/user/ogr_sql_dialect.html</a> if this is the case. But what I<br>
> can't understand is another side of the story, how to select just the<br>
> geometry field with OGR SQL dialect from some data sources.<br>
> <br>
> This command returns both time column and geometry from GML<br>
> ogrinfo -sql "select time from timetest" timetest.gml<br>
> I can select just the geometry from geometry column named as<br>
> "geometryProperty" by ogrinfo -sql "select geometryProperty from timetest"<br>
> timetest.gml<br>
> <br>
> The behavior is the same with SpatiaLite when I use "-dialect OGRSQL".<br>
> Geometry is selected automatically but it can be selected also by name. But<br>
> I cannot discover any way to select just the geometry field from shapefile.<br>
> Is it because geometry column in shapefile does not have any name as this<br>
> Python test seems to prove?<br>
> >>> from osgeo import ogr<br>
> >>> shp_ds = ogr.Open('timetest.shp')<br>
> >>> shp_lyr = shp_ds.GetLayer(0)<br>
> >>> shp_lyr.GetGeometryColumn()<br>
> <br>
> ''<br>
> <br>
> OpenJUMP JML driver does not report a name for geometry column either. Is it<br>
> rather an exception that drivers give names for the geometry columns? I<br>
> know I can select only the geometry column with the SQLite dialect, but I<br>
> am curious.<br>
<br>
Indeed, in OGR SQL, if the geometry column has a non-empty name (typically for <br>
database drivers), you can use it to select it. When the geometry column name <br>
is empty, such as in the shapefile or OpenJUMP driver, you can use <br>
_ogr_geometry_ (with leading and trainling underscore) to select it. This is <br>
likely undocumented.<br>
<br>
As far as I can see, this was introduced per<br>
<a href="https://github.com/OSGeo/gdal/commit/1519dcb00bd7802d2f701114f0e1cb9b524a547f" target="_blank">https://github.com/OSGeo/gdal/commit/1519dcb00bd7802d2f701114f0e1cb9b524a547f</a><br>
to fix an issue with comparison with empty string literals (this was before <br>
OGR enforced the difference between single-quoted string literals and double-<br>
quoted identifiers)<br>
<br>
Even<br>
<br>
-- <br>
Spatialys - Geospatial professional services<br>
<a href="http://www.spatialys.com" target="_blank">http://www.spatialys.com</a><br>
<br>
<br>
------------------------------<br>
<br>
Message: 2<br>
Date: Mon, 9 Nov 2020 15:35:11 +0000<br>
From: "Rahkonen Jukka (MML)" <<a href="mailto:jukka.rahkonen@maanmittauslaitos.fi" target="_blank">jukka.rahkonen@maanmittauslaitos.fi</a>><br>
To: Even Rouault <<a href="mailto:even.rouault@spatialys.com" target="_blank">even.rouault@spatialys.com</a>>,<br>
"<a href="mailto:gdal-dev@lists.osgeo.org" target="_blank">gdal-dev@lists.osgeo.org</a>" <<a href="mailto:gdal-dev@lists.osgeo.org" target="_blank">gdal-dev@lists.osgeo.org</a>><br>
Subject: Re: [gdal-dev] Selecting just the geometry column with OGR<br>
SQL<br>
Message-ID: <<a href="mailto:d59795ddfd994dee9f5eca836361324a@maanmittauslaitos.fi" target="_blank">d59795ddfd994dee9f5eca836361324a@maanmittauslaitos.fi</a>><br>
Content-Type: text/plain; charset="us-ascii"<br>
<br>
Even Rouault wrote: <br>
maanantai 9. marraskuuta 2020 17.10<br>
<br>
> Jukka,<br>
<br>
>> I believe that OGR SQL dialect adds the geometry column into SQL <br>
>> selection by default. Am I right with this? I volunteer to edit the <br>
>> documentation <a href="https://gdal.org/user/ogr_sql_dialect.html" target="_blank">
https://gdal.org/user/ogr_sql_dialect.html</a> if this is <br>
>> the case. But what I can't understand is another side of the story, <br>
>> how to select just the geometry field with OGR SQL dialect from some data sources.<br>
>> <br>
>> This command returns both time column and geometry from GML ogrinfo <br>
>> -sql "select time from timetest" timetest.gml I can select just the <br>
>> geometry from geometry column named as "geometryProperty" by ogrinfo <br>
>> -sql "select geometryProperty from timetest"<br>
>> timetest.gml<br>
>> <br>
>> The behavior is the same with SpatiaLite when I use "-dialect OGRSQL".<br>
>> Geometry is selected automatically but it can be selected also by <br>
>> name. But I cannot discover any way to select just the geometry field from shapefile.<br>
>> Is it because geometry column in shapefile does not have any name as <br>
>> this Python test seems to prove?<br>
>> >>> from osgeo import ogr<br>
>> >>> shp_ds = ogr.Open('timetest.shp')<br>
>> >>> shp_lyr = shp_ds.GetLayer(0)<br>
>> >>> shp_lyr.GetGeometryColumn()<br>
>> <br>
>> ''<br>
>> <br>
>> OpenJUMP JML driver does not report a name for geometry column either. <br>
>> Is it rather an exception that drivers give names for the geometry <br>
>> columns? I know I can select only the geometry column with the SQLite <br>
>> dialect, but I am curious.<br>
<br>
> Indeed, in OGR SQL, if the geometry column has a nonempty name (typically for <br>
> database drivers), you can use it to select it. When the geometry column name <br>
> is empty, such as in the shapefile or OpenJUMP driver, you can use <br>
> _ogr_geometry_ (with leading and trainling underscore) to select it. This is likely undocumented.<br>
<br>
> As far as I can see, this was introduced per <br>
> <a href="https://github.com/OSGeo/gdal/commit/1519dcb00bd7802d2f701114f0e1cb9b524a547f" target="_blank">
https://github.com/OSGeo/gdal/commit/1519dcb00bd7802d2f701114f0e1cb9b524a547f</a><br>
> to fix an issue with comparison with empty string literals (this was before OGR enforced the<br>
> difference between single-quoted string literals and double- quoted identifiers)<br>
<br>
I could not manage to guess the right syntax for Windows right ahead but solved the quiz easily
<br>
with your double-quote hint<br>
<br>
First trial:<br>
ogrinfo -sql "select _ogr_geometry_ from timetest" timetest.shp<br>
INFO: Open of `timetest.shp'<br>
using driver `ESRI Shapefile' successful.<br>
ERROR 1: SQL Expression Parsing Error: syntax error, unexpected $undefined. Occurred around :<br>
select _ogr_geometry_ from timetest<br>
<br>
Second trial:<br>
ogrinfo -sql "select \"_ogr_geometry_\" from timetest" timetest.shp <br>
<br>
-Jukka-<br>
<br>
> Even<br>
<br>
--<br>
Spatialys - Geospatial professional services <a href="http://www.spatialys.com" target="_blank">
http://www.spatialys.com</a><br>
<br>
<br>
------------------------------<br>
<br>
Message: 3<br>
Date: Mon, 9 Nov 2020 08:47:05 -0700 (MST)<br>
From: jratike80 <<a href="mailto:jukka.rahkonen@maanmittauslaitos.fi" target="_blank">jukka.rahkonen@maanmittauslaitos.fi</a>><br>
To: <a href="mailto:gdal-dev@lists.osgeo.org" target="_blank">gdal-dev@lists.osgeo.org</a><br>
Subject: Re: [gdal-dev] WCS_GDAL Driver Image SubSet Request<br>
Message-ID: <<a href="mailto:1604936825138-0.post@n6.nabble.com" target="_blank">1604936825138-0.post@n6.nabble.com</a>><br>
Content-Type: text/plain; charset=us-ascii<br>
<br>
Hi,<br>
<br>
Do you want to make is somehow impossible to request data beyound the<br>
BoundingBox that you defined in the XML file? If you just want to get data<br>
from the rectangle that you or your users define, gdal_translate with the<br>
regular -spat option should work out-of-the-box.<br>
<br>
-Jukka Rahkonen- <br>
<br>
<br>
juliermeopensourcedeveloper wrote<br>
> Hi all,<br>
> <br>
> Is it possible to define a bounding box in order to request a coverage<br>
> subset using the GDAL WCS Driver?<br>
> I have tried the code in [1]. I am trying to get information only for the<br>
> bbox_interest=( -58,-11,-57.75,-10.75) as shown below, however the subset<br>
> image or bbox_interest request has not been delivered. The West Longitude<br>
> reaches -50.17 and The Lower Latitude reaches -18.92 which are limits out<br>
> of the bbox_interest.<br>
> The image extent is:<br>
> <GetCoverageExtra><br>
> &BoundingBox=-63.280290728152664,-13.794296142085976,-53.624004749255874,-5.636644661059637<br>
> </GetCoverageExtra><br>
> Any comment on this issue will be greatly appreciated.<br>
> Thank you very much for your time in advance.<br>
> *CODE [1]:*<br>
> <WCS_GDAL><br>
> <ServiceURL><br>
> <a href="http://172.21.14.45:6060/geoserver/wcs?version=1.0.0&" target="_blank">
http://172.21.14.45:6060/geoserver/wcs?version=1.0.0&</a>;<br>
> </<br>
> ServiceURL><br>
> <CoverageName><br>
> cbers_deter_amazonia:CBERS-4_AWFI_170_111_21052020<br>
> </<br>
> CoverageName><br>
> <GetCoverageExtra><br>
> &BoundingBox=-58,-11,-57.75,-10.75<br>
> </GetCoverageExtra><br>
> <CoverageOffering><br>
> <description><br>
> Generated from GeoTIFF<br>
> </description><br>
> <name><br>
> cbers_deter_amazonia:CBERS-4_AWFI_170_111_21052020<br>
> </name><br>
> <label><br>
> CBERS-4_AWFI_170_111_21052020<br>
> </label><br>
> <lonLatEnvelope srsName="urn:ogc:def:crs:OGC:1.3:CRS84"><br>
> <pos><br>
> -58 -11<br>
> </pos><br>
> <pos><br>
> -57.75 -10.75<br>
> </pos><br>
> </lonLatEnvelope><br>
> <keywords><br>
> <keyword><br>
> WCS<br>
> </keyword><br>
> <keyword><br>
> GeoTIFF<br>
> </keyword><br>
> </keywords><br>
> <domainSet><br>
> <spatialDomain><br>
> <Envelope srsName="EPSG:4326"><br>
> <pos><br>
> -58 -11<br>
> </pos><br>
> <pos><br>
> -57.75 -10.75<br>
> </pos><br>
> </Envelope><br>
> <RectifiedGrid dimension="2" srsName="EPSG:4326"><br>
> <limits><br>
> <GridEnvelope><br>
> <low><br>
> 0 0<br>
> </low><br>
> <high><br>
> 0 0<br>
> </high><br>
> </GridEnvelope><br>
> </limits><br>
> <axisName><br>
> x<br>
> </axisName><br>
> <axisName><br>
> y<br>
> </axisName><br>
> <origin><br>
> <pos><br>
> -58 -10.75<br>
> </pos><br>
> </origin><br>
> <offsetVector><br>
> 0.005364603321609328 0.0<br>
> </offsetVector><br>
> <offsetVector><br>
> 0.0 -0.004532028600570188<br>
> </offsetVector><br>
> </RectifiedGrid><br>
> </spatialDomain><br>
> </domainSet><br>
> <rangeSet><br>
> <RangeSet><br>
> <name><br>
> CBERS-4_AWFI_170_111_21052020<br>
> </name><br>
> <label><br>
> CBERS-4_AWFI_170_111_21052020<br>
> </label><br>
> <axisDescription><br>
> <AxisDescription><br>
> <name><br>
> Band<br>
> </name><br>
> <label><br>
> Band<br>
> </label><br>
> <values><br>
> <interval><br>
> <min><br>
> 1<br>
> </min><br>
> <max><br>
> 4<br>
> </max><br>
> </interval><br>
> </values><br>
> </AxisDescription><br>
> </axisDescription><br>
> </RangeSet><br>
> </rangeSet><br>
> <supportedCRSs><br>
> <requestResponseCRSs><br>
> EPSG:4326<br>
> </requestResponseCRSs><br>
> </supportedCRSs><br>
> <supportedFormats nativeFormat="GeoTIFF"><br>
> <formats><br>
> GeoTIFF<br>
> </formats><br>
> <formats><br>
> GIF<br>
> </formats><br>
> <formats><br>
> JPEG<br>
> </formats><br>
> <formats><br>
> PNG<br>
> </formats><br>
> <formats><br>
> TIFF<br>
> </formats><br>
> </supportedFormats><br>
> <supportedInterpolations default="nearest neighbor"><br>
> <interpolationMethod><br>
> nearest neighbor<br>
> </interpolationMethod><br>
> <interpolationMethod><br>
> bilinear<br>
> </interpolationMethod><br>
> <interpolationMethod><br>
> bicubic<br>
> </interpolationMethod><br>
> </supportedInterpolations><br>
> </CoverageOffering><br>
> <PreferredFormat><br>
> GeoTIFF<br>
> </PreferredFormat><br>
> <BandCount><br>
> 4<br>
> </BandCount><br>
> <BandType><br>
> Byte<br>
> </BandType><br>
> </WCS_GDAL><br>
> _______________________________________________<br>
> gdal-dev mailing list<br>
<br>
> <a href="mailto:gdal-dev@.osgeo" target="_blank">gdal-dev@.osgeo</a><br>
<br>
> <a href="https://lists.osgeo.org/mailman/listinfo/gdal-dev" target="_blank">https://lists.osgeo.org/mailman/listinfo/gdal-dev</a><br>
<br>
<br>
<br>
<br>
<br>
--<br>
Sent from: <a href="http://osgeo-org.1560.x6.nabble.com/GDAL-Dev-f3742093.html" target="_blank">
http://osgeo-org.1560.x6.nabble.com/GDAL-Dev-f3742093.html</a><br>
<br>
<br>
------------------------------<br>
<br>
Subject: Digest Footer<br>
<br>
_______________________________________________<br>
gdal-dev mailing list<br>
<a href="mailto:gdal-dev@lists.osgeo.org" target="_blank">gdal-dev@lists.osgeo.org</a><br>
<a href="https://lists.osgeo.org/mailman/listinfo/gdal-dev" target="_blank">https://lists.osgeo.org/mailman/listinfo/gdal-dev</a><br>
<br>
------------------------------<br>
<br>
End of gdal-dev Digest, Vol 198, Issue 14<br>
*****************************************<o:p></o:p></p>
</blockquote>
</div>
</div>
</div>
</blockquote>
</div>
</div>
</body>
</html>