function [Q]=cir_quad(I,xc,yc,rc,t)
%求中心為（xc,yc），rc-t/2到rc+t/2圓環(huán)內(nèi)積分
It = I;
%%
% 1、計(jì)算點(diǎn)到圓心距離矩陣
for r = 1 : size(I,1)
    for c = 1 : size(I,2)
        Dist(r,c) = sqrt((r-yc)^2+(c-xc)^2);
        if abs(r-yc) > abs(c-xc)
            It(r,c) = 0;
        end
    end
end

%%
% 2、對環(huán)內(nèi)“1”值點(diǎn)求和
Q = sum(It(find((Dist<=rc+t/2)&(Dist>=rc-t/2))));

%%
% 3、歸一化
%Q = Q/(4*rc*t);