谷歌地球引擎GEE:一景遥感影像ee.Image数据的基本处理方法

  本文主要对GEE 中的ee.Image格式数据图层基本处理操作加以介绍。

  本文是谷歌地球引擎(Google Earth Engine,GEE)系列教学文章的第十篇,更多GEE 文章请参考专栏:GEE学习与应用。

  首先,导入# 谷歌地球引擎GEE的JavaScript代码基础知识与用法中提及的遥感影像。

js 复制代码
var landsat_5=ee.Image("LANDSAT/LT05/C01/T1_TOA/LT05_123032_20100605");

  这里需要注意:在# 谷歌地球引擎GEE的JavaScript代码基础知识与用法中也提到了,通过上述代码导入的这一景遥感影像是成像时间为2010年06月05日、Path与Row号分别为123与032(覆盖北京市)的Landsat 5 Collection 1 Tier 1 的大气表观反射率TOA Reflectance 产品;其成像时间是一个确定的值,因此是一景图像 ,此即ee.Image格式的数据。而若是导入很多景图像------例如将上述代码括号中的内容改为"LANDSAT/LT05/C01/T1_TOA",那么所导入的遥感影像就成了不限制成像时间的、全球的Landsat 5 Collection 1 Tier 1 大气表观反射率TOA Reflectance 产品,即从Landsat 5 卫星从升空成像一直到退役这个时间范围内的影像全部包括在内;因此在全球任意一个空间范围内,所导入的遥感影像都是很多时相所对应的很多景图像 ,相当于是好几景图像重叠、叠合在一起,此即ee.ImageCollection格式的数据。关于ee.ImageCollection格式数据的基本处理操作,我们将会在后期的博客中进行介绍。

  在# 谷歌地球引擎GEE:遥感影像可视化修改的代码方法中,我们介绍了在Map.addLayer()函数中对图层可视化参数进行配置的方法。

js 复制代码
Map.setCenter(116.36863, 39.961029,10);
Map.addLayer(landsat_5,{bands:["B4","B3","B2"]},"LANDSAT5");

  而上述代码有一个问题------我们必须按照一定顺序填写Map.addLayer()函数的参数,不能跳过。例如,如果在Map.addLayer()函数中,我们不想配置{bands:["B4","B3","B2"]}这个可视化参数,但还需要配置"LANDSAT5"这个地图中所显示图层名称的参数,那么就需要将可视化参数部分用{}留空。即以上代码就需要改为:

js 复制代码
Map.addLayer(landsat_5,{},"LANDSAT5");

  针对这种情况,我们可以引入参数名称,将Map.addLayer()函数的参数修改为字典的形式,这样就可以消除函数参数顺序的问题。不过这里要注意:不要忘记字典前后需要有大括号{}包围哦。

js 复制代码
Map.addLayer({visParams:{bands:["B4","B3","B2"]},eeObject:landsat_5,name:"LANDSAT5_New"});

  其中,visParams、eeObject、name就是Map.addLayer()函数的参数名称。函数的参数名称可以在编辑栏编写代码时,按下GEE 快捷键Ctrl+Space实现获取(有点类似于代码提示的功能);而在Windows 系统中,由于这一对快捷键同时是Windows 输入法中切换中英文输入的快捷键,因此由于快捷键冲突导致GEE 中无快捷键Ctrl+Space实现获取参数名称。针对这一问题,我个人比较习惯直接搜索对应的函数,在GEE官方文档中查看相关函数的详细解释,如下图所示。

  执行代码,可以看到这种基于字典形式的参数输入方式与# 谷歌地球引擎GEE:遥感影像可视化修改的代码方法中采用的参数输入方式效果是一致的。

  此外,# 谷歌地球引擎GEE:遥感影像可视化修改的代码方法中还介绍了将全部可视化参数首先存储于一个独立的参数中,随后在Map.addLayer()函数中调用这一独立参数的方法;这一点本文就不再赘述。

  在# 谷歌地球引擎GEE的JavaScript代码基础知识与用法中,我们通过Map.setCenter()函数,将交互式地图的显示位置自动定位到北京师范大学附近的位置;与此同时,Landsat 5 遥感影像的覆盖范围非常广泛,而上述自动定位的范围则比较小,因此可以通过ROI对遥感影像显示区域加以限制。

js 复制代码
Map.setCenter(116.36863, 39.961029,13);
//Map.addLayer(landsat_5,{bands:["B4","B3","B2"]},"LANDSAT5");
var roi=ee.Geometry.Point([116.36863, 39.961029]).buffer(2000);
Map.addLayer(roi);
Map.addLayer({visParams:{bands:["B4","B3","B2"]},eeObject:landsat_5.clip(roi),name:"LANDSAT5_New"});

  其中,首先基于ee.Geometry.Point()函数,通过一个给定的经、纬度,确定一个点类型的矢量要素,随后通过.buffer()函数构建缓冲区,作为ROI ;并将ROI 与经过.clip()函数裁剪后的遥感影像分别在地图中显示。其中,.buffer()函数的详细介绍请参考# 谷歌地球引擎GEE:绘制遥感影像的直方图、时间走势曲线图;.clip()函数的详细介绍请参考# 谷歌地球引擎GEE:基于矢量裁剪栅格图像。

相关推荐
西安小哥7 分钟前
端侧跑大模型,真正的门槛不在参数量:一次 596 秒到 10 秒的压缩路径复盘
前端·人工智能
陈随易28 分钟前
Cloudflare 出钱,Vite 团队挑人,44 个维护者分 30 万美元
前端·后端·程序员
恋猫de小郭1 小时前
Flutter THLC 的 WebView 滚动终于从 70 FPS 提升到 120 FPS
android·前端·flutter
严肃小白兔2 小时前
TanStack Router 更新构建工具链:Node.js 兼容性、路由分割与前端工程实践
前端·node.js·状态模式
IT_陈寒2 小时前
Vue这个奇葩的响应式问题让我加班到凌晨
前端·人工智能·后端
风骏时光牛马2 小时前
功能可用性自动化测试验证
前端
LucianaiB12 小时前
【元宝搜索开发者征文大赛】AI时代,学得慢就不用学?用元宝搜索打造一个教程保质期,确保每一个教程符合你的现状
前端·人工智能·ai编程
excel12 小时前
JavaScript 对象池:减少高频对象创建与 GC 压力
前端
陈随易14 小时前
纯 Go 写的 JS 引擎,跑分超过了 C 语言的 QuickJS
前端·后端·程序员
子兮曰16 小时前
Chrome 权限元素实战:144 与 151 的渐进增强写法
前端·javascript·chrome