dtiguicode.html

mri_toolbox是一个工具用来MRI. 来自于SourceForge, 我上传这个软件,希望能结识对医疗软件感兴趣的兄弟.

0093             e2=D(2,1);e2=e2*2;
0094             e3=D(1,1);e3=e3*2;
0095             lambda1(i,j)=e1;
0096             lambda2(i,j)=e2;
0097             lambda3(i,j)=e3;
0098             <span class="comment">%V=abs(V);       %get the x,y,z components of e3</span>
0099             largest=find(E==D(3,1));                            <span class="comment">%find the value of the largest eigenvalue in the original matrix</span>
0100             sz=length(largest);
0101             r(i,j)=V(1,largest(sz));                                <span class="comment">%in order to get its corresponding eigenvector components in x, y and z</span>
0102             g(i,j)=V(2,largest(sz));
0103             b(i,j)=V(3,largest(sz));
0104             trace=(e1+e2+e3)/3;
0105             <span class="comment">%FA=(sqrt(3/2)*sqrt((1/3)*((e1-e2)^2+(e2-e3)^2+(e3-e1)^2)))/(sqrt(e1^2+e2^2+e3^3)); %Another formula for FA</span>
0106             FA=(sqrt(3*((e1-trace)^2+(e2-trace)^2+(e3-trace)^2)))/(sqrt(2*(e1^2+e2^2+e3^2)));  <span class="comment">%Le Bihan 2001</span>
0107             tADC(i,j)=trace;
0108             FAmap(i,j)=FA;
0109         <span class="keyword">end</span>;
0110     <span class="keyword">end</span>;
0111 <span class="keyword">end</span>;
0112 
0113 
0114 
0115 
0116 rgb=cat(3,r,g,b);
0117 xyz=cat(3,r,(g*-1),b);
0118 FAmap3=cat(3,FAmap,FAmap,FAmap);                                <span class="comment">%needed to combine with the colormap, need 3D.</span>
0119 cm=abs(rgb).*FAmap3;                                            <span class="comment">%FA weighting for colormap.  Has absolute values for vector components</span>
0120 
0121 <span class="comment">%eval (['rgb' num2str(numslice) '=rgb;']);</span>
0122 <span class="comment">%nameofmap=['rgb' num2str(numslice)];</span>
0123 <span class="comment">%assignin('base',nameofmap,rgb);</span>
0124 
0125 
0126 FAmaps=cat(3,FAmaps,FAmap);
0127 tADCmaps=cat(3,tADCmaps,tADC);
0128 lambda1_maps=cat(3,lambda1_maps,lambda1);
0129 lambda2_maps=cat(3,lambda2_maps,lambda2);
0130 lambda3_maps=cat(3,lambda3_maps,lambda3);
0131 cmaps=cat(4,cmaps,cm);
0132 xyzmaps=cat(4,xyzmaps,xyz);
0133 
0134 
0135 
0136 
0137 
0138 
0139 <span class="comment">%eval (['xyz' num2str(numslice) '=rgb;']);                       %xyz is the same as cm (colormap) but without abs.  But, apparently, doesn't make a difference.</span>
0140 <span class="comment">%nameofmap=['xyz' num2str(numslice)];</span>
0141 <span class="comment">%assignin('base',nameofmap,rgb);</span>
0142 
0143 <span class="comment">%image (cm);axis image</span>
0144 <span class="comment">%eval (['FA' num2str(numslice) '=FAmap;']);</span>
0145 <span class="comment">%nameofmap=['FA' num2str(numslice)];</span>
0146 <span class="comment">%assignin('base',nameofmap,FAmap);</span>
0147 <span class="keyword">if</span> (doFAmap==1)
0148     number=num2str(numslice);                                   <span class="comment">%FOUR LINES FOR FILE WRITING, REPEAT FOR THE OTHER IMAGES</span>
0149     st1=<span class="string">'FA'</span>;st2=<span class="string">'.dcm'</span>;
0150     nameforfile=[st1 number st2];
0151     dicomwrite(FAmap,nameforfile);
0152     disp(<span class="string">'writing FA map to DICOM file'</span>)
0153 <span class="keyword">end</span>
0154 <span class="comment">%eval (['tADC' num2str(numslice) '=tADC;']);</span>
0155 <span class="comment">%nameofmap=['tADC' num2str(numslice)];</span>
0156 <span class="comment">%assignin('base',nameofmap,tADC);</span>
0157 <span class="keyword">if</span> (dotADC==1)
0158     number=num2str(numslice);                                   <span class="comment">%FOUR LINES FOR FILE WRITING, REPEAT FOR THE OTHER IMAGES</span>
0159     st1=<span class="string">'tADC'</span>;st2=<span class="string">'.dcm'</span>;
0160     nameforfile=[st1 number st2];
0161     dicomwrite(tADC,nameforfile);
0162     disp(<span class="string">'writing tADC map to DICOM file'</span>)
0163 <span class="keyword">end</span>
0164 <span class="comment">%eval (['lambda1' num2str(numslice) '=lambda1;']);</span>
0165 <span class="comment">%nameofmap=['lambda1_' num2str(numslice)];</span>
0166 <span class="comment">%assignin('base',nameofmap,lambda1);</span>
0167 <span class="keyword">if</span> (dolambda1==1)
0168     number=num2str(numslice);                                   <span class="comment">%FOUR LINES FOR FILE WRITING, REPEAT FOR THE OTHER IMAGES</span>
0169     st1=<span class="string">'lambda1_'</span>;st2=<span class="string">'.dcm'</span>;
0170     nameforfile=[st1 number st2];
0171     dicomwrite(lambda1,nameforfile);
0172     disp(<span class="string">'writing lambda1 map to DICOM file'</span>)
0173 <span class="keyword">end</span>
0174 <span class="comment">%eval (['lambda2' num2str(numslice) '=lambda2;']);</span>
0175 <span class="comment">%nameofmap=['lambda2_' num2str(numslice)];</span>
0176 <span class="comment">%assignin('base',nameofmap,lambda2);</span>
0177 <span class="keyword">if</span> (dolambda2==1)
0178     number=num2str(numslice);                                   <span class="comment">%FOUR LINES FOR FILE WRITING, REPEAT FOR THE OTHER IMAGES</span>
0179     st1=<span class="string">'lambda2_'</span>;st2=<span class="string">'.dcm'</span>;
0180     nameforfile=[st1 number st2];
0181     dicomwrite(lambda2,nameforfile); 
0182     disp(<span class="string">'writing lambda2 map to DICOM file'</span>)
0183 <span class="keyword">end</span>
0184 <span class="comment">%eval (['lambda3' num2str(numslice) '=lambda3;']);</span>
0185 <span class="comment">%nameofmap=['lambda3_' num2str(numslice)];</span>
0186 <span class="comment">%assignin('base',nameofmap,lambda3);</span>
0187 <span class="keyword">if</span> (dolambda3==1)
0188     number=num2str(numslice);                                   <span class="comment">%FOUR LINES FOR FILE WRITING, REPEAT FOR THE OTHER IMAGES</span>
0189     st1=<span class="string">'lambda3_'</span>;st2=<span class="string">'.dcm'</span>;
0190     nameforfile=[st1 number st2];
0191     dicomwrite(lambda3,nameforfile);
0192     disp(<span class="string">'writing lambda3 map to DICOM file'</span>)
0193 <span class="keyword">end</span>
0194 <span class="comment">%eval (['cm' num2str(numslice) '=cm;']);</span>
0195 <span class="comment">%nameofmap=['cm' num2str(numslice)];</span>
0196 <span class="comment">%assignin('base',nameofmap,cm);</span>
0197 <span class="keyword">if</span> (docm==1)
0198     number=num2str(numslice);                                   <span class="comment">%FOUR LINES FOR FILE WRITING, REPEAT FOR THE OTHER IMAGES</span>
0199     st1=<span class="string">'cm'</span>;st2=<span class="string">'.dcm'</span>;
0200     nameforfile=[st1 number st2];
0201     dicomwrite(cm,nameforfile); 
0202     disp(<span class="string">'writing colormap to DICOM file'</span>)
0203 <span class="keyword">end</span>
0204 
0205 mes=[num2str(numslice) <span class="string">' of '</span> num2str(slices) <span class="string">' finished'</span>];
0206 disp(mes)
0207 <span class="keyword">end</span>
0208  
0209 <span class="comment">%assignin('base','ten',ten)</span>
0210 clear FA FAmap tADC lambda1 lambda2 lambda3 cm FAmap3 rgb bo b1 b2 b3 b4 b5 b6 ADC101 ADC_101 ADC011;
0211 clear ADC01_1 ADC110 ADC_110 r g b xx yy zz xy xz yz trace ten tensorvol transform ADCe ADCm D V;
0212 assignin(<span class="string">'base'</span>,<span class="string">'tensorvol'</span>,0)
0213 
0214 
0215 assignin(<span class="string">'base'</span>,<span class="string">'FAmaps'</span>,FAmaps);
0216 assignin(<span class="string">'base'</span>,<span class="string">'tADCmaps'</span>,tADCmaps);
0217 assignin(<span class="string">'base'</span>,<span class="string">'lambda1_maps'</span>,lambda1_maps);
0218 assignin(<span class="string">'base'</span>,<span class="string">'lambda2_maps'</span>,lambda2_maps);
0219 assignin(<span class="string">'base'</span>,<span class="string">'lambda3_maps'</span>,lambda3_maps);
0220 assignin(<span class="string">'base'</span>,<span class="string">'cmaps'</span>,cmaps);
0221 assignin(<span class="string">'base'</span>,<span class="string">'xyzmaps'</span>,xyzmaps);
0222 
0223 toc</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> &copy; 2003</address>
</body>
</html>