У меня есть два растровых файла, один с разрешением 30 м и один с разрешением 100 м. Я хотел бы взять координаты каждого пикселя растра с разрешением 30 м, а затем извлечь значения для каждой из этих координат из растра с разрешением 100 м. Конечной целью является растр с разрешением 30 м со значениями из растра 100 м.
Я использую Python, и я пытался использовать модуль rasterio , но я пытаюсь понять, как извлечь точки широты / долготы (или строки / столбца) каждый пиксель.
Мой CRS: WGS84, EPSG 4326
Мое разрешение: res: (0.0002694945900000002, 0.00026949459)
Пока мой код:
import numpy as np
import skimage.transform as st
import rasterio
import pyproj
30m_file = rasterio.open(path_to_30m_file)
100m_file = rasterio.open(path_to_100m_file)
band1 = 30m_file.read(1)
print("band1:", band1)
#Access values based on georeferenced space
x, y = 30m_file.bounds.left + 0.003, 30m_file.bounds.top - 0.002
row, col = 30m_file.index(x, y)
print("row:", row)
print("col:", col)
print("example value:", band1[row, col])
#The logic would then follow:
# for every pixel lat/lon in 30m file
# open 100m file
# extract 100m value at those 30m pixel lat/lon
# keeping 30m resolution
Но мне мешают получить доступ к координатам и затем проанализировать их.