22

这个问题与这个问题有关:How to remove convexity defects in sudoku square

我试图nikie's answerMathematica to OpenCV-Python. 但我被困在程序的最后一步。

即我得到了正方形的所有交点,如下所示:

在此处输入图像描述

现在,我想将其转换为一个完美的正方形(450,450),如下所示:

在此处输入图像描述

(不要介意两个图像的亮度差异)。

问题: 我如何在 OpenCV-Python 中做到这一点?我正在使用cv2版本。

4

2 回答 2

35

除了 etarion 的建议,您还可以使用remap功能。我写了一个快速脚本来展示如何做到这一点。如您所见,这在 Python 中非常简单。这是测试图像:

扭曲的图像

这是翘曲后的结果:

扭曲的图像

这是代码:

import cv2
from scipy.interpolate import griddata
import numpy as np

grid_x, grid_y = np.mgrid[0:149:150j, 0:149:150j]
destination = np.array([[0,0], [0,49], [0,99], [0,149],
                  [49,0],[49,49],[49,99],[49,149],
                  [99,0],[99,49],[99,99],[99,149],
                  [149,0],[149,49],[149,99],[149,149]])
source = np.array([[22,22], [24,68], [26,116], [25,162],
                  [64,19],[65,64],[65,114],[64,159],
                  [107,16],[108,62],[108,111],[107,157],
                  [151,11],[151,58],[151,107],[151,156]])
grid_z = griddata(destination, source, (grid_x, grid_y), method='cubic')
map_x = np.append([], [ar[:,1] for ar in grid_z]).reshape(150,150)
map_y = np.append([], [ar[:,0] for ar in grid_z]).reshape(150,150)
map_x_32 = map_x.astype('float32')
map_y_32 = map_y.astype('float32')

orig = cv2.imread("tmp.png")
warped = cv2.remap(orig, map_x_32, map_y_32, cv2.INTER_CUBIC)
cv2.imwrite("warped.png", warped)

我想你可以谷歌并找到 griddata 的作用。简而言之,它进行插值,在这里我们使用它来将稀疏映射转换为密集映射,因为 cv2.remap 需要密集映射。我们只需要将值转换为 float32,因为 OpenCV 抱怨 float64 类型。请让我知道情况如何。

更新:如果您不想依赖 Scipy,一种方法是在代码中实现 2d 插值功能,例如,请参阅 Scipy 中 griddata 的源代码或更简单的类似http://inasafe.readthedocs .org/en/latest/_modules/engine/interpolation2d.html仅依赖于 numpy。不过,我建议为此使用 Scipy 或其他库,尽管我明白为什么只需要 CV2 和 numpy 对于这样的情况可能会更好。我想听听您的最终代码如何解决数独问题。

于 2012-04-29T18:34:39.350 回答
1

如果你有源点和端点(你只需要 4 个),你可以将它们插入 cv2.getPerspectiveTransform,并在 cv2.warpPerspective 中使用该结果。给你一个很好的平坦结果。

于 2017-10-02T03:36:54.720 回答