GDAL и WCS

Ответить
Andrey Zhukov
Гуру
Сообщения: 838
Зарегистрирован: 10 дек 2009, 23:24
Репутация: 169
Ваше звание: старик-гисовик
Откуда: Москва
Контактная информация:

GDAL и WCS

Сообщение Andrey Zhukov » 20 фев 2015, 16:25

Пытаемся забрать фрагмент (по subset) WCS через GDAL utils. Просто по url все работает, а с gdal начинается беда.
Запускаем gdal_translate -of GTIFF get5.xml wcs.tif
ошибка gdal_translate :
ERROR 1: Could not recognize crs EPSG:4326&version=2.0.1&request=GetCoverage&coverageId=srtm_full&subset=Lat(40,41)&subset=Long(40,41)&format=geotiff&Response_CRS=EPSG:4326

GDALOpen failed - 1
Could not recognize crs EPSG:4326&version=2.0.1&request=GetCoverage&coverageId=srtm_full&subset=Lat(40,41)&subset=Long(40,41)&format=geotiff&Response_CRS=EPSG:4326
gdalinfo тоже выдает ошибку.
ошибка gdalinfo :
ERROR 1: Could not recognize crs EPSG:4326&service=WCS&version=2.0.1&request=GetCoverage&coverageId=srtm_full&subset=Lat(40,41)&subset=Long(40,41)&format=geotiff&Response_CRS=EPSG:4326&CRS=EPSG:4326

gdalinfo failed - unable to open 'get5.xml'.
Пробовали и через vrt (по http://grasswiki.osgeo.org/wiki/WCS). Тоже ерунда получается: созданный vrt QGIS показывает как черный квадрат.
Пока что решили, что проще забирать нужные куски тупо по url.
Никто подобным не заморачивался? Может, мы что-то не так делаем?

Аватара пользователя
Дмитрий Барышников
Гуру
Сообщения: 2572
Зарегистрирован: 17 ноя 2009, 19:17
Репутация: 261
Откуда: Москва

Re: GDAL и WCS

Сообщение Дмитрий Барышников » 20 фев 2015, 19:00

Хотелось бы увидеть пример файла get5.xml

Andrey Zhukov
Гуру
Сообщения: 838
Зарегистрирован: 10 дек 2009, 23:24
Репутация: 169
Ваше звание: старик-гисовик
Откуда: Москва
Контактная информация:

Re: GDAL и WCS

Сообщение Andrey Zhukov » 24 фев 2015, 09:34

Тот файл был утерян, но шишки набивать продолжаем.

XML:

Код: Выделить всё

<WCS_GDAL>
  <ServiceURL>http://172.24.3.223:8080/geoserver/wcs?</ServiceURL>
  <CoverageName>sfdem</CoverageName>
  <GetCoverageExtra>%26service=WCS%26version=2.0.1%26request=GetCoverage%26coverageId=sfdem%26subset=E(600000,602000)%26subset=N(4920000,4922000)%26format=geotiff</GetCoverageExtra>
</WCS_GDAL>
Ответ:
ERROR 1: Could not recognize crs EPSG:26713&service=WCS&version=2.0.1&request=GetCoverage&coverageId=sfdem&subset=E(600000,602000)&subset=N(4920000,4922000)&format=geotiff

gdalinfo failed - unable to open 'get6.xml'.
Пробовали в xml добавлять crs, ошибка все та же.

Аватара пользователя
Дмитрий Барышников
Гуру
Сообщения: 2572
Зарегистрирован: 17 ноя 2009, 19:17
Репутация: 261
Откуда: Москва

Re: GDAL и WCS

Сообщение Дмитрий Барышников » 24 фев 2015, 17:48

А есть что-то публичное? А то так не проверить/отладить.
ERROR 1: Failed to connect to 172.24.3.223 port 8080: Network is unreachable

Andrey Zhukov
Гуру
Сообщения: 838
Зарегистрирован: 10 дек 2009, 23:24
Репутация: 169
Ваше звание: старик-гисовик
Откуда: Москва
Контактная информация:

Re: GDAL и WCS

Сообщение Andrey Zhukov » 24 фев 2015, 18:06

У нас наружу не смотрит, но нашел подобный:
http://172.24.3.223:8080
надо поменять на
http://map.ices.dk

Там такой же сервис болтается. Его мы тоже пробовали, ошибка та же.

Аватара пользователя
Дмитрий Барышников
Гуру
Сообщения: 2572
Зарегистрирован: 17 ноя 2009, 19:17
Репутация: 261
Откуда: Москва

Re: GDAL и WCS

Сообщение Дмитрий Барышников » 24 фев 2015, 20:13

Приведенный публичный адрес у меня отработал:
$ gdalinfo wcs.xml
Driver: WCS/OGC Web Coverage Service
Files: wcs.xml
Size is 634, 477
Coordinate System is:
PROJCS["NAD27 / UTM zone 13N",
GEOGCS["NAD27",
DATUM["North_American_Datum_1927",
SPHEROID["Clarke 1866",6378206.4,294.9786982139006,
AUTHORITY["EPSG","7008"]],
AUTHORITY["EPSG","6267"]],
PRIMEM["Greenwich",0],
UNIT["degree",0.0174532925199433],
AUTHORITY["EPSG","4267"]],
PROJECTION["Transverse_Mercator"],
PARAMETER["latitude_of_origin",0],
PARAMETER["central_meridian",-105],
PARAMETER["scale_factor",0.9996],
PARAMETER["false_easting",500000],
PARAMETER["false_northing",0],
UNIT["metre",1,
AUTHORITY["EPSG","9001"]],
AUTHORITY["EPSG","26713"]]
Origin = (589980.000000000000000,4928010.000000000000000)
Pixel Size = (30.000000000000000,-30.000000000000000)
Corner Coordinates:
Upper Left ( 589980.000, 4928010.000) (103d52' 5.32"W, 44d30' 6.31"N)
Lower Left ( 589980.000, 4913700.000) (103d52'14.26"W, 44d22'22.60"N)
Upper Right ( 609000.000, 4928010.000) (103d37'44.24"W, 44d29'56.87"N)
Lower Right ( 609000.000, 4913700.000) (103d37'55.07"W, 44d22'13.21"N)
Center ( 599490.000, 4920855.000) (103d44'59.72"W, 44d26' 9.97"N)
Band 1 Block=634x477 Type=Float32, ColorInterp=Undefined
Сам файл wcs.xml
<WCS_GDAL>
<ServiceURL>http://map.ices.dk/geoserver/wcs?</ServiceURL>
<CoverageName>sf:sfdem</CoverageName>
</WCS_GDAL>
После выполнения запроса GDAL лорписал в файл доп. информацию
<WCS_GDAL>
<ServiceURL>http://map.ices.dk/geoserver/wcs?</ServiceURL>
<CoverageName>sf:sfdem</CoverageName>
<CoverageOffering>
<description>Generated from sfdem</description>
<name>sf:sfdem</name>
<label>sfdem is a Tagged Image File Format with Geographic information</label>
<lonLatEnvelope srsName="urn:ogc:def:crs:OGC:1.3:CRS84">
<pos>-103.87108701853181 44.370187074132616</pos>
<pos>-103.62940739432703 44.5016011535299</pos>
</lonLatEnvelope>
<keywords>
<keyword>WCS</keyword>
<keyword>sfdem</keyword>
<keyword>sfdem</keyword>
</keywords>
<domainSet>
<spatialDomain>
<Envelope srsName="EPSG:26713">
<pos>589980.0 4913700.0</pos>
<pos>609000.0 4928010.0</pos>
</Envelope>
<RectifiedGrid dimension="2" srsName="EPSG:26713">
<limits>
<GridEnvelope>
<low>0 0</low>
<high>633 476</high>
</GridEnvelope>
</limits>
<axisName>E</axisName>
<axisName>N</axisName>
<origin>
<pos>589995.0 4927995.0</pos>
</origin>
<offsetVector>30.0 0.0</offsetVector>
<offsetVector>0.0 -30.0</offsetVector>
</RectifiedGrid>
</spatialDomain>
</domainSet>
<rangeSet>
<RangeSet>
<name>sfdem</name>
<label>sfdem is a Tagged Image File Format with Geographic information</label>
<axisDescription>
<AxisDescription>
<name>Band</name>
<label>Band</label>
<values>
<singleValue>1</singleValue>
</values>
</AxisDescription>
</axisDescription>
</RangeSet>
</rangeSet>
<supportedCRSs>
<requestResponseCRSs>EPSG:26713</requestResponseCRSs>
</supportedCRSs>
<supportedFormats nativeFormat="GeoTIFF">
<formats>ARCGRID</formats>
<formats>GeoTIFF</formats>
<formats>GIF</formats>
<formats>GTOPO30</formats>
<formats>IMAGEMOSAIC</formats>
<formats>JPEG</formats>
<formats>PNG</formats>
<formats>TIFF</formats>
</supportedFormats>
<supportedInterpolations default="nearest neighbor">
<interpolationMethod>nearest neighbor</interpolationMethod>
<interpolationMethod>bilinear</interpolationMethod>
<interpolationMethod>bicubic</interpolationMethod>
</supportedInterpolations>
</CoverageOffering>
<PreferredFormat>GeoTIFF</PreferredFormat>
<BandIdentifier>Band</BandIdentifier>
<BandCount>1</BandCount>
<BandType>Float32</BandType>
</WCS_GDAL>

Andrey Zhukov
Гуру
Сообщения: 838
Зарегистрирован: 10 дек 2009, 23:24
Репутация: 169
Ваше звание: старик-гисовик
Откуда: Москва
Контактная информация:

Re: GDAL и WCS

Сообщение Andrey Zhukov » 25 фев 2015, 09:27

Да, без дополнительных параметров все работает.
Но вот с этим параметром отваливается:

Код: Выделить всё

<GetCoverageExtra>%26service=WCS%26version=2.0.1%26request=GetCoverage%26coverageId=sfdem%26subset=E(600000,602000)%26subset=N(4920000,4922000)%26format=geotiff</GetCoverageExtra>

Andrey Zhukov
Гуру
Сообщения: 838
Зарегистрирован: 10 дек 2009, 23:24
Репутация: 169
Ваше звание: старик-гисовик
Откуда: Москва
Контактная информация:

Re: GDAL и WCS

Сообщение Andrey Zhukov » 25 фев 2015, 09:55

Нашли компромиссное решение: http://grasswiki.osgeo.org/wiki/WCS
Немного подпилили xml:

Код: Выделить всё

<WCS_GDAL>
  <ServiceURL>http://demo.opengeo.org/geoserver/wcs?</ServiceURL>
  <CoverageName>Img_Sample</CoverageName>
 <Version>1.0.0</Version>
</WCS_GDAL>
Через vrt прекрасно работает, в том числе запрос на нарезку.
Важный момент: написанный на gdal.org комментарий про Version по умолчанию неверный. Version надо указывать обязательно!

Ответить

Вернуться в «GDAL/OGR»

Кто сейчас на конференции

Сейчас этот форум просматривают: нет зарегистрированных пользователей и 22 гостя