cfar恒虚警检测

clear;
clc;
load T72_mag.mat
I=mobandata{150};
subplot（231）
imshow（I）;
title（‘原图像‘）
%求图像背景的均值和方差
[mn]=size（I）;
N=m*n;
s=0;
v=0;
for i=1:40
    for j=1:128
        s=s+I（ij）;
        v=v+I（ij）*I（ij）;
    end    
end

for i=80:128
    for j=1:128
        s=s+I（ij）;
        v=v+I（ij）*I（ij）;
    end
end

for i=1:128
    for j=1:44
        s=s+I（ij）;
        v=v+I（ij）*I（ij）;
    end
end

for i=1:128
    for j=86:128
        s=s+I（ij）;
        v=v+I（ij）*I（ij）;
    end
end
mean=s/（2*40*128+2*48*128）;
var=v/（2*40*128+48*128）-（s/（2*40*128+2*48*128））^2;
sigma=sqrt（var）;
T=raylinv（0.999sigma）;
for i=1:m
    for j=1:n
        if I（ij）>T
            I（ij）=1;
        else I（ij）=0;
        end
    end
end
subplot（232）
imshow（I）;
title（‘阈值判别后的图象‘）
for i=3:m-2
    for j=3:n-2
        count=0;
        for x=i-2:i+2
            for y=j-2:j+2
              count=count+I（xy）;
            end
        end
        if count<9