我已经编写了一些 Python(抱歉不是 c++)代码,这些代码将允许进行屏蔽自适应阈值处理。它不是很快,但它可以满足您的需求,并且您可以将它用作 C++ 代码的基础。它的工作原理如下:
- 将图像中的蒙版像素设置为零。
- 确定卷积块内每个像素的未屏蔽邻居的数量。
- 执行卷积,并通过块内未屏蔽的邻居数对其进行平均。这会产生像素邻域块内的平均值。
- 阈值,通过将图像与平均邻域值进行比较,
mean_conv
- 重新添加图像的被屏蔽(非阈值)部分。
图像显示,初始图像,蒙版,最终处理的图像。
代码如下:
import cv
import numpy
from scipy import signal
def thresh(a, b, max_value, C):
return max_value if a > b - C else 0
def mask(a,b):
return a if b > 100 else 0
def unmask(a,b,c):
return b if c > 100 else a
v_unmask = numpy.vectorize(unmask)
v_mask = numpy.vectorize(mask)
v_thresh = numpy.vectorize(thresh)
def block_size(size):
block = numpy.ones((size, size), dtype='d')
block[(size - 1 ) / 2, (size - 1 ) / 2] = 0
return block
def get_number_neighbours(mask,block):
'''returns number of unmasked neighbours of every element within block'''
mask = mask / 255.0
return signal.convolve2d(mask, block, mode='same', boundary='symm')
def masked_adaptive_threshold(image,mask,max_value,size,C):
'''thresholds only using the unmasked elements'''
block = block_size(size)
conv = signal.convolve2d(image, block, mode='same', boundary='symm')
mean_conv = conv / get_number_neighbours(mask,block)
return v_thresh(image, mean_conv, max_value,C)
image = cv.LoadImageM("image.png", cv.CV_LOAD_IMAGE_GRAYSCALE)
mask = cv.LoadImageM("mask.png", cv.CV_LOAD_IMAGE_GRAYSCALE)
#change the images to numpy arrays
original_image = numpy.asarray(image)
mask = numpy.asarray(mask)
# Masks the image, by removing all masked pixels.
# Elements for mask > 100, will be processed
image = v_mask(original_image, mask)
# convolution parameters, size and C are crucial. See discussion in link below.
image = masked_adaptive_threshold(image,mask,max_value=255,size=7,C=5)
# puts the original masked off region of the image back
image = v_unmask(original_image, image, mask)
#change to suitable type for opencv
image = image.astype(numpy.uint8)
#convert back to cvmat
image = cv.fromarray(image)
cv.ShowImage('image', image)
#cv.SaveImage('final.png',image)
cv.WaitKey(0)
写完这篇文章后,我发现this great link 有很多图像示例,解释得很好,我在上面的示例中使用了他们的文本图像。
注意。 scipy signal.convolve2d() 似乎不尊重 Numpy 掩码,因此上述解决方法是必要的。