引用代码：
functionkydys=ziwen _ fangiantu (ziwentuxiang ) ) ) ) ) )。
%函数功能计算指纹方向图
%函数参数指纹图像Zhiwentuxiang
%函数返回值指纹方向图kydys
W=4； %窗口大小(2W 1) ) 2W 1) ) ) ) )。
W=4；
%Sobel为xsdxhd/p x_fangxiang=[-1 0 1；
-2 0 2；
-1 0 1] )；
y_fangxiang=[1 2 1；
0 0 0；
-1 -2 -1]；
GX=Zeros(sizezhiwentuxiang；
kydys=Zeros(sizezhiwentuxiang；
Fori=2: sizezhiwentuxiang (1，1 )- 1
forj=2: sizezhiwentuxiang (1，2 )- 1
Sum_x=0； %j侧kndqz/p Sum_y=0； %i方kndqz/p for k=-1 : 1
for r=-1 : 1
sum_x=sum_xx_fangxiang(k2，r 2 ) *zhiwentuxiang ) Ik，j r )；
sum_y=sum_y_fangxiang(k2，r 2 ) *zhiwentuxiang ) Ik，j r )；
结束
结束
GX(I，j )=Sum_x；
结束
结束
fori=w 1: sizezhiwentuxiang (1，1 ) -清洁的触发鼠/p forj=w 1: sizezhiwentuxiang ) 1，2 ) -清洁的触发鼠/p Vx=0；
Vy=0；
for io=-W :清洁的触发鼠/p for jo=-W :清洁的触发鼠/pvy=vy2*GX(Iio，j jo ) ) gy ) Iio，j jo )；
结束
结束
if Vx ~=0
th=atan(vy/VX )/2；
kydys(I，j )=th；
else
kydys(I，j )=pi/2；
结束
结束
结束