Google Earth Engine谷歌地球引擎投影、参考坐标系及转换参数获取( 二 )


Google Earth Engine谷歌地球引擎投影、参考坐标系及转换参数获取

文章插图
 
  因此,我们需要对波段进行筛选 。首先,依据第二篇GEE教学博客(Google Earth Engine谷歌地球引擎GEE中JavaScript脚本语言代码基础规则与函数语句_疯狂学习GIS的博客-CSDN博客)中提及的GEE数据搜索方法,我们搜索并查看Landsat 8 Collection 1 Tier 1的大气表观反射率TOA Reflectance产品的波段信息 。
Google Earth Engine谷歌地球引擎投影、参考坐标系及转换参数获取

文章插图
 
  可以看到,其中B2、B3与B4波段分别为蓝、绿、红三种颜色对应的波段;我们就以这三个波段为例继续进行后续操作 。基于第六篇GEE教学博客(Google Earth Engine谷歌地球引擎GEE栅格代数与NDVI波段计算手动求取_疯狂学习GIS的博客-CSDN博客_gee波段计算)中介绍的.select()函数,将上述三个波段取出,并重新使用.projection()函数获取其投影信息 。
var band=landsat_5.select("B[2-4]");print("Projection and transformation information of this image:",band.projection());  执行代码,可以看到投影信息已经成功得到 。
Google Earth Engine谷歌地球引擎投影、参考坐标系及转换参数获取

文章插图
 
  通过右侧打印成功的投影信息我们还可以注意到,使用.projection()函数获取的投影信息包括crs与transform两个部分——这也再一次证明了波段之间除了投影坐标系的基准参考坐标系需要一致,还需要转换参数一致,才可以使用.projection()函数 。
  随后,我们还可以将当前图层投影信息的线性比例(单位为m)加以获取 。这里获取的比例即Nominal Pixel Size,在GEE官方手册中将其称为“当前图层的金字塔中,最底层的名义像素大小”;而一个图层的最底层金字塔其实就是该图层自身,因此这里求得的其实就是遥感影像自身的空间分辨率 。
print("Pixel size in meters:",band.projection().nominalScale());
Google Earth Engine谷歌地球引擎投影、参考坐标系及转换参数获取

文章插图
 
  可以看到,获得的结果为30m,亦即Landsat 8可见光波段的空间分辨率 。

【Google Earth Engine谷歌地球引擎投影、参考坐标系及转换参数获取】


推荐阅读