function?[imgfilter]=fbp2（rthetaxpimgrowsimgcolsfilter_type）;
[nk]=size（r）;
xp_offset=abs（min（xp））+1;
width=2^nextpow2（n）;
tmpimg=zeros（imgrowsimgcols）;
ramlak=2*[0:（width/2-1）width/2:-1:1]‘/width;
if（filter_type==0）
????filter=ones（width1）;
elseif（filter_type==1）
????filter=ramlak;
elseif（filter_type==2）
????Sinc?=?abs（sinc（2*（0:width-1）/2:-1:1））;
????temp?=?[?sinc（1:（width/2））?sinc（width/2:-1:1）];
????filter=ramlak.*temp‘;
elseif（filter_type==3）
????cosine=abs（cos（2*pi*（0:width-1）/（2*width）））;
????filter=ramlak.*cosine;
elseif（filter_type==4）
????hamming=0.54-0.46*cos（2*pi*（width-1）/width）;
????temp=[hamming（width/2:width）?hamming（1:widyh/2-1）];
????filter=ramlak.*temp‘;
elseif（filter_type=