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

目标：对中国的CH4（TROPOMI）提取异常值
做法：对每个省的CH4求中值，然后用该省的CH4影像减该中值，每个省都这样做，最终获得一张全国的影像。

最帅的代码：

var china = ee.FeatureCollection("users/peisipand/cn_province");

var image = ee.ImageCollection('COPERNICUS/S5P/OFFL/L3_CH4')
  .select('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 = image.reduceRegions({
    reducer: ee.Reducer.median(),
    collection: china,
    scale: 3000,
    tileScale:4
}).filter(ee.Filter.notNull(['median']));

var province_anomaly = median_value.map(function(province) {
  var median_value = ee.Image(ee.Number(province.get('median')));
  return ee.Image(image.clip(province).subtract(median_value));
});

//这一步非常关键，理论上上一步map完之后返回的就应该是ImageCollection类型，但是不加这个强制转换就会报错
var province_anomaly = ee.ImageCollection(province_anomaly) 
print(province_anomaly)  //可以打印出来看看状况

// Map.addLayer(province_anomaly.mean(),band_viz, 'China CH4')

Map.centerObject(china, 4);

Export.image.toDrive({
  image: province_anomaly.mean(),
  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.先导出到assets里，再从assets中下载图像
2.export中的region可以自己绘制一个polygon，而不是采用边界复杂的shp
3.也是最好用的方法：region:china.geometry().bounds()

province_anomaly = province_anomaly.mean().clip(china)
Export.image.toDrive({
  image: province_anomaly,
  description: 'CH4_anomaly_all_provinces',
  folder: 'CH4/four_provinces',
  region: china.geometry().bounds(),
  scale: 3000,
  crs: 'EPSG:4326',
  maxPixels: 1e13});