Вибачте, якщо наступне питання є дещо дурним, але я лише ДУЖЕ новий для всієї цієї ГІС речі.
Я намагаюся перетворити деякі проектовані зображення GeoTiff в WGS84 за допомогою gdal в python. Я знайшов публікацію, яка окреслює процес перетворення точок в прогнозованих 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 крок?
Спасибі!