1. 程式人生 > >python+opencv圖片旋轉矩形分割

python+opencv圖片旋轉矩形分割

有時候需要對有角度的矩形框內影象從原圖片中分割出來。這裡的程式思想是,先將圖片進行矩形角度的旋轉,使有角度的矩形處於水平狀態後,根據原來座標分割圖片。
參考:http://blog.csdn.net/yjl9122/article/details/71217872
修改原來的程式:

1.旋轉函式的輸入僅為矩形的四點座標
2.角度由公式計算出來
3.矩形四點pt1,pt2,pt3,pt4由txt檔案讀入
4.在旋轉程式中還處理了順時針和逆時針及出現矩形框翻轉的問題。

# -*- coding:utf-8 -*-
import cv2
from math import *
import
numpy as np import time,math import os import re '''旋轉影象並剪裁''' def rotate( img, # 圖片 pt1, pt2, pt3, pt4 ): print pt1,pt2,pt3,pt4 withRect = math.sqrt((pt4[0] - pt1[0]) ** 2 + (pt4[1] - pt1[1]) ** 2) # 矩形框的寬度 heightRect = math.sqrt((pt1[0] - pt2[0]) ** 2 + (pt1[1] - pt2[1]) **2
) print withRect,heightRect angle = acos((pt4[0] - pt1[0]) / withRect) * (180 / math.pi) # 矩形框旋轉角度 print angle if pt4[1]>pt1[1]: print "順時針旋轉" else: print "逆時針旋轉" angle=-angle height = img.shape[0] # 原始影象高度 width = img.shape[1] # 原始影象寬度 rotateMat = cv2.getRotationMatrix2D((width / 2
, height / 2), angle, 1) # 按angle角度旋轉影象 heightNew = int(width * fabs(sin(radians(angle))) + height * fabs(cos(radians(angle)))) widthNew = int(height * fabs(sin(radians(angle))) + width * fabs(cos(radians(angle)))) rotateMat[0, 2] += (widthNew - width) / 2 rotateMat[1, 2] += (heightNew - height) / 2 imgRotation = cv2.warpAffine(img, rotateMat, (widthNew, heightNew), borderValue=(255, 255, 255)) cv2.imshow('rotateImg2', imgRotation) cv2.waitKey(0) # 旋轉後圖像的四點座標 [[pt1[0]], [pt1[1]]] = np.dot(rotateMat, np.array([[pt1[0]], [pt1[1]], [1]])) [[pt3[0]], [pt3[1]]] = np.dot(rotateMat, np.array([[pt3[0]], [pt3[1]], [1]])) [[pt2[0]], [pt2[1]]] = np.dot(rotateMat, np.array([[pt2[0]], [pt2[1]], [1]])) [[pt4[0]], [pt4[1]]] = np.dot(rotateMat, np.array([[pt4[0]], [pt4[1]], [1]])) # 處理反轉的情況 if pt2[1]>pt4[1]: pt2[1],pt4[1]=pt4[1],pt2[1] if pt1[0]>pt3[0]: pt1[0],pt3[0]=pt3[0],pt1[0] imgOut = imgRotation[int(pt2[1]):int(pt4[1]), int(pt1[0]):int(pt3[0])] cv2.imshow("imgOut", imgOut) # 裁減得到的旋轉矩形框 cv2.waitKey(0) return imgRotation # rotated image # 根據四點畫原矩形 def drawRect(img,pt1,pt2,pt3,pt4,color,lineWidth): cv2.line(img, pt1, pt2, color, lineWidth) cv2.line(img, pt2, pt3, color, lineWidth) cv2.line(img, pt3, pt4, color, lineWidth) cv2.line(img, pt1, pt4, color, lineWidth) # 讀出檔案中的座標值 def ReadTxt(directory,imageName,last): fileTxt=directory+"//rawLabel//"+imageName[:7]+last # txt檔名 getTxt=open(fileTxt, 'r') # 開啟txt檔案 lines = getTxt.readlines() length=len(lines) for i in range(0,length,4): pt2=list(map(float,lines[i].split(' ')[:2])) pt1=list(map(float,lines[i+1].split(' ')[:2])) pt4=list(map(float,lines[i+2].split(' ')[:2])) pt3=list(map(float,re.split('\n| ',lines[i+3])[:2])) # float轉int pt2=list(map(int,pt2)) pt1=list(map(int,pt1)) pt4=list(map(int,pt4)) pt3=list(map(int,pt3)) imgSrc = cv2.imread(imageName) drawRect(imgSrc, tuple(pt1),tuple(pt2),tuple(pt3),tuple(pt4), (0, 0, 255), 2) cv2.imshow("img", imgSrc) cv2.waitKey(0) rotate(imgSrc,pt1,pt2,pt3,pt4) if __name__=="__main__": directory = "G://grasp//grapCode//trainImage//jpg//4" last = 'cneg.txt' imageName="pcd0247r.png" ReadTxt(directory,imageName,last)

原帶角度的矩形框:
這裡寫圖片描述

旋轉矩形框:
這裡寫圖片描述

分割:
這裡寫圖片描述