Виступ у підсумковій конференції
Нікоряка Віктора Васильовича
Код Python для обчислення NDWI
from osgeo import gdal
import numpy as np
def calculate_ndwi(green_band_path, nir_band_path, output_path):
"""
Обчислює NDWI (Normalized Difference Water Index) і зберігає результат у GeoTIFF.
Параметри:
----------
green_band_path : str
Шлях до файлу з зеленим каналом (Green).
nir_band_path : str
Шлях до файлу з ближнім інфрачервоним каналом (NIR).
output_path : str
Шлях для збереження результату у форматі GeoTIFF.
Формула:
--------
NDWI = (Green - NIR) / (Green + NIR)
Приклад використання:
---------------------
calculate_ndwi("green_band.tiff", "nir_band.tiff", "output_ndwi.tiff")
"""
try:
# Відкриття файлів каналів
green_ds = gdal.Open(green_band_path)
nir_ds = gdal.Open(nir_band_path)
if green_ds is None or nir_ds is None:
raise ValueError("Не вдалося відкрити один із вказаних файлів.")
# Зчитування даних з каналів
green_band = green_ds.GetRasterBand(1).ReadAsArray().astype(np.float32)
nir_band = nir_ds.GetRasterBand(1).ReadAsArray().astype(np.float32)
# Розрахунок NDWI
ndwi = (green_band - nir_band) / (green_band + nir_band)
ndwi[np.isinf(ndwi)] = np.nan # Уникнення помилок через ділення на нуль
# Створення GeoTIFF для результату
driver = gdal.GetDriverByName('GTiff')
output_ds = driver.Create(
output_path,
green_ds.RasterXSize,
green_ds.RasterYSize,
1,
gdal.GDT_Float32
)
output_ds.SetGeoTransform(green_ds.GetGeoTransform())
output_ds.SetProjection(green_ds.GetProjection())
output_ds.GetRasterBand(1).WriteArray(ndwi)
output_ds.FlushCache()
print(f"NDWI успішно збережено у: {output_path}")
except Exception as e:
print(f"Помилка під час обчислення NDWI: {e}")
finally:
# Закриття файлів
green_ds = None
nir_ds = None
output_ds = None
# Приклад виклику функції
calculate_ndwi("path_to_green_band.tiff", "path_to_nir_band.tiff", "output_ndwi.tiff")##Карта землекористування LULC
Графік змін землекористування
Порівняння середніх значень типів землекористування
