matpiv大作业.docVIP

MATPIV 1.6.1 阅读报告 第一部分 MATPIV 1.6.1的部分简介 在这一部分，我们将通过对例子的讲解来认识MATPIV，这些图像可以在MATPIV安装程序的子文件demo3下找到。它们是从一个水波实验中筛选出的。图像在figure2.1中被展示出来。 简单的方法 我们可以简单的执行第一次PIV计算通过改变demo3的地址录和执行下列命令： [x,y,u,v]=matpiv(’mpim1b.bmp’,’mpim1c.bmp’,64,0.0012,0.5,’single’); 这条命令将会使图像在64*64像素大小的转换窗口中单向转换并产生50%的重叠，单项传递指的是通过一次交互遍历所有的图像进行计算，为了达到更高的精度，最好使用多次交互。这些我们后来会谈。 结果包含四个以像素为计量单位的变量：x，y，u，v，结果可以通过下列命令被可视化： quiver (x ,y, u, v) 这样就都好了，但正常情况下我们更想知道以像cm和s这样工程单位表示出的斜率和位置。要做到这一点，我们需要知道在图像里，以物理坐标衡量每个像素有多大。另外，如果能掩盖切掉水波的特定区域更好。为了节省计算时间，没有必要计算图2.1中波峰以上真空区域的面积。 第二部分 definewoco.m文件的解读 definewoco文件是计算像素于坐标系中大小的文件。这个功能需要图像中有明显的坐标系点的痕迹，这些点或是需要被清楚的定位，或是有确定的距离间隔，点需要配置黑的或白的背景。 definewoco(woco.bmp,o);通过一个定义好的峰来假定点的位置，definewoco(woco.bmp,+);假定图像中有一个距离确定的网格，因此图像可以和假定的十字交叉线相互关联。definewoco(woco.bmp,x);假定图像中包含x形式的网格点。因此图像可以和假定的x方向的点相互关联。definewoco(woco.bmp,.);假定点的峰是没有被定义好的，举个例子，如果所求峰平坦并且相对比较大，输入的图像可以通过与高斯正态分布曲线相互关联来估算原图中每个点的中心位置。这个选项增加了允许使用者了解坐标系中点精确大小的功能。它在点非常多的时候很有帮助，比如说二十像素大小的区域。 用户然后要用鼠标左键标记出局部区域周围的每个坐标点。然后用户将被要求输入每个点的坐标。最后一步，为了将作图因素带来的误差放入计算考虑之内，用户需要在线性和非线性的作图功能做出选择，一般情况下，线性作图会有更好的效果。 最后结果会保存在当前目录下的文件中，文件命名为worldcoX.mat，X是由用户任意定义的数字，这样做是为了避免在同一个工作目录下会有两个或两个以上的坐标系同时存在的情况发生。 第三部分 程序注释： format long (长字节形式) if ischar (filename) （如果是文件名字符串则为真） A=imread(filename); （则定义变量A为指令“从文件中读取图像数据”） Else （如果是假） A=filename; （则定义变量A为文件名） End if isrgb(A) （如果图像中有三原色） A=double(rgb2gray(A)); （则将彩色图像转变为灰色图像，并转换为双精度的数据形式后赋值给A） Else （如果图像中没有三原色） A=double(A); （则将原数据转换为双精度的数据形式后赋给A） End my_ver=version; （定义my_ver为版本值） my_ver=str2num(my_ver(1:3)); （将my_ver字符串矩阵中的一至三行转换为数字矩阵，并重新赋给my_ver） if my_ver=6.5 （如果该矩阵的值大于或等于6.5为真） pixval on End if strcmp(typ,+)==1 （如果字符串typ和+采用比较字符串函数处理后所得值为真值） load articross.mat （加载articross.mat 文件） disp(....calculating....this may take a few seconds.) （显示“计算中，等一会。。。） b=A./max(A(:)); （令变量B等于方阵阵A对应行列式的值和方阵A的最大非零子阵对应的行列式的值相比） A=xcorrf2(b-mean(b(:)),kr-mean(kr(:)))./(std(kr(:))*std(b(:))*size(kr,1)*ay);(将两个矩阵进行平面傅里叶变换所得结果赋给A） [ax,ay]=size(A)； （将A的尺寸值赋给ax ay矩阵） [bx,by]=size(b); （将B的尺寸值赋