📄 dtiguicode.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 DTIguicode</title> <meta name="keywords" content="DTIguicode"> <meta name="description" content=""> <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> > DTIguicode.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>DTIguicode</h1><h2><a name="_name"></a>PURPOSE <a href="#_top"><img alt="^" border="0" src="../up.png"></a></h2><div class="box"><strong></strong></div><h2><a name="_synopsis"></a>SYNOPSIS <a href="#_top"><img alt="^" border="0" src="../up.png"></a></h2><div class="box"><strong>This is a script file. </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"></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)"><li><a href="DTIvol.html" class="code" title="function varargout = DTIvol(varargin)">DTIvol</a> DTIVOL M-file for DTIvol.fig</li></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 tic0002 0003 0004 0005 getDicomsgui <span class="comment">%uses my script to load an entire Dicom directory.</span>0006 0007 disp(<span class="string">'Calculating tensor...'</span>)0008 0009 noise=evalin(<span class="string">'base'</span>,<span class="string">'noise'</span>); <span class="comment">%noise threshold</span>0010 doFAmap=evalin(<span class="string">'base'</span>,<span class="string">'doFAmap'</span>); <span class="comment">%write DICOMs?</span>0011 dotADC=evalin(<span class="string">'base'</span>,<span class="string">'dotADC'</span>);0012 dolambda1=evalin(<span class="string">'base'</span>,<span class="string">'dolambda1'</span>);0013 dolambda2=evalin(<span class="string">'base'</span>,<span class="string">'dolambda2'</span>);0014 dolambda3=evalin(<span class="string">'base'</span>,<span class="string">'dolambda3'</span>);0015 docm=evalin(<span class="string">'base'</span>,<span class="string">'docm'</span>);0016 counter=evalin(<span class="string">'base'</span>,<span class="string">'counter'</span>); <span class="comment">%counter comes from getDicomsgui</span>0017 extension=evalin(<span class="string">'base'</span>,<span class="string">'extension'</span>);0018 0019 0020 slices=counter/7;0021 0022 0023 FAmaps=[];0024 tADCmaps=[];0025 lambda1_maps=[];0026 lambda2_maps=[];0027 lambda3_maps=[];0028 cmaps=[];0029 xyzmaps=[];0030 0031 0032 0033 0034 0035 0036 0037 warning off MATLAB:divideByZero <span class="comment">%supress matlab warning</span>0038 0039 0040 <span class="keyword">for</span> numslice=1:slices0041 0042 bo=tensorvol(:,:,numslice,1);<span class="comment">%bo=bo'; %transposes the images to correct for getDicoms</span>0043 b1=tensorvol(:,:,numslice,2);<span class="comment">%b1=b1';</span>0044 b2=tensorvol(:,:,numslice,3);<span class="comment">%b2=b2';</span>0045 b3=tensorvol(:,:,numslice,4);<span class="comment">%b3=b3';</span>0046 b4=tensorvol(:,:,numslice,5);<span class="comment">%b4=b4';</span>0047 b5=tensorvol(:,:,numslice,6);<span class="comment">%b5=b5';</span>0048 b6=tensorvol(:,:,numslice,7);<span class="comment">%b6=b6';</span>0049 0050 0051 ADC101=log(im2double(b1)./im2double(bo))./(-1000); <span class="comment">%Calculate ADC maps for each direction</span>0052 ADC_101=log(im2double(b2)./im2double(bo))./(-1000);0053 ADC011=log(im2double(b3)./im2double(bo))./(-1000);0054 ADC01_1=log(im2double(b4)./im2double(bo))./(-1000);0055 ADC110=log(im2double(b5)./im2double(bo))./(-1000);0056 ADC_110=log(im2double(b6)./im2double(bo))./(-1000);0057 0058 transform=[1 1 0 2 0 0;1 0 1 0 2 0;0 1 1 0 0 2;1 1 0 -2 0 0;1 0 1 0 -2 0;0 1 1 0 0 -2]; <span class="comment">%the transformation matrix to obtain tensor elements</span>0059 0060 xx=zeros(256); <span class="comment">%Create empty tensor elements</span>0061 yy=zeros(256);0062 zz=zeros(256);0063 xy=zeros(256);0064 xz=zeros(256);0065 yz=zeros(256);0066 lambda1=zeros(256); <span class="comment">%create empty images</span>0067 lambda2=zeros(256);0068 lambda3=zeros(256);0069 r=zeros(256);0070 g=zeros(256);0071 b=zeros(256);0072 FAmap=zeros(256);0073 tADC=zeros(256);0074 0075 <span class="keyword">for</span> i=1:256 <span class="comment">%Fill in tensor elements pixel by pixel.</span>0076 <span class="keyword">for</span> j=1:2560077 <span class="keyword">if</span> (bo(i,j)>noise) <span class="comment">%noise threshold. Change this to zero if you want the whole image to be calculated (takes longer)</span>0078 ADCm=[ADC110(i,j);ADC101(i,j);ADC011(i,j);ADC_110(i,j);ADC_101(i,j);ADC01_1(i,j)]; <span class="comment">%see Hedehus, http://www-radiology.stanford.edu/majh/</span>0079 ADCe=inv(transform)*ADCm;0080 xx(i,j)=ADCe(1,1); <span class="comment">%get tensor elements</span>0081 yy(i,j)=ADCe(2,1);0082 zz(i,j)=ADCe(3,1);0083 xy(i,j)=ADCe(4,1);0084 xz(i,j)=ADCe(5,1);0085 yz(i,j)=ADCe(6,1);0086 ten=[xx(i,j) xy(i,j) xz(i,j);xy(i,j) yy(i,j) yz(i,j);xz(i,j) yz(i,j) zz(i,j)]; <span class="comment">%the tensor itself</span>0087 [V,D]=eig(ten);0088 D=eig(ten); <span class="comment">%get eigenvalues</span>0089 D=abs(D); 0090 E=D; <span class="comment">%get the eigenvalues before re-ordering</span>0091 D=sort(D); <span class="comment">%sort the eigenvalues (upwards)</span>0092 e1=D(3,1);e1=e1*2;⌨️ 快捷键说明
复制代码
Ctrl + C
搜索代码
Ctrl + F
全屏模式
F11
切换主题
Ctrl + Shift + D
显示快捷键
?
增大字号
Ctrl + =
减小字号
Ctrl + -