1.H = cv2.getPerspectiveTransform(rect, transform_axes) 获得投射变化后的H矩阵

参数说明:rect表示原始的位置左上,右上,右下,左下, transform_axes表示变换后四个角的位置

2.cv2.warpPerspective(gray, H, (width, height)) 根据H获得变化后的图像

参数说明: gray表示输入的灰度图像, H表示变化矩阵,(width, height)表示变换后的图像大小
3. cv2.approxPloyDP(contour, length*0.02, True) # 计算轮廓的近似值

参数说明:contour表示输入的轮廓值, length表示轮廓的大小, True表示是否闭合

4. cv2.threshold(wrap, 0, 255, cv2.THRESH_BINARY_INV | cv2.THRESH_AUTO) # 进行图像的二值化操作
参数说明: wrap表示输入图片,0表示最小值,255表示最大值,THRESH_BINARY_INV | cv2.THRESH_AUTO表示方法,即大于最小阈值的为0, 小于最小阈值的为255

5. cv2.draw_contours(image, [contour], -1, 255, -1) # 对轮廓进行画图

参数说明:image表示输入图片. [contour] 表示单个轮廓, -1表示轮廓的index, 255表示颜色,-1表示对轮廓进行填充

6. cv2.bitwise(image, image, mask=mask)  # 与判断,如果两个都大于1才为1,否则都为0 

参数说明:image表示输入图片,image表示输出图片, mask=mask表示掩膜

7. cv2.countNoneZeros(mask) # 用于计算图像中不等于0的像素点

参数说明:mask表示经过图像与掩膜的与操作以后的图像

8.cv2.putText(image, text, org, font, font_scale, color, ticking) # 在图像上加上文字

参数说明:image表示输入图片,text表示文字,org表示文字位置,font表示字体的格式,font_scale表示字体的胖瘦,color表示颜色,ticking表示粗细

答题卡识别判卷, 需要做的就是将画圈圈的地方识别出来。

    

              原始图像                             涂了铅笔的比没涂白色部分要多,使用掩膜进行逐个遍历

第一部分:对图片进行透视变化,使得图片中的答题卡可以凸显出来

第一步:读入图片

第二步:使用cv2.gray() 进行灰度值的变化
第三步:使用cv2.GuassianBlur() 进行高斯滤波操作

第四步:使用cv2.Canny找出图片的边缘信息,为了下一步找出轮廓值做准备

第五步:使用cv2.findCountor() 找出图片的轮廓值,使用key = cv2.ContourArea对轮廓值进行排序

第六步:对轮廓值进行遍历,使用cv2.approxPolyDP(c, 0.02*lenght, True) 获得轮廓的近似值,如果近似值的维度为4,即为最外层答题卡的四个角的维度,则跳出循环

第七步:构造函数,对求得的轮廓值,根据其位置信息,将其变成左上,右上,右下,左上的形式的rect

第八步:计算width,height,计算出变换后的四个角的位置(transform_axes)

第九步:使用H = cv2.getPerspectiveTransform(rect, transform_axes) rect为变化前的位置,transform_axes为变化后的位置,以x,y表示,获得变化矩阵H 

第十步:使用cv2.warpPerspective(thresh, H, (width, height)) 获得最终的变化结果

 

第二步:根据正确答案和图卡的答案来判断正确率

第一步:根据上述获得的wrapped,进行二值化操作, 使用cv2.COLOR_BINARY_INV将图了颜色的地方变成白色

第二步:使用cv2.findContours进行轮廓检测

第三步:对获得的轮廓,使用外接矩形cv2.boundingRect的(x, y, w, h)对轮廓进行筛选,找出答案的轮廓
第四步:对轮廓进行从上到下的排序

第五步:分别对每一行进行遍历,np.arange(0, len(question_cnts, 5)), 遍历每一行的轮廓值,根据图像构造掩膜mask图像,使用cv2.draw_contours在掩膜上画出对应位置的白色轮廓, 使用cv2.bitwise_and 获得原始图像当前位置等于1的点,使用cv2.countNoneZero计算出不是零的点,我们发现涂了颜色的要比没涂颜色的要大

第六步:计算出颜色值最大的位置和Nonezeros的值,correct+= 1

第七步:使用cv2.draw_contours在正确位置上进行画图操作

第八步:使用cv2.putText将正确率打印在图像上

 

代码:

# 第一部分代码:

import cv2
import numpy as np
import matplotlib.pyplot as plt


# 正确答案
ANSWER_KEY = {0: 1, 1: 4, 2: 0, 3: 3, 4: 1}
def cv_show(name, image):
    cv2.imshow(name, image)
    cv2.waitKey(0)
    cv2.destroyAllWindows()

def four_point_transform(image, dotCnt):
    # 第七步:按照轮廓近似的位置,将结果变成左上,右上,右下,左下的位置
    rect = sort_dotCnt(dotCnt)
    (tl, tr, br, bl) = rect
    # 第八步:计算轮廓的宽和长,构造出变化后四个角的位置
    widthA = np.sqrt(((tl[0]-tr[0])**2) + ((tl[1]-tr[1])**2))
    widthB = np.sqrt(((bl[0]-br[0])**2) + ((bl[1] - br[1])**2))
    widthMax = max(int(widthB), int(widthA))

    heightA = np.sqrt(((tl[0]-bl[0])**2) + ((tl[1]-bl[1])**2))
    heightB = np.sqrt(((tr[0]-br[0])**2) + ((tr[1] - br[1])**2))
    heightMax = max(int(heightA), int(heightB))
    # 四个角的位置信息
    transform_axis = np.array([
        [0, 0],
        [widthMax-1, 0],
        [widthMax-1, heightMax-1],
        [0, heightMax-1],
    ], dtype='float32')

    # 第九步:使用cv2.getPerspectiveTransform构造变化矩阵H
    H = cv2.getPerspectiveTransform(rect, transform_axis)
    # 第十步:使用cv2.warpPerspective()获得变化后的矩阵
    wrapped = cv2.warpPerspective(image, H, (widthMax, heightMax))

    cv_show('wrapped', wrapped)

    return wrapped






def sort_dotCnt(kps):

    rect = np.zeros((4, 2), dtype='float32')
    s = kps.sum(axis=1)
    # 找出左上和右下
    rect[0] = kps[np.argmin(s)]
    rect[2] = kps[np.argmax(s)]
    # 找出右上和左下
    diff = np.diff(kps, axis=1)
    rect[1] = kps[np.argmin(diff)]
    rect[3] = kps[np.argmax(diff)]

    return rect

def sorted_contours(cnt, model='left-to-right'):

    if model == 'top-to-bottom':
        cnt = sorted(cnt, key=lambda x:cv2.boundingRect(x)[1])

    elif model == 'left-to-right':
        cnt = sorted(cnt, key=lambda x:cv2.boundingRect(x)[0])

    return cnt



# 第一部分:使用投射变化进行图像的转换, 将答题卡放在图片的正中央
# 第一步:读入图片
image = cv2.imread('images/test_01.png')
# 第二步:进行灰度化操作
gray = cv2.cvtColor(image, cv2.COLOR_BGR2GRAY)
cv_show('gray', gray)

# 第三步: 进行高斯变化
Guass = cv2.GaussianBlur(gray, (5, 5), 0)
# 第四步:使用cv2.Canny 找出图像的边缘部分
edges = cv2.Canny(Guass, 50, 200)
cv_show('edges', edges)
# 第五步:使用cv2.findCountours找出图像的轮廓, 并对轮廓进行排序
cnt = cv2.findContours(edges, cv2.RETR_EXTERNAL, cv2.CHAIN_APPROX_SIMPLE)[1]
ret = cv2.drawContours(image.copy(), cnt, -1, (0, 0, 255), 2)
cv_show('ret', ret)

dotCnt = None

if len(cnt) > 0:

    cnt = sorted(cnt, key=cv2.contourArea, reverse=True)
    # 第六步: 遍历轮廓,获得cv2.approxPolyDP获得轮廓的近似值, 如果approx=4,即为填图卡的轮廓近似
    for c in cnt:

        length = cv2.arcLength(c, True)
        # 获得轮廓的近似值
        approx = cv2.approxPolyDP(c, 0.02*length, True)
        # 如果轮廓的维度是4,即由4个点组成,就是一个方框,那就是外面的方框
        if len(approx) == 4:

            dotCnt = approx
            break


wrapped = four_point_transform(gray, dotCnt.reshape(4, 2))
cv_show('wrapped', wrapped)

            

        原始图片                            灰度变化后图片         边缘检测的图片                 检测到的轮廓图            最终透射变化结果

 

# 第二部分代码:

# 第二部分:对变化后的wrapped进行操作,构造mask根据有无填涂的特性,进行位置的计算

#第一步:对图像进行二值化操作
thresh = cv2.threshold(wrapped, 0, 255, cv2.THRESH_BINARY_INV | cv2.THRESH_OTSU)[1]
cv_show('thresh', thresh)

# 第二步:对图像进行轮廓检测
cnts = cv2.findContours(thresh, cv2.RETR_EXTERNAL, cv2.CHAIN_APPROX_SIMPLE)[1]

# 第三步:使用cv2.boudingRect将答案信息的轮廓进行筛选

questionCnts = []
for c in cnts:

    (x, y, w, h) = cv2.boundingRect(c)
    arc = w / float(h)

    # 根据实际情况找出合适的轮廓
    if w > 20 and h > 20 and arc >= 0.9 and arc<=1.1:
        questionCnts.append(c)

# 第四步:将轮廓进行从上到下的排序
questionCnts = sorted_contours(questionCnts, model='top-to-bottom')

# 第五步:对每一行的轮廓进行遍历,构造掩码,对每一行答案的轮廓进行遍历,使用cv2.countNoneZeros()计算出涂了颜色的答案,并计算最终的结果
correct = 0
for j, qC in enumerate(np.arange(0, len(questionCnts), 5)):

    questionCnt = questionCnts[qC:qC+5]
    # 进行排序
    questionCnt = sorted_contours(questionCnt)
    questionCnt = np.array(questionCnt)
    # 每个问题的答案
    bubbled = None
    for k, c in enumerate(questionCnt):
        # 根据轮廓值构造掩膜, 因为有涂答案的地方是白色的
        mask = np.zeros(thresh.shape, np.uint8)
        # 使用cv2.drawContours画出轮廓值
        cv2.drawContours(mask, [c], -1, 255, -1)
        cv_show('mask', mask)
        mask = cv2.bitwise_and(thresh, thresh, mask=mask)
        # 第六步:使用cv2.countNoneZero计算出白色最多的位置
        Nonezero = cv2.countNonZero(mask)
        if bubbled == None or Nonezero > bubbled[0]:
            bubbled = (Nonezero, k)

    question_answer = int(bubbled[1])
    answer_True = ANSWER_KEY[j]

    # 如果结果和正确结果相同,就+1
    if question_answer == answer_True:
        correct += 1
    # 第七步:对正确答案的轮廓进行画图操作
    cv2.drawContours(wrapped, [questionCnt[answer_True]], -1, 0, 2)

cv_show('wrapped', wrapped)

score = (correct / 5) * 100
# 第八步:使用cv2.putText将正确率打印在图上
cv2.putText(wrapped, 'correct_score:%d' %score, (10, 30), cv2.FONT_HERSHEY_SIMPLEX,
            0.8, 0, 2)
cv2.imshow('orginal', image)
cv2.imshow('wrapped', wrapped)

cv2.waitKey(0)
cv2.destroyAllWindows()

               

   Binary_Inv二值化的结果         draw_contours画出掩膜轮廓                      最终的结果图

 

posted on 2019-02-26 18:35  python我的最爱  阅读(3405)  评论(0编辑  收藏  举报