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