CN112731569B - 基于恒星辐射光谱及通量的全天空成像仪辐射定标方法 - Google Patents
基于恒星辐射光谱及通量的全天空成像仪辐射定标方法 Download PDFInfo
- Publication number
- CN112731569B CN112731569B CN202011547331.9A CN202011547331A CN112731569B CN 112731569 B CN112731569 B CN 112731569B CN 202011547331 A CN202011547331 A CN 202011547331A CN 112731569 B CN112731569 B CN 112731569B
- Authority
- CN
- China
- Prior art keywords
- star
- target
- flux
- determining
- radiation
- 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
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01W—METEOROLOGY
- G01W1/00—Meteorology
- G01W1/18—Testing or calibrating meteorological apparatus
Landscapes
- Environmental & Geological Engineering (AREA)
- Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Atmospheric Sciences (AREA)
- Biodiversity & Conservation Biology (AREA)
- Ecology (AREA)
- Environmental Sciences (AREA)
- Image Processing (AREA)
Abstract
本发明实施例公开了一种基于恒星辐射光谱及通量的全天空成像仪辐射定标方法。该辐射定标方法包括:选取全天空图像中的目标恒星后确定在极光观测站处全天空图像时刻所述目标恒星对应的天顶角,并根据天顶角获取在极光观测谱线上目标恒星对应的天顶角处的大气透射率;计算得到目标恒星在极光观测谱线上的光子辐射通量以及目标恒星对应的第一立体角,并根据光子辐射通量和第一立体角确定目标恒星的恒星辐射亮度;确定目标恒星对应的探测器输出信号值,并根据探测器输出信号值、大气透射率和恒星辐射亮度确定全天空成像仪的辐射系数。本发明实施例的技术方案,以实现不受地域局限快速准确对全天空成像仪进行辐射定标,全天空成像仪辐射定标更加灵活。
Description
技术领域
本发明实施例涉及辐射定标技术领域,尤其涉及一种基于恒星辐射光谱及通量的全天空成像仪辐射定标方法。
背景技术
全天空成像仪是一种重要的地基空间环境光学观测设备,对全天空成像仪进行辐射参数定标,能够获取全天空图像中极光结构在不同观测谱线上的真实激发强度,这既是对比不同观测设备数据中极光激发强度的前提,也是定量研究相关物理过程以及机制的基础。
全天空成像仪在出厂前,主要是通过实验室的定标实验对全天空成像仪的“辐射强度-图像数值”的辐射参数进行定标,而在全天空成像仪的长期使用过程中,由于仪器设备的维护和部件性能的老化等原因,原有的光学性能、探测器增益和偏置等性能将会发生不同程度的变化。因此,为了保证观测数据质量,则需要定期对全天空成像仪重新进行辐射定标。目前,全天空成像仪的辐射定标都是在光学实验室中进行的,使用的定标源是积分球、钨丝灯和朗伯平面等人造辐射源,对于安装在野外观测站点(比如南北极地区)的全天空极光成像仪而言,受站点条件的限制,在实验室重新定标非常困难。
发明内容
本发明实施例提供一种基于恒星辐射光谱及通量的全天空成像仪辐射定标方法,以实现不受地域局限快速准确对全天空成像仪进行辐射定标,全天空成像仪辐射定标更加灵活。
第一方面,本发明实施例提供了一种基于恒星辐射光谱及通量的全天空成像仪辐射定标方法,该辐射定标方法包括:
选取全天空图像中的目标恒星后确定在极光观测站处全天空图像时刻所述目标恒星对应的天顶角,并根据所述天顶角获取在极光观测谱线上所述目标恒星对应的所述天顶角处的大气透射率;
计算得到所述目标恒星在所述极光观测谱线上的光子辐射通量以及所述目标恒星对应的第一立体角,并根据所述光子辐射通量和所述第一立体角确定所述目标恒星的恒星辐射亮度;
确定所述目标恒星对应的探测器输出信号值,并根据所述探测器输出信号值、所述大气透射率和所述恒星辐射亮度确定全天空成像仪的辐射系数。
第二方面,本发明实施例还提供了一种基于恒星辐射光谱及通量的全天空成像仪辐射定标装置,该辐射定标装置包括:
大气透射率获取模块,用于选取全天空图像中的目标恒星后确定在极光观测站处全天空图像时刻所述目标恒星对应的天顶角,并根据所述天顶角获取在极光观测谱线上所述目标恒星对应的所述天顶角处的大气透射率;
恒星辐射亮度确定模块,用于计算得到所述目标恒星在所述极光观测谱线上的光子辐射通量以及所述目标恒星对应的第一立体角,并根据所述光子辐射通量和所述第一立体角确定所述目标恒星的恒星辐射亮度;
辐射定标模块,用于确定所述目标恒星对应的探测器输出信号值,并根据所述探测器输出信号值、所述大气透射率和所述恒星辐射亮度确定全天空成像仪的辐射系数。
第三方面,本发明实施例还提供了一种全天空成像仪,该全天空成像仪包括:
一个或多个处理器;
存储装置,用于存储多个程序,
当所述多个程序中的至少一个被所述一个或多个处理器执行时,使得所述一个或多个处理器实现本发明第一方面实施例所提供的一种基于恒星辐射光谱及通量的全天空成像仪辐射定标方法。
第四方面,本发明实施例还提供了一种计算机可读存储介质,其上存储有计算机程序,该程序被处理器执行时实现本发明第一方面实施例所提供的一种基于恒星辐射光谱及通量的全天空成像仪辐射定标方法。
本发明实施例的技术方案,通过选取全天空图像中的目标恒星后确定在极光观测站处全天空图像时刻所述目标恒星对应的天顶角,并根据所述天顶角获取在极光观测谱线上所述目标恒星对应的所述天顶角处的大气透射率;计算得到所述目标恒星在所述极光观测谱线上的光子辐射通量以及所述目标恒星对应的第一立体角,并根据所述光子辐射通量和所述第一立体角确定所述目标恒星的恒星辐射亮度;确定所述目标恒星对应的探测器输出信号值,并根据所述探测器输出信号值、所述大气透射率和所述恒星辐射亮度确定全天空成像仪的辐射系数。解决了当前处于野外观测站点的全天空极光成像仪进行辐射定标困难的问题,以实现不受地域局限快速准确对全天空成像仪进行辐射定标,全天空成像仪辐射定标更加灵活。
附图说明
图1是本发明实施例一提供的一种基于恒星辐射光谱及通量的全天空成像仪辐射定标方法的流程示意图;
图2A是本发明实施例二提供的一种基于恒星辐射光谱及通量的全天空成像仪辐射定标方法的流程示意图;
图2B是本发明实施例提供的黄河站拍摄的全天空图像以及全天空图像对应的星空图;
图2C是本发明实施例提供的模拟在极光观测谱线上恒星对应天顶角处的大气透射率的示意图;
图2D是本发明实施例提供的五车二星对应的光谱型曲线图;
图2E是本发明实施例提供的确定恒星在全天空图像中的位置原理示意图;
图2F是本发明实施例提供的全天空图像中五车二星对应的星点邻域示意图;
图2G是本发明实施例提供的目标恒星对应的第一立体角的原理示意图;
图2H是本发明实施例提供的利用每张全天空图像单独进行定标的误差的结果示意图;
图3是本发明实施例三提供的一种基于恒星辐射光谱及通量的全天空成像仪辐射定标装置的结构示意图;
图4是本发明实施例四提供的一种全天空成像仪的硬件结构示意图。
具体实施方式
为了使本发明的目的、技术方案和优点更加清楚,下面结合附图对本发明具体实施例作进一步的详细描述。可以理解的是,此处所描述的具体实施例仅仅用于解释本发明,而非对本发明的限定。
另外还需要说明的是,为了便于描述,附图中仅示出了与本发明相关的部分而非全部内容。在更加详细地讨论示例性实施例之前应当提到的是,一些示例性实施例被描述成作为流程图描绘的处理或方法。虽然流程图将各项操作(或步骤)描述成顺序的处理,但是其中的许多操作可以被并行地、并发地或者同时实施。此外,各项操作的顺序可以被重新安排。当其操作完成时所述处理可以被终止,但是还可以具有未包括在附图中的附加步骤。所述处理可以对应于方法、函数、规程、子例程、子程序等等。
实施例一
图1为本发明实施例一提供的一种基于恒星辐射光谱及通量的全天空成像仪辐射定标方法的流程示意图,本实施例可适用于对任意地点的全天空成像仪进行快速而准确的辐射定标的情况,该方法可以由基于恒星辐射光谱及通量的全天空成像仪辐射定标装置来执行,该装置可以通过软件和/或硬件的形式实现。具体包括如下步骤:
S110、选取全天空图像中的目标恒星后确定在极光观测站处全天空图像时刻所述目标恒星对应的天顶角,并根据所述天顶角获取在极光观测谱线上所述目标恒星对应的所述天顶角处的大气透射率。
其中,通过全天空成像仪获取一组全天空图像,一组全天空图像包括一张、两张或多张全天空图像,具体选取的全天空图像的数量可以由本领域技术人员根据实际辐射定标需求进行选择设置,本实施例对此不作任何限制。
进一步的,全天空图像中恒星清晰且恒星周围背景干净,即全天空图像中恒星周围无云且无极光遮挡,以作为本实施例所需的全天空图像。
在本实施例中,目标恒星可以选择全天空图像中较明显且天顶角较低的恒星,可选的,全天空图像中较明显可以设定为目标恒星中心对应的探测器信号输出值与恒星背景对应的探测器信号输出值的比值大于15,天顶角较低可以选择天顶角角度值小于65°,上述亮度或天顶角的具体数值并非对本实施例进行任何限制,而由本领域技术人员根据实际设置需要进行选择,本实施例在此仅对目标恒星的选择进行解释说明。
具体的,在选取全天空图像中的目标恒星后,利用星空软件Stellarium确定目标恒星在极光观测台站处与全天空图像拍摄时间对应的天顶角,即Stellarium中的地点设置为极光观测台站处的地理位置,时间设置为全天空图像拍摄时间。可以理解的是,天顶角的获取不限于利用星空软件Stellarium进行获取,也可以通过其他天顶角获取的计算机软件或设备进行获取,本实施例对天顶角的具体获取不作任何限制。
进一步的,基于天顶角利用基于大气传输软件MODTRAN的界面化软件PcModWin模拟,得到在极光观测谱线上所述目标恒星对应的所述天顶角处的大气透射率。可以理解的是,大气透射率的获取不限于利用基于大气传输软件MODTRAN的界面化软件PcModWin模拟进行获取,也可以通过其他大气透射率获取的计算机软件或设备进行获取,本实施例对大气透射率的具体获取不作任何限制。
S120、计算得到所述目标恒星在所述极光观测谱线上的光子辐射通量fphoton以及所述目标恒星对应的第一立体角wstar,并根据所述光子辐射通量和所述第一立体角确定所述目标恒星的恒星辐射亮度Istar。
具体的,计算得到所述目标恒星在所述极光观测谱线上的光子辐射通量,包括:确定所述极光观测谱线对应的滤波片波段以及所述滤波片波段内的通量零点zp(fv),并根据所述通量零点计算有效波长λeff处的有效波长辐射通量根据所述有效波长辐射通量以及所述目标恒星对应的光谱型确定所述目标恒星在所述极光观测谱线上的辐射通量f,并根据所述辐射通量确定所述光子辐射通量。
进一步的,在根据所述辐射通量确定所述光子辐射通量之前,还包括:按照公式:计算得到所述极光观测谱线处的光子能量,其中,e为光子能量;h为普朗克常数;c为真空中的光速;λ为极光观测谱线波长;相应的,根据所述辐射通量确定所述光子辐射通量,包括:根据所述辐射通量和所述光子能量按照公式:计算得到所述光子辐射通量;其中,fphoton为所述光子辐射通量;f为所述辐射通量;e为所述光子能量。
在上述实施例的基础上,在计算得到所述目标恒星在所述极光观测谱线上的光子辐射通量以及所述目标恒星对应的第一立体角之前,还包括:提取所述目标恒星在全天空图像中的星点邻域,并根据所述星点邻域确定所述目标恒星在全天空图像中对应的区域范围。相应的,计算得到所述目标恒星对应的第一立体角,包括:获取所述目标恒星在所述区域范围内各像素点到天顶的最近距离Rmin、最远距离Rmax以及所述最近距离对应的天顶角和所述最远距离对应的天顶角之间的圆环对应的第二立体角Δw;根据所述最近距离Rmin、最远距离Rmax以及第二立体角Δw确定每个像素点对应的第三立体角wpixel,并根据所述第三立体角wpixel确定所述目标恒星对应的第一立体角,即按照公式: 计算得到每个像素点对应的第三立体角wpixel;其中,wpixel为每个像素点对应的第三立体角;Δw为所述最近距离对应的天顶角和所述最远距离对应的天顶角之间的圆环对应的第二立体角;Rmin为所述目标恒星在所述区域范围内各像素点到天顶的最近距离;Rmax为所述目标恒星在所述区域范围内各像素点到天顶的最远距离;相应的,根据所述第三立体角wpixel确定所述目标恒星对应的第一立体角,包括:按照公式:wstar=wpixel·m,计算得到所述目标恒星对应的第一立体角;其中,wstar为所述目标恒星对应的第一立体角;m为所述目标恒星的区域范围内的像素个数。
S130、确定所述目标恒星对应的探测器输出信号值Nstar,并根据所述探测器输出信号值、所述大气透射率和所述恒星辐射亮度确定全天空成像仪的辐射系数。
其中,目标恒星对应的探测器输出信号值通过所述目标恒星在全天空图像中的星点邻域以及所述目标恒星在全天空图像中对应的区域范围确定。
在上述实施例的基础上,根据所述探测器输出信号值、所述大气透射率和所述恒星辐射亮度确定全天空成像仪的辐射系数,包括:按照公式:计算得到全天空成像仪的辐射系数;其中,K为所述辐射系数;Ttra为所述大气透射率;Istar为所述恒星辐射亮度;Nstar为所述目标恒星对应的探测器输出信号值。
本发明实施例的技术方案,通过选取全天空图像中的目标恒星后确定在极光观测站处全天空图像时刻所述目标恒星对应的天顶角,并根据所述天顶角获取在极光观测谱线上所述目标恒星对应的所述天顶角处的大气透射率;计算得到所述目标恒星在所述极光观测谱线上的光子辐射通量以及所述目标恒星对应的第一立体角,并根据所述光子辐射通量和所述第一立体角确定所述目标恒星的恒星辐射亮度;确定所述目标恒星对应的探测器输出信号值,并根据所述探测器输出信号值、所述大气透射率和所述恒星辐射亮度确定全天空成像仪的辐射系数。解决了当前处于野外观测站点的全天空极光成像仪进行辐射定标困难的问题,以实现不受地域局限快速准确对全天空成像仪进行辐射定标,全天空成像仪辐射定标更加灵活。
实施例二
图2A为本发明实施例二提供的一种基于恒星辐射光谱及通量的全天空成像仪辐射定标方法的流程示意图,本实施例以上述实施例为基础进行优化。
相应的,本实施例的方法具体包括:
S210、选取全天空图像中的目标恒星后确定在极光观测站处全天空图像时刻所述目标恒星对应的天顶角,并根据所述天顶角获取在极光观测谱线上所述目标恒星对应的所述天顶角处的大气透射率。
具体的,选取一组恒星清晰且恒星周围背景干净的全天空图像,并在选取后的全天空图像中选择较明显且天顶角较低的恒星,作为目标恒星,利用星空软件Stellarium确定目标恒星在极光观测台站处全天空图像时刻对应的天顶角。示例性的,以黄河站557.7nm波段的全天空极光成像仪为例,图2B是本发明实施例提供的黄河站拍摄的全天空图像以及全天空图像对应的星空图,如图2B所示,全天空图像的拍摄时间为2004年2月25日00:05:32UT时,即全天空图像时刻,图2B中对恒星Kochab北极二、Dubhe天枢、Capella五车二、Pollux北河三、Castor北河二以及Mirfak天船三等进行了标记,左图为黄河站全天空图像,右图为全天空图像对应的星空图。
进一步的,图2C是本发明实施例提供的模拟在极光观测谱线上恒星对应天顶角处的大气透射率的示意图,参见图2C,横坐标为波长微米(wavelength microns),纵坐标为透射率(Transmission),图2C中示出了利用基于大气传输软件MODTRAN的界面化软件PcModWin模拟,得到极光观测谱线上目标恒星对应天顶角为59.6428°处的大气透射率的模拟结果。
S220、确定所述极光观测谱线对应的滤波片波段以及所述滤波片波段内的通量零点,并根据所述通量零点计算有效波长处的有效波长辐射通量。
具体的,通过下表1可以确定所述极光观测谱线对应的滤波片波段以及所述滤波片波段内的通量零点,表中滤波片波段包括U、B、V、R、I、J、H、K、Kp、L,每个滤波片波段对应的有效波长λeff与滤波片波段对应的通量零点zp(fλ)一一对应,示例性的,通过查询表1可以得到U波段的有效波长λeff为0.366um,U波段的通量零点zp(fλ)为-0.152。
表1不同滤波片波段的有效波长和对应的通量零点
波段 | U | B | V | R | I | J | H | K | Kp | L |
λ<sub>eff</sub>/um | 0.366 | 0.438 | 0.545 | 0.641 | 0.798 | 1.22 | 1.63 | 2.19 | 2.12 | 3.45 |
zp(f<sub>λ</sub>) | -0.152 | -0.602 | 0.000 | 0.555 | 1.271 | 2.655 | 3.760 | 4.906 | 4.780 | 6.775 |
进一步的,按照下述公式计算得到有效波长处的有效波长辐射通量为:
S230、根据所述有效波长辐射通量以及所述目标恒星对应的光谱型确定所述目标恒星在所述极光观测谱线上的辐射通量,并根据所述辐射通量确定所述光子辐射通量。
根据恒星的光谱型,计算目标恒星在极光观测谱线上的辐射通量f为:
示例性的,图2D是本发明实施例提供的五车二星对应的光谱型曲线图,参见图2D,横坐标为波长(wavelength,单位),纵坐标为辐射通量相对值(relative value),例如V段有效波长为其对应的辐射通量相对值为1.02。
进一步的,在根据所述辐射通量确定所述光子辐射通量之前,还包括:按照公式:计算得到所述极光观测谱线处的光子能量e;其中,e为光子能量,单位为J;h为普朗克常数,取值为6.62607×10-34J·sec;c为真空中的光速,取值为2.9979×108m·sec-1;λ为极光观测谱线波长,单位为nm。
根据所述辐射通量和所述光子能量按照公式:计算得到所述目标恒星在极光观测谱线处的光子辐射通量;其中,fphoton为所述光子辐射通量,单位为photons·sec-1·m-2·nm-1;f为所述辐射通量;e为所述光子能量。
S240、提取所述目标恒星在全天空图像中的星点邻域。
图2E是本发明实施例提供的确定恒星在全天空图像中的位置原理示意图,参见图2E所示,恒星在全天空图像中的轨迹近似为一个圆形,假设该圆的圆心为C,半径为r,则对于任意时刻t,恒星在全天空图像坐标系下的坐标为:
其中,①xC和yC分别为圆心C的横坐标和纵坐标,可以利用任意三个时刻t0、t1和t2恒星在全天空图像中的位置坐标M(x0,y0)、N(x1,y1)和P(x2,y2)计算获得:
②r为圆C的半径,单位为像素:
β为t0时刻恒星与圆心C连线与X轴正方向的夹角:
图2F是本发明实施例提供的全天空图像中五车二星对应的星点邻域示意图。参见图2F所示,全天空图像时刻为2004年2月25日00:05:52时,黄河站557.7nm波段全天空图像中的五车二对应的星点邻域,以坐标(xt,yt)为中心,示例性的,在其周围11×11像素的范围内寻找亮度值最大的像素,将该像素作为目标恒星的实际位置,则目标恒星在全天空图像中的实际坐标满足以下要求:
其中,(xm,ym)为目标恒星在t时刻的实际位置坐标;N(xm,ym)和N(x,y)分别为全天空图像中坐标为(xm,ym)和(x,y)处像素对应的探测器输出信号值。
进一步的,确定目标恒星对应的星点邻域,大小为11×11像素:
其中,x和y分别为星点邻域内像素点的坐标。
S250、根据所述星点邻域确定所述目标恒星在全天空图像中对应的区域范围。
具体的,利用区域生长算法根据所述星点邻域确定目标恒星在全天空图像中对应的区域范围。
以目标恒星在t时刻的实际位置坐标(xm,ym)为星点中心,控制目标恒星区域范围向星点中心的上、下、左、右方向上进行生长,继续参见图2F所示,以星点中心(394,288)向星点中心的上、下、左、右方向上进行生长,得到浅灰色区域为目标恒星对应的区域范围。
可以理解的是,恒星区域范围在这4个方向上的生长方式类似,以右方向为例说明如下。在右方向上,控制恒星区域在星点中心所在的行上自星点中心向右进行生长,将满足生长规则的像素点判定为属于恒星区域范围。具体的生长规则为:
其中,N(x,y)为全天空图像中坐标(x,y)处的探测器输出信号值,g1为预设参数,为阈值,设置为星点邻域最外层像素点处探测器输出信号值的均值,即:
进一步的,控制恒星区域范围继续向星点中心的左上、左下、右上、右下方向上进行生长,恒星区域范围在这4个方向上的生长方式类似,以右上方向为例说明如下。在右上方向上,恒星区域从星点中心向上逐行进行生长,每一行内,自星点中心所在的列从左向右进行生长。对于每一个像素点,设置可信度P,并初始化为0。当该像素点的探测器输出信号值大于阈值T时,将该像素点的信号值N(x,y)与该像素点8邻域像素点的探测器输出信号值N(x+Δr,y+Δu)进行比较,当满足下述两式其中一个:
ds≥dnanf N(x,y)≤N(x+Δr,y+Δu)-g2 Δr∈[-1,1],Δu∈[-1,1];或者,
ds≤dnand N(x,y)≥N(x+Δr,y+Δu)+g2 Δr∈[-1,1],Δu∈[-1,1]
则,将可信度P加1。
其中,ds为像素点(x,y)到星点中心(xm,ym)的距离为:
dn为8邻域内的像素点(x+Δr,y+Δu)到星点中心的距离为:
对于像素点(x,y),遍历完其8邻域内像素点后,当可信度P大于可信度阈值PT时,即判定像素点(x,y)属于恒星区域范围。综上,最终确定的恒星区域范围如图2F中的浅灰色区域所示。
S260、获取所述目标恒星在所述区域范围内各像素点到天顶的最近距离、最远距离以及所述最近距离对应的天顶角和所述最远距离对应的天顶角之间的圆环对应的第二立体角。
对全天空极光成像仪进行几何定标,得到全天空图像中天顶的位置坐标以及成像半径与天顶角之间的关系。示例性的,对于黄河站557.7nm波段的全天空极光成像仪,全天空图像中的天顶位置坐标为(255.5925,254.8603),天顶角θ(单位度)与成像半径R(单位像素)之间的关系为:
θ=1.867×10(-6)×R3-4.250×10(-4)×R2+0.3544×R-1.062
图2G是本发明实施例提供的目标恒星对应的第一立体角的原理示意图。参见图2G所示,恒星区域范围内各像素点到天顶的最近距离为Rmin(图2G左图标示为rmin),对应的天顶角为θmin,对应的立体角为wmin;恒星区域范围内各像素点到天顶的最远距离为Rmax(图2G左图标示为rmax),对应的天顶角为θmax,对应的立体角为wmax。
wmin和wmax的计算公式为:
则,所述最近距离对应的天顶角θmin和所述最远距离θmax对应的天顶角之间的圆环对应的第二立体角Δw为:
Δw=wmax-wmin
S270、根据所述最近距离、最远距离以及第二立体角确定每个像素点对应的第三立体角,并根据所述第三立体角确定所述目标恒星对应的第一立体角。
具体的,在全天空图像中θmin与θmax相差很小,假设圆环中每个像素对应的立体角都相等,则每个像素对应的第三立体角wpixel为:
其中,wpixel为每个像素点对应的第三立体角;Δw为所述最近距离对应的天顶角和所述最远距离对应的天顶角之间的圆环对应的第二立体角;Rmin为所述目标恒星在所述区域范围内各像素点到天顶的最近距离;Rmax为所述目标恒星在所述区域范围内各像素点到天顶的最远距离;
则,目标恒星对应的第一立体角wstar为:
wstar=wpixel·m
其中,wstar为所述目标恒星对应的第一立体角,单位为sr;m为所述目标恒星的区域范围内的像素个数。
S280、根据所述光子辐射通量和所述第一立体角确定所述目标恒星的恒星辐射亮度。
S290、确定所述目标恒星对应的探测器输出信号值,并根据所述探测器输出信号值、所述大气透射率和所述恒星辐射亮度确定全天空成像仪的辐射系数。
具体的,确定所述目标恒星对应的探测器输出信号值Nstar(Nstar的单位为count)为:
背景的探测器输出信号值选取于图2F中环绕恒星区域范围的浅灰色区域的深色区域,即为目标恒星确定的区域范围向8邻域方向分别生长一个像素得到。其中,背景的探测器输出信号值Nbg的计算公式为:
其中,Sbg为恒星背景像素的坐标集合,定义为恒星区域范围内各像素点8邻域像素中不在恒星区域范围内的像素点,范围如图2F中的深灰色部分所示;n为背景像素的个数。
进一步的,根据所述探测器输出信号值、所述大气透射率和所述恒星辐射亮度计算全天空极光成像仪的辐射系数为
其中,K为辐射系数,单位为Rayleigh·nm-1·count-1;Ttra为大气透射率;Istar为所述恒星辐射亮度;Nstar为所述探测器输出信号值。
在上述实施例的基础上,图2H是本发明实施例提供的利用每张全天空图像单独进行定标的误差的结果示意图。参见图2H所示,横坐标为全天空图像序号(Image number),纵坐标为误差百分比(error percentage),对北极黄河站557.7nm波段的全天空极光成像仪拍摄全天空图像,利用2004年2月24日22:56:52UT-23:08:22UT时间段内70张全天空图像进行辐射定标,得到的平均辐射定标结果为1.0884Rayleigh·nm-1·count-1,相对于实验室的定标结果,误差为-0.23%,可以通过实验结果证明相较于现有技术的优势。
本发明实施例的技术方案,将恒星作为定标源,根据恒星的星等和光谱型计算恒星对应的光子辐射通量,再根据恒星在全天空图像中的像点确定恒星对应的立体角后计算恒星在极光观测波段的亮度,再确定恒星对应的探测器输出信号值,实现对全天空极光成像仪的辐射定标,进而实现对任意地点的全天空极光成像仪进行快速而准确的辐射定标。
实施例三
图3为本发明实施例三提供的一种基于恒星辐射光谱及通量的全天空成像仪辐射定标装置的结构示意图,本实施例可适用于对任意地点的全天空极光成像仪进行快速而准确的辐射定标的情况。
如图3所示,所述辐射定标装置包括:大气透射率获取模块310、恒星辐射亮度确定模块320和辐射定标模块330,其中:
大气透射率获取模块310,用于选取全天空图像中的目标恒星后确定在极光观测站处全天空图像时刻所述目标恒星对应的天顶角,并根据所述天顶角获取在极光观测谱线上所述目标恒星对应的所述天顶角处的大气透射率;
恒星辐射亮度确定模块320,用于计算得到所述目标恒星在所述极光观测谱线上的光子辐射通量以及所述目标恒星对应的第一立体角,并根据所述光子辐射通量和所述第一立体角确定所述目标恒星的恒星辐射亮度;
辐射定标模块330,用于确定所述目标恒星对应的探测器输出信号值,并根据所述探测器输出信号值、所述大气透射率和所述恒星辐射亮度确定全天空成像仪的辐射系数。
本实施例的基于恒星辐射光谱及通量的全天空成像仪辐射定标装置,通过选取全天空图像中的目标恒星后确定在极光观测站处全天空图像时刻所述目标恒星对应的天顶角,并获取在极光观测谱线上所述目标恒星对应的所述天顶角处的大气透射率;计算得到所述目标恒星在所述极光观测谱线上的光子辐射通量以及所述目标恒星对应的第一立体角,并根据所述光子辐射通量和所述第一立体角确定所述目标恒星的恒星辐射亮度;确定所述目标恒星对应的探测器输出信号值,并根据所述探测器输出信号值、所述大气透射率和所述恒星辐射亮度确定全天空成像仪的辐射系数。解决了当前处于野外观测站点的全天空极光成像仪进行辐射定标困难的问题,以实现不受地域局限快速准确对全天空成像仪进行辐射定标,全天空成像仪辐射定标更加灵活。
在上述各实施例的基础上,计算得到所述目标恒星在所述极光观测谱线上的光子辐射通量,包括:
确定所述极光观测谱线对应的滤波片波段以及所述滤波片波段内的通量零点,并根据所述通量零点计算有效波长处的有效波长辐射通量;
根据所述有效波长辐射通量以及所述目标恒星对应的光谱型确定所述目标恒星在所述极光观测谱线上的辐射通量,并根据所述辐射通量确定所述光子辐射通量。
在上述各实施例的基础上,在根据所述辐射通量确定所述光子辐射通量之前,还包括:
其中,e为光子能量;h为普朗克常数;c为真空中的光速;λ为极光观测谱线波长;
相应的,根据所述辐射通量确定所述光子辐射通量,包括:
其中,fphoton为所述光子辐射通量;f为所述辐射能通量;e为所述光子能量。
在上述各实施例的基础上,在计算得到所述目标恒星在所述极光观测谱线上的光子辐射通量以及所述目标恒星对应的第一立体角之前,还包括:
提取所述目标恒星在全天空图像中的星点邻域,并根据所述星点邻域确定所述目标恒星在全天空图像中对应的区域范围。
在上述各实施例的基础上,计算得到所述目标恒星对应的第一立体角,包括:
获取所述目标恒星在所述区域范围内各像素点到天顶的最近距离、最远距离以及所述最近距离对应的天顶角和所述最远距离对应的天顶角之间的圆环对应的第二立体角;
根据所述最近距离、最远距离以及第二立体角确定每个像素点对应的第三立体角,并根据所述第三立体角确定所述目标恒星对应的第一立体角。
在上述各实施例的基础上,根据所述最近距离、最远距离以及第二立体角确定每个像素点对应的第三立体角,包括:
其中,wpixel为每个像素点对应的第三立体角;Δw为所述最近距离对应的天顶角和所述最远距离对应的天顶角之间的圆环对应的第二立体角;Rmin为所述目标恒星在所述区域范围内各像素点到天顶的最近距离;Rmax为所述目标恒星在所述区域范围内各像素点到天顶的最远距离;
相应的,根据所述第三立体角确定所述目标恒星对应的第一立体角,包括:
按照公式:wstar=wpixel·m,计算得到所述目标恒星对应的第一立体角;
其中,wstar为所述目标恒星对应的第一立体角;m为所述目标恒星的区域范围内的像素个数。
在上述各实施例的基础上,根据所述探测器输出信号值、所述大气透射率和所述恒星辐射亮度确定全天空成像仪的辐射系数,包括:
其中,K为所述辐射系数;Ttra为所述大气透射率;Istar为所述恒星辐射亮度;Nstar为所述探测器输出信号值。
上述各实施例所提供的基于恒星辐射光谱及通量的全天空成像仪辐射定标装置可执行本发明任意实施例所提供的基于恒星辐射光谱及通量的全天空成像仪辐射定标方法,具备执行基于恒星辐射光谱及通量的全天空成像仪辐射定标方法相应的功能模块和有益效果。
实施例四
图4为本发明实施例四提供的一种全天空成像仪的结构示意图,如图4所示,该全天空成像仪包括处理器410、存储器420、输入装置430和输出装置440;全天空成像仪中处理器410的数量可以是一个或多个,图4中以一个处理器410为例;全天空成像仪中的处理器410、存储器420、输入装置430和输出装置440可以通过总线或其他方式连接,图4中以通过总线连接为例。
存储器420作为一种计算机可读存储介质,可用于存储软件程序、计算机可执行程序以及模块,如本发明实施例中的基于恒星辐射光谱及通量的全天空成像仪辐射定标方法对应的程序指令/模块(例如,基于恒星辐射光谱及通量的全天空成像仪辐射定标装置中的大气透射率获取模块310、恒星辐射亮度确定模块320和辐射定标模块330)。处理器410通过运行存储在存储器420中的软件程序、指令以及模块,从而执行全天空成像仪的各种功能应用以及数据处理,即实现上述的基于恒星辐射光谱及通量的全天空成像仪辐射定标方法。
存储器420可主要包括存储程序区和存储数据区,其中,存储程序区可存储操作***、至少一个功能所需的应用程序;存储数据区可存储根据终端的使用所创建的数据等。此外,存储器420可以包括高速随机存取存储器,还可以包括非易失性存储器,例如至少一个磁盘存储器件、闪存器件、或其他非易失性固态存储器件。在一些实例中,存储器420可进一步包括相对于处理器410远程设置的存储器,这些远程存储器可以通过网络连接至全天空成像仪。上述网络的实例包括但不限于互联网、企业内部网、局域网、移动通信网及其组合。
输入装置430可用于接收输入的数字或字符信息,以及产生与全天空成像仪的用户设置以及功能控制有关的键信号输入。输出装置440可包括显示屏等显示设备。
实施例五
本发明实施例五还提供一种包含计算机可执行指令的存储介质,所述计算机可执行指令在由计算机处理器执行时用于执行一种基于恒星辐射光谱及通量的全天空成像仪辐射定标方法,该辐射定标方法包括:
选取全天空图像中的目标恒星后确定在极光观测站处全天空图像时刻所述目标恒星对应的天顶角,并根据所述天顶角获取在极光观测谱线上所述目标恒星对应的所述天顶角处的大气透射率;
计算得到所述目标恒星在所述极光观测谱线上的光子辐射通量以及所述目标恒星对应的第一立体角,并根据所述光子辐射通量和所述第一立体角确定所述目标恒星的恒星辐射亮度;
确定所述目标恒星对应的探测器输出信号值,并根据所述探测器输出信号值、所述大气透射率和所述恒星辐射亮度确定全天空成像仪的辐射系数。
当然,本发明实施例所提供的一种包含计算机可执行指令的存储介质,其计算机可执行指令不限于如上所述的方法操作,还可以执行本发明任意实施例所提供的基于恒星辐射光谱及通量的全天空成像仪辐射定标方法中的相关操作。
通过以上关于实施方式的描述,所属领域的技术人员可以清楚地了解到,本发明可借助软件及必需的通用硬件来实现,当然也可以通过硬件实现,但很多情况下前者是更佳的实施方式。基于这样的理解,本发明的技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来,该计算机软件产品可以存储在计算机可读存储介质中,如计算机的软盘、只读存储器(Read-Only Memory,ROM)、随机存取存储器(RandomAccess Memory,RAM)、闪存(FLASH)、硬盘或光盘等,包括若干指令用以使得一台计算机设备(可以是个人计算机,服务器,或者网络设备等)执行本发明各个实施例所述的方法。
值得注意的是,上述基于恒星辐射光谱及通量的全天空成像仪辐射定标装置的实施例中,所包括的各个单元和模块只是按照功能逻辑进行划分的,但并不局限于上述的划分,只要能够实现相应的功能即可;另外,各功能单元的具体名称也只是为了便于相互区分,并不用于限制本发明的保护范围。
注意,上述仅为本发明的较佳实施例及所运用技术原理。本领域技术人员会理解,本发明不限于这里所述的特定实施例,对本领域技术人员来说能够进行各种明显的变化、重新调整和替代而不会脱离本发明的保护范围。因此,虽然通过以上实施例对本发明进行了较为详细的说明,但是本发明不仅仅限于以上实施例,在不脱离本发明构思的情况下,还可以包括更多其他等效实施例,而本发明的范围由所附的权利要求范围决定。
Claims (10)
1.一种基于恒星辐射光谱及通量的全天空成像仪辐射定标方法,其特征在于,包括:
选取全天空图像中的目标恒星后确定在极光观测站处全天空图像时刻所述目标恒星对应的天顶角,并根据所述天顶角获取在极光观测谱线上所述目标恒星对应的所述天顶角处的大气透射率;
计算得到所述目标恒星在所述极光观测谱线上的光子辐射通量以及所述目标恒星对应的第一立体角,并根据所述光子辐射通量和所述第一立体角确定所述目标恒星的恒星辐射亮度;
确定所述目标恒星对应的探测器输出信号值,并根据所述探测器输出信号值、所述大气透射率和所述恒星辐射亮度确定全天空成像仪的辐射系数。
2.根据权利要求1所述的方法,其特征在于,计算得到所述目标恒星在所述极光观测谱线上的光子辐射通量,包括:
确定所述极光观测谱线对应的滤波片波段以及所述滤波片波段内的通量零点,并根据所述通量零点计算有效波长处的有效波长辐射通量;
根据所述有效波长辐射通量以及所述目标恒星对应的光谱型确定所述目标恒星在所述极光观测谱线上的辐射通量,并根据所述辐射通量确定所述光子辐射通量。
4.根据权利要求1所述的方法,其特征在于,在计算得到所述目标恒星在所述极光观测谱线上的光子辐射通量以及所述目标恒星对应的第一立体角之前,还包括:
提取所述目标恒星在全天空图像中的星点邻域,并根据所述星点邻域确定所述目标恒星在全天空图像中对应的区域范围。
5.根据权利要求4所述的方法,其特征在于,计算得到所述目标恒星对应的第一立体角,包括:
获取所述目标恒星在所述区域范围内各像素点到天顶的最近距离、最远距离以及所述最近距离对应的天顶角和所述最远距离对应的天顶角之间的圆环对应的第二立体角;
根据所述最近距离、最远距离以及第二立体角确定每个像素点对应的第三立体角,并根据所述第三立体角确定所述目标恒星对应的第一立体角。
6.根据权利要求5所述的方法,其特征在于,根据所述最近距离、最远距离以及第二立体角确定每个像素点对应的第三立体角,包括:
其中,wpixel为每个像素点对应的第三立体角;Δw为所述最近距离对应的天顶角和所述最远距离对应的天顶角之间的圆环对应的第二立体角;Rmin为所述目标恒星在所述区域范围内各像素点到天顶的最近距离;Rmax为所述目标恒星在所述区域范围内各像素点到天顶的最远距离;
相应的,根据所述第三立体角确定所述目标恒星对应的第一立体角,包括:
按照公式:wstar=wpixel·m,计算得到所述目标恒星对应的第一立体角;
其中,wstar为所述目标恒星对应的第一立体角;m为所述目标恒星的区域范围内的像素个数。
8.一种基于恒星辐射光谱及通量的全天空成像仪辐射定标装置,其特征在于,包括:
大气透射率获取模块,用于选取全天空图像中的目标恒星后确定在极光观测站处全天空图像时刻所述目标恒星对应的天顶角,并根据所述天顶角获取在极光观测谱线上所述目标恒星对应的所述天顶角处的大气透射率;
恒星辐射亮度确定模块,用于计算得到所述目标恒星在所述极光观测谱线上的光子辐射通量以及所述目标恒星对应的第一立体角,并根据所述光子辐射通量和所述第一立体角确定所述目标恒星的恒星辐射亮度;
辐射定标模块,用于确定所述目标恒星对应的探测器输出信号值,并根据所述探测器输出信号值、所述大气透射率和所述恒星辐射亮度确定全天空成像仪的辐射系数。
9.一种全天空成像仪,其特征在于,所述全天空成像仪包括:
一个或多个处理器;
存储装置,用于存储一个或多个程序;
当所述一个或多个程序被所述一个或多个处理器执行,使得所述一个或多个处理器实现如权利要求1-7中任一项所述的基于恒星辐射光谱及通量的全天空成像仪辐射定标方法。
10.一种计算机可读存储介质,其上存储有计算机程序,其特征在于,该程序被处理器执行时实现如权利要求1-7中任一项所述的基于恒星辐射光谱及通量的全天空成像仪辐射定标方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011547331.9A CN112731569B (zh) | 2020-12-24 | 2020-12-24 | 基于恒星辐射光谱及通量的全天空成像仪辐射定标方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011547331.9A CN112731569B (zh) | 2020-12-24 | 2020-12-24 | 基于恒星辐射光谱及通量的全天空成像仪辐射定标方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112731569A CN112731569A (zh) | 2021-04-30 |
CN112731569B true CN112731569B (zh) | 2022-07-12 |
Family
ID=75605575
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011547331.9A Active CN112731569B (zh) | 2020-12-24 | 2020-12-24 | 基于恒星辐射光谱及通量的全天空成像仪辐射定标方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112731569B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112799144B (zh) * | 2020-12-30 | 2024-01-16 | 中国极地研究中心 | 一种基于梯度下降法的全天空成像仪的几何定标方法 |
CN115619716B (zh) * | 2022-09-20 | 2023-05-02 | 中国人民解放军61540部队 | 一种夜间极光移动参数计算方法、***、电子设备及介质 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102901516A (zh) * | 2012-09-29 | 2013-01-30 | 航天恒星科技有限公司 | 一种基于绝对辐射定标的多光谱影像辐射校正方法 |
CN105737980A (zh) * | 2016-02-14 | 2016-07-06 | 国家卫星气象中心 | 辐射定标方法及装置 |
CN109738973A (zh) * | 2019-02-14 | 2019-05-10 | 浙江中控太阳能技术有限公司 | 一种基于太阳方位的全天空成像仪标定方法 |
WO2019244510A1 (ja) * | 2018-06-19 | 2019-12-26 | 古野電気株式会社 | 雲観測装置、雲観測システム、雲観測方法、及びプログラム |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9565377B2 (en) * | 2013-04-30 | 2017-02-07 | International Business Machines Corporation | Multifunctional sky camera system for total sky imaging and spectral radiance measurement |
-
2020
- 2020-12-24 CN CN202011547331.9A patent/CN112731569B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102901516A (zh) * | 2012-09-29 | 2013-01-30 | 航天恒星科技有限公司 | 一种基于绝对辐射定标的多光谱影像辐射校正方法 |
CN105737980A (zh) * | 2016-02-14 | 2016-07-06 | 国家卫星气象中心 | 辐射定标方法及装置 |
WO2019244510A1 (ja) * | 2018-06-19 | 2019-12-26 | 古野電気株式会社 | 雲観測装置、雲観測システム、雲観測方法、及びプログラム |
CN109738973A (zh) * | 2019-02-14 | 2019-05-10 | 浙江中控太阳能技术有限公司 | 一种基于太阳方位的全天空成像仪标定方法 |
Non-Patent Citations (2)
Title |
---|
Calibration of the AIRS Microwave Instruments;Bjorn H. Lambrigtsen;《IEEE Transactions on Geoscience and Remote Sensing》;20031231;全文 * |
基于星点的全天空极光图像参数标定及其验证;张会朋 等;《地球物理学报》;20200229;全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN112731569A (zh) | 2021-04-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112731569B (zh) | 基于恒星辐射光谱及通量的全天空成像仪辐射定标方法 | |
Horvath et al. | Ground-based full-sky imaging polarimetry of rapidly changing skies and its use for polarimetric cloud detection | |
Cota et al. | PICASSO: an end-to-end image simulation tool for space and airborne imaging systems | |
CN107644416B (zh) | 一种基于地基云图的实时动态云量反演方法 | |
US10740609B1 (en) | System and method for space object detection in daytime sky images | |
CN107092752A (zh) | 一种基于光线追迹的光学相机仿真成像方法和*** | |
CN112258579A (zh) | 图像处理方法、装置、电子设备和计算机可读存储介质 | |
Jia et al. | High spatial resolution shortwave infrared imaging technology based on time delay and digital accumulation method | |
CN111340734B (zh) | 一种利用卷积神经网络模型的图像紫边校正方法 | |
KR102071879B1 (ko) | 하늘 영상을 이용한 미세먼지 측정 시스템 | |
CN111145351B (zh) | 考虑地物类型的Minnaert地形校正模型优化方法 | |
Ueda et al. | Underwater image synthesis from RGB-D images and its application to deep underwater image restoration | |
CN116519557A (zh) | 一种气溶胶光学厚度反演方法 | |
CN110689505B (zh) | 一种基于场景的星载遥感仪器自适应校正方法和*** | |
Kokka et al. | Flat-field calibration method for hyperspectral frame cameras | |
CN116452461A (zh) | 一种遥感数据处理方法、装置、电子设备及存储介质 | |
Zhou et al. | Photometric Calibration for the Beijing–Arizona Sky Survey and Mayall z-band Legacy Survey | |
CN117726550B (zh) | 一种多尺度门控注意力遥感图像去雾方法和*** | |
CN110070509A (zh) | 一种可见光遥感图像实时可视化***及可视化方法 | |
CN108225731B (zh) | 彩色面阵相机mtf测试方法及装置 | |
CN109141651A (zh) | 月基对地观测平台热红外传感器成像模拟方法 | |
Díaz et al. | Enhanced gap fraction extraction from hemispherical photography | |
CN110702228A (zh) | 一种航空高光谱影像的边缘辐射校正方法 | |
CA2784817C (en) | Filter setup learning for binary sensor | |
Yuan et al. | Tunable-liquid-crystal-filter-based low-light-level color night vision system and its image processing method |
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 |