function?[f_anglef_grad]?=?FingerprintOrient_convolution_f（fr1）
%%%%%%%%%%%%%?function?description?%%%%%%%%%%%%%%%%%%%%%%%%
%%%%input????f:Array?of?img??;??r1:window?of?average?that?control
%%%%relationship?between?the?pix?and?its?neighbors
%%%%
%%%%output???f_angle:array?that?is?to?show?Orientation?field???;
%%%%?????????f_grad?:array?that?is?to?show?Grad?field?
%%%%????????????????????????????@author??Faith?15.12.22
%%%%????????????????????????????????????v2.0????16.8.3
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%%%%%?讀入圖片取出接口?%%%%%%%%%%%
f1?=?f;
r?=?r1;

%%%%獲取相關圖片信息%%%%%%
f2?=?im2double（f1）;
%?f2?=?f1;
[heightwidth]?=?size（f1）;

EPI	=?57.29578;?????%%%%%360/2pi，角度和弧度的換算

%%%%%%%%初始化矩陣%%%%%%%%%%%%%%%%
f_temp?=?f2;
f_angle?=?zeros（heightwidth）;
Vx?=?zeros（heightwidth）;
Vy?=?zeros（heightwidth）;
f_grad?=??zeros（heightwidth）;
%%%%%%%%%%%%%%%%%%%%%%%%%soble算子%%%%%%%%%%%%%%%%%%%%%%%%%%
%?soble_x?=?[-1?0?1;?-2?0?2;?-1?0?1];
%?soble_y?=?[1?2?1;?0?0?0;?-1?-2?-1];

%%%%%%%%%%%%%%?整體像素遍歷?%%%%%%%%%%%%%%%

%%%%%%%向量法%%%%%%%%%%%%%%%
f_linex?=?f_temp（:）;
f_tempy?=?f_temp‘;
f_liney?=?f_tempy（:）;

soble_x?=?[-1?-2?-1?zeros（12*height?-3）?1?2?1];
soble_y?=?[1?2?1?zeros（12*width-3）?-1?-2?-1];

f_num?=?size（f_linex1）;
soblex_num?=?size（soble_x2）;
sobley_num?=?size（soble_y2）;

%%%%%%%%%%%%%%%%%%基礎乘加%%%%%%%%%%%%%%%%%
%?for?i?=?1:（f_num?-soblex_num?+?1）
%?????f_lineVx（i）?=??f_linex（i:i+length（soble_x）-1）‘?*?soble_x‘;
%?end
%?for?i?=?1:（f_num?-sobley_num?+?1）
%?????f_lineVy（i）?=??f_liney（i:i+length（soble_y）-1）‘?*?soble_y‘;
%?end
%?
%?
%?f_x?=?[zeros（1（soblex_num-1）/2）?f_lineVx?zeros（1（soblex_num-1）/2）];
%?f_y?=?[zeros（1（sobley_num-1）/2）?f_lineVy?zeros（1（sobley_num-1）/2）];
%%%%%%%%%%

%%%%%%%%%%%%%%%%?conv?卷積%%%%%%%%%%%%%%%%%%
tic

f_lineVx?=?-conv（f_linex‘soble_x）;
f_lineVy?=?-conv（f_liney‘soble_y）;

f_x?=?f_lineVx（1+（soblex_num-1）/2:?length（f_lineVx）-（soblex_num-1）/2）;
f_y?=?f_lineVy（1+（sobley_num-1）/2:?length（f_lineVy）-（sobley_num-1）/2）;

%%%%%%%%%%%%

%%%%%%%%%回復圖像相應大小%%%%%%%%%%
f_xx?=?reshape（f_xheightwidth）;
f_yy?=?reshape（f_ywidthheight）;
f_yy?=?f_yy‘;

toc
%%%%%%%%%%%%%%%%%%%%%%%%%

%%%%%%%%%%%%%%%%原遍歷法%%%%%%%%%%%%%%%%%%%%%%%%%%%
tic

for?x?=?1:height
????for?y?=?1:?width
?????????lvx?=?0;
?????????lvy?=?0;
?????????gradsum?=?0;
?????????num?=?0;
????????????%%%%%%%%%%%%%%%%%%%r?*?r?窗下的處理%%%%%%%%%%%%%%%%%%%
????????????for?i?=?-r:r
????????????????%%%%%%%%%%%%%%%%%%?邊界處理?%%%%%%%%%%%%%%%%%%%%%
????????????????if（x?+?i?-1??height?）
????????????????????continue;
????????????????end
????????????????????for?j?=?-r:r
????????????????????????if?（y?+?j?-1??width）
????????????????????????????continue;
????????????????????????end
????????????????????????
????????????????????????%%%%%%%%%%%%%%%?soble?算子?x方向?和?y方向?的偏導
????????????????????????%%%%%%%%%%%%%%%?x?-1?0?1???y???1??2??1
???????????????

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件????????4198??2016-08-03?12:43??FingerprintOrient_convolution_f.m