一分钟了解"用matlab计算信杂比SCR 以及 背景区域的标准差"

%求SCR和方差的方法

clear,clc,close all

imgName=input('请输入bmp文件名,要有单引号,要有后缀:');

img=imread(imgName);

figure,imshow(img,[])

leftupPointRow=input('请输入目标的左上角的行数(第二坐标):');

leftupPointCol=input('请输入目标的左上角的列数(第一坐标):');
rightdownPointRow=input('请输入目标的右下角的行数(第二坐标):');

rightdownPointCol=input('请输入目标的右下角的列数(第一坐标):');
backgroundAreaWidth=10;
targetVector=img(leftupPointRow:rightdownPointRow,leftupPointCol:rightdownPointCol)

targetSize=size(targetVector);

targetVector=double( targetVector(:) );

backgroundVector1=img(leftupPointRow-10:rightdownPointRow+10,leftupPointCol-backgroundAreaWidth:leftupPointCol-1)

backgroundVector2=img(leftupPointRow-10:rightdownPointRow+10,rightdownPointCol+1:rightdownPointCol+10)

backgroundVector3=img(leftupPointRow-10:leftupPointRow-1,leftupPointCol:rightdownPointCol)

backgroundVector4=img(rightdownPointRow+1:rightdownPointRow+10,leftupPointCol:rightdownPointCol)

backgroundVector=[backgroundVector1(:);backgroundVector2(:);backgroundVector3(:);backgroundVector4(:)];

backgroundVector=double(backgroundVector);
scr=(mean(targetVector)-mean(backgroundVector))/std(backgroundVector);
fprintf(['目标的尺寸是',num2str(targetSize(1)),'x',num2str(targetSize(2)),'\n']);

fprintf(['SCR=',num2str(scr),'\n']);

fprintf(['std(background)=',num2str(std(backgroundVector)),'\n']);