0

为了 100% 可重现,整个代码(和使用的图像)如下。在此处输入图像描述

以下是代码中的步骤。

  1. 我将轮廓 ( [[386, 330], [398, 320], [382, 300], [370, 310], [386, 330]]) 叠加在图像上。这很好用。

  2. 使用 theta (= 90) 旋转图像warpAffine

  3. 在完全相同大小和相同旋转功能的虚拟图像中旋转轮廓中的每个点。

  4. 找到旋转点的坐标。

  5. 在旋转图像上叠加旋转点。

我希望他们排队,但他们没有

可能是什么原因?

#! /usr/bin/env python                                                                                             
import numpy as np
import cv2
from shapely.geometry.polygon import LinearRing

theta = 90
image_path = 'image.tiff'
orig_image = cv2.imread(image_path)

wnd_str = 'unrotated'
cv2.namedWindow(wnd_str,cv2.WINDOW_NORMAL)
cv2.moveWindow(wnd_str, 0, 0)

rows,cols, channels = orig_image.shape
print 'Loaded image shape {}'.format(orig_image.shape)

blank_image = np.zeros((rows, cols, 3), np.uint8)
print 'Blank image shape {}'.format(blank_image.shape)

M = cv2.getRotationMatrix2D((rows/2, cols/2), theta, 1.0)
rot_image = cv2.warpAffine(orig_image, M, (rows, cols), flags=cv2.INTER_CUBIC+cv2.BORDER_CONSTANT)
print 'Rotated image shape {}'.format(rot_image.shape)

white = (255, 255, 255)

#contours overlayed on unrotated image                                                                             
pts = [[386, 330], [398, 320], [382, 300], [370, 310], [386, 330]]
poly_pts = []
for p in pts: poly_pts.append(p)
poly_pts = np.array(poly_pts[0:-1], np.int32)
poly_pts = poly_pts.reshape((-1,1,2))

cv2.polylines(orig_image, [poly_pts], True, white, 1)
cv2.imshow(wnd_str, orig_image)
cv2.waitKey(0)

#generate contours for the rotated image                                                                           
rot_poly_pts = []
for p in pts:
    x, y = p

    blank_image = np.zeros((rows, cols, 3), np.uint8)

    blank_image[x,y] = (255, 255, 255)
    blank_image_affine = cv2.warpAffine(blank_image, M, (rows, cols), flags=cv2.INTER_CUBIC+cv2.BORDER_CONSTANT)

    rotated_y, rotated_x, rotated_z = np.unravel_index(blank_image_affine.argmax(), blank_image_affine.shape)

    rot_poly_pts.append([rotated_x, rotated_y])

rot_poly_pts = np.array(rot_poly_pts[0:-1], np.int32)
rot_poly_pts = rot_poly_pts.reshape((-1,1,2))
cv2.polylines(rot_image, [rot_poly_pts], True, white, 1)

wnd_str = 'rotated {}'.format(theta)
cv2.namedWindow(wnd_str,cv2.WINDOW_NORMAL)
cv2.moveWindow(wnd_str, 0, 0)
cv2.imshow(wnd_str, rot_image)

cv2.waitKey(0)
cv2.destroyAllWindows()


  [1]: https://i.stack.imgur.com/pyxmq.png

这是具有对齐轮廓的原始图像 在此处输入图像描述

这是旋转后的图像和叠加。在此处输入图像描述

4

1 回答 1

1

wrapAffine生成两边都有裁剪像素的图像。此过程的一个有效替代方法是转置然后翻转。

dst = cv2.flip(img.transpose(1,0,2),0)

dst  = destination image;
img  = source image;

img.transpose(col,row,channel) = 1,0,2对应于这个顺序,因为我们只想转置行和列并保持通道不变。

cv2.flip= 以抵消转置操作期间引起的镜像效应。根据您的源图像(风景或肖像)0改变标志。1

于 2017-02-23T02:14:45.397 回答