アフィン座標を緯度/経度に変換するにはどうすればよいですか？

私はGISを初めて使用します。

私はgdal土地利用/土地被覆地図を読み込むために使用しており、特定の土地被覆タイプの緯度/経度を選択して、緯度/経度でのみ表現される別のデータセットにインデックスを付ける必要があります。Unfortuantelyは、私は具体的には、geotransformから私に与えられたx、y座標の形式を理解していないoriginXとoriginY、以下の：

geotransform = dataset.GetGeoTransform()
originX = geotransform[0]
originY = geotransform[3]

これらの値を印刷すると、のような座標が得られます(447466.693808, 4952570.40529)。これらは元の緯度と経度にどのように関係していますか？

編集：

ここに私が探していたものを手に入れた簡単なpythonの例があります：

srs = osr.SpatialReference()
srs.ImportFromWkt(dataset.GetProjection())

srsLatLong = srs.CloneGeogCS()
ct = osr.CoordinateTransformation(srs,srsLatLong)
print ct.TransformPoint(originX,originY)

盗難：tolatlong.py

gdal_translateは、次の方法を使用して、他の投影法（この場合はEPSG：4326が必要）にデータを再投影します。

gdal_translate -a_srs epsg:4326 srcfile new_file 

または、gdaltrasformを使用してポイントを変換できます（そして、Python（？）からもアクセスできると確信しています）

geotransformはhttps://gdal.org/user/raster_data_model.htmlに文書化されています。これは、データセットから（x、y）座標を直接取得し、線形変換を適用して（u、v）を取得するという考え方です。

u = a*x + b*y
v = c*x + d*y

（これを線形変換の定義とすることができます）、uにgeotransform [0]を、vにgeotransform [3]を追加して原点をシフトします。これにより、（x、y）の「アフィン変換」が得られます。回転、スケールの変更、スキューエラーの多少の修正、および既知の座標系と一致するようにデータ固有の座標（x、y）を再配置することを実際に意図しています。結果は、投影座標を生成することになっています。 これは単に、（経度、緯度）を取り、既知の座標に変換する数学的手順があることを意味します。これは「投影」と呼ばれます。「投影解除」はその逆です。したがって、必要な投影法がわかっている場合は、アフィン変換された（x、y）座標に適用します 緯度と経度を取得します。

ところで、定数a、b、c、dの値は、geotransform配列の1、2、4、および5エントリによって与えられます。

以下を使用できます。

import gdal, numpy as n
def coord(file):
    padfTransform = file.GetGeoTransform()
    indices = n.indices(file.ReadAsArray().shape)
    xp = padfTransform[0] + indices[1]*padfTransform[1] + indices[1]*padfTransform[2]   
    yp = padfTransform[3] + indices[0]*padfTransform[4] + indices[0]*padfTransform[5]  
    return xp,yp
file = gdal.Open('Yourfile.tif') # A GeoTiff file
x,y = coord(file)

coordは、すべてのピクセルの経度（x）と緯度（y）を返します。座標はピクセルの左隅にあることに注意してください