Относительный путь с zip в файле vrt (vsizip)?

Я пытаюсь создать файл vrt с большим количеством заархивированных файлов tiff. Я обнаружил, что эта xml-команда добавляет tiff-файл внутри zip-файла в vrt:

<SourceFilename relativeToVRT="1">/vsizip/c:/data/DSM_633_54_TIF_UTM32-ETRS89.zip/DSM_1km_6333_546.tif</SourceFilename>

Это прекрасно работает, но проблема в том, что это работает только с абсолютным путем. Поскольку я работаю в докере, я не могу использовать абсолютный путь, только относительный.

Если я просто использую относительный путь, он не может найти изображения (все они черные и имеют значения 1e [огромное число], когда я открываю его в QGis).

Есть ли способ использовать относительный путь при добавлении TIFF-изображения из ZIP-файла в VRT??


person Little geek    schedule 02.02.2017    source источник
comment
Только что нашел это: erouault.blogspot.dk/ 2011/12/ выглядит интересно   -  person Poul K. Sørensen    schedule 03.02.2017


Ответы (1)


Я думаю, вы можете использовать абсолютные пути:

Вот пример (обратите внимание на двойную косую черту).

pksorensen@swarm-master-4D58E13F-0:~$ docker -H 172.16.0.5:2375 run --rm -v $(docker -H 172.16.0.5:2375 volume create -d azurefile -o share=dhm-danmarks-hoejdemodel):/data  geodata/gdal gdalinfo /vsizip//data/DSM_633_54_TIF_UTM32-ETRS89.zip/DSM_1km_6333_546.tif
Driver: GTiff/GeoTIFF
Files: /vsizip//data/DSM_633_54_TIF_UTM32-ETRS89.zip/DSM_1km_6333_546.tif
Size is 2500, 2500
Coordinate System is:
PROJCS["ETRS89 / UTM zone 32N",
    GEOGCS["ETRS89",
        DATUM["European_Terrestrial_Reference_System_1989",
            SPHEROID["GRS 1980",6378137,298.257222101,
                AUTHORITY["EPSG","7019"]],
            TOWGS84[0,0,0,0,0,0,0],
            AUTHORITY["EPSG","6258"]],
        PRIMEM["Greenwich",0,
            AUTHORITY["EPSG","8901"]],
        UNIT["degree",0.0174532925199433,
            AUTHORITY["EPSG","9122"]],
        AUTHORITY["EPSG","4258"]],
    PROJECTION["Transverse_Mercator"],
    PARAMETER["latitude_of_origin",0],
    PARAMETER["central_meridian",9],
    PARAMETER["scale_factor",0.9996],
    PARAMETER["false_easting",500000],
    PARAMETER["false_northing",0],
    UNIT["metre",1,
        AUTHORITY["EPSG","9001"]],
    AXIS["Easting",EAST],
    AXIS["Northing",NORTH],
    AUTHORITY["EPSG","25832"]]
Origin = (546000.000000000000000,6334000.000000000000000)
Pixel Size = (0.400000000000000,-0.400000000000000)
Metadata:
  AREA_OR_POINT=Area
Image Structure Metadata:
  COMPRESSION=DEFLATE
  INTERLEAVE=BAND
Corner Coordinates:
Upper Left  (  546000.000, 6334000.000) (  9d45'36.85"E, 57d 8'49.01"N)
Lower Left  (  546000.000, 6333000.000) (  9d45'36.19"E, 57d 8'16.67"N)
Upper Right (  547000.000, 6334000.000) (  9d46'36.34"E, 57d 8'48.65"N)
Lower Right (  547000.000, 6333000.000) (  9d46'35.66"E, 57d 8'16.31"N)
Center      (  546500.000, 6333500.000) (  9d46' 6.26"E, 57d 8'32.66"N)
Band 1 Block=256x256 Type=Float32, ColorInterp=Gray
  NoData Value=-9999

Из-за этого это означает, что мы предполагаем, что сгенерированный файл vrt будет иметь одинаковые привязки тома для любого последующего вызова, но я думаю, что это должно быть хорошо :)

person Poul K. Sørensen    schedule 03.02.2017