@Aaron은 다음과 같이 물었습니다.
향상된 멀티 코어 및 멀티 스레드 작업을 위해 -multi 옵션을 사용하는 @wwnick의 대답에 대한 gdalwarp 버전을 찾고 싶습니다.
약간의 부인
gdalwarp비록 성능이 많이 나올 것이라고 확신하지는 않지만을 사용합니다 . 지금까지 나는 I / O 바인딩되었습니다.이 스크립트를 큰 래스터에서 많은 작은 부분으로 자르는 것은 CPU를 많이 사용하지 않는 것처럼 보이므로 병목 현상이 디스크에 기록되고 있다고 가정합니다. 타일 또는 이와 유사한 것을 동시에 다시 투영하려는 경우 변경 될 수 있습니다. 튜닝 팁이 있습니다 . 간단한 플레이는 나에게 아무런 개선을 가져 오지 않았으며 CPU는 결코 제한 요소가 아닌 것으로 보였습니다.
면책 조항을 제외 gdalwarp하고 래스터를 여러 개의 작은 타일로 분할하는 데 사용되는 스크립트가 있습니다 . 바닥 구분으로 인해 약간의 손실이있을 수 있지만 원하는 타일 수를 선택하여 처리 할 수 있습니다. 그것은 될 것입니다 n+1곳 n얻가 당신에 의해 분할 개수 tile_width와 tile_height변수.
import subprocess
import gdal
import sys
def gdalwarp(*args):
return subprocess.check_call(['gdalwarp'] + list(args))
src_path = sys.argv[1]
ds = gdal.Open(src_path)
try:
out_base = sys.argv[2]
except IndexError:
out_base = '/tmp/test_'
gt = ds.GetGeoTransform()
width_px = ds.RasterXSize
height_px = ds.RasterYSize
# Get coords for lower left corner
xmin = int(gt[0])
xmax = int(gt[0] + (gt[1] * width_px))
# get coords for upper right corner
if gt[5] > 0:
ymin = int(gt[3] - (gt[5] * height_px))
else:
ymin = int(gt[3] + (gt[5] * height_px))
ymax = int(gt[3])
# split height and width into four - i.e. this will produce 25 tiles
tile_width = (xmax - xmin) // 4
tile_height = (ymax - ymin) // 4
for x in range(xmin, xmax, tile_width):
for y in range(ymin, ymax, tile_height):
gdalwarp('-te', str(x), str(y), str(x + tile_width),
str(y + tile_height), '-multi', '-wo', 'NUM_THREADS=ALL_CPUS',
'-wm', '500', src_path, out_base + '{}_{}.tif'.format(x, y))
