使用Google Earth Engine制作作物时间序列变化曲线
2019-01-22 14:54
2651 查看
//Sentinel2数据制作时间序列NDVI,没有进行云量筛选,有待完善 //定义图表中不同作物类型的颜色 var COLOR = { DADOU: '7FFF00', GAOLIANG: 'ff0000', YUMI: '0000ff' }; //定义感兴趣区 var gaoliang = ee.Feature( ee.Geometry.Rectangle(123.269, 48.0495, 123.272, 48.0511), {label: 'Gaoliang'}); var yumi = ee.Feature( ee.Geometry.Rectangle(123.2892, 48.0599, 123.2907, 48.0582), {label: 'Yumi'}); var dadou = ee.Feature( ee.Geometry.Rectangle(123.2442, 48.0813, 123.2455, 48.0799), {label: 'Dadou'}); //将三种作物感兴趣区合并 var cropRegions = new ee.FeatureCollection([gaoliang,yumi,dadou]); //筛选s2数据 var s2= ee.ImageCollection('COPERNICUS/S2') .filterDate('2018-01-01', '2019-01-01') .filterBounds(cropRegions); //计算每副影像的NDVI并制作数据集 var ndvi = s2.map(function(image) { return image.select().addBands(image.normalizedDifference(['B8', 'B4']).select([0], ['NDVI'])); }); print(ndvi) //渲染NDVI显示颜色 var vis = {min: -0.2, max: 1, palette: [ 'FFFFFF', 'CE7E45', 'FCD163', '66A000', '207401', '056201', '004C00', '023B01', '012E01', '011301' ]}; Map.addLayer(ndvi, vis, 'NDVI'); Map.addLayer(cropRegions, {color: COLOR.GAOLIANG},'ROI'); // Map.addLayer(gaoliang, {color: COLOR.GAOLIANG}); // Map.addLayer(yumi, {color: COLOR.YUMI}); // Map.addLayer(dadou, {color: COLOR.DADOU}); //定义图表及样式 var ndviTimeSeries = ui.Chart.image.seriesByRegion({ imageCollection: ndvi, regions: cropRegions, reducer: ee.Reducer.mean(), band: 'NDVI', scale: 10, xProperty: 'system:time_start', seriesProperty: 'label' }); ndviTimeSeries.setChartType('ScatterChart'); ndviTimeSeries.setOptions({ title: 'Sentinel-2数据作物时间序列NDVI变化', vAxis: { title: 'NDVI' }, lineWidth: 1, pointSize: 4, series: { 0: {color: COLOR.GAOLIANG}, 1: {color: COLOR.YUMI}, 2: {color: COLOR.DADOU} } }); print(ndviTimeSeries);
相关文章推荐
- Facebook开源时间序列内存数据库Beringei,追求极致压缩率——如果是int根据大多数时间序列中的值与相邻数据点相比并没有显著的变化,只要使用XOR将当前值与先前值进行比较,然后存储发生变化的比特。最终,该算法将整个数据集至少压缩了90%
- Sentinel 2 Atmospheric Correction in Google Earth Engine(使用GEE进行哨兵2数据大气校正)
- 国产“google earth”制作依靠卫星影像图下的“虚拟地球仪” 时间:2009-05-13 19:38来源:互联网
- 使用Google Earth Engine(一):提取下载MODIS、Landsat点数据
- 使用 Google Earth + GPSBabel 生成 GPS NMEA 文件
- 使用 Eclipse 在 Google App Engine 上创建 mashup之第 1 部分:创建应用程序
- 使用RMAN的不完全恢复-基于时间/SCN/日志序列
- Google App Engine 使用指南~
- 使用 Eclipse + PyDev 开发 Google App Engine 程序
- 在C++中嵌入JavaScript——Google V8 JavaScript Engine使用体验
- 在.NET平台使用Google Earth COM API示例(01)
- 使用 Eclipse 在 Google App Engine 上创建 mashup
- 使用GAE(Google App Engine)一步一步开发自己的博客[1]:新手上路
- 使用jquery制作背景变化菜单
- Google App Engine 使用指南
- 使用 amcharts 和 highcharts 绘制多曲线时间趋势图的通用方法
- Android之使用AchartEngineActivity引擎绘制柱状图、曲线图
- 在.NET平台使用Google Earth COM API示例(05)
- Android之使用AchartEngineActivity引擎绘制柱状图、曲线图
- [转]使用Google App Engine Helper for Django