- 1、本文档共11页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
火焰高度和频率
一 火焰高度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)
您可能关注的文档
最近下载
- 《市级国土空间生态修复规划编制指南》.pdf
- 坠床跌倒护理持续质量改进ppt.pptx
- 英语新闻ppt(含讲稿).pptx VIP
- 感觉信息处理及自我调节能力检核表Sensoryprocessingandself-regulationchecklist,SPSRC.docx
- 做守纪知纪明纪的模范.pptx VIP
- 餐厅服务方案.docx VIP
- 苏教版三年级数学上册知识点整理.docx VIP
- 【微专题】2014-2024高考语文文言文挖空训练( I 卷、II 卷、甲卷、乙卷)教师版.doc VIP
- 公司安全事故隐患内部举报、报告奖励制度范本.docx
- 第一单元《中国共产党的领导》练习题试题含答案.docx VIP
文档评论(0)