谷歌地球使用教程-三维地图动画，谷歌地球模型

本文的主要内容如下。

G

中间

投影信息

和

参考坐标系

和

空间变换参数

请他们给你介绍一下。

本文是Google Earth Engine (GEE) 教育文章系列的第12 篇文章。

G

文章请参见GEE学习与应用专栏（GEE学习与应用_疯狂GIS学习博客-CSDN博客）。

在前11 篇文章中

G

我们还在我们的教学博客上提供详细信息。

G

各种代码规则和具体操作，但不涉及地理的重要部分——

投影

；这是因为我们是

G

系统内部进行各种地理操作时，会根据输入输出数据的属性自动调整投影信息，无需用户自行调整对应的地理坐标。系统、投影坐标系等，非常有用。但是，在某些情况下您可能需要指定，例如：

G

您可能需要自己执行投影操作，例如当选择特定图层的投影坐标系时，或者导入图像中的不同波段具有不同的投影信息时。本文介绍

G

如何获取图层投影信息、参考坐标系、投影坐标转换参数。

G

中重投影的介绍和操作将在下一篇博客中详细讲解。

首先，2023年08月3日，按照GEE教育博客第9期（Google Earth EngineGoogle地球引擎GEE_疯狂GIS学习博客-CSDN博客的JavaScript脚本语言代码基本规则及函数语句）中提到的导入遥感图像的方法导入图像。月亮，和

小路

123号

， 排

032号

Landsat 8 Collection 1 Tier 1 的大气视反射率TOA 反射率产品（覆盖北京）。根据该场景遥感影像的中心经纬度对地图进行缩放，并将遥感影像显示在地图上。

var landsat_5=ee.Image(\’LANDSAT/LC08/C01/T1_TOA/LC08_123032_20230803\’);Map.centerObject(landsat_5,8);Map.addLayer(landsat_5);

Map.centerObject()

该函数表示给定地理特征根据其中心经度和纬度的缩放比例。它的两个参数是作为参考的地理特征（在本文中是刚刚导入的场景的遥感图像）和比例因子。关于GEE中缩放因子的详细解释，请参见GEE教育博客第9部分（Google Earth EngineGoogle地球引擎GEE_疯狂GIS学习博客-CSDN博客中JavaScript脚本语言代码的基本规则和函数语句）。接下来我们看一下我们之前导入的场景。

陆地卫星8号

以遥感图像为例，我们将获取并打印其投影信息。

print(\’该图像：的投影和变换信息\’, landsat_5.projection());投影（）

用于获取图像投影信息的函数。但是，当我运行上面的代码时，出现以下错误。

您可以看到，您无法使用.projection()，因为Landsat 8 Collection 1 Tier 1 TOA 反射率产品图像中每个波段之间的投影信息不一致。

此函数检索该场景图像的投影信息。针对这种情况，首先打印该场景的图像并检查其波段信息。

print(landsat_5); 运行代码以获取该遥感图像的波段信息。

当您打印结果时，您将看到该场景的图像12

每个波段的投影坐标系一致，均为EPSG:32650

，即WGS 84

UTM区50N以下

坐标系。每个波段之间的投影坐标系似乎匹配，但为什么我仍然收到上述错误？

如果你看GEE官方文档，你会发现不仅每个波段之间的投影坐标系必须一致，而且每个波段对应图层的空间分辨率也必须一致。只有这样才能进行投影。信息一致，可以使用.projection()。

此函数检索该场景图像的投影信息。我们再看一下print()

该功能打印的遥感影像信息在其第7波段（即B8）。

）这肯定与其他乐队相矛盾。

展开波段信息将更清楚地显示第7 个波段（即B8）。

）以及——与其他频段的区别第7频段与其他频段的crs不一样。

一致，crs_transform

我们与其他乐队不同。其中，crs

（即坐标参考系，CRS）表示该波段的投影信息的基准参考坐标系，crs_transform

表示波段投影坐标系与基准参考坐标系之间的变换参数。由此看来，参考坐标系crs

是一致的，但是空间分辨率不一致，所以它的变换参数crs_transform

这种不匹配不可避免地会产生第七个带（即B8）。

）与其他波段的投影信息不同，导致上述误差。

因此，有必要对频带进行滤波。首先按照第二篇GEE教育博客（Google Earth EngineGoogle地球引擎GEE_疯狂GIS学习博客-CSDN博客JavaScript脚本语言代码的基本规则和函数语句）中描述的搜索GEE数据的方法来搜索波段信息。显示。 TOA Landsat 8 Collection 1 Tier 1 大气表观反射率反射率产品。

可以看到是B2。

,B3

含B4

这些带对应三种颜色：蓝色、绿色和红色。继续使用这三个频段作为示例。它基于GEE教育博客第6期（Google Earth Engine GEE栅格代数与手动计算NDVI波段计算_疯狂GIS学习博客-CSDN博客_Gee波段计算）中介绍的.select()。

在您的函数中，采用上面的三个带并重用.projection()。

获取其投影信息的函数。

var Band=landsat_5.select(\’B[2-4]\’);print(\’此图像：的投影和变换信息\’,band.projection()); 当我运行代码时，成功检索到投影信息。我明白这一点。获得。

您还可以通过使用.projection() 并在右侧打印成功的投影信息来验证。

函数得到的投影信息包括crs

带转换

——这两个部分再次证明，在使用.projection()之前，除了投影坐标系的基础参考系之外，变换参数必须在各个频段保持一致。

功能。接下来还可以获取当前图层投影信息的线性比例（单位：m）。这里获得的比率是标称像素大小。这在GEE官方手册中称为“当前层金字塔底层的标称像素大小”，层金字塔的底层实际上是层本身。这里得到的实际上是遥感图像本身的空间分辨率。

print(\’Pixel size inmeter:\’,band.projection().nominalScale()); 可以看到，得到的结果是30。

m，这是Landsat 8 在可见光波段的空间分辨率。