功能描述：
clear; 
clc; 
load Net 
 
Rt=100; 
St=115; 
 
Image2=imread('dst8\h10001.jpg'); 
Image2=imresize(Image2,[240,320]); 
h=initObj(Image2,0.01); 
 
oldsize=0; 
newsize=0; 
time=1; 
OldOutImage=zeros(123,163); 
OldOutImage1=zeros(240,320); 
OldDif=0; 
Dif=0; 
OldX=0; 
OldY=0; 
newX=0; 
newY=0; 
isfire=false; 
newcosb=0; 
oldcosb=0; 
%% 
tempfeature=[0;0;0]; 
%% 
for i=1:200 
if i<10 
  str=strcat('dst2\h1000',int2str(i),'.jpg'); 
elseif i<100 && i>=10 
    str=strcat('dst2\h100',int2str(i),'.jpg'); 
else 
    str=strcat('dst2\h10',int2str(i),'.jpg'); 
end 
     
Outimg=zeros(size(Image2(:,:,1))); 
Image1=imread(str); 
Image1=imresize(Image1,[240,320]); 
Image1(:,:,1) = medfilt2(Image1(:,:,1),[3 3]); 
Image1(:,:,2) = medfilt2(Image1(:,:,2),[3 3]); 
Image1(:,:,3) = medfilt2(Image1(:,:,3),[3 3]); 
outimg=getObj(Image1,h); 
outimg=~outimg; 
outimg=~outimg; 
se = strel('disk',2); 
outimg=imerode(outimg,se); 
ImageR=double(Image1(:,:,1)); 
ImageG=double(Image1(:,:,2)); 
ImageB=double(Image1(:,:,3)); 
HSVimg=double(rgb2hsv(Image1))*255; 
YCBCRimg=double(rgb2ycbcr(Image1)); 
oldsize=newsize; 
newsize=0; 
for m=1:size(ImageR,1) 
    for n=1:size(ImageR,2) 
        if outimg(m,n)==1 
            if ImageR(m,n)>Rt && YCBCRimg(m,n,3)>YCBCRimg(m,n,2) %&&  ImageR(m,n)>Image2(m,n,1) 
                if ImageR(m,n)>ImageG(m,n) && ImageG(m,n)>ImageB(m,n) 
                    if  HSVimg(m,n,2)>=((255-ImageR(m,n))*St/Rt) && HSVimg(m,n,2)>=90  
                        Outimg(m,n)=255; 
                        newsize=newsize+1; 
                    end 
                end 
            end 
        end 
    end 
end 
%% ͼ���� 
se = strel('disk',2,0); 
Outimg = imdilate(Outimg,se); 
se = strel('disk',2,0); 
Outimg = imerode(Outimg,se); 
Outimg = imerode(Outimg,se); 
 
% Label=bwlabel(Outimg);              
% area_num=regionprops(Label,'Area'); 
% if isempty(area_num) 
%     continue; 
% end 
%  len_area=length(area_num); 
%  area_bounding=regionprops(Label,'BoundingBox'); 
%  areasize=0; 
%  largest=1; 
%  for oo=1:len_area 
%      if area_bounding(oo).BoundingBox(3)*area_bounding(oo).BoundingBox(3)>areasize 
%          areasize=area_bounding(oo).BoundingBox(3)*area_bounding(oo).BoundingBox(4); 
%          largest=oo; 
%      end 
%  end 
% area_bounding(largest).BoundingBox(2)=round(area_bounding(largest).BoundingBox(2)); 
% area_bounding(largest).BoundingBox(4)=round(area_bounding(largest).BoundingBox(4)); 
% area_bounding(largest).BoundingBox(1)=round(area_bounding(largest).BoundingBox(1)); 
% area_bounding(largest).BoundingBox(3)=round(area_bounding(largest).BoundingBox(3)); 
%  tempimg=zeros(size(Outimg)); 
%  tempimg(area_bounding(largest).BoundingBox(2)+1:area_bounding(largest).BoundingBox(2)+area_bounding(largest).BoundingBox(4)-1,area_bounding(largest).BoundingBox(1)+1:area_bounding(largest).BoundingBox(1)+area_bounding(largest).BoundingBox(3)-1)= tempimg(area_bounding(largest).BoundingBox(2)+1:area_bounding(largest).BoundingBox(2)+area_bounding(largest).BoundingBox(4)-1,area_bounding(largest).BoundingBox(1)+1:area_bounding(largest).BoundingBox(1)+area_bounding(largest).BoundingBox(3)-1)+1; 
%  Outimg=Outimg.*tempimg; 
  
  
  
%% ������� 
[AY,BY,CY,DY]=dwt2(Image2(:,:,1),'db4'); 
AY=BY+CY+DY; 
[AN,BN,CN,DN]=dwt2(Image1(:,:,1),'db4'); 
AN=BN+CN+DN; 
DIN=AN-AY; 
Smok=zeros(size(DIN)); 
for k2=1:size(DIN,1) 
    for k3=1:size(DIN,2) 
        if DIN(k2,k3)<0 
        Smok(k2,k3)=255; 
        else 
            Smok(k2,k3)=0; 
        end 
    end 
end 
 
 
%% 
ImageW=double(Image1(:,:,1)).*(Outimg(:,:,1)/255); 
[A,B,C,D]=dwt2(ImageW,'db4'); 
A=B+C+D; 
Dif=xor(A,OldOutImage(:,:,1)); 
OldOutImage(:,:,1)=A; 
%% �������ƶȼ��� 
jiao=Outimg & OldOutImage1; 
bing=Outimg | OldOutImage1; 
jiao=size(find(jiao),1); 
bing=size(find(bing),1); 
if bing~=0; 
sim11=jiao/bing; 
else 
    sim11=0; 
end 
OldOutImage1(:,:,1)=Outimg(:,:,1); 
Dif=size(find(Dif),1); 
diff=abs(Dif-OldDif)/(OldDif+1); 
OldDif=Dif; 
%if abs(newsize-oldsize)/(oldsize+1)>0.2 && diff>0.2 
%% ������涯̬���� 
fireX=0; 
fireY=0; 
sumX=0; 
sumY=0; 
for k2=1:size(ImageR,1) 
    for k3=1:size(ImageR,2) 
         fireX=Outimg(k2,k3)*HSVimg(k2,k3,3)*k3+fireX; 
         fireY=Outimg(k2,k3)*HSVimg(k2,k3,3)*k2+fireY; 
         sumX=sumX+Outimg(k2,k3)*HSVimg(k2,k3,3); 
         sumY=sumY+Outimg(k2,k3)*HSVimg(k2,k3,3); 
    end 
end 
if sumX~=0 && sumY~=0 
fireX=fireX/sumX; 
fireY=fireY/sumY; 
else 
    fireX=0; 
    fireY=0; 
end 
 
newX=fireX; 
newY=fireY; 
OldX=newX; 
OldY=newY; 
%% ����Ƕȼ�� 
[Y2,X2]=find(Outimg(:,:,1)); 
[C,I]=min(X2); 
opintleft=[0+Y2(I),0+X2(I)]; 
[C,I]=max(X2); 
opintright=[Y2(I), X2(I)]; 
[C,I]=min(Y2); 
opinttop=[Y2(I),X2(I)]; 
[C,I]=max(Y2); 
opintdown=[Y2(I),X2(I)]; 
 
if ~isempty(opinttop) && ~isempty(opintdown) && ~isempty(opintright) 
AE=sqrt((opinttop(1)-opintright(1))^2+(opintright(2)-opinttop(2))^2); 
BE=sqrt((opintdown(1)-opintright(1))^2+(opintdown(2)-opintright(2))^2); 
CE=sqrt((opinttop(1)-opintdown(1))^2+(opintdown(2)-opinttop(2))^2); 
newcosb=-(AE^2+BE^2-CE^2)/(2*AE*BE); 
end 
cosb=acos(abs(newcosb-oldcosb))*180/pi;  
cosb=90-cosb; 
if cosb>30 
    isfire=true; 
end 
 
%% ���ֱ������� 
if ~isempty(opinttop) && ~isempty(opintdown) && ~isempty(opintright) 
bao=Outimg(opinttop(1):opintdown(1),opintleft(2):opintright(2),1); 
sizare=size(find(bao),1); 
sizeare=sizare/((opintdown(1)-opinttop(1))*(opintright(2)-opintleft(2))); 
else 
sizeare=0;     
end 
 
%% �����ж�     
if i>=21 
tempfeature=tempfeature+[diff;sim11;sizeare]; 
if mod(i,5)==0  
    tempfeature=real(tempfeature)./5; 
    if cosb~=0 && diff~=0 && sim11 ~=0 && sizeare~=0 
    XT = sim(net,tempfeature); 
    XT = full(compet(XT)); 
    if XT(1)==1 && newsize>=400 
         time=time+1  
     end 
    end 
    tempfeature=[0;0;0]; 
      
end 
 
if time>=4 
    disp('fire'); 
end 
end 
oldcosb=newcosb; 
 
figure(1);             
imshow(uint8(Image1)); 
drawnow;  
figure(2); 
imshow(uint8(Outimg)); 
drawnow;  
end

 

聯系：highspeedlogic

QQ ：1224848052

微信：HuangL1121

郵箱：1224848052@qq.com

網站：http://www.mat7lab.com/

網站：http://www.hslogic.com/

微信掃一掃：