基于辐射传递的大视场红外探测图像仿真方法技术

一种基于辐射传递的大视场红外探测图像仿真方法，方法的输入涉及两类图像，第一类为场景灰度图像，该类图像为大视场图像，对应的谱段与待仿真谱段相似但不相同；第二类为辐射亮度测量图像，该类图像为场景的多幅小视场图像，对应的谱段与待仿真谱段相同。一、将大视场场景灰度图像进行几何分辨率调整，调整至要求的分辨率；二、建立大视场场景灰度图像中的灰度与多幅小视场辐射亮度测量图像中的强度之间的线性辐射传递关系，将大视场场景灰度图像转换为辐射强度图像；三、基于点目标对应的探测信噪比指标计算噪声均方差，添加对应强度的高斯噪声完成噪声仿真；四、随机添加一定辐射强度的点目标并对辐射图像量化完成图像仿真。

【技术实现步骤摘要】
基于辐射传递的大视场红外探测图像仿真方法
本专利技术涉及一种基于辐射传递的大视场红外探测图像仿真方法，属于目标检测

技术介绍
云背景中红外点目标检测算法的设计、测试都离不开图像源，针对场景中的中层云、高层卷云或冰晶云等不同类型云场对目标检测算法性能的测试必不可少，因此大视场红外探测图像的快速仿真具有一定应用需求。目前，基于三维辐射传输计算的探测器辐射图像仿真是红外探测图像仿真的高精度方法，该方法通过逐次散射迭代进行多次散射源函数计算，一次性计算整个物理场的源函数，而后对3D物理场进行严格的三维辐射传输解算。但该方法与谱段以及云场景中的物质成分严格相关，计算速度较慢。
技术实现思路
本专利技术提出了基于辐射传递的大视场红外探测图像仿真方法，以避免大视场图像仿真复杂三维辐射传输中计算速度较慢的问题，为此本专利技术采用如下的技术方案：一种基于辐射传递的大视场红外探测图像仿真方法，包括：步骤一、以谱段Δλ作为图像仿真谱段，统计Δλ谱段对应的复杂云背景条件下的多幅小视场辐射亮度测量图像中的辐射亮度最大值Lmax和最小值Lmin；步骤二、选择相近谱段ΔλA的复杂云背景条件下的大视场灰度图像A，由于亮度为场景的固有属性，在大视场中场景辐射亮度的最大值和最小值同样分别为Lmax和Lmin；对图像进行双线性插值至仿真要求的分辨率d(m)，得到图像G；将图像G中像元对应的辐射亮度值转化为辐射强度值，具体方式为：I(W/sr)＝L(W/(sr·m2))×d2(m2)其中，I为辐射强度；步骤三、大视场灰度图像G中像素灰度范围为Gmin～Gmax，对应的强度范围为Imin～Imax，灰度与辐射强度之间的线性系数R为：其中，Imax和Imin分别为Lmax、Lmin对应的辐射强度；步骤四、对大视场灰度图像G中的像素灰度进行辐射强度赋值，具体方式为：F(x,y)＝R·[G(x,y)-Gmin]+Imin其中，G(x,y)为大场景图像中(x,y)处的像素灰度，F(x,y)为对应的辐射强度。步骤五、目标的探测性能一般用辐射强度IT时的探测信噪比SNRT描述，即为SNRT＝IT/σ，其中σ为探测器时间噪声均方差，噪声满足高斯分布；IT为目标辐射强度，则噪声均方差可描述为：σ＝IT/SNRT则添加噪声的图像表示为：FN(x,y)＝F(x,y)+N(x,y)其中，FN(x,y)为添加噪声后辐射强度分布图，N(x,y)为0均值，均方差为σ的高斯噪声；步骤六、根据目标辐射强度与数量要求，在添加噪声的辐射强度分布图像中随机位置处添加特定能量ITAR的点形态目标，则添加目标的图像表示为：FT(x,y)＝FN(x,y)+ITARδ(x-x′,y-y′)其中，当x＝x′,y＝y′时，δ(x-x′,y-y′)＝1；否则，δ(x-x′,y-y′)＝0；FT(x,y)为添加目标后辐射强度分布图，(x′,y′)为随机添加目标的位置，ITAR为目标的强度；步骤七、根据量化策略对添加目标后的辐射强度分布图进行量化，完成对图像的仿真。量化策略可通过辐射强度与辐射量化步长IΔ进行比值得到。具体来说，采用辐射量化步长IΔ进行量化仿真方法为：其中，DN(x,y)为输出图像，fix(·)为向最近的整数取整操作。所述大视场灰度图像的对应谱段与待仿真谱段相似。基于辐射传递的大视场红外探测图像仿真方法，包括；第一辐射强度分布图像生成模块，用于利用大视场灰度图像及多幅小视场辐射亮度测量图像获得待仿真谱段对应的第一辐射强度分布图像；第二辐射强度分布图像生成模块，用于根据目标的探测信噪比添加噪声到第一辐射强度仿真图像中，得到第二辐射强度分布图像；第三辐射强度分布图像生成模块，用于随机添加特定的能量点形态目标到所述第二辐射强度分布图像中，获得第三辐射强度分布图像；灰度仿真图像生成模块，用于对所述第三辐射强度分布图像进行量化，得到灰度仿真图像。本专利技术面向点目标检测算法性能测试应用的图像快速仿真需求，提出了基于辐射传递的大视场红外探测图像快速仿真方法。该方法避免了大视场图像仿真中的复杂三维辐射传输计算，基于小视场辐射亮度测量图像与大视场灰度图像建立辐射传递模型，基于点目标探测信噪比建立探测特性模型，实现了含有点目标的大视场红外探测图像的快速仿真。附图说明图1为本专利技术所述的基于辐射传递的图像仿真方法的流程示意图。图2为本专利技术所述的基于辐射传递的图像仿真方法中步骤11的流程示意图。图3为本专利技术所述的基于辐射传递的图像仿真方法中步骤12的流程示意图。图4(1)为实施例一的具体实施步骤中大视场图像源的示意图。图4(2)为实施例一的具体实施步骤中对大视场图像源进行分辨率调整后的示意图。图5(a)为实施例一中的6幅小视场MODIS图像中的一幅示意图。图5(b)为实施例一中的6幅小视场MODIS图像中的一幅示意图。图5(c)为实施例一中的6幅小视场MODIS图像中的一幅示意图。图5(d)为实施例一中的6幅小视场MODIS图像中的一幅示意图。图5(e)为实施例一中的6幅小视场MODIS图像中的一幅示意图。图5(f)为实施例一中的6幅小视场MODIS图像中的一幅示意图。图6为实施例一中6.535-6.895μm的射强度分布图像F0的示意图。图7为实施例一中第二辐射强度分布图像Fn的示意图。图8为实施例一中第三辐射强度分布图像Ft的示意图。图9为实施例一中量化后的灰度仿真图像DN的示意图。具体实施方式下面将结合本专利技术实施例中的附图，对本专利技术实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅是本专利技术一部分实施例，而不是全部的实施例。基于本专利技术中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其它实施例，都属于本专利技术保护的范围。结合图1所示，本具体实施方式提供的基于辐射传递的大视场红外探测图像仿真方法包括：步骤11、利用大视场灰度图像及多幅小视场辐射亮度测量图像获得待仿真谱段对应的第一辐射强度分布图像。步骤12、根据目标的探测信噪比添加噪声到第一辐射强度仿真图像中，得到第二辐射强度分布图像。步骤13、随机添加特定的能量点形态目标到所述第二辐射强度分布图像中，获得第三辐射强度分布图像。步骤14、对所述第三辐射强度分布图像进行量化，得到灰度仿真图像。可选的，结合图2所示，步骤11可以具体包括：步骤110、统计待仿真谱段对应的复杂云背景下的多幅小视场辐射亮度图像中的辐射亮度最大值和最小值。其中，若待仿真谱段为Δλ，统计Δλ谱段对应的复杂云背景条件下的多幅小视场辐射亮度测量图像中的辐射亮度最大值Lmax和最小值Lmin。步骤111、对所述预定大视场灰度图像A进行双线性插值至仿真要求的分辨率d后得到图像G，并将所述图像G中像元对应的辐射亮度值转化为辐射强度值。其中，大视场灰度图像A根据待仿真谱段确定，即选择与Δλ相近的谱段ΔλA的复杂云背景条件下的大视场灰度图像A，由于辐射亮度值为场景的固有属性，在大视场中场景辐射亮度的最大值和最小值同样分别为Lmax和Lmin；对图像进行双线性插值至仿真要求的分辨率d(单位：m)，得到图像G；将图像G中像元对应的辐射亮度值转化为辐射强度值，具体方式为：I＝L×d2其中，I为辐射强度，单位为(W/sr)；L为辐射亮度值，单位为(W/(sr·m2))。步骤本文档来自技高网...

【技术保护点】
一种基于辐射传递的图像仿真方法，其特征在于，步骤一、以谱段Δλ作为图像仿真谱段，统计Δλ谱段对应的复杂云背景条件下的多幅小视场辐射亮度测量图像中的辐射亮度最大值Lmax和最小值Lmin；步骤二、选择相近谱段ΔλA的复杂云背景条件下的大视场灰度图像A，由于亮度为场景的固有属性，在大视场中场景辐射亮度的最大值和最小值同样分别为Lmax和Lmin；对图像进行双线性插值至仿真要求的分辨率d(m)，得到图像G；将图像G中像元对应的辐射亮度值转化为辐射强度值，具体方式为：I(W/sr)＝L(W/(sr·m

【技术特征摘要】
1.一种基于辐射传递的图像仿真方法，其特征在于，步骤一、以谱段Δλ作为图像仿真谱段，统计Δλ谱段对应的复杂云背景条件下的多幅小视场辐射亮度测量图像中的辐射亮度最大值Lmax和最小值Lmin；步骤二、选择相近谱段ΔλA的复杂云背景条件下的大视场灰度图像A，由于亮度为场景的固有属性，在大视场中场景辐射亮度的最大值和最小值同样分别为Lmax和Lmin；对图像进行双线性插值至仿真要求的分辨率d(m)，得到图像G；将图像G中像元对应的辐射亮度值转化为辐射强度值，具体方式为：I(W/sr)＝L(W/(sr·m2))×d2(m2)其中，I为辐射强度；步骤三、大视场灰度图像G中像素灰度范围为Gmin～Gmax，对应的强度范围为Imin～Imax，灰度与辐射强度之间的线性系数R为：其中，Imax和Imin分别为Lmax、Lmin对应的辐射强度；步骤四、对大视场灰度图像G中的像素灰度进行辐射强度赋值，具体方式为：F(x,y)＝R·[G(x,y)-Gmin]+Imin其中，G(x,y)为大场景图像中(x,y)处的像素灰度，F(x,y)为对应的辐射强度。步骤五、目标的探测性能一般用辐射强度IT时的探测信噪比SNRT描述，即为SNRT＝IT/σ，其中σ为探测器时间噪声均方差，噪声满足高斯分布；IT为目标辐射强度，则噪声均方差可描述为：σ＝IT/SNRT则添加噪声的图像表示为：FN(x,y)＝F(x,y)+N(x,y)其中，FN(x,y)为添加噪声后辐射强度分布图，N(x,y)为0均值，均方差为σ的高斯噪...

【专利技术属性】
技术研发人员：侯晴宇，巩晋南，张伟，王志鹏，
申请(专利权)人：哈尔滨工业大学，
类型：发明
国别省市：黑龙江,23