python - 为 3D 体积绘制 seaborn 热图动画时出错
问题描述
试图可视化两个体积img_3D和mask_3D之间的互相关,使用Seaborn heatmap和来自 Matplotlib 的动画将 3D 互相关结果可视化为 2D 图像的渐进动画,但我遇到了一个错误,请你帮忙告诉我如何摆脱这个错误,并正确地可视化热图?
提前致谢。
Traceback (most recent call last):
File "C:\Users\User\AppData\Local\Programs\Python\Python37\lib\tkinter\__init__.py", line 1705, in __call__
return self.func(*args)
File "C:\Users\User\AppData\Local\Programs\Python\Python37\lib\site-packages\matplotlib\backends\_backend_tk.py", line 259, in resize
self.draw()
File "C:\Users\User\AppData\Local\Programs\Python\Python37\lib\site-packages\matplotlib\backends\backend_tkagg.py", line 9, in draw
super(FigureCanvasTkAgg, self).draw()
File "C:\Users\User\AppData\Local\Programs\Python\Python37\lib\site-packages\matplotlib\backends\backend_agg.py", line 392, in draw
else nullcontext()):
File "C:\Users\User\AppData\Local\Programs\Python\Python37\lib\contextlib.py", line 112, in __enter__
return next(self.gen)
File "C:\Users\User\AppData\Local\Programs\Python\Python37\lib\site-packages\matplotlib\backend_bases.py", line 2788, in _wait_cursor_for_draw_cm
self.set_cursor(self._lastCursor)
File "C:\Users\User\AppData\Local\Programs\Python\Python37\lib\site-packages\matplotlib\backends\_backend_tk.py", line 544, in set_cursor
window.configure(cursor=cursord[cursor])
File "C:\Users\User\AppData\Local\Programs\Python\Python37\lib\tkinter\__init__.py", line 1485, in configure
return self._configure('configure', cnf, kw)
File "C:\Users\User\AppData\Local\Programs\Python\Python37\lib\tkinter\__init__.py", line 1476, in _configure
self.tk.call(_flatten((self._w, cmd)) + self._options(cnf))
_tkinter.TclError: invalid command name "."
使用的代码是:
# Import Libraries
#====================================
import numpy as np
np.random.seed(0)
import matplotlib.pyplot as plt
from matplotlib.animation import FuncAnimation
import nibabel as nib
from scipy.signal import correlate
import seaborn as sns
sns.set()
#===================================
img = np.load('img.npy')
act = np.load('act.npy')
# Mode : 'full', 'valid', 'same'
result = correlate(img, act,mode='same')
print(img.shape, act.shape, result.shape)
def updatefig(sl):
for sl in range(result.shape[2]):
print(sl,' / ',result.shape[2])
sns.heatmap(result[...,sl],cbar=False)
ax.set_title("frame {}".format(sl))
# Note that using time.sleep does *not* work here!
plt.pause(0.1)
fig, ax = plt.subplots()
ani = FuncAnimation(fig, updatefig, frames=range(result.shape[2]), interval=5, blit=True)
plt.show()
解决方案
检查此代码:
import numpy as np
np.random.seed(0)
import matplotlib.pyplot as plt
from matplotlib.animation import FuncAnimation
from scipy.signal import correlate
import seaborn as sns
sns.set()
img = np.load('img.npy')
act = np.load('act.npy')
result = correlate(img, act, mode = 'same')
def updatefig(sl):
ax.cla()
print(sl + 1, ' / ', result.shape[2])
sns.heatmap(result[..., sl], cbar = False)
ax.set_title("frame {}".format(sl + 1))
ax.axis('off')
fig, ax = plt.subplots()
ani = FuncAnimation(fig, updatefig, frames = result.shape[2], interval = 5)
plt.show()
这给了我这个动画(我将下面报告的动画减半以将文件大小减少到 2 MB 以下,上面的代码重现了所有 40 帧):
编辑
为了给热图添加一个固定的颜色条,检查这个代码:
import numpy as np
np.random.seed(0)
import matplotlib.pyplot as plt
from matplotlib.animation import FuncAnimation
from scipy.signal import correlate
import seaborn as sns
sns.set()
img = np.load('img.npy')
act = np.load('act.npy')
result = correlate(img, act, mode = 'same')
def updatefig(sl):
ax.cla()
print(sl + 1, ' / ', result.shape[2])
sns.heatmap(result[..., sl],
ax = ax,
cbar = True,
cbar_ax = cbar_ax,
vmin = result.min(),
vmax = result.max())
ax.set_title("frame {}".format(sl + 1))
ax.axis('off')
grid_kws = {'width_ratios': (0.9, 0.05), 'wspace': 0.2}
fig, (ax, cbar_ax) = plt.subplots(1, 2, gridspec_kw = grid_kws, figsize = (10, 8))
ani = FuncAnimation(fig, updatefig, frames = result.shape[2], interval = 5)
plt.show()
产生这个动画(剪成前一个):
推荐阅读
- windows - 获取 USB_BANDWIDTH_INFO 结构的问题
- github-api - PyGithub - 如何获取特定分支的所有提交列表?
- flutter - 如何在 Flutter 中仅将边框/阴影添加到底部?
- python - 在 Python 中:如果拆分字符串与列值匹配,如何比较并生成新列?
- linux - Yocto - 未应用补丁
- c# - 最大长度验证阻止我在 mvc6 中编写超过最大长度的内容
- python - 在 FTP 发送/接收软件的 while 循环中尝试/排除块
- aframe - AR.js 场景中可点击对象的位置错误
- python - 为什么我不能总结熊猫的价值?
- mysql - Sequelize 模型导出函数