GEE栅格插值(以TROPOMI的NO2数据为例)

GEE插值


GEE提供了IDW和Kriging两种插值

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

// Define an area to perform interpolation over.
var yunnan = shp.filter(ee.Filter.eq('OBJECTID', 25));

Map.addLayer(ee.Image().paint(yunnan, 0, 2), {}, 'yunnanborder');


var NO2 = ee.ImageCollection('COPERNICUS/S5P/OFFL/L3_NO2').select('tropospheric_NO2_column_number_density').filterDate('2019-01-06','2019-01-08').mean().rename('NO2');

// Sample the methane composite to generate a FeatureCollection.
var samples = NO2.addBands(ee.Image.pixelLonLat())
  .sample({region: yunnan, numPixels: 1500,
    scale:1000, projection: 'EPSG:4326'})
  .map(function(sample) {
    var lat = sample.get('latitude');
    var lon = sample.get('longitude');
    var NO2 = sample.get('NO2');
    return ee.Feature(ee.Geometry.Point([lon, lat]), {NO2: NO2});
  });
  
Map.addLayer(samples);


  // get mean
var mean_value = samples.reduceColumns({
reducer: 'mean',
selectors: ['NO2']
});

// get stdev
var sd_value = samples.reduceColumns({
reducer: 'stdDev',
selectors: ['NO2']
});

print(mean_value,sd_value)

// Do the interpolation, valid to 70 kilometers.
var interpolated = samples.inverseDistance({
  range: 7e4,
  propertyName: 'NO2',
  mean: mean_value.get('mean'),
  stdDev: sd_value.get('stdDev'),
  gamma: 0.3});

  
// var band_viz = {
//   min: 0,
//   max: 0.0002,
//   palette: ['black', 'blue', 'purple', 'cyan', 'green', 'yellow', 'red']
// };

// Map.centerObject(yunnan, 7);
// Map.addLayer(NO2, band_viz, 'NO2');
// Map.addLayer(interpolated, band_viz, 'NO2 Interpolated');


  Export.image.toDrive({
      image: interpolated.clip(yunnan),
      description: 'IDW-test',
      crs: 'EPSG:4326',
      scale: 10000,
      region: yunnan,
      maxPixels: 1e13,
      folder: 'IDW'
    })


结果如下
插值前:
在这里插入图片描
插值后:

在这里插入图片描述

参考的是GEE官方提供的文档

评论 6
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值