%----------------------------------------------------------
%angular_filter_bank
%precomputes?angular?filter?bank?and?also?generates?a?C
%header?file?with?filter?coefficients.?The?matlab?variable
%angf?is?stored?in?file?‘a(chǎn)ngular_filters.mat‘?and?the?C?
%variables?are?stored?in?‘a(chǎn)ngular_fiter.h‘.
%There?are?TSTEPS?number?of?raised?cosine?filter?that?span?the?
%interval?[0180]
%usage
%?angular_filter_bank（BWfname）
%BW????-?bandwidth?（radians）
%fname?-?file?name
%Contact:
%???ssc5@eng.buffalo.edu
%???www.eng.buffalo.edu/~ssc5
%References:
%B.G.Sherlock?and??D.M.Monro?and??K.Millard“Fingerprint?Enhancement?by
%directional?Fourier?Filtering“IEEE?Visual?Image?Signal?Processing
%141（2）?pp.?87--94?1994
%----------------------------------------------------------
function?angular_filter_bank（BWfname）
close?all;
%---------------
%parameters
%---------------
FFTN????=???32;
TSTEPS??=???12;?%15?degrees?interval
DELTAT??=???pi/TSTEPS;
%---------------
%precompute
%---------------
[xy]???=???meshgrid（-FFTN/2:FFTN/2-1-FFTN/2:FFTN/2-1）;
r???????=???sqrt（x.^2+y.^2）;
th??????=???atan2（yx）;
th（th<0）=???th（th<0）+2*pi;??%unsigned
filter??=???[];

%-------------------------
%generate?the?filters
%-------------------------
for?t0??=???0:DELTAT:（TSTEPS-1）*DELTAT
?????t1?????=?t0+pi;????????????????????????????????%for?the?other?lobe
?????%-----------------
?????%first?lobe
?????%-----------------
?????d??????????=?angular_distance（tht0）;
?????msk????????=?1+cos（d*pi/BW）;?
?????msk（d>BW）??=?0;
?????rmsk???????=?msk;??????????????????????????????%save?first?lobe

?????%-----------------
?????%second?lobe
?????%-----------------
?????d??????????=?angular_distance（tht1）;
?????msk????????=?1+cos（d*pi/BW）;?
?????msk（d>BW）??=?0;
?????rmsk???????=?（rmsk+msk）;

?????imagesc（rmsk）;pause;
?????rmsk???=?transpose（rmsk）;
?????filter?=?[filterrmsk（:）];
end;
?????%-----------------
?????%save?the?filters
?????%-----------------
?????angf??=?filter;
?????eval（sprintf（‘save?%s?angf‘fname））;

%----------------------
%write?to?a?C?file
%----------------------
fp?=?fopen（sprintf（‘%s.h‘fname）‘w‘）;
fprintf（fp‘{\n‘）;
for?i?=?1:size（filter2）
????i
????k?=?1;
????fprintf（fp‘{‘）;
????for?j?=?1:size（filter1）
????????fprintf（fp‘%f‘filter（ji））;
????????if（k?==?32）?k=0;?fprintf（fp‘\n‘）;?end;
????????k?=?k+1;
????end;
????fprintf（fp‘}\n‘）;
end;
fprintf（fp‘};\n‘）;
fclose（fp）;
%end?function?radial_filter_bank

%-----------------------------------------
%angular_distance
%computes?angular?distance-acute?angle?
%-----------------------------------------
function?d?=?angular_distance（tht0）
????d?=?abs（th-t0）;
????d?=?min（d2*pi-d）;
%end?function?angular_distance

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----

?????文件?????166272??2005-06-30?00:41??圖像增強(qiáng)源代碼\1_1.tif

?????文件?????166272??2005-06-30?00:41??圖像增強(qiáng)源代碼\1_2.tif

?????文件??????98488??2005-07-09?15:26??圖像增強(qiáng)源代碼\angular_filters_pi_2.mat

?????文件??????98488??2005-07-09?15:26??圖像增強(qiáng)源代碼\angular_filters_pi_4.mat

?????文件???????2774??2005-06-30?00:41??圖像增強(qiáng)源代碼\angular_filter_bank.m

?????文件????????414??2005-07-15?22:23??圖像增強(qiáng)源代碼\complex_gradient.m

?????文件???????1115??2005-07-17?20:19??圖像增強(qiáng)源代碼\compute_coherence.m

?????文件???????1019??2005-07-31?00:53??圖像增強(qiáng)源代碼\compute_region_mask.m

?????文件??????11678??2005-07-31?00:53??圖像增強(qiáng)源代碼\fft_enhance_cubs.m

?????文件????????574??2005-07-16?16:27??圖像增強(qiáng)源代碼\gabor_kernel.m

?????文件??????10832??2005-07-31?00:54??圖像增強(qiáng)源代碼\hong_enhancement.m

?????文件????????278??2005-07-17?20:29??圖像增強(qiáng)源代碼\imscale.m

?????文件???????2394??2005-06-30?00:41??圖像增強(qiáng)源代碼\orientation_image_rao.m

?????文件???????1026??2005-06-30?00:41??圖像增強(qiáng)源代碼\otsu_threshold.m

?????文件????????841??2005-07-09?17:48??圖像增強(qiáng)源代碼\pad_image.m

?????文件???????2625??2005-07-17?20:30??圖像增強(qiáng)源代碼\pseudo_matched_filter.m

?????文件?????131256??2005-01-31?12:45??圖像增強(qiáng)源代碼\radial_filters.mat

?????文件???????2143??2005-07-17?20:31??圖像增強(qiáng)源代碼\radial_filter_bank.m

?????文件???????2523??2005-07-31?00:54??圖像增強(qiáng)源代碼\segment_print.m

?????文件???????2954??2005-07-17?20:33??圖像增強(qiáng)源代碼\smoothen_frequency_image.m

?????文件???????1097??2005-07-17?20:33??圖像增強(qiáng)源代碼\smoothen_orientation_image.m

?????文件????????844??2005-07-17?20:34??圖像增強(qiáng)源代碼\view_orientation_image.m

?????目錄??????????0??2009-06-05?11:20??圖像增強(qiáng)源代碼

-----------?---------??----------?-----??----

???????????????705907????????????????????23