作业标准1：辐射定标及波段运算
铁路外部环境安全隐患变化卫星监测技术
特别策划·铁路科技保安全铁路外部环境安全隐患变化卫星监测技术杨移超（中国铁路北京局集团有限公司北京高铁工务段，北京100070）摘要：为保障铁路运营安全，科学、高效、定期开展铁路周边安全环境整治工作，需大力提升技防手段，其中卫星监测技术可以快速准确获取铁路外部环境安全隐患的地理位置和属性特征，有效降低人力、物力投入。
因此应大力推广铁路外部环境安全隐患变化卫星监测技术应用，同时要及时跟进作业流程规范化梳理，依托京沪高铁廊坊段卫星监测技术应用，详细论述技术流程、数据预处理、隐患影像特征、项目内容，形成一整套满足工程应用的作业指导。
相比较传统的现场筛查安全隐患变化，卫星监测技术获取的安全隐患问题库更精细具体，规范化安全隐患变化卫星监测作业流程可以更好地提升监测精度和工作效率。
关键词：铁路外部环境；安全隐患；卫星监测；数据预处理；影像特征；京沪高铁；铁路安全中图分类号：U298 文献标识码：A 文章编号：1001-683X（2023）10-0151-07 DOI：10.19549/j.issn.1001-683x.2023.01.10.0030 引言卫星监测技术在土地利用、水资源、林业、采矿等行业均已系统化运用［1-4］。
随着陆地卫星观测体系不断健全完善，覆盖能力不断提升，为全天候、全要素、全流程监测提供了重要数据基础。
近年来，随着高速铁路建设快速发展，铁路营业里程不断增长，铁路外部环境安全隐患时刻威胁着列车运行安全。
铁路外部环境安全隐患按问题属性主要分为18个类别，包括违法施工、危险物品、上跨并行等。
在这些问题类别中，有学者研究部分安全隐患如硬飘浮物、轻飘浮物、堆放隐患等适合采用卫星监测手段进行检测调查［5-6］。
铁路外部环境安全隐患造成的列车事故通常由外力因素触发，具有随机性和偶发性，精确排查出安全隐患的地理位置及自身属性后，综合考虑隐患自身属性、环境特征和空间关系可以进行安全风险评价。
铁基金项目：中国国家铁路集团有限公司科技研究开发计划项目（N2020T005）作者简介：杨移超（1983—），男，高级工程师。
辐射定标教程
在新版本ENVI中（4.5版本以及更新），没有单独设立ASTER和MODIS定标的工具。
对于ASTER L1A/L1B和MODIS 02级数据，在打开数据时会自动完成对数据的定标。
如图1所示打开ASTER L1B的结果，在波段列表中，自动读取各个波段的中心波长信息，并按照波段范围信息（VNIR、SWIR、TIR）分组波段。
其中VNIR、SWIR自动定标为辐射亮度，单位是：W/m2/sr/μm；TIR数据定标为大气表观温度值，单位：开尔文。
打开其中一个数据，浏览像元值，可以看到已经定标为浮点型的辐射亮度值。
图1 ASTER L1B数据如图2为打开MODIS 02级1km数据，其中250米和500米的波段经过重采样为1km 加入这个数据集中。
ENVI根据各个波段的中心波长信息定标为三个类型数据：反射率数据（Reflectance）、辐射亮度值数据（Radiance）和发射率数据（Emissive）。
其中反射率和发射率为0~1无单位值，辐射亮度值单位是：W/m2/μm/sr。
图2 MODIS 02级数据如果打开原始的ASTER和MODIS的DN值数据，可以在ENVI主菜单中选择File->preferences，切换到Miscellaneous面板，将Auto-Correct ASTER/MODIS项设置为NO。
3.2 Landsat数据定标ENVI4.7版本改进了Landsat数据定标的功能，对于Landsat4/5数据可以手动选择以下两种定标公式：式中:•QCAL为原始量化的DN值•LMINλ为QCAL = 0时的辐射亮度值•LMAXλ为QCAL = QCALMAX时的辐射亮度值注：LMINλ和LMAXλ的值取自Chander, Markham, and Helder (2009)的研究成果。
•QCALMIN是最小量化定标像素值（与LMINλ类似）。
取值如下：1：LPGS产品1：04 April 2004之后的NLAPS产品0：04 April 2004之前的NLAPS产品注：如果没有元数据信息，QCALMIN取默认值1（TM和ETM+)）或者0 (MSS)。
辐射定标
线形定标公式
定标过程一般采取线形公式进行转换： 定标过程一般采取线形公式进行转换： L = a*DN + b a(gain)、b(offset)通常可以从遥感数据头文件读出 a(gain)、b(offset)通常可以从遥感数据头文件读出 L
线形区域
DN
二、辐射纠正—反射率 辐射纠正— 的计算
电磁辐射与辐射源——地物 电磁辐射与辐射源——地物——传感器的几何 地物——传感器的几何 关系 水平地面的假设 山地辐射纠正 辐亮度向反射率的转换
定标参数的确定
• 定标公式针对何种真实物理量 反射率？ 反射率？ • a是负数时，定标前后图象视觉相反 是负数时， ？ • 定标参数的确定都是对波段的波长积分 响 应 函 数 λ
遥感数据的星上辐射定标
DN值的影响因素 DN值的影响因素 参考光源 地面定标测量 遥感数据的辐射定标——地表辐亮度的计算 遥感数据的辐射定标 地表辐亮度的计算
公式的GAINS/BIASES可从头文件中获取
%1=%1*0.00398-0.0100 %2=%2*0.00964-0.0232 %3=%3*0.00540-0.0078 %4=%4*0.01043-0.0193 %5=%5*0.00235-0.0080 %6=%6*0.05516+1.2378 %7=%7*0.00154-0.0040 TM6的 TM6的GAINS/BIASES 是不变的，为：0.055158 / 1.2378 是不变的，
式中除10为将单位由 式中除10为将单位由 mW /(cm ⋅ Sr ⋅ µm ) 转化为 W /( M ⋅ Sr ⋅ µm )
2
2
实际操作
计算出各波段的反射辐亮度 计算出各波段的行星反射率 对比地物反射辐亮度与行星反射率 请考查计算后各典型地物反射率的特点， 请考查计算后各典型地物反射率的特点，并 分析原因
遥感原理与应用第6章-遥感作业(1)
第六章遥感图像辐射校正名词解释：辐射定标、绝对定标、相对定标、辐射校正、大气校正、图像增强、累积直方图、直方图匹配、NDVI、图像融合1、辐射定标：是指传感器探测值的标定过程方法，用以确定传感器入口处的准确辐射值。
2、绝对定标：建立传感器测量的数字信号与对应的辐射能量之间的数量关系，对目标作定量的描述，得到目标的辐射绝对值。
3、相对定标：又称传感器探测元件归一化，是为了校正传感器中各个探测元件响应度差异而对卫星传感器测量到的原始亮度值进行归一化的一种处理过程。
最终得到的是目标中某一点辐射亮度与其他点的相对值。
4、辐射校正：是指消除或改正遥感图像成像过程中附加在传感器输出的辐射能量中的各种噪声的过程。
5、大气校正：是指消除大气对阳光和来自目标的辐射产生的吸收和散射影响的过程。
6、图像增强：为了特定目的，突出遥感图像中的某些信息，削弱或除去某些不需要的信息，使图像更易判读。
7、累积直方图：以累积分布函数为纵坐标，图像灰度为横坐标得到的直方图称为累积直方图。
8、直方图匹配：是通过非线性变换使得一个图像的直方图与另一个图像直方图类似。
也称生物量指标变化，可使植9、NDVI：归一化差分植被指数。
NDVI=B7−B5B7+B5被从水和土中分离出来。
10、图像融合：是指将多源遥感图像按照一定的算法，在规定的地理坐标系，生成新的图像的过程。
问答题：1．根据辐射传输方程，指出传感器接收的能量包含哪几方面，辐射误差及辐射误差纠正内容是什么。
根据辐射传输方程，传感器接收的电磁波能量包含三部分：1）太阳经大气衰减后照射到地面，经地面发射后又经过大气的二次衰减进入传感器的能量；2）大气散射、反射和辐射的能量；3）地面本身辐射的能量经过大气后进入传感器的能量。
辐射误差包括：1）传感器本身的性能引起的辐射误差；2）大气的散射和吸收引起的辐射误差；3）地形影响和光照条件的变化引起的辐射误差。
辐射误差纠正的内容是传感器辐射定标和辐射误差校正等。
辐射定标步骤
辐射定标步骤
辐射定标是指通过比较已知源的辐射强度和待测源的辐射强度来确定待测源的辐射水平的过程。
辐射定标的步骤如下：
1. 准备辐射源：选择已知辐射强度的辐射源作为定标源，例如放射性同位素或已经经过精确测定的辐射源。
2. 测量定标源的辐射强度：使用辐射测量仪器（例如辐射计或辐射剂量仪）对定标源进行辐射强度的测量。
确保使用的仪器具有良好的准确性和稳定性。
3. 确定仪器的响应：根据已知辐射强度和仪器测量值的对比，计算出仪器的响应因子。
这个响应因子可以用来将仪器的测量值转换为辐射强度。
4. 测量待测源的辐射强度：使用同一台仪器测量待测源的辐射强度。
确保在测量中采取适当的防护和安全措施。
5. 计算待测源的辐射水平：根据仪器的响应因子和待测源的测量值，计算出待测源的辐射强度。
可以使用之前计算得到的响应因子将仪器的测量值转换为相应的辐射水平。
通过以上步骤，可以将待测源的辐射水平与已知源进行对比，从而确定待测源的辐射水平。
这个过程对于辐射监测、核安全等领域具有重要意义。
辐射定标
具体的定标公式：
LMAX LMIN L LMIN ( )(QCAL QCALMIN) QCALMAX QCALMIN
式中，QCAL为原始量化的DN值，为QCAL=0 时的辐射亮度值。为QCAL=QCALMAX时的辐 射亮度值，QCALMAX/QCALMIN是最大/最小 量化定标像素值。一般情况下，QCALMIN=0。 所以公式又可简化为 L LMIN LMAX LMIN QCAL
splitfilenaห้องสมุดไป่ตู้me()
getMaxDN ()
getMinDN()
调用后返回存放DN值数 组地址
getDNvalu e() 调用后返回存放DN值数 组地址 getGrayVa lue())
图1 定标算法函数间逻辑关系
功能1：打开原始文件 （1）操作->打开原始数据->选择“pp_01.txt”，打开
图2 打开文件“pp_01.txt”
（2）显示测试原始图像和定标信息
图3 原始数据图像显示。
如果对原始数据文件打开并定标写入成功，则会 弹出下面的提示框：
图4 提示信息
功能2：定标 步骤：定标->线性定标->选择“file_01.txt”，打开； 下图是定标后的结果显示：
图7 测试文件定标结果
QCALMAX
其中增益值gain=
LMAX LMIN QCALMAX
偏移量offset= LMIN
定义类：class Calibration
showPictu re() 调用后返回特 征值 调用后返回特征 QCALMAX GetChara cter() 调用后返回特征 QCALMIN readfile() 调用后返回特征 文件名
(完整)基于单窗算法反演地表温度的ENVI操作教程
单窗算法反演地表温度教程1.1 算法原理1.1.1单窗算法单窗算法（MW 算法）是覃志豪于2001年提出的针对TM 数据只有一个热红外波段的地面温度反演算法。
经过众多学者验证，单窗算法具有很高的反演精度，且同样适用于ETM+和landsat 8数据。
公式如下：6666666666/)))1(()1((C T D T D C D C b D C a T a sensor s -++--+--=式中，LST 为地表温度（K ），T sensor 是传感器上的亮度温度（K ），T a 是大气平均温度（K ）；a 、b 为参考系数，当地表温度为0-70℃时，a = -67.355351，b = 0.458606；C 、D为中间变量，计算公式为：式中，为地表比辐射率，为地面到传感器的大气总透射率。
因此单窗算法反演地表温度的关键是计算得到亮度温度T senso 、地表比辐射率、大气透射率和大气平均作用温度T a 。
1.1.2参数计算1.1.2.1辐射亮温计算利用Planck 公式将图像像元对应传感器辐射强度值转换为对应的亮度温度值。
公式如下式中，T senso 为亮度温度值；影像预处理后得到的光谱辐射值，λL 单位为，K1 、K2为常量，可由数据头文件获取。
)/(2m sr m w μ⋅⋅计算图像辐射亮温之前，需采用辐射定标参数将像元灰度值DN转换为热辐射强度值，公式如下：式中，M L 为增益参数，A L 为偏移参数，该参数可直接在影像通文件数据中获取，且ENVI 软件中已经集成，不需要自己在查找。
1.1.2.2地表比辐射率计算根据覃志豪针对TM 影像提出的混合像元分解法来确定区域地表福辐射率。
对于城市区域，我们简单的将其分为水体、自然表面和建筑表面三种，因此针对混合像元尺度上的地表比辐射率通过下式来估算：式中，为混合像元的地表比辐射率；P V 为植被覆盖率；R V 为植被的温度比率；R M 为建筑表面的温度比率；V 表示植被法地表比辐射率，m 表示建筑表面的地表比辐射率；d表示辐射校正项。
实验六-辐射定标
波段使用高增益参数，4在太阳高度角（elevation angle） 低于45度也使用高增益，反之使用低增益
7
一. 辐射定标
2. Land sat8 数据 定标：
8
一. 辐射定标——方法一
• 针对特定的卫星，如Landsat • 要求文件含有头文件信息 L******_MTL.txt（如果没
float（bx）*a+b c.为bx指定对应波段
17
一. 辐射定标——方法二
例如：TM5 2008年的数据，对第一波段
18
一. 辐射定标——方法二
例如：TM5 2008年的数据，对第一波段 辐射定标后的结果分析：图像数值的改变，由原来的 DN值整数变为有小数点的辐射亮度值
19
一. 辐射定标
⑤对每个波段进行辐射定标后，使用layer stacking再次 将定标后的各波段单个文件合并成一个文件（建议修 改头文件中的波段名称） 完成DNRadiance的转换
有则需要单个波段校正）
• Basic Tools—Preprocessing—Calibration Utilities—Landsat Calibration
9
一. 辐射定标——方法一 • 例1：Landsat5 TM • 输入原文件的MTL信息，注意选择X6的txt
10
一. 辐射定标——方法一 • 例1：Landsat5 TM • 自动读出了文件的时间 • 选择Radiance！！
11
一. 辐射定标——方法一 • 例1：Landsat5 TM • 得到了辐亮度图像
12
一. 辐射定标——方法一 • 例2：Landsat7 ETM+ • 没有集成的头文件，需要一个波段一个波段
辐射定标公式
辐射定标公式
辐射定标公式是用于计算辐射量的公式，通常用于实验室和工业领域的辐射测量。
辐射量可以通过对某种射线的能量或强度进行测量来确定，而辐射定标公式可以用来将这些量值转换为实际的辐射量。
其中最常见的辐射定标公式是剂量当量公式，即：
D = H * Q
其中，D表示剂量当量，H表示测量到的放射线能量或强度，Q表示该射线的放射性当量系数。
该公式用于计算剂量当量的单位是西弗（Sievert），它是国际单位制（SI）中用于表达辐射剂量和效应的标准计量单位。
另一个常见的辐射定标公式是能谱测量所用的峰面积公式，即：
A = K * I * F
其中，A表示测量到的峰面积，K表示测量所用的探测器常数，I 表示射线的强度，F表示测量的效率因子。
该公式适用于估算不同能量下的射线强度和对应的峰面积。
除了这些公式，还有很多其他的辐射定标公式，例如用于测量辐射源活度的公式等。
这些公式都是用于描述不同类型和目的的辐射测量，可以根据需要灵活运用。
ENVI
ENVI+IDL 上机操作一、图像数据的输入与输出（选择、读取、显示、保存）二、辐射定标（波段运算）三、地表温度的反演（单窗算法）一、图像的输入输出1、图像数据的选择与读取（以jpg文件为例）1）jpg文件的选择a、直接获取jpg文件路径语法：Result = FILEPATH( Filename [, ROOT_DIR=string][, SUBDIRECTORY=string/string_array] [, /TERMINAL] [, /TMP] )b、通过对话框打开jpg文件语法：Result = DIALOG_PICKFILE( [,DEFAULT_EXTENSION=string] [, /DIRECTORY] [,DIALOG_PARENT=widget_id] [,DISPLAY_NAME=string] [, FILE=string] [, FILTER=string/string array] [,/FIX_FILTER] [, GET_PATH=variable] [, GROUP=widget_id] [, /MULTIPLE_FILES] [, /MUST_EXIST] [, /OVERWRITE_PROMPT] [, PATH=string] [, /READ | , /WRITE] [, RESOURCE_NAME=string][,TITLE=string] )2）jpg文件的读取（二、三为数组）语法：result=read_image(fname)（借助于前面打开图像操作）语法：Result= DIALOG_READ_IMAGE ( [Filename] [, BLUE=variable][, DIALOG_PARENT=widget_id] [, FILE=variable] [,FILTER_TYPE=string] [, /FIX_FILTER] [, GET_PATH=variable] [, GREEN=variable] [, IMAGE=variable] [, PATH=string][, QUERY=variable] [, RED=variable] [ ,TITLE=string] )（打开读取一步到位）3)jpg文件的显示语法：TV, Image [, Position] [, /CENTIMETERS | , /INCHES] [, /ORDER] [, TRUE={1 | 2 | 3}] [, /WORDS] [, XSIZE=value] [, YSIZE=value]常用代码：window,1,title='彩色图像显示',$xsize=(size(data,/dimension))[2],$ysize=(size(data,/dimension))[2]tv,data,/true或者sz=size(data)ns=sz[2]&nl=sz[2]window,xsize=ns,ysize=nltv,data,true=14)图像文件的保存（写入）主要有两种情况一种是基于ENVI的二次开发envi_enter_data,envi_write_envi_file还有一种是基于IDL语法：WRITE_IMAGE, Filename, Format, Data [, Red, Green, Blue] [, /APPEND]例：IDL>fn=dialog_pickfile(title=’图像文件保存为：’) IDL>write_image,fn,’jpeg’,img二、辐射定标（波段运算）辐射定标公式为：L i=Gain i*DN i+Bias i其中L i为第i波段的辐射亮度；DN i为第i波段的灰度值；Gain i和Bias i 为第i波段定标系数可以使用IDL中的readu命令读入到内存中，然后再调用。
辐射定标及波段运算遥感实验二
测绘与海洋信息学院《遥感原理与应用技术A》实验报告实验名称：遥感图像的辐射定标姓名：学号：班级：指导教师：日期：2018-4-8地理信息系统实验室2017-2018学年第二学期一、实验目的与任务（1）了解辐射定标的原理；（2）使用ENVI软件自带的定标工具定标（3）学习使用波段运算进行辐射定标。
二、试验设备与数据设备：遥感图像处理系统数据：焦作2004年3-7和4-8数据三、辐射定标原理及目的目的：消除传感器本身的误差，确定传感器入口处的准确辐射值。
原理：辐射定标是将传感器记录的电压或数字量化值DN转化为绝对辐射亮度的过程，或者转化为与地表反射率、表面温度等物理量有关的相对值的处理过程。
四、实验过程辐射定标的结果可以是表观辐亮度（L），也可以是表观反射率(ρ) 一般有两种方式：第一种：利用计算公式，在ENVI中利用band math（波段运算）计算辐亮度或者反射率；第二种：利用ENVI自带的对TM的定标工具，进行定标，获取辐亮度或者反射率。
第一种方法：利用计算公式，通过ENVI的波段运算进行定标：1）计算表观辐亮度的公式：radiance=（（lmax-lmin）/（qcalmax-qcalmin））*（qcal-qcalmin）+lmin其中：radiance –表观辐亮度qcal-----DN（也就是影像数据本身）；lmax 和lmin是从参数表中查询；qcalmax 是DN值的最大值，对于TM是8bit来说，qcalmax=255；Qcalmin 是DN值的最小值，一般为0。
所以上面的公式针对TM数据可以简写成：radiance=（（lmax-lmin）/qcalmax））*qcal + lmin 即：上面的这个公式还可以进一步简化为：两个公式比较，可以看出，公式的中字母的对应关系。
注意：公式中需要的数据，可以通过后面的表格中查询获取2）表观反射率的计算ρ =π*L*d2/（ESUN*cos（θ））其中ρ为表观反射率；L为上一步计算出来的表观辐亮度；d为日地距离，这个数据通过下面的表格中获取；ESUN为大气层外的太阳辐射,也可以说是传感器接收处的太阳辐射;θ为太阳天顶角。
辐射定标的方法
辐射定标的方法一、辐射定标基础。
1.1 辐射定标是啥。
辐射定标啊，就像是给测量辐射的仪器找个标准的“尺子”。
咱们知道，仪器测量辐射的时候，它得知道自己测出来的值到底对应多少真实的辐射量，这就需要定标了。
这就好比一个人要称重东西，得先确定秤准不准是一个道理。
这是个很关键的事儿，如果定标没做好，那后面测量出来的数据就可能是“乱弹琴”，完全不靠谱。
1.2 为啥要辐射定标。
这原因可不少呢。
不同的仪器啊，就像不同的人看东西有偏差一样，它们对辐射的感知和测量可能会有很大的差异。
如果不定标，那数据就没有可比性。
比如说，一个研究团队用不同的设备在不同地方测量太阳辐射，要是没定标，那这些数据凑在一起就像一盘散沙，根本没法用。
而且，准确的辐射定标能让我们对环境啊、气象啊等很多方面有更精确的了解，就像给我们的研究装上了精准的导航仪。
二、辐射定标方法分类。
2.1 实验室定标。
这是一种比较传统、比较靠谱的方法。
在实验室里啊，各种条件都能控制得很好，就像在温室里种花一样。
可以给仪器提供已知的、稳定的辐射源，然后根据仪器的测量结果来确定定标参数。
这就像是给学生做小测验一样，题目都是已知答案的，就看仪器这个“学生”能不能答对。
不过呢，这种方法也有缺点，实验室的环境毕竟和实际的使用环境不太一样，就像在教室里学游泳，到了真正的泳池里可能还是会手忙脚乱。
2.2 现场定标。
现场定标就不一样了，它是在仪器实际工作的环境里进行定标。
这就好比是在战场上训练士兵，很考验仪器的“实战能力”。
现场定标可以利用一些自然的辐射源，比如说太阳。
太阳这个大光源啊，到处都有，利用它来定标很方便。
但是呢，现场的环境复杂多变，干扰因素很多，就像在集市里听人说话，周围嘈杂得很，可能会影响定标的准确性。
2.3 交叉定标。
交叉定标有点像“借东风”的感觉。
它是利用已经定标好的仪器来对其他仪器进行定标。
比如说，有个很厉害的、大家都信得过的仪器已经定标好了，那就可以用它来给其他类似的仪器定标。
辐射定标
L = F0 cosθ s R(θ s,ϕs ;θv,ϕv )
辐亮度向反射率的转换
反射率： 反射率：地物反射亮度与相同光照条件下的标 准板反射辐亮度之比。 准板反射辐亮度之比。 标准板反射辐亮度为： 标准板反射辐亮度为：
L0 = F0
π
cos θ s
故地物反射率为： 故地物反射率为：
πL(θ s , ϕ s ;θ v , ϕ v ) R(θ s , ϕ s ;θ v , ϕ v ) = F0 cosθ s
TM图像的辐射定标
辐射公式及单位
公式： 公式：
L = DN * G + B
单位： 以前为： 单位： Land 5 以前为： mW
/(cm ⋅ Sr ⋅ µm )
2
2
Land 7 以后为： W /( M ⋅ Sr ⋅ µm ) 以后为： 二者之比为1 二者之比为1：10
广州981222辐的辐射定标公式
电磁辐射与辐射源——地物 电磁辐射与辐射源——地物 ——传感器的几何关系 ——传感器的几何关系
入射辐照度下入射方向是余弦关系 对于水平地面，入射辐照度为： 对于水平地面，入射辐照度为：
F = F0 cosθ s
对单位面积，反射辐亮度与观测射方向也是余弦关系， 对单位面积，反射辐亮度与观测射方向也是余弦关系， 对于水平地面，单位反射辐亮充度为： 对于水平地面，单位反射辐亮充度为：
L = F0 cosθ s cosθv R(θ s,ϕs ;θv,ϕv )
对足够大的地面面积， 对足够大的地面面积，反射辐亮度仅与二向反射率有 如地面为朗伯体且水平， 关；如地面为朗伯体且水平，则任何观测方向上其反 射辐亮度均相等。对于传感器来说， 射辐亮度均相等。对于传感器来说，一般地物均充满 整个瞬时视场，故反射亮度与观测方向不是余弦关系， 整个瞬时视场，故反射亮度与观测方向不是余弦关系， 其反射辐亮度为： 其反射辐亮度为：
辐射定标实习报告
辐射定标一、实验数据源描述：（一）实验数据来源：地点：青岛区域时间：2009年7月15号搭载传感器：Landsat TM 5传感器特征：含有7个波段，其中第6波段为热红外波段，空间分辨率是120m ，其他波段空间分辨率是30m 。
辐射分辨率是8bit ，重访周期是16天，幅宽是185km 。
TM 各波段及其作用：TM Band 1蓝绿波段 0.45-0.52 30 用于水体穿透,分辨土壤植被 Band 2 绿色波段 0.52-0.60 30 分辨植被Band 3 红色波段 0.63-0.69 30 处于叶绿素吸收区域,用于观测道路/裸露土壤/植被种类效果很好Band 4 近红外波段 0.76-0.90 30 用于估算生物量, 尽管这个波段可以从植被中区分出水体，分辨潮湿土壤，但对道路辨认效果不如TM3Band 5 中红外波段 1.55-1.75 30 用于分辨道路/裸露土壤/水, 它在不同植被之间有好的对比度,并且有较好的穿透大气、云雾的能力。
Band 6 热红外波段 10.40-12.50 120 感应发出热辐射的目标。
Band 7中红外波段 2.08-2.35 30 对于岩石/矿物的分辨很有用, 也可用于辨识植被覆盖和湿润土壤。
（二）青岛 - 自然环境1、地质地貌：青岛为海滨丘陵城市，地势东高西低，南北两侧隆起，中间低凹，其中山地约占全市总面积的15.5％、丘陵占25.1％、平原占37.7％、洼地占21.7％。
全市海岸分为岬湾相间的山基岩岸、山地港湾泥质粉砂岸及基岩砂砾质海岸等3种基本类型。
浅海海底则有水下浅滩、现代水下三角洲及海冲蚀平原等。
青岛所处大地构造位置为新华夏隆起带次级构造单元——胶南隆起区东北缘和胶莱凹陷区中南部。
区内缺失整个古生界地层及部分中生界地层，但白垩系青山组火山岩层发育充分，在本市出露十分广泛。
岩浆岩以元古代胶南期月季山式片麻状花岗岩及中生代燕山晚期的艾山式花岗闪长岩和崂山式花岗岩为主。
辐射定标(像元亮度值,辐射亮度亮温)、表观反射率、地表反射率、反照率、比辐射率
环境一号卫星光学数据绝对定标环境一号卫星光学数据的遥感器校正分为绝对定标和相对辐射定标。
对目标作定量的描述，得到目标的辐射绝对值。
要建立传感器测量的数字信号与对应的辐射能量之间的数量关系，即定标系数，在卫星发射前后都要进行。
卫星发射前的绝对定标是在地面实验室或实验场，用传感器观测辐射亮度值已知的标准辐射源以获得定标数据。
卫星发射后，定标数据主要采用敦煌外场测量数据，此值一般在图像头文件信息中可以读取。
以下两表为敦煌场地测定的绝对定标数据。
表HJ 1A/B星绝对辐射定标系数（DN/W⋅m-2⋅sr-1⋅μm-1）利用绝对定标系数将DN值图像转换为辐亮度图像的公式为：L=DN/coe式中coe为绝对定标系数，转换后辐亮度单位为W⋅m-2⋅sr-1⋅μm-1。
由于以上定标系数为敦煌场采用单点法对中等反射率目标（戈壁）测定的结果，因此对于太阳反射光谱波段，建议针对中等反射率地物采用上面提供的绝对辐射定标系数。
对于HJ1B的红外相机，近红外波段绝对定标系数为4.2857，短波红外波段绝对定标系数为18.5579。
定标公式同前。
HJ-1B红外相机热红外通道绝对辐射定标系数为：增益53.473，单位：DN/（W⋅m-2⋅sr-1⋅μm-1）；截距26.965，单位：DN。
利用绝对定标系数将DN值图像转换为辐亮度图像的公式为L=（DN-b）/coe，式中coe为绝对定标系数的增益，b为截距，转换后辐亮度单位为W⋅m-2⋅sr-1⋅μm-1。
HJ1B红外相机中红外波段则条带较为严重，不利于定量化应用。
遥感数字图像遥感数字图像是以数字形式记录的二维遥感信息，即其内容是通过遥感手段获得的，通常是地物不同波段的电磁波谱信息。
其中的像素值称为亮度值（或称为灰度值、DN值）。
遥感概念DN值（Digital Number ）是遥感影像像元亮度值，记录的地物的灰度值。
无单位，是一个整数值，值大小与传感器的辐射分辨率、地物发射率、大气透过率和散射率等有关。
辐射定标原理与方法
辐射定标原理与方法嗨，朋友！今天咱们来聊聊一个超酷的话题——辐射定标。
你可能会想，这是啥玩意儿？听我慢慢道来。
我有个朋友叫小李，他在一家做遥感图像的公司上班。
有一次，我去他公司玩，看到那些花花绿绿的遥感图像，觉得特别神奇。
我就问他：“小李啊，这些图像怎么这么准确呢？看起来就像真的从天上俯瞰大地一样。
”小李笑着跟我说：“这可多亏了辐射定标啊。
”我当时就懵了，辐射定标？这听起来就像是来自外太空的术语。
那到底什么是辐射定标呢？简单来说，辐射定标就像是给测量辐射的仪器做个精准的尺子。
你想啊，假如你要量东西，尺子不准，那量出来的结果肯定不对呀。
辐射定标就是要让测量辐射的仪器知道，多少辐射量对应的是多少数值。
这就好比我们要知道一杯水有多少毫升，得先有个标准的刻度杯一样。
辐射定标有它的原理呢。
想象一下，辐射就像一群调皮的小粒子在到处跑。
仪器要捕捉它们，就像我们用网去捞鱼一样。
但是每个网的大小、疏密不一样，捞到鱼的数量和种类可能就不一样。
仪器也一样，不同的仪器对辐射的感应是不同的。
所以我们得建立一个标准，让所有的仪器都按照这个标准来衡量辐射量。
这就像是让所有捞鱼的网都按照同样的标准制作，这样捞到鱼的数量才有可比性。
那怎么进行辐射定标呢？这里面的方法可不少。
有一种叫实验室定标。
这就像是运动员在训练基地进行特训一样。
在实验室里，我们可以控制各种条件，比如温度、湿度、光照等等。
把仪器放在这样一个稳定的环境里，然后用已知辐射量的光源去照射它。
就好比给运动员一个标准的训练环境和明确的训练目标一样。
仪器这个时候就会产生一个对应的响应值。
通过多次这样的实验，我们就能建立起辐射量和仪器响应值之间的关系。
这关系就像是一把钥匙，能让我们准确地解读仪器测量到的数值。
还有一种定标方法叫现场定标。
这有点像在实际比赛场地对运动员进行最后的调整。
现场定标就是把仪器带到实际测量的地方，比如说在一片田野或者一座山上。
在这个真实的环境里，利用一些已经被精确测量过辐射量的目标物来对仪器进行定标。
遥感科学-第九章-光学传感器定标
中国科学院大学 姜小光
2017年12月20日
主要内容
辐射定标的概念、目的与意义 辐射定标的分类与基本原理 发射前定标（实验室定标、场地定标） 发射后定标（星上定标、场地定标、交叉定标）
辐射定标的概念
辐射定标： 定义：建立遥感器记录的无量纲的数字信号（DN值） 与对应目标的辐射能量之间的定量关系，即得到各波 段的定标系数。
Li=ai*DNi+bi
Li:辐射能量；DNi：灰度值；ai：增益； bi：偏置
遥感数据通常记录的是灰度值（DN），定标就是要将 灰度值转换为物理量，如反射率、辐射亮度的过程。
辐射定标的概念
辐射定标（续）： 相对辐射定标：校正探测器不同像元的响应度的过程 （如均匀性校正、同一地方不同时间图像输出值之间 的相对关系的确定等）
1.发射前定标（续）： （1）实验室定标：
辐射定标分类
积分球： 一个内壁涂有白色漫反射材料的空腔球体，球壁上开 一个或几个窗孔，用作进光孔和放置光接收器件的接 收孔。积分球的内壁是良好的球面，球内壁上涂以理 想的漫反射材料（氧化镁或硫酸钡），氧化镁涂层在 可见光谱范围内的光谱反射比都在99%以上，这样，进 入积分球的光经过内壁涂层多次反射，在内壁上形成 均匀照度。
辐射定标分类
1.发射前定标（续）： （2）外场地定标（续） 地面定标： 在卫星发射前，在室外直接使用太阳光替代实验室光 源，可以避免室内光源在光谱形状、量级上与实际光 源存在差异，同时也是对实验室定标的有效检验。 航空飞行定标： 研制相应的机载模拟器或测试载荷，在载荷上天前， 利用航空平台进行载荷的定标。开展航空飞行测试， 除了可在实际工况下验证载荷的工程指标及实际应用 效能，还有利于改进数据处理方法，发展更有针对性 的参数反演算法和应用模型。
ENVI上机指导(二)
ENVI 上机指导（二）——基于6S 模型的遥感数据大气纠正图像预处理流程分为几个步骤：图像的配准、重采样、定标以及大气纠正。
其中，图像配准是做所有工作的前提，是图像的几何纠正。
在进行定性或者定量遥感时都要进行图像配准来确定我们所要研究的目标。
在进行完图像配准（几何纠正）之后，为了使得输出图像的配置与输入图像相对应，因此要进行重采样。
定标以及大气纠正则是进行辐射量纠正的重要的过程，是进行定量遥感不可少的步骤。
辐射定标和大气纠正都属于图像的预处理，辐射定标的目的是把图像上的DN 值转为辐亮度，大气纠正的目的是消除或减少大气对图像的干扰，得到地物的真实反射率。
一、 数据：使用的数据为广东省汕头市的ETM+影像，成像时间为2001年11月22日，2:28:18.000（格林威治时间）。
波段列表：TM1：0.45~0.52 (um) 蓝波段 30m TM2：0.52~0.60 (um) 绿波段 30m TM3：0.63~0.69 (um) 红波段 30m TM4：0.76~0.90 (um) 近红外波段 30m TM5：1.55~1.75 (um) 短波红外波段 30m TM7：2.08~2.35 (um) 短波红外波段 30m二、 辐射定标1. 首先对图像进行辐射定标，将图像的DN 值转化为辐亮度，offset DN gain L +=* （1）其中，λλλλλλmin max min max min max DN DN L L e fullDNrang L L gain −−=−=，，λmin Loffset =每个角标中含有λ的参数表示波段不同则取值不同，具体参数可从卫星影像的头文件中得到。
L 是某个波段光谱辐射亮度；gain 为增量校正系数，offset 为校正偏差量，DN 是图像灰度值，DNmax和DNmin为遥感器最大和最小灰度值Lmax， Lmin分别为最大和最小灰度值所相应的辐射亮度。
2.在ENVI中进行辐射定标——波段运算1）选择 Basic ToolsÆBand Math.2）在标签为 “Enter an expression:” 的文本框内，输入数学表达式。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
一总述
1 遥感图像处理的目的
遥感的目的是为了获得地物的几何属性和物理属性.但是由于受到大气,目标,传感器等诸多因素的影响,原始的遥感影像中除了有目标地物的信息以外还包含有大气,传感器的运行状态等信息,如果我们只是利用原始的遥感影像,将不能提取出所感兴趣的有效信息, 所以为了实现遥感的最终目的,提取所需的信息,我们必须对遥感影像进行处理.
2 ENVI简介
目前已经开发了一些进行遥感图像处理的软件,例如ENVI,PCI,ERDAS等.现在就简单介绍一下ENVI.
ENVI是由美国RSI公司开发的一套功能齐全的遥感图像处理系统，是处理、分析并显示多光谱数据、高光谱数据和雷达数据的高级工具。
其完全是由IDL开发，方便灵活，可扩展性强,并可用IDL进行二次开发。
现在最高版是4.7版本的.
我们来大概熟悉一下
ENVI的主菜单:
可以看出ENVI的主菜单中主要有以下一些工具: 基本工具,分类,空间变换,滤波,波谱工具,制图工具,矢量工具,地形分析,雷达工具
来看一下主菜单中的FILE菜单, 通过选择Open Image File可以打开ENVI 图像文件或其它已知格式的二进制图像文件。
ENVI 自动地识别和读取下列类型的文件：TIFF、GeoTIFF、GIF、JPEG、BMP、SRF、HDF、PDS、MAS-50、NLAPS、RADARSAT 和A VHRR 。
数据仍保留它原有格式，必要的信息从数据头文件中读取。
ENVI也直接读取其它几种文件类型（参见“O pen External File”）。
注意:
若你得到“File does not appear to be a valid Radarsat file” 这样一个错误消息，使用File > Open External File 来选择正确的数据类型。
当ENVI 第一次打开一个文件，它需要关于文件特征的特定信息。
通常，这些信息存储在与图像文件同名的一个独立的文本头文件，但是文件扩展名为.hdr 。
若文件打开时没有找到ENVI头文件，你必须在Header Information 对话框中输入一些基本的参数.
另外一些数据格式没有.hdr 文件也能自动打开。
这些格式包括：TIFF、GeoTIFF、GIF、JPEG、BMP、SRF、HDF、PDS、MAS-50、NLAPS、RADARSAT 和A VHRR 。
(ENVI 头文件中含有丰富的信息,例如:
ENVI
description = {
Create New File Result [Tue Oct 19 15:47:45 2004]}
samples = 2000
lines = 2000 (图像的大小)
bands = 7
header offset = 0
file type = ENVI Standard
data type = 4
interleave = bsq (存储方式)
sensor type = Unknown
byte order = 0
x start = 3252
y start = 4828
(坐标投影情况)
map info = {Transverse Mercator, 1.0000, 1.0000, 428019.2564, 4457909.5476,
3.0000000000e+001, 3.0000000000e+001, , units=Meters}
projection info = {3, 0.0, 0.0, 0.000000, 67036062.030306, 500000.0, 0.0, 1.000000, Transverse Mercator, units=Meters}
wavelength units = nm (图像波段的情况)
wavelength = {
482.000000, 564.000000, 659.000000, 837.000000, 1649.000000,
10000.000000, 2219.000000}
关于ENVI的一些基本知识,我们就介绍到这里,如果想了解更多的,请参考用户手册和ENVI 中的HELP.
下面是关于ENVI的一些具体应用.
二图像的预处理
拿到一幅原始图像,我们先要进行辐射定标,目的是把图像上的DN值转为辐亮度或者是反射率.另外通过大气纠正,我们可以消除一些大气的干扰.
1 可以采用简单的波段运算
例如: 我们把2004年四月份的TM图像第3波段的DN值转化为表观反射率。
运用公式：radiance=（（lmax-lmin）/（qcalmax-qcmin））*（qcal-qcamin）+lmin
qcal-----DN
high gain ：b3-------lmax=264，lmin=-1.17；
ρ=π*L*d2/（ESUN*cos（θ））
L----radiance θ=42.43。
d=0.9909 π=3.1415
ESUN：b3：1554
我们在ENVI中运用band math进行相关的处理:
对于波段三: L3=1.039880*b3-1.17
ρ3=π*L3*d2/（ESUN*cos（θ））=3.1415* (1.039880*b3-1.17)*
0.9909^2/(1554*(cos42.43))= 3.1415* (1.039880*b3-1.17)* 0.9909^2/(1554*0.7381)
运算之前
实施运算
Basic tools > band math
输入运算式:
运算之后:
对于其它波段只要知道相关的参数,可以用同样的方法作简单的定标. 其中关于ETM图像的一些参数如下:
2用ENVI中的工具
(1) spectral> FLAASH
FLAASH相当于ENVI的插件,需要格外安装,其是根据类似于MODTRAN4的方法进行大气纠正,并得到地面的反射率.
(2)简单的经验线性法(不考虑交叉辐射)
此种方法比较简单,精度一般.需要已知图像上各波段地表反射率的最高值与最低值,通过建立反射率与图像中所选取的感兴趣区域的DN值的线性关系来进行简单的大气纠正.
spectral> preprocessing> Calibration Utilities>Empirical Line
compute factors and calibrate(计算参数,进行校准)
calibrate using existing factors(用已有的参数进行校准)
我们选择计算参数进行校准
过程如下:
选择compute factors and calibrate
,
点击上面的import spectra
打开一个ROI,选择上面IMPOT的from roi from input file,选好后点击OK
点击APPL Y
点击下一个import spectra 选择from ascii file
选好后,得到
把图像中亮(暗)的地方与高(低)反射率一一对应,此中RED对应C2(高),GREEN对应C3(低),所以:
选择OK,如果出现出错的信息,不用管,继续,下图为所建立的线性关系.
上一条线为截距,下一条线为斜率.
除了上面的方法还可以使用其它的纠正模型进行大气定标和大气纠正.
三NDVI的计算
(1)使用BAND MA TH
将经过了辐射定标(或者是定标+大气纠正)的图像,运用工具中的band math
采用公式：NDVI=（b1-b2）/（b1+b2）b1---近红外的反射率b2---红光的反射率
例如:
(2)使用IDL
IDL简介
IDL(interactive data language)交互式数据语言是进行数据分析,可视化表现和应用开发的软件工具.其特性包括: 高级图像处理能力,交互式二维和三维图形技术,面对对象的编程方式,opengl图形加速,量化可视化表现,集成数学与统计学算法,灵活的数据输入输出方式,
跨平台图形用户界面工具包,连接ODBC兼容数据库存及多种程序连接工具等.
以下是idl的主界面:
IDL程序的特点:
(a)分隔符为”,”,而非空格; (b) 不分大小写; (c)变量无需事先申明
IDL的符号(需注意的)
●$
–作为一行的第一个字符时，返回到操作系统下，如：$ dir
–作为一行的最后一个字符时，相当于一行未写完，换行。
●；后面是注释
●@ 批作业如：@test
关于IDL的知识,大家可以找相应的资料看
下面是简单利用一下IDL做波段运算
A 运用IDL编写波段运算的公式,再在BAND MATH中调用
例如:
选择ENVI中的BAND MATH,
选择相应的波段进行运算.
B 直接运用IDL编程,计算NDVI
例如:。