GEE重投影——NICFI数据集重投影到WGS84坐标

简介

重投影

遥感影像的重投影是将一幅影像从一个地理坐标系统转换到另一个地理坐标系统的过程。在遥感应用中,重投影非常重要,因为不同的遥感影像可能来自于不同的遥感传感器或不同的地理坐标系统,而在进行数据分析和集成时,需要保证所有影像在同一地理坐标系统下。

重投影的过程涉及到将原始影像上的像素点重新映射到目标坐标系统上的像素点。这一过程需要使用到原始影像的地理信息和两个坐标系统之间的转换参数。

在重投影过程中,首先需要确定原始影像的地理参考信息,包括投影方式(如经纬度投影或UTM投影)、地理坐标系统(如WGS84坐标系统)、像素分辨率等。然后,确定目标坐标系统,包括投影方式和目标坐标参考信息。

接下来,需要确定原始坐标系统和目标坐标系统之间的转换参数,这些参数通常包括原始坐标系的投影参数、地理转换参数和尺度因子。这些参数可以通过地理信息系统(GIS)软件或专门的重投影工具来计算和提供。

一旦确定了转换参数,就可以进行重投影。这个过程涉及到将原始影像上的每个像素点的地理坐标转换为目标坐标系统下的地理坐标,并根据目标坐标系统的投影方式重新计算像素的位置。最后,根据目标坐标系统的分辨率和像素值插值方法,可以将原始影像像素的值重新映射到目标坐标系统的像素上。

重投影的结果是一个经过转换的新影像,它的地理坐标和分辨率与目标坐标系统相匹配。通过重投影,可以将来自不同传感器或不同坐标系统的遥感影像集成在一起,进行数据分析、特征提取和地理空间分析等任务。

WGS84:espg4326

WGS84(World Geodetic System 1984)是由美国国防部和国家地理空间情报局(NGA)共同制定的地理坐标系统。它被广泛用于全球定位系统(GPS)中,以及地图、航海、航空、测量和GIS(地理信息系统)等领域。WGS84基于地球的椭球体模型,以提供更准确和统一的地理位置表示。

WGS84坐标使用经度和纬度来表示一个地点的位置。经度测量地点相对于本初子午线(通常是格林尼治子午线)的东西方向,以度为单位,范围从-180°到+180°。负值表示位于本初子午线的西侧,正值表示位于本初子午线的东侧。

纬度测量地点相对于地球赤道的北南方向,同样以度为单位,范围从-90°到+90°。负值表示位于赤道的南侧,正值表示位于赤道的北侧。

WGS84坐标的原点是地球的质心,它是基于大量的测量数据和全球地理参考框架(例如大地水准面、地球形状和引力场)来定义的。

WGS84坐标系统在全球范围内都是通用的,允许不同国家和地区之间进行位置的准确比较和数据交换。它被广泛应用于GPS设备、地图制作、导航系统、测量和制图等领域。

需要注意的是,由于地球的形状和引力场的复杂性,WGS84坐标系统仍然有一定的精度误差。因此,在需要更精确的测量和定位应用中,可能需要使用更高精度的大地测量系统或区域特定的坐标系统。

结果

代码:

var nicfi = ee.ImageCollection("projects/planet-nicfi/assets/basemaps/americas"),
    geometry = 
    /* color: #d63000 */
    /* shown: false */
    ee.Geometry.Polygon(
        [[[-77.75868232024801, -4.95936558796817],
          [-77.25331122649801, -6.991610669222523],
          [-75.53944403899801, -9.817082479248054],
          [-73.07850653899801, -10.617171980297915],
          [-72.85877997649801, -9.29705977121429],
          [-73.5839357588108, -7.294011116424077],
          [-71.2987795088108, -4.168022106966874],
          [-71.2768068525608, -3.422595504484709],
          [-72.7709474775608, -1.8862060891271266],
          [-74.6605959150608, -0.3924031878546906],
          [-75.4296388838108, -0.7659185993699265],
          [-77.5170412275608, -2.8302194271787755],
          [-77.7367677900608, -4.4747643892937585]]]),
    someGeometry = 
    /* color: #bf04c2 */
    /* shown: false */
    ee.Geometry.Point([-73.24629103109629, -3.7319791101813076]),
    Iquitos = 
    /* color: #98ff00 */
    /* shown: false */
    ee.Geometry.Point([-73.25989842477257, -3.7680075489480154]),
    hyperion = 
    /* color: #0b4a8b */
    /* shown: false */
    /* displayProperties: [
      {
        "type": "rectangle"
      }
    ] */
    ee.Geometry.Polygon(
        [[[-76.48946704975148, -5.6377560469052685],
          [-76.48946704975148, -6.691790662938007],
          [-75.97534122577687, -6.691790662938007],
          [-75.97534122577687, -5.6377560469052685]]], null, false),
    cloud = 
    /* color: #ffffff */
    /* shown: false */
    ee.FeatureCollection(
        [ee.Feature(
            ee.Geometry.Point([-76.11586732529143, -5.800144605786184]),
            {
              "landcover": 1,
              "system:index": "0"
            }),
        ee.Feature(
            ee.Geometry.Point([-76.11539525650481, -5.802236688126257]),
            {
              "landcover": 1,
              "system:index": "1"
            }),
        ee.Feature(
            ee.Geometry.Point([-76.11586732529143, -5.803218683282618]),
            {
              "landcover": 1,
              "system:index": "2"
            }),
        ee.Feature(
            ee.Geometry.Point([-76.11934346817473, -5.812632941303244]),
            {
              "landcover": 1,
              "system:index": "3"
            }),
        ee.Feature(
            ee.Geometry.Point([-76.11824912689666, -5.812483509865488]),
            {
              "landcover": 1,
              "system:index": "4"
            }),
        ee.Feature(
            ee.Geometry.Point([-76.10711259506682, -5.8166675751229695]),
            {
              "landcover": 1,
              "system:index": "5"
            }),
        ee.Feature(
            ee.Geometry.Point([-76.10698384903411, -5.814746324531442]),
            {
              "landcover": 1,
              "system:index": "6"
            }),
        ee.Feature(
            ee.Geometry.Point([-76.08117026947478, -5.808576931062461]),
            {
              "landcover": 1,
              "system:index": "7"
            }),
        ee.Feature(
            ee.Geometry.Point([-76.08132047317962, -5.809900474457839]),
            {
              "landcover": 1,
              "system:index": "8"
            }),
        ee.Feature(
            ee.Geometry.Point([-76.08014030121306, -5.809089270810307]),
            {
              "landcover": 1,
              "system:index": "9"
            }),
        ee.Feature(
            ee.Geometry.Point([-76.07417792059695, -5.782639073843204]),
            {
              "landcover": 1,
              "system:index": "10"
            }),
        ee.Feature(
            ee.Geometry.Point([-76.07016533591067, -5.782703119279921]),
            {
              "landcover": 1,
              "system:index": "11"
            }),
        ee.Feature(
            ee.Geometry.Point([-76.06954306341922, -5.784389646504057]),
            {
              "landcover": 1,
              "system:index": "12"
            }),
        ee.Feature(
            ee.Geometry.Point([-76.0654875633887, -5.781657042905717]),
            {
              "landcover": 1,
              "system:index": "13"
            }),
        ee.Feature(
            ee.Geometry.Point([-76.06628149725711, -5.780333433293494]),
            {
              "landcover": 1,
              "system:index": "14"
            }),
        ee.Feature(
            ee.Geometry.Point([-76.06881350256717, -5.780034553275225]),
            {
              "landcover": 1,
              "system:index": "15"
            }),
        ee.Feature(
            ee.Geometry.Point([-76.07250422217166, -5.780290736157687]),
            {
              "landcover": 1,
              "system:index": "16"
            }),
        ee.Feature(
            ee.Geometry.Point([-76.07025116659915, -5.786972798704372]),
            {
              "landcover": 1,
              "system:index": "17"
            }),
        ee.Feature(
            ee.Geometry.Point([-76.06950014807498, -5.785393021074328]),
            {
              "landcover": 1,
              "system:index": "18"
            }),
        ee.Feature(
            ee.Geometry.Point([-76.06973618246829, -5.7897267248277835]),
            {
              "landcover": 1,
              "system:index": "19"
            }),
        ee.Feature(
            ee.Geometry.Point([-76.06645315863406, -5.791263794037819]),
            {
              "landcover": 1,
              "system:index": "20"
            }),
        ee.Feature(
            ee.Geometry.Point([-76.07177466131961, -5.793483997735426]),
            {
              "landcover": 1,
              "system:index": "21"
            }),
        ee.Feature(
            ee.Geometry.Point([-76.07362002112185, -5.791797497646509]),
            {
              "landcover": 1,
              "system:index": "22"
            }),
        ee.Feature(
            ee.Geometry.Point([-76.08857506691616, -5.721113412209441]),
            {
              "landcover": 1,
              "system:index": "23"
            }),
        ee.Feature(
            ee.Geometry.Point([-76.08698719917935, -5.720835851880167]),
            {
              "landcover": 1,
              "system:index": "24"
            }),
        ee.Feature(
            ee.Geometry.Point([-76.08610743462246, -5.721775286297117]),
            {
              "landcover": 1,
              "system:index": "25"
            }),
        ee.Feature(
            ee.Geometry.Point([-76.08788842140834, -5.725575709775938]),
            {
              "landcover": 1,
              "system:index": "26"
            }),
        ee.Feature(
            ee.Geometry.Point([-76.0762798207919, -5.725661112260444]),
            {
              "landcover": 1,
              "system:index": "27"
            }),
        ee.Feature(
            ee.Geometry.Point([-76.05264655366109, -5.76371047574705]),
            {
              "landcover": 1,
              "system:index": "28"
            }),
        ee.Feature(
            ee.Geometry.Point([-76.0544489981191, -5.76181039437343]),
            {
              "landcover": 1,
              "system:index": "29"
            }),
        ee.Feature(
            ee.Geometry.Point([-76.04882708802388, -5.758672043800184]),
            {
              "landcover": 1,
              "system:index": "30"
            }),
        ee.Feature(
            ee.Geometry.Point([-76.0508655668752, -5.756195510021284]),
            {
              "landcover": 1,
              "system:index": "31"
            }),
        ee.Feature(
            ee.Geometry.Point([-76.05367652192281, -5.756366305800796]),
            {
              "landcover": 1,
              "system:index": "32"
            }),
        ee.Feature(
            ee.Geometry.Point([-76.05906239762471, -5.753377372250828]),
            {
              "landcover": 1,
              "system:index": "33"
            }),
        ee.Feature(
            ee.Geometry.Point([-76.04554009163768, -5.685020694030229]),
            {
              "landcover": 1,
              "system:index": "34"
            }),
        ee.Feature(
            ee.Geometry.Point([-76.0450036498347, -5.681369467776653]),
            {
              "landcover": 1,
              "system:index": "35"
            }),
        ee.Feature(
            ee.Geometry.Point([-76.04787897789866, -5.681775160725996]),
            {
              "landcover": 1,
              "system:index": "36"
            }),
        ee.Feature(
            ee.Geometry.Point([-76.04762148583323, -5.6825011368693765]),
            {
              "landcover": 1,
              "system:index": "37"
            }),
        ee.Feature(
            ee.Geometry.Point([-76.04457449639231, -5.683824973362626]),
            {
              "landcover": 1,
              "system:index": "38"
            }),
        ee.Feature(
            ee.Geometry.Point([-76.04532551491648, -5.68794592876145]),
            {
              "landcover": 1,
              "system:index": "39"
            }),
        ee.Feature(
            ee.Geometry.Point([-76.10173342611019, -5.948838098929204]),
            {
              "landcover": 1,
              "system:index": "40"
            }),
        ee.Feature(
            ee.Geometry.Point([-76.0989439287347, -5.950139966650249]),
            {
              "landcover": 1,
              "system:index": "41"
            }),
        ee.Feature(
            ee.Geometry.Point([-76.10115406896297, -5.949585072588717]),
            {
              "landcover": 1,
              "system:index": "42"
            }),
        ee.Feature(
            ee.Geometry.Point([-76.10658286000911, -5.950118624581319]),
            {
              "landcover": 1,
              "system:index": "43"
            }),
        ee.Feature(
            ee.Geometry.Point([-76.10752699758235, -5.9477283076184495]),
            {
              "landcover": 1,
              "system:index": "44"
            }),
        ee.Feature(
            ee.Geometry.Point([-76.10638974096004, -5.9452312689816935]),
            {
              "landcover": 1,
              "system:index": "45"
            }),
        ee.Feature(
            ee.Geometry.Point([-76.10467312719051, -5.951697935444226]),
            {
              "landcover": 1,
              "system:index": "46"
            }),
        ee.Feature(
            ee.Geometry.Point([-76.10711930181209, -5.9512497530929345]),
            {
              "landcover": 1,
              "system:index": "47"
            }),
        ee.Feature(
            ee.Geometry.Point([-76.10638974096004, -5.9525942990503315]),
            {
              "landcover": 1,
              "system:index": "48"
            }),
        ee.Feature(
            ee.Geometry.Point([-76.10714075948421, -5.952231485385892]),
            {
              "landcover": 1,
              "system:index": "49"
            }),
        ee.Feature(
            ee.Geometry.Point([-76.10896466161434, -5.953725422466508]),
            {
              "landcover": 1,
              "system:index": "50"
            }),
        ee.Feature(
            ee.Geometry.Point([-76.12544415380184, -5.979057707656797]),
        
  • 61
    点赞
  • 41
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 0
    评论
以下是在 Google Earth Engine 中导入 Landsat 5 数据集的代码示例: ```javascript // 定义时间范围 var startDate = ee.Date('1984-01-01'); var endDate = ee.Date('2012-05-05'); // 导入Landsat 5表面反射率数据集 var collection = ee.ImageCollection('LANDSAT/LT05/C01/T1_SR') .filterDate(startDate, endDate) .filterBounds(geometry); // 选择需要的波段 var bands = ['B1', 'B2', 'B3', 'B4', 'B5', 'B6', 'B7']; // 选择云掩膜(可选) var maskClouds = function(image) { var qa = image.select('pixel_qa'); var cloud = qa.bitwiseAnd(1 << 5).neq(0); return image.updateMask(cloud.not()); }; // 应用云掩膜 var maskedCollection = collection.map(maskClouds); // 打印导入的图像集合信息 print(maskedCollection); // 选择特定图像进行进一步分析 var image = maskedCollection.first(); print(image); ``` 在上述代码中,我们首先定义了一个时间范围,然后使用 `ee.ImageCollection()` 函数导入 Landsat 5 表面反射率数据集。通过 `filterDate()` 方法可以筛选出指定时间范围内的图像。我们还可以使用 `filterBounds()` 方法来限定导入数据的地理范围,其中 `geometry` 是表示地理范围的几何图形。 接下来,我们选择需要的波段,可以根据自己的需求修改 `bands` 数组。如果需要进行云掩膜处理,可以定义 `maskClouds` 函数并应用于图像集合,通过 `map()` 方法应用云掩膜。 最后,我们可以打印导入的图像集合信息,并选择特定的图像进行进一步的分析。在示例中,我们选择了第一张图像进行打印。 请注意,在使用上述代码之前,需要在代码中指定感兴趣区域的几何图形,可以使用 `ee.Geometry()` 函数来定义几何图形。

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

此星光明

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值