CN114663490B - 一种基于月球方位的全天空成像仪标定方法 - Google Patents

一种基于月球方位的全天空成像仪标定方法 Download PDF

Info

Publication number
CN114663490B
CN114663490B CN202210221753.XA CN202210221753A CN114663490B CN 114663490 B CN114663490 B CN 114663490B CN 202210221753 A CN202210221753 A CN 202210221753A CN 114663490 B CN114663490 B CN 114663490B
Authority
CN
China
Prior art keywords
moon
sky imager
lunar
calculating
shooting
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
CN202210221753.XA
Other languages
English (en)
Other versions
CN114663490A (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.)
North China Electric Power University
China Three Gorges Corp
Original Assignee
North China Electric Power University
China Three Gorges Corp
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 North China Electric Power University, China Three Gorges Corp filed Critical North China Electric Power University
Priority to CN202210221753.XA priority Critical patent/CN114663490B/zh
Publication of CN114663490A publication Critical patent/CN114663490A/zh
Application granted granted Critical
Publication of CN114663490B publication Critical patent/CN114663490B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00—Image analysis
    • G06T7/60—Analysis of geometric attributes
    • G06T7/66—Analysis of geometric attributes of image moments or centre of gravity
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01W—METEOROLOGY
    • G01W1/00—Meteorology
    • G01W1/18—Testing or calibrating meteorological apparatus
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00—Image analysis
    • G06T7/70—Determining position or orientation of objects or cameras
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00—Image analysis
    • G06T7/80—Analysis of captured images to determine intrinsic or extrinsic camera parameters, i.e. camera calibration

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • General Physics & Mathematics (AREA)
  • Environmental & Geological Engineering (AREA)
  • Geometry (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Atmospheric Sciences (AREA)
  • Biodiversity & Conservation Biology (AREA)
  • Ecology (AREA)
  • Environmental Sciences (AREA)
  • Image Processing (AREA)
  • Studio Devices (AREA)

Abstract

本发明适用于涉及气象学仪表的校准领域，尤其是涉及一种基于月球方位的全天空成像仪标定方法，包括以下步骤：S1：固定并安装全天空成像仪，并记录全天空成像仪所在安装位置经纬度；S2：设置全天空成像仪拍摄间隔自动采集获取一系列第一地基云图，并记录每张第一地基云图的拍摄日期和拍摄时间；S3：图像识别获取第一地基云图中月球的质心坐标；S4：根据拍摄日期、拍摄时间及安装位置经纬度，计算各拍摄时间下的月球高度角H和月球方位角L；S5：根据月球质心坐标、月球高度角和月球方位角，对全天空成像仪进行参数拟合，从而完成对全天空成像仪的标定。本发明具有标定简单快捷、标定精度高、标定成本低的技术特点。

Description

一种基于月球方位的全天空成像仪标定方法
技术领域
本发明涉及气象学仪表的校准领域，尤其是涉及一种基于月球方位的全天空成像仪标定方法。
背景技术
云团运动是造成短期光伏发电波动的主要原因，瞬时云变导致强烈的时空变化将显著影响馈入电网的太阳能电力稳定性。对未来指定时间段云团运动方向和速度的准确预测将大幅降低这一影响。地基云图观测是进行云量自动监测和预报的有效手段，该方法通过安装在地面的全天空成像仪获取局地天空的全部信息。
由于全天空成像仪采用视场角为180度的鱼眼镜头，其所摄图像具有较为严重的径向畸变。因此云团高精度预测的实现需建立在对全天空成像仪准确标定的基础上。
目前已存在多种全天空成像仪的标定方法，不同标定方法各有其优缺点。
如申请号为CN201410424349.8的中国专利文献公开了一种鱼眼镜头标定后鱼眼图像校正的方法，通过自制的标定板实现鱼眼摄像机模型参数的标定。该方法受棋盘格标定板加工工艺的影响精度较差，标定板与相机距离较近提高了误差影响的灵敏度，同时标定板视野覆盖率小需进行多次拍摄易造成误差积累，标定过程操作繁琐且标定结果准确率低。
申请号为CN201910114490.0的中国专利文献公开了一种基于太阳方位的全天空成像仪标定方法，通过确定图像太阳质心坐标与实际的太阳高度角、方位角完成对全天空成像仪的标定。该方法计算精度较低，首先因为太阳极亮，相机容易过曝光，为质心的识别带来困难；同时如果所在观测地区为北半球，一年中太阳运动的高度角跨度小，在冬季尤为显著，导致数据涵盖范围局限，拟合所得标定参数准确度低。
发明内容
为解决上述方法存在的不足，本发明公开了一种基于月球方位的全天空成像仪标定方法，能较好地解决由鱼眼镜头畸变而造成的标定范围局限以及标定过程复杂、计算量大等问题，也能解决标定参数准确度低的问题，能够显著提高了全天空成像仪的标定效率和精度。
一种基于月球方位的全天空成像仪标定方法，包括以下步骤：
S1：固定并安装所述全天空成像仪，记录所述全天空成像仪所在安装位置经纬度及时区；
S2：设置所述全天空成像仪拍摄间隔自动采集获取一系列第一地基云图，并记录每张所述第一地基云图的拍摄日期和拍摄时间；
S3：图像识别获取所述第一地基云图中月球的质心坐标；
S31：将所述全天空成像仪拍摄RGB格式的第一地基云图转换为HSV格式的第二地基云图；
S32：对所述第二地基云图进行空间分层，获取所述第二地基云图中V层图像；
S33：对所述V层图像进行阈值为T的二值化处理，获取V层图像灰度的第三地基云图；
S34：判断第一地基云图中的月球区域：
对所述第三地基云图进行圆形区域查找，设定所述圆形区域的半径为r个像素值，r∈(a，b)，查找到的圆形区域即为月球；
S35：计算所述圆形区域质心点坐标，即为所述月球质心坐标，并选取下一时间节点的第一地基云图进行所述S31至S35，直至所有第一地基云图完成图像识别。
S4：根据所述拍摄日期、拍摄时间、所述安装位置经纬度和时区，计算各拍摄时间下的月球高度角H和月球方位角L；
S5：根据所述月球质心坐标、所述月球高度角H和所述月球方位角L，对所述全天空成像仪进行参数拟合，从而完成对所述全天空成像仪的标定。
进一步地，所述步骤S1中，记录所述全天空成像仪所在安装位置经纬度及时区具体为：设定所述经纬度的经度为Lon，东经为正，西经为负；设定所述经纬度的纬度为Lat，北纬为正，南纬为负；设定所述时区为tz，东时区为正，西时区为负。
进一步地，所述步骤S4包括以下步骤：
S41：计算所述日期的儒略星历世纪T和儒略世纪JC：
所述拍摄时间为Y年M月D日hh时mm分ss秒，将世界时转换为儒略日JD：
JD＝[365.25×(Y+4716)]+[30.60001×(M+1)]+DT-1524.5；
其中，DT＝D+(hh-tz+(mm+ss)/60)/60/24，ΔT表示地球时与国际标准时间的差值；当M＝1或2时，Y＝Y-1，M＝M+12；
将所述儒略日JD转换为儒略星历日JDE:
JDE＝JD+ΔT/86400；
则，
T＝(JDE-2451545)/36525；
JC＝(JD-2451545)/36525；
S42：计算所述月球的开普勒元素，其中：
月球平黄经LM：
LM＝218.3164591+481267.88134236×T-0.0013268×T2+T3/538841-T4/65194000；
月球平近点角LW：
LW＝134.9634114+477198.8676313×T-0.0089970×T2+T3/69699-T4/14712000；
月球经度参数F：
F＝93.2720993+483202.0175273×T-0.0034029×T2+T3/352600-T4/863310000；
月球黄道倾角E：
E＝84381.448-4680.93×U-1.55×U2+1999.25×U3-51.38×U4-249.67×U5-39.05×U6+7.12×U7+27.87×U8+5.79×U9+2.45×U10；
其中，U＝JDE/10；
S43：计算所述月球与所述太阳的角距离D：
D＝297.8502042+445267.1115168×T-0.0016300×T2+T3/545868-T4/113065000；
S45：计算所述太阳的平近角点SW：
SW＝357.5291092+35999.0502909×T-0.0001536×T2+T3/24490000；
S46：计算考虑行星摄动月球黄经周期项∑I和黄纬周期项∑b：
∑I＝∑I0+3958×sinA1+1962×sin(LM-F)+318×sinA2；
∑b＝∑b0-2235×sinLM+382×sinA3+175×sin(A1-F)+175×sin(A1+F)+127×sin(LM-LW)-115×sin(LM+LW)；
其中，∑I0和∑b0为未考虑行星摄动下的月球黄经周期项和黄纬周期项；
A1、A2、A3分别用来处理金星摄动、木星摄动以及地球扁率摄动影响：
A1＝119.75+131.849×T；
A2＝53.09+479264.290×T；
A3＝313.45+481266.484×T；
S47：计算所述月球对应地心的高度角H0和方位角L0：
H0＝LM+∑I/1000000；
L0＝∑b/1000000；
S48：计算所述月球赤经A和格林威治平均恒星时V：
A＝arctan2×((sinH0×cosE-tanL0×sinE)/cosH0)；
V＝280.46061837+360.98564736629×(JD-2451545)+0.000387933×JC2-JC3/37810000；
S49：计算所述安装位置的月球高度角H和方位角L：
H＝arcsin(sinLat×sinLon+cosLat×cosLot×cosX)；
L＝arctan(2×sinX/(cosX×sinLat-tanLon×cosLat))；
其中，X＝V+Lon-A表示所述位置地方时角。
进一步地，所述步骤S5包括以下步骤：
S51:根据多张所述第一地基云图对应的所述月球质心坐标，计算每张所述第一地基云图月球质心点与所述第一地基云图图像中心点的像素距离R和两者连线的方向向量l；根据多张所述第一地基云图拍摄时刻，计算相应所述月球高度角H和所述月球方位角L；
S52:根据所述多个像素距离R和所述月球高度角H，拟合得到所述全天空成像仪高度角的多项式参数p1、p2、p3：
H＝p1×R2+p2×R+p3；
S53:记录所述多张第一地基云图中方位角L为0度和90度的所述方向向量为l1和l2，所述全天空成像仪内任意方向与l1和l2的夹角表示为q1和q2，则所述全天空成像仪任意方向方位角L0：
进一步地，步骤S34中，设置圆形区域查找的查找灵敏度为S。
进一步地，步骤S34中，若对所述圆形区域无查找结果则将灵敏度S增加一个梯度设定值进行重复搜索，直至查找到相应圆形区域。
进一步地，步骤S33中，设定所述二值化处理的阈值T初始设置为0.99。
采用本发明的一种基于月球方位的全天空成像仪标定方法，相对于现有技术，至少具有以下有益效果：
(1)本发明无需标定板等外部设备，标定过程便捷简单，避免了由于引入棋盘格所致的标定视野范围局限及其加工工艺差异对全天空成像仪标定结果误差的影响，达到了标定精准、操作简单、成本低廉的技术效果；
(2)本发明不涉及复杂的图像变换等数学过程，计算结果迅速，同时避免多步图像投影变换造成的误差积累效应，达到了高效快捷的技术效果；
(3)本发明采用月球作为标定参照物，相比于以太阳作为参照物的地理标定方法，月球亮度适中不会引起相机过曝光，同时观测时段为夜晚，月球与周围天空亮度差异明显，避免了其他光源对识别的影响，使计算结果更为精准；北半球月球高度角全年多数时间较太阳具有更广运动跨度范围，标定数据覆盖更为全面，拟合结果更为准确，达到了大幅降低识别误差、提高识别精度的技术效果。
(4)采用本发明的月球识别方法，识别方法简单，能够快速有效地识别月球，进一步提高了标定的效率。
附图说明
为了更清楚地说明本发明实施例的技术方案，下面将对本发明实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面所描述的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1是本发明实施例的一种基于月球方位的全天空成像仪标定方法的流程示意图。
具体实施方式
以下的说明提供了许多不同的实施例、或是例子，用来实施本发明的不同特征。以下特定例子所描述的元件和排列方式，仅用来精简的表达本发明，其仅作为例子，而并非用以限制本发明。
一种基于月球方位的全天空成像仪标定方法，如图1所示，包括以下步骤：
S1：固定并安装所述全天空成像仪，记录所述全天空成像仪所在安装位置经纬度及时区；
设定所述经纬度的经度为Lon，东经为正，西经为负；设定所述经纬度的纬度为Lat，北纬为正，南纬为负；设定所述时区为tz，东时区为正，西时区为负。
S2：设置所述全天空成像仪拍摄间隔自动采集获取一系列第一地基云图，并记录每张所述第一地基云图的拍摄日期和拍摄时间；
本发明的全天空成像仪在多日夜晚时段采集获取各第一地基云图的拍摄时间为Y年M月D日hh时mm分ss秒；
S3：图像识别获取所述第一地基云图中月球的质心坐标；
S31：将所述全天空成像仪拍摄RGB格式的第一地基云图转换为HSV格式的第二地基云图；其中，R为红、G为绿、B为蓝色通道，H为色调、S为饱和度、V为明度；
S32：对所述第二地基云图进行空间分层，获取所述第二地基云图中V层图像，也就是获取第二地基云图的明度图像；
S33：对所述V层图像进行阈值为T的二值化处理，获取V层图像灰度的第三地基云图；
作为优选，设定所述二值化处理的阈值T初始设置为0.99，即当V层图像的明度大于0.99即设置为1，小于0.99设置为0；
S34：判断第一地基云图中的月球区域：
对所述第三地基云图进行圆形区域查找，设定所述圆形区域的半径为r个像素，r∈(a，b)，查找到的圆形区域即为月球；例如，设定r∈(5，10)，即设定了查找的圆形区域的半径为5-10个像素。
作为优选，设置圆形区域查找的灵敏度为S，本发明中的灵敏度是指灵敏度指圆形区域查找结果与圆形的匹配程度，灵敏度越高，查找区域与圆形吻合度越高。
若对所述圆形区域无查找结果则将灵敏度S增加一个梯度设定值进行重复搜索，直至查找到相应圆形区域，例如梯度设定值可以是0.02等。
S35：计算所述圆形区域质心点坐标，即为所述月球质心坐标，并选取下一时间节点的第一地基云图进行所述S31至S35，直至所有第一地基云图完成图像识别。
该步骤中，查找出的圆形区域的圆心即为坐标，相应地，圆心坐标即为质心点坐标。
S4：根据所述拍摄日期、拍摄时间、所述安装位置经纬度和时区，计算各拍摄时间下的月球高度角H和月球方位角L，包括以下步骤：
S41：计算所述日期的儒略星历世纪T和儒略世纪JC：
所述拍摄时间为Y年M月D日hh时mm分ss秒，将世界时转换为儒略日JD：
JD＝[365.25×(Y+4716)]+[30.60001×(M+1)]+DT-1524.5；
其中，DT＝D+(hh-tz+(mm+ss)/60)/60/24，DT为十进制时间日期，ΔT表示地球时与国际标准时间的差值，由观测得到，例如全天空成像仪位于东八区时区，则ΔT＝70秒；
上式中，当M＝1或2时，Y＝Y-1，M＝M+12；
将所述儒略日JD转换为儒略星历日JDE：
JDE＝JD+ΔT/86400；
则，
T＝(JDE-2451545)/36525；
JC＝(JD-2451545)/36525；
S42：计算所述月球的开普勒元素，其中：
月球平黄经LM：
LM＝218.3164591+481267.88134236×T-0.0013268×T2+T3/538841-T4/65194000；
月球平近点角LW：
LW＝134.9634114+477198.8676313×T-0.0089970×T2+T3/69699-T4/14712000；
月球经度参数F：
F＝93.2720993+483202.0175273×T-0.0034029×T2+T3/352600-T4/863310000；
月球黄道倾角E：
E＝84381.448-4680.93×U-1.55×U2+1999.25×U3-51.38×U4-249.67×U5-39.05×U6+7.12×U7+27.87×U8+5.79×U9+2.45×U10；
其中，U＝JDE/10；
S43：计算所述月球与所述太阳的角距离D：
D＝297.8502042+445267.1115168×T-0.0016300×T2+T3/545868-T4/113065000；
S45：计算所述太阳的平近角点SW：
SW＝357.5291092+35999.0502909×T-0.0001536×T2+T3/24490000；
S46：计算考虑行星摄动月球黄经周期项∑I和黄纬周期项∑b：
∑I＝∑I0+3958×sinA1+1962×sin(LM-F)+318×sinA2；
∑b＝∑b0-2235×sinLM+382×sinA3+175×sin(A1-F)+175×sin(A1+F)+127×sin(LM-LW)-115×sin(LM+LW)；
其中，∑I0和∑b0为未考虑行星摄动下的月球黄经周期项和黄纬周期项，由月球经度和纬度周期组合系数所得；A1、A2、A3分别用来处理金星摄动、木星摄动以及地球扁率摄动影响：
A1＝119.75+131.849×T；
A2＝53.09+479264.290×T；
A3＝313.45+481266.484×T；
S47：计算所述月球对应地心的高度角H0和方位角L0：
H0＝LM+∑I/1000000；
L0＝∑b/1000000；
S48：计算所述月球赤经A和格林威治平均恒星时V：
A＝arctan2×((sinH0×cosE-tanL0×sinE)/cosH0)；
V＝280.46061837+360.98564736629×(JD-2451545)+0.000387933×JC2-JC3/37810000；
S49：计算所述安装位置的月球高度角H和方位角L：
H＝arcsin(sinLat×sinLon+cosLat×cosLot×cosX)；
L＝arctan(2×sinX/(cosX×sinLat-tanLon×cosLat))；
其中，X＝V+Lon-A表示所述位置地方时角。
S5：根据所述月球质心坐标、所述月球高度角H和所述月球方位角L，对所述全天空成像仪进行参数拟合，从而完成对所述全天空成像仪的标定。
S51:根据多张所述第一地基云图对应的所述月球质心坐标，计算每张所述第一地基云图月球质心点与所述第一地基云图图像中心点的像素距离R和两者连线的方向向量l；根据多张所述第一地基云图拍摄时刻，计算相应所述月球高度角H和所述月球方位角L；
S52:根据所述多个像素距离R和所述月球高度角H，拟合得到所述全天空成像仪高度角的多项式参数p1、p2、p3：
H＝p1×R2+p2×R+p3；
S53:记录所述多张第一地基云图中方位角L为0度和90度的所述方向向量为l1和l2，所述全天空成像仪内任意方向与l1和l2的夹角表示为q1和q2，则所述全天空成像仪任意方向方位角L0：
其中，p1、p2、p3、q1、q2为全天空成像仪标定模型所需参数，通过多张第一地基云图及其对应数据进行参数拟合，可得到上述各标定参数，由此完成对全天空成像仪的标定。
以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明的保护范围之内。

Claims (5)

1.一种基于月球方位的全天空成像仪标定方法，其特征在于，包括以下步骤：
S1：固定并安装所述全天空成像仪，记录所述全天空成像仪所在安装位置经纬度及时区；
S2：设置所述全天空成像仪拍摄间隔自动采集获取一系列第一地基云图，并记录每张所述第一地基云图的拍摄日期和拍摄时间；
S3：图像识别获取所述第一地基云图中月球的质心坐标；
S31：将所述全天空成像仪拍摄RGB格式的第一地基云图转换为HSV格式的第二地基云图;
S32：对所述第二地基云图进行空间分层，获取所述第二地基云图中V层图像；
S33：对所述V层图像进行阈值为T的二值化处理，获取V层图像灰度的第三地基云图；
S34：判断第一地基云图中的月球区域：
对所述第三地基云图进行圆形区域查找，设定所述圆形区域的半径为r个像素值，，查找到的圆形区域即为月球；
S35：计算所述圆形区域质心点坐标，即为所述月球质心坐标，并选取下一时间节点的第一地基云图进行所述S31至S35，直至所有第一地基云图完成图像识别；
S4：根据所述拍摄日期、拍摄时间、所述安装位置经纬度和时区，计算各拍摄时间下的月球高度角H和月球方位角L；
S5：根据所述月球质心坐标、所述月球高度角H和所述月球方位角L，对所述全天空成像仪进行参数拟合，从而完成对所述全天空成像仪的标定；
所述步骤S4包括以下步骤：
S41：计算所述日期的儒略星历世纪T和儒略世纪JC：
所述拍摄时间为Y年M月D日hh时mm分ss秒，将世界时转换为儒略日JD：
；
其中，，/>表示地球时与国际标准时间的差值；当M=1或2时，Y=Y-1，M=M+12；
将所述儒略日JD转换为儒略星历日JDE:
；
则，
；
；
S42：计算所述月球的开普勒元素，其中：
月球平黄经LM：
；
月球平近点角LW：
；
月球经度参数F：
；
月球黄道倾角E：
其中，U=JDE/10；
S43：计算所述月球与所述太阳的角距离D：
；
S45：计算所述太阳的平近角点SW：
；
S46：计算考虑行星摄动月球黄经周期项和黄纬周期项/>：
其中，和/>为未考虑行星摄动下的月球黄经周期项和黄纬周期项；A1、A2、A3分别用来处理金星摄动、木星摄动以及地球扁率摄动影响：
S47：计算所述月球对应地心的高度角H0和方位角L0：
S48：计算所述月球赤经A和格林威治平均恒星时V：
S49：计算所述安装位置的月球高度角H和方位角L：
其中，表示所述位置地方时角；
所述步骤S5包括以下步骤：
S51:根据多张所述第一地基云图对应的所述月球质心坐标，计算每张所述第一地基云图月球质心点与所述第一地基云图图像中心点的像素距离R和两者连线的方向向量l；根据多张所述第一地基云图拍摄时刻，计算相应所述月球高度角H和所述月球方位角L；
S52:根据所述多个像素距离R和所述月球高度角H，拟合得到所述全天空成像仪高度角的多项式参数p1、p2、p3：
；
S53:记录所述多张第一地基云图中方位角L为0度和90度的所述方向向量为l1和l2，所述全天空成像仪内任意方向与l1和l2的夹角表示为q1和q2，则所述全天空成像仪任意方向方位角L0：
。
2.根据权利要求1所述的一种基于月球方位的全天空成像仪标定方法，其特征在于，所述步骤S1中，记录所述全天空成像仪所在安装位置经纬度及时区具体为：设定所述经纬度的经度为Lon，东经为正，西经为负；设定所述经纬度的纬度为Lat，北纬为正，南纬为负；设定所述时区为tz，东时区为正，西时区为负。
3.根据权利要求1所述的一种基于月球方位的全天空成像仪标定方法，其特征在于，步骤S34中，设置圆形区域查找的查找灵敏度为S。
4.根据权利要求1所述的一种基于月球方位的全天空成像仪标定方法，其特征在于，步骤S34中，若对所述圆形区域无查找结果则将灵敏度S增加一个梯度设定值进行重复搜索，直至查找到相应圆形区域。
5.根据权利要求1所述的一种基于月球方位的全天空成像仪标定方法，其特征在于，步骤S33中，设定所述二值化处理的阈值T初始设置为0.99。
CN202210221753.XA 2022-03-07 2022-03-07 一种基于月球方位的全天空成像仪标定方法 Active CN114663490B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210221753.XA CN114663490B (zh) 2022-03-07 2022-03-07 一种基于月球方位的全天空成像仪标定方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210221753.XA CN114663490B (zh) 2022-03-07 2022-03-07 一种基于月球方位的全天空成像仪标定方法

Publications (2)

Publication Number Publication Date
CN114663490A CN114663490A (zh) 2022-06-24
CN114663490B true CN114663490B (zh) 2024-04-05

Family

ID=82028840

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210221753.XA Active CN114663490B (zh) 2022-03-07 2022-03-07 一种基于月球方位的全天空成像仪标定方法

Country Status (1)

Country Link
CN (1) CN114663490B (zh)

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103968859A (zh) * 2014-04-30 2014-08-06 中国科学院长春光学精密机械与物理研究所 一种超大视场紫外临边成像仪的几何定标方法
KR20170125716A (ko) * 2016-05-04 2017-11-15 임재형 대상물의 위치 정보 측정 장치 그리고 그 제어 방법
CN109579874A (zh) * 2018-12-14 2019-04-05 天津津航技术物理研究所 一种光电平台现场标北方法
CN109738973A (zh) * 2019-02-14 2019-05-10 浙江中控太阳能技术有限公司 一种基于太阳方位的全天空成像仪标定方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103968859A (zh) * 2014-04-30 2014-08-06 中国科学院长春光学精密机械与物理研究所 一种超大视场紫外临边成像仪的几何定标方法
KR20170125716A (ko) * 2016-05-04 2017-11-15 임재형 대상물의 위치 정보 측정 장치 그리고 그 제어 방법
CN109579874A (zh) * 2018-12-14 2019-04-05 天津津航技术物理研究所 一种光电平台现场标北方法
CN109738973A (zh) * 2019-02-14 2019-05-10 浙江中控太阳能技术有限公司 一种基于太阳方位的全天空成像仪标定方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
基于太阳形状的直接辐射精细测量;牛一森;中国优秀硕士论文电子期刊网;20220115;全文 *

Also Published As

Publication number Publication date
CN114663490A (zh) 2022-06-24

Similar Documents

Publication Publication Date Title
CN110244766B (zh) 一种光伏电站无人机巡检航线的规划方法及系统
Rodríguez-Benítez et al. Assessment of new solar radiation nowcasting methods based on sky-camera and satellite imagery
KR101890673B1 (ko) 일사량 추정을 위해 전천사진으로부터 운량을 계산하는 방법, 상기 계산한 운량을 이용한 태양광 발전량 예측 장치
AU2017261849A1 (en) "Solar power forecasting"
CN107644416B (zh) 一种基于地基云图的实时动态云量反演方法
CN104913780B (zh) 集成gnss和ccd天顶筒的高精度垂线偏差快速测量方法
CN106600589B (zh) 一种基于海岸线区域遥感图的星载光谱仪指向配准方法
CN107710098A (zh) 基于制图模型来控制太阳能追踪器的定向的方法
CN108509844B (zh) 基于地基云图的vper云团识别方法
JP2009282406A (ja) 定点観測装置、及び定点観測方法
CN116182805A (zh) 基于遥感图像的国土测绘方法
CN109738973B (zh) 一种基于太阳方位的全天空成像仪标定方法
Spurný et al. Atmospheric trajectory and heliocentric orbit of the Ejby meteorite fall in Denmark on February 6, 2016
CN110689505B (zh) 一种基于场景的星载遥感仪器自适应校正方法和系统
CN102752504A (zh) 一种宽视场线阵ccd相机的相对辐射校正方法
CN113936031A (zh) 一种基于机器视觉的云影轨迹预测方法
Li et al. Hybrid cloud detection algorithm based on intelligent scene recognition
CN114663490B (zh) 一种基于月球方位的全天空成像仪标定方法
Savoy et al. Geo-referencing and stereo calibration of ground-based whole sky imagers using the sun trajectory
CN112085260B (zh) 一种基于云图预测的光伏功率超短期预测方法
Liu et al. High-spatial-resolution nighttime light dataset acquisition based on volunteered passenger aircraft remote sensing
CN113674391A (zh) 一种升降轨成像双星地物表观辐亮度差异计算方法
Barnie et al. Volcanic plume height monitoring using calibrated web cameras at the Icelandic Meteorological Office: system overview and first application during the 2021 Fagradalsfjall eruption
CN109727217B (zh) 基于改进Criminisi算法的地基云图修复方法
JP2791926B2 (ja) 受光量測定方法、日照条件測定方法、太陽電池の設置方法および太陽電池システム

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