Differenziazione del tipo di ghiaccio usando i dati satellitari

Vorrei classificare l'area glaciale per categorie di neve (se presente) e ghiaccio, ma ciò che è più importante: tra ghiaccio vecchio e ghiaccio fresco. Hanno proprietà diverse che è possibile riconoscere sul campo, ma puoi farlo con i dati satellitari? (preferibilmente Landsat a causa di una risoluzione spaziale di 30 / 15m)

— adamczi
fonte

Quali caratteristiche hanno il ghiaccio vecchio e nuovo nel campo?

— Aaron

1) La neve fresca è molto meno compattata rispetto al vecchio ghiaccio del ghiacciaio (diventa ghiaccio per compattazione). Pertanto, ciò può essere in qualche modo collegato alla riflettanza IR, che viene assorbita dall'acqua. 2) Anche la neve fresca ha un'albedo anche fino a quasi il 100%, ma la neve vecchia può arrivare fino a ~ 40% (nessuna classificazione rigorosa ovviamente). Vorrei fare uso di IR, poiché la composizione True Color non è così utile come mi piacerebbe che fosse.

— Adamczi,

Sembra un semplice problema di classificazione delle immagini. È necessario iniziare con i dati di allenamento, che possono essere raccolti sul campo o selezionando abilmente i pixel dalle immagini.

— Aaron

Penso che la strada da percorrere qui sia un algoritmo di classificazione supervisionato come Maximum Likelihood, Random Forests, ecc. Che utilizza tutte le bande spettrali disponibili. Conosci questi metodi? Non sono sicuro di cosa intendi per "composizione IR". Ti riferisci alla creazione di immagini composite come composito a falsi colori (ovvero NIR, R, G)? In tal caso, sei molto limitato nelle applicazioni di tali prodotti.

— Aaron

@adamczi prova a lavorare su google-earth-engine. saranno disponibili algoritmi di classificazione supervisionata e dati SAR (il tuo upload o il cloud di Google).

— csheth,

Risposte:

Per questo dovrai usare i dati delle microonde. I dati ottici non li tagliano. Se vuoi ancora passare all'ottica, dimmi quale metodologia hai seguito. Molto dipende anche dalla topografia, LULC della tua zona. La classificazione dei dati a microonde non è di per sé semplice, dovrai consultare molta letteratura e scegliere una metodologia adatta a te. Si prega di consultare la metodologia che ho seguito nella mia tesi di laurea specialistica: http://www.iirs.gov.in/iirs/sites/default/files/StudentThesis/Sanjay_MTech_2013-15.pdf

Si prega di chiedere se avete qualche domanda dopo aver esaminato la letteratura.

— coprire con una volta
fonte

Ecco un esempio che può aiutarti a iniziare su Google-Earth-Engine utilizzando la banda C di Sentinel-1:

var pt = ee.Geometry.Point(96.7868, 29.31409);

// Filter collection around point. Also read up on Sentinel-1's 
// polarization 
var collection = ee.ImageCollection('COPERNICUS/S1_GRD').filterBounds(pt)
.filter(ee.Filter.listContains('transmitterReceiverPolarisation', 'VV'))
.select('VV');

// select an appropriate date
var beforesnow = collection.filterDate('2016-11-01', '2016-12-01').mosaic();
var aftersnow = collection.filterDate('2017-02-01', '2017-03-01').mosaic();

// bands for Sentinel-2
var bands = ['B2', 'B3', 'B4'];

// Some Sentinel-2 images for reference
var S2 = ee.ImageCollection('COPERNICUS/S2').filterBounds(pt)
.select(bands);
var S2before = S2.filterDate('2016-10-01', '2016-11-30').mosaic();
var S2after = S2.filterDate('2017-01-01', '2017-02-01').mosaic();

Map.addLayer(S2before, {bands: ['B4', 'B3', 'B2'], min: 300,max: 5000}, 'S2 Before');
Map.addLayer(S2after, {bands: ['B4', 'B3', 'B2'], min:873,max: 12522}, 'S2 After');

Map.centerObject(pt, 13);

// you may change the min, max later when tinkering with the layers tab in // the map
Map.addLayer(beforesnow, {min:-30,max:0}, 'Before snow');
Map.addLayer(aftersnow, {min:-30,max:0}, 'After snow');

//Some information on the Sentinel-1 collection
print('Collection: ', collection);

Dovrai classificare l'immagine usando algoritmi di classificazione supervisionati menzionati qui: https://developers.google.com/earth-engine/classification

Altre informazioni sull'uso di Sentinel-1 https://developers.google.com/earth-engine/sentinel1

Su Google Earth Engine and Glaciers: http://www.geo.uzh.ch/~mzemp/share/scratch/msc/MSc.Thesis_NoahZeltner_UsingGoogleEarthEngineForGlobalGlacierChangeAssessment.pdf

Su SAR e zone del ghiacciaio: http://www.sciencedirect.com/science/article/pii/S0034425713001703

— csheth
fonte