Pracuję nad kodem statystycznym, który generuje różne wykresy za pomocą Matlaba. Typy wykresów wahają się od prostych wykresów kołowych i słupkowych po trójwymiarowe siatki histogramów.

Teraz chcielibyśmy mieć ładne GUI z oprogramowaniem. Mamy prototyp GUI Matlaba, ale GUI Matlaba ma wiele problemów, więc chcielibyśmy przejść do bardziej niezawodnego GUI. Moim najlepszym rozwiązaniem wydaje się PySide + matplotlib, ale jak dotąd nie znalazłem sposobu na wykreślenie siatki 3D. Kod matlab używa contourslice. Nie ma podobnego połączenia w matplotlib. Czy więc ktoś może mi pomóc dowiedzieć się, w jaki sposób uzyskać taki wykres z matplotlibem? Jak dotąd moim jedynym pomysłem jest wykreślenie 6 powierzchni w celu utworzenia sześcianu.

W rzeczywistości mile widziane są również zalecenia dotyczące innych kombinacji bibliotek GUI / Graph. Kod stats to C ++, więc Python to tylko jedna z wielu opcji. Sądząc po niektórych odpowiedziach, które widziałem na StackOverflow, matplotlib może być niedopuszczalnie powolny w 3D. Może R będzie lepszy?

Oto kod matlab:

clf
xlo = -1.800000e+01; 
xhi = 1.000000e+01; 
ylo = 1.000000e+01; 
yhi = 3.000000e+01; 
zlo = -1.000000e+03; 
zhi = 1.000000e+03; 
X=zeros(16,16,16);
Y=zeros(16,16,16);
Z=zeros(16,16,16);
V=zeros(16,16,16);
% fill in X, Y, Z, and V  huge amount of text
xt = [-1.800000e+01:2.800000e-01:1.000000e+01];
yt = [1.000000e+01:2.000000e-01:3.000000e+01];
zt = [-1.000000e+03:2.000000e+01:1.000000e+03];
isoval = -1.428280e+01;
h = patch(isosurface(X,Y,Z,V,isoval),... 
'FaceColor', 'blue', ... 
'EdgeColor', 'none', ... 
'AmbientStrength', 0.2, ... 
'SpecularStrength', 0.7, ... 
'DiffuseStrength', 0.4);
isonormals(X,Y,Z,V,h);
patch(isocaps(X,Y,Z,V,isoval), ...
'FaceColor', 'interp', ... 
'EdgeColor', 'none'); 
axis([xlo xhi ylo yhi zlo zhi])
daspect([2.800000e+01,2.000000e+01,2.000000e+03])
set(gca,'linewidth',2)
set(gca,'fontweight','bold')
set(gca,'fontsize',12)
grid on
box on
colormap('default'); colorbar

view(3) 
set(gcf,'Renderer','zbuffer')
lighting phong
cin = 'n';
if (cin == 'y')
xin = questdlg('Axis to slide through ?', 'Axis', 'X', 'Y', 'Z', 'X');
xin = lower(xin);
for i = 1 : 101
if gcf ~= plotFigure
return
end
if (xin == 'y')
h = contourslice(X,Y,Z,V,xt(i),[],[],101);
elseif (xin == 'x')
h = contourslice(X,Y,Z,V,[],yt(i),[],101);
elseif (xin == 'z')
h = contourslice(X,Y,Z,V,[],[],zt(i),101);
end
axis([-1.8000e+01  1.0000e+01  1.0000e+01  3.0000e+01 -1.0000e+03  1.0000e+03 -8.6774e+01  4.2066e+02])
set(gca,'linewidth',2)
set(gca,'fontweight','bold')
set(gca,'fontsize',12)
grid on
box on
view(3)
set(h, 'Linewidth', 10)
colorbar
pause(0.3)
if (i < 101)
clf
end
end
end