Автоматизация оцифровки данных тепловой карты с использованием OpenCV - PullRequest
Новая
       7

Автоматизация оцифровки данных тепловой карты с использованием OpenCV

0 голосов
/ 02 апреля 2020

Это продолжение моего предыдущего вопроса здесь

Я пытался преобразовать данные цвета в тепловой карте в значения RGB.

исходное изображение

На приведенном ниже изображении слева представлен субплот, представленный на панели D исходного изображения. Это имеет 6 х 6 клеток (6 строк и 6 столбцов). Справа мы видим бинаризованное изображение с белым цветом, выделенным в ячейке, по которой щелкают после запуска приведенного ниже кода. Вход для запуска кода - изображение ниже. Вывод (mean = [ 27.72 26.83 144.17]) представляет собой среднее значение цвета BGR в ячейке, которая выделена белым цветом на правом изображении ниже.

enter image description here

Действительно хорошее решение в качестве ответа на мой предыдущий вопрос было указано следующее 

import cv2
import numpy as np

# print pixel value on click
def mouse_callback(event, x, y, flags, params):
    if event == cv2.EVENT_LBUTTONDOWN:
        # get specified color
        row = y
        column = x
        color = image[row, column]
        print('color = ', color)
        # calculate range
        thr = 20  # ± color range
        up_thr = color + thr
        up_thr[up_thr < color] = 255
        down_thr = color - thr
        down_thr[down_thr > color] = 0
        # find points in range
        img_thr = cv2.inRange(image, down_thr, up_thr)  # accepted range
        height, width, _ = image.shape
        left_bound = x - (x % round(width/6))
        right_bound = left_bound + round(width/6)
        up_bound = y - (y % round(height/6))
        down_bound = up_bound + round(height/6)
        img_rect = np.zeros((height, width), np.uint8)  # bounded by rectangle
        cv2.rectangle(img_rect, (left_bound, up_bound), (right_bound, down_bound), (255,255,255), -1)
        img_thr = cv2.bitwise_and(img_thr, img_rect)
        # get points around specified point
        img_spec = np.zeros((height, width), np.uint8)  # specified mask
        last_img_spec = np.copy(img_spec)
        img_spec[row, column] = 255
        kernel = np.ones((3,3), np.uint8)  # dilation structuring element
        while cv2.bitwise_xor(img_spec, last_img_spec).any():
            last_img_spec = np.copy(img_spec)
            img_spec = cv2.dilate(img_spec, kernel)
            img_spec = cv2.bitwise_and(img_spec, img_thr)
            cv2.imshow('mask', img_spec)
            cv2.waitKey(10)
        avg = cv2.mean(image, img_spec)[:3]
        print('mean = ', np.around(np.array(avg), 2))
        global avg_table
        avg_table[:, 5 - int(x / (width/6)), 5 - int(y / (height/6))] = avg
        print(avg_table)

# average value of each cell in 6x6 matrix
avg_table = np.zeros((3, 6, 6))

# create window and callback
winname = 'img'
cv2.namedWindow(winname)
cv2.setMouseCallback(winname, mouse_callback)

# read & display image
image = cv2.imread('ip.jpg', 1)
image = image[3:62, 2:118]  # crop the image to 6x6 cells
cv2.imshow(winname, image)
cv2.waitKey()  # press any key to exit
cv2.destroyAllWindows()

В настоящее время требуется нажать на каждую ячейку, представленную на изображении, чтобы получить среднее значение. Полученное таким образом среднее значение сохраняется в переменной avg_table.

. Хотелось бы узнать, есть ли способ автоматизировать этот процесс, т.е. избегать щелчков и сохранять значения непосредственно в avg_table.

Любые предложения будут действительно полезны

...