本发明属于光谱学、遥感和图像识别的交叉技术领域,更具体地,涉及一种图谱关联的分谱段估计目标亮温的方法。
背景技术:
使用遥感红外数据进行温度的反演以获得目标的亮温对于地面监测有非常重要的意义,尤其是在对水文、生态、环境等方面。
目前利用遥感技术对地表温度信息的提取主要是利用热红外波段,但热红外波段数据受大气、背景环境影响较大,反演精度较低。随着科技的不断进步,已有学者提出了一些很好的利用中红外或短波红外来识别和反演温度的方法:程洁等提出利用中红外波段提取土壤发射率的算法,利用提取出土壤的发射率反演地表温度;短波红外谱段(1.3~3.0μm),对高温目标有一定的敏感性,对森林火、活火山等高温目标的识别精度较高,能够清楚的显示火点的位置和大小,因此,短波红外谱段常被用于捕捉高温信息。
以上地面遥感技术对于各类目标,均采用相同波段的光谱进行温度的反演,所得到的亮温精度不高。而且,以上地面遥感技术获得是大区域的地表遥感数据,反演得到的是大区域面积的平均温度,难以得到各小局部区域的精确反演结果。
技术实现要素:
针对现有技术的缺陷和改进需求,本发明提供了一种图谱关联的分谱段估计目标亮温的方法,其目的在于,精确估计目标各局部区域的亮温。
为实现上述目的,本发明提供了一种图谱关联的分谱段估计目标亮温的方法,包括:
(1)凝视测量目标的红外宽光谱及对应的长波红外图像;
(2)根据长波红外图像识别目标的类型,若为小目标,则转入步骤(5);若为大目标,则转入步骤(6);
(3)根据红外宽光谱的形态及长波红外图像识别待测对象的类型,并根据识别结果确定进行温度反演的目标谱段;
不同类型待测对象,其红外宽光谱的形态会有所不同;当待测对象为高温物体时,在光谱的短波波段和/或中波波段会出现波峰,使用短波波段对这类待测对象进行温度反演的精度更高;待测对象反射太阳光的能量强时,在光谱的短波波段和/或中波波段也会出现波峰,此时应使用长波波段或全波段进行温度的反演;中、低温待测对象在短波波段和中波波段均不会出现波峰,对这类目标,应使用长波波段进行温度反演;本发明根据红外宽光谱的形态及长波红外图像能够准确识别待测对象的类型,针对不同类型的待测对象使用不同的谱段进行温度反演,能够有效提高亮温估计的准确度;
(4)通过温度反演分别得到目标谱段内每一个波长所对应的亮温,并获得所有波长对应的亮温的平均值作为待测对象的亮温;
(5)将目标作为待测对象并执行步骤(3)~(4),以估计目标的亮温,处理结束;
(6)将目标划分为多个局部区域,分别以各局部区域为待测对象并执行步骤(3)~(4),从而估计目标各局部区域的亮温,处理结束。
本发明所提供的图谱关联的分谱段估计目标亮温的方法,对于小目标,直接估计其亮温,对于大目标,则先将其划分为多个局部区域并分别估计各局部区域的亮温,由此能够获得大目标的亮温分布;由于大目标各局部区域的亮温可能不一致,相比于传统的只能获得大区域面积平均温度的亮温估计方法,本发明具有更重要的实际意义。
进一步地,步骤(2)中,根据长波红外图像识别目标的类型,包括:
获得目标在红外图像中所占像素数的比值,若比值大于预设的尺寸阈值,则将目标识别为大目标;否则,将目标识别为小目标。
进一步地,步骤(3)包括:
获得红外宽光谱对应的波段b0,并分别获得波段b0中的短波波段和中波波段;
若短波波段和中波波段均没有波峰,则识别待测对象为中温对象或低温对象,并确定波段b0或者波段b0中的长波波段为目标谱段;
若短波波段和/或中波波段有波峰,且待测对象在长波红外图像中为亮目标,则识别待测对象为高温对象,并确定出现波峰的波峰波段为目标谱段;
若短波波段和/或中波波段有波峰,且待测对象在长波红外图像中不是亮目标,则识别待测对象为反射太阳光能量较多的对象,并确定波段b0中的长波波段为目标谱段。
更进一步地,波峰波段的获取方式为:
根据红外宽光谱计算中间波段内的平均dn值
将波段范围b1向其邻域扩大一倍,以得到波峰波段为:
其中,中间波段为波段范围b0中的短波波段或中波波段,λi和λj分别为波段范围b1的下界和上界,且λi<λj,λi′和λj′分别为波段b的下界和上界,且λi′<λj′,λk为中间变量。
更进一步地,判断待测对象在长波红外图像中是否为亮目标的方式为:
遍历长波红外图像以获得灰度直方图;
根据灰度直方图分别获得灰度值最高的前α%像素的平均灰度值
若
其中,α和β均为预设的百分比阈值,vth为预设的灰度值阈值。
进一步地,步骤(4)中,通过温度反演分别得到目标谱段内每一个波长所对应的亮温,包括:
(41)对于目标谱段内的每一个波长λ,分别获得黑体在温度t1和t2下与波长λ相对应的dn值s(λ,t1)和s(λ,t2),以计算与波长λ相对应的仪器响应函数k(λ)和仪器辐射偏置b(λ);
(42)根据红外宽光谱获得待测对象与波长λ相对应的dn值s(λ,t),并根据仪器响应函数k(λ)、仪器辐射偏置b(λ)以及dn值s(λ,t)计算校正后的辐射亮度g(λ,t);
(43)从辐射亮度g(λ,t)中去除大气程辐射gpath(λ)以得到待测对象在像方的辐射亮度g′(λ,t′);
(44)根据辐射亮度g′(λ,t′)和大气透过率τ(λ)得到待测对象在物方的辐射亮度g″(λ,t″),并根据辐射亮度g″(λ,t″)计算待测对象的亮温t″;
其中,大气程辐射gpath(λ)和大气透过率τ(λ)均与波长λ相对应。
总体而言,通过本发明所构思的以上技术方案,能够取得以下有益效果:
(1)本发明所提供的图谱关联的分谱段估计目标亮温的方法,同时获得目标的红外宽光谱和长波红外图像,能够根据红外宽光谱的形态及长波红外图像能够准确识别待测对象的类型。
(2)本发明所提供的图谱关联的分谱段估计目标亮温的方法,针对不同类型的待测对象使用不同的谱段进行温度反演,能够有效提高亮温估计的准确度。
(3)本发明所提供的图谱关联的分谱段估计目标亮温的方法,对于小目标,直接分谱段估计其亮温;对于大目标,则先将其划分为多个局部区域并分别分谱段估计各局部区域的亮温,由此能够获得大目标的亮温分布;由于大目标各局部区域的亮温可能不一致,相比于传统的只能获得大区域面积平均温度的亮温估计方法,本发明具有更重要的实际意义。
附图说明
图1为本发明实施例提供的图谱关联的分谱段估计目标亮温的方法流程图;
图2(a)为本发明实施例提供的凝视测量海水亮带的红外宽光谱和长波红外图像的示意图;
图2(b)为本发明实施例提供的凝视测量防浪堤的红外宽光谱和长波红外图像的示意图;
图3(a)为本发明实施例提供的海水亮带长波红外图像;
图3(b)为本发明实施例提供的防浪堤长波红外图像;
图4为本发明实施例提供的海水亮带和防浪堤的红外宽光谱;
图5为本发明实施例提供的飞机各部位的长波红外图像;(a)~(d)分别为飞机机头、机翼、尾喷管及尾焰的长波红外图像;
图6为本发明实施例提供的飞机各部位的红外宽光谱;
图7为本发明实施例提供的飞机各部位的光谱辐射校正曲线;(a)~(c)分别为飞机机体、机翼及尾喷管在6.5~9.5um波段内的关铺辐射校正曲线;(d)为飞机尾焰在3.2~4.2um波段内的光谱辐射校正曲线;
图8为本发明实施例提供的飞机各部位的亮温估计结果;(a)~(c)分别为飞机机体、机翼及尾喷管的亮温估计结果;(d)为飞机尾焰在3.2~4.2um波段内的亮温估计结果;(e)为飞机尾焰在6.5~9.5um波段内的亮温估计结果;
图9(a)为本发明实施例提供的防浪堤的光谱辐射校正曲线;
图9(b)为本发明实施例提供的防浪堤的亮温估计结果;
图10(a)为本发明实施例提供的海水亮带的光谱辐射校正曲线;
图10(b)为本发明实施例提供的海水亮带在9~11um波段内的亮温估计结果;
图10(c)为本发明实施例提供的海水亮带在3.2~4.2um波段内的亮温估计结果。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。此外,下面所描述的本发明各个实施方式中所涉及到的技术特征只要彼此之间未构成冲突就可以相互组合。
本发明提供了一种图谱关联的分谱段估计目标亮温的方法,其整体思路在于:同时获得目标的红外宽光谱及对应的长波红外图像;根据长波红外图像识别目标类型;对于小目标,则直接分谱段估计其亮温;对于大目标,则将其划分为多个局部区域,并分别分谱段估计各局部区域的亮温,以得到目标的亮温分布;在分谱段估计亮温的过程中,结合红外宽光谱的形态和长波红外图像识别待测对象的类型,并根据识别结果确定更为具体的目标谱段,在目标谱段内进行温度反演以提高亮温估计的准确度。
本发明所提供的图谱关联的分谱段估计目标亮温的方法,如图1所示,包括:
(1)凝视测量目标的红外宽光谱及对应的长波红外图像;
(2)根据长波红外图像识别目标的类型,若为小目标,则转入步骤(5);若为大目标,则转入步骤(6);
在一个可选的实施方式中,步骤(2)中,根据长波红外图像识别目标的类型,包括:
获得目标在红外图像中所占像素数的比值,若比值大于预设的尺寸阈值,则将目标识别为大目标;否则,将目标识别为小目标;
(3)根据红外宽光谱的形态及长波红外图像识别待测对象的类型,并根据识别结果确定进行温度反演的目标谱段;
根据维恩位移定律,黑体辐射光谱分布曲线极大值对应的波长λ随黑体温度的升高,向短波长方向移动;
不同类型的待测对象,其红外宽光谱的形态会有所不同;当待测对象为高温物体时,在光谱的短波波段和/或中波波段会出现波峰,使用短波波段对这类待测对象进行温度反演的精度更高;待测对象反射太阳光的能量强时,在光谱的短波波段和/或中波波段也会出现波峰,此时应使用长波波段或全波段进行温度的反演;中、低温待测对象在短波波段和中波波段均不会出现波峰,对这类目标,应使用长波波段进行温度反演;本发明根据红外宽光谱的形态及长波红外图像能够准确识别待测对象的类型,针对不同类型的待测对象使用不同的谱段进行温度反演,能够有效提高亮温估计的准确度;
在一个可选的实施方式中,步骤(3)具体包括:
获得红外宽光谱对应的波段b0,并分别获得波段b0中的短波波段和中波波段;
若短波波段和中波波段均没有波峰,则识别待测对象为中温对象或低温对象,并确定波段b0或者波段b0中的长波波段为目标谱段;
若短波波段和/或所述中波波段有波峰,且待测对象在长波红外图像中为亮目标,则识别待测对象为高温对象,并确定出现波峰的波峰波段为目标谱段;
若短波波段和/或所述中波波段有波峰,且待测对象在长波红外图像中不是亮目标,则识别待测对象为反射太阳光能量较多的对象,并确定波段b0中的长波波段为目标谱段;
在本实施例中,波峰波段的获取方式为:
根据红外宽光谱计算中间波段内的平均dn值
将波段范围b1向其邻域扩大一倍,以得到波峰波段为:
其中,中间波段为波段范围b1中的短波波段或中波波段,λi和λj分别为波段范围b1的下界和上界,且λi<λj,λi′和λj′分别为波段b的下界和上界,且λi′<λj′,λk为中间变量;
判断待测对象在长波红外图像中是否为亮目标的方式为:
遍历长波红外图像以获得灰度直方图;
根据灰度直方图分别获得灰度值最高的前α%像素的平均灰度值
若
其中,α和β均为预设的百分比阈值,vth为预设的灰度值阈值;
(4)通过温度反演分别得到目标谱段内每一个波长所对应的亮温,并获得所有波长对应的亮温的平均值作为待测对象的亮温;
在一个可选的实施方式中,步骤(4)中,通过温度反演分别得到目标谱段内每一个波长所对应的亮温,具体包括:
(41)对于目标谱段内的每一个波长λ,分别获得黑体在温度t1和t2下与波长λ相对应的dn值s(λ,t1)和s(λ,t2),以计算与波长λ相对应的仪器响应函数k(λ)和仪器辐射偏置b(λ);
(42)根据红外宽光谱获得待测对象与波长λ相对应的dn值s(λ,t),并根据仪器响应函数k(λ)、仪器辐射偏置b(λ)以及dn值s(λ,t)计算校正后的辐射亮度g(λ,t);
(43)从辐射亮度g(λ,t)中去除大气程辐射gpath(λ)以得到待测对象在像方的辐射亮度g′(λ,t′);
(44)根据辐射亮度g′(λ,t′)和大气透过率τ(λ)得到待测对象在物方的辐射亮度g″(λ,t″),并根据辐射亮度g″(λ,t″)计算待测对象的亮温t″;
其中,大气程辐射gpath(λ)和大气透过率τ(λ)均与波长λ相对应;
(5)将目标作为待测对象并执行步骤(3)~(4),以估计目标的亮温,处理结束;
(6)将目标划分为多个局部区域,分别以各局部区域为待测对象并执行步骤(3)~(4),从而估计目标各局部区域的亮温,处理结束。
本发明所提供的图谱关联的分谱段估计目标亮温的方法,对于小目标,直接估计其亮温,对于大目标,则先将其划分为多个局部区域并分别估计各局部区域的亮温,由此能够获得大目标的亮温分布;由于大目标各局部区域的亮温可能不一致,相比于传统的只能获得大区域面积平均温度的亮温估计方法,本发明具有更重要的实际意义。
应用实例
使用华中科技大学张天序教授课题组研制的图谱协同遥感探测系统跟踪和凝视相结合的方式测量多组飞机不同部位、海水亮带与背景1.7~20um红外宽光谱数据和8~12um长波红外图像,以及面源黑体光谱数据。在亮温估计是,将实测1.7~20um和2~12um的宽光谱分为短波波段(2~3um)、中波波段(3~5um)和长波波段(8~12um),针对不同类型(如高温、低温、反射太阳光强)目标温度的反演使用不同谱段光谱。
图谱协同遥感探测系统凝视测量海水亮带与防浪堤的红外宽光谱和长波红外图像,数据采集示意图分别如图2(a)和图2(b)所示,两类各采集光谱样本100个,红外图像各100帧,每个光谱样本在1.7~20um波长上的dn值数量n=1545,所采集到的海水亮带和防浪堤的长波红外图像分别如图3(a)和图3(b)所示,所采集到的海水亮带和防浪堤的红外宽光谱如图4所示;对高温目标,测量机场飞机尾焰光谱样本100个,红外图像100帧,每个光谱样本在2~12um波长上的dn值数量n=1042;飞机机头、机翼、尾喷管及尾焰的长波红外图像分别如图5(a)~图5(d)所示,飞机各部位的红外宽光谱如图6所示。
从图3和图5可知,对实测光谱对应的红外图像,遍历图像获得灰度直方图,计算灰度直方图中高灰度值前α%且高于图像光斑灰度值的像素数n和低灰度值前β%的像素数m,在本应用实例中,α=10,β=10,分别计算这n和m个像素数的平均灰度值
从图4和图6可知,红外宽光谱方面,对于海面亮带,反射太阳光能量造成的波峰出现在光谱短波波段和中波波段,短波波段的光谱dn值的平均值
band=λi-k,j+k
对于海面上的防浪堤,在光谱短波波段没有波峰,且在长波红外图像中不是亮目标,因此判定防浪堤为中、低温物体。
飞机尾焰在红外图像中很亮,且在红外光谱中波波段有波峰出现,中波波段的光谱dn得到平均值
根据上述分析,可以识别各对象(即目标或目标的局部区域)的类型,对不同目标或同一目标的不同部件选择合适谱段范围进行温度的反演。由于飞机尾焰是高温目标且在短波波段有波峰出现,所以选择用波峰出现的波峰波段进行温度反演;海水亮带为低温物体但在短波波段有波峰出现,短波反射太阳光能量多,所以选择长波波段进行温度的反演;防浪堤属于中、低温目标,且反射太阳光能量少,所以选择长波红外波段进行温度反演。
上述目标与背景的表观温度是由实测目标红外宽光谱经辐射定标和大气校正后得到的辐射亮度谱,再经普朗克公式计算所得。辐射校正使用两点校正公式,大气校正用的大气透过率和大气程辐射由modtran根据当时测量条件计算所得。测量条件如表1所示。
表1测量条件
利用实测黑体光谱数据对目标与背景光谱进行辐射校正,两点校正公式为s(λ,t1)=k(λ)h(λ,t1)+b(λ)
s(λ,t2)=k(λ)h(λ,t2)+b(λ)
其中,s(λ,t1)和s(λ,t2)分别表示两个不同温度的实测黑体dn值,k(λ)表示仪器响应函数,h(λ,t1)和h(λ,t2)分别表示这两个不同温度的普朗克数据,b(λ)表示仪器辐射偏置,
对观测到的目标dn值s(λ,t)有,
s(λ,t)=so(λ,t)+sp(λ,t)
so(λ,t)=τ(λ)*starget(λ,t)
其中so(λ,t)是实际目标进入到图谱协同设备的dn值,starget(λ,t)是目标在物方的dn值,τ(λ)为当时测量条件下的大气透过率,计算得到校正后的辐射亮度g(λ,t)
减去大气程辐射gpath(λ)后的目标在像方的辐射亮度为g′(λ,t′),其像方亮温为t′。
g′(λ,t′)=g(λ,t)-gpath(λ)
g′(λ,t′)=τ(λ)*g″(λ,t″)
g′(λ,t′)除大气透过率τ(λ)后得到目标在物方辐射亮度g″(λ,t″),经普朗克公式计算可得目标物方亮温为t″;
其中,c1为普朗克第一辐射常数,c1=3.74l8×10-16(w*m2);c2为普朗克第二辐射常数c2=1.438786×10-2(m*k)。
飞机机体、机翼及尾喷管在6.5~9.5um波段内以及飞机尾焰在3.2~4.2um波段内的光谱辐射校正曲线分别如图7(a)、7(b)、7(c)和7(d)所示;对应波段内,飞机机体、机翼及尾喷管的亮温估计结果分别如图8(a)~(c)所示,飞机尾焰在3.2~4.2um波段和6.5~9.5um波段内的亮温估计结果分别如图8(d)和图8(e)所示。根据图8所示亮温估计结果可知,如果使用长波波段来估计高温目标(尾焰)的亮温误差较短、中波有波峰的波段大。从图8还可得知,对于飞机的不同部位,其亮温不同,因此通过测量同一目标的不同局部小区域可以得到该目标的温度分布。
防浪堤的光谱辐射校正曲线和亮温估计结果分别如图9(a)和图9(b)所示;海水亮带的光谱辐射校正曲线、9~11um波段内的亮温估计结果以及3.2~4.2um波段内的亮温估计结果分别如图10(a)~图10(c)所示。根据图10所示的亮温估计结果可知,如果使用短、中波波段来估计海水亮带的亮温,误差较长波波段大。
总的来说,根据本发明所提供的图谱关联的分谱段估计目标亮温的方法,通过识别待测对象类型,并对不同类型的目标使用不同波段进行亮温的估计能够有效提高亮温估计的准确度。
本领域的技术人员容易理解,以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。