粒子图像测速技术.docVIP

热灾害实验技术粒子图像测速技术 实验目的：测量流场的速度分布 实验装置：摄像头，片光源，计算机 实验原理与方法: 互相关算法: 畛(亦At) 畛(亦At) Q nftiao) 根据DPIV互相关算法理论，过一采样窗口，同一特定位置，可顺序 获取两幅数字图像。粒子的位移可从一幅图像到另一幅相对应的图像 经互相关计算获得。两幅顺序获取图像中相同位置的两个查问窗口函 数f(m, n)和g(m, n),从信号系统的观点出发，g(m, n)可以看作 是f(m, n)经线性转换后叠加以噪声而成，如图所示。 (叫 n)G(uf v)f (眄 n) s (叫 n)G(uf v) S(u, v) bf(u,v) d(a,n) D(u,v) 1、 预先在流场屮撒入一定浓度和大小的示踪粒子，用脉冲激光片光 源照明二维流场，在垂直片光的方向上摄像，获得两次或者多次曝光 的粒子图像。 2、 对采取的图像使用Photoshop进行处理使之变成长宽为256x256 的灰度图像，使用matbb和互相关算法编程，对图像进行处理，得 到二维流场屮速度矢量分布图。 实验结果： 隔两帧得到的灰度图： 得到速度矢量图如图所示: 得到速度矢量图如图所示: 源程序： function [XNumz YNumz AutoFlag]=g(al,a2z OutFile) fidl = fopen(ml.jpg, 1r *); fid2 = fopen(a2 ? jpg 1, 1r *); fid3 = fopen(OutFile ? jpgz 1w); XL=256; YL=25 6; frewind(fidl); Imagel=fread(fidl,[XL,YL]); frewind(fid2); Image2=fread(fid2z [XL,YL]); WinSizeX=64; WinSizeY=64; DeltaX=l6; DeltaY=l6; XNum = (XL-WinSizeX)/DeltaX; YNum = (YL-WinSizeY)/DeltaY; m=3.4968e+006 +4.7425e-01Oi; frewind(fid3); for jl=l:YNum+1; for il=l:XNum+1; CorrWinl=Image1((il-1)*DeltaX+1:(il-1)*DeltaX+WinSizeXz (jl-1)*Del taY+1:(jl-1)*DeltaY+WinSizeY); CorrWin2=Image2((il-1)*DeltaX+l:(il-1)*DeltaX+WinSizeXz (jl-1)*Del taY+1:(jl-1)*DeltaY+WinSizeY); Fl=fft2(CorrWinl); F2=fft2(CorrWin2); S=F2 ? *conj (Fl); Pks = ifft2 (S); Pks=f ftshift(Pks); Pks=Pks/a; if (il==2) SURF(abs(Pks)); pause end [tempiz maxx]=max(Pks); [temp2,maxy]=max(tempi); if temp2==0.0 || abs(maxx(maxy)-WinSizeX/2-1) WinSizeX/3.0 || abs(maxy-WinSizeY/2-1)WinSizeY/3.0 VX(il, jl)=0; VY(ilz jl)=0; Corr(ilz jl)=0; else TX=maxx(maxy); TY=maxy; Corr (ilz jl)=temp2; % 用高斯拟合峰值的精确位置 yO=Pks (TX-1, TY) ; %横向 yl=Pks(TX,TY); y2=Pks(TX+1,TY); CX1=(y0-y2)/(2.0*(y0+y2-2.0*y1)); yO=Pks (TX.TY-l) ; %纵向 yl=Pks(TX,TY); y2=Pks(TX,TY+1); CY1=(y0-y2)/(2 ?0*(y0+y2-2.0*yl)); yO=Pks (TX-1ZTY-1) ; %第…对角线 yl=Pks(TX,TY); y2=Pks(TX+1ZTY+1); CX2= (y0-y2)/(2.828* (y0+y2-2.0*y1)); CY2= (y0-y2)/(2.828*(y0+y2-2.0*y1)); yO=Pks (TX+1, TY-1) ; %笫-:对角线 yl=Pks(TX,TY); y2=Pks(TX-1,TY+1); CX3=-(yO-y2)/(2.828*(yO+y2-2.0*y1)); CY3=(y0-y2)/(2.828*(y0+y2-2.0*yl));