bmiafig08.m

来自「beamlet变化的工具箱」· M 代码 · 共 69 行

function BMIAfig08(OpenNewWindow,HaveTitle,LoadData)

if nargin<2
    HaveTitle=1;
end
if nargin<1
    OpenNewWindow=1;
end
if OpenNewWindow
    figure
end

% BMIAfig08: 
%clc; close all;
fprintf('\n');
disp('Figure 8')
disp('Illustrating the nearly identical sensitivity of')
disp('the maximum beam statistic and the maximum beamlet statistic.')
disp('Eleven different Signal-to-noise Ratios (SNR) are tested.')
disp('The lines with \Diamond, \oplus, \ast, \Box and \times')
disp('indicate: ')
disp(' o the pixel SNR of the underlying line object; ')
disp(' o the median of maximum raw pixel intensities; ')
disp(' o the threshold of detectability for .05 false-alarm rate; ')
disp(' o the median of the maximum standardized beamlet coefficients; and ')
disp(' o the median GLRT score. ')
disp('All the medians refer to results of 100 Monte Carlo simulations. ')
disp('The size of the test image is $n=128$. ')
disp('The underlying line segment is given in the paper.')
fprintf('\n');
if ~LoadData
disp('Warning: it takes a while to run the simulations.')
%k = input('Would you like to use a saved simulation result (0-no/1-yes)?');
end
fprintf('\n');

if LoadData, 
  disp('You need to download the file "fig08_128_100.mat" on the web.')
  disp('To download, visit "http://www.isye.gatech.edu/~beamlab/data/".')
  disp('Install the data file at [Your BMIALab Directory]/Applications/MaxBeamletCoeff/.')
  disp('Now waiting...');
  pause;
  disp('Loading the data...');

  fig08draw; 
else%if k==0, 
  fig08run;
  fig08draw; 
end;

fprintf('\n');
disp('For more information, visit subdirectory [Applications/MaxBeamletCoeff].');
disp('*** END Figure 8 ***');


%
% Part of BeamLab Version:200
% Built:Friday,23-Aug-2002 00:00:00
% This is Copyrighted Material
% For Copying permissions see COPYING.m
% Comments? e-mail [e-mail]
%
%
% Part of BeamLab Version:200
% Built:Saturday,14-Sep-2002 00:00:00
% This is Copyrighted Material
% For Copying permissions see COPYING.m
% Comments? e-mail [e-mail]
%