Hα全日面太阳图像云污染的去除

朱海波，杨云飞，邓辉，季凯帆

(昆明理工大学云南省计算机技术应用重点实验室，云南 昆明 650500)

Hα全日面太阳图像云污染的去除

朱海波，杨云飞，邓辉，季凯帆

(昆明理工大学云南省计算机技术应用重点实验室，云南 昆明 650500)

Hα全日面太阳图像对太阳物理研究有着重要的意义，但其观测过程中可能受天气的影响，导致观测图像含有云层覆盖的污染。主要通过对这种Hα全日面太阳图像的云污染进行处理，从而得到较为清晰的太阳活动细节。具体算法为:首先将含云污染的图像扣除模板图像的标准临边昏暗轮廓，然后对其进行中值滤波获得云层图像，最后在原始图像上扣除云层图像得到修复后的全日面图像。实验证明此方法可以有效地去除Hα全日面太阳图像上的云污染，并对太阳活动区域影响较小。

Hα全日面太阳图像；云污染的去除；平均径向轮廓曲线；中值滤波

CN53－1189/P ISSN1672－7673

Hα全日面太阳观测站点众多，观测历史悠久，数据资料极为丰富。例如全球日震网(Global Oscillation Network Group，GONG)[1]就有6个站点[2]:Big Bear太阳天文台、Mauna Loa天文台、Learmonth太阳天文台、Udaipur太阳天文台、Teide天文台和Cerro Tololo美洲洲际太阳天文台。Hα全日面太阳图像的研究，有助于了解太阳活动，对太阳物理研究有着重要的意义。然而在Hα全日面观测过程中，可能受到地球大气云层的影响，观测的Hα全日面图像中会有云污染，遮盖日面上的活动现象。本文主要针对云污染的Hα全日面太阳图像进行处理，消除其对太阳活动区域的影响，从而将被污染的区域及其太阳活动清晰地呈现出来，方便后续的太阳物理研究。

在云污染的去除上本文借鉴了卫星遥感图像中云污染去除[3]的方法，即将云层类比为一个不均匀的“云透镜”。当光线通过云层的时候，不同区域的透射率是不一样的。如果能得到这个透射率的变化，并假设其为线性的，就可以通过除法将其扣除，获得改正后的图像。通常Hα全日面图像上太阳的亮度呈现临边昏暗[4]现象。图1(a)是一张标准的Hα全日面图像(观测站点为Learmonth天文台，当地观测时间为2013年1月1日4点5分34秒)；图1(b)是图1(a)的平均径向轮廓曲线，表现了其临边昏暗的现象，其中横坐标表示Hα全日面太阳图像上的点到图像上太阳中心的距离，纵坐标是以日面中心归一化后所对应的强度。图1(c)是一张被云污染的Hα全日面图像(观测站点也为Learmonth天文台，当地观测时间为2013年1月1日2点45分34秒)；图1(d)是图1(c)的平均径向轮廓曲线，显然与临边昏暗曲线有较大的区别。本文所述的平均径向轮廓曲线通过在极坐标下，对半径尺度上所有点的强度值取中值得到的，这是获得临边昏暗曲线较为常用的一种做法，不但能反映全日面像的基本特征，同时还降低对太阳活动区域的影响。

一般来讲，同一架仪器获得的全日面太阳像的临边昏暗背景在一定时间内是相对稳定的，图1 (b)和(d)的区别明显是有云污染导致。通过比较被云污染的图像和标准临边昏暗背景，就可以获得这个“云透镜”的透射率的分布。

1 云污染Hα图像的修复处理方法

本文采用的方法主要包括如下几部分:首先选择一副没有云污染的图像并获得标准的临边昏暗轮廓；其次将被污染图像扣除这个标准轮廓[5－6]，得到一个比较“平坦”的太阳圆面；然后对其进行中值滤波消除太阳各种活动现象的影响，得到云层图像；最后用原始图像扣除云层图像得到修复后的全日面图像。

图1 (a)Hα全日面太阳图像；(b)(a)图的平均径向轮廓曲线；(c)含云污染的Hα全日面太阳图像；(d)(c)图的平均径向轮廓曲线Fig.1 (a)An Hα full-disk solar image.(b)The direction-averaged radial profile of the image in(a).(c)An Hα full-disk solar image affected by cloud extinction.(d)The direction-averaged radial profile of the image in(c)

首先，为判断一幅图像是否被云污染，或者为得到较为标准的临边昏暗图像背景，必须获得全日面像的径向轮廓曲线。

本文采用类似美国大熊湖天文台的径向轮廓曲线测量方式，基本步骤如下:

步骤1:求得图像中太阳圆面的圆心和半径。Hα全日面图像中天空背景和太阳圆面明显地存在两个灰度等级，基于此点，首先通过设定阙值，将图像进行二值化；然后对此图太阳圆面的边缘点进行圆拟合，得到圆心和半径；

步骤2:以太阳圆心为原点，将图像从笛卡尔坐标系转换到极坐标系；

步骤3:在极坐标系下，将一定角度范围内的径向轮廓进行中值处理，就得到在该角度范围内的平均径向轮廓值。图1(b)和(d)即为360°范围内的平均径向轮廓值，也就是整个太阳的平均径向轮廓。

如果将图像按90°分割成4个象限，就可以得到不同象限的平均径向轮廓曲线。在没有云污染的良好观测条件下，这4条曲线是彼此重合的。因此它们可以用于判断Hα图像是否存在云污染。具体如下:

(1)观察4个象限的平均径向轮廓曲线是否符合临边昏暗特性。如果4条曲线都很好地呈现临边昏暗性质，则图像每个象限区域内都不存在云污染；否则，比如呈现临边增亮或者径向轮廓变化极为不光滑，则可能存在云污染，导致图像质量降低；

(2)在4条平均径向轮廓曲线都具有临边昏暗性质的情况下，利用它们的残差值判断这4条曲线的重合度。如果它们基本重合，说明图像4个象限都不存在云污染，图像质量较好；否则，图像至少一个象限整体受云覆盖。

根据这两个准则，不但可以筛选出被云污染的Hα全日面像，同时也可以确定哪些图像拥有较好的质量，可以作为标准的临边昏暗背景。

在找到模板图像后，可以求得模板图像的360°平均径向轮廓曲线并进行高阶曲线拟合，以彻底消除背景图像上太阳活动的影响。对于在极坐标系下，用这个平均轮廓代替所有角度方向上的径向曲线，得到一个光滑的曲面，然后再将此曲面转换回笛卡尔坐标系，就可以得到一个标准的太阳临边昏暗背景图像。

用原始图像扣除模板图像的临边昏暗背景图像，得到“平坦”的太阳表面，但在这个表面上既有云层信息也有太阳的活动区域。通过对此图像进行中值滤波(采用“田”形滤波方法)，可以过滤掉上面的太阳活动信息，得到原始图像的云层图像。中值滤波窗口选择要大于太阳圆面的活动区域，这样可以在滤波过程中比较好地去除太阳圆面中的光斑和暗点，同时也不能过于太大，太大会使云层信息扩散到其它非云污染区域。本文采用太阳圆面上最大活动区域的两倍作为中值滤波窗口。最后将原始图和云层图像进行归一化，用原始图像扣除云层图像，得到修复后的Hα全日面太阳图像。

2 实验结果及分析

本文的实验数据来源于GONG网站。图2(a)为模板图像(观测站点为大熊湖太阳天文台，当地观测时间为2013年1月1日晚上23点32分54秒)，由图也可以看出这是一张不含云污染的Hα图像。图2(b)为模板图像的4个象限的平均径向轮廓曲线，4条曲线都呈现临边昏暗性质且基本重合，符合评价标准，这说明模板图像不含云信息，质量较好。图2(c)为模板图像的平均径向轮廓曲线，图2(d)为通过模板图像的平均径向轮廓曲线求得的模板临边昏暗背景图像。

图2 (a)模板Hα图像；(b)4条径向轮廓曲线；(c)标准临边昏暗轮廓；(d)模板临边昏暗图像Fig.2 (a)An Hα full-disk solar image used as a template.(b)The direction-averaged radial profiles in the four quadrants of the image in(a).(c)The limb-darkening profile of the template in(a).(d)An image showing the profile in(c)

图3 (a)原始图像；(b)扣除临边昏暗后的图像；(c)云层图像；(d)结果图像；(e)原始图像的三维展示；(f)原始图像的4条径向轮廓曲线；(g)结果图像的三维展示；(h)结果图像的4条径向轮廓曲线Fig.3 (a)A raw Hα full-disk solar image affected by cloud extinction.(b)The result image after the removal of the limb-darkening profile from the image in(a).(c)The image of cloud extinction obtained by performing median filtering on the image in(b).(d)The result image after the removal of the cloud image in(c)from the raw image in(a).(e)A three-dimensional plot of the raw image.(f)The direction-averaged radial profiles in the four quandrants of the raw image.(g)A three-dimensional plot of the result image.(h)The direction-averaged radial profiles in the four quandrants of the result image

处理过程以及结果如图3。图3显示了一个云污染图像的去除过程。图3(a)为实验中的原始图像(观测站点为大熊湖太阳天文台，当地观测时间为2013年1月1日晚上23点26分54秒)，显然是一张含有云污染的Hα全日面太阳图像。首先将原始图像图3(a)扣除模板图像的标准临边昏暗背景图像图2(d)，得到图像如图3(b)。对此图像作中值滤波得到云层图像，见图3(c)。原始图像扣除云层图像得到修复后的Hα图像，见图3(d)。

图3(e)为原始图像的三维展示，图像上“凹”进去的部分正是原始图像中被云遮盖的部分；图3(g)为复原后图像的三维展示，可以看出图3(e)中“凹”进去的部分也得到很好的修复，说明原始图像上的云污染得到很好的去除。图3(f)为原始图像的4条径向轮廓曲线，4条曲线变化不光滑且彼此不重合，不符合评价标准，说明图像上含有云污染，图像质量差；图3(h)为结果图像的4条径向轮廓曲线，4条曲线重合且呈现出临边昏暗现象，符合评价标准，说明原始图像中的云污染得到较好地去除。对比原始图像图3(a)和结果图像图3(d)，也可以看出原始图像中的云污染在修复后的Hα图像中已不复存在。综上所述，本文所述方法有效地去除了Hα全日面太阳图像中的云污染。

还可以通过比较处理前后图像的细节部分来观察Hα图像中太阳活动区域在处理前后的变化，图4为太阳圆面中活动频繁的区域在处理前后的变化，其中图4(a)和(c)为原始图像中太阳活动频繁的局部区域，图4(b)和(d)是与图4(a)和(c)相对应的区域在修复后的全日面图像的显示。

图4 原图与结果图像的局部对比Fig.4 Comparisons between the raw image and the result image in Fig.3 at some locations

通过观察图4可以发现，原始图像中的云污染得到有效去除，被云遮盖的部分清晰地呈现出来；太阳圆面上的光斑和暗点也完整地保留下来，活动区域受到的影响很小。可以看出，云污染的去除对太阳圆面的活动区域影响相对较少。

为了进一步验证本算法的通用性，对其它时间其它天文台的云污染较为严重的Hα全日面太阳图像进行了处理。图5(a)为大熊湖太阳天文台于当地时间2011年7月24日观测的图像，图5(d)为此图像经过本文方法处理的结果；图5(b)为Learmonth天文台于当地时间2013年1月1日观测的图像，图5(e)为处理结果；图5(c)为当地时间2012年12月31日Learmonth天文台观测的图像，图5(f)为处理结果。

图5 Hα图像云去除Fig.5 A set of Hα solar images affected by cloud extinction and corresponding result images after the cloud removal

图5原始图像中的云污染都得到了很好的去除，而且对太阳活动区域的影响较小，处理结果较为理想。这也验证了本文所述方法在Hα全日面太阳图像云污染的去除中具有较好的通用性。

3 总 结

本文基于平均径向轮廓曲线，采用中值滤波方法，去除Hα全日面太阳图像上的云污染。比较实验中处理前后图像的三维展示，处理前后图像4个象限的径向轮廓曲线以及原始图像和结果图像发现:Hα图像上的云污染有效地去除了，被云遮盖部分清晰地显示。对比处理前后图像的细节部分，可以看到被云遮盖部分的太阳活动能观察到了，太阳表面的谱斑、暗条等依然很好地存在，对太阳表面的活动区域影响较小。对不同观测器、不同时间以及云污染程度不同的图像进行了处理，得到的结果验证了该方法在去除Hα图像云污染上的有效性。

本方法也存在一些不足:首先处理过程中需要其他图像充当模板图像进行辅助处理；其次由于中值滤波是一种非线性平滑技术，这就导致修复后的Hα全日面图像整体辨识度降低，也会使结果图像中太阳活动的强度和范围缩小；第三，本文的云污染去除原理是基于“云透镜”的线性表现，这显然是一种近似。虽然文中的算法表现出较好的结果，但必须强调的是用这种恢复后的图像做高精度的测量工作是不合适的，因为被云污染的图像到底是一种降质的图像，而图像处理只能在有限的条件下尽可能地改善图像质量，但不能从根本上得到精确的结果。只有在实在无法控制观测条件的时候，图像质量恢复才是一种可以考虑的选择。

致谢：感谢GONG项目提供的实验图像。

参考文献:

[1] Leibacher J W.The Global Oscillation Network Group(GONG)Project[J].Advances in Space Research，1999，24(2):173－176.

[2] Frank H，George F，Suzanne F，et al.The Global Oscillation Network Group Site Survey II. Results[J].Solar Physics，1994，152(2):351－379.

[3] Li M，Liew S C，Kwoh L K.Automated production of cloud-free and cloud shaodow-free image mosaics from cloudy satellite imagery[C]//International Archives of Photogrammetry，International congress for photogrammetry and remote sensing，ISPRS XXth congress 20th，2004，25(B3): 754－758.

[4] Neckel H，Labs D.Solar limb darkening 1986－1990(λλ303 to 1099 nm)[J].Solar Physics，1994，153(1-2):91－114.

[5] Zharkova V V，Ipson S S，Zharkov S I，et al.A full-disk image standardisation of the synoptic solar observations at the Meudon Observatory[J].Solar Physics，2003，214(1):89－105.

[6] Denker C，Johannesson A，Marquette W，et al.Synoptic Hα full-disk observations of the sun from Big Bear Solar Observatory I.instrumentation，image Processing，data products，and first results[J].Solar Physics，1999，184(1):87－102.

Removal of Cloud Extinction for Hα Full-disk Solar Images

Zhu Haibo，Yang Yunfei，Deng Hui，Ji Kaifan
(Key Laboratory of Applications of Computer Technologies of the Yunnan Province，Kumming University of Science and Technology of Kumming，Kunming 650500，China，Email:jkf＠cnlab.net)

Hα full-disk solar images are important in research of solar science.Observations leading to such images may be affected by weather so that they show appreciable cloud extinction.This paper focuses on an approach to remove cloud extinction in Hα full-disk solar images.The approach can make details of solar activities easily recognizable in result images.We first evaluate and select a template image of Hα solar surface brightness.The template is constructed as an image modeling the Hα solar full disk with an isotropic surface brightness distribution following the direction-averaged Hα solar radial profile of a(cloud-extinction free)real image.We next subtract the template from a raw Hα full-disk solar image affected by cloud extinction，and perform median filtering on the difference image to obtain an image of cloud extinction.We finally subtract the raw image by the image of cloud extinction.The subtraction result is an Hα full-disk solar image with cloud extinction effectively removed.We have conducted some experiments，which validate our approach and show that it has little effect on image details of solar activities.Particularly，our experiments demonstrate that the performance of the approach is stable against patterns of cloud extinction.We also discuss shortcomings of our approach.

Hα full-disk solar image；Removal of cloud extinction；Direction-averaged radial profile；Median filtering

P182

A

1672－7673(2014)02－0134－06

2013－07－08；

2013－07－21

朱海波，男，硕士.研究方向:大规模数据处理.Email:498815638＠qq.com

季凯帆，男，研究员.研究方向:天文科学数据处理.Email:jkf＠cnlab.net