Red Edge Position Enmap

In "Fundamental of geological and environmental remote sensing R.K.Vincent, Prentice Hall (1997)" cap 7 viene indicato che lo spostamento del Red Edge della vegetazione puo' essere un indicatore diagnostico per il tipo di mineralogia del suolo quando questo non possa essere visibile a causa della copertura vegetale

 

 

pag 197 del libro 

Questo spostamento e' indicativamente verso il blu per zone a maggiore mineralizzazione (anche se puo' accadere l'opposto)

Gavorrano

Ho provato con lo script sottostante a verificare cosa accadeva nei dintorni dei depositi minerari della Toscana Meridionale

 Questi sono gli spettri di alcuni punti nell'intorno

se si effettua la normalizzazione fondamentalmente non si vede nessun spostamento del red edg

 

 

import h5py

import numpy as np

import rasterio

from rasterio.transform import from_origin

file_path = 'gavorrano.h5'

output_tiff = 'mappa_red_edge_enmap.tif'

def get_band(f, band_number):

band_name = f"bands/band_{band_number:03d}"

return f[band_name][:].astype(float)

with h5py.File(file_path, 'r') as f:

r670 = get_band(f, 42)

r700 = get_band(f, 52)

r740 = get_band(f, 62)

r780 = get_band(f, 75)

re_reflectance = (r670 + r780) / 2

rep = 700 + 40 * ((re_reflectance - r700) / (r740 - r700))

rep = np.where((rep > 680) & (rep < 760), rep, -999)

rows, cols = rep.shape

transform = from_origin(0, 0, 30, 30) # (West, North, x_size, y_size)

crs = 'EPSG:4326' # 32632 per UTM 32N

with rasterio.open(

output_tiff,

'w',

driver='GTiff',

height=rows,

width=cols,

count=1,

dtype=rep.dtype,

crs=crs,

transform=transform,

nodata=-999

) as dst:

dst.write(rep, 1)