다음 질문이 다소 어리 석다면 사과하지만, 나는이 전체 GIS에 아주 새로운 것입니다.
파이썬에서 gdal을 사용하여 투영 된 geoTiff 이미지를 WGS84로 변환하려고합니다. 다음과 비슷한 것을 사용하여 투영 된 GeoTiff 내의 점을 변환하는 프로세스를 개략적으로 설명하는 게시물을 찾았습니다.
from osgeo import osr, gdal
# get the existing coordinate system
ds = gdal.Open('path/to/file')
old_cs= osr.SpatialReference()
old_cs.ImportFromWkt(ds.GetProjectionRef())
# create the new coordinate system
wgs84_wkt = """
GEOGCS["WGS 84",
DATUM["WGS_1984",
SPHEROID["WGS 84",6378137,298.257223563,
AUTHORITY["EPSG","7030"]],
AUTHORITY["EPSG","6326"]],
PRIMEM["Greenwich",0,
AUTHORITY["EPSG","8901"]],
UNIT["degree",0.01745329251994328,
AUTHORITY["EPSG","9122"]],
AUTHORITY["EPSG","4326"]]"""
new_cs = osr.SpatialReference()
new_cs .ImportFromWkt(wgs84_wkt)
# create a transform object to convert between coordinate systems
transform = osr.CoordinateTransformation(old_cs,new_cs)
#get the point to transform, pixel (0,0) in this case
width = ds.RasterXSize
height = ds.RasterYSize
gt = ds.GetGeoTransform()
minx = gt[0]
miny = gt[3] + width*gt[4] + height*gt[5]
#get the coordinates in lat long
latlong = transform.TransformPoint(x,y)
내 질문은,이 점들을 변환하고 새로운 WGS84 GeoTiff 파일을 만들고 싶다면 이것이 가장 좋은 방법입니까? 1 단계의 작업과 같이 수행 할 기능이 있습니까?
감사!