GEE中影像镶嵌的mosaic函数的应用

本文通过实验验证了Google Earth Engine (GEE) 中影像镶嵌时选取影像值的原则。实验使用了两个常量影像,值分别为1和0,结果显示镶嵌后的影像值为0,即列表中后一个影像的值。

 GEE中影像镶嵌一般会取覆盖在最上边的影像值作为镶嵌影像的值,那么那个才是最上边得影像呢?

var imageCollection=ee.ImageCollection([image1,image2])
var image3=imageCollection.mosaic()

那么image3是取image1的值还是取image2的值呢,为此,我做了一个实验测试,创建了南宁市的常量影像image1和image2,其影像值分别为1和0,测试镶嵌结果image3的值是1还是0,代码如下:

//创建常量影像
var image1=ee.Image(1).clip(nanning.first().geometry()).toInt()
var image2=ee.Image(0).clip(nanning.first().geometry()).toInt()
var imageCollection=ee.ImageCollection([image1,image2])
print(imageCollection)
//镶嵌
var image3=imageCollection.mosaic()
//加载,inspector查看影像值
Map.addLayer(image1,{},'image1')
Map.addLayer(image2,{},'image2')
Map.addLayer(image3,{},'image3')
Map.centerObject(nanning.first().geometry())

最终发现镶嵌结果影像的值为0,也就是列表中索引较大位置的影像值,因而实际镶嵌过程中是取最新的影像值</

我的研究区是阿拉尔市,使用GEE代码预处理遥感地图影像,遇到的问题是导出10m的研究区影像是一块块的,以下是我的代码,我出现了什么问题?该如何修改// 1. 加载矢量边界 var region = table; // 可视化边界 Map.centerObject(region, 10); Map.addLayer(region, {}, 'Region Boundary'); // 2. 加载Sentinel-2影像集合 var collection = ee.ImageCollection("COPERNICUS/S2_SR_HARMONIZED") .filterBounds(region) // 过滤边界范围 .filterDate('2022-01-01', '2022-05-31') // 时间范围 .filter(ee.Filter.lt('CLOUDY_PIXEL_PERCENTAGE', 10)); // 云量过滤 print(collection) // 3. 影像合成:取中位数影像 var image = collection.median().clip(region); // 4. 显示所有波段 var visParamsAllBands = { min: 0, max: 3000, gamma: 1.4 }; Map.addLayer(image, visParamsAllBands, 'All Bands'); // 5. 单独显示部分常用波段组合 var visParamsRGB = {bands: ['B4', 'B3', 'B2'], min: 0, max: 3000}; Map.addLayer(image, visParamsRGB, 'True Color (RGB)'); var visParamsFalseColor = {bands: ['B8', 'B4', 'B3'], min: 0, max: 3000}; Map.addLayer(image, visParamsFalseColor, 'False Color (NIR)'); // 6. 计算NDVI(归一化植被指数) var ndvi = image.normalizedDifference(['B8', 'B4']).rename('NDVI'); var ndviParams = {min: -1, max: 1, palette: ['blue', 'white', 'green']}; Map.addLayer(ndvi, ndviParams, 'NDVI'); // 7. 导出所有波段影像到Google云盘 Export.image.toDrive({ image: image, description: 'Sentinel2_AllBands_2022', scale: 10, region: region.geometry().bounds(), crs: 'EPSG:4326', // 坐标系 fileFormat: 'GeoTIFF', folder: 'Sentinel2_Exports', maxPixels:1e13 }); // 8. 导出NDVI影像到Google云盘 Export.image.toDrive({ image: ndvi, description: 'Sentinel2_NDVI_2022', scale: 10, region: region.geometry().bounds(), crs: 'EPSG:4326', fileFormat: 'GeoTIFF', folder: 'Sentinel2_Exports' });
04-01
评论 5
成就一亿技术人!
拼手气红包6.0元
还能输入1000个字符
 
红包 添加红包
表情包 插入表情
 条评论被折叠 查看
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

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

余额充值