📄 plot3d.html
字号:
<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN" "http://www.w3.org/TR/REC-html40/loose.dtd"><html><head> <title>Description of plot3d</title> <meta name="keywords" content="plot3d"> <meta name="description" content="PLOT3D - produces an image of a 3D object"> <meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1"> <meta name="generator" content="m2html © 2003 Guillaume Flandin"> <meta name="robots" content="index, follow"> <link type="text/css" rel="stylesheet" href="../m2html.css"></head><body><a name="_top"></a><div><a href="../index.html">Home</a> > <a href="index.html">mri_toolbox</a> > plot3d.m</div><!--<table width="100%"><tr><td align="left"><a href="../index.html"><img alt="<" border="0" src="../left.png"> Master index</a></td><td align="right"><a href="index.html">Index for mri_toolbox <img alt=">" border="0" src="../right.png"></a></td></tr></table>--><h1>plot3d</h1><h2><a name="_name"></a>PURPOSE <a href="#_top"><img alt="^" border="0" src="../up.png"></a></h2><div class="box"><strong>PLOT3D - produces an image of a 3D object</strong></div><h2><a name="_synopsis"></a>SYNOPSIS <a href="#_top"><img alt="^" border="0" src="../up.png"></a></h2><div class="box"><strong>function [b,d]=plot3d(a,alfa,beta); </strong></div><h2><a name="_description"></a>DESCRIPTION <a href="#_top"><img alt="^" border="0" src="../up.png"></a></h2><div class="fragment"><pre class="comment"> PLOT3D - produces an image of a 3D object
To use: [b,d]=plot3d(a,alfa,beta);
This function produces an image of a 3D object
defined by matrix a(l,m,n) in terms of voxels.
The image is a view after rotating the object
by angles alfa and beta (in degrees).
b is the image and d is its distance to the viewer matrix
The first figure depicts the object using only its gray level values
The second image depicts the object using some lighting effect
rotate3d may be used for reorientation but the obtained image is planar
Kindly be patient the prog. is very slow!
July-9-1998 Dr. H Azhari
Example:
x=zeros(40,40,40);
x(10:30,10:30,10:30)=ones(21,21,21);
x(15:25,15:25,:)=zeros(11,11,40);
plot3d(x,30,20);
colormap gray</pre></div><!-- crossreference --><h2><a name="_cross"></a>CROSS-REFERENCE INFORMATION <a href="#_top"><img alt="^" border="0" src="../up.png"></a></h2>This function calls:<ul style="list-style-image:url(../matlabicon.gif)"></ul>This function is called by:<ul style="list-style-image:url(../matlabicon.gif)"></ul><!-- crossreference --><h2><a name="_source"></a>SOURCE CODE <a href="#_top"><img alt="^" border="0" src="../up.png"></a></h2><div class="fragment"><pre>0001 <a name="_sub0" href="#_subfunctions" class="code">function [b,d]=plot3d(a,alfa,beta);</a>0002 <span class="comment">% PLOT3D - produces an image of a 3D object</span>0003 <span class="comment">%</span>0004 <span class="comment">% To use: [b,d]=plot3d(a,alfa,beta);</span>0005 <span class="comment">%</span>0006 <span class="comment">% This function produces an image of a 3D object</span>0007 <span class="comment">% defined by matrix a(l,m,n) in terms of voxels.</span>0008 <span class="comment">% The image is a view after rotating the object</span>0009 <span class="comment">% by angles alfa and beta (in degrees).</span>0010 <span class="comment">% b is the image and d is its distance to the viewer matrix</span>0011 <span class="comment">% The first figure depicts the object using only its gray level values</span>0012 <span class="comment">% The second image depicts the object using some lighting effect</span>0013 <span class="comment">% rotate3d may be used for reorientation but the obtained image is planar</span>0014 <span class="comment">% Kindly be patient the prog. is very slow!</span>0015 <span class="comment">% July-9-1998 Dr. H Azhari</span>0016 <span class="comment">%</span>0017 <span class="comment">% Example:</span>0018 <span class="comment">% x=zeros(40,40,40);</span>0019 <span class="comment">% x(10:30,10:30,10:30)=ones(21,21,21);</span>0020 <span class="comment">% x(15:25,15:25,:)=zeros(11,11,40);</span>0021 <span class="comment">% plot3d(x,30,20);</span>0022 <span class="comment">% colormap gray</span>0023 0024 0025 time0=clock;0026 <span class="comment">%-----------------</span>0027 <span class="comment">% Intial values</span>0028 [l,m,n]=size(a);0029 M=max(size(a));0030 M1=round(sqrt(2)*M)+2;0031 b=zeros(M1,M1);0032 d=-M1*ones(M1,M1);0033 alfa=alfa*pi/180;0034 beta=beta*pi/180;0035 T=zeros(3,3); <span class="comment">% Rotation matrix</span>0036 T(1,1)=cos(alfa);0037 T(1,2)=sin(alfa);0038 T(1,3)=0;0039 T(2,1)=-cos(beta)*sin(alfa);0040 T(2,2)=cos(beta)*cos(alfa);0041 T(2,3)=sin(beta);0042 T(3,1)=sin(beta)*sin(alfa);0043 T(3,2)=-sin(beta)*cos(alfa);0044 T(3,3)=cos(beta);0045 xyz=[1,1,1]';0046 <span class="comment">%-------------------------</span>0047 <span class="keyword">for</span> i=1:l,0048 <span class="keyword">for</span> j=1:m,0049 <span class="keyword">for</span> k=1:n,0050 xyz=round(T*[i-l/2,j-m/2,k-n/2]');0051 x1=round(M1/2+xyz(1)+1);0052 y1=round(M1/2+xyz(2)+1);0053 z1=round(M1/2+xyz(3)+1);0054 <span class="keyword">if</span> a(i,j,k)>0,0055 <span class="keyword">if</span> d(x1,y1)<z1,0056 b(x1,y1)=a(i,j,k);0057 d(x1,y1)=z1;0058 <span class="keyword">end</span>0059 <span class="keyword">end</span>0060 0061 0062 <span class="keyword">end</span>0063 <span class="keyword">end</span>0064 <span class="keyword">end</span>0065 0066 imagesc(b)0067 figure0068 [fx,fy]=gradient(d);0069 h=(b.*(d.^1.5+mean(mean(d/2))*(cos((atan(fy)+atan(fx))/2)).^2));0070 pcolor(flipud(h));shading <span class="string">'interp'</span> ;0071 <span class="comment">%pcolor(flipud((b.*d).^1.5));</span>0072 <span class="comment">%shading 'interp'</span>0073 0074 <span class="comment">%------------------------------</span>0075 Elapsed_Time=etime(clock,time0)</pre></div><hr><address>Generated on Fri 21-May-2004 12:38:21 by <strong><a href="http://www.artefact.tk/software/matlab/m2html/">m2html</a></strong> © 2003</address></body></html>⌨️ 快捷键说明
复制代码
Ctrl + C
搜索代码
Ctrl + F
全屏模式
F11
切换主题
Ctrl + Shift + D
显示快捷键
?
增大字号
Ctrl + =
减小字号
Ctrl + -