Порог Otsu: AttributeError: у объекта 'AxesSubplot' нет атрибута 'ravel' - PullRequest
Новая
       51

Порог Otsu: AttributeError: у объекта 'AxesSubplot' нет атрибута 'ravel'

1 голос
/ 12 марта 2019

Я загрузил отличные файлы (они также были конвертированы из .pack CT сканирований).Моя цель состоит в том, чтобы использовать алгоритм Treashold Otsu, чтобы замаскировать его от фона и сравнить два изображения.Когда я пытаюсь построить график, я получаю сообщение об ошибке 

AttributeError: 'AxesSubplot' object has no attribute 'ravel'

Ниже приведен код и прикрепленный скриншот.

import SimpleITK as sitk
import matplotlib.pyplot as plt
import numpy as np
from skimage.filters import threshold_otsu

#THRESHOLD OTSU
img = sitk.GetArrayFromImage(sitk.ReadImage("\\\\x.x.x.x/users/ddff/python/nifts/prr_ipsi.nii"))
print(img.shape)
thresh = threshold_otsu(img.flatten())
#thresh = thresh.reshape(img.shape)
binary = img <= thresh

#I can plot this image slice fine
plt.imshow(img[20,:,:])

fig, axes = plt.subplots(ncols=1)
ax = axes.ravel()
ax[0] = plt.subplot(1, 3, 1)
ax[1] = plt.subplot(1, 3, 2)
ax[2] = plt.subplot(1, 3, 3, sharex=ax[0], sharey=ax[0])

ax[0].imshow(img[20,:,:], cmap=plt.cm.gray)
ax[0].set_title('Original Breast Delineation')
ax[0].axis('off')

ax[1].hist(thresh, bins=256)
ax[1].set_title('Histogram ')
ax[1].axvline(thresh, color='r')

ax[2].imshow(binary[20,:,:], cmap=plt.cm.gray)
ax[2].set_title('Thresholded')
ax[2].axis('off')

plt.show()[enter image description here][1]

1 Ответ

2 голосов
/ 12 марта 2019

axes - это всего лишь одна цифра с 1 столбцом, поэтому нет ничего для ravel или flatten. Это будет работать, если у вас есть несколько подзаговоров. Тем не менее, вы можете сделать следующее без ravel, если у вас есть только одна строка или один столбец. 

fig, ax = plt.subplots(ncols=3, sharex=True, sharey=True)

ax[0].imshow(img[20,:,:], cmap=plt.cm.gray)
ax[0].set_title('Original Breast Delineation')
ax[0].axis('off')

ax[1].hist(thresh, bins=256)
ax[1].set_title('Histogram ')
ax[1].axvline(thresh, color='r')

ax[2].imshow(binary[20,:,:], cmap=plt.cm.gray)
ax[2].set_title('Thresholded')
ax[2].axis('off')

Если вам нужна двумерная матрица экземпляров подзаговоров, вы можете воспользоваться предложением Thomas Kühn. 

fig, ax = plt.subplots(ncols=3, sharex=True, sharey=True, squeeze=False)

и затем вы можете получить доступ к участкам как 

ax[0][0].imshow()
ax[0][1].imshow()
......
...