MatlabのようなJulia3D表面プロット、対数スケールのカラーマップ

Aug 16 2020

次の図のように3D表面プロットをプロットしようとしています(出典:ウィキペディア):

これが私の試みjuliaです;

x  = [-2:0.05:2;]
y = [-1:0.05:3;]
z = (1 .-x').^2 .+ 100 .*(y.-x'.^2).^2

minZ = minimum(z[:]);  
maxZ = maximum(z[:]);

c =  minZ .+ (maxZ-minZ).*log.(1 .+z .- minZ) ./ log(1+maxZ-minZ)
Plots.plot(x,y,z,st=:surface,color=cgrad(:jet,c),
xlabel = "x",ylabel="y",zlabel="f(x,y)")

これが私の質問です:

  1. 対数目盛が機能しない場合でも、ジュリアプロットのカラーマップをMatlabのように表示するにはどうすればよいですか。
  2. Matlabプロットの方位角と標高のビュー(カメラ)は(-30,30)ですが、juliaプロットカメラオプションはまったく機能しません。camera = (-30,30)。負の値はcameraオプションでは機能しないようです。ビュー(xおよびyオプション)をmatlabのように表示するにはどうすればよいですか?

参考のためにMatlabcodel:

clear all;
close all;
hfig = figure(1);

s = 0.05;
X = [-2 : s : 2+s];
Y = [-1 : s : 3+s];
[X, Y] = meshgrid(X, Y);
Z = (1-X).^2 + 100*(Y-X.^2).^2;

% Use log scale of Z for a batter usage of the color spectrum  
minZ = min(Z(:));  
maxZ = max(Z(:));
C = minZ + (maxZ-minZ).*log(1+Z-minZ)./log(1+maxZ-minZ); 


colormap(jet);
surf(X, Y, Z, C, 'EdgeColor', 'none', 'LineStyle', 'none');
 
axis([-2, 2, -1, 3, 0, 2500]);
xlabel('x', 'fontsize', 18);
ylabel('y', 'fontsize', 18);
zlabel('f', 'fontsize', 18);

ありがとうございました

回答

1 PrzemyslawSzufel Aug 16 2020 at 01:36

多分これは役に立つでしょう?

pyplot()
Plots.plot(x,y,z,st=:surface,color=Colors.diverging_palette(20,300,40, logscale=true,  wcolor=colorant"red",dcolor1=colorant"green", dcolor2=colorant"red", b=0.1,d1=1.0,d2=1.0),
       xlabel = "x",ylabel="y",zlabel="f(x,y)",camera=(-30,30))