Làm cách nào để chuyển đổi tọa độ affine thành lat / lng?

Tôi cực kỳ mới đối với GIS.

Tôi đang sử dụng gdalđể đọc trong bản đồ sử dụng đất / bản đồ và tôi cần chọn lat / lng của một số loại che phủ đất nhất định để lập chỉ mục thành một tập dữ liệu khác nhau chỉ được thể hiện bằng lat / lng. Thật không may, tôi không hiểu dạng tọa độ x và y được cung cấp cho tôi từ geotransform, cụ thể là originXvà originYbên dưới:

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

In những giá trị này cho tôi tọa độ như thế nào (447466.693808, 4952570.40529). Làm thế nào để những điều này liên quan đến vĩ độ và kinh độ ban đầu?

Biên tập:

Đây là một ví dụ python đơn giản cho tôi những gì tôi đang tìm kiếm:

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

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

Bị đánh cắp từ: tolatlong.py

gdal_translate sẽ định hướng lại dữ liệu của bạn từ bất kỳ phép chiếu nào sang bất kỳ hình thức nào khác (trong trường hợp này bạn muốn EPSG: 4326) bằng cách sử dụng:

gdal_translate -a_srs epsg:4326 srcfile new_file 

hoặc bạn có thể sử dụng gdaltrasform để chuyển đổi các điểm (và tôi chắc chắn bạn cũng có thể truy cập từ Python (?))

Geotransform được ghi lại tại https://gdal.org/user/raster_data_model.html . Ý tưởng là bạn lấy tọa độ (x, y) trực tiếp từ tập dữ liệu, áp dụng phép biến đổi tuyến tính để lấy (u, v) với

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

(bạn có thể coi đây là định nghĩa của phép biến đổi tuyến tính), sau đó chuyển nguồn gốc bằng cách thêm geotransform [0] vào u và geotransform [3] thành v. Nó thực sự có ý định xoay, thay đổi tỷ lệ, có lẽ sửa một chút cho một số lỗi sai và đặt lại tọa độ cụ thể của dữ liệu (x, y) để khớp với hệ tọa độ đã biết. Kết quả được cho là tạo ra tọa độ dự kiến. Điều này đơn giản có nghĩa là có một thủ tục toán học lấy (kinh độ, vĩ độ) và biến chúng thành tọa độ đã biết: đây được gọi là "phép chiếu". "Không dự đoán" là làm ngược lại; vì vậy, nếu bạn biết phép chiếu nào là cần thiết, bạn áp dụng điều đó cho tọa độ biến đổi affine (x, y) để có được vĩ độ và kinh độ.

Nhân tiện, các giá trị của các hằng số a, b, c, d được cho bởi các mục 1, 2, 4 và 5 trong mảng geotransform.

Bạn có thể sử dụng như sau:

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)

tọa độ sẽ trả về kinh độ (x) và vĩ độ (y) của tất cả các pixel. Hãy nhớ rằng tọa độ là góc bên trái của pixel