UPDATE: We seem to have found where the trouble lies, but sadly it is a problem that is built into GDAL's handling of WCS requests and it doesn't seem like we can bypass it. Instead we are going to write our own WCS client to fetch the GeoTiff data and then load in into GDAL from memory. See "EDIT 4" at the bottom for details.
The trouble I'm having is in writing the xml string which I send to the server. I either get an error message from the server when I try to fetch the Dataset src or it passes but I've apparently fetched all the data on the server (which can be seen on the size of the raster) despite wanting just a piece of it.
using(OSGeo.GDAL.Dataset src = OSGeo.GDAL.Gdal.Open(xml, OSGeo.GDAL.Access.GA_ReadOnly))
{
int rasters = src.RasterCount;
for(int i = 1; i <= rasters; ++i)
{
using(OSGeo.GDAL.Band raster = src.GetRasterBand(i))
{
//Yet to be written code
}
}
}
I've been unable to find a list of which attributes that can be sent to the server. The attributes I've tried are all attributes which I've seen other people try while I've been googling for help.
Some attributes I've seen different alternatives for, and I've tried each. For projection I've tried SRS, CRS and Projection. For image size I've tried WIDTH and HEIGHT and also SizeX and SizeY. For bounding box I've tried BoundingBox and BBOX.
This is one way that I've attempted to construct the xml string:
string xml = "<WCS_GDAL>";
xml += "<ServiceURL>http://example.com/geoserver/wcs?</ServiceURL>";
xml += "<VERSION>1.1.0</VERSION>";*
xml += "<REQUEST>GetCoverage</REQUEST>";
xml += "<CoverageName>coverage_name</CoverageName>";
xml += "<SRS>EPSG:3006</SRS>";
xml += "<CRS>EPSG:3006</CRS>";
xml += "<Projection>EPSG:3006</Projection>";
xml += "<GetCoverageExtra>";
{
xml += "&CRS=EPSG:3006";
xml += "&BBOX=372500.0,6164500.0,373000.0,6165000.0" + ",EPSG:3006";**
xml += "&WIDTH=200&HEIGHT=200";
xml += "&FORMAT=GTiff";
}
xml += "</GetCoverageExtra>";
xml += "</WCS_GDAL>";
*I've also tried 1.0.0.
**1.1.0 requires a projection to be part of the bounding box while for 1.0.0 GDAL doesn't allow more than 4 entries in the list.
This is another way that I tried to do it. In this case by affecting the final ServiceURL directly. GDAL will automatically append version and request type before the first ampersand, so I had to include them in the XML so that it doesn't append the wrong ones. GDAL will also cast an exception if CoverageName is not included in the XML.
string link = "http://example.com/geoserver/wcs?";
//link += "&VERSION=1.1.0";
//link += "&REQUEST=GetCoverage";
link += "&coverage=coverage_name";
link += "&CRS=EPSG:3006";
link += "&BBOX=372500.0,6164500.0,373000.0,6165000.0,EPSG:3006";
link += "&WIDTH=200";
link += "&HEIGHT=200";
link += "&FORMAT=GTiff";
string xml = "<WCS_GDAL>";
xml += "<ServiceURL>";
xml += link;
xml += "</ServiceURL>";
xml += "<VERSION>1.1.0</VERSION>";*
xml += "<REQUEST>GetCoverage</REQUEST>";
xml += "<CoverageName>coverage_name</CoverageName>";
xml += "</WCS_GDAL>";
If "<REQUEST>GetCoverage</REQUEST>" isn't included it seems that GDAL defaults to DescribeCoverage. (EDIT: This is incorrect. GDAL doesn't allow you to do DescribeCoverage. GDAL automatically does DescribeCoverage before it does GetCoverage.) And if version is not included GDAL seems to default to 1.0.0.
Depending on what attributes I include, what values I've given them and what version I specify I get a few different Java exceptions from the server. Sometimes it says java.util.ArrayList cannot be cast to org.geotools.geometry.GeneralEnvelope (with version 1.0.0) and sometimes it says java.lang.ClassCastException: [Ljava.lang.String; cannot be cast to java.lang.String.
The exception that I think managed to go the furthest and which I can't seem to get past is org.geotools.coverage.grid.InvalidGridGeometryException: Unspecified image's size..
When I use the version that directly affects the link I can bypass the problem of GDAL not allowing 5 entries for the bounding box with version 1.0.0, which makes the server not cast the ArrayList to GeneralEnvelope cast exception. However, it seems that it then completely ignores the bounding box and returns all the data on the server. If I switch version to 1.1.0 then I get the unspecified image's size exception.
There should be a version 2.0.0 of WCS, but considering the problems I'm having and the fact that I have yet to see any examples of someone else using it has made me unwilling to try it right now. EDIT: Scratch that, GDAL doesn't seem to support it... At least that's what the exceptions tells me.
In short: I need help with finding a way to fetch a GeoTiff of a small area via WCS.
EDIT: This is the XML which I'm currently trying to use, as asked for by iant:
<WCS_GDAL>
<ServiceURL>http://example.com/geoserver/wcs?</ServiceURL>
<SERVICE>WCS</SERVICE>
<VERSION>1.1.0</VERSION>
<REQUEST>GetCoverage</REQUEST>
<CoverageName>coverage_name</CoverageName>
<CRS>urn:ogc:def:crs:EPSG:3006</CRS>
<GetCoverageExtra>&BoundingBox=372500.0,6164500.0,373000.0,6165000.0,urn:ogc:def:crs:EPSG:3006&WIDTH=200&HEIGHT=200&RESX=200&RESY=200&SizeX=200&SizeY=200&FORMAT=image/GeoTiff</GetCoverageExtra>
</WCS_GDAL>
As you can see, I've resorted to using a lot of redundancy in hope that one of the attributes will be right. This XML isn't written like this btw. I do as I showed above where I build them in a string. Though some of the attributes and values are different since I wrote that.
EDIT 2: iant actually asked for the generated URL string. Here it is:
/geoserver/wcs?SERVICE=WCS&VERSION=1.1.0&REQUEST=GetCoverage&IDENTIFIER=coverage_name&FORMAT=image%2Ftiff&BOUNDINGBOX=200025,7699925,200075,7699975,urn:ogc:def:crs:EPSG::3006&RangeSubset=contents
It seems like GDAL has messed up the URL big time. It's supposed to include everything in the GetCoverageExtra tags straight of, from what I hear. Yet it doesn't have the WIDTH & HEIGHT... It also messed up the values in the bounding box for some reason.
EDIT 3: After URL-escaping the ampersands in the GetCoverageExtra tag the URL is more complete:
/geoserver/wcs?SERVICE=WCS&VERSION=1.1.0&REQUEST=GetCoverage&IDENTIFIER=swe%3Anh_riks_Sweref_99_TM_geotiff&FORMAT=image%2Ftiff&BOUNDINGBOX=200025,7699925,200075,7699975,urn:ogc:def:crs:EPSG::3006&RangeSubset=contents&BOUNDINGBOX=372500.0,6164500.0,373000.0,6165000.0,urn:ogc:def:crs:EPSG:3006&WIDTH=200
Something is still off though. There's two bounding boxes! Both the one which is wrong, which seems like some kind of default value, and the one which I wrote. WIDTH is included, but not HEIGHT, however I found that this was my own fault because I screwed up the ampersand before HEIGHT. After fixing it I get the following exception from the server: java.util.ArrayList cannot be cast to net.opengis.ows11.BoundingBoxType. The default value seems to have double colons between EPSG and 3006, but even after adding another colon I get the same error. I also tried to remove the ".0" from the bounding box values.
EDIT 4 (final?): By looking at the GDAL source code we seem to have found what GDAL is doing wrong. When opening a Dataset via WCS GDAL first automatically sends a DescribeCoverage to the server and uses the returned data to set a bounding box. This means that the bounding box is set to cover all the data in the coverage I'm asking for. This bounding box is then sent with the GetCoverage request regardless of what data you try to give GDAL and the server seems to be unable to handle receiving two bounding boxes in a single request. So there seems to be no way of getting a specified area in WCS via GDAL. My colleague is writing a separate class to fetch specific areas from the server without GDAL and then load that data into GDAL from memory. It seems to be the only way for now, other than compiling our own version of GDAL.
