Question

嗨我对图像进行FFT操作后的幅度谱现在我想画一个掩模。然而，使用霍夫线等，我得到太多的线，掩码失败。以下是我需要创建的图像和蒙版的链接

https://www.dropbox.com/s/vgbmf3ju50uf821/fftimg4.jpg?dl=0 https://www.dropbox.com/s/ktl0yghan9t868o/image4.jpg?dl=0

提前致谢

这是代码

import cv2
import numpy as np
from matplotlib import pyplot as plt
from scipy.misc import imshow, imsave, imread ,imsave
import mahotas


scaler = .5
img = imread("color2.jpeg")[:,:,:3]
mask = imread("image4.jpg")[:,:,:3]

img = cv2.resize(img, (0,0), fx=scaler, fy=scaler)
imggray = np.mean(img, -1)
imfft = np.fft.fft2(imggray)
mags = np.abs(np.fft.fftshift(imfft))
angles = np.angle(np.fft.fftshift(imfft))
visual = np.log(mags)
visual2 = (visual - visual.min()) / (visual.max() - visual.min())*255
cv2.imshow('Visual 2', visual2.astype(np.uint8))

imsave('fftimg4.jpg',visual2)

height,width,depth = img.shape
masking = np.zeros((height,width))

vis = visual2.astype(np.uint8)
edges = cv2.Canny(vis,50,180)
cv2.imshow('Gradient', edges)

lines = cv2.HoughLines(edges,1,np.pi/180,200,200,10)
for rho,theta in lines[0]:
   a = np.cos(theta)
   b = np.sin(theta)
   x0 = a*rho
   y0 = b*rho
   x1 = int(x0 + 1000*(-b))
   y1 = int(y0 + 1000*(a))
   x2 = int(x0 - 1000*(-b))
   y2 = int(y0 - 1000*(a))

   cv2.line(masking,(x1,y1),(x2,y2),(255,255,255),1)

   cv2.imshow('HoughLines', masking)

   mask = (np.mean(mask,-1) > 20)
   visual[mask] = np.mean(visual)
   newmagsshift = np.exp(visual)
   newffts = newmagsshift * np.exp(1j*angles)
   newfft = np.fft.ifftshift(newffts)
   imrev = np.fft.ifft2(newfft)
   newim2 = 255 - np.abs(imrev).astype(np.uint8)
   imsave("fftimg2.jpg",  newim2  )

   cv2.imshow('Image without Lines', newim2.astype(np.uint8))

cv2.waitKey(0)

Answer 1

首先，请澄清一下。在fft中你至少有三个垂直的＆＃39;线条和一堆水平＆＃39;，这是正确的吗？另外，您真的需要自动执行此操作，还是可以手动完成此过滤？您的霍夫线算法是否为您提供了一系列线的方程式？

似乎你可以过滤掉那些经过零的线，因为这是你所期望的。通过零并且也最接近预期角度（例如，在您的示例中为-5度）。

如何在FFT幅度谱上绘制掩模线

1 个答案: