v3d_show.m

MATLAB三维显示地震数据, 开放小工具。 这个开放MABLAB三维绘图小工具, 非常好用, 可进行SLICE, CONTOUR, 等交互功能.

function v3d_show(x,y,z,v,options)
% Standard view after loading of data

% Adaptation of function "v3d_show" by Robert Barsch; the original
% version is available at The Matlab Central File Exchange, File ID 2255.
% http://www.mathworks.com/matlabcentral/fileexchange/loadFile.do?objectId=2255&objectType=file
%
% Modified by E. Rietsch: October 15, 2006

global V3D_HANDLES

figure_handle=V3D_HANDLES.figure_handle;

set(figure_handle,'Tag','V3D:FIGURE');
figure(figure_handle)

%	Delete all objects of current axes
set(figure_handle,'HandleVisibility','on')
axis_handle=cla;
set(axis_handle,'parent',figure_handle,'Tag','V3D:AXES');
setappdata(figure_handle,'AxisHandle',axis_handle)

%	Proportional axes
if options.equal_axes
   axis equal;
end

% 3D-View
myview(3);

%	Axis labels
xlabel(info2label(options.xinfo));
ylabel(info2label(options.yinfo));
zlabel(info2label(options.zinfo));

   %	Set direction of the z-axis
set(axis_handle,'ZDir',options.zdir,'Layer','top','Box','on')

% Limits setzen
set(axis_handle,'xlim',[min(x) max(x)],'ylim',[min(y) max(y)],'zlim',[min(z) max(z)]);  

% Userdaten des V3D-Fensters setzen
% x,y,z,v 
ud.x=x; % x Vektor
ud.y=y; 
ud.z=z;
ud.v=v; % v Data matrix
%ud.options=options;

userdata.xmin=min(x);
userdata.xmax=max(x);
userdata.ymin=min(y);
userdata.ymax=max(y);
userdata.zmin=min(z);
userdata.zmax=max(z);
setappdata(figure_handle,'userdata',userdata)
setappdata(figure_handle,'options',options)

% Colormap setzen
ud.cmap.name='jet'; % Colormap jet
ud.cmap.reverse=0; % Color direction not reverse
ud.cmap.brighten=0; % normal brightness
ud.cmap.log=0; % nicht logarithmic

% Alphadaten
ud.alphadata=1; % Alpha data available
ud.alphav=ones(length(y),length(x),length(z)); % Alphadaten everywhere = 1.0 (completely visible)

% Userdaten setzen
set(figure_handle,'userdata',ud);

% Auslesen von Minimum und Maximum der Daten
cmin=min(v(:));
cmax=max(v(:));

% Testen ob Maximum ungleich Minimum ansonsten Fehler
if (cmax == cmin) 
    error('Global minimum of volume ist equal to global maximum!')
 end

%	Parameters
method_list={'linear','nearest','cubic'};
lighting_list={'none','flat','gouraud','phong'};

% Userdata des Objektes setzen -> ben鰐igt f黵s auslesen in Slicemen