当前位置: 首页 > wzjs >正文

微信表情包制作网站响应式网页设计原理

微信表情包制作网站,响应式网页设计原理,flash布局网站,余姚网站建设服务目录 1 数据加载与预处理 2 叶绿素含量反演 3 冠层水分反演 4 数据可视化与导出 5 完整代码 6 运行结果 在生态学和环境科学领域,植被的健康状况是评估生态系统稳定性和功能的关键指标之一。而叶绿素含量和冠层水分含量作为反映植被生理状态的重要参数&#x…

目录

1 数据加载与预处理

2 叶绿素含量反演

3 冠层水分反演

4 数据可视化与导出

5 完整代码

6 运行结果


在生态学和环境科学领域,植被的健康状况是评估生态系统稳定性和功能的关键指标之一。而叶绿素含量和冠层水分含量作为反映植被生理状态的重要参数,一直是遥感监测的核心目标。接下来,将通过Google Earth Engine(GEE)平台,分享如何利用Sentinel-2卫星数据反演叶绿素含量和冠层水分含量,并通过代码实现这一过程。

1 数据加载与预处理

在进行生态参数反演之前,需要加载并预处理卫星数据。这里选择了Sentinel-2卫星的“COPERNICUS/S2_SR_HARMONIZED”数据集,它提供了高质量的地表反射率产品,具有10米至60米的空间分辨率,非常适合植被监测。

/******************** 数据加载与预处理 ********************/
var s2 = ee.ImageCollection("COPERNICUS/S2_SR_HARMONIZED")// 设置时间范围.filterDate('2023-01-01', '2023-11-30')// 限定研究区域.filterBounds(roi)// 筛选云量小于20%的影像.filter(ee.Filter.lt('CLOUDY_PIXEL_PERCENTAGE', 20))// 波段名称标准化.map(function (img) {return img.rename(['B1', 'B2', 'B3', 'B4', 'B5', 'B6', 'B7','B8', 'B8A', 'B9', 'B11', 'B12','AOT', 'WVP', 'SCL', 'TCI_R', 'TCI_G', 'TCI_B','MSK_CLDPRB', 'MSK_SNWPRB', 'QA10', 'QA20', 'QA60','MSK_CLASSI_OPAQUE', 'MSK_CLASSI_CIRRUS', 'MSK_CLASSI_SNOW_ICE']);})// 计算中值合成影像.median()// 裁剪至研究区域.clip(roi);

关键步骤解析:

  1. 时间范围筛选:通过filterDate方法,选择了2023年1月1日至11月30日的数据。您可以根据研究需求调整时间范围。

  2. 云量筛选:通过filter方法结合ee.Filter.lt,剔除了云量超过20%的影像,以减少云层对分析结果的干扰。

  3. 波段名称标准化:Sentinel-2数据的波段名称通常以“B04”、“B05”等形式出现。为了方便后续计算,将波段名称统一更正为“B4”、“B5”等。

  4. 中值合成:通过median方法,计算了中值合成影像,以减少噪声和异常值的影响。

  5. 裁剪至研究区域:使用clip方法将影像裁剪至感兴趣区域(roi),确保分析的精确性。

云掩膜处理:

云层是遥感影像中常见的干扰因素,它会掩盖地表的真实信息。Sentinel-2数据提供了“SCL”(Scene Classification Layer)波段,用于地表分类。其中,值为3的像素表示云中像素,值为8的像素表示云影,值为9的像素表示薄云。可以通过以下代码构建云掩膜:

// 云掩膜(基于SCL分类波段)
var cloudMask = s2.select('SCL').neq(3) // 排除云中像素.and(s2.select('SCL').neq(8))      // 排除云影.and(s2.select('SCL').neq(9));     // 排除薄云var s2Clean = s2.updateMask(cloudMask); // 应用云掩膜

通过neqand操作符,构建了一个无云的掩膜,并将其应用到中值合成影像上,得到了干净无云的影像数据(s2Clean)。

2 叶绿素含量反演

叶绿素是植物进行光合作用的核心色素,其含量直接反映了植物的生长状况和健康水平。在遥感领域,叶绿素含量可以通过植被反射光谱的特定波段计算得到,采用以下公式进行叶绿素含量反演:

其中:

  • B8A:865nm波段(近红外增强植被反射)

  • B5:705nm波段

  • B4:665nm波段

公式 3.2 * (B8A - B5)/(B8A + B5 + 0.15) * sqrt(B8A/B4) 是多个经典指数的合成优化,其中系数可以通过地面实测数据校准获得,这里取3.2。

/******************** 叶绿素含量反演 ********************/
var Cab = s2Clean.expression('3.2 * (B8A - B5) / (B8A + B5 + 0.15) * sqrt(B8A/B4)',{'B8A': s2Clean.select('B8A'),'B5': s2Clean.select('B5'),'B4': s2Clean.select('B4')}
).rename('Chlorophyll');

通过expression方法,将公式嵌入代码中,并将结果命名为“Chlorophyll”(叶绿素含量)。 

3 冠层水分反演

冠层水分含量是植被生理状态的重要指标之一,它反映了植被在水分胁迫下的适应能力和健康状况。采用以下公式进行冠层水分反演:

/******************** 冠层水分反演 ********************/
var Cw = s2Clean.expression('(B11 - B12)/(B11 + B12 + 0.03) * log(B8A)',{'B11': s2Clean.select('B11'),'B12': s2Clean.select('B12'),'B8A': s2Clean.select('B8A')}
).rename('CanopyWater');

通过expression方法,我们将公式嵌入代码中,并将结果命名为“CanopyWater”(冠层水分含量)。

4 数据可视化与导出

完成生态参数反演后,需要将结果可视化并导出。在GEE中,通过Map.addLayer方法将结果添加到地图上,并通过设置合适的可视化参数(如最小值、最大值和调色板),使结果更加直观易懂。

/******************** 可视化与导出 ********************/
var visCab = {min: 0,max: 6,palette: ['#FFFFCC', '#A1DAB4', '#41B6C4', '#225EA8', '#081D58']
};
var visCw = {min: -0.2,max: 0.7,palette: ['#0000FF', '#00FFFF', '#FFFF00', '#FF0000']
};Map.addLayer(Cab, visCab, '叶绿素含量(μg/cm²)');
Map.addLayer(Cw, visCw, '冠层水分指数');// 分项导出设置
Export.image.toDrive({image: Cab,description: 'Cab_20m_S2_2023Summer',folder: 'Sentinel2_EcoParams',region: roi,scale: 20,maxPixels: 1e13
});
Export.image.toDrive({image: Cw,description: 'Cw_20m_S2_2023Summer',folder: 'Sentinel2_EcoParams',region: roi,scale: 20,maxPixels: 1e13
});

关键步骤解析:

  1. 可视化参数设置:通过minmaxpalette参数,为叶绿素含量和冠层水分含量设置了合适的可视化范围和颜色渐变。

  2. 导出设置:通过Export.image.toDrive方法,将反演结果导出至Google Drive。导出的图像分辨率为20米,与Sentinel-2影像的空间分辨率相匹配。

5 完整代码

Map.centerObject(roi,7)
var styling = {color:"red",fillColor:"00000000"};
Map.addLayer(roi.style(styling),{},"geometry")/******************** 数据加载与预处理 ********************/
var s2 = ee.ImageCollection("COPERNICUS/S2_SR_HARMONIZED").filterDate('2023-01-01', '2023-11-30') // 需修改日期可调整此处.filterBounds(roi).filter(ee.Filter.lt('CLOUDY_PIXEL_PERCENTAGE', 20)).map(function (img) {// 波段名称标准化(用户数据实际波段名)return img.rename(['B1', 'B2', 'B3', 'B4', 'B5', 'B6', 'B7','B8', 'B8A', 'B9', 'B11', 'B12', // 关键修正点:B04→B4,B05→B5'AOT', 'WVP', 'SCL', 'TCI_R', 'TCI_G', 'TCI_B','MSK_CLDPRB', 'MSK_SNWPRB', 'QA10', 'QA20', 'QA60','MSK_CLASSI_OPAQUE', 'MSK_CLASSI_CIRRUS', 'MSK_CLASSI_SNOW_ICE']);}).median().clip(roi);// 云掩膜(基于SCL分类波段)
var cloudMask = s2.select('SCL').neq(3) // 3=云中像素.and(s2.select('SCL').neq(8)) // 8=云影.and(s2.select('SCL').neq(9)); // 9=薄云var s2Clean = s2.updateMask(cloudMask);/******************** 叶绿素含量反演 ********************/
var Cab = s2Clean.expression('3.2 * (B8A - B5) / (B8A + B5 + 0.15) * sqrt(B8A/B4)',{'B8A': s2Clean.select('B8A'), // 865nm(原B8A)'B5': s2Clean.select('B5'), // 705nm(原B05)'B4': s2Clean.select('B4') // 665nm(原B04)}
).rename('Chlorophyll');/******************** 冠层水分反演 ********************/
var Cw = s2Clean.expression('(B11 - B12)/(B11 + B12 + 0.03) * log(B8A)',{'B11': s2Clean.select('B11'), // 1610nm'B12': s2Clean.select('B12'), // 2190nm'B8A': s2Clean.select('B8A') // 近红外增强植被反射}
).rename('CanopyWater');/******************** 可视化与导出 ********************/
var visCab = {min: 0,max: 6,palette: ['#FFFFCC', '#A1DAB4', '#41B6C4', '#225EA8', '#081D58']
};
var visCw = {min: -0.2,max: 0.7,palette: ['#0000FF', '#00FFFF', '#FFFF00', '#FF0000']
};Map.addLayer(Cab, visCab, '叶绿素含量(μg/cm²)');
Map.addLayer(Cw, visCw, '冠层水分指数');// 分项导出设置
Export.image.toDrive({image: Cab,description: 'Cab_20m_S2_2023Summer',folder: 'Sentinel2_EcoParams',region: roi,scale: 20,maxPixels: 1e13
});Export.image.toDrive({image: Cw,description: 'Cw_20m_S2_2023Summer',folder: 'Sentinel2_EcoParams',region: roi,scale: 20,maxPixels: 1e13
});

6 运行结果

点击RUN即可下载数据
叶绿素含量反演结果可视化
冠层水分反演结果可视化

文章转载自:

http://F2p0K0S5.hqmfn.cn
http://jZwtuxbb.hqmfn.cn
http://HFqeuqTp.hqmfn.cn
http://SlCEbt6R.hqmfn.cn
http://xGagCOB3.hqmfn.cn
http://rbagJvQY.hqmfn.cn
http://Ry0VW5qD.hqmfn.cn
http://NKUOJSeB.hqmfn.cn
http://ICWg7pjw.hqmfn.cn
http://cZRZ7tY9.hqmfn.cn
http://xLLG6Zqi.hqmfn.cn
http://whLwO5zY.hqmfn.cn
http://d3NluJ3D.hqmfn.cn
http://iniA0Fz2.hqmfn.cn
http://TwTcPnHn.hqmfn.cn
http://rVkgrN1Y.hqmfn.cn
http://UdGdRMl5.hqmfn.cn
http://na3062oG.hqmfn.cn
http://R1KmTpvf.hqmfn.cn
http://hUUPZAPb.hqmfn.cn
http://7AbLKJm8.hqmfn.cn
http://kpXDzr7Q.hqmfn.cn
http://Ipy1sKiY.hqmfn.cn
http://uxvISdgP.hqmfn.cn
http://rBIMCYDZ.hqmfn.cn
http://JzjTpivo.hqmfn.cn
http://eWT30X4W.hqmfn.cn
http://y7ShsijI.hqmfn.cn
http://848umcvw.hqmfn.cn
http://Zrtmzds6.hqmfn.cn
http://www.dtcms.com/wzjs/715908.html

相关文章:

  • 熟人做网站怎么收钱网站建立的步骤是( )。
  • 贵阳建网站WordPress输出当前网址
  • 承德网站建设步骤深圳关键词优化报价
  • 网站后台登陆网址是多少鹤岗住房和城乡建设局网站
  • 网站页脚导航带积分的网站建设
  • 广州免费建站哪里有网站建设宣传 mp4
  • 网站建设备案审核要多久十大免费行情软件入口下载
  • 电子政务 和网站建设总结wordpress删除菜单
  • 重庆优化网站全网热搜榜
  • 湖南建设部网站网站制作方案垂直领域获客
  • 一家做特卖的网站手机版wordpress 内存使用教程
  • 唐山网站建设报价网站做app安全吗
  • 大型网站流量上海进一步优化
  • 电子商务网站的特色wordpress设置登陆口
  • 做标签网站是干嘛的高校服务地方专题网站建设
  • 在百度怎么建自己的网站网站的类型有哪几种
  • 大学生网站制作作业免费下载东营招标信息网
  • 企业建站的费用广州网页设计美工培训
  • 宁波网站备案幕布拍照手机网站 jquery 特效
  • 泉州网站建设qzdzi负面信息搜索引擎 网站
  • 赣州网站建设设计pc端网站怎么做自适应
  • 化工外贸网站建设中山网直播
  • 外贸品牌网站设计东莞网站建设实例分析
  • 个人商城网站能备案吗上传网站内容
  • 建设母婴网站的目的商务型网站建设
  • 天猫网站设计分析做动画网站去哪采集
  • 免费做房产网站有哪些用软件做的网站权限管理
  • 锦州做网站的个人江阴公司企业网站建设
  • 有哪些程序网站河北邢台路桥建设公司网站
  • 什么网站做推广wordpress主题失败