[发明专利]基于行星际和地磁参数的极光卵边界位置预测方法

说明书

技术领域

本发明属于空间物理领域，涉及极光物理，具体涉及一种基于行星际和地磁参数的极光卵边界位置预测方法。

背景技术

极光卵是极光粒子沉降以磁极为中心在地球南北极地区形成的椭圆带状区域，是地球上来自太阳的能量粒子的影响区域。极光卵能从全局尺度上对空间天气进行表征和预测，所以一个好的极光卵模型能有效地帮助理解空间天气并预测其对极区和次极区的影响。

现有的对极光卵大小和位置的研究工作分为两大类：一种是定性分析，如讨论行星际磁场(IMF)的By、Bz分量对极光卵边界位置的影响；另一种是单变量回归分析，回归量取某个地磁指数，如极光电集流指数(AE)，地磁扰动指数(KP)等。

可事实上，一方面，极光活动是受所有行星际参数共同影响的，所以传统的一元回归分析或定量分析的方法不能充分表征极光卵的变化。另一方面，从物理源头来说，这些地磁活动指数(AE，Kp等)与极光是同一个层面的东西：它们都受太阳风-磁层-电离层耦合作用的影响。换句话说，这些地磁活动指数仅仅能从一个侧面来反映极光卵的属性，因为它们都受同一个“源”的影响，但并非因果关系。此外，以往研究还有一个明显缺点是将地磁指数按大小划分为几个区间，每个区间里只取其均值来做回归分析。这样就会出现一些非常相近的值被划分到了两个区间里，导致回归结果完全不同；此外，仅仅一个均值也不足以表征整个区间。

卫星图像是同时确定整个(或绝大部分)极光卵边界的唯一途径。可是从海量的图像中人工确定极光卵边界是非常繁琐的，几乎不太可能完成。

发明内容

本发明的目的在于提供一种基于行星际和地磁参数的UVI图像极光卵边界位置建模和预测方法，以克服上述现有技术存在的缺陷，本发明可用于对极光卵边界(包括高纬极向边界和低纬赤道向边界)位置预测，可以有效地帮助理解空间天气并预测其对极区和次极区的影响。

为达到上述目的，本发明采用如下技术方案：

基于行星际和地磁参数的极光卵边界位置预测方法，包括以下步骤：

步骤1：选取Polar卫星拍摄到的UVI极光图像，并对UVI极光图像预处理，以消除噪声点、移除无效背景像素；

步骤2：对预处理后的UVI极光图像中的极光卵边界进行提取；

步骤3：利用NASAOMNI数据，并根据每幅极光图像发生的UT时间，得到每个极光卵边界点对应的行星际和地磁参数；

步骤4：对步骤2获得的极光卵边界数据结合步骤3获得的行星际地磁参数值进行筛选，去除无效边界点，构建有效数据集；

步骤5：根据步骤4得到的有效数据集，用一元线性回归模型，对赤道向和极向极光卵边界位置进行建模，因变量取赤道向或极向极光卵边界的地磁纬度值，自变量取行星际和地磁参数；

步骤6：根据步骤5建立的极光卵边界位置模型，对任意给定的行星际和地磁条件，能够预测出此条件下的赤道向和极向极光卵边界位置。

进一步地，步骤1中对UVI极光图像预处理的方法为：

步骤1.1：移除UVI极光图像中纬度低于50度的区域；

步骤1.2：掩膜：根据UVI极光图像外边界呈椭圆状的特点，用直径为228像素的圆去截取经过步骤1.1处理后的UVI极光图像，并从中间把宽度限制到200像素，最后图像大小为228*200；

步骤1.3：清除经步骤1.2处理后的UVI极光图像中像素值为负值的点；

步骤1.4：将经过步骤1.3处理后的UVI极光图像中面积小于50个像素，极光强度大于整幅图像灰度均值加3倍方差的区域的像素值重置为整幅图像像素的灰度均值。

进一步地，步骤2中对UVI极光图像中极光卵边界进行提取的方法为：

步骤2.1：利用包含空间信息的模糊c均值聚类方法将UVI极光图像聚成6类；

步骤2.2：根据聚类结果把UVI极光图像分割为极光卵区和背景区域；

步骤2.3：对聚类后出现的与周围区域类别不同的区域，用其周围邻域的类别信息对其进行平滑处理；

步骤2.4：获取极光卵边界并将其映射到磁地方时-地磁纬度坐标系下。

进一步地，步骤2.2中根据聚类结果把UVI极光图像分割为极光卵区和背景区域的准则为：将强度最小的2个或3个聚类簇视为背景区域，其他则为极光卵区。

进一步地，步骤3中所述的行星际参数包括行星际磁场三分量Bx、By、Bz和太阳风动压Pdyn、太阳风速度Vp、太阳风密度Np；地磁参数包括极光电集流指数AE、极盖指数PC。