<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=UTF-8">
</head>
<body>
<div class="moz-cite-prefix">It's already some time I worked with
the WCS driver and then most of the time went into wondering the
small differences between servers.</div>
<div class="moz-cite-prefix"><br>
</div>
<div class="moz-cite-prefix">With
<a class="moz-txt-link-freetext" href="https://demo.geo-solutions.it/geoserver/wcs">https://demo.geo-solutions.it/geoserver/wcs</a> I can do this
workflow:</div>
<div class="moz-cite-prefix"><br>
</div>
<div class="moz-cite-prefix">1) gdalinfo
"WCS:<a class="moz-txt-link-freetext" href="https://demo.geo-solutions.it/geoserver/wcs">https://demo.geo-solutions.it/geoserver/wcs</a>"</div>
<div class="moz-cite-prefix"><br>
</div>
<div class="moz-cite-prefix">..takes some time as it makes calls and
then it stores into $HOME/.gdal/wcs_cache/</div>
<div class="moz-cite-prefix">(note the version added to the URL)<br>
</div>
<div class="moz-cite-prefix"><br>
</div>
<div class="moz-cite-prefix">- an entry into db:
KfSWH=<a class="moz-txt-link-freetext" href="https://demo.geo-solutions.it/geoserver/wcs?version=2.0.1">https://demo.geo-solutions.it/geoserver/wcs?version=2.0.1</a></div>
<div class="moz-cite-prefix">- two files: KfSWH.xml and
KfSWH.aux.xml, the first is the capabilities XML, the 2nd is GDAL
aux XML</div>
<div class="moz-cite-prefix"><br>
</div>
<div class="moz-cite-prefix">and prints gdalinfo output, which
contains for example the subdatasets.</div>
<div class="moz-cite-prefix">You'll get the same output, but faster,
if you issue the above command again</div>
<div class="moz-cite-prefix"><br>
</div>
<div class="moz-cite-prefix">From the subdataset list I can now pick
one coverage.<br>
</div>
<div class="moz-cite-prefix"><br>
</div>
<div class="moz-cite-prefix">2) gdalinfo
"WCS:<a class="moz-txt-link-freetext" href="https://demo.geo-solutions.it/geoserver/wcs?version=2.0.1&coverage=nurc__Img_Sample">https://demo.geo-solutions.it/geoserver/wcs?version=2.0.1&coverage=nurc__Img_Sample</a>"</div>
<div class="moz-cite-prefix"><br>
</div>
<div class="moz-cite-prefix">.. takes again some time as it makes
calls and then it stores new data into the cache</div>
<div class="moz-cite-prefix"><br>
</div>
<div class="moz-cite-prefix">-
cJLrY=<a class="moz-txt-link-freetext" href="https://demo.geo-solutions.it/geoserver/wcs?version=2.0.1&coverage=nurc__Img_Sample">https://demo.geo-solutions.it/geoserver/wcs?version=2.0.1&coverage=nurc__Img_Sample</a></div>
<div class="moz-cite-prefix">- three files: cJLrY.xml,
cJLrY.xml.aux.xml, cJLrY.DC.xml, the first is the WCS_GDAL XML</div>
<div class="moz-cite-prefix"><br>
</div>
<div class="moz-cite-prefix">and prints</div>
<div class="moz-cite-prefix"><br>
</div>
<div class="moz-cite-prefix">ERROR 1: HTTP error code : 404<br>
ERROR 1: InvalidSubsetting: Empty intersection after subsetting</div>
<div class="moz-cite-prefix"><br>
</div>
<div class="moz-cite-prefix">so it probably needs some special
switches, which are described in the documentation.</div>
<div class="moz-cite-prefix"><br>
</div>
<div class="moz-cite-prefix">Best,</div>
<div class="moz-cite-prefix"><br>
</div>
<div class="moz-cite-prefix">Ari<br>
</div>
<div class="moz-cite-prefix"><br>
</div>
<div class="moz-cite-prefix">Rahkonen Jukka (MML) kirjoitti
11.11.2020 klo 23.55:<br>
</div>
<blockquote type="cite"
cite="mid:77b40009b5f046418ab7a464b50f8504@maanmittauslaitos.fi">
<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]-->
<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 style="mso-fareast-language:EN-US"
lang="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 style="mso-fareast-language:EN-US"
lang="EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><span style="mso-fareast-language:EN-US"
lang="EN-US">-Jukka Rahkonen-<o:p></o:p></span></p>
<p class="MsoNormal"><span style="mso-fareast-language:EN-US"
lang="EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><b>Lähettäjä:</b> 1520 gis
<a class="moz-txt-link-rfc2396E" href="mailto:juliermeopensourcedeveloper@gmail.com"><juliermeopensourcedeveloper@gmail.com></a>
<br>
<b>Lähetetty:</b> keskiviikko 11. marraskuuta 2020 16.55<br>
<b>Vastaanottaja:</b> Rahkonen Jukka (MML)
<a class="moz-txt-link-rfc2396E" href="mailto:jukka.rahkonen@maanmittauslaitos.fi"><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"><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"
moz-do-not-send="true">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" moz-do-not-send="true">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&"
moz-do-not-send="true">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"
moz-do-not-send="true">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"
moz-do-not-send="true">
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&"
moz-do-not-send="true">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"
moz-do-not-send="true">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" moz-do-not-send="true">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" moz-do-not-send="true">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" moz-do-not-send="true">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" moz-do-not-send="true">gdal-dev@lists.osgeo.org</a>;
Rahkonen Jukka (MML) <<a
href="mailto:jukka.rahkonen@maanmittauslaitos.fi"
target="_blank" moz-do-not-send="true">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" moz-do-not-send="true">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" moz-do-not-send="true">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" moz-do-not-send="true">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" moz-do-not-send="true">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" moz-do-not-send="true">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"
moz-do-not-send="true">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" moz-do-not-send="true">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" moz-do-not-send="true">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" moz-do-not-send="true">
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" moz-do-not-send="true">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" moz-do-not-send="true">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" moz-do-not-send="true">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" moz-do-not-send="true">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" moz-do-not-send="true">
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" moz-do-not-send="true">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" moz-do-not-send="true">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" moz-do-not-send="true">even.rouault@spatialys.com</a>><br>
To: <a href="mailto:gdal-dev@lists.osgeo.org"
target="_blank" moz-do-not-send="true">gdal-dev@lists.osgeo.org</a><br>
Cc: "Rahkonen Jukka (MML)" <<a
href="mailto:jukka.rahkonen@maanmittauslaitos.fi"
target="_blank" moz-do-not-send="true">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" moz-do-not-send="true">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" moz-do-not-send="true">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" moz-do-not-send="true">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"
moz-do-not-send="true">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" moz-do-not-send="true">jukka.rahkonen@maanmittauslaitos.fi</a>><br>
To: Even Rouault <<a
href="mailto:even.rouault@spatialys.com"
target="_blank" moz-do-not-send="true">even.rouault@spatialys.com</a>>,<br>
"<a href="mailto:gdal-dev@lists.osgeo.org"
target="_blank" moz-do-not-send="true">gdal-dev@lists.osgeo.org</a>"
<<a href="mailto:gdal-dev@lists.osgeo.org"
target="_blank" moz-do-not-send="true">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" moz-do-not-send="true">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" moz-do-not-send="true">
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" moz-do-not-send="true">
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"
moz-do-not-send="true">
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" moz-do-not-send="true">jukka.rahkonen@maanmittauslaitos.fi</a>><br>
To: <a href="mailto:gdal-dev@lists.osgeo.org"
target="_blank" moz-do-not-send="true">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" moz-do-not-send="true">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" moz-do-not-send="true">
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" moz-do-not-send="true">gdal-dev@.osgeo</a><br>
<br>
> <a
href="https://lists.osgeo.org/mailman/listinfo/gdal-dev"
target="_blank" moz-do-not-send="true">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" moz-do-not-send="true">
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" moz-do-not-send="true">gdal-dev@lists.osgeo.org</a><br>
<a
href="https://lists.osgeo.org/mailman/listinfo/gdal-dev"
target="_blank" moz-do-not-send="true">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>
<br>
<fieldset class="mimeAttachmentHeader"></fieldset>
<pre class="moz-quote-pre" wrap="">_______________________________________________
gdal-dev mailing list
<a class="moz-txt-link-abbreviated" href="mailto:gdal-dev@lists.osgeo.org">gdal-dev@lists.osgeo.org</a>
<a class="moz-txt-link-freetext" href="https://lists.osgeo.org/mailman/listinfo/gdal-dev">https://lists.osgeo.org/mailman/listinfo/gdal-dev</a></pre>
</blockquote>
</body>
</html>