作者:林科院-姜怡
一、数据下载
HJ卫星的数据可以在-陆地观测卫星数据服务平台-上进行下载
http://218.247.138.119:7777/DSSPlatform/productSearch.html
二、数据读取
下载插件openhjccd
地址:https://www.ixxin.cn/2018/01/06/envi53openhjccd/
将插件装在envi拓展工具安装目录下:D:\ENVI\ENVI53\extensions
打开ENVI,选择File-open as-HJ CCD,即可打开环境卫星数据
但是如果要进一步进行预处理的话需要做以下工作:
1、将刚刚读取的环境卫星数据重新另存为envi文件的格式,选择file-save as等
2、定标系数的更改
重新打开刚刚存储的环境卫星数据,即:file-open或者直接拖入。更改图片的定标信息:
操作步骤:右击图像,选择View Metadata,选择Spectral,选择edit metadata,更改Data gain values与Data offset values中的数据,定标系数的更改值可以在中国资源卫星数据中心进行下载,找到与自己数据匹配的年份与对应的传感器,http://www.cresda.com/CN/Downloads/dbcs/6709.shtml
值得注意的是,由于HJ卫星的某些原因,我们要根据具体情况更改定标系数,如果公式是 ,那么就直接按照中国资源卫星数据中心中的系数直接填写,就不必更改;如果公式中是DN/gain,那么记得将gain取倒数。
三、辐射定标
在Toolbox工具箱中,双击Radiometric Correction/Radiometric Calibration工具。在File Selection对话框中,选择要处理的数据,接下来是参数设置,类型选择radiance,点击Apply Flaash Settings ,以便今后的进行大气纠正,选择储存路径即可。
四、构建波谱响应函数(参考文献:ENVI-IDL技术殿堂的博客)
在中国资源卫星数据中心可以下载,数据是以四列数字形式提供,第一列表示波长,在envi中,使用波谱曲线来描述波谱响应函数,也就是以波长作为x轴,波谱响应值作为y轴,存储格式为ENVI波谱库文件(.sli)。下面介绍ENVI中的波谱响应函数的制作。
(1) 在传统界面中,选择Window->Start New Plot Window,ENVI Plot Window窗口中,选择File->Input Data->ASCII,如图2所示,自动将第一列作为X轴,后面4列作为Y轴。单击OK。
(2) 如下图所示,生成了4条曲线。选择Edit->Data Parameters,更改每一条曲线的名称:b1,b2,b3,b4,便于区分。
(3) 选择File->Save Plot As->Spectral Library,将波谱曲线保存为波谱库文件。
五、大气纠正
(1) 工具箱/Radiometric Correction/Atmospheric Correction Module/FLAASH Atmospheric Correction打开FLAASH大气校正模块;
(2) 点击Input Radiance Image,选择BIL格式的环境小卫星数据,在Radiance Scale Factors面板中选择Use single scale factor for all bands,单击OK;
注:定标后的单位是W⋅m−2⋅sr−1⋅μm−1,与FLAASH要求的单位(μW)/(cm2*nm*sr)相差10倍关系,因此在Radiance Scale Factors中输入10缩放系数。
(3) 设置输出文件及路径设置;
(4) 传感器基本信息设置:
• 成像中心点经纬度FLAASH自动从影像中获取。
• 传感器高度(Sensor Altitude):650km
• 像元大小(pixel Size):30m
• 成像区域平均高度可以通过统计DEM数据获取(envi中有自带的900M DEM数据)
• 成像时间:从数据头文件中读取(.XML);
(5) 大气模型和气溶胶模型,根据经纬度和影像区域选择
(6) 气溶胶反演方法选择None(缺少短波红外),能见度设置为40km。
(7) 单击Multispectral Setting按钮,在Filter Function File 导入光谱响应曲线“环境1A星CCD2光谱响应.sli”,单击OK;
(8) 单击Advanced Settings,在高级设置中,Tile Size 默认的是Cash size 的大小,手动改为50-100Mb(根据内存大小设定),单击OK;
(9) 设置好后,在大气校正模块面板中,单击Apply。
(10)大气校正完成后,检查大气校正的结果,分别加载校正前后的图像(选择CIR假彩色方式加载,可以更好的识别植被),查看典型地物的大气校正前后的光谱曲线。
校正前:
校正后:
3、参考文献:
ENVI下环境小卫星图像处理——ENVI-IDL技术殿堂的博客
ENVI下环境小卫星图像处理——环境小卫星CCD数据(资源一号、三号)的大气校正(ENVI5.0及之后版本)