基于三维随机球仿真的雪粒径谱分布特征参数优化方法，本发明专利技术属于遥感与地理信息系统领域。本发明专利技术是要解决现有方法成像过程中将三维降维至二维图像时不可避免地存在信息损失并且提取的参数不够精确的问题，而提供了基于三维随机球仿真的雪粒径谱分布特征参数优化方法。一、获得二值化的积雪切片图像；二、利用二值圆拟合算法对二值化的积雪切片图像进行填充；三、根据步骤二中的填充圆列表ResultList对二维模拟切片图像的二维谱分布参数进行拟合与估计；四、统计二维模拟切片图像上的模粒径。本发明专利技术主要应用于卫星遥感对地观测领域。
【技术实现步骤摘要】
本专利技术属于遥感与地理信息系统领域，主要应用于卫星遥感对地观测。
技术介绍
地表雪盖对于全球水文循环、大气过程以及全球气候变化均具有重要影响，多年来卫星遥感技术作为重要的工具应用于大尺度的雪盖参数监测和反演，而具有全天时和全天候监测能力的微波遥感技术对于雪深和雪水当量的监测更是其他手段所无法替代。微波遥感辐射传输模型是微波遥感反演地表参数的理论基础，而研究表明积雪粒子的粒径是辐射传输模型最敏感的参数之一，但这一参数的测量一直以来具有很大的主观性和不确定性。实测雪粒径的方法有多种，一个简单办法是将颗粒放到有毫米刻度的盘上，通过比较颗粒尺寸与盘上刻度线的间隔，可估计平均尺寸，这种目测估计的方法具有一定的主观性；另外的测量方法是筛分法，即选择一系列不同筛孔直径的标准筛，按照孔径从小到大依次摞起，固定在振筛机上，选择适当的模式及时长，自动振动即可实现筛分，这种测量方法易于操作，但操作过程中容易破坏雪层的固有结构和粒径，使得测量结果不够准确。为了准确模拟雪层的微波散射和辐射特征需要明确知道粒径的统计分布特征，而以上两种方法不易做到。雪粒子的显微图像采集与处理的方法以其采集过程准确，测量结果客观，不破坏雪层结构等优势为许多研究者所采用，如立体测量法，图像采集与处理法等，此类方法通过体视显微镜拍摄雪粒子的图像，并通过数字图像处理技术(二值化，圆形膨胀拟合)获得粒径尺寸的范围和分布，而由此获得的粒径所服从的概论分布对于随机分布的雪粒子的微波(散射、辐射)特性的准确计算和模拟尤为重要。但雪粒子的显微图像采集与处理的方法同时存在以下问题:首先，由于积雪切片图像获取本身的步骤较多，实验室控制要求比较严格，所以获取的图像十分有限；其次，获取图像过程中需要对雪层进行切割处理，而图像也是基于切割的断面采集的，我们实际上最终是为了获得雪层中三维粒径的分布，所以由二维图像推测三维特征的过程中本身就存在损失和较大误差，获得的参数未必能够体现雪层中雪粒径分布的三维特征。
技术实现思路
本专利技术是要解决现有方法成像过程中将三维降维至二维图像时不可避免地存在信息损失并且提取的参数自然不够精确的问题，而提供了。按以下步骤实现: 一、二值化的积雪切片图像的获得；二、利用二值圆拟合算法对二值化的积雪切片图像进行填充，得到不同半径的红色圆盘的统计列表；三、根据步骤二中的不同半径的红色圆盘的统计列表对二维模拟切片图像的二维谱分布参数进行拟合与估计；四、统计二维模拟切片图像上的模粒径，并进一步和原始的二维谱分布参数进行误差分析统计输出三维优化雪粒径。专利技术效果:(I)尽量从较少的采样图像中获取更多的信息，能够弥补三维到二维时的信息损失，使得提取的雪粒径分布参数更加准确真实，从而避免了需要更多的二维切片图像所带来的累计误差，总体精度可提高90%以上；(2)利用初始的雪粒径参数进行计算机模拟得到三维随机球拟合图像，模拟切片过程，多次随机产生切片图像，避免了由于切片图像制作过程切割位置不同所带来的误差，使得获得的优化参数获 取更加准确，为积雪微波辐射传输模型的研究以及地表积雪参数的准确提取提供了技术保障。附图说明图1是本专利技术的整体流程图；图2是具体实施方式一中的积雪样本的二值化切片图像；图3是具体实施方式二中的输入的图像进行圆填充流程图；图4是具体实施方式三中的二值图像圆拟合算法获得的粒径拟合结果；图5是具体实施方式三中的二维谱分布参数的估计流程图；图6是具体实施方式三中的三维切片未必能够恰好切至雪粒子的直径处的演示图；图7是具体实施方式三中的二维雪粒径分布函数图；其中，am = 0.38, P = 9, Q =2，ftot = 0.4 ；图8是具体实施方式四中的三维随机球仿真(3DRBS)优化算法的实现流程图；图9是具体实施方式四中的单位体积内(Icm3)三维模拟雪粒径图；图10是具体实施方式四中的模拟三维雪粒子的在0.5cm处的切片图像；图11是具体实施方式四中的模拟三维雪粒子的在0.8cm处的切片图像；图12是具体实施方式四中的迭代的分布曲线变化过程，虚线是初始的分布曲线，实线为迭代曲线。具体实施例方式具体实施方式一:本实施方式的按以下步骤实现:一、二值化的积雪切片图像的获得；二、利用二值圆拟合算法对二值化的积雪切片图像进行填充，得到不同半径的红色圆盘的统计列表；三、根据步骤二中的不同半径的红色圆盘的统计列表对二维模拟切片图像的二维谱分布参数进行拟合与估计；四、统计二维模拟切片图像上的模粒径，并进一步和原始的二维谱分布参数进行误差分析统计输出三维优化雪粒径。本实施方式效果:(I)尽量从较少的采样图像中获取更多的信息，能够弥补三维到二维时的信息损失，使得提取的雪粒径分布参数更加准确真实，从而避免了需要更多的二维切片图像所带来的累计误差，总体精度可提高90%以上；(2)利用初始的雪粒径参数进行计算机模拟得到三维随机球拟合图像，模拟切片过程，多次随机产生切片图像，避免了由于切片图像制作过程切割位置不同所带来的误差，使得获得的优化参数获取更加准确，为积雪微波辐射传输模型的研究以及地表积雪参数的准确提取提供了技术保障。具体实施方式二:本实施方式与具体实施方式一不同的是:步骤二中利用二值圆拟合算法对二值化的积雪切片图像进行填充具体步骤如下:步骤A:设定圆直径的阈值参数aT，对于雪粒径处理，该范围指定在0.2mm_2.5mm ；步骤B:在 二值图像的白色区域(雪粒子)放入尽可能大直径的圆，获得图像中的填充圆列表ResultList，进一步得到雪粒径的大小及位置数据；B1:对于图像中的每一个白像素(对应数组值为I),初始化圆的半径r = O ;B2:判断以该像元为圆心、r为半径的圆范围内是否包含值为O的像素或者图像边界，如果包含，r = r+l，继续判断；否则，转步骤B3 ；B3:判断r是否大于0，如果是，将该圆的数据(包括:中心位置和半径)存储在填充圆列表List，转入步骤BI，计算其他的白像素；否则，直接转步骤BI ；B4:对List中的所有圆按照半径从大到小排序；B5:初始化填充圆列表ResultList ；B6:判断List中是否还包含数据，如果是，执行步骤B8 ;否则，跳入B8输出ResultList ；B7:删除List中最大元素x ;找到List中与x相交叠的所有圆，放入列表TList之中；将TList的所有元素从List中删除；将X加入到ResultList ；B8:输出填充圆列表ResultList ；步骤C:在二值图像的白色区域(雪粒子)继续利用尽可能大的圆填充剩余部分的白色区域；步骤D:检查是否剩余的白色区域能够利用阈值范围内的圆继续填充，如果能，转到步骤C，否则，结束。其它步骤及参数与具体实施方式一相同。具体实施方式三:本实施方式与具体实施方式一或二不同的是:步骤三中对二维模拟切片图像的二维谱分布参数进行拟合与估计具体步骤如下:雪粒径按照ga_a分布进行统计:n (a) =Kpp exp (_K2aQ) (I)其中，a表示雪粒径；雪粒子的占空比ft(lt可以用下式表示为:本文档来自技高网...
【技术保护点】
基于三维随机球仿真的雪粒径谱分布特征参数优化方法，其特征在于基于三维随机球仿真的雪粒径谱分布特征参数优化方法按以下步骤实现：一、二值化的积雪切片图像的获得；二、利用二值圆拟合算法对二值化的积雪切片图像进行填充，得到填充圆列表ResultList；三、根据步骤二中的填充圆列表ResultList对二维模拟切片图像的二维谱分布参数进行拟合与估计；四、统计二维模拟切片图像上的模粒径，并进一步和原始的二维谱分布参数进行误差分析统计输出三维优化雪粒径。
【技术特征摘要】
1.基于三维随机球仿真的雪粒径谱分布特征参数优化方法，其特征在于基于三维随机球仿真的雪粒径谱分布特征参数优化方法按以下步骤实现: 一、二值化的积雪切片图像的获得； 二、利用二值圆拟合算法对二值化的积雪切片图像进行填充，得到填充圆列表ResultList ； 三、根据步骤二中的填充圆列表ResultList对二维模拟切片图像的二维谱分布参数进行拟合与估计； 四、统计二维模拟切片图像上的模粒径，并进一步和原始的二维谱分布参数进行误差分析统计输出三维优化雪粒径。2.根据权利要求1所述的基于三维随机球仿真的雪粒径谱分布特征参数优化方法，其特征在于步骤二中利用二值圆拟合算法对二值化的积雪切片图像进行填充，具体步骤如下: 步骤A:设定圆直径的阈值参数aT，对于雪粒径处理，该范围指定在0.2mm-2.5mm ； 步骤B:在二值图像的白色区域(雪粒子)放入尽可能大直径的圆，获得图像中填充圆列表ResultList，进一步得到雪粒径的大小及位置数据； B1:对于图像中的每一个白像素(对应数组值为I)，初始化圆的半径r = O ; B2:判断以该像元为圆心、r为半径的圆范围内是否包含值为O的像素或者图像边界，如果包含，r = r+Ι，继续判断；否则，转步骤B3 ； B3:判断!■是否大于O，如果是，将该圆的数据(包括:中心位置和半径)存储在填充圆列表List，转入步骤B I，计算其他的白像素；否则，直接转步骤BI ； B4:对List中的所有圆按照半径从大到小排序； B5:初始化填充圆列表ResultList ； B6:判断List中是否还包含数据，如果是，执行步骤B8 ;否则，跳入B8输出ResultList ； B7:删除List中最大元素X ;找到List中与x相交叠的所有圆,放入列表TList之中；将TList的所有元素从List中删除；将X加入到ResultList ； B8:输出填充圆列表ResultList ； 步骤C:在二值图像的白色区域(雪粒子)继续利用尽可能大的圆填充剩余部分的白色区域； 步骤D:检查是否剩余的白色区域能够利用阈值范围内的圆继续填充，如果能，转到步骤C，否则，结束。3.根...
【专利技术属性】
技术研发人员：李晓峰，赵凯，武黎黎，
申请(专利权)人：中国科学院东北地理与农业生态研究所，
类型：发明
国别省市：
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。