2017-04-12 15 views
0

でnumpyの配列位置にGEOTIFFのxy共ords変換私は以下のリンクに記載されているように、私はnumpyの配列に読み込まれたGeoTIFFのファイルがありますパイソン

Working with TIFFs (import, export) in Python using numpy

そのGEOTIFF配列のサイズを私は(465,465)を持っており、gdalinfoを使ってファイルのメタデータを取得しています.WGS84をCRSとして使用しています。

私がファイルとしたいのは、QGISとGdalinfoに表示されているx y lat lon座標を、インポートされたnumpy配列のポイントの実際の位置に変換することです。 GDALのドキュメントから

gt = ds.GetGeoTransform()

:あなたはあなたがそれを得ることができる開かれたGDALのデータセットから、geotransformを使用する必要が

答えて

0

アフィン変換が返さ6つの係数で構成されていby ピクセル/ライン座標を ジオリファレンス空間にマップするGDALDataset :: GetGeoTransform()は、次の関係を使用します。

あなたのラスタが回転していない場合

Xgeo = GT(0) + Xpixel*GT(1) + Yline*GT(2) 
Ygeo = GT(3) + Xpixel*GT(4) + Yline*GT(5) 

http://www.gdal.org/gdal_datamodel.html、その非常に単純な、ただ決議により、原点と除算を引きます。

このようなアフィン変換を実行できるライブラリがいくつかあります。たとえば、適切な名前の 'affine'ライブラリです。その使用方法の詳細については、

http://www.perrygeo.com/python-affine-transforms.html