Dzień dobry,
to mój pierwszy raz na tym forum, więc jeśli popełnię jakieś faux pas to z góry przepraszam.
Z GISu jestem absolutnym amatorem, więc pewnie popełniam jakąś głupotę, że nie chce działać, ale do brzegu.
Próbuję za pomocą gdalwarp z poziomu konsoli ściągnąć na dysk kawałek mapy topograficznej z geoportal.
W pierwszym kroku za pomocą gdal_translate tworzę sobie plik xml z danymi warstwy:
Kod:
gdal_translate "WMTS:https://mapy.geoportal.gov.pl/wss/service/WMTS/guest/wmts/TOPO?request=getCapabilities&service=WMTS" c:\_\tmp\geopTopo.xml -of WMTS
Powstaje plik o następującej treści:
Kod:
<GDAL_WMTS>
<GetCapabilitiesUrl>https://mapy.geoportal.gov.pl/wss/service/WMTS/guest/wmts/TOPO?request=getCapabilities&service=WMTS</GetCapabilitiesUrl>
<Layer>MAPA TOPOGRAFICZNA</Layer>
<Style>default</Style>
<TileMatrixSet>EPSG:2180</TileMatrixSet>
<InfoFormat>text/html</InfoFormat>
<DataWindow>
<UpperLeftX>139521.1561254875</UpperLeftX>
<UpperLeftY>790426.826687225</UpperLeftY>
<LowerRightX>880075.247647375</LowerRightX>
<LowerRightY>134423.4188497</LowerRightY>
</DataWindow>
<BandsCount>3</BandsCount> <------------ tu było 4 ale zmieniłem na 3, nie potrzebuję kanału alpha
<DataType>Byte</DataType>
<Cache /> <------------------------------ ten fragment usunąłem bo czasami jest przyczyną innego błędu ale tu zadam tylko jedno pytanie
<UnsafeSSL>true</UnsafeSSL>
<ZeroBlockHttpCodes>204,404</ZeroBlockHttpCodes>
<ZeroBlockOnServerException>true</ZeroBlockOnServerException>
</GDAL_WMTS>
Teraz korzystając z informacji powyżej chcę ściągnąć kawałek mapy za pomocą gdalwarp:
Kod:
gdalwarp -overwrite C:\_\tmp\geopTopo.xml C:\_\tmp\test5.jpeg -of jpeg -te 7437567 5536742 7441385 5539164 -te_srs EPSG:2178 -t_srs EPSG:2178 -ts 1236 746
i tu pojawia się problem. Dostaję taką informację zwrotną:
Kod:
Creating output file that is 1236P x 746L.
0ERROR 1: GDALWMS: Unable to download block 326, 410.
URL: Recv failure: Connection was reset
HTTP status code: 0, error: Recv failure: Connection was reset.
Add the HTTP status code to <ZeroBlockHttpCodes> to ignore this error (see http://www.gdal.org/frmt_wms.html).
ERROR 1: GDAL_WMS>, band 1: IReadBlock failed at X offset 326, Y offset 410: GDALWMS: Unable to download block 326, 410.
URL: Recv failure: Connection was reset
HTTP status code: 0, error: Recv failure: Connection was reset.
Add the HTTP status code to <ZeroBlockHttpCodes> to ignore this error (see http://www.gdal.org/frmt_wms.html).
ERROR 1: IReadBlock failed at X offset 0, Y offset 0: GDAL_WMS>, band 1: IReadBlock failed at X offset 326, Y offset 410: GDALWMS: Unable to download block 326, 410.
URL: Recv failure: Connection was reset
HTTP status code: 0, error: Recv failure: Connection was reset.
Add the HTTP status code to <ZeroBlockHttpCodes> to ignore this error (see http://www.gdal.org/frmt_wms.html).
Co robię źle?
Dodam, że jeżeli dokładnie te same kroki wykonam dla ortofotomapy to wszystko działa bez zarzutu.
gdalwarp - nie działa dla mapy topograficznej z geoportal
