大学本科毕业设计（论文）
end %m大于n-h退出循环； end %h大于n/2退出循环；
avg_k=avg_k+k; %求出变差平均值放入avg_k； end
从像素矩阵中提取i行A行阵，A1，A2，A3……Ai 读图I1 NO A阵列数是否为偶数 YES 向0取偶整数 整数变程h≤n/2 YES 求变差函数值行阵 求平均值变差函数值矩阵 绘 图
图 3-1 水平方向上流程图
NO 3.2 变差函数在垂直方向的流程图及程序设计
取像素矩阵列方向上的矩阵设计流程图如图3-2所示，程序与水平方向相似，
15
大学本科毕业设计（论文）
这里不再赘述。
从像素矩阵中提取i列A列阵，A1，A2，A3……Ai NO A阵行数是否为偶数 YES 向0取偶整数 NO 整数变程h≤n/2 YES 求变差函数值行阵 求平均值变差函数值矩阵 绘 图 图 3-2 垂直方向上流程图
读图I1 3.3 变差函数在任意方向的流程图及程序设计
当取像素矩阵I1任意方向上的矩阵时设计程序就会比较复杂。此时我们选取3个变量分别为start_line,n_rank,AA。选取目标矩阵A窗口的起始行为I1矩阵的start_line行，n_rank为选取的列数，AA为选取的角度。
在选取A阵时需要注意3点：
16
大学本科毕业设计（论文）
(1)目标矩阵A窗口的长度的确定。n_rank为给出I1像素阵m*n中n方向上（即x轴上）的列数（如图3-3所示，I1为任意给出的一方阵），也就是选取的目标矩阵（旋转某角度后）A窗口的列数。因为以任意角度划线取出一个矩阵（假设放到一行阵中），列数均和a线、x轴的两行矩阵的列数相同，虽然a线比x轴距离长，但是在矩阵中对应的列数是相同的，比如a轴为[1 2 3 4 5 6 7 8 9]，x轴为[1 1 1 1 1 1 1 1]，列数相同。而对于b线上的矩阵为[1 1 1 2 2 3 4 4 5]，没有画到线的列，行数向0取整。因此A窗口阵的列数与任何方向上的列数都相同。
(2)确定x轴长度求出y轴长度。n_rank为给出I1像素阵x轴方向上的列数（如图3-3所示），根据公式n_line=n_rank*tand(AA)，求出y轴上的行数，加上起始行start_line，就能确定出在像素矩阵中I1的坐标（位置）。
(3)角度AA的方向是以x轴为边向下（顺时针）旋转AA角度（如图3-3中a线为45?），并且为锐角。
0 1 1 1 1 1 1 1 1 1 2 2 2 2 2 2 2 2 2 3 3 3 3 3 3 3 3 3 4 4 4 4 4 4 4 4 4 5 5 5 5 5 5 5 5 5 6 6 6 6 6 6 6 6 6 7 7 7 7 7 7 7 7 7 8 8 8 8 8 8 8 8 8 9 9 9 9 9 9 9 9 9 Y 轴 X轴 b a 图3-3 I1像素阵
流程图如图3-4。
17
大学本科毕业设计（论文）
读图 I1
确定起始行为第start_line行第1列 给定n_rank列（X轴方向），AA角 用公式n_line=n_rank*tan(AA)，确定出像素矩阵I1 中对应的像素点，并且存入A阵中，取i行A阵， 记为A1，A2，A3……Ai NO A阵行数是否为偶数 YES 向0取偶整数 NO 整数变程h≤n/2 YES 求变差函数值行阵 求平均值变差函数值矩阵 绘 图 图 3-4 任意方向上流程图
18
大学本科毕业设计（论文）
此处仅就A阵窗口的选取程序作以解释，变差函数程序不再赘述： %n_line=n_rank*tand(AA);
n_rank=ceil(n_rank/2)*2; %对列数向正无穷圆整；
avg_k=zeros(1,n_rank/2); %初始化 avg_k矩阵； for z=start_line:start_line+9; %A阵总共取得行数为10行(可自选行数)； for j=1:n_rank; %I1的列j从1取到n_rank； i=z+fix(j*tand(AA))-1; %求出向0取整；
A(z,j)=I1(i,j,1); %从I1中取出的像素存入A阵中； end %当i＞n_rank时退出循环； end %当z＞10时退出循环；
上面叙述了选取的A阵窗口是从像素矩阵的左上角划线到右下角的方向（从X轴起顺时针方向旋转AA角度），而在本实验中，为了全方位的检测出纹理图像的周期性，我们还选择了从左下角向右上角（从X轴起逆时针旋转AA角度）划线选取A阵窗口。程序设计中只需将i=z+fix(j*tand(AA))-1改为i=z+fix(j*tand(360-AA))-1，流程图相同。
3.4 本章小结
本章主要讲到在各个方向上变差函数的流程图设计以及程序设计的相关内容。
（1）在取像素矩阵行和列方向上的矩阵时主要介绍变差函数的计算方法以及程序设计，以及变量的选取[13]。
（2）在取像素矩阵任意方向上的矩阵时主要弄清3.3节中的3个注意点，然后运用3.1中的变差程序进行变差计算。
19