一种盐碱化土壤含盐量的在线测量方法技术

一种盐碱化土壤含盐量的在线测量方法，本发明专利技术涉及一种盐碱化土壤含盐量的在线测量方法，它为了解决现有盐碱化土壤含盐量的测量方法采样费时费力，含盐量的在线测量精度较低的问题。在线测量方法：一、通过土壤开裂实验获取开裂土壤样本照片，然后提取纹理特征，再测量土壤中盐分含量，建立土壤含盐量反演模型；二、在线测量野外土壤含盐量，计算野外土壤归一化灰度共生矩阵和相应拟合关系最好的纹理特征值，得到土壤含盐量的初步在线测量结果；三、通过土壤样本含盐量的最大值和最小值，修正土壤含盐量测量结果。本发明专利技术基于盐碱土含盐量对土壤干缩开裂程度的影响，建立开裂样本的纹理特征量与样本含盐量的定量模型，提高测量精度。

【技术实现步骤摘要】
一种盐碱化土壤含盐量的在线测量方法
本专利技术涉及一种盐碱化土壤含盐量的在线测量方法，尤其涉及室内盐碱化土壤含盐量反演模型的建立及通过在线测量盐碱土纹理特征实现土壤含盐量的在线测量方法。
技术介绍
土壤盐碱化是成为世界性的土地退化过程，它严重破坏土壤性状并且在很大程度上降低了农作物的产量，同时对生态环境也有较大的危害。因此，快速准确的测量土壤含盐量水平对于确定土壤盐碱化程度、指导盐碱化土壤的改良具有重要的现实意义。目前，传统的含盐量测量方法主要针对野外采样、室内测量土壤内各盐离子的总含量来确定的，这种土壤含盐量的测量方法费时费力且不能实时反应含盐量的野外真实状况，遥感手段虽然具有大面积实时的特点，但是遥感数据本身的空间分辨率和光谱分辨率有限，使其局限于大面积的含盐量程度的确定和制图应用，而小尺度的单点含盐量测量精度却很差。松嫩平原的苏打型盐碱土是一种典型的盐碱化土壤，其土壤粘粒含量较高具有很强的膨胀收缩特性。因此在降水过后，土壤表面随水分的蒸发而出现的干缩开裂现象非常普遍。截至目前，已有很多研究关注于土壤裂纹特征的识别和测量以及土壤理化参数、环境条件等对开裂程度的影响等方面。盐分含量作为盐碱化土壤主要的化学特性，对土壤的开裂过程和裂纹特征也有很大程度的影响。
技术实现思路
本专利技术的目的是为了解决现有盐碱化土壤含盐量的测量方法采样费时费力，含盐量的在线测量精度较低的问题，而提供一种盐碱化土壤含盐量的在线测量方法。本专利技术盐碱化土壤含盐量的在线测量方法按下列步骤实现：一、建立土壤含盐量反演模型1、土壤开裂实验选择不同盐分含量的多组土壤样本，对所有的土壤样本进行风干、研磨处理并过2mm孔径的筛子，然后将每组过筛后的土壤样本分为土壤样本A和土壤样本B，将其中的土壤样本A进行干缩开裂实验，所述的干缩开裂实验是先将所有的土壤样本A分别配置成含水量相同的饱和泥浆，把每组配置好的饱和泥浆分别倒入50厘米×50厘米×3厘米大小的木质样本盒中，抹平泥浆后放置在实验室中进行自然干燥以产生干缩裂纹，直至样本盒中样品的质量不再减少时，得到干燥开裂的土壤样本A；将数码相机安装在固定的实验平台上，数码相机距离地面高度为1m，以数码相机镜头在地面上的投影为中心确定50厘米×50厘米的正方形区域以保证所有干燥开裂的土壤样本A拍照过程的几何畸变影响相同，将50厘米×50厘米的黑白网格定标板覆盖在干燥开裂的土壤样本A表面，然后使用数码相机对干燥开裂的土壤样本A进行拍照，获取开裂土壤样本照片；2、提取纹理特征对每张开裂土壤样本照片进行几何畸变校正，对几何畸变校正后的样本图像进行剪切以获取样本的彩色裂纹区域图像，然后将每张彩色裂纹区域图像转化为灰度图，再根据灰度直方图选取阈值进而进行二值化处理，取反操作后则将每个开裂土壤样本的裂纹图像分割为未开裂的块区(用黑色象元表示)和收缩区域(包括裂纹部分和样本盒内四周的收缩部分，用白色象元表示)；计算每个开裂土壤样本二值化图像的灰度共生矩阵，将距离步长设置为40，4个方向分别为0°、45°、90°和135°，然后将四个方向得到的归一化灰度共生矩阵进行算数平均从而以去除方向的影响，根据每个开裂土壤样本灰度共生矩阵计算各统计纹理特征值，各统计纹理特征值的计算公式如下：1)对比度：2)角二阶矩：3)熵：4)一致性：5)相关性：6)集群荫：7)最大概率：Maxprobability＝max{p(i,j)}，(7)8)和平均：px+y(i)表示x行和y列同时出现灰度级i的概率值；9)和熵：10)和方差：11)相关信息特征1：12)相关信息特征2：Inforofcorrelation2＝{1-exp[-2*(HXY2-HXY)]}1/2，(12)其中：p(i,j)是归一化灰度共生矩阵(i,j)位置元素值，即图像灰度级i，j同时出现的概率值，Ng是共生矩阵灰度级数，i、j均为正整数，μx，μy为px(i)，py(j)均值，σx，σy为px(i)，py(j)方差，n为灰度级；3、测量土壤盐分将每个土壤样本B配置成土壤浸提液，分别测量土壤样本B中的Ca2+、Mg2+、Cl-、CO2-和HCO3-含量，计算每个土壤样本B中Ca2+、Mg2+、Cl-、CO2-和HCO3-含量的总和，即为土壤样本B的含盐量；4、建立反演模型利用土壤样本B的含盐量及计算得到的各统计纹理特征值建立线性模型，其中线性模型中的横坐标(自变量)为纹理特征量，纵坐标(因变量)为含盐量，选择与土壤样本B含盐量拟合关系最好的纹理特征量，建立土壤含盐量野外在线测量的标准线性反演模型；二、在线测量野外土壤含盐量针对待测野外盐碱化土壤的在线测量，将数码相机放置在移动平台上，调整平台高度使数码相机本文档来自技高网...

【技术保护点】
一种盐碱化土壤含盐量的在线测量方法，其特征在于盐碱化土壤含盐量的在线测量方法是按下列步骤实现：一、建立土壤含盐量反演模型1、土壤开裂实验选择不同盐分含量的多组土壤样本，对所有的土壤样本进行风干、研磨处理并过2mm孔径的筛子，然后将每组过筛后的土壤样本分为土壤样本A和土壤样本B，将其中的土壤样本A进行干缩开裂实验，所述的干缩开裂实验是先将所有的土壤样本A分别配置成含水量相同的饱和泥浆，把每组配置好的饱和泥浆分别倒入50厘米×50厘米×3厘米大小的木质样本盒中，抹平泥浆后放置在实验室中进行自然干燥以产生干缩裂纹，直至样本盒中样品的质量不再减少时，得到干燥开裂的土壤样本A；将数码相机安装在固定的实验平台上，数码相机距离地面高度为1m，以数码相机镜头在地面上的投影为中心确定50厘米×50厘米的正方形区域以保证所有干燥开裂的土壤样本A拍照过程的几何畸变影响相同，将50厘米×50厘米的黑白网格定标板覆盖在干燥开裂的土壤样本A表面，然后使用数码相机对干燥开裂的土壤样本A进行拍照，获取开裂土壤样本照片；2、提取纹理特征对每张开裂土壤样本照片进行几何畸变校正，对几何畸变校正后的样本图像进行剪切以获取样本的彩色裂纹区域图像，然后将每张彩色裂纹区域图像转化为灰度图，再根据灰度直方图选取阈值进而进行二值化处理，取反操作后则将每个开裂土壤样本的裂纹图像分割为未开裂的块区和收缩区域；计算每个开裂土壤样本二值化图像的灰度共生矩阵，将距离步长设置为40，4个方向分别为0°、45°、90°和135°，然后将四个方向得到的归一化灰度共生矩阵进行算数平均从而以去除方向的影响，根据每个开裂土壤样本灰度共生矩阵计算各统计纹理特征值，各统计纹理特征值的计算公式如下：1)对比度：contra=ΣnNg-1n2{Σi=1NgΣj=1Ngp(i,j)},|i-j|=n---(1)]]>2)角二阶矩：ASM=ΣiΣj{p(i,j)}2,---(2)]]>3)熵：Entropy=-Σi=1NgΣj=1Ngp(i,j)log(p(i,j)),---(3)]]>4)一致性：Homo=Σi=1NgΣj=1Ng11+(i-j)2p(i,j),---(4)]]>5)相关性：Correlation=Σi=1NgΣj=1Ng(ij)p(i,j)-μxμyσxσy,---(5)]]>6)集群荫：Clustershade=Σi=1NgΣj=1Ng((i-μi)+(j-μj))3p(i,j),---(6)]]>7)最大概率：Maxprobability＝max{p(i,j)}，   (7)8)和平均：SumAverage=Σi=22Ngipx+y(i),---(8)]]>px+y(i)表示x行和y列同时出现灰度级i的概率值；9)和熵：SumEntropy=-Σi=22Ngpx+y(i)log{px+y(i)},---(9)]]>10)和方差：SumVariance=Σi=22Ng(i-SumEntropy)2px+y(i),---(10)]]>11)相关信息特征1：Inforofcorrelation1=HXY-HXY1max(HX,HY),---(11)]]>12)相关信息特征2：Inforofcorrelation2＝{1‑exp[‑2*(HXY2‑HXY)]}1/2，   (12)其中：p(i,j)是归一化灰度共生矩阵(i,j)位置元素值，即图像灰度级i，j同时出现的概率值，Ng是共生矩阵灰度级数，i、j均为正整数，px(i)=Σj=1Ngp(i,j),---(13)]]>py(j)=Σi=1Ngp(i,j),---(14)]]>px-y(k)=Σi=1NgΣj=1Ngp(i,j),i+j=k,k=2,3,...2Ng,---(15)]]>HXY=-Σi=1NgΣj=1Ngp(i,j)log{p(i,j)},---(16)]]>HXY1=-Σi=1NgΣj=1Ngp(i,j)log{px(i)py(j)},---(17)]]>HXY2=-Σi=1NgΣj=1Ngpx(i)py(j)log{px(i)py(j)},---(18)]]>μx，μy为px(i)，py(j)均值，σx，σy为px(i)，py(j)方差，n为灰度级；3、测量土壤盐分将每个土壤样本B配置成土壤浸提液，分别测量土壤样本B中的Ca2+、Mg2+、Cl‑、CO2‑和HCO3‑...

【技术特征摘要】
1.一种盐碱化土壤含盐量的在线测量方法，其特征在于盐碱化土壤含盐量的在线测量方法是按下列步骤实现：一、建立土壤含盐量反演模型1、土壤开裂实验选择不同盐分含量的多组土壤样本，对所有的土壤样本进行风干、研磨处理并过2mm孔径的筛子，然后将每组过筛后的土壤样本分为土壤样本A和土壤样本B，将其中的土壤样本A进行干缩开裂实验，所述的干缩开裂实验是先将所有的土壤样本A分别配置成含水量相同的饱和泥浆，把每组配置好的饱和泥浆分别倒入50厘米×50厘米×3厘米大小的木质样本盒中，抹平泥浆后放置在实验室中进行自然干燥以产生干缩裂纹，直至样本盒中样品的质量不再减少时，得到干燥开裂的土壤样本A；将数码相机安装在固定的实验平台上，数码相机距离地面高度为1m，以数码相机镜头在地面上的投影为中心确定50厘米×50厘米的正方形区域以保证所有干燥开裂的土壤样本A拍照过程的几何畸变影响相同，将50厘米×50厘米的黑白网格定标板覆盖在干燥开裂的土壤样本A表面，然后使用数码相机对干燥开裂的土壤样本A进行拍照，获取开裂土壤样本照片；2、提取纹理特征对每张开裂土壤样本照片进行几何畸变校正，对几何畸变校正后的样本图像进行剪切以获取样本的彩色裂纹区域图像，然后将每张彩色裂纹区域图像转化为灰度图，再根据灰度直方图选取阈值进而进行二值化处理，取反操作后则将每个开裂土壤样本的裂纹图像分割为未开裂的块区和收缩区域；计算每个开裂土壤样本二值化图像的灰度共生矩阵，将距离步长设置为40，4个方向分别为0°、45°、90°和135°，然后将四个方向得到的归一化灰度共生矩阵进行算数平均从而以去除方向的影响，根据每个开裂土壤样本灰度共生矩阵计算各统计纹理特征值，各统计纹理特征值的计算公式如下：1)对比度：2)角二阶矩：3)熵：4)一致性：5)相关性：6)集群荫：7)最大概率：Maxprobability＝max{p(i,j)}，(7)8)和平均：px+y(i)表示x行和y列同时出现灰度级i的概率值；9)和熵：10)和方差：11)相关信息特征1：12)相关信息特征2：Inforofcorrelation2＝{1-exp[-2*(HXY2-HXY)]}1/2，(12)其中：p(i,j)是归一化灰度共生矩阵(i,j)位置元素值，即图像灰度级i，j同时出现的概率值，N...

【专利技术属性】
技术研发人员：任建华，赵凯，李晓洁，
申请(专利权)人：中国科学院东北地理与农业生态研究所，
类型：发明
国别省市：黑龙江;23