GEE区域统计(求各省的中值并做进一步处理)

更新时间:2023-07-17 23:32:33 阅读: 评论:0

GEE区域统计(求各省的中值并做进⼀步处理)
⽬标:对中国的CH4(TROPOMI)提取异常值
做法:对每个省的CH4求中值,然后⽤该省的CH4影像减该中值,每个省都这样做,最终获得⼀张全国的影像。
最帅的代码:
var china = ee.FeatureCollection("urs/peisipand/cn_province");
心理课的收获与感悟var image = ee.ImageCollection('COPERNICUS/S5P/OFFL/L3_CH4')
意大利薄饼.lect('CH4_column_volume_mixing_ratio_dry_air')
.filterDate('2019-03-01','2021-04-14').mean().rename('CH4').clip(china);
var band_viz ={
min:0,
max:30,
palette:['black','blue','purple','cyan','green','yellow','red']
};
是造句
var median_value = duceRegions({
reducer: dian(),
collection: china,席珍
scale:3000,
tileScale:4
}).filter(Null(['median']));
var province_anomaly = median_value.map(function(province){
var median_value = ee.Image(ee.('median')));赞叹
return ee.Image(image.clip(province).subtract(median_value));
相逢是首歌
});
//这⼀步⾮常关键,理论上上⼀步map完之后返回的就应该是ImageCollection类型,但是不加这个强制转换就会报错
var province_anomaly = ee.ImageCollection(province_anomaly)
print(province_anomaly)//可以打印出来看看状况
小学英语质量分析// Map.addLayer(an(),band_viz, 'China CH4')
Drive({
image: an(),
哈尔滨医科大学研究生院
description:'CH4_anomaly_all_provinces',
folder:'CH4/four_provinces',
region: china,
scale:3000,
crs:'EPSG:4326',
maxPixels:1e15});
⼀般帅的代码:
思路reduceRegions完得到featureCollection之后,reduceToImage()将⽮量转为栅格2,像素值取median属性,然后栅格1减栅格2,但是获得的栅格2分辨率是11千⽶,⽽栅格1的分辨率是3km,在省份的边缘可能存在⼀些问题
这两段代码均可运⾏,并能加载到图层,但是导出图像时,⼜会遇到错误。
Error creating or submitting task Request payload size exceeds the limit: 4194304
解决⽅法有以下⼏种:
1.先导出到asts⾥,再从asts中下载图像
3.也是最好⽤的⽅法:ry().bounds()
province_anomaly = an().clip(china)
Drive({
image: province_anomaly,
description:'CH4_anomaly_all_provinces',
folder:'CH4/four_provinces',
region: ry().bounds(),
scale:3000,
crs:'EPSG:4326',
maxPixels:1e13});

本文发布于:2023-07-17 23:32:33,感谢您对本站的认可!

本文链接:https://www.wtabcd.cn/fanwen/fan/82/1102107.html

版权声明:本站内容均来自互联网,仅供演示用,请勿用于商业和其他非法用途。如果侵犯了您的权益请与我们联系,我们将在24小时内删除。

标签:栅格   代码   收获
相关文章
留言与评论(共有 0 条评论)
   
验证码:
推荐文章
排行榜
Copyright ©2019-2022 Comsenz Inc.Powered by © 专利检索| 网站地图