参考文献 :         基于 HJ 卫星的北京市气溶胶光学厚度遥感反演      冉凡坤

深蓝算法反演气溶胶光学厚度的依据>>>

在亮像元地表,蓝光波段的地表反射率较小,只有可见光波段的 1/2~1/4,可以通过对蓝光波段地表反射率的模拟实现对气溶胶光学厚度的反演。同时,深蓝波段的分子散射和吸收都远远的低于其它波段,气溶胶垂直廓线的影响对深蓝波段气溶胶光学厚度反演的影响甚微。

反演算法>>>
1. 计算卫星在  412、490 和 670nm 的表观反射率;

2. 对卫星数据进行高程修正,数据基于数字高程;

3. 基于蓝波段的反射率变化和吸收气溶胶指数对云进行提取;

4. 建立412nm、490nm 和 670nm 的地表反射率库,地表反射率取的是晴朗无云条件下的地表反射率;

5. 基于地表反射率和辐射传输模型 6S,通过建立查找表得到气溶胶光学厚度;

敏感度测试>>>

当地表反射率较小时(< 0.2),卫星观测到的表观反射率与气溶胶光学厚度有着很好的线性相关性,且敏感度很大。

但随着地表反射率的增大(>0.2),卫星观测到的表观反射率与气溶胶光学厚度的相关性几乎为 0。

因此当利用深蓝算法进行气溶胶光学厚度反演时,必须保证地表反射率不高于 0.1。

蓝波段地表反射率的分布比较集中,80%介于 0~0.1 之间,60%小于 0.05,只有 10%高于 0.2    =====》   蓝波段反演AOD度敏感度非常高

参考文献 :         基于深蓝算法的 HJ-1CCD数据快速大气校正模型        王中挺 等

原理>>>> 蓝波段大气反射较强,而地表反射较弱。将清晰天的地表反射率(即地表反射率库)代入地表反射率公式,反演气溶胶。

备注>>>> 当地表反射率<0.1,反射率随AOD有较好的线性关系;

当地表反射率>0.2,基本无变化。                 =====》去云阈值的设立依据

MODIS地表反射率产品>>>>包括250  m单日合成、500 m与1 km单日合成、250 m 8日合成、500 m 8日合成、CMG网格全球日合成(分辨率为0.05°)等5种形式

为消除地表反射率的短期变化,选择8日合成产品进行地表反射率库的构建。该产品采用正弦投影,将全球分为10°×10°的小块存放。

将MODIS的8日合成地表反射率产品进行波段提取、重投影、重采样、图像剪切和镶嵌等处理得到地表反射率库。

&&&对地表反射率数据进行综合统计,分析蓝波段(MODIS第3波段)地表反射率在不同值域区间的分布情况。

&&&得到蓝波段地表反射率累积百分比。====》蓝波段地表反射率整体较小。

====》在地表反射率小于0.1时,表观反射率随着AOD有着较好的线性关系;

在地表反射率大于0.2时,表观反射率随着AOD基本没有变化。

因此,当地表反射率小于0.1时,本算法能够较好的运行。

查找表>>>>1152个观测几何,12个太阳天顶角(0—66°,间隔6°)、6个观测天顶角(0—35°,间隔为7°)、16个太阳与卫星之间的相对方位角(0—180°,间隔12°);

大气气溶胶模式参数假设为大陆型气溶胶,并设立6个大气气溶胶光学厚度值(0.0001、0.25、0.5、1、1.5和1.95);

波段自定义为HJ-1  CCD相机第一波段的响应函数;暂时不考虑地表反射,反射率设为0。

????如果不使用HJ-1CCD数据,如何确定自定义波段的响应函数????

HJ-1CCD观测几何的转换>>>.

辐射定标及云去除>>>>首先,从CCD数据的辅助xml文件中读取辐射定标系数,将DN值转换为表观辐亮度;

然后,根据上一步生成的太阳天顶角,,将获得的表观辐亮度转换为表观反射率。

云的去除则是通过阈值实现。研究表明(Acker-man  等,2010),在红光波段卫星观测到的云的反射率一般大于0.2====》选取0.2作为云去除阈值。

深蓝算法反演AOD入门记录(一)相关推荐

  1. 使用MODIS Level 1B 1KM 数据反演AOD实验流程

    1.实验思路 实验的整体思路: 第一部分是MODIS数据预处理工作. 数据预处理工作主要包括对 MODIS LIB数据进行去除蝴蝶结效应.几何校正.波段合成.裁剪等操作 第二部分为城市地区AOD及PM ...

  2. 算法竞赛从入门到进阶pdf_好书送不停 | 算法竞赛入门到进阶

    本书内容: 本书是算法竞赛的入门和进阶教材,包括算法思路.模板代码.知识体系.赛事相关等内容.本书把竞赛常用的知识点和竞赛题结合起来,讲解清晰.透彻,帮助初学者建立自信心,快速从实际问题入手,模仿经典 ...

  3. jieba入门记录——nltk中文语料处理

    jieba入门记录--nltk中文语料处理 环境:pycharm(Anaconda) 1.官网下载jieba:https://pypi.org/project/jieba/ 2.将下载好的压缩包解压到 ...

  4. 【玩转数据系列十三】机器学习算法基于信用卡消费记录做信用评分

    原文链接 机器学习算法基于信用卡消费记录做信用评分 背景 如果你是做互联网金融的,那么一定听说过评分卡.评分卡是信用风险评估领域常用的建模方法,评分卡并不简单对应于某一种机器学习算法,而是一种通用的建 ...

  5. 【Envi】基于单窗算法的地表温度反演实验操作记录

    文章目录 比前言还前 前言 参考博客 技术流程 1.数据预处理 1.1数据获取 1.2辐射亮度温度 1.2.1热红外波段辐射定标与亮度 1.2.2辐射亮度温度计算 1.3地表比辐射率 1.3.1多光谱 ...

  6. 算法与数据结构入门一篇就搞定

    算法与数据结构开篇 你真的会数据结构吗? 公司开发一个客服电话系统,小菜需要完成客户排队模块的开发,经过三次修改: 第一次:小菜使用了数据库设计了一张客户排队表,并且设置了一个自动增长的整型id字段, ...

  7. 算法基础、算法比赛快速入门(java)

    想用Java快速入门算法?这篇文章你得看! 提示:本文章适合想要入门算法,并且想 "快速" 达到一定成果的同学们阅读~ 文章非常非常非常长(可能是你见过最长的算法基础篇章)!!! ...

  8. 字符串算法:从入门到劝退

    大家好,我是黄人和,一只小菜鸡.本文长期更新,记得收藏点赞关注三连 前几日王烨程问我KMP的缘故,今天中午他又在问我这个.之所以很久没有回答他,一来是近来事务繁忙,未尝闲暇,二来是我也对字符串算法有点 ...

  9. python dfs算法_python从入门到进阶,就靠这份书单了

    自2018年3月起,在全国计算机二级考试中加入了"Python语言程序设计"科目. 从 2018 年起,浙江省信息技术教材将不会再使用晦涩难懂的 VB 语言,而是改学更简单易懂的 ...

最新文章

  1. Android入门:DDMS视图
  2. .NET Core 实例接口代理转发
  3. Python保存最后N个元素
  4. 高级cmd攻击命令_一步一步学习DVWA渗透测试(Command Injection命令行注入)-第七次课...
  5. c 调用openoffice word转pdf_批量 Word 转 PDF 方法
  6. 前端学习(2703):重读vue电商网站24之配置axios
  7. python绘图颜色深浅代表数值_画图理解Python的深浅拷贝
  8. Lintcode 138. 子数组之和 [Easy]
  9. Flexsim在固定资源类中没有分拣传送带?
  10. Linux课程设计:Linux系统下多进程的创建与通信
  11. android微信支付+指纹支付密码错误,为什么微信付款无法指纹支付密码?如何开启指纹支付?...
  12. 铁道部售票网站怎么啦?
  13. 为Win7关机按钮菜单增加“休眠”选项
  14. mysql实验训练2 数据查询操作_实验训练2:数据查询操作.doc
  15. win10安装SQL Server 2017 遇到的一大堆问题(+SSMS)【最后用虚拟机解决】
  16. java学习体会论文_关于学习java的心得体会
  17. Java之BigInteger的用法
  18. [网络工程师]-防火墙-防火墙技术
  19. 转换成大写日期(中文)
  20. Android筑基——BroadcastReceiver 的动态注册、发送和接收过程(基于api21)

热门文章

  1. 做了一个app,返回三国武将的排序网页,记录一下,省得以后找不到
  2. [指针八]有关指针的面试题
  3. 执著如泪,是滴入心中的破碎
  4. Linux学习-02-Linux的安装
  5. (20201209已解决)从window访问wsl地址
  6. 用计算机画函数,用计算机画函数图像 优质课教案设计
  7. ILM(Index Lifecycle Management)
  8. 三星pay显示无法连接服务器,三星智付Samsung pay正式开启:服务器宕机 三星紧急发公告...
  9. c语言点餐怎么变量相加营业额,c语言点餐系统软件.docx
  10. 潍坊学院计算机系宿舍,关于表彰潍坊学院第十届大学生宿舍文化节先进集体和先进个人的决定...