Coregistrazione immagini satellitari con AROSIC

Nonostante le immagini satellitare siano gia' georiferite per fare analisi multitemporale a livello pixel serve che le immagini si sovrappongano a livello subpixel

Per Enmap si puo' usare Arosics  https://github.com/GFZ/arosics 

Si installa tramite 

pip install arosics

in seguito all'interno di QGis  si installa il plugin Coregistratio. Il plugin non crea nessuna icona e deve richiamato tramite Processing Toolbox

Si prendono i due geotiff caratterizzata da Spectral_Image caricandoli in QGis e visualizzzando in entrambi a scala di grigio la sola banda 55 (o comunque una banda IR che sia sufficiente luminosa e definita)

A questo punto si puo' lanciare Automated Global Co-Registration per una prima analisi

Si devono inserire l'immagine di reference e quella di target e si lancia

 

Algorithm started at: 2026-11-01T20:57:09

Algorithm 'Automated global Co-Registration' starting…

Input parameters:

{ 'ALIGN_GRIDS' : True, 'IMG_REF' : '/home/linnocenti/Downloads/dims_op_oc_oc-en_704147776_1/ENMAP.HSI.L2A/ENMAP01-____L2A-DT0000087498_20240814T104352Z_033_V010506_20261028T222244Z/reference.TIF', 'INPUT' : '/home/linnocenti/Downloads/dims_op_oc_oc-en_704147776_1/ENMAP.HSI.L2A/ENMAP01-____L2A-DT0000099945_20241103T104352Z_033_V010506_20261028T221806Z/target.TIF', 'MATCHING_WINDOW_CENTER' : None, 'MATCHING_WINDOW_SIZE' : 256, 'MATCH_GSD' : True, 'MAX_SHIFT' : 5, 'OUTPUT' : 'TEMPORARY_OUTPUT', 'RESAMPLING' : 2 }

Image to image Co-Registration:

Processing file: /home/linnocenti/Downloads/dims_op_oc_oc-en_704147776_1/ENMAP.HSI.L2A/ENMAP01-____L2A-DT0000099945_20241103T104352Z_033_V010506_20261028T221806Z/target.TIF

Perform automatic subpixel co-registration with AROSICS...

Calculating footprint polygon and actual data corner coordinates for reference image...

Bounding box of calculated footprint for reference image:

(666825.0, 4851015.0, 702795.0, 4886055.0)

Calculating footprint polygon and actual data corner coordinates for image to be shifted...

Bounding box of calculated footprint for image to be shifted:

(666105.0, 4847745.0, 702045.0, 4882845.0)

Matching window position (X,Y): 684449.5912947934/4866898.006505515

Detected integer shifts (X/Y): 1/1

Detected subpixel shifts (X/Y): -0.4115450978279114/0.25640174746513367

Calculated total shifts in fft pixel units (X/Y): 0.5884549021720886/1.2564017474651337

Calculated total shifts in reference pixel units (X/Y): 0.5884549021720886/1.2564017474651337

Calculated total shifts in target pixel units (X/Y): 0.5884549021720886/1.2564017474651337

Calculated map shifts (X,Y): 17.65364706516266/-37.69205242395401

Calculated absolute shift vector length in map units: 41.62141360683733

Calculated angle of shift vector in degrees from North: 334.9032417536094

Original map info: ['UTM', 1.0, 1.0, 666075.0, 4882875.0, 30.0, 30.0, 32, 'North', 'WGS-84']

Updated map info: ['UTM', 1.0, 1.0, '666092.6536470652', '4882837.307947576', 30.0, 30.0, 32, 'North', 'WGS-84']

Image similarity within the matching window (SSIM before/after correction): 0.7379 => 0.7687

Estimated reliability of the calculated shifts: 84.9%

Correcting geometric shifts...

Writing GeoArray of size (1171, 1199, 224) to /tmp/processing_nJnJJo/5c2dac6e43944e51904b71996b75fb54/OUTPUT.tif.

 

L'aumento del parametro SSIM da 0.738 to 0.769 indica che la traslazione migliora la coregistrazione. Reliability non e' ottimale ma comunque accettabile. Lo shift proposto e' +17.65 m E, −37.69 m N quindi superiore ad 1 px in N

La coregistrazione e' ottimale solo nella Custom matching windows size e oiu' e' grande migliore e' il risultato. Cliccando si puo' scegliere dove e' centrata la finestra, altrimenti viene scelta al centro dell'immagine