2016-04-24 21:26:12 -0500 received badge ● Enthusiast 2016-04-21 14:59:49 -0500 asked a question how to compute image gradient using Sobel operator Hi everyone, I'm trying to get an orientation map from a fingerprint image but I don't know if the way I compute the gradient in each direction(x and y) is good. That's my python code Dx = cv2.Sobel(norm_gray_img,cv2.CV_16S,1,0,3)#in x direction Dy = cv2.Sobel(norm_gray_img,cv2.CV_16S,0,1,3)#in y direction  I also see in some question the use of "cv2.convertScaleAbs" but I don't know why it is used. all the code are: def local_vx_vy(Dx,Dy,point,img_size,block_size=(17,17)): #point = (x,y) x on col navigation, y on row navigation # img_size = (row, col) size of the image row is the nombre # of lignes and col the number of column x,y = point row, col = img_size h,w = block_size ofi, ofj = w/2,h/2 Gxy = Gxx = Gyy = 0 for j in range(-ofj,ofj+1): for i in range(-ofi,ofi+1): if 0 <= x+i and x+i <= col-1 and 0 <= y+j and y+j <= row-1: #print Dx[y+j][x+i], Dy[y+j][x+i], "dx,dy" Gxy += (Dx[y+j][x+i])*(Dy[y+j][x+i]) Gxx += Dx[y+j][x+i]**2 Gyy += Dy[y+j][x+i]**2 vx = (2*Gxy)/(w*h) vy = (Gxx*Gyy)/(w*h) return (vy,vx) def compute_orientation(norm_gray_img,block_size=(17,17)): h,w = block_size row,col = norm_gray_img.shape ofj, ofi = h/2, w/2 print norm_gray_img Dx = cv2.Sobel(norm_gray_img,cv2.CV_16S,1,0,3) Dx = cv2.convertScaleAbs(Dx).astype(np.float) Dy = cv2.Sobel(norm_gray_img,cv2.CV_16S,0,1,3) Dy = cv2.convertScaleAbs(Dy).astype(np.float) Vy = norm_gray_img[0:row-h+1:w, 0:col-h+1:h] Vx = Vy.copy() angle = Vx.copy() r = c = -1 for j in range(0,row-h+1,h): r += 1 for i in range(0,col-w+1,w): c += 1 Vy[r][c],Vx[r][c] = local_vx_vy(Dx,Dy,(i,j),norm_gray_img.shape,block_size)# c = -1 theta = 0.5*(np.arctan2(Vy,Vx) + M_PI) Vx, Vy = cv2.blur(np.cos(theta),(5,5))/25, cv2.blur(np.sin(theta),(5,5))/25 theta = 0.5*(np.arctan2(Vy,Vx) + M_PI) return angle  with the code above and the following plotting function the result is not what I expect  def plot_point(point, angle, length, ax): ''' point - Tuple (x, y) angle - Angle you want your end point at in degrees. length - Length of the line you want to plot. Will plot the line on a 10 x 10 plot. ''' # unpack the first point x, y = point # find the end point #print angle endx = x + length endy = length*math.tan(angle)+y ax.plot([x, endx], [y,endy],color='red ') 2016-04-19 13:51:17 -0500 commented answer Fingerprint orientation map with this I can't see how it is possible to get an orientation map 2016-04-19 13:33:34 -0500 asked a question implementation of an orientation map for fingerprint enhancement I'm implementing a function to get an orientation map of a fingerprint image using opencv and python but something is wrong I don't know what this is my code def compute_real_orientation(array_I, w=17, h=17, low_pass_filter=cv2.blur,filter_size=(5,5), blur_size=(5,5),**kwargs): row, col = array_I.shape array_I = array_I.astype(np.float) Ox = array_I[0:row-h+1:h,0:col-w+1:w].copy() Ox[:] = 0.0 Vx = Ox.copy() Vy = Vx.copy() Oy = Ox.copy() angle = Vx.copy()#array to contain all the 17*17 blocks's orientatons c = r = -1 for i in xrange(0, row-h+1, h): r+=1 for j in xrange(0, col-w+1, w): c+=1 Dx = cv2.Sobel(array_I[i:i+h,j:j+w],-1,1,0)#gradient component x for a 17*17block Dy = cv2.Sobel(array_I[i:i+h,j:j+w],-1,0,1)#gradient component y for 17*17 block for k in range(0,h): for l in range(0,w): Vy[r][c] += 2*(Dx[k][l])*(Dy[k][l]) Vx[r][c] += ((Dx[k][l])**2 - (Dy[k][l])**2) #get the orientation angle for the given 16*16 block angle[r][c] = 0.5*(math.atan(Vy[r][c]/Vx[r][c])) c = -1 #smoothing process of the whole array angle row, col = angle.shape for i in range(0, row): for j in range(0, col): Ox[i][j] = math.cos(2*angle[i][j]) Oy[i][j] = math.sin(2*angle[i][j]) Ox = low_pass_filter(Ox, blur_size) Oy = low_pass_filter(Oy, blur_size) for i in range(0, row): for j in range(0, col): #take the final orientation of all 17*17 blocks angle[i][j] = 0.5*math.atan(Oy[i][j]/Ox[i][j]) return angle  I'm implementing the following algorithm algorithm at 2.4 Orientation Image section but my code is not working properly, I don't get the right orientation map. Can any one help me troubleshooting this? B.R