<div dir="ltr">Dear Jukka Rahkonen,<div><br></div><div>Thank you very much for your reply.</div><div><br></div><div>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.</div><div><br></div><div>Extract of <a href="http://172.21.14.45:8181/geoserver/ows/wcs?request=GetCapabilities">http://172.21.14.45:8181/geoserver/ows/wcs?request=GetCapabilities</a> is shown in [1].</div><div><br></div><div>Regarding gdal_translate GDAL_WCS driver Image subset extract, looking at <a href="https://gdal.org/drivers/raster/wcs.html#examples">https://gdal.org/drivers/raster/wcs.html#examples</a>, I tried the following:<br></div><div><br></div><div>Similar to GDAL_WMTS <<span class="gmail-n" style="color:rgb(51,51,51);font-family:"Source Code Pro","Andale Mono WT","Andale Mono","Lucida Console","Lucida Sans Typewriter","DejaVu Sans Mono","Bitstream Vera Sans Mono","Liberation Mono","Nimbus Mono L",Monaco,"Courier New",Courier,monospace;font-size:12px;box-sizing:border-box">gdal_translate</span><span style="background-color:rgb(248,248,248);color:rgb(51,51,51);font-family:"Source Code Pro","Andale Mono WT","Andale Mono","Lucida Console","Lucida Sans Typewriter","DejaVu Sans Mono","Bitstream Vera Sans Mono","Liberation Mono","Nimbus Mono L",Monaco,"Courier New",Courier,monospace;font-size:12px"> </span><span class="gmail-s2" style="font-family:"Source Code Pro","Andale Mono WT","Andale Mono","Lucida Console","Lucida Sans Typewriter","DejaVu Sans Mono","Bitstream Vera Sans Mono","Liberation Mono","Nimbus Mono L",Monaco,"Courier New",Courier,monospace;font-size:12px;box-sizing:border-box;color:rgb(187,136,68)">"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 style="background-color:rgb(248,248,248);color:rgb(51,51,51);font-family:"Source Code Pro","Andale Mono WT","Andale Mono","Lucida Console","Lucida Sans Typewriter","DejaVu Sans Mono","Bitstream Vera Sans Mono","Liberation Mono","Nimbus Mono L",Monaco,"Courier New",Courier,monospace;font-size:12px"> </span><span class="gmail-n" style="color:rgb(51,51,51);font-family:"Source Code Pro","Andale Mono WT","Andale Mono","Lucida Console","Lucida Sans Typewriter","DejaVu Sans Mono","Bitstream Vera Sans Mono","Liberation Mono","Nimbus Mono L",Monaco,"Courier New",Courier,monospace;font-size:12px;box-sizing:border-box">wmts</span><span class="gmail-o" style="color:rgb(51,51,51);font-family:"Source Code Pro","Andale Mono WT","Andale Mono","Lucida Console","Lucida Sans Typewriter","DejaVu Sans Mono","Bitstream Vera Sans Mono","Liberation Mono","Nimbus Mono L",Monaco,"Courier New",Courier,monospace;font-size:12px;box-sizing:border-box;font-weight:bold">.</span><span class="gmail-n" style="color:rgb(51,51,51);font-family:"Source Code Pro","Andale Mono WT","Andale Mono","Lucida Console","Lucida Sans Typewriter","DejaVu Sans Mono","Bitstream Vera Sans Mono","Liberation Mono","Nimbus Mono L",Monaco,"Courier New",Courier,monospace;font-size:12px;box-sizing:border-box">xml</span><span style="background-color:rgb(248,248,248);color:rgb(51,51,51);font-family:"Source Code Pro","Andale Mono WT","Andale Mono","Lucida Console","Lucida Sans Typewriter","DejaVu Sans Mono","Bitstream Vera Sans Mono","Liberation Mono","Nimbus Mono L",Monaco,"Courier New",Courier,monospace;font-size:12px"> </span><span class="gmail-o" style="color:rgb(51,51,51);font-family:"Source Code Pro","Andale Mono WT","Andale Mono","Lucida Console","Lucida Sans Typewriter","DejaVu Sans Mono","Bitstream Vera Sans Mono","Liberation Mono","Nimbus Mono L",Monaco,"Courier New",Courier,monospace;font-size:12px;box-sizing:border-box;font-weight:bold">-</span><span class="gmail-n" style="color:rgb(51,51,51);font-family:"Source Code Pro","Andale Mono WT","Andale Mono","Lucida Console","Lucida Sans Typewriter","DejaVu Sans Mono","Bitstream Vera Sans Mono","Liberation Mono","Nimbus Mono L",Monaco,"Courier New",Courier,monospace;font-size:12px;box-sizing:border-box">of</span><span style="background-color:rgb(248,248,248);color:rgb(51,51,51);font-family:"Source Code Pro","Andale Mono WT","Andale Mono","Lucida Console","Lucida Sans Typewriter","DejaVu Sans Mono","Bitstream Vera Sans Mono","Liberation Mono","Nimbus Mono L",Monaco,"Courier New",Courier,monospace;font-size:12px"> </span><span class="gmail-n" style="color:rgb(51,51,51);font-family:"Source Code Pro","Andale Mono WT","Andale Mono","Lucida Console","Lucida Sans Typewriter","DejaVu Sans Mono","Bitstream Vera Sans Mono","Liberation Mono","Nimbus Mono L",Monaco,"Courier New",Courier,monospace;font-size:12px;box-sizing:border-box">WMTS>,</span></div><div>I tried running:</div><div>(1)</div><div> 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">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 </div><div><br></div><div><b>The response:</b> I did not get the xml as wmts.xml showed above</div><div><br></div><div>Note: Including the -of option (WCS), the outcome fails</div><div><br></div><div>(2)</div><div>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">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<br></div><div><br></div><div><b>The response:</b> There is no -spat option for gdal_translate</div><div><br></div><div><b>gdal_translate --help</b></div><div><br></div><div>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<br></div><div><br></div><div><br></div><div><br></div><div>[1] WCS GetCapabilities request </div><div><br></div><div><div class="gmail-line" style="color:rgb(0,0,0);font-family:monospace;font-size:13px"><span class="gmail-html-tag"><wcs:CoverageSummary></span></div><div class="gmail-opened" style="margin-left:1em;color:rgb(0,0,0);font-family:monospace;font-size:13px"><div class="gmail-line"><span class="gmail-html-tag"><ows:Title></span>CBERS-4_AWFI_170_111_21052020<span class="gmail-html-tag"></ows:Title></span></div><div class="gmail-line"><span class="gmail-html-tag"><ows:Abstract></span>Generated from GeoTIFF<span class="gmail-html-tag"></ows:Abstract></span></div><div class="gmail-folder" id="gmail-folder28"><div class="gmail-line"><span class="gmail-folder-button gmail-fold" style="display:inline-block;width:10px;background-position:initial;background-size:initial;background-repeat:initial;background-origin:initial;background-clip:initial;background-color:initial;vertical-align:bottom;height:10px"></span><span class="gmail-html-tag"><ows:Keywords></span></div><div class="gmail-opened" style="margin-left:1em"><div class="gmail-line"><span class="gmail-html-tag"><ows:Keyword></span>deter-amazonia-CBERS-4_AWFI_170_111_21052020<span class="gmail-html-tag"></ows:Keyword></span></div><div class="gmail-line"><span class="gmail-html-tag"><ows:Keyword></span>WCS<span class="gmail-html-tag"></ows:Keyword></span></div><div class="gmail-line"><span class="gmail-html-tag"><ows:Keyword></span>GeoTIFF<span class="gmail-html-tag"></ows:Keyword></span></div></div><div class="gmail-line"><span class="gmail-html-tag"></ows:Keywords></span></div></div><div class="gmail-line"><span class="gmail-html-tag"><wcs:CoverageId></span>deter-amazonia__CBERS-4_AWFI_170_111_21052020<span class="gmail-html-tag"></wcs:CoverageId></span></div><div class="gmail-line"><span class="gmail-html-tag"><wcs:CoverageSubtype></span>RectifiedGridCoverage<span class="gmail-html-tag"></wcs:CoverageSubtype></span></div><div class="gmail-folder" id="gmail-folder29"><div class="gmail-line"><span class="gmail-folder-button gmail-fold" style="display:inline-block;width:10px;background-position:initial;background-size:initial;background-repeat:initial;background-origin:initial;background-clip:initial;background-color:initial;vertical-align:bottom;height:10px"></span><span class="gmail-html-tag"><ows:WGS84BoundingBox></span></div><div class="gmail-opened" style="margin-left:1em"><div class="gmail-line"><span class="gmail-html-tag"><ows:LowerCorner></span>-63.280290728152664 -13.794296142085976<span class="gmail-html-tag"></ows:LowerCorner></span></div><div class="gmail-line"><span class="gmail-html-tag"><ows:UpperCorner></span>-53.624004749255874 -5.636644661059637<span class="gmail-html-tag"></ows:UpperCorner></span></div></div><div class="gmail-line"><span class="gmail-html-tag"></ows:WGS84BoundingBox></span></div></div><div class="gmail-folder" id="gmail-folder30"><div class="gmail-line"><span class="gmail-folder-button gmail-fold" style="display:inline-block;width:10px;background-position:initial;background-size:initial;background-repeat:initial;background-origin:initial;background-clip:initial;background-color:initial;vertical-align:bottom;height:10px"></span><span class="gmail-html-tag"><ows:BoundingBox<span class="gmail-html-attribute"> <span class="gmail-html-attribute-name">crs</span>="<span class="gmail-html-attribute-value"><a href="http://www.opengis.net/def/crs/EPSG/0/EPSG:4326">http://www.opengis.net/def/crs/EPSG/0/EPSG:4326</a></span>"</span>></span></div><div class="gmail-opened" style="margin-left:1em"><div class="gmail-line"><span class="gmail-html-tag"><ows:LowerCorner></span>-63.280290728152664 -13.794296142085976<span class="gmail-html-tag"></ows:LowerCorner></span></div><div class="gmail-line"><span class="gmail-html-tag"><ows:UpperCorner></span>-53.624004749255874 -5.636644661059637<span class="gmail-html-tag"></ows:UpperCorner></span></div></div><div class="gmail-line"><span class="gmail-html-tag"></ows:BoundingBox></span></div></div></div><div class="gmail-line" style="color:rgb(0,0,0);font-family:monospace;font-size:13px"><span class="gmail-html-tag"></wcs:CoverageSummary></span></div></div><div><br></div><div><br></div><div>Any comment will be very appreciated.</div><div><br></div><div>GDAL <a href="https://gdal.org/drivers/raster/wmts.html">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">https://gdal.org/drivers/raster/wcs.html</a> and could not find an example pointing to -spat option. I got an -spat option <br><div>to<a href="https://gdal.org/programs/gdal_grid.html?highlight=spat#cmdoption-gdal_grid-spat"> gdal_grid</a> though. </div><div><pre style="box-sizing:border-box;font-family:"Source Code Pro","Andale Mono WT","Andale Mono","Lucida Console","Lucida Sans Typewriter","DejaVu Sans Mono","Bitstream Vera Sans Mono","Liberation Mono","Nimbus Mono L",Monaco,"Courier New",Courier,monospace;font-size:12px;background-color:rgb(248,248,248);margin-top:0px;margin-bottom:0px;padding:12px;overflow:auto;line-height:1.5;color:rgb(51,51,51)"><span class="gmail-n" style="box-sizing:border-box">gdal_translate</span> <span class="gmail-s2" style="box-sizing:border-box;color:rgb(187,136,68)">"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 class="gmail-n" style="box-sizing:border-box">wmts</span><span class="gmail-o" style="box-sizing:border-box;font-weight:bold">.</span><span class="gmail-n" style="box-sizing:border-box">xml</span> <span class="gmail-o" style="box-sizing:border-box;font-weight:bold">-</span><span class="gmail-n" style="box-sizing:border-box">of</span> <span class="gmail-n" style="box-sizing:border-box">WMTS</span></pre><pre style="box-sizing:border-box;font-family:"Source Code Pro","Andale Mono WT","Andale Mono","Lucida Console","Lucida Sans Typewriter","DejaVu Sans Mono","Bitstream Vera Sans Mono","Liberation Mono","Nimbus Mono L",Monaco,"Courier New",Courier,monospace;font-size:12px;background-color:rgb(248,248,248);margin-top:0px;margin-bottom:0px;padding:12px;overflow:auto;line-height:1.5;color:rgb(51,51,51)"><span class="gmail-n" style="box-sizing:border-box">gdal_translate "WCS:<a href="http://www.dpi.inpe.br/fipcerrado-geoserver/ows/wcs?request=GetCapabilities">http://www.dpi.inpe.br/fipcerrado-geoserver/ows/wcs?request=GetCapabilities</a>" \</span></pre><pre style="box-sizing:border-box;font-family:"Source Code Pro","Andale Mono WT","Andale Mono","Lucida Console","Lucida Sans Typewriter","DejaVu Sans Mono","Bitstream Vera Sans Mono","Liberation Mono","Nimbus Mono L",Monaco,"Courier New",Courier,monospace;font-size:12px;background-color:rgb(248,248,248);margin-top:0px;margin-bottom:0px;padding:12px;overflow:auto;line-height:1.5;color:rgb(51,51,51)"><span class="gmail-n" style="box-sizing:border-box">-oo CoverageName=deter-cerrado:CBERS-4_AWFI_161_099_13072018 CBERS-4_AWFI_161_099_13072018.xml -of WCS<br><br></span></pre><pre style="box-sizing:border-box;font-family:"Source Code Pro","Andale Mono WT","Andale Mono","Lucida Console","Lucida Sans Typewriter","DejaVu Sans Mono","Bitstream Vera Sans Mono","Liberation Mono","Nimbus Mono L",Monaco,"Courier New",Courier,monospace;font-size:12px;background-color:rgb(248,248,248);margin-top:0px;margin-bottom:0px;padding:12px;overflow:auto;line-height:1.5;color:rgb(51,51,51)"><span class="gmail-n" style="box-sizing:border-box"><br></span></pre></div></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Mon, Nov 9, 2020 at 5:00 PM <<a href="mailto:gdal-dev-request@lists.osgeo.org">gdal-dev-request@lists.osgeo.org</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">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" rel="noreferrer" 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: <1659214.2tRo4lr4Oo@even-i700><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" rel="noreferrer" 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" rel="noreferrer" 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" rel="noreferrer" 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" rel="noreferrer" 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" rel="noreferrer" 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" rel="noreferrer" 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&" rel="noreferrer" 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>
> gdal-dev@.osgeo<br>
<br>
> <a href="https://lists.osgeo.org/mailman/listinfo/gdal-dev" rel="noreferrer" 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" rel="noreferrer" 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" rel="noreferrer" 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>
*****************************************<br>
</blockquote></div>