微信公众号搜"智元新知"关注
微信扫一扫可直接关注哦!

Google Earth Engine:波段值直方图不相等

如何解决Google Earth Engine:波段值直方图不相等

我正在Google Earth Engine中对Landsat-8图像进行云遮罩,然后创建不同年份的EVI像素值的直方图。对于我创建的每个层,我使用相同的几何体裁剪输出,但是直方图中的像素数彼此之间有很大差异。下面是其中一些直方图的图片,这些直方图已导出到excel,并显示为折线图,以便于查看。

Screenshot of my graph values

下面是我的代码摘录,整个内容由重复多次但相同年份的同一代码组成,而不是每个部分中的“ 2018”。我已经尝试过使用各种可用的Landsat-8集合,并根据需要编辑了代码,但是它们看起来都不正确,或者我无法使其正常工作(我对GEE还是很陌生)。我还确定我的金字塔政策,规模和区域在整个过程中都是相同的。

为什么会这样,我该如何解决?我了解到,由于使用了云遮罩,每个输出中的总像素会略有变化,但是在这种规模上不会产生任何效果,因为只有很小的区域被遮罩了。

var forest = broADLEAF.merge(MIXED).filterBounds(geometry4);
////////////////CLOUD MASK////////////////
function maskL8sr(image) {
  var cloudShadowBitMask = 1 << 3;
  var cloudsBitMask = 1 << 5;
  var qa = image.select('pixel_qa');
  var mask = qa.bitwiseAnd(cloudShadowBitMask).eq(0)
      .and(qa.bitwiseAnd(cloudsBitMask).eq(0));

  return image.updateMask(mask).divide(10000)
      .select("B[0-9]*")
      .copyProperties(image,["system:time_start"]);
}
//////////////2018/////////////////////////
var collection18 = ee.ImageCollection('LANDSAT/LC08/C01/T1_SR')
    .filterDate('2018-06-01','2018-08-31')
    .filterBounds(geometry4)
    .sort('CLOUD_COVER',false)
    .map(maskL8sr);

var composite18 = collection18.median();

// display the results.
Map.addLayer(composite18,{bands: ['B4','B3','B2'],min: 0,max: 0.3});
//////////////////////////////////////////////////////////
var evi2018 = composite18.expression(
  '2.5 * ((NIR - RED) / (NIR + 6 * RED - 7.5 * BLUE + 1))',{
    'NIR': composite18.select('B5'),'RED': composite18.select('B4'),'BLUE': composite18.select('B2')
});

var evi2018real = evi2018.clip(forest);

var histogram5 = ui.Chart.image.histogram({
  image: evi2018real,region: geometry4,scale: 30,maxPixels: 1e13})
    .setSeriesNames(['18']);

print(histogram5);
Map.addLayer(evi2018real,{bands: ['constant'],max: 0.75},"2018");

版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。

相关推荐


Selenium Web驱动程序和Java。元素在(x,y)点处不可单击。其他元素将获得点击?
Python-如何使用点“。” 访问字典成员?
Java 字符串是不可变的。到底是什么意思?
Java中的“ final”关键字如何工作?(我仍然可以修改对象。)
“loop:”在Java代码中。这是什么,为什么要编译?
java.lang.ClassNotFoundException:sun.jdbc.odbc.JdbcOdbcDriver发生异常。为什么?
这是用Java进行XML解析的最佳库。
Java的PriorityQueue的内置迭代器不会以任何特定顺序遍历数据结构。为什么?
如何在Java中聆听按键时移动图像。
Java“Program to an interface”。这是什么意思?