CN112799144B - 一种基于梯度下降法的全天空成像仪的几何定标方法 - Google Patents

一种基于梯度下降法的全天空成像仪的几何定标方法 Download PDF

Info

Publication number
CN112799144B
CN112799144B CN202011614882.2A CN202011614882A CN112799144B CN 112799144 B CN112799144 B CN 112799144B CN 202011614882 A CN202011614882 A CN 202011614882A CN 112799144 B CN112799144 B CN 112799144B
Authority
CN
China
Prior art keywords
star point
real
zenith
sky
azimuth
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
Application number
CN202011614882.2A
Other languages
English (en)
Other versions
CN112799144A (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.)
POLAR RESEARCH INSTITUTE OF CHINA
Xidian University
Original Assignee
POLAR RESEARCH INSTITUTE OF CHINA
Xidian University
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 POLAR RESEARCH INSTITUTE OF CHINA, Xidian University filed Critical POLAR RESEARCH INSTITUTE OF CHINA
Priority to CN202011614882.2A priority Critical patent/CN112799144B/zh
Publication of CN112799144A publication Critical patent/CN112799144A/zh
Application granted granted Critical
Publication of CN112799144B publication Critical patent/CN112799144B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V8/00—Prospecting or detecting by optical means
    • G01V8/10—Detecting, e.g. by using light barriers
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V13/00—Manufacturing, calibrating, cleaning, or repairing instruments or devices covered by groups G01V1/00 – G01V11/00
    • Y—GENERAL 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
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation

Landscapes

  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • Engineering & Computer Science (AREA)
  • Manufacturing & Machinery (AREA)
  • Image Processing (AREA)

Abstract

本发明实施例公开了一种基于梯度下降法的全天空成像仪的几何定标方法。该方法包括：选取全天空图像中的目标星点，并确定目标星点在全天空图像坐标系下的星点位置坐标以及在极光观测站处全天空图像时刻目标星点的真实方位角和真实天顶角；基于真实方位角、全天空图像中的天顶位置和地理方位构建误差函数，并根据误差函数确定真实天顶位置和真实地理方位；根据星点位置坐标和真实天顶位置确定目标星点的成像半径，并根据真实天顶角和成像半径拟合得到全天空图像中成像半径与真实天顶角之间的关系。本发明实施例的技术方案，以实现不受地域局限快速准确对全天空成像仪进行几何定标，全天空成像仪几何定标更加灵活。

Description

一种基于梯度下降法的全天空成像仪的几何定标方法
技术领域
本发明实施例涉及几何定标技术领域，尤其涉及一种基于梯度下降法的全天空成像仪的几何定标方法。
背景技术
全天空成像仪是一种重要的地基空间环境光学观测设备，准确的“物像”光学几何对应关系是获取极光结构的实际尺度大小、确定全天空极光图像到地理、地磁坐标系投影变换以及实现全天空成像仪与卫星联合观测的基础，这对于研究相关的空间物理过程和机制是极其重要的。
传统的全天空极光成像仪几何定标方法需要在实验室条件下，需要借助专门的定标设备，通过旋转全天空相机对特定的标志物进行拍照，拟合得到成像半径与天顶角之间关系，上述几何定标方法不能确定全天空图像中天顶的像素点位置以及地理方位信息。此外，在全天空成像仪的实际使用过程中，由于成像固件(CCD等)、光路镜头更新、设备维修，全天空成像仪的“物像”关系均会发生变化，对于安装在野外观测站点(比如南北极地区)的全天空极光成像仪而言，采用实验室重新定标是一件非常困难的事情。另一方面，已有的基于星点的全天空极光成像仪几何定标方法中，需要手动解算天顶在全天空图像中的位置，过程复杂，利用单颗星点确定地理方位信息，存在偶然误差，拟合成像半径与天顶角的关系时所使用的星点数量较少，难以保证较高的精度。
发明内容
本发明实施例提供一种基于梯度下降法的全天空成像仪的几何定标方法，以实现不受地域局限快速准确对全天空成像仪进行几何定标，全天空成像仪几何定标更加灵活。
第一方面，本发明实施例提供了一种基于梯度下降法的全天空成像仪的几何定标方法，该几何定标方法包括：
选取全天空图像中的目标星点，并确定所述目标星点在全天空图像坐标系下的星点位置坐标以及在极光观测站处全天空图像时刻所述目标星点的真实方位角和真实天顶角；
基于所述真实方位角、所述全天空图像中的天顶位置和地理方位构建误差函数，并根据所述误差函数确定真实天顶位置和真实地理方位；
根据所述星点位置坐标和所述真实天顶位置确定所述目标星点的成像半径，并根据所述真实天顶角和所述成像半径拟合得到所述全天空图像中成像半径与真实天顶角之间的关系。
第二方面，本发明实施例还提供了一种基于梯度下降法的全天空成像仪的几何定标装置，该几何定标装置包括：
真实数据确定模块，用于选取全天空图像中的目标星点，并确定所述目标星点在全天空图像坐标系下的星点位置坐标以及在极光观测站处全天空图像时刻所述目标星点的真实方位角和真实天顶角；
误差函数构建模块，用于基于所述真实方位角、所述全天空图像中的天顶位置和地理方位构建误差函数；
几何定标确定模块，用于根据所述误差函数确定真实天顶位置和真实地理方位，并根据所述星点位置坐标和所述真实天顶位置确定所述目标星点的成像半径，并根据所述真实天顶角和所述成像半径拟合得到所述全天空图像中成像半径与真实天顶角之间的关系。
第三方面，本发明实施例还提供了一种全天空成像仪，该全天空成像仪包括：
一个或多个处理器；
存储装置，用于存储多个程序，
当所述多个程序中的至少一个被所述一个或多个处理器执行时，使得所述一个或多个处理器实现本发明第一方面实施例所提供的基于梯度下降法的全天空成像仪的几何定标方法。
第四方面，本发明实施例还提供了一种计算机可读存储介质，其上存储有计算机程序，该程序被处理器执行时实现本发明第一方面实施例所提供的基于梯度下降法的全天空成像仪的几何定标方法。
本发明实施例的技术方案，通过选取全天空图像中的目标星点，并确定所述目标星点在全天空图像坐标系下的星点位置坐标以及在极光观测站处全天空图像时刻所述目标星点的真实方位角和真实天顶角；基于所述真实方位角、所述全天空图像中的天顶位置和地理方位构建误差函数，并根据所述误差函数确定真实天顶位置和真实地理方位；根据所述星点位置坐标和所述真实天顶位置确定所述目标星点的成像半径，并根据所述真实天顶角和所述成像半径拟合得到所述全天空图像中成像半径与真实天顶角之间的关系。解决了现有技术中安装在野外观测站点的全天空成像仪采用实验室重新定标困难且过程复杂精度差的问题，以实现不受地域局限快速准确对全天空成像仪进行几何定标，全天空成像仪几何定标更加灵活。
附图说明
图1是本发明实施例一提供的一种基于梯度下降法的全天空成像仪的几何定标方法的流程图；
图2A是本发明实施例二提供的一种基于梯度下降法的全天空成像仪的几何定标方法的流程图；
图2B是本发明实施例提供的示例性的全天空图像以及黄河站处与全天空图像拍摄时间相对应的星空图的示意图；
图2C是本发明实施例提供的目标星点的星点误差项确定原理示意图；
图3是本发明实施例三提供的一种基于梯度下降法的全天空成像仪的几何定标装置的结构图；
图4是本发明实施例四提供的一种全天空成像仪的硬件结构示意图。
具体实施方式
为了使本发明的目的、技术方案和优点更加清楚，下面结合附图对本发明具体实施例作进一步的详细描述。可以理解的是，此处所描述的具体实施例仅仅用于解释本发明，而非对本发明的限定。
另外还需要说明的是，为了便于描述，附图中仅示出了与本发明相关的部分而非全部内容。在更加详细地讨论示例性实施例之前应当提到的是，一些示例性实施例被描述成作为流程图描绘的处理或方法。虽然流程图将各项操作(或步骤)描述成顺序的处理，但是其中的许多操作可以被并行地、并发地或者同时实施。此外，各项操作的顺序可以被重新安排。当其操作完成时所述处理可以被终止，但是还可以具有未包括在附图中的附加步骤。所述处理可以对应于方法、函数、规程、子例程、子程序等等。
实施例一
图1为本发明实施例一提供的一种基于梯度下降法的全天空成像仪的几何定标方法的流程图，本实施例可适用于对任意地点的全天空成像仪进行快速而准确的几何定标的情况，该方法可以由基于梯度下降法的全天空成像仪的几何定标装置来执行，该装置可以通过软件和/或硬件的形式实现。具体包括如下步骤：
S110、选取全天空图像中的目标星点，并确定所述目标星点在全天空图像坐标系下的星点位置坐标以及在极光观测站处全天空图像时刻所述目标星点的真实方位角和真实天顶角。
其中，通过全天空成像仪获取一组全天空图像，一组全天空图像包括一张、两张或多张全天空图像，具体选取的全天空图像的数量可以由本领域技术人员根据实际几何定标需求进行选择设置，本实施例对此不作任何限制。
进一步的，全天空图像中星点清晰，即选取多张星点清晰的全天空图像，以便于根据选取的全天空图像确定目标星点。
具体的，选取全天空图像中的目标星点，包括：选取多张全天空图像，并将标记出的每张全天空图像中的星点作为所述目标星点。
可以理解的是，每张全天空图像中可以标记多个目标星点，同时，每个目标星点可以被标记于多张全天空图像中，也可以只被标记于某一张全天空图像中，本实施例对目标星点的具体存在的全天空图像不作任何限制。
在上述实施例的基础上，在全天空图像中建立坐标系，进而确定目标星点在全天空图像坐标系下的星点坐标位置。进一步的，在选取全天空图像中目标星点后，利用星空软件Stellarium确定目标恒星在极光观测台站处与全天空图像拍摄时间对应的星空图，进而得到全天空图像中目标星点的真实方位角和真实天顶角，即Stellarium中的地点设置为极光观测台站处的地理位置，时间设置为全天空图像拍摄时间。可以理解的是，多个目标星点则可以得到全天空图像中各个目标星点的真实方位角和真实天顶角。
需要说明的是，方位角和天顶角的获取不限于利用星空软件Stellarium进行获取，也可以通过其他方位角和天顶角获取的计算机软件或设备进行获取，本实施例对方位角和天顶角的具体获取不作任何限制。
S120、基于所述真实方位角、所述全天空图像中的天顶位置和地理方位构建误差函数，并根据所述误差函数确定真实天顶位置和真实地理方位。
其中，误差函数用于基于假设全天空图像中的天顶位置和地理方位，推算各目标星点的方位角，并与各星点的真实方位角进行比较，进而确定真实天顶位置和真实地理方位。
在上述实施例的基础上，基于所述真实方位角、所述全天空图像中的天顶位置和地理方位构建误差函数，包括：基于所述真实方位角、所述全天空图像中的天顶位置和地理方位计算所述目标星点的星点误差项；根据所述星点误差项确定误差函数。
进一步的，基于所述真实方位角、所述全天空图像中的天顶位置和地理方位计算所述目标星点的星点误差项，包括：按照公式：
计算得到所述目标星点的星点误差项；其中，errori为第i个目标星点的星点误差项；∠SiZN中Si为第i个目标星点、Z为天顶、ZN为天顶Z指向正北方向的射线；Azi为第i个目标星点的所述真实方位角。根据所述星点误差项确定误差函数，包括：按照公式：/>构建误差函数；其中，error为误差函数；errori为第i个目标星点的星点误差项；m为目标星点个数。
在上述实施例的基础上，根据所述误差函数确定真实天顶位置和真实地理方位，包括：利用梯度下降法对天顶位置和地理方位进行更新，求解所述误差函数的最小值，并将所述最小值对应的天顶位置和地理方位确定为所述真实天顶位置和所述真实地理方位。
具体的，利用梯度下降法求解误差函数的最小值，进而得到全天空图像中的所述真实天顶位置和所述真实地理方位。
S130、根据所述星点位置坐标和所述真实天顶位置确定所述目标星点的成像半径，并根据所述真实天顶角和所述成像半径拟合得到所述全天空图像中成像半径与真实天顶角之间的关系。
具体的，根据所述星点位置坐标和所述真实天顶位置确定所述目标星点的成像半径，包括：按照公式：确定所述目标星点的成像半径；其中，(xi，yi)为第i个目标星点的星点位置坐标；/>为所述真实天顶位置。
进一步的，根据多个目标星点的所述真实天顶角和所述成像半径拟合得到所述全天空图像中成像半径与真实天顶角之间的关系。
本发明实施例的技术方案，通过选取全天空图像中的目标星点，并确定所述目标星点在全天空图像坐标系下的星点位置坐标以及在极光观测站处全天空图像时刻所述目标星点的真实方位角和真实天顶角；基于所述真实方位角、所述全天空图像中的天顶位置和地理方位构建所述误差函数，并根据所述误差函数确定真实天顶位置和真实地理方位；根据所述星点位置坐标和所述真实天顶位置确定所述目标星点的成像半径，并根据所述真实天顶角和所述成像半径拟合得到所述全天空图像中成像半径与真实天顶角之间的关系。解决了现有技术中安装在野外观测站点的全天空成像仪采用实验室重新定标困难且过程复杂精度差的问题，以实现不受地域局限快速准确对全天空成像仪进行几何定标，全天空成像仪几何定标更加灵活。
实施例二
图2A为本发明实施例二提供的一种基于梯度下降法的全天空成像仪的几何定标方法的流程图，本实施例以上述实施例为基础进行优化。
相应的，本实施例的方法具体包括：
S210、选取全天空图像中的目标星点，并确定所述目标星点在全天空图像坐标系下的星点位置坐标以及在极光观测站处全天空图像时刻所述目标星点的真实方位角和真实天顶角。
具体的，选择多张星点清晰的全天空图像，图2B是本发明实施例提供的示例性的全天空图像以及黄河站处与全天空图像拍摄时间对应的星空图的示意图，参见图2B中的左图，在全天空图像中建立坐标系xoy，坐标原点o在全天空图像的正中心处，x轴水平向右，y轴垂直向下，标记在全天空图像中的目标星点全天空图像坐标系下的星点位置坐标为(xi，yi)，其中i为星点序号。
进一步的，继续参见图2B中的右图，右图为利用星空软件Stellarium确定极光观测台站处全天空图像时刻对应的星空图，通过该星空图可以确定各目标星点的真实方位角Azi和真实天顶角Zai。
S220、基于所述真实方位角、所述全天空图像中的天顶位置和地理方位计算所述目标星点的星点误差项。
具体的，图2C是本发明实施例提供的目标星点的星点误差项确定原理示意图，参见图2C，圆o与全天空图像的边沿相切，圆心在坐标原点o处，半径为255.5像素，对应方程为：
x2+y2＝255.52
假设全天空图像中天顶的像素点位置为Z(xz，yz)，由天顶指向正北方向的射线ZN与圆o相交于点N，射线oN与x正向夹角为θ，则N点的坐标为(255.5cos(θ)，-255.5sin(θ))。
对于第i个目标星点，计算角∠SiZN的大小为：
基于所述真实方位角、所述全天空图像中的天顶位置和地理方位计算所述目标星点的星点误差项，包括：
按照公式：计算得到所述目标星点的星点误差项；
其中，errori为第i个目标星点的星点误差项；∠SiZN中Si为第i个目标星点、Z为天顶、ZN为天顶Z指向正北方向的射线；Azi为第i个目标星点的所述真实方位角。
S230、根据所述星点误差项确定误差函数。
在上述实施例的基础上，根据所述星点误差项确定误差函数，包括：按照公式：构建误差函数；其中，error为误差函数；errori为第i个目标星点的星点误差项；m为目标星点个数。
S240、利用梯度下降法求解所述误差函数的最小值，并将所述最小值对应的天顶位置和地理方位确定为所述真实天顶位置和所述真实地理方位。
具体的，误差函数一阶连续可微，求解误差函数关于参数xZ、yZ和θ的导数，分别和error′θ；
则沿梯度的反方向更新各参数，求解误差函数的最小值为：
θZ (k+1)＝θZ (k)-λ2·error′θ(k)
其中，k为当前迭代次数，λ1为xZ和yZ的更新步长，λ2为θ的更新步长。
进一步的，经过K次迭代后，误差函数关于和θ(K)的导数均为0，误差函数取得最小值。此时，全天空图像中天顶的坐标为/>N的坐标为(255.5cos(θ(K))，-255.5sin(θ(K)))，ZN与/>正方向的夹角β为：
由于全天空图像中的天顶位置和正北方位都具有唯一性，因此，误差函数取最小值，上述过程求得的真实天顶位置即为与天顶指向正北的方位ZN即为真实地理方位。
S250、根据所述星点位置坐标和所述真实天顶位置确定所述目标星点的成像半径，并根据所述真实天顶角和所述成像半径拟合得到所述全天空图像中成像半径与真实天顶角之间的关系。
具体的，计算第i个星点的成像半径ri，具体的公式如下：
其中，(xi，yi)为第i个目标星点在全天空图像坐标系中的星点位置坐标，为经过K次迭代后得到的真实天顶位置坐标。
进一步的，第i个星点的天顶角为Zai，则根据多颗星点可以拟合得到所述全天空图像中真实天顶角与成像半径之间的关系。
示例性的，本发明实施例对Kochab北极二、Mizar开阳、Dubhe天枢、Vega织女一、Capella五车二、Procyon南河三、Betelgeuse参宿四、Arcturus大角以及Jupiter木星进行了标记，共获得129颗星点，其中，部分星点对应的全天空图像获取时间以及全天空图像中星点的像素点坐标、方位角和高度角等信息如下表1所示。
表1部分星点对应的全天空图像获取时间以及全天空图像中星点的像素点坐标、方位角和高度角等信息
本发明实施例的目的是提供一种基于梯度下降法的全天空成像仪的几何定标方法，利用全天空图像中星点的位置信息，构建关于全天空图像中天顶位置以及地理方位的误差函数，并利用梯度下降法求解误差函数的最小值，实现对全天空成像仪的几何定标。利用本发明实施例提供的方法可以对任意地点的全天空成像仪进行快速、精确、自动的几何定标，并得到全天空图像中的天顶位置、地理方位以及成像半径与天顶角之间的关系。
实施例三
图3为本发明实施例三提供的一种基于梯度下降法的全天空成像仪的几何定标装置的结构图，本实施例可适用于对任意地点的全天空成像仪进行快速而准确的几何定标的情况。
如图3所示，所述几何定标装置包括：真实数据确定模块310、误差函数构建模块320和几何定标确定模块330，其中：
真实数据确定模块310，用于选取全天空图像中的目标星点，并确定所述目标星点在全天空图像坐标系下的星点位置坐标以及在极光观测站处全天空图像时刻所述目标星点的真实方位角和真实天顶角；
误差函数构建模块320，用于基于所述真实方位角、所述全天空图像中的天顶位置和地理方位构建误差函数；
几何定标确定模块330，用于根据所述误差函数确定真实天顶位置和真实地理方位，并根据所述星点位置坐标和所述真实天顶位置确定所述目标星点的成像半径，并根据所述真实天顶角和所述成像半径拟合得到所述全天空图像中成像半径与真实天顶角之间的关系。
本实施例的基于梯度下降法的全天空成像仪的几何定标装置，通过选取全天空图像中的目标星点，并确定所述目标星点在全天空图像坐标系下的星点位置坐标以及在极光观测站处全天空图像时刻所述目标星点的真实方位角和真实天顶角；基于所述真实方位角、所述全天空图像中的天顶位置和地理方位构建误差函数，并根据所述误差函数确定真实天顶位置和真实地理方位；根据所述星点位置坐标和所述真实天顶位置确定所述目标星点的成像半径，并根据所述真实天顶角和所述成像半径拟合得到所述全天空图像中成像半径与真实天顶角之间的关系。解决了现有技术中安装在野外观测站点的全天空成像仪采用实验室重新定标困难且过程复杂精度差的问题，以实现不受地域局限快速准确对全天空成像仪进行几何定标，全天空成像仪几何定标更加灵活。
在上述各实施例的基础上，选取全天空图像中的目标星点，包括：
选取多张全天空图像，并将标记出的每张全天空图像中的星点作为所述目标星点。
在上述各实施例的基础上，基于所述真实方位角、所述全天空图像中的天顶位置和地理方位构建误差函数，包括：
基于所述真实方位角、所述全天空图像中的天顶位置和地理方位计算所述目标星点的星点误差项；
根据所述星点误差项确定所述误差函数。
在上述各实施例的基础上，基于所述真实方位角、所述全天空图像中的天顶位置和地理方位计算所述目标星点的星点误差项，包括：
按照公式：计算得到所述目标星点的星点误差项；
其中，errori为第i个目标星点的星点误差项；∠SiZN中Si为第i个目标星点、Z为天顶、ZN为天顶Z指向正北方向的射线；Azi为第i个目标星点的所述真实方位角。
在上述各实施例的基础上，根据所述星点误差项确定所述误差函数，包括：
按照公式：构建所述误差函数；
其中，error为误差函数；errori为第i个目标星点的星点误差项；m为目标星点个数。
在上述各实施例的基础上，根据所述误差函数确定真实天顶位置和真实地理方位，包括：
利用梯度下降法求解所述误差函数的最小值，并将所述最小值对应的天顶位置和地理方位确定为所述真实天顶位置和所述真实地理方位。
在上述各实施例的基础上，根据所述星点位置坐标和所述真实天顶位置确定所述目标星点的成像半径，包括：
按照公式：确定所述目标星点的成像半径；
其中，(xi，yi)为第i个目标星点的星点位置坐标；为所述真实天顶位置。
上述各实施例所提供的基于梯度下降法的全天空成像仪的几何定标装置可执行本发明任意实施例所提供的基于梯度下降法的全天空成像仪的几何定标方法，具备执行基于梯度下降法的全天空成像仪的几何定标方法相应的功能模块和有益效果。
实施例四
图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 (6)

1.一种基于梯度下降法的全天空成像仪的几何定标方法，其特征在于，包括：
选取全天空图像中的目标星点，并确定所述目标星点在全天空图像坐标系下的星点位置坐标以及在极光观测站处全天空图像时刻所述目标星点的真实方位角和真实天顶角；
基于所述真实方位角、所述全天空图像中的天顶位置和地理方位构建误差函数，并根据所述误差函数确定真实天顶位置和真实地理方位；
基于所述真实方位角、所述全天空图像中的天顶位置和地理方位构建误差函数，包括：
基于所述真实方位角、所述全天空图像中的天顶位置和地理方位计算所述目标星点的星点误差项；
根据所述星点误差项确定所述误差函数；
基于所述真实方位角、所述全天空图像中的天顶位置和地理方位计算所述目标星点的星点误差项，包括：
按照公式：计算得到所述目标星点的星点误差项；
其中，errori为第i个目标星点的星点误差项；∠SiZN中Si为第i个目标星点、Z为天顶、ZN为天顶Z指向正北方向的射线；Azi为第i个目标星点的所述真实方位角；
根据所述星点误差项确定误差函数，包括：
按照公式：构建所述误差函数；
其中，error为误差函数；errori为第i个目标星点的星点误差项；m为目标星点个数；
根据所述误差函数确定真实天顶位置和真实地理方位，包括：
利用梯度下降法求解所述误差函数的最小值，并将所述最小值对应的天顶位置和地理方位确定为所述真实天顶位置和所述真实地理方位；
根据所述星点位置坐标和所述真实天顶位置确定所述目标星点的成像半径，并根据所述真实天顶角和所述成像半径拟合得到所述全天空图像中成像半径与真实天顶角之间的关系。
2.根据权利要求1所述的几何定标方法，其特征在于，选取全天空图像中的目标星点，包括：
选取多张全天空图像，并将标记出的每张全天空图像中的星点作为所述目标星点。
3.根据权利要求1所述的几何定标方法，其特征在于，根据所述星点位置坐标和所述真实天顶位置确定所述目标星点的成像半径，包括：
按照公式：确定所述目标星点的成像半径；
其中，(xi,yi)为第i个目标星点的星点位置坐标；为所述真实天顶位置；K为当前迭代次数。
4.一种基于梯度下降法的全天空成像仪的几何定标装置，其特征在于，包括：
真实数据确定模块，用于选取全天空图像中的目标星点，并确定所述目标星点在全天空图像坐标系下的星点位置坐标以及在极光观测站处全天空图像时刻所述目标星点的真实方位角和真实天顶角；
误差函数构建模块，用于基于所述真实方位角、所述全天空图像中的天顶位置和地理方位构建误差函数；
基于所述真实方位角、所述全天空图像中的天顶位置和地理方位构建误差函数，包括：
基于所述真实方位角、所述全天空图像中的天顶位置和地理方位计算所述目标星点的星点误差项；
根据所述星点误差项确定所述误差函数；
基于所述真实方位角、所述全天空图像中的天顶位置和地理方位计算所述目标星点的星点误差项，包括：
按照公式：计算得到所述目标星点的星点误差项；
其中，errori为第i个目标星点的星点误差项；∠SiZN中Si为第i个目标星点、Z为天顶、ZN为天顶Z指向正北方向的射线；Azi为第i个目标星点的所述真实方位角；
根据所述星点误差项确定误差函数，包括：
按照公式：构建所述误差函数；
其中，error为误差函数；errori为第i个目标星点的星点误差项；m为目标星点个数；
根据所述误差函数确定真实天顶位置和真实地理方位，包括：
利用梯度下降法求解所述误差函数的最小值，并将所述最小值对应的天顶位置和地理方位确定为所述真实天顶位置和所述真实地理方位；
几何定标确定模块，用于根据所述误差函数确定真实天顶位置和真实地理方位，并根据所述星点位置坐标和所述真实天顶位置确定所述目标星点的成像半径，并根据所述真实天顶角和所述成像半径拟合得到所述全天空图像中成像半径与真实天顶角之间的关系。
5.一种全天空成像仪，其特征在于，所述全天空成像仪包括：
一个或多个处理器；
存储装置，用于存储一个或多个程序；
当所述一个或多个程序被所述一个或多个处理器执行，使得所述一个或多个处理器实现如权利要求1-3中任一项所述的基于梯度下降法的全天空成像仪的几何定标方法。
6.一种计算机可读存储介质，其上存储有计算机程序，其特征在于，该程序被处理器执行时实现如权利要求1-3中任一项所述的基于梯度下降法的全天空成像仪的几何定标方法。
CN202011614882.2A 2020-12-30 2020-12-30 一种基于梯度下降法的全天空成像仪的几何定标方法 Active CN112799144B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011614882.2A CN112799144B (zh) 2020-12-30 2020-12-30 一种基于梯度下降法的全天空成像仪的几何定标方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011614882.2A CN112799144B (zh) 2020-12-30 2020-12-30 一种基于梯度下降法的全天空成像仪的几何定标方法

Publications (2)

Publication Number Publication Date
CN112799144A CN112799144A (zh) 2021-05-14
CN112799144B true CN112799144B (zh) 2024-01-16

Family

ID=75804634

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011614882.2A Active CN112799144B (zh) 2020-12-30 2020-12-30 一种基于梯度下降法的全天空成像仪的几何定标方法

Country Status (1)

Country Link
CN (1) CN112799144B (zh)

Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103544679A (zh) * 2012-07-09 2014-01-29 中国电力科学研究院 一种全天空成像仪地基云图畸变校正方法
CN103605887A (zh) * 2013-11-13 2014-02-26 北京航空航天大学 一种在轨卫星成像波浪水面太阳耀光时空分布仿真模型
CN203464930U (zh) * 2013-09-10 2014-03-05 国家电网公司 全天空成像仪
WO2018055619A1 (en) * 2016-09-25 2018-03-29 Israel Aerospace Industries Ltd. Celestial compass and method of calibrating
CN107917880A (zh) * 2017-11-06 2018-04-17 中国科学院寒区旱区环境与工程研究所 一种基于地基云图的云底高度反演方法
JP2019060754A (ja) * 2017-09-27 2019-04-18 国立研究開発法人情報通信研究機構 光学画像を用いた雲の高度および風速の計測方法
CN109738973A (zh) * 2019-02-14 2019-05-10 浙江中控太阳能技术有限公司 一种基于太阳方位的全天空成像仪标定方法
CN110764153A (zh) * 2019-11-15 2020-02-07 上海航天测控通信研究所 一种星载微波成像仪热镜背瓣在轨误差校正系统及方法
JP2020071854A (ja) * 2018-10-31 2020-05-07 大日本印刷株式会社 画像処理装置、画像処理方法及びプログラム
CN111536965A (zh) * 2020-05-25 2020-08-14 中国科学院长春光学精密机械与物理研究所 一种以全天空成像为背景的星图显示方法
CN112731569A (zh) * 2020-12-24 2021-04-30 中国极地研究中心 基于恒星辐射光谱及通量的全天空成像仪辐射定标方法

Patent Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103544679A (zh) * 2012-07-09 2014-01-29 中国电力科学研究院 一种全天空成像仪地基云图畸变校正方法
CN203464930U (zh) * 2013-09-10 2014-03-05 国家电网公司 全天空成像仪
CN103605887A (zh) * 2013-11-13 2014-02-26 北京航空航天大学 一种在轨卫星成像波浪水面太阳耀光时空分布仿真模型
WO2018055619A1 (en) * 2016-09-25 2018-03-29 Israel Aerospace Industries Ltd. Celestial compass and method of calibrating
JP2019060754A (ja) * 2017-09-27 2019-04-18 国立研究開発法人情報通信研究機構 光学画像を用いた雲の高度および風速の計測方法
CN107917880A (zh) * 2017-11-06 2018-04-17 中国科学院寒区旱区环境与工程研究所 一种基于地基云图的云底高度反演方法
JP2020071854A (ja) * 2018-10-31 2020-05-07 大日本印刷株式会社 画像処理装置、画像処理方法及びプログラム
CN109738973A (zh) * 2019-02-14 2019-05-10 浙江中控太阳能技术有限公司 一种基于太阳方位的全天空成像仪标定方法
CN110764153A (zh) * 2019-11-15 2020-02-07 上海航天测控通信研究所 一种星载微波成像仪热镜背瓣在轨误差校正系统及方法
CN111536965A (zh) * 2020-05-25 2020-08-14 中国科学院长春光学精密机械与物理研究所 一种以全天空成像为背景的星图显示方法
CN112731569A (zh) * 2020-12-24 2021-04-30 中国极地研究中心 基于恒星辐射光谱及通量的全天空成像仪辐射定标方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
Calibration and Validation of an All-Sky Imager;HUO Juan 等;ATMOSPHERIC AND OCEANIC SCIENCE LETTERS;第2卷(第4期);第220-223页 *
ORION software tool for the geometrical calibration of all-sky cameras;Juan Carlos Antuña-Sa´nchez 等;PLOS ONE;全文 *
利用全天空可见光图像反演天空辐亮度;霍娟 等;气象学报;第68卷(第06期);第800-807页 *
基于星点的全天空极光图像参数标定及其验证;张会朋 等;地球物理学报;第63卷(第02期);第401-411页 *

Also Published As

Publication number Publication date
CN112799144A (zh) 2021-05-14

Similar Documents

Publication Publication Date Title
CN107767422B (zh) 一种鱼眼镜头的校正方法、装置及便携式终端
CN108765328B (zh) 一种高精度多特征平面模板及其畸变优化和标定方法
CN108592950B (zh) 一种单目相机和惯性测量单元相对安装角标定方法
CN107192376B (zh) 基于帧间连续性的无人机多帧图像目标定位校正方法
CN106709944B (zh) 卫星遥感图像配准方法
CN106971408A (zh) 一种基于时空转换思想的摄像机标定方法
Zhou et al. Novel autonomous on-orbit calibration method for star sensors
CN114001756B (zh) 一种小视场星敏感器外场地面寻星方法
CN113137958B (zh) 一种rtk主机的放样控制方法、系统及存储介质
CN111598930B (zh) 彩色点云生成方法、装置及终端设备
CN112799144B (zh) 一种基于梯度下降法的全天空成像仪的几何定标方法
CN110411449B (zh) 一种航空侦察载荷目标定位方法、系统及终端设备
CN116309798A (zh) 一种无人机成像定位方法
CN110779517A (zh) 激光雷达的数据处理方法及装置、存储介质和计算机终端
CN115797468A (zh) 一种鱼眼摄像头安装高度的自动校正方法、装置以及设备
CN112085751B (zh) 一种基于云图阴影匹配算法的云层高度估算方法
CN109919998B (zh) 卫星姿态确定方法、装置和终端设备
CN114322996A (zh) 一种多传感器融合定位系统的位姿优化方法和装置
CN112053402A (zh) 利用偏振地理信息数据库获取航向角的方法
CN113657160B (zh) 一种船舶关联方法、船舶关联装置及电子设备
Kang et al. Positioning Errors of Objects Measured by Convolution Neural Network in Unmanned Aerial Vehicle Images
CN109269477A (zh) 一种视觉定位方法、装置、设备和存储介质
CN117406185B (zh) 雷达与相机间的外参标定方法、装置、设备及存储介质
CN111010494B (zh) 一种带有地理编码的光学卫星视频稳像方法及系统
CN117745818A (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