Fast Radial Symmetry Transform/快速径向对称变换（代码）

<span style="font-size:18px;">function [ filtered ] = frst2d( image, radii, alpha, stdFactor, mode )
%FRST2D Applies Fast radial symmetry transform to image
%   Check paper Loy, G., & Zelinsky, A. (2002). A fast radial symmetry transform for detecting points of interest. Computer Vision?ECCV 2002.
%   Three different modes are available: 'bright', 'dark', 'both'

% Example:
% I = imread('iris.jpg');
% Id = rgb2gray(I);
% IdSmall = imresize(Id, 0.25);
% 
% f = frst2d(IdSmall, 13, 2, 0.1, 'dark');
% 
% subplot(2,1,1);
% imshow(IdSmall, []);
% subplot(2,1,2);
% imshow(f,[]);

bright = false;
dark = false;
if (strcmp(mode, 'bright'))
    bright = true;     %寻找亮区域的对称
elseif (strcmp(mode, 'dark'))
    dark = true;
elseif (strcmp(mode, 'both'))
    bright = true;                   %寻找暗区域的对称
    dark = true;
else
    error('invali