一种云微粒子数据区域提取方法

allin2023-03-17  134



1.本发明属于图像处理及气象领域,具体涉及对机载云微粒子探测设备获取的数据进行数据划分,并基于聚类方法和同区域搜索方法对云微粒子区域进行确定和提取。


背景技术:

2.气候变化影响着云的形成和发展,全球范围内的人工影响天气工作正在努力增加降水并减轻冰雹。为了更有效的进行人工影响天气作业,需要有关特定地区和季节特定的云属性的详细信息。为了知道何时何地播种云,了解特定的气溶胶和云特性至关重要,所以需要利用探头仪器获得液态水含量、凝结核、液滴大小、液滴数量、气溶胶数量,质量浓度等信息。
3.目前,国内外气象部门和科研机构都已经在利用机载云粒子成像探头对云层内部的微物理结构进行直接观测。相比于利用雷达等遥感设备对云进行观测,云粒子成像仪获取的云信息更加准确和直观,并且机载探测设备所记录的粒子图像是云微物理研究的基础,
4.然而,由于云微粒子探测设备获取到的云微粒子原始数据中存在破碎粒子、虚假目标、粒子遮挡等问题,如果不对这些问题进行处理直接进行数据分析会使得云微粒子分类的准确性以及利用云微粒子分类结果来进行气象研究的可靠性降低。而当前对于粒子区域的提取主要是基于最小外接矩形、形态学处理、模板匹配等方法来实现的,这些方法对于处理云微粒子原始数据中存在的问题效果并不理想,存在破碎粒子、虚假目标干扰严重、粒子遮挡无法处理的情况。


技术实现要素:

5.本发明针对上述现有技术的不足,提出了云微粒子数据区域提取方法,包括以下步骤:
6.一种云微粒子数据区域提取方法,主要包括以下步骤:
7.步骤1:对云微粒子探测设备采集的原始云微粒子图像数据进行数据划分,得到相应数据划分结果;
8.步骤2:对步骤1所得的每个dl
n,m
中的所有图像数据进行图像聚类方法处理,得到不同批次、数据块的dl
n,m
数据集的处理后的包含像素块的图像数据结果sf
n,m,j
,其中变量n表示批次编号,m表示数据块编号,j表示相同批次和数据块内的图像数据编号,j∈[1,j],变量j为dl
n,m
中图像数据数量;
[0009]
步骤3:对所有图像数据sf
n,m,j
进行像素块同区域搜索:首先根据像素块直方图统计分布特征对像素块属性进行初步确认,将部分像素块属性标记为中心像素块或附属像素块,然后对初步确认的中心像素块和附属像素块分别进行中心像素块同区域搜索和附属像素块同区域搜索,对初步确认的像素块属性进行进一步确认,并对未确认属性的像素块进行属性确认;
[0010]
步骤4:遍历所有图像数据sf
n,m,j
中的所有中心像素块,以中心像素块为核心,以归属于该中心像素块的附属像素块为扩展区域,两者共同作为单个云微粒子区域,若某个附属像素块同时归属多个中心像素块,重复将其归属于不同中心像素块,组成不同单个云微粒子区域,最后以单个云微粒子区域为对象进行云微粒子提取。
[0011]
进一步的,步骤1对云微粒子探测设备采集的原始云微粒子图像数据进行数据划分具体为:
[0012]
首先对云微粒子探测设备采集的原始云微粒子图像数据dl
all
进行批次划分:按照云微粒子探测设备采集数据的时间段不同将dl
all
划分为n个批次,dln为第n个批次的原始云微粒子图像数据,n∈[1,n],然后对每个dln进行数据块划分:按照云微粒子探测设备采集数据时所处海拔局度不同将dln划分为m个数据块,dl
n,m
为第n个批次中的第m个数据块,m∈[1,m]。
[0013]
进一步的,步骤3:对所有图像数据sf
n,m,j
进行像素块同区域搜索包括:
[0014]
步骤3.1:按照从左到右,从上到下的规则遍历图像数据sf
n,m,j
中的所有像素块,并对像素块进行编号得到sf
n,m,j
(q),q表示图像中像素块编号,q∈[1,q],变量q为图像中的像素块数量;
[0015]
步骤3.2:对步骤3.1所得的像素块sf
n,m,j
(q)生成其对应的直方图h
n,m,j,q
,将满足相同n、m、j条件下的h
n,m,j,q
分为一组,得到n
·m·
j组;
[0016]
步骤3.3:对步骤3.2所得的n
·m·
j组直方图h
n,m,j,q
,设置阈值gv,统计每个直方图中灰度值大于等于gv的区域中的最大纵坐标值和小于gv的区域中的最大纵坐标值,并计算前者与后者的比值rh

n,m,j,q

[0017]
步骤3.4:对步骤3.3所得的n
·m·
j组比值数据集进行分析,设置阈值th

,统计每组比值数据集里小于阈值的比值rh

n,m,j,q

,q

表示比值数据集中满足rh

n,m,j,q
<th

条件的筛选出的比值对应的图像中像素块编号;
[0018]
步骤3.5:对步骤3.4所得的n
·m·
j组筛选出的像素块,按照逐组处理的方式,通过每组像素块在图像数据中的空间分布,将这些筛选出的像素块设为中心像素块,并将其8邻域像素块设为该中心像素块的附属像素块,如果出现某附属像素块同时属于多个中心像素块,将其随机分配给其中一个中心像素块;按照从左到右,从上到下的规则,对所有中心像素块进行步骤3.6至步骤3.8的处理;
[0019]
步骤3.6:对中心像素块进行同区域搜索,搜索规则为:从中心像素块上方开始,按照顺时针顺序,逐一搜索中心像素块8邻域像素块;
[0020]
步骤3.7:对附属像素块进行同区域搜索,搜索规则为:从附属像素块上方开始,按照顺时针顺序,逐一搜索附属像素块8邻域像素块;
[0021]
步骤3.8:重复步骤3.7对所有轮次附属像素块进行附属像素块同区域搜索直到搜索8邻域没有附属像素块。
[0022]
进一步的,所述步骤3.6具体为:
[0023]
对搜索到的某一邻域像素块,首先判断其属性是否为中心像素块,如果是,执行步骤3.6.1,如果不是,判断其是否为其他中心像素块的附属像素块,如果是,执行步骤3.6.2,如果不是,执行步骤3.6.3;当搜索完8邻域范围结束;
[0024]
步骤3.6.1:更改该邻域像素块属性为其依附的中心像素块的第1轮附属像素块,
且将其所依附的中心像素块的编号和轮次数附加到该像素块编号之后作为该附属像素块的标记号,继续搜索;
[0025]
步骤3.6.2:设置阈值th

,如果该邻域像素块的rh

n,m,j,q
满足条件rh

n,m,j,q
<th

,在其原有属性基础上增加一项属性为当前中心像素的第1轮附属像素块,且将其当前附属的中心像素块编号和轮次数增加到该像素块已有标记号后面作为该附属像素块新增标记号;如果该邻域像素块的rh

n,m,j,q
不满足条件rh

n,m,j,q
<th

,跳过该像素块,继续搜索;
[0026]
步骤3.6.3:基于设置的阈值th

,如果邻域像素块的rh

n,m,j,q
满足条件rh

n,m,j,q
<th

,更改该邻域像素块属性为其依附的中心像素块的第1轮附属像素块,且将其附属的中心像素块编号和轮次数附加到该像素块编号之后作为该附属像素块的标记号,继续搜索;如果邻域像素块的rh

n,m,j,q
不满足条件rh

n,m,j,q
<th

,跳过该像素块,继续搜索。
[0027]
进一步的,所述步骤3.7具体为:
[0028]
对搜索到的某一邻域像素块,首先判断该邻域像素块编号是否与此附属像素块标记号中的附属的中心像素块编号相同,如果相同,跳过这些邻域像素块,如果不相同,执行步骤3.7.1;
[0029]
步骤3.7.1:设置阈值th

,如果邻域像素块的rh

n,m,j,q
满足条件rh

n,m,j,q
<th

,将其设置为第e轮附属像素块,且将其附属的前一轮附属像素块标记号作为该附属像素块标记号,并且将轮次数更新为该轮的轮次数,继续搜索;如果邻域像素块的rh

n,m,j,q
不满足条件rh

n,m,j,q
<th

,跳过该像素块,继续搜索,当搜索完8邻域范围时结束搜索,其中,变量e表示附属像素块是从步骤3.5得到的中心像素块开始,第e次8邻域搜索得到的附属像素块,e∈[2,e],变量e为从中心像素块开始,搜索完所有附属像素块的轮次。
[0030]
本发明与传统的云微粒子图像提取及分类方法相比,存在以下优点,从而解决了相应技术问题:
[0031]
1、对云微粒子数据进行了基于直方图分布比例的像素块筛选,初步筛选了云微粒子数据聚类结果中的中心像素块,确保了云微粒子数据区域提取的可靠性。
[0032]
2、对云微粒子数据进行了基于中心像素块的同区域搜索,提高了云微粒子数据区域提取的准确性。
[0033]
3、对云微粒子数据进行了基于附属像素块的同区域搜索,提高了云微粒子数据区域提取的完整性和准确性。
附图说明
[0034]
图1为一种微粒子数据区域提取方法流程图

具体实施方式
[0035]
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整的描述,方法流程图如图1所示,包括以下步骤:
[0036]
一种云微粒子数据区域提取方法,主要包括以下步骤:
[0037]
步骤1:对云微粒子探测设备采集的原始云微粒子图像数据进行数据划分,得到相应数据划分结果;
[0038]
步骤2:对步骤1所得的每个dl
n,m
中的所有图像数据进行图像聚类方法处理,得到
不同批次、数据块的dl
n,m
数据集的处理后的包含像素块的图像数据结果sf
n,m,j
,其中变量n表示批次编号,m表示数据块编号,j表示相同批次和数据块内的图像数据编号,j∈[1,j],变量j为dl
n,m
中图像数据数量;
[0039]
步骤3:对所有图像数据sf
n,m,j
进行像素块同区域搜索:首先根据像素块直方图统计分布特征对像素块属性进行初步确认,将部分像素块属性标记为中心像素块或附属像素块,然后对初步确认的中心像素块和附属像素块分别进行中心像素块同区域搜索和附属像素块同区域搜索,对初步确认的像素块属性进行进一步确认,并对未确认属性的像素块进行属性确认;
[0040]
步骤4:遍历所有图像数据sf
n,m,j
中的所有中心像素块,以中心像素块为核心,以归属于该中心像素块的附属像素块为扩展区域,两者共同作为单个云微粒子区域,若某个附属像素块同时归属多个中心像素块,重复将其归属于不同中心像素块,组成不同单个云微粒子区域,最后以单个云微粒子区域为对象进行云微粒子提取。
[0041]
进一步的,步骤1对云微粒子探测设备采集的原始云微粒子图像数据进行数据划分具体为:
[0042]
首先对云微粒子探测设备采集的原始云微粒子图像数据dl
all
进行批次划分:按照云微粒子探测设备采集数据的时间段不同将dl
all
划分为n个批次,dln为第n个批次的原始云微粒子图像数据,n∈[1,n],然后对每个dln进行数据块划分:按照云微粒子探测设备采集数据时所处海拔高度不同将dln划分为m个数据块,dl
n,m
为第n个批次中的第m个数据块,m∈[1,m]。
[0043]
进一步的,步骤3:对所有图像数据sf
n,m,j
进行像素块同区域搜索包括:
[0044]
步骤3.1:按照从左到右,从上到下的规则遍历图像数据sf
n,m,j
中的所有像素块,并对像素块进行编号得到sf
n,m,j
(q),q表示图像中像素块编号,q∈[1,q],变量q为图像中的像素块数量;
[0045]
步骤3.2:对步骤3.1所得的像素块sf
n,m,j
(q)生成其对应的直方图h
n,m,j,q
,将满足相同n、m、j条件下的h
n,m,j,q
分为一组,得到n
·m·
j组;
[0046]
步骤3.3:对步骤3.2所得的n
·m·
j组直方图h
n,m,j,q
,设置阈值gv,统计每个直方图中灰度值大于等于gv的区域中的最大纵坐标值和小于gv的区域中的最大纵坐标值,并计算前者与后者的比值rh

n,m,j,q
,gv∈[100,120];
[0047]
步骤3.4:对步骤3.3所得的n
·m·
j组比值数据集进行分析,设置阈值th

,统计每组比值数据集里小于阈值的比值rh

n,m,j,q

,q

表示比值数据集中满足rh

n,m,j,q
<th

条件的筛选出的比值对应的图像中像素块编号,th

∈[0.15,0.25];
[0048]
步骤3.5:对步骤3.4所得的n
·m·
j组筛选出的像素块,按照逐组处理的方式,通过每组像素块在图像数据中的空间分布,将这些筛选出的像素块设为中心像素块,并将其8邻域像素块设为该中心像素块的附属像素块,如果出现某附属像素块同时属于多个中心像素块,将其随机分配给其中一个中心像素块;按照从左到右,从上到下的规则,对所有中心像素块进行步骤3.6至步骤3.8的处理;;
[0049]
步骤3.6:对中心像素块进行同区域搜索,搜索规则为:从中心像素块上方开始,按照顺时针顺序,逐一搜索中心像素块8邻域像素块;对搜索到的某一邻域像素块,首先判断其属性是否为中心像素块,如果是,执行步骤3.6.1,如果不是,判断其是否为其他中心像素
块的附属像素块,如果是,执行步骤3.6.2,如果不是,执行步骤3.6.3;当搜索完8邻域范围结束;
[0050]
步骤3.7:对附属像素块进行同区域搜索,搜索规则为:从附属像素块上方开始,按照顺时针顺序,逐一搜索附属像素块8邻域像素块;对搜索到的某一邻域像素块,首先判断该邻域像素块编号是否与此附属像素块标记号中的附属的中心像素块编号相同,如果相同,跳过这些邻域像素块,如果不相同,执行步骤3.7.1;
[0051]
步骤3.8:重复步骤3.7对所有轮次附属像素块进行附属像素块同区域搜索直到搜索8邻域没有附属像素块。
[0052]
进一步的,所述步骤3.6具体为:
[0053]
步骤3.6.1:更改该邻域像素块属性为其依附的中心像素块的第1轮附属像素块,且将其所依附的中心像素块的编号和轮次数附加到该像素块编号之后作为该附属像素块的标记号,继续搜索;
[0054]
步骤3.6.2:设置阈值th

,如果该邻域像素块的rh

n,m,j,q
满足条件rh

n,m,j,q
<th

,在其原有属性基础上增加一项属性为当前中心像素的第1轮附属像素块,且将其当前附属的中心像素块编号和轮次数增加到该像素块已有标记号后面作为该附属像素块新增标记号;如果该邻域像素块的rh

n,m,j,q
不满足条件rh

n,m,j,q
<th

,跳过该像素块,继续搜索,th

∈[1.5,2];
[0055]
步骤3.6.3:基于设置的阈值th

,如果邻域像素块的rh

n,m,j,q
满足条件rh

n,m,j,q
<th

,更改该邻域像素块属性为其依附的中心像素块的第1轮附属像素块,且将其附属的中心像素块编号和轮次数附加到该像素块编号之后作为该附属像素块的标记号,继续搜索;如果邻域像素块的rh

n,m,j,q
不满足条件rh

n,m,j,q
<th

,跳过该像素块,继续搜索。
[0056]
进一步的,所述步骤3.7具体为:
[0057]
步骤3.7.1:设置阈值th

,如果邻域像素块的rh

n,m,j,q
满足条件rh

n,m,j,q
<th

,将其设置为第e轮附属像素块,且将其附属的前一轮附属像素块标记号作为该附属像素块标记号,并且将轮次数更新为该轮的轮次数,继续搜索;如果邻域像素块的rh

n,m,j,q
不满足条件rh

n,m,j,q
<th

,跳过该像素块,继续搜索,当搜索完8邻域范围时结束搜索,其中,变量e表示附属像素块是从步骤3.5得到的中心像素块开始,第e次8邻域搜索得到的附属像素块,e∈[2,e],变量e为从中心像素块开始,搜索完所有附属像素块的轮次。
[0058]
显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的其他实施例,都属于本发明保护的范围。

技术特征:
1.一种云微粒子数据区域提取方法,主要包括以下步骤:步骤1:对云微粒子探测设备采集的原始云微粒子图像数据进行数据划分,得到相应数据划分结果dl
n,m
;步骤2:对数据划分结果dl
n,m
中的所有图像数据进行图像聚类方法处理,得到不同批次、数据块的dl
n,m
数据集的处理后的包含像素块的图像数据结果sf
n,m,j
,其中变量n表示批次编号,m表示数据块编号,j表示相同批次和数据块内的图像数据编号,j∈[1,j],变量j为dl
n,m
中图像数据数量;步骤3:对所有图像数据sf
n,m,j
进行像素块同区域搜索:首先根据像素块直方图统计分布特征对像素块属性进行初步确认,将部分像素块属性标记为中心像素块或附属像素块,然后对初步确认的中心像素块和附属像素块分别进行中心像素块同区域搜索和附属像素块同区域搜索,对初步确认的像素块属性进行进一步确认,并对未确认属性的像素块进行属性确认;步骤4:遍历所有图像数据sf
n,m,j
中的所有中心像素块,以中心像素块为核心,以归属于该中心像素块的附属像素块为扩展区域,两者共同作为单个云微粒子区域,若某个附属像素块同时归属多个中心像素块,重复将其归属于不同中心像素块,组成不同单个云微粒子区域,最后以单个云微粒子区域为对象进行云微粒子提取。2.如权利要求1所述的一种云微粒子数据区域提取方法,其特征在于步骤1对云微粒子探测设备采集的原始云微粒子图像数据进行数据划分具体为:首先对云微粒子探测设备采集的原始云微粒子图像数据dl
all
进行批次划分:按照云微粒子探测设备采集数据的时间段不同将dl
all
划分为n个批次,dl
n
为第n个批次的原始云微粒子图像数据,n∈[1,n],然后对每个dl
n
进行数据块划分:按照云微粒子探测设备采集数据时所处海拔高度不同将dl
n
划分为m个数据块,dl
n,m
为第n个批次中的第m个数据块,m∈[1,m]。3.如权利要求1所述的一种云微粒子数据区域提取方法,其特征在于步骤3对所有图像数据sf
n,m,j
进行像素块同区域搜索具体为:步骤3.1:按照从左到右,从上到下的规则遍历图像数据sf
n,m,j
中的所有像素块,并对像素块进行编号得到sf
n,m,j
(q),q表示图像中像素块编号,q∈[1,q],变量q为图像中的像素块数量;步骤3.2:对步骤3.1所得的像素块sf
n,m,j
(q)生成其对应的直方图h
n,m,j,q
,将满足相同n、m、j条件下的h
n,m,j,q
分为一组,得到n
·
m
·
j组;步骤3.3:对步骤3.2所得的n
·
m
·
j组直方图h
n,m,j,q
,设置阈值gv,统计每个直方图中灰度值大于等于gv的区域中的最大纵坐标值和小于gv的区域中的最大纵坐标值,并计算前者与后者的比值rh

n,m,j,q
;步骤3.4:对步骤3.3所得的n
·
m
·
j组比值数据集进行分析,设置阈值th

,统计每组比值数据集里小于阈值的比值rh

n,m,j,q

,q

表示比值数据集中满足rh

n,m,j,q
<th

条件的筛选出的比值对应的图像中像素块编号;步骤3.5:对步骤3.4所得的n
·
m
·
j组筛选出的像素块,按照逐组处理的方式,通过每组像素块在图像数据中的空间分布,将这些筛选出的像素块设为中心像素块,并设置该中心像素块的附属像素块;按照从左到右,从上到下的规则,对所有中心像素块进行步骤3.6
至步骤3.8的处理;步骤3.6:对中心像素块进行同区域搜索,搜索规则为:从中心像素块上方开始,按照顺时针顺序,逐一搜索中心像素块8邻域像素块;步骤3.7:对附属像素块进行同区域搜索,搜索规则为:从附属像素块上方开始,按照顺时针顺序,逐一搜索附属像素块8邻域像素块;步骤3.8:重复步骤3.7对所有轮次附属像素块进行附属像素块同区域搜索直到搜索8邻域没有附属像素块。4.如权利要求3所述的一种云微粒子数据区域提取方法,其特征在于所述步骤3.6具体为:对搜索到的某一邻域像素块,首先判断其属性是否为中心像素块,如果是,执行步骤3.6.1,如果不是,判断其是否为其他中心像素块的附属像素块,如果是,执行步骤3.6.2,如果不是,执行步骤3.6.3;当搜索完8邻域范围结束;步骤3.6.1:更改该邻域像素块属性为其依附的中心像素块的第1轮附属像素块,且将其所依附的中心像素块的编号和轮次数附加到该像素块编号之后作为该附属像素块的标记号,继续搜索;步骤3.6.2:设置阈值th

,如果该邻域像素块的rh

n,m,j,q
满足条件rh

n,m,j,q
<th

,在其原有属性基础上增加一项属性为当前中心像素的第1轮附属像素块,且将其当前附属的中心像素块编号和轮次数增加到该像素块已有标记号后面作为该附属像素块新增标记号;如果该邻域像素块的rh

n,m,j,q
不满足条件rh

n,m,j,q
<th

,跳过该像素块,继续搜索;步骤3.6.3:基于设置的阈值th

,如果邻域像素块的rh

n,m,j,q
满足条件rh

n,m,j,q
<th

,更改该邻域像素块属性为其依附的中心像素块的第1轮附属像素块,且将其附属的中心像素块编号和轮次数附加到该像素块编号之后作为该附属像素块的标记号,继续搜索;如果邻域像素块的rh

n,m,j,q
不满足条件rh

n,m,j,q
<th

,跳过该像素块,继续搜索。5.如权利要求3所述的一种云微粒子数据区域提取方法,其特征在于所述步骤3.7具体为:对搜索到的某一邻域像素块,首先判断该邻域像素块编号是否与此附属像素块标记号中的附属的中心像素块编号相同,如果相同,跳过这些邻域像素块,如果不相同,执行步骤3.7.1;步骤3.7.1:设置阈值th

,如果邻域像素块的rh

n,m,j,q
满足条件rh

n,m,j,q
<th

,将其设置为第e轮附属像素块,且将其附属的前一轮附属像素块标记号作为该附属像素块标记号,并且将轮次数更新为该轮的轮次数,继续搜索;如果邻域像素块的rh

n,m,j,q
不满足条件rh

n,m,j,q
<th

,跳过该像素块,继续搜索,当搜索完8邻域范围时结束搜索,其中,变量e表示附属像素块是从步骤3.5得到的中心像素块开始,第e次8邻域搜索得到的附属像素块,e∈[2,e],变量e为从中心像素块开始,搜索完所有附属像素块的轮次。6.如权利要求3所述的一种云微粒子数据区域提取方法,其特征在于所述步骤3.3:设置阈值gv,其中阈值gv∈[100,120]。7.如权利要求3所述的一种云微粒子数据区域提取方法,其特征在于所述步骤3.4:设置阈值th

,其中阈值th

∈[0.15,0.25]。8.如权利要求4所述的一种云微粒子数据区域提取方法,其特征在于所述步骤3.6.2:
设置阈值th

,其中阈值th

∈[1.5,2]。9.如权利要求3所述的一种云微粒子数据区域提取方法,其特征在于所述步骤3.5设置该中心像素块的附属像素块具体为,将该中心像素块的8邻域像素块设为该中心像素块的附属像素块。10.如权利要求9所述的一种云微粒子数据区域提取方法,其特征在于将该中心像素块8邻域像素块设为该中心像素块的附属像素块,如果出现某附属像素块同时属于多个中心像素块,将其随机分配给其中一个中心像素块。

技术总结
本发明涉及一种云微粒子数据区域提取方法,主要包括对云微粒子数据进行数据集筛选和聚类,然后对所有图像数据中的像素块属性进行初步确认,在对初步确认的中心像素块和附属像素块分别进行中心像素块同区域搜索和附属像素块同区域搜索后,以中心像素块为核心,以归属于该中心像素块的附属像素块为扩展区域,两者共同作为单个云微粒子区域,以单个云微粒子区域为对象进行云微粒子提取,通过该方法提高了云微粒子数据区域提取的完整性和准确性。了云微粒子数据区域提取的完整性和准确性。了云微粒子数据区域提取的完整性和准确性。


技术研发人员:张无瑕 刘说 张宇航 杨玲
受保护的技术使用者:成都信息工程大学
技术研发日:2022.04.11
技术公布日:2022/7/5
转载请注明原文地址: https://www.8miu.com/read-6394.html

最新回复(0)