[发明专利]土壤裂纹参数在线测量系统及采用该系统实现土壤裂纹参数提取方法

权利要求书

1.土壤裂纹参数在线测量系统，其特征在于，它包括三脚架(1)、数码相机(2)、笔记本电脑(3)、环形灯(4)、数字光度计(5)、取样框(6)、黑白比色板(7)和黑白格网定标板(8)，

数码相机(2)固定在三脚架(1)上，所述数码相机(2)的摄像区域为其正下方，环形灯(4)用于给数码相机(2)的摄像区域提供照明，黑白比色板(7)用于对数码相机进行白平衡处理和灰度值拉伸处理，数字光度计(5)用于测量每次拍摄时的背景环境光照强度，取样框(6)放置在数码相机(2)的摄像区域内，用于拍摄照片时确定标准的裂纹区域尺寸，数码相机(2)通过数据连接线将拍摄的照片传送给笔记本电脑(3)。

2.采用权利要求1所述的土壤裂纹参数在线测量系统实现土壤裂纹参数提取方法，其特征在于，所述方法为：

步骤一、将取样框(6)放置在土壤上，使得待检测土壤裂纹区域位于所述取样框(6)内，通过三脚架(1)将数码相机(2)固定在取样框(6)的正上方，使得数码相机(2)的摄像区域覆盖取样框(6)，并且取样框(6)的中心与数码相机(2)镜头中心正下方的投影点重合，采用固定环形灯(4)对取样框(6)所在区域进行照明；

步骤二、利用数字光度计(5)测量被测土壤样本上方的环境光照强度，根据所述光照强度设置数码相机(2)的ISO、光圈尺寸和曝光时间的参数，保证不同环境下光照对相片的影响相同，利用黑白比色板(7)对数码相机(2)进行白平衡处理和灰度值拉伸处理，将对比度和色差对照片的影响减到最小；

步骤三、控制数码相机(2)连续采集图像信息，每次土壤裂纹拍照后，在取样框(6)内放置黑白格网定标板(8)后再次拍照，用于对土壤裂纹照片进行几何畸变校正；

步骤四、数码相机(2)对土壤裂纹照片进行几何定标，对标定后的照片按取样框(6)尺寸进行裁剪，用于提取取样框(6)内土壤裂纹参数；

笔记本电脑(3)内部嵌入有土壤裂纹参数提取软件，该土壤裂纹参数提取软件对土壤裂纹照片进行图像处理操作以获得土壤裂纹参数的提取为：

利用图像处理方法，将裁剪后的裂纹照片进行灰度图转换和二值处理，计算代表裂纹区域的像素个数，提取裁剪后的取样框(6)内总像素个数并按取样框(6)尺寸比例计算土壤裂纹面积；

根据取样框(6)尺寸计算取样框面积，结合提取的土壤裂纹面积计算土壤面裂纹率；

利用图像处理方法设置纹理提取灰度级、方向和步长，计算图像的灰度共生矩阵，进而根据计算的灰度共生矩阵提取取样框(6)内开裂土壤的各个统计纹理特征量参数，公式如下：

1)对比度：

contra=ΣnNg-1n2{Σi=1NgΣj=1Ngp(i,j)},|i-j|=n,]]>

2)角二阶矩：

ASM=ΣiΣj{p(i,j)}2,]]>

3)熵：

Entropy=-Σi=1NgΣj=1Ngp(i,j)log(p(i,j)),]]>

4)一致性：

Homo=Σi=1NgΣj=1Ng11+(i-j)2p(i,j),]]>

5)相关性：

Correlation=Σi=1NgΣj=1Ng(ij)p(i,j)-μxμyσxσy,]]>

6)集群荫：

Clustershade=Σi=1NgΣj=1Ng((i-μi)+(j-μj))3p(i,j),]]>

7)集群突出性：

Clustershade=Σi=1NgΣj=1Ng((i-μi)+(j-μj))4p(i,j),]]>

8)最大概率：

Maxprobability＝max{p(i,j)}，

9)和平均：

SumAverage=Σi=22Ngipx+y(i),]]>

10)和熵：

SumEntropy=-Σi=22Ngpx+y(i)log{px+y(i)},]]>

11)和方差：

SumVariance=Σi=22Ng(i-SumEntropy)2px+y(i),]]>

12)相关信息特征1：

Inforofcorrelation1=HXY-HXY1max(HX,HY),]]>

13)相关信息特征2：

Inforofcorrelation2＝{1-exp[-2*(HXY2-HXY)]}^1/2，

其中：

p(i,j)是归一化灰度共生矩阵(i,j)位置元素值，Ng是共生矩阵灰度级数，i、j均为正整数，

px(i)=Σj=1Ngp(i,j),]]>

py(j)=Σi=1Ngp(i,j),]]>

px-y(k)=Σi=1NgΣj=1Ngp(i,j),i+j=k,k=2,3,....2Ng,]]>

px-y(k)=Σi=1NgΣj=1Ngp(i,j),|i-j=k|,k=0,1,2,...Ng-1,]]>

HXY=-Σi=1NgΣj=1Ngp(i,j)log{p(i,j)},]]>

HXY1=-Σi=1NgΣj=1Ngp(i,j)log{px(i)py(j)},]]>

HXY2=-Σi=1NgΣj=1Ngpx(i)py(j)log{px(i)py(j)},]]>

μ_x，μ_y为p_x(i)，p_y(j)均值，σ_x，σ_y为p_x(i)，p_y(j)方差，n为灰度级，

选择各纹理特征量4个方向0°、45°、90°和135°的均值作为最终的各纹理特征量。