CN104680167B - 基于深度学习的极光卵位置确定方法 - Google Patents
基于深度学习的极光卵位置确定方法 Download PDFInfo
- Publication number
- CN104680167B CN104680167B CN201510102348.6A CN201510102348A CN104680167B CN 104680167 B CN104680167 B CN 104680167B CN 201510102348 A CN201510102348 A CN 201510102348A CN 104680167 B CN104680167 B CN 104680167B
- Authority
- CN
- China
- Prior art keywords
- image
- mtd
- aurora
- mrow
- mtr
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 32
- 238000013135 deep learning Methods 0.000 title claims abstract description 24
- 239000005441 aurora Substances 0.000 claims abstract description 114
- 230000011218 segmentation Effects 0.000 claims abstract description 31
- 238000003709 image segmentation Methods 0.000 claims abstract description 9
- 238000006243 chemical reaction Methods 0.000 claims abstract description 4
- 210000004681 ovum Anatomy 0.000 claims description 47
- 102000002322 Egg Proteins Human genes 0.000 claims description 37
- 108010000912 Egg Proteins Proteins 0.000 claims description 37
- 238000012549 training Methods 0.000 claims description 15
- 238000004422 calculation algorithm Methods 0.000 claims description 13
- 238000004364 calculation method Methods 0.000 claims description 2
- 238000003384 imaging method Methods 0.000 claims description 2
- 238000003064 k means clustering Methods 0.000 claims 1
- 230000000694 effects Effects 0.000 abstract description 8
- 238000011160 research Methods 0.000 abstract description 4
- 230000005389 magnetism Effects 0.000 abstract 3
- 238000010276 construction Methods 0.000 abstract 1
- 235000013601 eggs Nutrition 0.000 description 40
- 230000006870 function Effects 0.000 description 10
- 238000004088 simulation Methods 0.000 description 6
- 230000008859 change Effects 0.000 description 5
- 238000012360 testing method Methods 0.000 description 5
- 238000010168 coupling process Methods 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 230000004907 flux Effects 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 230000009471 action Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 230000005358 geomagnetic field Effects 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 238000005293 physical law Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000004062 sedimentation Methods 0.000 description 1
- 238000011524 similarity measure Methods 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Landscapes
- Image Analysis (AREA)
Abstract
本发明公开了一种深度学习的极光卵位置确定方法，主要解决现有极光卵位置确定方法中存在的不准确的问题，其步骤是：(1)从紫外极光图像中挑选完整的极光卵图像构成极光卵图像集；(2)利用形状信息和最大相似性区域合并准则对原始图像做图像分割；(3)利用图像的地磁坐标信息，将分割得到的图像进行坐标转换；(4)根据转换后图像得到极光卵赤道向和极向的边界点坐标；(5)利用极光卵图像集中图像的拍摄时间，构建极光卵对应时刻的地磁物理参数数据库；(6)将极光卵边界点坐标和地磁物理参数输入深度学习网络对极光卵位置进行确定。本发明提高了极光卵位置的确定精确，可用于研究地磁物理参数对极光活动的影响。
Description
技术领域
本发明属于图像处理技术领域，特别涉及极光卵位置的确定方法，可用于研究极光卵位置与地磁物理参数之间的关系，了解日地耦合过程，认知空间气候变化规律。
背景技术
了解地球周围环境的一些物理过程如太阳风中的等离子体、地球磁场和地球大气层共同作用产生的极光现象，对于地球空间科学的研究具有重要的意义。极光卵是在地球南北极地区，受到太阳风影响的高能粒子沉降所产生的环形带状发光区域。极光卵的极向边界和赤道向边界是重要的地球物理参数，其与太阳风、地磁活动有着密切的关系。对极光卵的边界位置进行确定，有助于了解日地耦合过程，认知空间气候变化规律。极光卵的边界位置是受多种地磁物理因素共同影响的，传统方法将极光卵边界位置表示为一种或几种地磁活动指数的函数，具有很大的局限性。
在早期的研究工作中人们将极光卵的位置表示为单变量的函数。1967年Feldstein提出基于原始数据的模型，将极光卵出现的可能性描绘成了Q指数的函数。Starkov将极向、赤道向和弥散极光卵边界表示为AL指数的函数，Zhang和Paxton使用了Epstein方程来计算电子通量和平均能通量，并将极光卵位置表示为Kp指数的函数，Carbary也将极光卵的边界位置用Kp指数来表示。Milan的研究表明在磁暴期间极光卵半径随环电流增强而增加。Lukianova统计了极盖区域边界与行星际磁场的关系，指出极盖区域边界在正午-午夜方向受行星际磁场z方向分量控制，在晨昏线方向受行星际磁场y方向分量影响。
由于极光卵的边界位置不是一种地磁活动指数能决定的，因而上述这些将极光卵边界位置表示为单一一种地磁活动指数的函数的方法具有局限性，目前研究的趋势是将极光卵的位置和大小用多种地磁物理参数来表征。Yang利用多元回归的方法将极光卵的位置表征为多种物理参数的函数，但多元回归的方法受回归方程假设的影响很大，二次多项式形式的回归方程并不能准确地反应客观的物理规律。
发明内容
本发明的目的在于克服上述方法的不足，提出了基于深度学习的极光卵位置确定方法，以客观准确地反应各个物理参数对不同极光现象的影响，分析空间气候变化规律。
实现本发明目的技术思路是：利用基于形状信息反馈和最大相似性区域合并准则的极光卵分割算法对原始极光卵图像做图像分割，然后对分割结果图像利用地磁时和地磁纬度信息转换图像坐标，得到在地磁时和地磁纬度坐标下的极光卵图像并取极光卵赤道向和极向的边界点坐标，最后训练深度学习网络，将极光卵赤道向和极向边界点坐标和地磁物理参数对应通过训练好的深度学习网络预测极光卵位置。其具体实现方案如下：
本发明基于深度学习的极光卵位置确定方法，包括如下步骤：
(1)选取Polar卫星携带的紫外极光成像仪所拍摄的紫外极光图像，并从中挑选完整的极光卵图像构成极光卵图像集；
(2)对极光卵图像集中的图像做图像分割：
(2a)利用分水岭算法对极光卵图像集中的图像进行初始分割，得到初始分割结果F；
(2b)利用K-means方法按灰度值大小将极光卵图像集中的图像像素点聚为三类，这三类聚类中心的灰度值按从小到大依次记为C1,C2和C3，并将位于第二灰度值C2附近的像素点和整幅图像边界上的像素点标记为背景，将位于第三灰度值C3附近的像素点标记为目标；
(2c)利用标记信息，将初始分割结果F中的所有区域划分为包含背景标记的区域集合Mb、包含目标标记的区域集合Mo和未包含任何标记的区域集合Mn，以包含背景标记的区域集合Mb和包含目标标记的区域集合Mo为初始区域，应用最大相似性区域对初始区域和其相邻区域进行合并，得到初步分割结果Ii；
(2d)应用最小二乘椭圆拟合算法将初步分割结果Ii中的极光卵边界拟合成椭圆，再用最大相似性区域合并算法对椭圆区域和其相邻区域进行合并，得到最终的分割结果Io；
(3)利用分割结果Io中极光卵区域像素点的地磁纬度Mlat值和地磁时Mlt值，将分割结果Io转换到Mlt-Mlat坐标系下，得到分割结果Io在Mlt-Mlat坐标系下的图像IM；
(4)对极光卵图像集中的所有图像经过图像分割和转换坐标系后，以图像中心为起点均匀引出24条射线，得到射线与极光卵内边界和外边界的交点，记录所有图像得到的交点的地磁纬度Mlat值并构成图像交点集合MI；
(5)根据极光卵图像集中图像的拍摄时间，从数据库中获得对应时刻的32个不同地磁物理参数的数值构成地磁物理参数集PI，将地磁物理参数集PI和图像交点集合MI按照图像的拍摄时间相对应，删除其中的无效数据，得到地磁物理参数样本集P和图像交点样本集M；
(6)根据参数样本集P和图像交点样本集M，训练得到深度学习网络。
本发明与现有技术相比具有如下优点：
1)本发明由于在图像分割阶段充分利用极光卵的形状信息，提出了基于极光卵形状信息反馈和最大相似性区域合并准则的极光卵分割，能够得到既完整又相对准确的极光卵边界；
2)本发明通过训练深度学习网络，利用深度学习的方法建立利用多种地磁物理参数确定极光卵边界位置的确定方法，能够有效找到很多的地磁活动指数与极光卵边界位置的内在联系，因而所构建的深度学习网络更能反映太阳风与地磁场的相互作用的规律。
本发明的技术过程和效果可结合以下附图详细说明。
附图说明
图1是本发明的实现流程图；
图2是本发明的极光卵图像分割流程图；
图3是本发明的极光卵位置确定流程图；
图4是本发明的极光卵边界点位置结果示意图。
具体实施方式
参照图1，本发明基于深度学习的极光卵位置确定方法，包括如下步骤：
步骤一，获取极光卵图像集。
从Polar卫星携带的紫外极光成像仪所拍摄的紫外极光图像中，选取1996年12月至1997年2月的极光卵图像，并从中挑选极光卵区域比较完整的图像，随着Polar卫星的运动导致紫外极光成像仪的拍摄角度发生变化，使得其拍摄的图像存在极光卵区域的缺损，剔除这些不完整的极光卵图像，最终得到具有完整型极光卵的图像构成极光卵图像集。
步骤二，对极光卵图像集中的图像做图像分割。
参照图2，本步骤的具体实现如下：
2.1)对极光卵图像集中的极光卵图像I，将其与Sobel算子进行卷积，得到极光卵图像I的水平梯度图像Gx和竖直梯度图像Gy
利用水平梯度图像Gx和竖直梯度图像Gy得到包含图像边缘信息的梯度模值图像G:
2.2)对梯度模值图像G直接输入分水岭算法，得到初始分割结果F；
2.3)用K-mean方法对极光卵图像I进行聚类：
2.3a)设置K-means聚类类别数为K＝3；
2.3b)选取初始聚类中心为灰度值为0、140和220的像素点，分别记为无效成像区域聚类中心C1、背景区域聚类中心C2和极光卵区域聚类中心C3；
2.3c)计算其余像素点到各聚类中心的欧氏距离，找出与其距离最近的聚类中心，并将其归类到该聚类中心所属的类别，迭代计算得到极光卵图像I所有像素点的归类结果；
2.3d)对归类后得到的三个类别簇，分别计算三个类别簇内像素点的灰度值均值，并把距离均值最近的像素点替代原有的聚类中心C1、C2和C3；
2.3f)重复步骤2.3.3)和2.3.4)直到聚类中心C1、C2和C3不再发生变化，得到极光卵图像I的聚类结果；
2.4)利用极光卵图像I的聚类结果，将位于聚类中心C2附近的像素点和整幅图像边界上的像素点标记为背景，将位于聚类中心C3附近的像素点标记为目标；
2.5)利用标记信息，将初始分割结果F中的所有区域划分到不同的集合：若区域内包含背景标记，则划分到包含背景标记的区域集合Mb中，若区域内包含目标标记，则划分到包含目标标记的区域集合Mo中，若区域未包含任何标记，则划分到未包含任何标记的区域集合Mn中；
2.6)对上述这三个区域集合内的区域进行区域合并：
2.6a)对包含目标标记的区域集合Mo中的区域O∈Mo，将其邻域的集合记作r为邻域的个数且
2.6b)选取集合中的一个区域Ai，将区域Ai的邻域集合记作k为区域Ai邻域的个数；
2.6c)遍历搜索区域Ai的邻域，若区域Ai满足：的条件，则将区域Ai合并到区域O，其中，ρ为相似性度量函数；
2.6d)对区域O的所有邻域{Ai}i＝1,2,...,r重复执行步骤2.6b)和2.6c)直到没有满足合并条件的区域为止，将合并到区域O的区域添加到包含目标标记的区域集合Mo内；
2.6e)对未包含任何标记的区域集合Mn中的区域N∈Mn，其邻域的集合记作p为邻域个数；
2.6f)选取集合中的一个区域Hi，对区域Hi的邻域集合记作q为区域Hi的邻域个数；
2.6g)遍历搜索区域Hi邻域，若区域Hi满足：的条件，其中，ρ为相似性度量函数，则将区域H1并入区域N；
2.6h)对区域N的所有邻域{Hi}i＝1,2,...,p重复执行步骤2.6f)和2.6g)直到没有满足合并条件的区域为止，将合并到区域N的区域添加到未包含任何标记的区域集合Mn；
2.6i)对包含目标标记的区域集合MO和未包含任何标记的区域集合Mn内的所有区域，重复执行步骤2.6a)到2.6h)，直到没有满足合并条件的区域为止，得到初始分割图像Ii；
2.7)应用最小二乘椭圆拟合算法拟合极光卵边界：
2.7a)提取初步分割结果Ii中极光卵的所有边界点的坐标值(xi,yi)i＝1,2,...n，n为边界点的个数，对所有坐标点进行归一化后得到坐标向量：
其中T表示转置；
2.7b)定义椭圆方程：
aX2+bXY+cY2+dX+eY+f＝0，
其中{a b c d e f}为椭圆方程各项的参数，将坐标向量代入椭圆方程求解出这些参数{a b c d e f}；
2.7c)利用上述椭圆参数{a b c d e f}，得到极光卵边界拟合的椭圆形状；
2.8)将极光卵边界拟合的椭圆形状标记为目标区域；
2.9)应用最大相似性区域合并算法对目标区域和其相邻区域进行合并，得到最终的分割结果图像Io。
步骤三，分割图像坐标转换。
利用分割结果图像Io中极光卵区域像素点的地磁纬度Mlat值和地磁时Mlt值，将分割结果Io转换到Mlt-Mlat坐标系下，对于转换坐标后缺少的像素点，取其临近位置的灰度值作为像素点灰度值，得到分割结果Io在Mlt-Mlat坐标系下的图像IM。
步骤四，选取极光卵边界点坐标。
以上述图像IM的中心为起点均匀引出24条射线，其分别对应24个地磁时Mlt值，24条射线分别与极光卵边界有交点，得到射线与极光卵内边界和外边界的交点，并把交点的地磁纬度Mlat值构成交点集合MI。
步骤五，选取物理参数。
根据极光卵图像集中图像的拍摄时间，从OMNI数据库中获得对应时刻的32个不同地磁物理参数的数值，包括太阳风压、太阳风速等参数，构成地磁物理参数集PI，将地磁物理参数集PI和图像交点集合MI按照图像的拍摄时间相对应，删除其中由于卫星拍摄角度变化导致的无效数据，得到由4432组样本数据构成的地磁物理参数样本集P和图像交点样本集M。
步骤六，根据地磁物理参数样本集P和图像交点样本集M，训练得到深度学习网络。
参考图3，本步骤具体操作如下：
6.1)将两层受限玻尔兹曼机RBM网络和一层径向基RBF网络进行组合，并将地磁物理参数样本集P作为该组合网络的输入，将图像交点样本集M作为该组合网络的输出；
6.2)设置组合网络的训练误差为ε＝4，并对其进行训练，即修改该组合网络中间层之间的连接权值，使得该组合网络的输出结果与图像交点样本集M之间的输出误差越来越小，直到输出误差小于设定的误差ε时结束训练；
6.3)训练结束后，固定该组合网络的中间层之间的连接权值，则该组合网络就为表示地磁物理参数和极光卵边界位置关系的深度学习网络。
本发明的效果通过一下仿真进一步说明：
1.仿真实验数据：
本发明所使用的数据来自Polar卫星携带的紫外极光成像仪所拍摄的1996年12月至1997年2月间的紫外极光图像。为了验证本发明的有效性，在紫外极光图像上进行了仿真实验。
2.仿真参数
实验中对利用深度学习网络确定极光卵边界位置的准确性，主要的评价准则是平均绝对误差：
其中，为第j个测试样本在第i个磁地方时处实际的地磁纬度值，为第j个测试样本在第i磁地方时处通过网络学习得到的地磁纬度值，k为测试样本总数。
2.仿真内容和结果
从地磁物理参数样本集P和图像交点样本集M随机选取3500组数据作为训练集，设置深度学习网络的隐含层节点数为64，训练误差设置为4，得到深度学习网络，并利用该网络对剩下的932组数据进行测试，进行100次测试，得到利用深度学习网络得到的极光卵边界点位置。
将上述仿真结果进行可视化表示，将通过网络确定的极光卵边界点位置和实际的极光卵边界点位置，同时在极光卵图像上进行表示，结果如图4所示，图4中圆圈代表极光卵内边界点的实际位置，加号代表利用该发明确定的极光卵内边界点的位置，正方形代表极光卵外边界点的实际位置，三角形代表利用该发明确定的极光卵外边界点的位置。
3.仿真结果分析，将利用深度学习网络得到的极光卵边界点位置与实际的极光卵边界点位置做对比，得到极光卵内边界点位置的平均绝对误差为1.8，外边界点平均绝对误差为1.6，可以看出本发明对极光卵边界点位置的确定有着较高的精度。
Claims (4)
1.一种基于深度学习的极光卵位置确定方法，包括如下步骤：
(1)选取Polar卫星携带的紫外极光成像仪所拍摄的紫外极光图像，并从中挑选完整的极光卵图像构成极光卵图像集；
(2)对极光卵图像集中的图像做图像分割：
(2a)利用分水岭算法对极光卵图像集中的图像进行初始分割，得到初始分割结果F；
(2b)利用K-means方法按灰度值大小将极光卵图像集中的图像像素点聚为三类，这三类聚类中心的灰度值按从小到大依次记为C1,C2和C3，并将位于第二灰度值C2附近的像素点和整幅图像边界上的像素点标记为背景，将位于第三灰度值C3附近的像素点标记为目标；
所述利用K-means方法按灰度值大小将极光卵图像集中的图像像素点聚为三类，其实现如下：
(2b1)设置K-means聚类方法聚类的类别数为K＝3；
(2b2)选取初始聚类中心为灰度值为0、140和220的像素点，分别记为无效成像区域聚类中心C1、背景区域聚类中心C2和极光卵区域聚类中心C3；
(2b3)计算其余像素点与各聚类中心的欧氏距离，将其归类到距离最近的类别，迭代计算得到图像所有像素点的归类结果；
(2b4)利用均值最小化准则对聚类中心C1、C2和C3进行更新；
(2b5)重复步骤(2b3)和(2b4)直到聚类中心不再发生变化，得到图像聚类结果；
(2c)利用标记信息，将初始分割结果F中的所有区域划分为包含背景标记的区域集合Mb、包含目标标记的区域集合Mo和未包含任何标记的区域集合Mn，以包含背景标记的区域集合Mb和包含目标标记的区域集合Mo为初始区域，应用最大相似性区域对初始区域和其相邻区域进行合并，得到初步分割结果Ii；
(2d)应用最小二乘椭圆拟合算法将初步分割结果Ii中的极光卵边界拟合成椭圆，再用最大相似性区域合并算法对椭圆区域和其相邻区域进行合并，得到最终的分割结果Io；
(3)利用分割结果Io中极光卵区域像素点的地磁纬度Mlat值和地磁时Mlt值，将分割结果Io转换到Mlt-Mlat坐标系下，得到分割结果Io在Mlt-Mlat坐标系下的图像IM；
(4)对极光卵图像集中的所有图像经过图像分割和转换坐标系后，以图像中心为起点均匀引出24条射线，得到射线与极光卵内边界和外边界的交点，记录所有图像得到的交点的地磁纬度Mlat值并构成图像交点集合MI；
(5)根据极光卵图像集中图像的拍摄时间，从数据库中获得对应时刻的32个不同地磁物理参数的数值构成地磁物理参数集PI，将地磁物理参数集PI和图像交点集合MI按照图像的拍摄时间相对应，删除其中的无效数据，得到地磁物理参数样本集P和图像交点样本集M；
(6)根据参数样本集P和图像交点样本集M，训练得到深度学习网络。
2.根据权利要求1所述的方法，其中步骤(2a)所述的利用分水岭算法对极光卵图像集中的图像进行初始分割，得到初始分割结果F，按如下步骤进行：
(2a1)对极光卵图像集中的图像I与Sobel算子进行卷积，得到图像的水平梯度图像Gx和竖直梯度图像Gy：
<mfenced open = "" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>G</mi>
<mi>x</mi>
</msub>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>1</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<mn>2</mn>
</mrow>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>2</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>1</mn>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>&CircleTimes;</mo>
<mi>I</mi>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>G</mi>
<mi>y</mi>
</msub>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<mn>2</mn>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>1</mn>
</mtd>
<mtd>
<mn>2</mn>
</mtd>
<mtd>
<mn>1</mn>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>&CircleTimes;</mo>
<mi>I</mi>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
(2a2)由水平梯度图像Gx和竖直梯度图像Gy得到包含图像边缘信息的图像G；
<mrow>
<mi>G</mi>
<mo>=</mo>
<msqrt>
<mrow>
<msubsup>
<mi>G</mi>
<mi>x</mi>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>G</mi>
<mi>y</mi>
<mn>2</mn>
</msubsup>
</mrow>
</msqrt>
</mrow>
(2a3)将图像G作为分水岭算法的输入，得到初始分割结果F。
3.根据权利要求1所述的方法，其中步骤(2d)所述的应用最小二乘椭圆拟合算法将初步分割结果Ii中的极光卵边界拟合成椭圆，按如下步骤进行：
(2d1)提取初步分割结果Ii中极光卵的所有边界点的坐标值(xi,yi)i＝1,2,...n，n为边界点的个数，对所有坐标点进行归一化后得到坐标向量：
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>X</mi>
<mo>=</mo>
<msup>
<mrow>
<mo>&lsqb;</mo>
<msub>
<mi>x</mi>
<mn>1</mn>
</msub>
<mo>,</mo>
<msub>
<mi>x</mi>
<mn>2</mn>
</msub>
<mo>,</mo>
<mo>...</mo>
<msub>
<mi>x</mi>
<mi>n</mi>
</msub>
<mo>&rsqb;</mo>
</mrow>
<mi>T</mi>
</msup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>Y</mi>
<mo>=</mo>
<msup>
<mrow>
<mo>&lsqb;</mo>
<msub>
<mi>y</mi>
<mn>1</mn>
</msub>
<mo>,</mo>
<msub>
<mi>y</mi>
<mn>2</mn>
</msub>
<mo>,</mo>
<mo>...</mo>
<msub>
<mi>y</mi>
<mi>n</mi>
</msub>
<mo>&rsqb;</mo>
</mrow>
<mi>T</mi>
</msup>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
其中T表示转置；
(2d2)定义椭圆方程：
aX2+bXY+cY2+dX+eY+f＝0
其中{a b c d e f}为实数，将坐标向量代入椭圆方程求解椭圆参数{a b c d e f}；
(2d3)利用得到的椭圆参数{a b c d e f}，得到极光卵边界拟合的椭圆。
4.根据权利要求1所述的方法，其中步骤(6)所述的根据参数样本集P和图像交点样本集M，训练得到深度学习网络，按如下步骤进行：
(6a)将两层受限玻尔兹曼机RBM网络和一层径向基RBF网络进行组合，并将地磁物理参数样本集P作为该组合网络的输入，将图像交点样本集M作为该组合网络的输出；
(6b)设置组合网络的训练误差为ε＝4，并对其进行训练，即修改该组合网络中间层之间的连接权值，使得该组合网络的输出结果与图像交点样本集M之间的输出误差越来越小，直到输出误差小于设定的误差ε时结束训练；
(6c)训练结束后，固定该组合网络的中间层之间的连接权值，则该组合网络就为表示地磁物理参数和极光卵边界位置关系的深度学习网络。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510102348.6A CN104680167B (zh) | 2015-03-09 | 2015-03-09 | 基于深度学习的极光卵位置确定方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510102348.6A CN104680167B (zh) | 2015-03-09 | 2015-03-09 | 基于深度学习的极光卵位置确定方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104680167A CN104680167A (zh) | 2015-06-03 |
| CN104680167B true CN104680167B (zh) | 2018-03-20 |
Family
ID=53315184
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510102348.6A Active CN104680167B (zh) | 2015-03-09 | 2015-03-09 | 基于深度学习的极光卵位置确定方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104680167B (zh) |
Families Citing this family (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105118047A (zh) * | 2015-07-15 | 2015-12-02 | 陕西师范大学 | 基于行星际和地磁参数的极光卵边界位置预测方法 |
| CN105069798A (zh) * | 2015-08-12 | 2015-11-18 | 西安电子科技大学 | 基于图割优化的局部主动轮廓图像分割方法 |
| CN107186708B (zh) * | 2017-04-25 | 2020-05-12 | 珠海智卓投资管理有限公司 | 基于深度学习图像分割技术的手眼伺服机器人抓取系统及方法 |
| CN108347643B (zh) * | 2018-03-05 | 2020-09-15 | 成都索贝数码科技股份有限公司 | 一种基于深度学习的字幕叠加截图的实现方法 |
| CN110188612B (zh) * | 2019-04-28 | 2023-02-10 | 西安电子科技大学 | 基于生成式对抗网络的极光卵强度图像建模方法 |
| CN113344032B (zh) * | 2021-05-13 | 2021-12-28 | 中国科学院国家空间科学中心 | 一种基于相似度算法的地磁Kp指数预报方法及系统 |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102930273A (zh) * | 2012-10-15 | 2013-02-13 | 西安电子科技大学 | 基于亮度自适应水平集的极光卵分割方法 |
-
2015
- 2015-03-09 CN CN201510102348.6A patent/CN104680167B/zh active Active
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102930273A (zh) * | 2012-10-15 | 2013-02-13 | 西安电子科技大学 | 基于亮度自适应水平集的极光卵分割方法 |
Non-Patent Citations (3)
| Title |
|---|
| A shape-initialized and intensity-adaptive level set method for auroral oval segmentation;Xi Yang等;《Information Science》;20140325;794-807 * |
| An Automatic MSRM Method with a Feedback Based on Shape Information for Auroral Oval Segmentation;Hui Liu等;《Intelligence Science and Big Data Engineering》;20131231;748-755 * |
| Z.-J Hu等.Synoptic distribution of dayside aurora- Multiple-wavelength all-sky observation at Yellow River Station in Ny-Ålesund, Svalbard.《Journal of Atmospheric and Solar-Terrestrial Physics》.2009, * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN104680167A (zh) | 2015-06-03 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104680167B (zh) | 基于深度学习的极光卵位置确定方法 | |
| Lu et al. | Cultivated land information extraction in UAV imagery based on deep convolutional neural network and transfer learning | |
| CN106599922B (zh) | 用于大规模数据标定的迁移学习方法及系统 | |
| CN108038445B (zh) | 一种基于多视角深度学习框架的sar自动目标识别方法 | |
| CN111985376A (zh) | 一种基于深度学习的遥感影像舰船轮廓提取方法 | |
| CN110132263B (zh) | 一种基于表示学习的星图识别方法 | |
| CN109029363A (zh) | 一种基于深度学习的目标测距方法 | |
| CN110334645B (zh) | 一种基于深度学习的月球撞击坑识别方法 | |
| CN106709421B (zh) | 一种基于变换域特征和cnn的细胞图像识别分类方法 | |
| CN109284779A (zh) | 基于深度全卷积网络的物体检测方法 | |
| CN114332348B (zh) | 一种融合激光雷达与图像数据的轨道三维重建方法 | |
| CN114612835A (zh) | 一种基于YOLOv5网络的无人机目标检测模型 | |
| CN112348758B (zh) | 一种光学遥感图像数据增强方法及目标识别方法 | |
| CN108776777A (zh) | 一种基于Faster RCNN的遥感影像对象间空间关系的识别方法 | |
| Herrault et al. | A comparative study of geometric transformation models for the historical" Map of France" registration | |
| CN112329559A (zh) | 一种基于深度卷积神经网络的宅基地目标检测方法 | |
| CN117765313A (zh) | 一种基于空间连续视图对比学习的多源遥感图像分类方法 | |
| CN104268581A (zh) | 一种基于综合像元级和亚像元级空间相关性特征的遥感亚像元制图方法 | |
| Xu et al. | A co-training approach to the classification of local climate zones with multi-source data | |
| CN115311502A (zh) | 基于多尺度双流架构的遥感图像小样本场景分类方法 | |
| CN106991411A (zh) | 基于深度形状先验的遥感图像目标精细化提取方法 | |
| CN105654122A (zh) | 基于核函数匹配的空间金字塔物体识别方法 | |
| CN107045722A (zh) | 融合静态信息与动态信息的视频序列分割方法 | |
| CN115690627A (zh) | 一种航拍图像旋转目标的检测方法及系统 | |
| Yildiz et al. | Amstertime: A visual place recognition benchmark dataset for severe domain shift |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |