2

大きな GeoTiff ファイル (1.4GB、緯度 4 度 x 経度 8 度) がたくさんあります。

それぞれを 1 度の緯度 x 1 度の長さのタイルにスライスする必要があります (それぞれにわずかな超過があります)。

GeoTiff ファイルの 1 つの gdalinfo 出力を次に示します。

Driver: GTiff/GeoTIFF
Files: index.tif
Size is 38401, 19201
Coordinate System is:
GEOGCS["WGS 84",
    DATUM["WGS_1984",
        SPHEROID["WGS 84",6378137,298.257223563,
            AUTHORITY["EPSG","7030"]],
        AUTHORITY["EPSG","6326"]],
    PRIMEM["Greenwich",0],
    UNIT["degree",0.0174532925199433],
    AUTHORITY["EPSG","4326"]]
Origin = (-112.000104166666674,56.000104166666667)
Pixel Size = (0.000208333333333,-0.000208333333333)
Metadata:
  AREA_OR_POINT=Area
Image Structure Metadata:
  INTERLEAVE=BAND
Corner Coordinates:
Upper Left  (-112.0001042,  56.0001042) (112d 0' 0.38"W, 56d 0' 0.37"N)
Lower Left  (-112.0001042,  51.9998958) (112d 0' 0.38"W, 51d59'59.62"N)
Upper Right (-103.9998958,  56.0001042) (103d59'59.62"W, 56d 0' 0.37"N)
Lower Right (-103.9998958,  51.9998958) (103d59'59.62"W, 51d59'59.62"N)
Center      (-108.0000000,  54.0000000) (108d 0' 0.00"W, 54d 0' 0.00"N)
Band 1 Block=38401x1 Type=Int16, ColorInterp=Gray
  NoData Value=-32767

サイズ以外は全て同じにしたいです。したがって、それぞれを 1 度 x 1 度の 32 個のタイルにスライスする必要があります。

たとえば
、-112.0001042、56.0001042 から -110.9998958、54.9998958 (基本的に -112, 56 から -111, 55) をカバーします。

gdal_retile.py や gdal_grid などのユーティリティが表示されますが、進歩していません。これに適したツール/コマンドラインオプションは何ですか?

4

1 に答える 1

3

以下のように、Python ループで gdal_translate を使用できます。

import os

for ulx in range(-112,-103):
    for uly in range(51,56):
        os.system('gdal_translate -projwin '+str(ulx)+' '+str(uly)+' '+str(ulx+1)+' '+str(uly+1)+' index.tif index_tile_ulx_'+str(ulx)+'_uly_'+str(uly)+'.tif')
于 2013-11-26T14:18:07.743 に答える