python - geotiff 画像から緯度と経度を計算する

Question

geotiff各ブロックを通常のブロック (たとえば、50 x 50 ピクセル) に分割した画像から、各ブロックの平均緯度と平均経度を計算するにはどうすればよいですか。

入力データは、たとえばhttp://eoimages.gsfc.nasa.gov/images/imagerecords/57000/57752/land_shallow_topo_2048.tifからダウンロードした画像のみです。

これは、次のように gdal を使用して Python で開くことができます。

import gdal
geotiff = gdal.Open ('land_shallow_topo_2048.tif')
colum_numbers,row_numbers,band_numbers=geotiff.RasterXSize,
                                       geotiff.RasterYSize,geotiff.RasterCount
print (colum_numbers,row_numbers,band_numbers)
2048 1024 3