高分四号数据计算表观反射率
前期准备:查找关于计算表观反射率的文献,总结出以上的几点:
(1)需要偏移和增益数据,在中国资源环境卫星获取。
需要注意的是,这个定标参数每年都会变化,需要根据数据日期选择,如下图所示。
(2)计算表观反射率的公式,那些参数需要根据自己使用的数据计算出,不可以直接套用我的计算出来的结果。
(3)所有的操作都是在ENVI上的band math的工具上进行的,因此,操作相对来说比较简单。
下面开始进行实践:
使用ENVI5.3软件完成,如上图所示:使用国产的小插件打开原始数据的PMS文件,打开之后如图所示。
查看数据的波段数量,这里需要注意的是,可以看出有5个波段,但是你用到什么波段,就对那个单波段处理,可以节省一大部分时间。
第一个波段是全色的,千万不要搞错。
在本次实验中,我只用到红色波段和近红外波段,因此,我只是用最后两个波段,B3是红包段,B4是近红外。
对红波段进行求辐射亮度值的操作,由表中可知,OFFSET的值都是0,因此,这个值就不加了。
所以,计算算法就是下图所示部分。
一定要选中红色波段输出,如下图所示。
命名是有所讲究的,至少让自己或者他人一下子就看懂就行。
输出后图像,上图所示。
同样的对近红外波段是同样的操作处理,但是近红外的Gains与红波段不是一样的,需要看清楚。
同样的,波段选择和命名要注意。
处理结果,如上图所示。
接下来计算表观反射率:输入计算公式(注:这个公式是自己根据影像数据算出来的,不是通用的),选择的波段是上个步骤的辐射亮度值。
波段选择和命名要注意,TOA是表观反射率的缩写!下图是输出结果。
近红外波段也是同样操作。
目前的结果都是单波段,可以使用Layerstacking 工具进行波段合成,使其变成两个波段的数据。
在工具栏里搜索即可。
选择两个波段,输出即可。本次实验只是针对两个波段进行,可以使用同样的操作对绿、蓝两个波段进行处理。
高分四号数据计算表观反射率相关推荐
- 【技术类】【ArcGIS对国产卫星的支持2:高分一号卫星】篇5、实时表观反射率产品生产
表观反射率就是指大气层顶的反射率,辐射定标的结果之一,又称视反射率.英文表示为:apparent reflectance (=地表反射率+大气反射率) 分发商提供的数据基本都是灰度值(DN值),必须转 ...
- 遥感中的表观反射率、反射率、反照率等之间的区别与联系
遥感反射率的定义:地物表面反射能量与到达地物表面的入射能量的比值. 遥感表观反射率的定义:地物表面反射能量与近地表太阳入射能量的比值. 大气校正就是将辐射亮度或者表观反射率转换为地表实际反射率,目的是 ...
- 辐射亮度与表观反射率的转换
文章目录 DN值 辐射亮度 表观反射率 转换公式 亮度温度 其他资料 补充相关资料 DN值 就是影像灰度值 辐射亮度 略 表观反射率 大气层顶的反射率 转换公式 DN值转辐射亮度: Lλ=Gain∗D ...
- 计算GF1和GF6的表观反射率和投影问题
记录一下怕自己忘了 主要是之前一种用classic做出来的一直没有坐标信息 很烦 于是想了此下策 前提是要装了中国卫星的插件比较方便 首先用ENVI将需要的波段转成envi标准格式 不用classic ...
- 反射率(表观反射率、地表反射率)
反射率是物体表面所能反射的辐射量和它所接受的辐射量的比值,一般在[0,1]范围,有的时候为了储存方便而扩大一定的倍数,如放大一万倍[0,10000].一些材料可以通过他们的反射光谱来识别,因此,为了更 ...
- IDL调用6S文件生成数据,及表观反射率有关二维图像
;;=============================================================================================== ;; ...
- 利用Python和GDAL把MODIS的HDF格式表观反射率产品MOD09A1转换为TIFF格式并重投影到与LANDSAT8-OLI一致
解释 1.MOD09A1是表观反射率产品,是HDF格式并且投影和常用的Landsat数据不一致.为了统一反演时使用,想要将该产品转换为TIFF格式并重投影.特别是要批量的进行转换和重投影. 2.采用P ...
- RSD处理高分5号高光谱(GF5 AHSI)数据(一)——正射校正和大气校正
RSD对高光谱数据处理增加了一种层类型,在原来的SWATH.GRID基础上增加了一种HYPER,HYPER大部分结构与GRID基本相同,主要区别是GRID处理数据通道上限是40个,HYPER增加到40 ...
- [RS] 基础概念区分:DN-辐射率-反射率
[反射率] 是指任何物体表面反射阳光的能力.比如说某物体的反射率是45%,这意思是说,此物体表面所接受到的太阳辐射中,有45%被反射了出去 在遥感中,有好几种反射率:地表反射率.表观反射率等等 &qu ...
- 实习一:6S大气校正
********************************************************************************************* 以下内容只用 ...
最新文章
- 【万级并发】电商库存扣减如何设计?不超卖!
- mac 中 Mach-O
- 八杯水微博程序演示及下载
- Java for LeetCode 042 Trapping Rain Water
- 丹琦女神新作:对比学习,简单到只需要Dropout两下
- 轻量级开源内容管理系统-YzmCMS v6.0源码
- 28线性空间02——坐标、坐标变换与基变换、过度矩阵
- LintCode Coins in a Line II
- linux 正向代理与反向代理
- 用计算机刻盘,电脑怎么刻盘_怎么将电脑桌面的文件刻录到cd光盘
- 龙华大浪注塑加工中模具原理及组成
- 郑明秋什么版本的MySQL_mysql数据库实用教程教学课件作者郑明秋代码数据库脚本代码9787568250825.docx...
- 第1期技术: DQN算法原理及实现过程
- 同步六进制加法计数电路(D触发器)
- lidatr2rosbag工具包,KITTI数据转换ros bag文件
- 尼尔机械纪元房顶机器人_尼尔机械纪元辅助机器人是什么 辅助机器人解析攻略...
- 虚拟机(一)虚拟机安装Mysql
- STM32关于使用定时器来实现串口通信的整活实验
- 转:目标50行代码之内完成3d编辑器功能
- python爬虫三大解析数据方法:bs4 及爬小说网案例
热门文章
- 计算机系统的组成和应用,计算机系统组成与应用领域的考点
- cron每两个小时_crontab 每隔 1 小时 2 小时的执行 job 写法
- python字符串对齐.format_python字符串格式化之.format
- java单元测试模拟输入_java – 单元测试:在定义模拟行为后调用...
- postman添加cookie_有了这款IDEA插件,再也不需要postman了
- android多线程计时器,Android 计时器Timer用法
- mysql 8.X主从复制
- 自动驾驶_产品象限图_清洁机器人场景
- Tensorflow训练渐渐变慢,迭代一段时间卡死
- 主干网络系列(1) -ResNet: Deep Residual Learning for Image Recognition