Google Earth Engine(GEE)——sentinel-5P 臭氧O3数据逐日数据分析

还是关于之前逐日数据的问题,那天我没有改数据就直接用了年度合成,因为我担心这些数据肯定不是逐日更新的结果果然是这样,所以我今天在有的同学的追问下再次进行了尝试。她的要求是让进行2020年12月4日到8日的数据。代码和结果如下

代码:

//矢量边界
var hh= ee.FeatureCollection("users/bqt2000204051/CHINA")
//这里把之前的年改称了天数,
var days = ee.List.sequence(4, 8);
days.evaluate(function(r) 
{
  r.map(function(r) {
  var start = ee.Date.fromYMD(2020, 11, r);
  var stop =  ee.Date.fromYMD(2020, 12, r);
  var scol= ee.ImageCollection("COPERNICUS/S5P/OFFL/L3_O3")
          .filterBounds(hh)
          .filterDate(start,stop)
          .select("O3_column_number_density")

          
;
print("scol",scol)
  var ndvi_before =scol.qualityMosaic("O3_column_number_density").clip(hh);
  ndvi_before=ndvi_before.set({name:ee.String(ee.Number(r).int())})
  Export.image.toDrive({
    image: ndvi_before.select("SO3_column_number_density"),
    region:hh,
    scale:1000,
 
  • 0
    点赞
  • 10
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 4
    评论
以下是使用GEEGoogle Earth Engine)进行Sentinel-1数据反演土壤湿度的代码示例: 1.导入Sentinel-1数据: ``` var sentinel1 = ee.ImageCollection('COPERNICUS/S1_GRD') .filterBounds(geometry) .filterDate('2018-01-01', '2021-12-31') .filter(ee.Filter.listContains('transmitterReceiverPolarisation', 'VV')) .filter(ee.Filter.listContains('transmitterReceiverPolarisation', 'VH')) .filter(ee.Filter.eq('instrumentMode', 'IW')) .filter(ee.Filter.eq('orbitProperties_pass', 'DESCENDING')) .select(['VV', 'VH']) ``` 2.计算土壤湿度: ``` var s1 = ee.Image(sentinel1.first()).clip(geometry) var vh = s1.select('VH') var vv = s1.select('VV') var b = 0.7 //地表粗糙度系数 var alpha = 1 //极化反照率系数 var eps = 0.015 //介电常数 //计算表面散射系数 var gamma_naught_vh = ee.Image.constant(10).multiply(vh).log10().multiply(10).subtract(83) var gamma_naught_vv = ee.Image.constant(10).multiply(vv).log10().multiply(10).subtract(83) //计算比例系数 var rho = gamma_naught_vh.divide(gamma_naught_vv) //计算土壤湿度 var soil_moisture = ee.Image.constant(1.5).multiply(rho.divide(Math.sin(Math.toRadians(40))).add(1)).log10().multiply(-10/b).divide(alpha*eps) ``` 3.可视化土壤湿度: ``` var vizParams = { min: -1, max: 1, palette: ['blue', 'white', 'green'] } Map.addLayer(soil_moisture, vizParams, 'soil moisture') ``` 注意:以上代码仅为示例,具体使用时需要根据实际情况进行修改。同时,需要注意在使用GEE时要了解相关的API和限制条件。

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 4
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

此星光明

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

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

余额充值