端を発したGetGeoTransformにマイクToewsポインタ、私はサポートされている任意の地理参照ラスター(と思う)の世界ファイルを構築する小さなGDAL Pythonスクリプトを作成するために管理しました。完全なコードは次のとおりgdal-makeworld.py
です。基本的なビットは次のとおりです。
geotransform = dataset.GetGeoTransform()
if geotransform is not None:
x, x_size, x_rot, y, y_rot, y_size = geotransform
# correct for centre vs corner of pixel
x = x_size/2+x
y = y_size/2+y
world_file.write('%s\n' % x_size)
world_file.write('%s\n' % x_rot)
world_file.write('%s\n' % y_rot)
world_file.write('%s\n' % y_size)
world_file.write('%s\n' % x)
world_file.write('%s\n' % y)
world_file.close()
''' geotransform tuple key:
[0] /* top left x */
[1] /* w-e pixel resolution */
[2] /* rotation, 0 if image is "north up" */
[3] /* top left y */
[4] /* rotation, 0 if image is "north up" */
[5] /* n-s pixel resolution */
'''
gdalcopyproj.py
私が出発点として使用した執筆のためのSchuyler Erleにさらに感謝します。
ハットチップ@AlisterH、「ピクセルの中心とコーナーの両方に対して正しい」、2019-05-30