Diferenciación del tipo de hielo mediante datos de satélite

Question

Me gustaría clasificar la zona glaciar por categorías de nieve (si la hay) y hielo, pero lo más importante: entre hielo viejo y hielo fresco. Tienen diferentes propiedades que son posibles de reconocer en el campo, pero ¿se puede hacer esto con datos de satélite? (preferiblemente Landsat por su resolución espacial de 30/15m)

Answer 1

Para ello, tendrás que utilizar datos de microondas. Los datos ópticos no son suficientes. Si quieres seguir con los datos ópticos, dime qué metodología has seguido. También depende en gran medida de la topografía y de la vegetación de la zona. La clasificación de los datos de microondas no es en sí misma sencilla, tendrás que consultar mucha literatura y elegir la metodología que más te convenga. Por favor, vea la metodología que seguí en mi tesis de maestría: http://www.iirs.gov.in/iirs/sites/default/files/StudentThesis/Sanjay_MTech_2013-15.pdf

Por favor, pregunte si tiene alguna duda después de revisar la literatura.

Answer 2

Aquí tienes un ejemplo que te puede ayudar a iniciarte en Google-Earth-Engine utilizando la banda C de 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);

Tendrá que clasificar la imagen utilizando los algoritmos de clasificación supervisada mencionados aquí: https://developers.google.com/earth-engine/classification

Más información sobre el uso de Sentinel-1 https://developers.google.com/earth-engine/sentinel1

Sobre el motor de Google Earth y los glaciares: http://www.geo.uzh.ch/~mzemp/share/scratch/msc/MSc.Thesis_NoahZeltner_UsingGoogleEarthEngineForGlobalGlacierChangeAssessment.pdf

Sobre el SAR y las zonas glaciares: http://www.sciencedirect.com/science/article/pii/S0034425713001703