CN116612199A - 一种深空探测彩色相机在轨联合颜色定标方法 - Google Patents

一种深空探测彩色相机在轨联合颜色定标方法 Download PDF

Info

Publication number
CN116612199A
CN116612199A CN202310892257.1A CN202310892257A CN116612199A CN 116612199 A CN116612199 A CN 116612199A CN 202310892257 A CN202310892257 A CN 202310892257A CN 116612199 A CN116612199 A CN 116612199A
Authority
CN
China
Prior art keywords
color
orbit
camera
deep space
chromatic aberration
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
Application number
CN202310892257.1A
Other languages
English (en)
Other versions
CN116612199B (zh
Inventor
赵汝进
刘赛男
曾骏哲
乔丽
马跃博
余国彬
刘恩海
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Institute of Optics and Electronics of CAS
Original Assignee
Institute of Optics and Electronics of CAS
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Institute of Optics and Electronics of CAS filed Critical Institute of Optics and Electronics of CAS
Priority to CN202310892257.1A priority Critical patent/CN116612199B/zh
Publication of CN116612199A publication Critical patent/CN116612199A/zh
Application granted granted Critical
Publication of CN116612199B publication Critical patent/CN116612199B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/80Analysis of captured images to determine intrinsic or extrinsic camera parameters, i.e. camera calibration
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00Image enhancement or restoration
    • G06T5/77Retouching; Inpainting; Scratch removal
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/90Determination of colour characteristics
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/30Subject of image; Context of image processing
    • G06T2207/30244Camera pose
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02ATECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/30Assessment of water resources

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Spectrometry And Color Measurement (AREA)

Abstract

本发明公开了一种深空探测彩色相机在轨联合颜色定标方法,涉及颜色定标技术领域,该方法利用在轨光谱数据进行在轨颜色定标,结合标准人眼视觉色彩响应函数,反演目标颜色真值。通过最小二乘数据拟合以及色差迭代优化方法求解出最优定标系数矩阵,使得相机的观测效果与标准人眼观察者相一致。由于该定标方法为结合实时光谱数据进行解算,使其对光照环境具有自适应性。本发明解决了在轨无色标板、光照环境复杂等情况所导致的严重偏色的问题,提高了颜色定标的准确性和在轨适应性。

Description

一种深空探测彩色相机在轨联合颜色定标方法
技术领域
本发明涉及颜色定标技术领域,尤其涉及一种深空探测彩色相机在轨联合颜色定标方法。
背景技术
深空探测成像技术是人类探索脱离地球引力场的太阳系以外未知空间的重要技术途径,尤其是深空探测彩色相机,具有获取标准人眼观察者图像,根据其颜色、纹理信息勘察地形地质构造、重建天体三维形貌等一系列重大的科学和工程意义。但由于深空环境的特殊性、相机自身颜色响应特性与人眼标准观察者不一致等原因,会导致图像出现色差。颜色定标是解决色差问题最科学有效的技术途径。然而,目前该领域正面临着地面定标系数在轨适应性差且无法使用在轨色标板的技术难点。
随着火箭运载能力和航天技术的不断发展,人造卫星已经从近地卫星逐渐迈向深空探测卫星,但目前国内的主要在轨定标方法大部分仍然只适用于近地卫星,比如高光谱微纳卫星中的高光谱相机,利用地面人工靶标进行在轨颜色定标;“高分二号”卫星通过敦煌定标场人工靶标和周边均匀地物进行在轨辐射定标;“资源一号”卫星采用地面多灰阶人工靶标在轨绝对辐射定标;“天绘一号”卫星以海洋为暗场、阿尔及利亚的沙漠为亮场,对搭载的高分辨相机、多光谱相机、三线阵相机进行了在轨相对辐射定标。以上的在轨定标方法均是基于地球戈壁、干盐、湖、沙漠、海洋、湖泊、云层、以及人工靶标等参考目标的场地进行定标,这些场地的辐射、颜色、几何信息均可通过地面实地测量获得先验信息,而对于深空探测而言,目标的相关信息是未知且没有先验的,因此目前大部分的近地卫星在轨定标方法并不能满足深空探测的需求。我国深空探测主要包括“嫦娥系列”卫星,但“嫦娥系列”卫星上的光学相机载荷仅进行了地面定标试验,没有在轨定标试验方法,例如文献“赵汝进,刘恩海,王进,余国彬.嫦娥三号卫星相机地形图像的彩色校正方法[J]. 宇航学报,2016,37(03):341-347,2011.”中设计了一种颜色定标与白平衡相结合的两步颜色校正方法,经地面验证达到了满意的效果,但其在轨适应性并不理想。发射深空探测卫星最多的国家是美国,其也通过大量在轨定标提高在轨测量精准度。比如火星表面“毅力号”、“好奇号”火星车均搭载有在轨色彩定标板,其中“毅力号”搭载了三个色标板;美国探月工程“阿波罗17号”拍摄月球表面图像是也是通过在轨色标板进行颜色定标。然而,对于非停留在星球表面的探测器而言,比如环绕卫星上的光学载荷,由于难以搭载在轨色标板而没有进行在轨颜色定标。比如火星勘测轨道器上的火星彩色成像仪,无在轨颜色定标,通过7个窄带滤光片(258~718nm)进行彩色复原;火星轨道“维京号”彩色相机,无在轨颜色定标,通过6个宽带滤光片(350~700nm)进行彩色复原。
综上所述,国内外的深空探测卫星中,环绕器上的星载相机几乎没有在轨颜色定标方法。大部分现有在轨定标方法均是近地卫星以地表已知先验信息进行在轨定标,或是着陆器搭载在轨色标板进行定标。深空探测的在轨颜色定标技术成为当前的航天技术难点之一。
因此,需要提出一种方法以提高深空彩色相机颜色定标的准确性和在轨适应性。
发明内容
为解决在轨无色标板以及光照环境变化复杂等情况导致的偏色问题,本发明根据在轨实际情况提出了一种深空探测彩色相机在轨联合颜色定标方法,该方法结合在轨光谱数据进行在轨定标,使得相机的观测效果与标准人眼观察者相一致。
为达到上述目的,本发明采用的技术方案如下:
一种深空探测彩色相机在轨联合颜色定标方法,包括如下步骤:
S1、利用在轨的深空探测彩色相机拍摄1张星体表面图像,并获取其相机系统sRGB颜色空间下的色度值N组,进而归一化为待校正色品值;
S2、当深空探测彩色相机在轨对某一目标区域成像时,利用在轨探测器搭载的光谱仪获得相应的目标光谱数据;
S3、根据S2所获得的目标光谱数据的光谱曲线,结合标准人眼视觉色彩响应函数,解算XYZ颜色空间下的标准人眼响应色度值;
S4、根据颜色空间转换公式,将S3所求得的人眼响应色度值转换为sRGB颜色空间下的色度值,进而归一化为标准色品值;
S5、将S1所获得的待校正色品值的数据与S4所求得的标准色品值的数据做数据拟合,利用最小二乘法初步求取颜色定标系数矩阵;
S6、利用CIE Lab色差公式求解色差,得到色差权重矩阵,并通过色差权重矩阵重新优化颜色定标系数矩阵,直至得到满意的色差结果,得到最优的颜色定标系数矩阵;
S7、将S6得到的颜色定标系数矩阵的定标系数作用于深空探测彩色相机在轨拍摄图像,反演出图像真实色彩。
有益效果:
本发明根据在轨实际情况提出了一种在轨联合颜色定标方法,该方法利用在轨光谱数据进行在轨定标,结合标准人眼响应函数,反演目标真实色彩,通过最小二乘数据拟合以及色差迭代优化方法求解出最优定标系数矩阵,使得相机的观测效果与人眼标准观察者相一致,解决了在轨无色标板、光照环境复杂等情况导致的色差大的问题,提高了颜色定标的准确性和在轨应用适应性。
附图说明
图1为本发明的一种深空探测彩色相机在轨联合颜色定标方法流程图。
具体实施方式
以下结合具体实施例对本发明作进一步详细描述。
如图1所示,本发明的一种深空探测彩色相机在轨联合颜色定标方法,具体包括如下步骤:
S1、利用在轨的深空探测彩色相机拍摄1张星体表面图像,并获取其相机系统sRGB颜色空间下的色度值N组,对应N个像素,进而归一化为待校正色品值,以R通道为例,色度值归一化为色品值的计算方法均为:
其中,r为R通道色品值,R、G、B分别为三个通道的色度值;
S2、当深空探测彩色相机在轨对某一目标区域成像时,在轨探测器其他载荷——光谱仪也会对该区域进行探测,利用光谱仪获得N个像素相对应的目标光谱数据。
为波长,星体表面的反射光线的辐射光谱为/>,k为目标表面发出的总光功率与人眼入瞳接收到光功率的比例,则“标准人眼”接收到的光线的辐射光谱/>为:
S3、根据S2所获得的目标光谱数据的N组光谱曲线,结合标准人眼视觉色彩响应函数,解算XYZ颜色空间下的标准人眼响应色度值。
CIE光谱三刺激值是CIE以百位正常视觉者,用CIE规定的红、绿、蓝三原色光,对等能光谱色从380~780nm所进行的专门性颜色混合匹配实验得到的,反映了“标准人眼”视觉的色彩响应特征。
设CIE三刺激值曲线函数为、/>、/>,则目标表面反射光线在“标准人眼”中的X、Y、Z响应/>、/>、/>分别为:
,
,
,
S4、根据XYZ空间与sRGB 空间之间的标准转换关系式,将S3所求得的人眼响应色度值转换为sRGB颜色空间下的色度值,进而归一化为标准色品值。
“标准人眼”在sRGB颜色空间下的响应函数、/>为:
分别对、/>、/>求积分,得到人眼对目标的三色响应值/>,归一化可求得标准色品值/>
S5、图像真实色彩应无限接近“标准人眼”响应色彩,因此将S1所获得的待校正色品值的数据与S4所求得的标准色品值的数据做数据拟合,利用最小二乘法初步求取颜色定标系数矩阵。
在最小二乘法误差约束下,可解算在色彩空间下颜色定标系数矩阵
,
其中,,为N组标准色品值;/>对应为第N组三通道标准色品值;
,为N组待校正色品值;/>对应为第N组三通道待校正色品值;
,为颜色定标系数矩阵,/>对应为R、G、B三通道定标系数。
将得到的颜色定标系数矩阵作用于待校正色品值。
S6、利用CIE Lab色差公式求解色差,得到色差权重矩阵,并通过色差权重矩阵重新优化颜色定标系数矩阵,直至得到满意的色差结果,得到最优的颜色定标系数矩阵,包括:
先将RGB值转化为标准的XYZ色度系统下的三刺激值x、y、z,再将x、y、z数据转换到Lab颜色空间得到L、a、b值,其中L为明度因数,a和b为两个色度因数,根据以下公式计算色差
,
其中,为标准颜色与校正后颜色的明度差,/>和/>为色度差。
将S5得到的颜色定标系数矩阵作用于待校正色品值后求得其色差,得到色差权重矩阵,并利用下式重新优化颜色定标系数矩阵:
利用以上方法即可求得最优颜色定标系数矩阵C。
S7、将S6得到的定标系数矩阵作用于深空探测彩色相机在轨拍摄图像:
由此反演出图像真实色彩,即为校正后符合人眼标准颜色响应的图像。
本发明未详述部分属于本技术领域的公知技术。以上所述仅为本发明的具体实例而已,并不用于以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。

Claims (1)

1.一种深空探测彩色相机在轨联合颜色定标方法,其特征在于,包括如下步骤:
S1、利用在轨的深空探测彩色相机拍摄1张星体表面图像,并获取其相机系统sRGB颜色空间下的色度值N组,进而归一化为待校正色品值;
S2、当深空探测彩色相机在轨对某一目标区域成像时,利用在轨探测器搭载的光谱仪获得相应的目标光谱数据;
S3、根据S2所获得的目标光谱数据的光谱曲线,结合标准人眼视觉色彩响应函数,解算XYZ颜色空间下的标准人眼响应色度值;
S4、根据颜色空间转换公式,将S3所求得的人眼响应色度值转换为sRGB颜色空间下的色度值,进而归一化为标准色品值;
S5、将S1所获得的待校正色品值的数据与S4所求得的标准色品值的数据做数据拟合,利用最小二乘法初步求取颜色定标系数矩阵;
S6、利用CIE Lab色差公式求解色差,得到色差权重矩阵,并通过色差权重矩阵重新优化颜色定标系数矩阵,直至得到满意的色差结果,得到最优的颜色定标系数矩阵;
S7、将S6得到的颜色定标系数矩阵作用于深空探测彩色相机在轨拍摄图像,反演出图像真实色彩。
CN202310892257.1A 2023-07-20 2023-07-20 一种深空探测彩色相机在轨联合颜色定标方法 Active CN116612199B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202310892257.1A CN116612199B (zh) 2023-07-20 2023-07-20 一种深空探测彩色相机在轨联合颜色定标方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202310892257.1A CN116612199B (zh) 2023-07-20 2023-07-20 一种深空探测彩色相机在轨联合颜色定标方法

Publications (2)

Publication Number Publication Date
CN116612199A true CN116612199A (zh) 2023-08-18
CN116612199B CN116612199B (zh) 2024-02-06

Family

ID=87683992

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202310892257.1A Active CN116612199B (zh) 2023-07-20 2023-07-20 一种深空探测彩色相机在轨联合颜色定标方法

Country Status (1)

Country Link
CN (1) CN116612199B (zh)

Citations (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2006217169A (ja) * 2005-02-02 2006-08-17 Fuji Photo Film Co Ltd 画像処理装置及びカラー撮像装置
US20080240558A1 (en) * 2007-03-19 2008-10-02 Sti Medical Systems, Llc Method of automated image color calibration
CN103079076A (zh) * 2013-01-22 2013-05-01 无锡鸿图微电子技术有限公司 自适应伽玛校正曲线的色彩校正矩阵的产生方法及装置
CN103905803A (zh) * 2014-03-18 2014-07-02 中国科学院国家天文台 一种图像的颜色校正方法及装置
CN109218698A (zh) * 2018-10-19 2019-01-15 浙江大学 一种高容错性的彩色数码相机颜色校正方法
CN109523475A (zh) * 2018-10-29 2019-03-26 中国资源卫星应用中心 一种基于色彩复原及增强的遥感图像处理方法
CN110047042A (zh) * 2019-03-11 2019-07-23 北京空间飞行器总体设计部 一种空间相对导航目标图像局部增强方法及系统
CN110726536A (zh) * 2019-09-25 2020-01-24 宁波永新光学股份有限公司 一种彩色数码反射显微镜颜色校正方法
CN110751611A (zh) * 2019-11-01 2020-02-04 深圳市博盛医疗科技有限公司 一种内窥镜系统的显示器色彩的校准方法及内窥镜系统
CN110830778A (zh) * 2018-08-14 2020-02-21 杭州海康威视数字技术股份有限公司 成像设备色彩校正方法、装置、电子设备及存储介质
CN111047539A (zh) * 2019-12-27 2020-04-21 上海工程技术大学 一种基于光谱反射率重建的织物图像颜色校准算法
US20220003988A1 (en) * 2018-11-12 2022-01-06 Uea Enterprises Ltd System and method for generating a colour filter for modifying the spectral response of a vision system
CN115942128A (zh) * 2022-12-12 2023-04-07 大连理工大学 一种基于异构平台的isp系统设计与实现方法
CN115965541A (zh) * 2022-11-02 2023-04-14 深圳思谋信息科技有限公司 颜色校正方法、装置、电子设备和存储介质

Patent Citations (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2006217169A (ja) * 2005-02-02 2006-08-17 Fuji Photo Film Co Ltd 画像処理装置及びカラー撮像装置
US20080240558A1 (en) * 2007-03-19 2008-10-02 Sti Medical Systems, Llc Method of automated image color calibration
CN103079076A (zh) * 2013-01-22 2013-05-01 无锡鸿图微电子技术有限公司 自适应伽玛校正曲线的色彩校正矩阵的产生方法及装置
CN103905803A (zh) * 2014-03-18 2014-07-02 中国科学院国家天文台 一种图像的颜色校正方法及装置
CN110830778A (zh) * 2018-08-14 2020-02-21 杭州海康威视数字技术股份有限公司 成像设备色彩校正方法、装置、电子设备及存储介质
CN109218698A (zh) * 2018-10-19 2019-01-15 浙江大学 一种高容错性的彩色数码相机颜色校正方法
CN109523475A (zh) * 2018-10-29 2019-03-26 中国资源卫星应用中心 一种基于色彩复原及增强的遥感图像处理方法
US20220003988A1 (en) * 2018-11-12 2022-01-06 Uea Enterprises Ltd System and method for generating a colour filter for modifying the spectral response of a vision system
CN110047042A (zh) * 2019-03-11 2019-07-23 北京空间飞行器总体设计部 一种空间相对导航目标图像局部增强方法及系统
CN110726536A (zh) * 2019-09-25 2020-01-24 宁波永新光学股份有限公司 一种彩色数码反射显微镜颜色校正方法
CN110751611A (zh) * 2019-11-01 2020-02-04 深圳市博盛医疗科技有限公司 一种内窥镜系统的显示器色彩的校准方法及内窥镜系统
CN111047539A (zh) * 2019-12-27 2020-04-21 上海工程技术大学 一种基于光谱反射率重建的织物图像颜色校准算法
CN115965541A (zh) * 2022-11-02 2023-04-14 深圳思谋信息科技有限公司 颜色校正方法、装置、电子设备和存储介质
CN115942128A (zh) * 2022-12-12 2023-04-07 大连理工大学 一种基于异构平台的isp系统设计与实现方法

Non-Patent Citations (6)

* Cited by examiner, † Cited by third party
Title
E.J. SPEYERER等: "Pre-flight and On-orbit Geometric Calibration of the Lunar Reconnaissance Orbiter Camera", 《SPACE SCI REV》, pages 357 - 392 *
NEKTARIOS A. VALOUS等: "Colour calibration of a laboratory computer vision system for quality evaluation of pre-sliced hams", 《MEAT SCIENCE》, pages 132 - 141 *
付小康等: "自适应加权根多项式回归颜色校正算法研究", 《电子测量技术》, vol. 46, no. 8, pages 45 - 50 *
余国彬等: "天问一号中分辨率相机的设计与验证", 《中国科学》, vol. 52, no. 3, pages 1 - 8 *
赵汝进等: "嫦娥三号卫星相机地形图像的彩色校正方法", 《宇航学报》, vol. 37, no. 3, pages 341 - 347 *
黄翔: "基于遗传算法的结构色逆向设计研究", 《中国优秀硕士学位论文全文数据库 信息科技辑》, no. 2023, pages 138 - 1858 *

Also Published As

Publication number Publication date
CN116612199B (zh) 2024-02-06

Similar Documents

Publication Publication Date Title
Helfenstein et al. Galileo photometry of asteroid 243 Ida
Suzuki et al. Initial inflight calibration for Hayabusa2 optical navigation camera (ONC) for science observations of asteroid Ryugu
Tyson et al. A deep CCD survey of 12 high-latitude fields
CN106895851B (zh) 一种光学遥感卫星多ccd多相机统一处理的传感器校正方法
Hanner et al. 1.1. 4 Pioneer 10 Observations of Zodiacal Light Brightness near the Ecliptic: Changes with Heliocentric Distance
Ren et al. A method and results of color calibration for the Chang'e-3 terrain camera and panoramic camera
CN104977024A (zh) 一种日盲紫外遥感相机绝对辐射定标系数在轨修正方法
Dalessandro et al. GeMS/GSAOI photometric and astrometric performance in dense stellar fields
Zhou et al. Photometric Calibration for the Beijing–Arizona Sky Survey and Mayall z-band Legacy Survey
Monnier et al. Diffraction-limited near-IR imaging at Keck reveals asymmetric, time-variable nebula around carbon star CIT 6
CN116612199B (zh) 一种深空探测彩色相机在轨联合颜色定标方法
Zacharias et al. Accurate optical positions of extragalactic radio reference frame sources
Bowen Absolute radiometric calibration of the IKONOS sensor using radiometrically characterized stellar sources
Smith et al. The MVACS surface stereo imager on Mars Polar Lander
CN114894737A (zh) 一种基于红外图像的光谱反射率重建方法
Bonanos et al. Reanalysis of Very Large Telescope Data for M83 with Image Subtraction—Ninefold Increase in Number of Cepheids
Anche et al. Simulation of high-contrast polarimetric observations of debris disks with the Roman Coronagraph Instrument
Clem et al. Secondary Standard Stars for the u′ g′ r′ i′ z′ Photometric System
Payne et al. Satellite photometric error determination
CN109556715A (zh) 一种航空高光谱多航带影像辐射畸变校正方法
Smith et al. Local ugriz Standard Stars in the Chandra Deep Field South
Boyle et al. CCD photometry of the M 67 cluster in the Vilnius photometric system
Zheng et al. The SAGE photometric survey: technical description
Liu et al. Deep-space exploration camera on-orbit white balance calibration based on correlated color temperature interpolation
Tereschenko SPECTROPHOTOMETRIC STANDARDS 8m-10m. II. THE EQUATORIAL ZONA FROM 0h TO 12h.

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