尋ねた@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))
