小编Gre*_*nes的帖子

我试图覆盖两个图像.第一个是512x512 numpy阵列(来自CT图像).第二个也是512x512 numpy数组,但我只对值大于0的像素(功能图像)感兴趣.

为此,我试图创建一个蒙版数组.

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

# Both images are loaded from a dicom. Both are numpy arrays of (512,512) 
Image1 = readimage(path)
Image2 = readimage(path)
# Create image 2 mask
mask = ma.masked_where(Image2>0, Image2)
Image2_mask = ma.masked_array(Image2,mask)

# Plot images
plt.figure(dpi=300)
y, x = np.mgrid[1:513,1:513]
plt.axes().set_aspect('equal', 'datalim')
plt.set_cmap(plt.gray())
plt.pcolormesh(x, y, Image1,cmap='gray')
plt.pcolormesh(x, y, Image2_mask,cmap='jet')
plt.axis([x.min(), x.max(), y.min(), y.max()])
plt.colorbar()
plt.show()

此代码不显示任何叠加.我做错了什么？有没有正确的方法？我来自Matlab环境,我对phyton很新手.

谢谢

我有一个512x512x512 numpy阵列.是否有任何有效的方法来执行均值过滤器,其中每个数组值都被所有3x3x3本地值替换？

我们正在寻找与scipy.ndimage.filters.median_filter中位数相似的但与平均值相关的东西.