深圳市天文台全天云图监测系统的设计与实现
梅林1, 吴翔2, 周作超2     
1. 深圳市国家气候观象台(深圳市天文台), 广东 深圳 518040;
2. 深圳市一思佛科技有限公司, 广东 深圳 518040
摘要: 全天云图监测对天文观测至关重要,目前大多数天文台站安装了实时云图监测相机,但共同的难题是如何在光线变化较大的情况下进行正确曝光。为解决这一问题,深圳市天文台自行开发了一套拍摄控制系统,根据相机可用的拍摄参数,计算得到曝光值索引曲线,采用0.3作为曝光值误差,在亮度变化较大的各种环境下进行正确曝光拍摄全天云图,并采用最大类间方差自适应阈值算法,对云进行有效识别。经多次测试,系统在亮度变化较大较快的环境下能拍摄到正确曝光的云图,为天文观测计划提供了极大的便利。
关键词: 全天相机     自动拍摄     曝光参数     云识别     图像处理    
Design and Implementation of the All-Sky Cloud Monitoring System at Shenzhen Astronomical Observatory
Mei Lin1, Wu Xiang2, Zhou Zuochao2     
1. Shenzhen National Climate Observatory(Shenzhen Astronomical Observatory), Shenzhen 518040, China;
2. Shenzhen YSF Technology Co., Ltd., Shenzhen 518040, China
Abstract: All-sky cloud monitoring is of key importance for astronomical observation. Most astronomical observatories in the world have all-sky cloud cameras now. But one common concern is to get accurate exposure in a vast range of illumination conditions. To solve this problem, Shenzhen Astronomical Observatory has developed independently an automatic exposure control system. By using the parameters of the all-sky camera, this system calculates an index curve of exposure value, with error of 0.3, to adjust the exposure of each picture that the camera captures. And based on the pictures with accurate exposure, cloud can be effectively detected by the maximum interclass variance adaptive threshold selection method. It has been proved in practice that this system works well in varying illumination changes, obtaining accurate exposure even in low-light circumstances, and therefore providing substantial convenience to observation planning.
Key words: All-sky camera     Automatic exposure     Exposure parameters     Cloud detection     Image processing    
1 全天云图监测的现状与不足

全天云图观测面临一个关键问题,如何在复杂变化的环境以及光线变化较大的情况下进行正确曝光。目前相机自带的自动测光程序无法较好解决该问题。以佳能单反相机为例,它提供了一款控制软件EOS Utility,支持通过计算机对相机进行快门速度、光圈、感光度、白平衡等曝光参数的控制,可进行等间隔连续拍摄,并对拍摄的图片按顺序存储到指定目录。当相机设置为自动拍摄模式时,控制软件可根据相机自带的测光模式对曝光参数进行自动调整,一般光线充足的白天能获得正确曝光的图片,但是在日落、日出、天气突变等天空亮度变化较大的环境下,导致曝光过度或曝光不足,尤其在夜晚,控制程序更不能获得正确的曝光图像。图 1图 2为使用相机自带测光程序拍摄的曝光过度和曝光不足的云图序列。

http://gdlp01.c-wss.com/gds/7/0300026207/01/eu-3-6-0-w-im-zh.pdf

图 1 使用相机自带的自动测光程序,拍摄得到的2015年12月29日黎明日出阶段的云图序列,出现曝光过度现象 Figure 1 Overexposured images captured by the camera built-in program in the dawn of December 29, 2015
图 2 使用相机自带的自动测光程序,拍摄得到的2015年12月29日黄昏日落阶段的云图序列,出现曝光不足现象 Figure 2 Underexposured images captured by the camera built-in program in the evening of December 29, 2015

为弥补相机自动测光程序的不足,目前国内各天文台主要采用以下3种修正方法。

(1) 预设固定曝光参数法,即采用预先设定好的曝光参数,分别用于白天和夜晚。该方法可以获得白天和夜晚的正确曝光,但是在黎明日出以及黄昏日落的时间段不能获得正确的曝光。中国科学院国家天文台兴隆观测站的全天云图监测相机就采用该方法。

http://www.xinglong-naoc.org/weather/yuntu.jhtml

(2) 预测天光法,即根据地理经纬度、海拔、日期时间推算太阳月亮的位置,得到一个天光亮度值,设定一条曝光曲线。这个方法只对晴朗的天空有效,如果阴天有云或者有其他光害情况,则无法正确曝光。

(3) 外部测光法,即使用测光仪,根据环境亮度的变化,测量天空亮度值,自动调整相机的曝光参数,准确曝光。云南天文台丽江天文观测站自行开发了一款控制程序,采用实时天光亮度设置拍摄参数,解决了全天24小时正确曝光拍摄的问题[1]

2 深圳市天文台全天云图监测系统的设计与实现

深圳市天文台自行开发了一套拍摄控制系统,根据相机可用的拍摄参数,计算得到曝光值索引曲线,采用0.3作为曝光值误差,在亮度变化较大的环境下进行正确曝光拍摄全天云图,并采用最大类间方差自适应阈值算法,对云进行有效识别。

2.1 系统构成

深圳市天文台全天云图监测系统(图 3)由全天相机、拍摄控制系统、云识别系统3部分构成。其中,全天相机负责对全天云图进行实时拍摄,拍摄控制系统对全天相机的拍摄参数进行自动调整,云识别系统对拍摄得到的全天云图进行云判别,最终监测结果通过深圳市天文台网站、深圳市综合气象监测业务平台、深圳天文微信以及移动APP等服务渠道对外发布。图 3为深圳市天文台全天云图监测系统结构图。

http://www.szmb.gov.cn/yljc/

图 3 深圳市天文台全天云图监测系统图示 Figure 3 Diagram of the all-sky cloud monitoring system at Shenzhen astronomical observatory
2.2 全天相机的安装

深圳市天文台位于广东省深圳市大鹏新区南澳半岛,距离市区80 km,地理坐标为北纬22°28.94′,东经114°33.36′,海拔高度170 m,其东、南、西三面环海,周围无遮挡,全天相机的安装位置选择在天文台最高点户外水池上。拍摄相机采用佳能EOS 5D Mark II单反相机,镜头采用佳能EF 8~15 mm F4.0变焦鱼眼镜头,并将拍摄焦距固定在8 mm,可实现对角线180°的覆盖。为保证相机能够在全天候环境下工作,相机和一台用于控制的微型电脑被安装在一个防护箱内。防护箱的顶部装有一个直径为20 cm的半球玻璃罩。玻璃罩和侧面密封防水,底部留有通风口。箱体内安装了一个螺旋导流风扇,可对相机、微型电脑进行散热,保证箱体内外环境温度一致,防止箱体过热或玻璃罩结露。拍摄图像通过光纤传输到远端服务器。图 4为安装于室外的全天相机。

图 4 安装于室外的全天相机 Figure 4 The all-sky cloud monitoring camera stands outdoors
2.3 根据曝光值设置对应的拍摄参数

为了获得正确的曝光全天云图,在相机可用光圈、快门、感光度3个参数的基础上,根据曝光值计算公式[2],进行曝光值预计算,公式为

https://en.wikipedia.org/wiki/Exposure_value#Tabulated_exposure_values

EV=log2(F2/t)-log2(ISO/100),

其中,EV为曝光值;F为光圈值;t为快门速度;ISO为感光度。将相机所有可用参数进行组合,计算出曝光值并按大小排序,可得到如表 1的曝光参数以及图 5的曝光值索引曲线。同时,将曝光值索引曲线相邻取值相减,即可得出曝光值的误差索引图,如图 6,将曝光值误差允许范围取值0.3。

表 1 曝光值索引与曝光参数 Table 1 The index of exposure value and corresponding exposure parameters
曝光索引 光圈 快门 感光度 曝光值
0 2.8 1 100 2.97
1 4 1 200 3
2 4 1/2 200 4
3 5.6 1/2 200 4.7
图 5 曝光值索引曲线图 Figure 5 The index curve of exposure value
图 6 曝光值误差索引图 Figure 6 The index of the exposure value errors

拍摄时,先使用预定参数拍摄一张照片,并计算这张照片的曝光值,假设参考图片曝光值为EVref,则误差量为

ΔEV=log2(EV/EVref).

若误差值在允许范围内,则此照片直接作为拍摄结果进行保存;若误差超出范围,则需根据曝光值索引曲线查找正确的曝光参数进行重拍。

经多次测试,该方法在不同环境下以及亮度变化较大较快的阶段均能拍摄到正确曝光的云图。图 7图 8为夜晚以及白天得到的正确曝光的云图,图 9图 12为黎明日出、黄昏日落过渡阶段以及白天、夜晚得到的正确曝光云图序列。

图 7 根据曝光值索引曲线,拍摄得到的2017年6月11日20时16分的正确曝光云图 Figure 7 Image with accurate exposure adjusted by index curve of exposure value on 20:16, June 11, 2017
图 8 根据曝光值索引曲线,拍摄得到的2017年6月3日16时36分的正确曝光云图 Figure 8 Image with accurate exposure adjusted by index curve of exposure value on 16:36, June 3, 2017
图 9 根据曝光值索引曲线和曝光值误差,拍摄得到的2017年7月12日03时至09时黎明日出过渡阶段的正确曝光云图 Figure 9 Images with accurate exposure adjusted by index curve of exposure value with errors from 3:00 to 9:00, July 12, 2017
图 10 根据曝光值索引曲线和曝光值误差,拍摄得到的2017年7月12日16时至21时黄昏日落过渡阶段的正确曝光云图 Figure 10 Images with accurate exposure adjusted by index curve of exposure value with errors, from 16:00 to 21:00, July 12, 2017
图 11 根据曝光值索引曲线和曝光值误差,拍摄得到的2017年7月12日10时至15时之间阳光直射时间段的正确曝光云图 Figure 11 Images with accurate exposure adjusted by index curve of exposure value with errors from 10:00 to 15:00, July 12, 2017
图 12 根据曝光值索引曲线和曝光值误差,拍摄得到的2017年1月23日00时至06时之间正确曝光云图 Figure 12 Images with accurate exposure adjusted by index curve of exposure value with errors, from 0: 00 to 6: 00, January 23, 2017
2.4 全天云图中云识别

已有不少文献阐述对全天相机拍摄到的云图进行云自动识别的方法。如文[3]介绍了一种数字云量白天观测处理方法,选择合适的云判别的辐射比阈值,可对全天云量进行正确判断。文[4]提出了利用图像直方图理论直接计算和阈值分割两种计算云量的方法,均可对白天云图进行正确的云量判断。本文采用文[5]提出的最大类间方差自适应阈值算法对白天云进行识别。

由于全天云图为鱼眼镜头拍摄,在执行云判别前,需对云图进行预处理:先将原始云图裁剪成1 900 × 1 800像素的文件,再排除云图中的非天空部分,排除后剩下的像素被分成若干小方格,最后对每个小方格进行云判断,判断流程如图 13

图 13 最大类间方差自适应阈值算法识别云流程 Figure 13 Flowchart of cloud detection by maximum interclass variance adaptive threshold selection method

图 14为具体事例:

图 14 (a) 将原始云图(2 785 × 1 856)裁剪为处理图片(1 900 × 1 800);(b)排除非天空部分(粉红色区域);(c)通过最大类间方差法计算整个图像的自适应阈值T;(d)本例中计算出自适应阈值T=167,通过自适应阈值分割后,白色部分为云 Figure 14 (a) The processed image (1900 × 1800) reduced from the original cloud image (2785 × 1856); (b) The image substracted non-sky parts (pink areas); (c) Calculating the adaptive threshold T over the image with the maximum interclass variance method; (d) The adaptive threshold in this example is T=167. The white regions in the image are clouds after separating with the adaptive threshold

逐个像素判断是否为云。先判断蓝红波比值是否符合云标准,符合直接判断为云。如不符合,再通过像素点灰度值与自适应阈值T进行比较,小于阈值则为云。得出为云的像素点占比例72.93%,判断为8分云。

云判别完成后,系统获取图片信息中的拍摄时间、感光度、光圈值、曝光时间,并连同云量值以及文件名存入数据库,每10 min执行一次。

3 结语

深圳市天文台全天云图监测系统上线后运行稳定,各时间段以及各种天气环境下均能拍摄曝光正确的全天云图,为天文观测提供了有效保障,同时也为观测数据提供了质量控制的参考依据。在流星雨等重要天象事件期间,系统还可根据不同需求,设置拍摄间隔时间,以记录重要天象信息,图 15为监测到的火流星。此外,深圳市天文台还装有天空成像仪,可提供白天低分辨率云图、云量和云高信息,未来,可将二者进行参照对比,进行数据标校。

图 15 2016年12月15日晚监测到的火流星 Figure 15 A meteor captured by the all-sky cloud monitoring system on December 15, 2016
参考文献
[1] 彭焕文, 辛玉新, 和寿圣, 等. 丽江天文观测站全天相机介绍[J]. 天文研究与技术, 2015, 12(1): 89–95
Peng Huanwen, Xin Yuxin, He Shousheng, et al. An introduction to the all-sky camera at the YNAO Lijiang astronomical station[J]. Astronomical Research & Technology, 2015, 12(1): 89–95.
[2] Jacobson R E, Ray S F, Atteridge G G, et al. The manual of photography:photographic and digital imaging[M]. 9th ed. [S.l.]: Focal Pres, 2000: 318.
[3] 施洋, 姚永强, 刘立勇. 天文选址数字云量白天观测处理方法[J]. 天文研究与技术——国家天文台台刊, 2008, 5(4): 415–419
Shi Yang, Yao Yongqiang, Liu Liyong. A method for observing and counting day-time cloud amounts in an astronomical site survey[J]. Astronomical Research & Technology——Publications of National Astronomical Observatories of China, 2008, 5(4): 415–419.
[4] 杨健, 沈彦燕, 宋志刚. 基于图像处理技术的地基云图云量的识别[J]. 气象水文海洋仪器, 2009(3): 42–45
Yang Jian, Shen Yanyan, Song Zhigang. The recognition of ground nephogram cloudage based on image processing technology[J]. Meteorological, Hydrological and Marine Instruments, 2009(3): 42–45.
[5] 杨俊, 吕伟涛, 马颖, 等. 基于自适应阈值的地基云自动检测方法[J]. 应用气象学报, 2009, 20(6): 713–721
Yang Jun, Lv Weitao, Ma Ying, et al. An automatic groung-based cloud detection method based on adaptive threshold[J]. Journal of Applied Meteorological Science, 2009, 20(6): 713–721. DOI: 10.11898/1001-7313.20090609
由中国科学院国家天文台主办。
0

文章信息

梅林, 吴翔, 周作超
Mei Lin, Wu Xiang, Zhou Zuochao
深圳市天文台全天云图监测系统的设计与实现
Design and Implementation of the All-Sky Cloud Monitoring System at Shenzhen Astronomical Observatory
天文研究与技术, 2018, 15(2): 232-239.
Astronomical Research and Technology, 2018, 15(2): 232-239.
收稿日期: 2017-08-18
修订日期: 2017-09-04

工作空间