相关文章推荐
严肃的面包  ·  google ...·  11 月前    · 
严肃的面包  ·  Google Earth ...·  11 月前    · 
严肃的面包  ·  Google Earth ...·  11 月前    · 
严肃的面包  ·  google earth map ...·  11 月前    · 
Map = geemap.Map()
# 加载影像数据
naip2004_2012 = ee.ImageCollection('USDA/NAIP/DOQQ') \
    .filterBounds(ee.Geometry.Point(-71.09, 42.4)) \
    .filterDate('2004-07-01', '2012-12-31') \
    .select(['R', 'G', 'B'])
# 最大值合成
composite = naip2004_2012.max()
Map.addLayer(composite, {}, 'max value composite')
Map.setCenter(-71.09, 42.4, 14)

结果如下:
在这里插入图片描述

2、影像镶嵌Mosaic

  • ee.ImageCollection.mosaic()
    在这里插入图片描述
    注意:镶嵌mosaic更主要体现的不同区域影像的合成操作
    需要同时拼接四个不同的 DOQQ,但位置不同

python代码如下

# 加载跨4幅影像的数据
naip2012 = ee.ImageCollection('USDA/NAIP/DOQQ') \
    .filterBounds(ee.Geometry.Rectangle([-71.18, 42.35, -71.09, 42.41])) \
    .filterDate('2012-01-01', '2012-12-31')
# 影像镶嵌mosaic
mosaic = naip2012.mosaic()
Map.addLayer(naip2012, {}, 'naip2012')
Map.setCenter(-71.12532, 42.3712, 12)

结果如下:
在这里插入图片描述

3、进一步挖掘Mosaic的特性

重点:mosaic()方法根据它们在集合中的顺序(最后在顶部)合成重叠图像。要控制镶嵌(或合成)中像素的来源,可使用图像掩膜mask。
例如,以下使用光谱阈值进行镶嵌新的图像数据

python代码如下:

# 加载影像,显示
naip = ee.Image('USDA/NAIP/DOQQ/m_4207148_nw_19_1_20120710')
Map.setCenter(-71.0915, 42.3443, 12)
Map.addLayer(naip, {}, 'naip')
# 生成NDVI、NDWI指数
ndvi = naip.normalizedDifference(['N', 'R']).rename('ndvi')
ndwi = naip.normalizedDifference(['G', 'N']).rename('ndwi')
# 利用上述2个指数,采用不同阈值条件,探测区域裸地分布
# 裸地1,植被指数低于0.2,同时水体指数低于0.3的区域
bare1 = ndvi.lt(0.2).And(ndwi.lt(0.3))
# 裸地2,植被指数低于0.2,同时水体指数低于0.8的区域。可见其敏感性比上一阈值较低
bare2 = ndvi.lt(0.2).And(ndwi.lt(0.8))
# 显示参数
vis_ndvi = {'min':-1, 'max':1, 'palette':['red', 'green']}
vis_ndwi = {'min':0.5, 'max':1, 'palette':['gray', 'blue']}
# 加载显示
Map.addLayer(ndvi, vis_ndvi, 'ndvi')
Map.addLayer(ndwi, vis_ndwi, 'ndwi')
Map.addLayer(bare1, {}, 'bare1')
Map.addLayer(bare2, {}, 'bare2')

显示结果:
在这里插入图片描述

# 掩膜并镶嵌可视化影像,可以看到最后的图层位于最上层(这是Mosaic的特点)
mosaic = ee.ImageCollection([
    ndwi.updateMask(ndwi.gte(0.5)).visualize(**vis_ndwi), 
    ndvi.updateMask(ndvi.gte(0.2)).visualize(**vis_ndvi),
    bare2.updateMask(bare2.And(bare1.Not())).visualize(**{'palette':['gray']}),
    bare1.updateMask(bare1).visualize(**{'palette':['white']})
]).mosaic()
Map.addLayer(mosaic, {}, 'visualization mosaic')

结果如下:
在这里插入图片描述

  1. 分析了影像合成与影像镶嵌的两种不同影响合成方法,其中影像合成composite更多的表现在同一区域不同时间影像的组合,多使用的合成方法有最大化合成(MVC,max value composite)、中值合成(median)、均值合成(mean)等合成方法,在具体工作中可根据实际需要采用不同的合成方法;镶嵌(mosaic)更多的表现在不同区域上的多幅影像的合成。
  2. 镶嵌的特点是最后的图层位于最上层,这也是Mosaic的特点之一。
  • https://developers.google.com/earth-engine/apidocs/ee-image-visualize?hl=en
  • https://github.com/giswqs/earthengine-py-notebooks/blob/master/ImageCollection/mosaicking.ipynb
  • https://developers.google.com/earth-engine/guides/ic_composite_mosaic?hl=en
如何在gee上绘制影像研究区影像过境的条带) 当对大型场景(研究区超过单景影像范围时)进行拍摄时,单张照片或像可能无法完全捕捉到整个场景。这时候就需要采用影像拼接技术,将多张照片或像拼接在一起,以形成一个更大的、完整的像。 在进行影像拼接时,需要将每张照片或像的像素位置映射到一个共同的坐标系。这就需要使用到行列号。行列号是用来表示每个像素位置的坐标系,其“行”代表垂直方向上的位置,而“列”代表水平方向上的位置。通过将每个像素的行列号映射到一个共同的坐标系,就可以将多张照片或像拼接在一起,以形成一个更大的、完整的像。 总之,行列号是影像处理非常重要的概念,它可以帮助我们描述和操作的每个像素位置,而影像拼接则是一种常用的技术,可以将多张照片或像拼接在一起,以形成一个更大的、完整的像。     此次的内容是遥感镶嵌。     在遥感像处理,我们有时为了获取更大范围的地面像,通常需要将多幅遥感像拼成一幅像,这就需要使用镶嵌遥感影像进行拼接操作。     我们需要进行镶嵌的遥感卫星影像数据有两种情况,南北方向和东西方方向的,所以我们在选择数据时尽量选择成像条件相似(同一轨道、同一时间)的像。     通过使用PCI软件镶嵌模块手动对长沙地区影像(南... var imageCollection=ee.ImageCollection([image1,image2]) var image3=imageCollection.mosaic() 那么image3是取image1的值还是取image2的值呢,为此,我做了一个实验测试,创建了南宁市的常量影像image1和image2,其影像值分别为1和0,测试镶嵌结果image3的值是1还是0,代码如下: //创建...
本篇主要介绍如何加载单幅Sentinel2影像,以及合成长时间序列感兴趣区域的影像集合,自己在写代码的过程,出现了一些问题,比如导出的影像会以多幅tif影像出现,导出的影像在Arcgis加载全为黑色,也看了一些博主的博客,有的解释很简单,在这我会详细解释一下,让大家少走弯路(被我走完了)。 遇到问题的同学可以直接跳转到第三部分。 目录1.单幅影像合成2.特定长时间序列影像合成3.导出为多幅影像以及导出影像为黑色的问题 1.单幅影像合成 此代码只对感兴趣区域,并辅以时间范围进行影像筛选,使用了经典的Sen
遥感生态指数(Remote Sensing Ecological Index)是一种基于遥感技术对生态系统进行监测和评估的指数,它通过对地表反射光谱数据的分析,反映出生态系统的植被覆盖程度、植被生长状况、土地利用类型、土地覆盖变化等生态环境信息。 遥感生态指数通常采用遥感影像数据来计算,包括多光谱影像、高光谱影像、雷达影像等。常用的遥感生态指数包括植被指数(Vegetation Index)、归一化植被指数(Normalized Difference Vegetation Index,NDVI)、土地覆盖指数(Land Cover Index)等。 其,植被指数可以用来反映植被的生长状况,归一化植被指数则是常用的植被指数,可以通过计算近红外波段和红光波段反射率之间的比值来计算;土地覆盖指数则可以用来反映不同土地利用类型的空间分布情况,如城市、耕地、林地、草地等。
merge简单说就是合并,即数据集合的合并mosaic即拼接,这两个命令常用作拼接不完整的影像时进行使用。 可以简单的说:b4=b1.merge(b2).merge(b3).mosaic() b4返回的是b3的值,如果把merge命令b2和b3位置进行调换时,mosaic()b4返回的就是b2的值了。 所以在写代码时,弄懂每个简单的命令很重要。 参考代码如下所示: https://code.earthengine.google.com/2a4fed8cb7c1b5fac9cf1ebb04f7669a
以下是一些关于 GEE 遥感影像处理的教程视频课程,供您参考: 1. 《Google Earth Engine基础入门教程》:https://www.bilibili.com/video/BV1xJ411J7tG 2. 《GEE遥感影像处理实战》:https://www.bilibili.com/video/BV1kE411T7r6 3. 《Google Earth Engine遥感应用》:https://www.bilibili.com/video/BV1cE411P7k5 4. 《Google Earth Engine遥感影像处理实战》:https://www.bilibili.com/video/BV1gE411V7rE 5. 《GEE遥感影像数据处理与分析》:https://www.bilibili.com/video/BV1pJ411c7mD 这些视频课程涵盖了 GEE 的基础知识、遥感影像数据的处理方法、常见的遥感影像分析应用等内容,可以帮助您更好地了解和掌握 GEE 遥感影像处理技术。