【图像检测】基于直方图实现火灾检测matlab代码
1 简介
直方图是表达一副图像灰度级分布情况的关于统计信息的灰度级函数。它反映了灰度级与出现这种灰度
的概率间的对应关系,其横坐标是灰度级,纵坐标是该灰度级出现的频度。对于高对比度图像的直方
图,其灰度级的分布范围宽且较均匀,而对比度低的图像,其直方图分布区间窄且集中于灰度级的某一
部分,因此直方图可以表征图像的一些重要特征信息。直方图均衡化处理是空域处理技术的基础。它是
通过灰度值的累积函数把原图像的直方图变成均匀分布的直方图,即把相对集中的灰度级分布进行扩
展,通过增加像素灰度级分布的动态范围,使对比度得到增加[76],达到增强图像整体对比度的效果。a
2 部分代码
clear all; close all; clc;
img=double(imread('flame2.jpg'));
imshow(img,[]);
[m n]=size(img);
img=sqrt(img); %伽马校正
% HOG特征提取
%下面是求边缘
fy=[-1 0 1]; %定义竖直模板
fx=fy'; %定义水平模板
Iy=imfilter(img,fy,'replicate'); %竖直边缘
Ix=imfilter(img,fx,'replicate'); %水平边缘
Ied=sqrt(Ix.^2+Iy.^2); %边缘强度
Iphase=Iy./Ix; %边缘斜率,有些为inf,-inf,nan,其中nan需要再处理一下
%下面是求cell
step=16; %step*step个像素作为一个单元
orient=9; %方向直方图的方向个数
jiao=360/orient; %每个方向包含的角度数
Cell=cell(1,1); %所有的角度直方图,cell是可以动态增加的,所以先设了一个
ii=1;
jj=1;
for i=1:step:m %如果处理的m/step不是整数,最好是i=1:step:m-step
ii=1;
for j=1:step:n %注释同上