CN112070663A - 一种多探测器拼接型相机辐射定标联合数据处理方法 - Google Patents
一种多探测器拼接型相机辐射定标联合数据处理方法 Download PDFInfo
- Publication number
- CN112070663A CN112070663A CN202010724073.0A CN202010724073A CN112070663A CN 112070663 A CN112070663 A CN 112070663A CN 202010724073 A CN202010724073 A CN 202010724073A CN 112070663 A CN112070663 A CN 112070663A
- Authority
- CN
- China
- Prior art keywords
- detector
- image
- radiance
- pixel
- grade
- 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.)
- Granted
Links
- 238000003672 processing method Methods 0.000 title claims abstract description 22
- 238000012937 correction Methods 0.000 claims abstract description 24
- 238000000034 method Methods 0.000 claims abstract description 23
- 238000012545 processing Methods 0.000 claims abstract description 10
- 239000000126 substance Substances 0.000 claims description 10
- 230000005855 radiation Effects 0.000 description 7
- 238000010586 diagram Methods 0.000 description 5
- 238000004364 calculation method Methods 0.000 description 3
- 238000003384 imaging method Methods 0.000 description 3
- 230000007547 defect Effects 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000000903 blocking effect Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T3/00—Geometric image transformations in the plane of the image
- G06T3/40—Scaling of whole images or parts thereof, e.g. expanding or contracting
- G06T3/4038—Image mosaicing, e.g. composing plane images from plane sub-images
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C25/00—Manufacturing, calibrating, cleaning, or repairing instruments or devices referred to in the other groups of this subclass
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/80—Geometric correction
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10032—Satellite or aerial image; Remote sensing
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Manufacturing & Machinery (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Photometry And Measurement Of Optical Pulse Characteristics (AREA)
- Image Processing (AREA)
Abstract
本发明公开了一种多探测器拼接型相机辐射定标联合数据处理方法,所述方法包括如下步骤:计算每个探测器各档辐亮度下图像平均灰度值;采用最小二乘法拟合出第一直线,根据第一直线得到每个探测器各像元的增益拟合系数和偏置拟合系数;得到Q个探测器各档辐亮度下的图像平均灰度值;采用最小二乘法拟合出第二直线,根据第二直线得到每个探测器的加权调整系数;根据每个探测器各像元的拟合系数和每个探测器的加权调整系数得到Q个探测器各像元的最终相对辐射定标校正公式,根据Q个探测器各像元的最终相对辐射定标校正公式得到辐射定标处理图像。本发明通过辐射定标联合数据处理方法,能够获取高精度辐射定标处理结果。
Description
技术领域
本发明属于多探测器拼接型相机技术领域,尤其涉及一种多探测器拼接型相机辐射定标联合数据处理方法。
背景技术
光学遥感相机主要任务是获取地面目标的几何信息和辐射信息,相机输出图像的灰度包含地物目标的辐射信息,而图像灰度幅值、信噪比和对比度将直接影响到对遥感图像的分析、判读和解译。相机辐射性能的全面测试对于相机的研制和应用至关重要。
为满足宽视场、高分辨率成像的需求,多镜头多探测器拼接型相机是研制的热点和目前的趋势。不同的拼接方法对数据处理的几何精度和辐射精度具有不同的技术难点,多镜头多探测器相机一般是外视场拼接成像方式,通过外视场拼接最终输出拼接后大幅面图像,各探测器和镜头的不一致性,带来的缺陷是采用传统辐射定标校正方法,各探测器内像元的一致性会处理较好,但各探测器间会存在偏差,导致最终拼接图像具有明显分块的效应。
发明内容
本发明解决的技术问题是:克服现有技术的不足,提供了一种探测器拼接型相机辐射定标联合数据处理方法,针对通过视场拼接最终输出拼接后大幅面响应不一致性问题,通过辐射定标联合数据处理方法,能够获取高精度辐射定标处理结果。
本发明目的通过以下技术方案予以实现:一种多探测器拼接型相机辐射定标联合数据处理方法,所述方法包括如下步骤:步骤一:多探测器拼接型相机包括Q个探测器,每个探测器像元数为m×n,采集Q个探测器在积分球输出L档辐亮度下的图像数据;其中,每个探测器获得N幅图像数据,每幅图像数据大小为m×n;m为每个探测器像元的总行数,n为每个探测器像元的总列数;步骤二:计算每个探测器各档辐亮度下图像平均灰度值;步骤三:根据每个探测器各档辐亮度下图像平均灰度值和每个探测器中每个像元在各档辐亮度下图像灰度DN(i,j,k,l)采用最小二乘法拟合出第一直线,根据第一直线得到每个探测器各像元的增益拟合系数和偏置拟合系数A(i,j,k)、B(i,j,k);步骤四:根据步骤二中的每个探测器各档辐亮度下图像平均灰度值得到Q个探测器各档辐亮度下的图像平均灰度值;步骤五:根据Q个探测器各档辐亮度下的图像平均灰度值和每个探测器各档辐亮度下图像平均灰度值采用最小二乘法拟合出第二直线,根据第二直线得到每个探测器的加权调整系数Ck、Dk;步骤六:根据每个探测器各像元的拟合系数A(i,j,k)、B(i,j,k)和每个探测器的加权调整系数Ck、Dk得到Q个探测器各像元的最终相对辐射定标校正公式,根据Q个探测器各像元的最终相对辐射定标校正公式得到辐射定标处理图像。
上述多探测器拼接型相机辐射定标联合数据处理方法中,在步骤二中,每个探测器各档辐亮度下图像平均灰度值通过以下公式得到:
其中,为第k个探测器在第l档辐亮度下图像的平均灰度值,i=1、2、······、m,i表示每个探测器像元的行数,j=1、2、······、n,j表示每个探测器像元的列数,k=1、2、······、Q,k表示探测器的序号,l=1、2、······、L,l表示第l档辐亮度。
上述多探测器拼接型相机辐射定标联合数据处理方法中,在步骤三中,第一直线为:
其中,表示第k个探测器的第i行j列像元在第l档辐亮度下的图像灰度观测值,表示第k个探测器的第i行j列像元在第l档辐亮度下的图像灰度估计值,A(i,j,k)表示第k个探测器的第i行j列像元的增益拟合系数,B(i,j,k)表示第k个探测器的第i行j列像元的偏置拟合系数。
上述多探测器拼接型相机辐射定标联合数据处理方法中,在步骤四中,Q个探测器各档辐亮度下图像平均灰度值通过以下公式得到:
上述多探测器拼接型相机辐射定标联合数据处理方法中,在步骤五中,第二直线为:
上述多探测器拼接型相机辐射定标联合数据处理方法中,在步骤六中,Q个探测器各像元的最终相对辐射定标校正公式为:
其中,a(i,j,k)=CkA(i,j,k),a(i,j,k)为第k个探测器第i行j列像元的增益校正参数,b(i,j,k)=CkB(i,j,k)+Dk,b(i,j,k)为第k个探测器第i行j列像元的偏置校正参数;为第k个探测器第i行j列像元的图像灰度观测值,为第k个探测器第i行j列像元的图像灰度校正后估计值。
本发明与现有技术相比具有如下有益效果:
本发明对多探测器拼接型相机辐射定标联合数据处理,最终输出拼接后高精度辐射定标处理图像,获取高精度辐射定标处理结果。图像的质量会相应提升,对高精度光学遥感相机具有重要意义。
附图说明
通过阅读下文优选实施方式的详细描述,各种其他的优点和益处对于本领域普通技术人员将变得清楚明了。附图仅用于示出优选实施方式的目的,而并不认为是对本发明的限制。而且在整个附图中,用相同的参考符号表示相同的部件。在附图中:
图1是本发明实施例提供的探测器拼接型相机辐射定标联合数据处理方法的流程图;
图2是本发明实施例提供的获取的待处理图像数据示意图;
图3是本发明实施例提供的4个探测器拼接原理示意图;
图4是本发明实施例提供的4个探测器原始图像示意图;
图5是本发明实施例提供的4个探测器原始图像单探测器校正处理示意图;
图6是本发明实施例提供的4个探测器联合校正处理处理示意图。
具体实施方式
下面将参照附图更详细地描述本公开的示例性实施例。虽然附图中显示了本公开的示例性实施例,然而应当理解,可以以各种形式实现本公开而不应被这里阐述的实施例所限制。相反,提供这些实施例是为了能够更透彻地理解本公开,并且能够将本公开的范围完整的传达给本领域的技术人员。需要说明的是,在不冲突的情况下,本发明中的实施例及实施例中的特征可以相互组合。下面将参考附图并结合实施例来详细说明本发明。
图1是本发明实施例提供的探测器拼接型相机辐射定标联合数据处理方法的流程图。如图1所示,该探测器拼接型相机辐射定标联合数据处理方法包括如下步骤:
第一步,完成Q个探测器在积分球输出L档辐亮度下的图像数据采集。Q个探测器,每个探测器获得N幅图像,每幅图像大小为m×n。获取的待处理图像数据示意如图2。多探测器拼接型相机含Q个探测器,每个探测器像元数为m×n。
第二步,计算每个探测器各等效辐亮度下图像平均灰度值。
其中,i=1,2,······,m表示探测器像元的行数,j=1,2,······,n表示探测器像元的列数,k=1,2,······,Q表示第k个探测器,l=1,2,······,L表示第l档辐亮度。
重复上述方法,计算出每个探测器各等效辐亮度下图像平均灰度值。
第三步,计算每个探测器各像元的拟合系数A(i,j,k)、B(i,j,k)。其中i=1,2,······,m表示探测器像元的行数,j=1,2,······,n表示探测器像元的列数,k=1,2,······,Q表示第k个探测器;A(i,j,k)表示第k个探测器的第i行j列像元的增益拟合系数,B(i,j,k)表示第k个探测器的第i行j列像元的偏置拟合系数。
具体计算方法如下。
(2)重复步骤(1),计算出第k个探测器的所有像元的拟合系数。
(3)重复步骤(1)、(2),计算出Q探测器的所有像元的拟合系数。
第四步,计算Q个探测器各等效辐亮度下整体的图像平均灰度值。
第五步,计算各探测器的加权调整系数Ck、Dk。
具体计算方法如下。
至此,计算出了第k个探测器的加权调整系数Ck、Dk。
重复上述方法,计算出Q探测器的所有加权调整系数。
第六步,获取Q个探测器各像素点的最终相对辐射定标校正公式
其中,a(i,j,k)=CkA(i,j,k)表示第k个探测器第i行j列像元的增益校正参数,b(i,j,k)=CkB(i,j,k)+Dk表示第k个探测器第i行j列像元的偏置校正参数;表示第k个探测器第i行j列像元的图像灰度观测值,表示第k个探测器第i行j列像元的图像灰度校正后估计值。
多探测器拼接型相机含4个探测器,每个探测器像元数为3000×4000,其拼接成像原理如图3所示,技术方案的步骤如下:
第一步,完成4个探测器8档辐亮度图像数据采集,其中一档亮度原始图像如图4所示。各档图像数据采集状态对应的等效辐亮度为Ln(n=1,2,…,8),在等效辐亮度Ln下像元输出灰度值DN可表示为DN(i,j,k,l),其中i=1,2,…,3000,j=1,2,…,4000表示像元位置,p=1,2,…,8表示第p级辐亮度,k=1,2,…,4表示第k个探测器。
第二步,使用最小二乘法拟合出每个探测器各像元在8档辐亮度下的线性响应曲线。每个探测器各像元的拟合系数A(i,j,k)、B(i,j,k),每个像元相对辐射定标的基本公式为。校正后如图5所示。
第三步,计算各探测器等效辐亮度l下整幅图像的平均灰度值。第k个探测器的计算公式为
第四步,使用最小二乘法拟合出每个探测器在8档辐亮度下的平均灰度曲线,各探测器的拟合系数Ck、Dk,每个探测器平均灰度调整公式
第五步,利用最终校正公式,对每个探测器的每个像元进行校正处理,校正后图像如图6所示。最终校正公式为
本实施例以单探测器不同辐亮度下的平均灰度为观测值,采用最小二乘法拟合求取单探测器每个像元的相对定标系数,然后将所有探测器不同辐亮度下平均灰度为观测值,联合多探测器数据再次采用最小二乘法拟合求取各探测器之间相对定标调整加权系数,最终获取联合校正处理关系。
本发明虽然已以较佳实施例公开如上,但其并不是用来限定本发明,任何本领域技术人员在不脱离本发明的精神和范围内,都可以利用上述揭示的方法和技术内容对本发明技术方案做出可能的变动和修改,因此,凡是未脱离本发明技术方案的内容,依据本发明的技术实质对以上实施例所作的任何简单修改、等同变化及修饰,均属于本发明技术方案的保护范围。
Claims (6)
1.一种多探测器拼接型相机辐射定标联合数据处理方法,其特征在于,所述方法包括如下步骤:
步骤一:多探测器拼接型相机包括Q个探测器,每个探测器像元数为m×n,采集Q个探测器在积分球输出L档辐亮度下的图像数据;其中,每个探测器获得N幅图像数据,每幅图像数据大小为m×n;m为每个探测器像元的总行数,n为每个探测器像元的总列数;
步骤二:计算每个探测器各档辐亮度下图像平均灰度值;
步骤三:根据每个探测器各档辐亮度下图像平均灰度值和每个探测器中每个像元在各档辐亮度下图像灰度采用最小二乘法拟合出第一直线,根据第一直线得到每个探测器各像元的增益拟合系数和偏置拟合系数;
步骤四:根据步骤二中的每个探测器各档辐亮度下图像平均灰度值得到Q个探测器各档辐亮度下的图像平均灰度值;
步骤五:根据Q个探测器各档辐亮度下的图像平均灰度值和每个探测器各档辐亮度下图像平均灰度值采用最小二乘法拟合出第二直线,根据第二直线得到每个探测器的加权调整系数;
步骤六:根据每个探测器各像元的拟合系数和每个探测器的加权调整系数得到Q个探测器各像元的最终相对辐射定标校正公式,根据Q个探测器各像元的最终相对辐射定标校正公式得到辐射定标处理图像。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010724073.0A CN112070663B (zh) | 2020-07-24 | 2020-07-24 | 一种多探测器拼接型相机辐射定标联合数据处理方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010724073.0A CN112070663B (zh) | 2020-07-24 | 2020-07-24 | 一种多探测器拼接型相机辐射定标联合数据处理方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112070663A true CN112070663A (zh) | 2020-12-11 |
CN112070663B CN112070663B (zh) | 2024-03-29 |
Family
ID=73656967
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010724073.0A Active CN112070663B (zh) | 2020-07-24 | 2020-07-24 | 一种多探测器拼接型相机辐射定标联合数据处理方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112070663B (zh) |
Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20080088719A1 (en) * | 2005-04-29 | 2008-04-17 | Eliezer Jacob | Digital camera with non-uniform image resolution |
CN101226639A (zh) * | 2008-01-29 | 2008-07-23 | 航天东方红卫星有限公司 | 一种星载tdiccd相机的相对辐射校正方法 |
CN103247029A (zh) * | 2013-03-26 | 2013-08-14 | 中国科学院上海技术物理研究所 | 一种用于拼接式探测器生成的高光谱图像几何配准方法 |
CN104077740A (zh) * | 2013-03-29 | 2014-10-01 | 中国科学院国家天文台 | 月球遥感图像的灰度均衡处理方法 |
CN104851106A (zh) * | 2015-06-03 | 2015-08-19 | 西北大学 | 一种ccd原始图像的处理方法 |
CN104977024A (zh) * | 2015-05-28 | 2015-10-14 | 北京空间机电研究所 | 一种日盲紫外遥感相机绝对辐射定标系数在轨修正方法 |
CN106643669A (zh) * | 2016-11-22 | 2017-05-10 | 北京空间机电研究所 | 一种多镜头多探测器航空相机单中心投影转换方法 |
CN107063296A (zh) * | 2016-11-17 | 2017-08-18 | 许昌学院 | 一种卫星遥感传感器在轨辐射定标方法 |
US20180035033A1 (en) * | 2015-03-13 | 2018-02-01 | Aqueti Incorporated | Multi-array camera imaging system and method therefor |
CN109737987A (zh) * | 2018-12-29 | 2019-05-10 | 中国科学院长春光学精密机械与物理研究所 | 一种多光合一大口径空间相机在轨星上红外辐射定标系统 |
CN110120077A (zh) * | 2019-05-06 | 2019-08-13 | 航天东方红卫星有限公司 | 一种基于卫星姿态调整的面阵相机在轨相对辐射定标方法 |
-
2020
- 2020-07-24 CN CN202010724073.0A patent/CN112070663B/zh active Active
Patent Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20080088719A1 (en) * | 2005-04-29 | 2008-04-17 | Eliezer Jacob | Digital camera with non-uniform image resolution |
CN101226639A (zh) * | 2008-01-29 | 2008-07-23 | 航天东方红卫星有限公司 | 一种星载tdiccd相机的相对辐射校正方法 |
CN103247029A (zh) * | 2013-03-26 | 2013-08-14 | 中国科学院上海技术物理研究所 | 一种用于拼接式探测器生成的高光谱图像几何配准方法 |
CN104077740A (zh) * | 2013-03-29 | 2014-10-01 | 中国科学院国家天文台 | 月球遥感图像的灰度均衡处理方法 |
US20180035033A1 (en) * | 2015-03-13 | 2018-02-01 | Aqueti Incorporated | Multi-array camera imaging system and method therefor |
CN104977024A (zh) * | 2015-05-28 | 2015-10-14 | 北京空间机电研究所 | 一种日盲紫外遥感相机绝对辐射定标系数在轨修正方法 |
CN104851106A (zh) * | 2015-06-03 | 2015-08-19 | 西北大学 | 一种ccd原始图像的处理方法 |
CN107063296A (zh) * | 2016-11-17 | 2017-08-18 | 许昌学院 | 一种卫星遥感传感器在轨辐射定标方法 |
CN106643669A (zh) * | 2016-11-22 | 2017-05-10 | 北京空间机电研究所 | 一种多镜头多探测器航空相机单中心投影转换方法 |
CN109737987A (zh) * | 2018-12-29 | 2019-05-10 | 中国科学院长春光学精密机械与物理研究所 | 一种多光合一大口径空间相机在轨星上红外辐射定标系统 |
CN110120077A (zh) * | 2019-05-06 | 2019-08-13 | 航天东方红卫星有限公司 | 一种基于卫星姿态调整的面阵相机在轨相对辐射定标方法 |
Also Published As
Publication number | Publication date |
---|---|
CN112070663B (zh) | 2024-03-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102300057B (zh) | 线阵ccd像元响应不一致性校正方法 | |
CN102042878A (zh) | 一种去除温漂的红外非均匀性校正方法 | |
CN109459135A (zh) | 一种ccd成像光谱仪图像校正方法 | |
CN112529807B (zh) | 卫星影像的相对辐射校正方法和装置 | |
CN110009688A (zh) | 一种红外遥感影像相对辐射定标方法、系统及遥感平台 | |
CN107093196B (zh) | 视频卫星面阵相机在轨相对辐射定标方法 | |
CN104406699A (zh) | 基于自适应红外图像校正算法红外热像仪 | |
CN102752504B (zh) | 一种宽视场线阵ccd相机的相对辐射校正方法 | |
CN105869129B (zh) | 针对非均匀校正后的热红外图像剩余非均匀噪声去除方法 | |
CN101635782A (zh) | 一种基于动态时间延迟积分的图像获取方法及器件 | |
US7796153B1 (en) | Equalization system and method for an imaging sensor | |
CN114125431A (zh) | 用于静止轨道光学大面阵相机的非均匀性定标校正方法 | |
CN101420520A (zh) | 一种多摄像机阵列的联合标定方法 | |
CN116612080B (zh) | 基于高光谱图像波谱曲线的变化检测方法 | |
CN112070663A (zh) | 一种多探测器拼接型相机辐射定标联合数据处理方法 | |
US20050194539A1 (en) | Method for improving measurement accuracy of infrared imaging radiometers | |
CN113608246B (zh) | 一种光学卫星tdiccd积分级数和增益设置方法及成像装置 | |
CN114859099A (zh) | 一种无暗像元的可见光线阵探测器的暗电流实时估计方法 | |
CN113936065A (zh) | 基于恒星源的遥感相机辐射定标方法、系统、装置及介质 | |
CN110910442B (zh) | 一种基于kernel-free图像复原的高速移动对象机器视觉尺寸检测方法 | |
CN114494039A (zh) | 一种水下高光谱推扫图像几何校正的方法 | |
CN113970374A (zh) | 一种分焦平面偏振探测系统定标方法 | |
CN112529810A (zh) | 一种面阵凝视相机的探测信噪比提升方法 | |
CN110827362A (zh) | 一种基于多项式相机响应函数与渐晕效应补偿的光度标定方法 | |
CN113570523B (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 |