注册 登录
  • 注册时,本站名称为:RGB空格3S博客,注意中间的空格。
  • 2018.10.14本人利用业余时间接遥感数据处理与编程小活,欢迎骚扰,QQ:853060844
  • 2018.8.14本站全面接入google广告
  • 本人闲暇时间(周末及其他闲暇时间)接遥感图像处理与IDL编程小活,欢迎骚扰,qq:853060844
  • 2017.2.14今天收到45条恶意评论(全是外文),故评论时请填写必要信息,匿名评论全部拉黑,迫不得已而为之
  • 2017.1.27,2017年春节及至,我谨代表本人祝大家新春快乐,本人年终总结文章请访问:2016年终总结
  • 为防止恶意转载,本站全面禁止复制,并添加图片水印:RGB 3S博客www.ixxin.cn。
  • 本站正式更名为RGB 3S博客,本站将撤消所有非3S内容,其将转移到新博客江湖时代

国产高光谱传感器HSI去条带几何配准以及大气校正(最完整预处理)

技术文档 admin 2519次浏览 已收录 0个评论
[隐藏]

数据简介

转载请在文章开头注明出处。
环境一号卫星HJ-1A星上搭载的HSI超光谱成像仪能完成对地刈宽为50公里、地面像元分辨率为100米、光谱范围(0.45-0.95µm)和115个光谱谱段的推扫成像。
操作软件ENVI5.3.1:https://www.ixxin.cn/2016/12/14/envi-5-3-sp1pojie/

预处理步骤

读取影像

ENVI5.3.1直接open HSI数据的.H5文件即可。
open
点击ImageData,双击Banddata即可加载到右边。
open2
同时获取以下属性信息:
以下信息用于FLAASH:
Scalefactor:100(定标值系数)
MapInformation→SceneCenterLat(纬度) SceneCenterLong(经度)
MapInformation→PixelSpacing(像元大小)
ImageAttributes→SceneStartTime(数据获取时间)
以下信息用于几何配准:
MapInformation:
Projection(投影)
Utm-zone(带号)
ProductUpperLeftX(左上角X坐标)
ProductUpperLeftY(左上角Y坐标)
单击Open Rasters按钮,打开图像文件。
点击File——save as 保存为ENVI标准格式。

几何配准

需要在Classic中进行。
加载上一步保存的影像。
右键数据,选择Edit Header。
点击Edit Attributs。
选择Map info
依次选择投影信息,椭球体:WGS84,填写带号,填写左上角坐标,填写像元大小100.(所有参数上一步均已获得)
mapinfo

条带去除

HSI前面波段数据质量很差,有很多条带,如图:
HSI条带
第一步,将条带旋转到垂直状态:
用到工具:rotate/flip data
角度为:-14.39
旋转
第二步去除条带:
用到的工具:spear vertical removal
参数设置为100
去条带
去条带后对比:
对比
第三步,旋转到原始状态
用到工具:rotate/flip data
角度为:14.39
去条带

FLAASH大气校正

第一步,数据转换
FLAASH需要BIL存储方式
在Toolbox工具箱中,双击Raster Management→Convert Interleave工具,选择刚打开的图像文件,在Convert File Parameters中,设置Output Interleave:BIL,选择输出文件名“HSI_BIL.dat”和路径,单击OK按钮。
第二步:FLAASH大气校正
(1)在Toolbox工具箱中,双击Radiometric Correction/Atmospheric Correction Module /FLAASH Atmospheric Correction工具,单击Input Radiance Image按钮,选择上一步保存的“HSI_BIL.dat”图像。
(2)在打开的Enter ASCII Filename containing wavelengths and fwhm对话框中,选择含有波段信息的txt文件HJ1A-HSI-BandInfo.txt,并在Input ASCII Filename对话框中手动输入Wavelength和FWHM的行数,如图所示。
提示:HJ1A-HIS数据的文件夹中包含一个文本文件,如本示例中的“HJ1A-HSI-2-65-A2-20110820-L20000597997_BandInfo.txt”,其中包括了每个波段的波谱范围和波段宽度。“HJ1A-HSI-BandInfo.txt”就是基于这个文件生成的,即第一列和第二例的平均值作为中心波长(wavelenth),第三列作为波段宽度(FWHM)。“HJ1A-HSI-BandInfo.txt”下载地址:http://pan.baidu.com/s/1o64sVzG
波长
(3)在出现的Radiance Scale Factors对话中选中Use single scale factor for all bands,缩放系数填写:1000。
提示:下载的HJ1A-HIS数据已经经过了定标,为了减少储存空间而扩大了100倍。单位是100* W/(m2 *µm * sr)。
(4)单击Output Reflectance File按钮,选择输出文件名和路径。
(5)单击Output Directory for FLAASH Files按钮,选择文件输出目录。
(6)传感器与图像信息
lLat:40.4616° Lon:116.195°
l传感器类型(Sensor Type):UNKNOWN-HIS
l传感器飞行高度Sensor Altitude:650
l图像区域平均海拔(Ground Elevation):0.5
l成像时间(Flight Date):2011-08-20 Flight Time:03:03:47.96。
(7)大气参数
l模型(Atmospheric Model):Mid-Latitude Summer
提示:根据图像中心纬度40°和8月份成像从帮助表中查询获取。
l执行水汽反演(Water Retrieval):Yes。
lWater Absorption Feature:820nm
l气溶胶模型(Aerosol Model):Urban
l气溶胶反演(Aerosol Retrieval):None
提示:HJ1A-HIS数据缺少短波红外波段,采用能见度值参与气溶胶计算。
(8) 其他按照默认:
• 初始能见度(Initial Visibility):30。
• Spectral Polishing:yes
• Width(number of bands):5
• Wavelength Recalibration:No
(9) 单击高光谱设置(Hyperspectral Settings)按钮,选择默认设置。
(10)单击Advanced Settings按钮,打开高级设置面板。修改以下几个参数,其他按照默认:Tile Size:100
(11)参数如图所示。单击Apply按钮,执行FLAASH处理。
FLAASH
结果:
结果

后语

参考文章:http://blog.sina.com.cn/s/blog_764b1e9d0102uxde.html
打赏直接扫描二维码,打赏名单请看:https://www.ixxin.cn/about.html


xxin blog , 版权所有丨如未注明 , 均为原创丨本网站采用BY-NC-SA协议进行授权 , 转载请注明国产高光谱传感器HSI去条带几何配准以及大气校正(最完整预处理)
喜欢 (2)
支付宝[853060844@qq.com]
分享 (0)
admin
关于作者:
坐标山科大遥感系小鲜肉一枚。
发表我的评论
取消评论
表情 贴图 加粗 删除线 居中 斜体 签到

Hi,您需要填写昵称和邮箱!

  • 昵称 (必填)
  • 邮箱 (必填)
  • 网址