火焰高度和频率.docx

一 火焰高度1截取火焰视频，要有火焰图像背景；2 matlab处理（total）获得火焰二值图；并生成average、combine和output三个文件。其中output为Tecplot后处理文件。3用Tecplot打开output。去mesh点contour. 4去除火焰多余的部分。点击活动按钮出现再点击more,出现移除备选框中火焰多余的部分。只保留0.5,0.9，0.95三个高度的火焰。然后记录火焰概率为0.5时的火焰最高最低点像素并相减得出火焰高的像素值。该像素与标度像素值（火焰图像中作为标尺的图像）作比较则可得出实际火焰平均高度。附,高度处理程序fire=mmreader(t.avi);nframes=get(fire,numberofframes);background=imread(0.jpg);for k=1:nframessingleframe=read(fire,k);pic=imsubtract(singleframe,background);p=rgb2gray(pic);level=0.08;i=im2bw(p,level);str= [num2str(k) .jpg];imwrite(i,str,jpg);end[height,width]=size(i);ishow=zeros(height,width);for k=1:nframesstr= [num2str(k) .jpg];i=imread(str);t=double(i);ishow=ishow+t;endishow=ishow./nframes;imwrite(ishow,combine.png);ishow=ishow./255;imwrite(ishow,average.png);fid=fopen(output.dat,wt);fprintf(fid,ZONE I=%d J=%d\n,width,height);for i=1:heightfor j=1:widthfprintf(fid,%14.8e\t%14.8e\t%14.8e\n,j,i,ishow(height+1-i,j));endendfclose(fid);二火焰频率1把待分析视频‘t.avi’和程序文件‘Beta0915V1’放入同一工作文件夹。2 matlab打开‘Beta0915V1’。3 点击运行该m文件。选取目标区域。（该程序中以第六帧为样本帧做相关性计算）点击鼠标左键并按住选取区域选中后双击。开始计算。、4 计算完成后工作目录中会生成一个文件。将里面的数据拷到origin中。在第一列插入序号数据。点击A(x)，全选中第一列，后点右键得出：再将该行数除以该视频的帧率得出视频的时间序列：点右键，选择set column values. 出现：点击Col（A）。选择col(A).后下面出现将其处以帧率（本例为25帧/s）.然后点ok。第一列即变为：5 作图。选中两列点击右下角中的获得：6 FFT转换。出现：点击ok.数据表中出现FFT转化的数据和曲线。曲线图中第一张：双击y轴数据，出现：更改起始值（From）和终止值（To）：起始值一般为0，终止值适当选择直到发现有除零点以外的峰值出现，若出现峰值即为火焰振荡频率。本例中未发现。程序：%Beta0915V1%程序功能：计算火焰图像的相关系数%读取视频%以第k帧为参考%在图像中截去干扰部分%第k帧GRB输出%保存第k帧GRAY输出%各帧转化为灰度图片（不输出）%计算所有帧与第k帧相关系数r%并输出成文本文件%Beta0915V1 %计算图像火焰部分的相关系数mv=mmreader(t.AVI);nframes=get(mv,NumberOfFrames);nframesget(mv)k= 6 ; %以第k张图片作为参考imwrite(read(mv,k),[RGB_,num2str(k),.png])RGB=imread([RGB_,num2str(k),.png]);[RGB rect]=imcrop(RGB); %裁剪出有效图片范围,rect,[xmin ymin width height]A = rgb2gray(RGB); %转化为灰度图片imwrite(A,[GRAY_,num2str(k),.png]);for k=1:600 B=rgb2gray(read(mv,k)); B=imcrop(B,rect); r(k)=corr2(A,B); 1000-k %指示计算进度endfid = fopen(fft.txt, wt);fprintf(fid, %8.6f\n, r);fclose(fid)