MODIS 质量控制图层 QC

MODIS产品质量控制文件使用方法_mctk如何进行批处理-CSDN博客icon-default.png?t=N7T8https://blog.csdn.net/tk20190411/article/details/106183163?ops_request_misc=%257B%2522request%255Fid%2522%253A%2522170195808116800211580116%2522%252C%2522scm%2522%253A%252220140713.130102334.pc%255Fblog.%2522%257D&request_id=170195808116800211580116&biz_id=0&utm_medium=distribute.pc_search_result.none-task-blog-2~blog~first_rank_ecpm_v1~rank_v31_ecpm-2-106183163-null-null.nonecase&utm_term=%E8%B4%A8%E9%87%8F%E6%8E%A7%E5%88%B6&spm=1018.2226.3001.4450

基于python和google earth engine的MODIS数据质量控制————以MOD11A1为例-CSDN博客icon-default.png?t=N7T8https://blog.csdn.net/VictoriaLy/article/details/130259754主要参考了以上2篇blog。

第一篇的代码正确,可以直接用。

第二篇的代码有1个地方需要修改,见下图。

修改原因:in_Array里的nan值,如果不处理,直接进行np.array(in_Array,dtype = np.uint16),nan会被转为0。在接下来筛选gooddata时,0(10进制)转为二进制是00,nan值被当作了gooddata,这是不对的。

另外,还要搞懂QC图层,每个像素值的含义:

python代码里,把十进制转为16位2进制之后,不足16位的是向左边补零直到满足16位。例如:11(十进制)转为二进制:1011,16位2进制就是:0000000000001011。

MCTK里转出数据时,选择double precision,导出的数据类型是float 64.但是,在提取QC掩码的时候,使用的值实际上是整型的。所以MCTK里也写了,先用double precision导出数据,接着就转为uint 16.  所以,float 64虽然有很多小数位,但都是0,即实际只用到了整数位。

  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 1
    评论

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

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值