아핀 좌표를 위도 / 경도로 변환하려면 어떻게해야합니까?

저는 GIS를 처음 접했습니다.

나는 gdal토지 사용 / 토지 표지 맵을 읽는 데 사용 하고 있으며 특정 토지 표지 유형의 위도 / 경도를 선택하여 다른 데이터 세트로 색인화하여 위도 / 경도로만 표현해야합니다. 불행히도, 나는 지형 변환에서 나에게 주어진 x 및 y 좌표의 형태, 특히 originX및 originY아래를 이해하지 못합니다 .

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

이 값을 인쇄하면과 같은 좌표가 표시 (447466.693808, 4952570.40529)됩니다. 이것들은 원래 위도와 경도와 어떤 관련이 있습니까?

편집하다:

여기 내가 찾던 것을 얻은 간단한 파이썬 예제가 있습니다.

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 을 사용 하여 점을 변환 할 수 있습니다 (그리고 파이썬에서도 액세스 할 수 있다고 확신합니까?)

지리 변환은 https://gdal.org/user/raster_data_model.html에 문서화되어 있습니다 . 아이디어는 데이터 세트에서 (x, y) 좌표를 직선으로 가져 와서 선형 변환을 적용하여 (u, v)를 얻는 것입니다.

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

( 이를 선형 변환 의 정의 로 사용할 수 있습니다 ) 그런 다음 geotransform [0]을 u에 추가하고 geotransform [3]을 v에 추가하여 원점을 이동하십시오. 그러면 (x, y)의 "병렬 변환"이 제공됩니다. 실제로 회전, 스케일 변경, 약간의 오차에 대해 약간 수정하고 알려진 좌표계와 일치하도록 데이터 특정 좌표 (x, y)를 다시 배치하려고합니다. 결과적으로 투영 된 좌표 가 생성 됩니다. 이것은 단순히 수학적 절차 (경도, 위도)를 알려진 좌표로 바꾸는 것을 의미합니다. 이것을 "투영"이라고합니다. "투영"은 그 반대입니다. 따라서 어떤 투영이 필요한지 아는 경우 아핀 변환 된 (x, y) 좌표에 적용합니다 위도와 경도를 얻을 수 있습니다.

그런데 상수 a, b, c, d의 값은 지오 변환 배열의 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)를 반환합니다. 좌표는 픽셀의 왼쪽 모서리에 있습니다.