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});