CN115752759B - 一种非均匀性校正单个偏移参数的校正算法 - Google Patents
一种非均匀性校正单个偏移参数的校正算法 Download PDFInfo
- Publication number
- CN115752759B CN115752759B CN202211511725.8A CN202211511725A CN115752759B CN 115752759 B CN115752759 B CN 115752759B CN 202211511725 A CN202211511725 A CN 202211511725A CN 115752759 B CN115752759 B CN 115752759B
- Authority
- CN
- China
- Prior art keywords
- ppij
- correction
- pixel
- coefficient
- storage model
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000012937 correction Methods 0.000 title claims abstract description 129
- 238000003860 storage Methods 0.000 claims abstract description 119
- 230000004044 response Effects 0.000 claims abstract description 95
- 238000000605 extraction Methods 0.000 claims abstract description 61
- 239000002131 composite material Substances 0.000 claims abstract description 27
- 230000005855 radiation Effects 0.000 claims abstract description 27
- 238000000034 method Methods 0.000 claims abstract description 21
- 230000004907 flux Effects 0.000 claims abstract description 18
- 230000008569 process Effects 0.000 claims abstract description 18
- 238000012545 processing Methods 0.000 claims description 9
- 239000000284 extract Substances 0.000 claims description 5
- 125000004122 cyclic group Chemical group 0.000 claims description 2
- 230000003287 optical effect Effects 0.000 claims description 2
- 238000004590 computer program Methods 0.000 claims 2
- 238000001514 detection method Methods 0.000 abstract description 5
- 238000003491 array Methods 0.000 abstract description 2
- 238000012790 confirmation Methods 0.000 abstract description 2
- 238000013461 design Methods 0.000 abstract description 2
- 238000004364 calculation method Methods 0.000 description 13
- 230000006870 function Effects 0.000 description 3
- 238000012795 verification Methods 0.000 description 3
- 101001121408 Homo sapiens L-amino-acid oxidase Proteins 0.000 description 2
- 102100026388 L-amino-acid oxidase Human genes 0.000 description 2
- 230000002776 aggregation Effects 0.000 description 1
- 238000004220 aggregation Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000003331 infrared imaging Methods 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000001052 transient effect Effects 0.000 description 1
Landscapes
- Transforming Light Signals Into Electric Signals (AREA)
Abstract
本发明涉及红外焦平面阵列非均匀性校正设计技术领域,且公开了一种非均匀性校正单个偏移参数的校正算法,利用标准校正参考源提供给非制冷红外焦平面阵列均匀入射辐射通量,测量每个探测单元输出响应,并按照线性假设条件对非制冷红外焦平面阵列的响应曲线进行拟合输出非均匀校正参数,通过构造像元响应模型中单偏移校正参数存储模型、把单偏移参数存储模型聚合生成复合校正存储模型,根据像元和采集数据地址获取复合校正存储模型,从复合校正存储模型中提取单偏移校正参数存储模型,从单偏移校正参数存储模型中提取单个校正偏移参数,整个提取过程经过若干次的识别确认,有效地保证了非制冷红外焦平面阵列上任一探测单元校正结果的准确性。
Description
技术领域
本发明涉及红外焦平面阵列非均匀性校正技术领域,具体为一种非均匀性校正单个偏移参数的校正算法。
背景技术
红外焦平面阵列广泛应用于红外跟踪和制导等领域,理想情况下,同一焦平面的各个像元对一个均匀且稳定的辐射目标的响应是一致的。但是由于空间噪声和瞬态随机噪声的影响,这些像元的响应值之间会有不同程度的差别。这种差别主要是由于各像元之间的非均匀性造成的,可以通过将每个像元的响应曲线校正到一个标准曲线上的方法来减弱这样的差别,也就是所谓的非均匀性校正方法。
改善探测器信号非均匀性一般来说有两种方法:一种改善探测器生产工艺、读出电路设计来和探测器耦合方法,以降低输出信号的非均匀性;另一种是采用后续图像处理电路对探测器在相同辐射强度下不同的响应值进行补偿的方法。相比而言,后一种方案效果显著并更易实现,深受各研究单位的重视。
发明内容
为了解决非制冷型红外焦平面阵列探测器在同一均匀辐射条件下各探测单元响应输出不同的技术问题,本发明提供一种非均匀性校正单个偏移参数的校正算法。
一种非均匀性校正单个偏移参数的校正算法,所述校正算法包括如下步骤:
步骤1,通过采集标准校正参考源在低温和高温情况下非制冷红外焦平面阵列上任一个像元的响应值,计算得到隶属于像元的校正系数;
步骤2,构造像元响应模型中单偏移校正参数存储模型,该单偏移校正参数存储模型包括:校正响应增益系数存储模型={隶属于像元的采集数据地址,校正响应增益系数,隶属于校正响应增益系数的提取参数},校正响应偏置系数存储模型={隶属于像元的采集数据地址,校正响应偏置系数,隶属于校正响应偏置系数的提取参数};
聚合单偏移校正参数存储模型生成复合校正存储模型={隶属于像元的采集数据地址,隶属于复合校正系数的提取参数,校正响应增益系数存储模型,校正响应偏置系数存储模型};
步骤3,获取非制冷红外焦平面阵列上隶属于像元的采集数据地址、与在辐射通量ΦK下的响应值RvK;
调用与像元的采集数据地址对应的复合校正存储模型,在隶属于复合校正系数的提取参数的参与作用下,若采集数据地址识别确认成功,则从复合校正存储模型中提取出单偏移校正参数存储模型;
在隶属于单偏移参数的提取参数的参与作用下,若采集数据地址识别确认成功,则从单偏移校正参数存储模型中提取出校正响应增益系数和校正响应偏置系数,据此输出响应值RvK的校正值ERv(ΦK)。
优选的,所述校正响应增益系数存储模型为存储模型SMG,其生成步骤如下:
设定存储模型SMG的存储系数Gr,然后执行如下的计算过程:
基于存储系数Gr的第一提取参数GR;
基于存储系数Gr的第二提取参数Gh;
基于存储系数Gr的第三提取参数GV;
输出存储模型SMG为{隶属于像元的采集数据地址,校正响应增益系数,第一提取参数GR,第三提取参数GV}。
优选的,所述校正响应偏置系数存储模型为存储模型SMOf,其具体生成步骤如下:
设定存储模型SMOf的存储系数Or,然后执行如下的计算过程:
基于存储系数Orppij的第一提取参数OR;
基于存储系数Orppij的第二提取参数Oh;
基于存储系数Orppij的第三提取参数OV;
输出存储模型SMOfppij为{隶属于像元的采集数据地址,校正响应偏置系数,第一提取参数OR,第三提取参数OV}。
优选的,所述复合校正存储模型为存储模型SM,其具体生成步骤如下:
根据存储模型SMG的第一提取参数GR和存储模型SMOf的第一提取参数OR,计算像元的复合校正偏移参数的第一提取参数R=GhGR+OhOR;
根据存储模型SMG的第三提取参数GV和存储模型SMOf的第三提取参数OV,计算像元的复合校正偏移参数的第三提取参数V=GV+OV;
输出存储模型SM为{隶属于像元的采集数据地址,第一提取参数R,第三提取参数V,存储模型SMG,存储模型SMOf}。
优选的,提取单偏移校正参数存储模型的具体方法如下:
计算像元的复合校正偏移参数的第一识别系数M1,其具体计算结果如下:
M1=e(V,P);
像元的复合校正偏移参数的第二识别系数M2,其具体计算结果如下:
M2=e(2Q,P)e(T,R)e(Q,2xP);
其中,x为像元的校正偏移参数的存储标记;
若第一识别系数M1与第二识别系数M2相等,则提取存储模型SMG和存储模型SMOf。
与现有技术相比,本发明具备以下有益的技术效果:
本发明:利用标准校正参考源提供给非制冷红外焦平面阵列均匀入射辐射通量,测量每个探测单元输出响应,并按照线性假设条件对非制冷红外焦平面阵列的响应曲线进行拟合,根据拟合曲线计算探测器各个单元在入射辐照度变化后非均匀校正参数,通过构造像元响应模型中单偏移校正参数存储模型、把单偏移参数存储模型聚合生成复合校正存储模型,根据像元和采集数据地址获取复合校正存储模型,从复合校正存储模型中提取单偏移校正参数存储模型,从单偏移校正参数存储模型中提取单个校正偏移参数,整个提取过程经过若干次的识别确认,有效地保证了非制冷红外焦平面阵列上任一探测单元校正结果的准确性。
附图说明
图1为数字图像处理模块的结构图。
具体实施方式
实施例1:
红外成像系统包括非制冷焦平面阵列型红外探测器、模拟信号处理模块、数字图像处理模块、电源控制模块和监视器,其中,数字图像处理模块的结构如图1所示,数字图像处理模块使用FPGA芯片作为核心处理器,FPGA芯片一方面通过视频帧存取控制口与用于存储经AD转换后数字图像数据的SD RAM进行数据交互、另一方面通过视频帧存取控制口与用于存储非均匀校正数据和硬件程序的SRAM进行数据交互;并且FPGA芯片一方面采用UART接口与运行有非均匀性校正算法的PC机进行数据交互、另一方面采用DA接口把标准视频图像信号传输至监视器;
数字图像处理模块还包括电源管理电路部分,FPGA芯片通过电源管理接口与电源控制模块进行数据交互;
在PC机上运行有非均匀性校正单个偏移参数的校正算法,该校正算法包括如下具体实施步骤:
步骤1,在光路中插入均匀辐射的标准校正参考源,依次采集标准校正参考源在低温TL和高温TH情况下非制冷红外焦平面阵列上任一个像元的响应值,具体实施步骤如下:
步骤1-1,设定在低温TL时标准校正参考源的辐射通量为ΦL、在高温TH时标准校正参考源的辐射通量为ΦH;
步骤1-2,设定非制冷红外焦平面阵列上第i行第j列像元为ppij,并且像元ppij的采集数据地址标记为IDppij∈{0,1}n,n表示固定长度;
步骤1-3,构造非制冷红外焦平面阵列上像元的响应模型,该响应模型基于如下假设条件:非制冷红外焦平面阵列上任一个像元在一定积分时间内对暗电流和背景噪声的响应为固定值,非制冷红外焦平面阵列上任一个像元的响应值与辐射源辐射波长无关,仅与光强度相关,且与辐照强度呈线性关系;
基于上述假设条件的像元响应模型为:
Rvij(Φ)=Gij×Φ+Ofij (1)
式(1)中,Φ为均匀辐射源辐射通量;
Rvij(Φ)为辐射通量为Φ时,非制冷红外焦平面阵列上第i行第j列像元的响应值;
Gij为非制冷红外焦平面阵列上第i行第j列像元的响应增益系数;
Ofij为非制冷红外焦平面阵列上第i行第j列像元的响应偏置系数;
步骤1-4,获取非制冷红外焦平面阵列上第i行第j列像元ppij在辐射通量ΦL下的响应值RvL-ppij、在辐射通量ΦH下的响应值RvH-ppij,根据步骤1-3的像元响应模型获得像元ppij的响应方程,该响应方程具体如下:
RvL-ppij=Gppij×ppij(ΦL)+Ofppij (2)
RvH-ppij=Gppij×ppij(ΦH)+Ofppij (3)
式(2)和(3)中,Gppij为非制冷红外焦平面阵列上第i行第j列像元ppij的校正响应增益系数,Ofppij为非制冷红外焦平面阵列上第i行第j列像元ppij的校正响应偏置系数;
步骤1-5,输出隶属于平面探测器上第i行第j列像元ppij的单个偏移参数的校正结果,其具体的实施步骤如下:
通过求解像元ppij的响应方程(2)和(3),输出像元ppij的单个偏移参数-校正响应增益系数Gppij,其具体计算过程如下:
Gppij=(RvH-ppij-RvL-ppij)/(ppij(ΦH)-ppij(ΦL)) (4)
通过求解像元ppij的响应方程(2)和(3),输出像元ppij的单个偏移参数-校正响应偏置系数Ofppij,其具体计算过程如下:
Ofppij={RvH-ppijppij(ΦL)-RvL-ppijppij(ΦH)}/(ppij(ΦH)-ppij(ΦL)) (5)
实施例2:
首先构造像元响应模型中单个偏移参数的存储模型、然后聚合生成复合校正偏移参数的存储模型,其具体实施步骤如下:
步骤2-1,定义公开参数PPs={GT,G1,G2,e,n1,p,P,Q,H0,H1,H2},其中的各个参数具体含义如下:
GT是n1阶循环群,G1和G2是阶为p的乘法群,G1和G2的生成元分别为P和Q,e表示双线性映射:G1×G2→GT,双线性映射e满足下述条件:
(1)对于任意a,b∈Zn1={0,1,2,…,n1-1},e(Pa,Qb)=e(P,Q)ab;
(2)存在u∈G1,v∈G2,使得e(u,v)≠0,其中0是单位元;
(3)对于所有的u∈G1,v∈G2,可以计算得到e(u,v);
H0表示散列函数:{0,1}*→Zp *={0,1,2,…,p-1},*表示任意长度;
H1表示散列函数:{0,1}*×G1→G2;
H2表示散列函数:{0,1}*→G1;
步骤2-2,首先设定像元ppij的校正偏移参数的存储标记xppij∈Zp *={0,1,2,…,p-1},然后计算像元ppij的校正偏移参数的第一提取系数Pppij=xppijP和第二提取系数Tppij=H2(P),并且计算像元ppij的校正偏移参数的提取标记Qppij=H1(IDppij||Pppij);
其中,||表示字符串连接;
步骤2-3,根据像元ppij的采集数据地址IDppij,生成像元ppij的校正响应增益系数Gppij的存储模型SMGppij,其具体生成步骤如下:
首先设定存储模型SMGppij的存储系数Grppij∈Zp *={0,1,2,…,p-1},然后执行如下的计算过程:
基于存储系数Grppij的第一提取参数GRppij=GrppijP;
基于存储系数Grppij的第二提取参数Ghppij=H0(IDppij||Gppij||Pppij||GRppij);
基于存储系数Grppij的第三提取参数GVppij=Qppij+GhppijGrppijTppij+xppijQ;
输出存储模型SMGppij为{ppij,IDppij,Gppij,GRppij,GVppij};
步骤2-4,根据像元ppij的采集数据地址IDppij,生成像元ppij的校正响应偏置系数Ofppij的存储模型SMOfppij,其具体生成步骤如下:
首先设定存储模型SMOfppij的存储系数Orppij∈Zp *={0,1,2,…,p-1},然后执行如下的计算过程:
基于存储系数Orppij的第一提取参数ORppij=OrppijP;
基于存储系数Orppij的第二提取参数Ohppij=H0(IDppij||Ofppij||Pppij||ORppij);
基于存储系数Orppij的第三提取参数OVppij=Qppij+OhppijOrppijTppij+xppijQ;
输出存储模型SMOfppij为{ppij,IDppij,Ofppij,ORppij,OVppij};
步骤2-5,生成像元ppij的复合校正偏移参数的存储模型SMppij,其具体生成步骤如下:
根据存储模型SMGppij的第一提取参数GRppij和存储模型SMOfppij的第一提取参数ORppij,计算像元ppij的复合校正偏移参数的第一提取参数Rppij=GhppijGRppij+OhppijORppij;
根据存储模型SMGppij的第三提取参数GVppij和存储模型SMOfppij的第三提取参数OVppij,计算像元ppij的复合校正偏移参数的第三提取参数Vppij=GVppij+OVppij;
输出存储模型SMppij为{ppij,IDppij,Rppij,Vppij,存储模型SMGppij,存储模型SMOfppij},该存储模型SMppij存储在SRAM中。
实施例3:
采用存储模型SMppij对非制冷红外焦平面阵列上第i行第j列像元ppij在任一辐射通量ΦK(ΦL≤ΦK≤ΦH)下的响应值RvK-ppij进行校正后输出,其具体的实施步骤如下:
获取非制冷红外焦平面阵列上第i行第j列像元ppij在任一辐射通量ΦK(ΦL≤ΦK≤ΦH)下的响应值RvK-ppij,同步地FPGA芯片通过时序控制获得像元ppij的采集数据地址IDppij;
实施例4:
根据像元ppij和采集数据地址IDppij,FPGA芯片从SRAM中读取得到对应的存储模型SMppij={ppij,IDppij,Rppij,Vppij,存储模型SMGppij,存储模型SMOfppij}、并从中提取存储模型SMGppij和存储模型SMOfppij,其具体的实施步骤如下:
根据已知的信息(公开参数PPs、像元ppij、采集数据地址IDppij、第一提取系数Pppij、第二提取系数Tppij、提取标记Qppij),FPGA芯片具体执行下述计算过程:
像元ppij的复合校正偏移参数的第一识别系数M1,其具体计算结果如下:
M1=e(Vppij,P);
像元ppij的复合校正偏移参数的第二识别系数M2,其具体计算结果如下:
M2=e(2Qppij,P)e(Tppij,Rppij)e(Q,2Pppij);
若第一识别系数M1与第二识别系数M2相等,则提取存储模型SMGppij和存储模型SMOfppij;
其中,第一识别系数M1与第二识别系数M2的方案可行性验证过程如下:
M1=e(Vppij,P)
=e(GVppij+OVppij,P)
=e(GhppijGrppijTppij+OhppijOrppijTppij+2Qppij+2xppijQ,P)
=e((GhppijGrppij+OhppijOrppij)Tppij,P)e(2Qppij,P)e(2xppijQ,P)
=e(Tppij,(GhppijGrppij+OhppijOrppij)P)e(2Qppij,P)e(2xppijQ,P)
=e(Tppij,Rppij)e(2Qppij,P)e(2xppijQ,P)
=e(Tppij,Rppij)e(2Qppij,P)e(Q,2Pppij)=M2。
实施例5:
FPGA芯片根据存储模型SMGppij,计算校正响应增益系数Gppij的第一识别系数GM1和第二识别系数GM2,具体执行下述计算过程:
第一识别系数GM1=e(GVppij,P);
第二识别系数GM2=e(Qppij,P)e(Tppij,GhppijGRppij)e(Q,Pppij);
若校正响应增益系数Gppij的第一识别系数GM1与第二识别系数GM2相等,则从存储模型SMGppij中提取校正响应增益系数Gppij;
其中,第一识别系数GM1与第二识别系数GM2的方案可行性验证过程如下:
GM1=e(GVppij,P)
=e(Qppij+GhppijGrppijTppij+xppijQ,P)
=e(Qppij,P)e(GhppijGrppijTppij,P)e(xppijQ,P)
=e(Qppij,P)e(Tppij,GhppijGrppijP)e(Q,xppijP)
=e(Qppij,P)e(Tppij,GhppijGRppij)e(Q,xppijP)=GM2;
FPGA芯片把校正响应增益系数Gppij传输至乘法器。
实施例6:
FPGA芯片根据存储模型SMOfppij,,计算校正响应偏置系数Ofppij的第一识别系数OM1和第二识别系数OM2,具体执行下述计算过程:
第一识别系数OM1=e(OVppij,P);
第二识别系数OM2=e(Qppij,P)e(Tppij,OhppijORppij)e(Q,Pppij);
若校正响应偏置系数Ofppij的第一识别系数OM1与第二识别系数OM2相等,则从存储模型存储模型SMOfppij中提取校正响应偏置系数Ofppij;
其中,第一识别系数OM1与第二识别系数OM2的方案可行性验证过程如下:
OM1=e(OVppij,P)
=e(Qppij+OhppijOrppijTppij+xppijQ,P)
=e(Qppij,P)e(OhppijOrppijTppij,P)e(xppijQ,P)
=e(Qppij,P)e(Tppij,OhppijOrppijP)e(Q,xppijP)
=e(Qppij,P)e(Tppij,OhppijORppij)e(Q,xppijP)=OM2;
FPGA芯片把校正响应偏置系数Ofppij传输至加法器;
实施例7:
根据校正响应增益系数Gppij和校正响应偏置系数Ofppij,输出非制冷红外焦平面阵列上像元ppij在任一辐射通量ΦK(ΦL≤ΦK≤ΦH)下响应值RvK-ppij的校正值ERv-ppij(ΦK)=Gppij×Rvppij(ΦK)+Ofppij,至此完成红外图像数据的非均匀性校正处理。
Claims (2)
1.一种非均匀性校正单个偏移参数的校正算法,其特征在于,所述校正算法包括如下步骤:
步骤一,在光路中插入均匀辐射的标准校正参考源,通过采集标准校正参考源在低温和高温情况下非制冷红外焦平面阵列上任一个像元的响应值,计算得到隶属于像元的校正系数,其具体执行过程为:
定义在低温TL时标准校正参考源的辐射通量为ΦL、在高温TH时标准校正参考源的辐射通量为ΦH;
定义非制冷红外焦平面阵列上第i行第j列像元为ppij,并且像元ppij的采集数据地址标记为IDppij∈{0,1}n,n表示固定长度;
构造非制冷红外焦平面阵列上像元的响应模型为:
Rvij(Φ)=Gij×Φ+Ofij (1)
其中,Φ为均匀辐射源辐射通量;Rvij(Φ)为辐射通量为Φ时,非制冷红外焦平面阵列上第i行第j列像元的响应值;Gij为非制冷红外焦平面阵列上第i行第j列像元的响应增益系数;Ofij为非制冷红外焦平面阵列上第i行第j列像元的响应偏置系数;
获取非制冷红外焦平面阵列上第i行第j列像元ppij在辐射通量ΦL下的响应值RvL-ppij、在辐射通量ΦH下的响应值RvH-ppij,根据像元响应模型获得像元ppij的响应方程,该响应方程具体如下:
RvL-ppij=Gppij×ppij(ΦL)+Ofppij (2)
RvH-ppij=Gppij×ppij(ΦH)+Ofppij (3)
其中,Gppij为非制冷红外焦平面阵列上第i行第j列像元ppij的校正响应增益系数,Ofppij为非制冷红外焦平面阵列上第i行第j列像元ppij的校正响应偏置系数;
通过求解像元ppij的响应方程(2)和(3),输出像元ppij的校正响应增益系数Gppij和校正响应偏置系数Ofppij:
Gppij=(RvH-ppij-RvL-ppij)/(ppij(ΦH)-ppij(ΦL))(4)
Ofppij={RvH-ppijppij(ΦL)-RvL-ppijppij(ΦH)}/(ppij(ΦH)-ppij(ΦL))(5)
步骤二,构造像元响应模型中单个偏移参数的存储模型、然后聚合生成复合校正偏移参数的存储模型,其具体执行过程为:
定义公开参数PPs={GT,G1,G2,e,n1,p,P,Q,H0,H1,H2},其中的各个参数具体含义如下:
GT是n1阶循环群,G1和G2是阶为p的乘法群,G1和G2的生成元分别为P和Q,e表示双线性映射:G1×G2→GT;
H0表示散列函数:{0,1}*→Zp *={0,1,2,…,p-1},*表示任意长度;
H1表示散列函数:{0,1}*×G1→G2;
H2表示散列函数:{0,1}*→G1;
定义像元ppij的校正偏移参数的存储标记xppij∈Zp *,读取像元ppij的采集数据地址IDppij,据此计算像元ppij的校正偏移参数的第一提取系数Pppij、第二提取系数Tppij和提取标记Qppij:
Pppij=xppijP;
Tppij=H2(P);
Qppij=H1(IDppij||Pppij);
其中,||表示字符串连接;
生成像元ppij的校正响应增益系数Gppij的存储模型SMGppij,其具体生成过程如下:
获取像元ppij的采集数据地址IDppij;
定义存储模型SMGppij的存储系数Grppij∈Zp *;
计算存储模型SMGppij的第一提取参数GRppij、第二提取参数Ghppij、第三提取参数GVppij:
GRppij=GrppijP;
Ghppij=H0(IDppij||Gppij||Pppij||GRppij);
GVppij=Qppij+GhppijGrppijTppij+xppijQ;
输出存储模型SMGppij为{ppij,IDppij,Gppij,GRppij,GVppij};
生成像元ppij的校正响应偏置系数Ofppij的存储模型SMOfppij,其具体生成过程如下:
获取像元ppij的采集数据地址IDppij;
定义存储模型SMOfppij的存储系数Orppij∈Zp *;
计算基于存储系数Orppij的第一提取参数ORppij、第二提取参数Ohppij和第三提取参数OVppij:
ORppij=OrppijP;
Ohppij=H0(IDppij||Ofppij||Pppij||ORppij);
OVppij=Qppij+OhppijOrppijTppij+xppijQ;
输出存储模型SMOfppij为{ppij,IDppij,Ofppij,ORppij,OVppij};
聚合生成像元ppij的复合校正偏移参数的存储模型SMppij,其具体生成过程如下:
根据存储模型SMGppij的第一提取参数GRppij和存储模型SMOfppij的第一提取参数ORppij,计算像元ppij的复合校正偏移参数的第一提取参数Rppij=GhppijGRppij+OhppijORppij;
根据存储模型SMGppij的第三提取参数GVppij和存储模型SMOfppij的第三提取参数OVppij,计算像元ppij的复合校正偏移参数的第三提取参数Vppij=GVppij+OVppij;
输出存储模型SMppij为{ppij,IDppij,Rppij,Vppij,存储模型SMGppij,存储模型SMOfppij},该存储模型SMppij存储在SRAM中;
获取非制冷红外焦平面阵列上第i行第j列像元ppij在任一辐射通量ΦK下的响应值RvK-ppij,同步地FPGA芯片通过时序控制获得像元ppij的采集数据地址IDppij;其中,ΦL≤ΦK≤ΦH;
根据像元ppij和采集数据地址IDppij,FPGA芯片从SRAM中读取存储模型SMppij={ppij,IDppij,Rppij,Vppij,存储模型SMGppij,存储模型SMOfppij};
之后结合公开参数PPs、像元ppij、采集数据地址IDppij、第一提取系数Pppij、第二提取系数Tppij及提取标记Qppij,FPGA芯片计算像元ppij的复合校正偏移参数的第一识别系数M1和第二识别系数M2:
M1=e(Vppij,P);
M2=e(2Qppij,P)e(Tppij,Rppij)e(Q,2Pppij);
若第一识别系数M1与第二识别系数M2相等,则提取存储模型SMGppij和存储模型SMOfppij;
FPGA芯片根据存储模型SMGppij,计算校正响应增益系数Gppij的第一识别系数GM1和第二识别系数GM2:
GM1=e(GVppij,P);
GM2=e(Qppij,P)e(Tppij,GhppijGRppij)e(Q,Pppij);
若校正响应增益系数Gppij的第一识别系数GM1与第二识别系数GM2相等,则FPGA芯片从存储模型SMGppij中提取校正响应增益系数Gppij,并把校正响应增益系数Gppij传输至乘法器;
FPGA芯片根据存储模型SMOfppij,计算校正响应偏置系数Ofppij的第一识别系数OM1和第二识别系数OM2:
OM1=e(OVppij,P);
OM2=e(Qppij,P)e(Tppij,OhppijORppij)e(Q,Pppij);
若校正响应偏置系数Ofppij的第一识别系数OM1与第二识别系数OM2相等,则FPGA芯片从存储模型SMOfppij中提取校正响应偏置系数Ofppij,并把校正响应偏置系数Ofppij传输至加法器;
根据校正响应增益系数Gppij和校正响应偏置系数Ofppij,输出非制冷红外焦平面阵列上像元ppij在辐射通量ΦK下响应值RvK-ppij的校正值ERv-ppij(ΦK):
ERv-ppij(ΦK)=Gppij×Rvppij(ΦK)+Ofppij;
至此,红外图像数据的非均匀性校正处理程序执行完毕。
2.一种PC机,该PC机上存储有计算机程序,其特征在于,所述计算机程序被执行时以实现权利要求1中所述的一种非均匀性校正单个偏移参数的校正算法的步骤。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211511725.8A CN115752759B (zh) | 2022-11-29 | 2022-11-29 | 一种非均匀性校正单个偏移参数的校正算法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211511725.8A CN115752759B (zh) | 2022-11-29 | 2022-11-29 | 一种非均匀性校正单个偏移参数的校正算法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115752759A CN115752759A (zh) | 2023-03-07 |
CN115752759B true CN115752759B (zh) | 2024-08-23 |
Family
ID=85340314
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202211511725.8A Active CN115752759B (zh) | 2022-11-29 | 2022-11-29 | 一种非均匀性校正单个偏移参数的校正算法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115752759B (zh) |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102609923A (zh) * | 2012-03-09 | 2012-07-25 | 深圳职业技术学院 | 红外图像处理方法及装置 |
CN104251742A (zh) * | 2014-09-30 | 2014-12-31 | 成都市晶林科技有限公司 | 基于两点线性和目标、环境的二元非线性红外探测器非均匀性校正方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6433333B1 (en) * | 2000-03-03 | 2002-08-13 | Drs Sensors & Targeting Systems, Inc. | Infrared sensor temperature compensated response and offset correction |
-
2022
- 2022-11-29 CN CN202211511725.8A patent/CN115752759B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102609923A (zh) * | 2012-03-09 | 2012-07-25 | 深圳职业技术学院 | 红外图像处理方法及装置 |
CN104251742A (zh) * | 2014-09-30 | 2014-12-31 | 成都市晶林科技有限公司 | 基于两点线性和目标、环境的二元非线性红外探测器非均匀性校正方法 |
Also Published As
Publication number | Publication date |
---|---|
CN115752759A (zh) | 2023-03-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN116055908B (zh) | 一种非制冷红外焦平面阵列非均匀性校正方法 | |
CN101361360B (zh) | 用于在使用光学黑色像素及束缚像素的成像器中设定黑色电平的方法及设备 | |
CN105865723A (zh) | 气体泄漏检测非均匀校正方法及气体泄漏检测装置 | |
CN101635782B (zh) | 一种基于动态时间延迟积分的图像获取方法及器件 | |
CN113709382B (zh) | 一种图像传感器、图像处理方法及存储介质 | |
CN114529830B (zh) | 基于混合卷积网络遥感图像时空融合方法 | |
CN105466563A (zh) | 一种两点校正红外热像仪的非均匀性的模块及方法 | |
CN104219511A (zh) | 色彩校正系统及方法 | |
CN204286600U (zh) | 一种两点校正红外热像仪的非均匀性的模块 | |
CN115752759B (zh) | 一种非均匀性校正单个偏移参数的校正算法 | |
CA2451463C (en) | Method and apparatus for readout of compound microbolometer arrays | |
Iwert et al. | The challenge of highly curved monolithic imaging detectors | |
CN106374224A (zh) | 电磁波成像系统及天线阵列信号校正方法 | |
CN107566763B (zh) | 大尺寸图像传感器及其图像校正方法 | |
CN113259558B (zh) | 一种光场相机的无损全聚焦方法及装置 | |
CN110392186B (zh) | 减少雾霾影响的成像装置及成像方法 | |
CN117451190A (zh) | 一种深度学习的散焦散射波前传感方法 | |
CN109949363A (zh) | 一种适用于太赫兹仿生复眼成像系统的目标识别定位系统 | |
WO2020071104A1 (ja) | 固体撮像装置と駆動方法および電子機器 | |
CN117522685A (zh) | 基于地面图像辅助的遥感图像超分辨重建方法及装置 | |
CN110487412A (zh) | 红外高光谱图像非均匀性校正方法、装置和计算机设备 | |
CN101852650A (zh) | 提高红外热像仪测温均匀性的装置及方法 | |
CN109118460A (zh) | 一种分光偏振光谱信息同步处理方法及系统 | |
CN105046674A (zh) | 一种多元并扫红外ccd图像的非均匀化校正方法 | |
CN108491852A (zh) | 一种面阵红外焦平面的盲元检测方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |