Dos áreas diferentes para la misma imagen en el motor de Google Earth

Question

IMPORTACIONES:

var polygon = /* color: #d63000 */ee.Geometry.Polygon(
        [[[-119.33364769821117, 46.05123532178373],
          [-119.3233620672313, 45.869732769408905],
          [-119.04111088542663, 45.873079023065166],
          [-119.0396574679861, 46.045448840018565]]]),
    landsat8 = ee.ImageCollection("LANDSAT/LC8_L1T_TOA_FMASK");

PRIMER CÓDIGO: En el primer código, se filtra la colección y se aplica el mosaico a pesar de tratarse de una sola imagen.

var PIXEL_SCALE = 30; // Meters. Resolution of most Landsat 8 bands
var PIXEL_AREA = PIXEL_SCALE * PIXEL_SCALE; // Square meters.

// Fmask classification values
var FMASK_CLEAR_GROUND = 0;
var FMASK_WATER = 1;
var FMASK_CLOUD_SHADOW = 2;
var FMASK_SNOW = 3;
var FMASK_CLOUD = 4;

var mosaic = landsat8
  .filterBounds(polygon)
  .filterDate('2016-08-01', '2016-08-30')
  .mosaic();

// Update the mask on our mosaic to mask cloud and cloud shadow pixels
var fmask = mosaic.select('fmask');
var cloudMask = fmask.neq(FMASK_CLOUD).and(fmask.neq(FMASK_CLOUD_SHADOW));
var maskedMosaic = mosaic.updateMask(cloudMask);

Map.addLayer(fmask, {min:0, max:4, palette:'green, blue, black, cyan, 
white'}, 'Fmask');
Map.addLayer(maskedMosaic.select('B4'), {min:0, max:0.5, palette:'yellow, 
green'}, 'Masked NIR');
Map.setCenter(-119.34, 45.97, 8);

// Calculate the number of pixels of each classification in our polygon
var regionCoverHistogram = mosaic.select('fmask')
 .reduceRegion(ee.Reducer.frequencyHistogram(), polygon, PIXEL_SCALE);
print('Fmask class pixel count within region', regionCoverHistogram);

var waterPixelCount = ee.Dictionary(regionCoverHistogram.get('fmask'))
  .get(FMASK_WATER.toString());

var waterArea = ee.Number(waterPixelCount).multiply(PIXEL_AREA);
print('Water Area (sq meters) in region', waterArea);

SEGUNDO CÓDIGO:

Mientras que en el segundo código, estoy tomando directamente el ID de imagen de la misma imagen y ejecutando el mismo algoritmo.

var PIXEL_SCALE = 30; // Meters. Resolution of most Landsat 8 bands
var PIXEL_AREA = PIXEL_SCALE * PIXEL_SCALE; // Square meters.

// Fmask classification values
var FMASK_CLEAR_GROUND = 0;
var FMASK_WATER = 1;
var FMASK_CLOUD_SHADOW = 2;
var FMASK_SNOW = 3;
var FMASK_CLOUD = 4;

var mosaic  = ee.Image('LANDSAT/LC8_L1T_TOA_FMASK/LC80440282016227LGN00');

// Update the mask on our mosaic to mask cloud and cloud shadow pixels
var fmask = mosaic.select('fmask');
var cloudMask = fmask.neq(FMASK_CLOUD).and(fmask.neq(FMASK_CLOUD_SHADOW));
var maskedMosaic = mosaic.updateMask(cloudMask);

Map.addLayer(fmask, {min:0, max:4, palette:'green, blue, black, cyan, 
white'}, 'Fmask');
Map.addLayer(maskedMosaic.select('B4'), {min:0, max:0.5, palette:'yellow, 
green'}, 'Masked NIR');
Map.setCenter(-119.34, 45.97, 8);

//
// Calculating Region Cover Statistics
//

// Calculate the number of pixels of each classification in our polygon
var regionCoverHistogram = mosaic.select('fmask')
 .reduceRegion(ee.Reducer.frequencyHistogram(), polygon, PIXEL_SCALE);
print('Fmask class pixel count within region', regionCoverHistogram);

var waterPixelCount = ee.Dictionary(regionCoverHistogram.get('fmask'))
.get(FMASK_WATER.toString());

var waterArea = ee.Number(waterPixelCount).multiply(PIXEL_AREA);
print('Water Area (sq meters) in region', waterArea);

Obtengo dos áreas de agua diferentes para la misma imagen. ¿Por qué?

Answer 1

Está obteniendo respuestas diferentes porque el primer código está utilizando la proyección por defecto del mosaico (WGS84) y el segundo está utilizando la proyección por defecto de la imagen (EPSG:32628). La reproyección está alterando significativamente la banda QA.

Siempre debe especificar un crs a reduceRegion para saber en qué proyección se está realizando el cálculo.