绘制高光谱3d立方体

结果如图

代码

import scipy.io as sio
import spectral
import numpy as np
import wx

data = sio.loadmat('Dataset/Indian_pines_corrected.mat')['indian_pines_corrected']
print(data.shape)
vmax = data.max(axis=(0, 1))
vmin = data.min(axis=(0, 1))
vv = (data - vmin) / (vmax - vmin)*255
vv = vv.astype(np.uint8)
app = wx.App()
spectral.settings.WX_GL_DEPTH_SIZE=16
spectral.view_cube(vv, bands=[29, 19, 9])
app.MainLoop()