avw_view_old.html

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

0120         <span class="string">'   set(AVWVIEW.handles.Ssag,''Value'',SagSlice);'</span>,<span class="keyword">...</span>
0121         <span class="string">'   clear Ssag; '</span>,<span class="keyword">...</span>
0122         <span class="string">'   set(AVWVIEW.gui,''UserData'',AVWVIEW);'</span>,<span class="keyword">...</span>
0123         <span class="string">'end; '</span>,<span class="keyword">...</span>
0124         <span class="string">'if ishandle(AVWVIEW.handles.Hc) &amp; AVWVIEW.handles.Hc, '</span>,<span class="keyword">...</span>
0125         <span class="string">'   Scor = squeeze(AVWVIEW.avw.img(:,CorSlice,:));'</span>,<span class="keyword">...</span>
0126         <span class="string">'   set(AVWVIEW.handles.Hc,''CData'',Scor); drawnow;'</span>,<span class="keyword">...</span>
0127         <span class="string">'   set(AVWVIEW.handles.Tcor,''String'',num2str(CorSlice));'</span>,<span class="keyword">...</span>
0128         <span class="string">'   set(AVWVIEW.handles.Scor,''Value'',CorSlice);'</span>,<span class="keyword">...</span>
0129         <span class="string">'   clear Scor; '</span>,<span class="keyword">...</span>
0130         <span class="string">'   set(AVWVIEW.gui,''UserData'',AVWVIEW);'</span>,<span class="keyword">...</span>
0131         <span class="string">'end; '</span>,<span class="keyword">...</span>
0132         <span class="string">'clear currentpoint imgvalue AxiSlice CorSlice SagSlice AVWVIEW;'</span>));
0133     
0134     <span class="keyword">if</span> zdim &gt; 1,
0135         slider_step(1) = 1/(zdim);
0136         slider_step(2) = 1/(zdim);
0137         G.Saxi = uicontrol(<span class="string">'Parent'</span>,GUI,<span class="string">'Style'</span>,<span class="string">'slider'</span>,<span class="string">'Units'</span>,<span class="string">'Normalized'</span>, Font, <span class="keyword">...</span>
0138             <span class="string">'Position'</span>,[.45 .56 .03 .40], <span class="string">'HorizontalAlignment'</span>, <span class="string">'center'</span>,<span class="keyword">...</span>
0139             <span class="string">'BusyAction'</span>,<span class="string">'queue'</span>,<span class="keyword">...</span>
0140             <span class="string">'Min'</span>,1,<span class="string">'Max'</span>,zdim,<span class="string">'SliderStep'</span>,slider_step,<span class="string">'Value'</span>,AxiSlice,<span class="keyword">...</span>
0141             <span class="string">'Callback'</span>,strcat(<span class="string">'AVWVIEW = get(gcbf,''Userdata''); '</span>,<span class="keyword">...</span>
0142             <span class="string">'AxiSlice = round(get(AVWVIEW.handles.Saxi,''Value''));'</span>,<span class="keyword">...</span>
0143             <span class="string">'set(AVWVIEW.handles.Saxi,''Value'',AxiSlice);'</span>,<span class="keyword">...</span>
0144             <span class="string">'Saxi = squeeze(AVWVIEW.avw.img(:,:,AxiSlice));'</span>,<span class="keyword">...</span>
0145             <span class="string">'set(AVWVIEW.handles.Ha,''CData'',Saxi); drawnow;'</span>,<span class="keyword">...</span>
0146             <span class="string">'set(AVWVIEW.handles.Taxi,''String'',num2str(AxiSlice));'</span>,<span class="keyword">...</span>
0147             <span class="string">'CorSlice = round(get(AVWVIEW.handles.Scor,''Value''));'</span>,<span class="keyword">...</span>
0148             <span class="string">'SagSlice = round(get(AVWVIEW.handles.Ssag,''Value''));'</span>,<span class="keyword">...</span>
0149             <span class="string">'imgvalue = AVWVIEW.avw.img(SagSlice,CorSlice,AxiSlice); '</span>,<span class="keyword">...</span>
0150             <span class="string">'set(AVWVIEW.handles.imval,''String'',sprintf(''%8.2f'',imgvalue));'</span>,<span class="keyword">...</span>
0151             <span class="string">'set(AVWVIEW.gui,''UserData'',AVWVIEW);'</span>,<span class="keyword">...</span>
0152             <span class="string">'clear imgvalue Saxi AxiSlice CorSlice SagSlice AVWVIEW;'</span>));
0153     <span class="keyword">end</span>
0154     G.Taxi = uicontrol(<span class="string">'Parent'</span>,GUI,<span class="string">'Style'</span>,<span class="string">'text'</span>,<span class="string">'Units'</span>,<span class="string">'Normalized'</span>, Font, <span class="keyword">...</span>
0155         <span class="string">'Position'</span>,[.45 .51 .03 .05], <span class="string">'HorizontalAlignment'</span>, <span class="string">'center'</span>,<span class="keyword">...</span>
0156         <span class="string">'BusyAction'</span>,<span class="string">'queue'</span>,<span class="keyword">...</span>
0157         <span class="string">'String'</span>,num2str(AxiSlice));
0158 <span class="keyword">end</span>
0159 
0160 <span class="comment">%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%</span>
0161 <span class="comment">%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%</span>
0162 
0163 <span class="comment">% Coronal Slice</span>
0164 <span class="keyword">if</span> xdim &gt; 1 &amp; zdim &gt; 1,
0165     
0166     [x,z] = meshgrid(1:xdim,1:zdim);
0167     Xcor = x'; clear x;
0168     Zcor = z'; clear z;
0169     Ycor = zeros(xdim,zdim);
0170     
0171     subplot(<span class="string">'position'</span>,[0.55 0.56 0.4 0.4])
0172     colormap(<span class="string">'gray'</span>);
0173     Scor = squeeze(avw.img(:,CorSlice,:));
0174     
0175     G.Hc = surf(Xcor,Ycor,Zcor,Scor,<span class="string">'EdgeColor'</span>,<span class="string">'none'</span>);
0176     axis square, daspect([1,1,1]);
0177     xlabel(<span class="string">'(Left &lt;&lt;) X (&gt;&gt; Right)'</span>)
0178     ylabel(<span class="string">'Y'</span>)
0179     zlabel(<span class="string">'Z'</span>)
0180     title(<span class="string">'Coronal'</span>)
0181     view([0,0]);
0182     
0183     <span class="comment">% This callback navigates with left click</span>
0184     set(G.Hc,<span class="string">'ButtonDownFcn'</span>,<span class="keyword">...</span>
0185         strcat(<span class="string">'AVWVIEW = get(gcbf,''Userdata''); '</span>,<span class="keyword">...</span>
0186         <span class="string">'currentpoint = get(get(AVWVIEW.handles.Hc,''Parent''),''CurrentPoint''); '</span>,<span class="keyword">...</span>
0187         <span class="string">'SagSlice = round(currentpoint(2,1)); '</span>,<span class="keyword">...</span>
0188         <span class="string">'AxiSlice = round(currentpoint(2,3)); '</span>,<span class="keyword">...</span>
0189         <span class="string">'CorSlice = round(str2num(get(AVWVIEW.handles.Tcor,''String''))); '</span>,<span class="keyword">...</span>
0190         <span class="string">'imgvalue = AVWVIEW.avw.img(SagSlice,CorSlice,AxiSlice); '</span>,<span class="keyword">...</span>
0191         <span class="string">'set(AVWVIEW.handles.imval,''String'',sprintf(''%8.2f'',imgvalue));'</span>,<span class="keyword">...</span>
0192         <span class="string">'set(AVWVIEW.gui,''UserData'',AVWVIEW);'</span>,<span class="keyword">...</span>
0193         <span class="string">'if ishandle(AVWVIEW.handles.Hs) &amp; AVWVIEW.handles.Hs, '</span>,<span class="keyword">...</span>
0194         <span class="string">'   Ssag = squeeze(AVWVIEW.avw.img(SagSlice,:,:));'</span>,<span class="keyword">...</span>
0195         <span class="string">'   set(AVWVIEW.handles.Hs,''CData'',Ssag); drawnow;'</span>,<span class="keyword">...</span>
0196         <span class="string">'   set(AVWVIEW.handles.Tsag,''String'',num2str(SagSlice));'</span>,<span class="keyword">...</span>
0197         <span class="string">'   set(AVWVIEW.handles.Ssag,''Value'',SagSlice);'</span>,<span class="keyword">...</span>
0198         <span class="string">'   clear Ssag; '</span>,<span class="keyword">...</span>
0199         <span class="string">'   set(AVWVIEW.gui,''UserData'',AVWVIEW);'</span>,<span class="keyword">...</span>
0200         <span class="string">'end; '</span>,<span class="keyword">...</span>
0201         <span class="string">'if ishandle(AVWVIEW.handles.Ha) &amp; AVWVIEW.handles.Ha, '</span>,<span class="keyword">...</span>
0202         <span class="string">'   Saxi = squeeze(AVWVIEW.avw.img(:,:,AxiSlice));'</span>,<span class="keyword">...</span>
0203         <span class="string">'   set(AVWVIEW.handles.Ha,''CData'',Saxi); drawnow;'</span>,<span class="keyword">...</span>
0204         <span class="string">'   set(AVWVIEW.handles.Taxi,''String'',num2str(AxiSlice));'</span>,<span class="keyword">...</span>
0205         <span class="string">'   set(AVWVIEW.handles.Saxi,''Value'',AxiSlice);'</span>,<span class="keyword">...</span>
0206         <span class="string">'   clear Saxi; '</span>,<span class="keyword">...</span>
0207         <span class="string">'   set(AVWVIEW.gui,''UserData'',AVWVIEW);'</span>,<span class="keyword">...</span>
0208         <span class="string">'end; '</span>,<span class="keyword">...</span>
0209         <span class="string">'clear currentpoint imgvalue AxiSlice CorSlice SagSlice AVWVIEW;'</span>));
0210     
0211     <span class="keyword">if</span> ydim &gt; 1,
0212         slider_step(1) = 1/(ydim);
0213         slider_step(2) = 1/(ydim);
0214         G.Scor = uicontrol(<span class="string">'Parent'</span>,GUI,<span class="string">'Style'</span>,<span class="string">'slider'</span>,<span class="string">'Units'</span>,<span class="string">'Normalized'</span>, Font, <span class="keyword">...</span>
0215             <span class="string">'Position'</span>,[.95 .56 .03 .40], <span class="string">'HorizontalAlignment'</span>, <span class="string">'center'</span>,<span class="keyword">...</span>
0216             <span class="string">'BusyAction'</span>,<span class="string">'queue'</span>,<span class="keyword">...</span>
0217             <span class="string">'Min'</span>,1,<span class="string">'Max'</span>,ydim,<span class="string">'SliderStep'</span>,slider_step,<span class="string">'Value'</span>,CorSlice,<span class="keyword">...</span>
0218             <span class="string">'Callback'</span>,strcat(<span class="string">'AVWVIEW = get(gcbf,''Userdata''); '</span>,<span class="keyword">...</span>
0219             <span class="string">'CorSlice = round(get(AVWVIEW.handles.Scor,''Value''));'</span>,<span class="keyword">...</span>
0220             <span class="string">'set(AVWVIEW.handles.Scor,''Value'',CorSlice);'</span>,<span class="keyword">...</span>
0221             <span class="string">'Scor = squeeze(AVWVIEW.avw.img(:,CorSlice,:));'</span>,<span class="keyword">...</span>
0222             <span class="string">'set(AVWVIEW.handles.Hc,''CData'',Scor); drawnow;'</span>,<span class="keyword">...</span>
0223             <span class="string">'set(AVWVIEW.handles.Tcor,''String'',num2str(CorSlice));'</span>,<span class="keyword">...</span>
0224             <span class="string">'AxiSlice = round(get(AVWVIEW.handles.Saxi,''Value''));'</span>,<span class="keyword">...</span>
0225             <span class="string">'SagSlice = round(get(AVWVIEW.handles.Ssag,''Value''));'</span>,<span class="keyword">...</span>
0226             <span class="string">'imgvalue = AVWVIEW.avw.img(SagSlice,CorSlice,AxiSlice); '</span>,<span class="keyword">...</span>
0227             <span class="string">'set(AVWVIEW.handles.imval,''String'',sprintf(''%8.2f'',imgvalue));'</span>,<span class="keyword">...</span>
0228             <span class="string">'set(AVWVIEW.gui,''UserData'',AVWVIEW);'</span>,<span class="keyword">...</span>
0229             <span class="string">'clear imgvalue Scor AxiSlice CorSlice SagSlice AVWVIEW;'</span>));
0230     <span class="keyword">end</span>
0231     G.Tcor = uicontrol(<span class="string">'Parent'</span>,GUI,<span class="string">'Style'</span>,<span class="string">'text'</span>,<span class="string">'Units'</span>,<span class="string">'Normalized'</span>, Font, <span class="keyword">...</span>
0232         <span class="string">'Position'</span>,[.95 .51 .03 .05], <span class="string">'HorizontalAlignment'</span>, <span class="string">'center'</span>,<span class="keyword">...</span>
0233         <span class="string">'BusyAction'</span>,<span class="string">'queue'</span>,<span class="keyword">...</span>
0234         <span class="string">'String'</span>,num2str(CorSlice));
0235 <span class="keyword">end</span>
0236 
0237 <span class="comment">%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%</span>
0238 <span class="comment">%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%</span>
0239 
0240 <span class="comment">% Sagittal Slice</span>
0241 <span class="keyword">if</span> ydim &gt; 1 &amp; ydim &gt; 1,
0242         
0243     [y,z] = meshgrid(1:ydim,1:zdim);
0244     Ysag = y'; clear y;
0245     Zsag = z'; clear z;
0246     Xsag = zeros(ydim,zdim);
0247     
0248     subplot(<span class="string">'position'</span>,[0.05 0.06 0.4 0.4])
0249     colormap(<span class="string">'gray'</span>);
0250     Ssag = squeeze(avw.img(SagSlice,:,:));
0251     
0252     G.Hs = surf(Xsag,Ysag,Zsag,Ssag,<span class="string">'EdgeColor'</span>,<span class="string">'none'</span>);
0253     axis square, daspect([1,1,1]);
0254     xlabel(<span class="string">'(Left &lt;&lt;) X (&gt;&gt; Right)'</span>)
0255     ylabel(<span class="string">'Y'</span>)
0256     zlabel(<span class="string">'Z'</span>)
0257     title(<span class="string">'Sagittal'</span>)
0258     view([90,0]);
0259     
0260     <span class="comment">% This callback navigates with mouse click</span>
0261     set(G.Hs,<span class="string">'ButtonDownFcn'</span>,<span class="keyword">...</span>
0262         strcat(<span class="string">'AVWVIEW = get(gcbf,''Userdata''); '</span>,<span class="keyword">...</span>
0263         <span class="string">'currentpoint = get(get(AVWVIEW.handles.Hs,''Parent''),''CurrentPoint''); '</span>,<span class="keyword">...</span>
0264         <span class="string">'CorSlice = round(currentpoint(1,2)); '</span>,<span class="keyword">...</span>
0265         <span class="string">'AxiSlice = round(currentpoint(1,3)); '</span>,<span class="keyword">...</span>
0266         <span class="string">'SagSlice = round(str2num(get(AVWVIEW.handles.Tsag,''String'')));'</span>,<span class="keyword">...</span>
0267         <span class="string">'imgvalue = AVWVIEW.avw.img(SagSlice,CorSlice,AxiSlice); '</span>,<span class="keyword">...</span>
0268         <span class="string">'set(AVWVIEW.handles.imval,''String'',sprintf(''%8.2f'',imgvalue));'</span>,<span class="keyword">...</span>
0269         <span class="string">'set(AVWVIEW.gui,''UserData'',AVWVIEW);'</span>,<span class="keyword">...</span>
0270         <span class="string">'if ishandle(AVWVIEW.handles.Hc) &amp; AVWVIEW.handles.Hc, '</span>,<span class="keyword">...</span>
0271         <span class="string">'   Scor = squeeze(AVWVIEW.avw.img(:,CorSlice,:));'</span>,<span class="keyword">...</span>
0272         <span class="string">'   set(AVWVIEW.handles.Hc,''CData'',Scor); drawnow;'</span>,<span class="keyword">...</span>
0273         <span class="string">'   set(AVWVIEW.handles.Tcor,''String'',num2str(CorSlice));'</span>,<span class="keyword">...</span>
0274         <span class="string">'   set(AVWVIEW.handles.Scor,''Value'',CorSlice);'</span>,<span class="keyword">...</span>
0275         <span class="string">'   clear Scor; '</span>,<span class="keyword">...</span>
0276         <span class="string">'   set(AVWVIEW.gui,''UserData'',AVWVIEW);'</span>,<span class="keyword">...</span>
0277         <span class="string">'end; '</span>,<span class="keyword">...</span>
0278         <span class="string">'if ishandle(AVWVIEW.handles.Ha) &amp; AVWVIEW.handles.Ha, '</span>,<span class="keyword">...</span>
0279         <span class="string">'   Saxi = squeeze(AVWVIEW.avw.img(:,:,AxiSlice));'</span>,<span class="keyword">...</span>
0280         <span class="string">'   set(AVWVIEW.handles.Ha,''CData'',Saxi); drawnow;'</span>,<span class="keyword">...</span>
0281         <span class="string">'   set(AVWVIEW.handles.Taxi,''String'',num2str(AxiSlice));'</span>,<span class="keyword">...</span>
0282         <span class="string">'   set(AVWVIEW.handles.Saxi,''Value'',AxiSlice);'</span>,<span class="keyword">...</span>
0283         <span class="string">'   clear Saxi; '</span>,<span class="keyword">...</span>
0284         <span class="string">'   set(AVWVIEW.gui,''UserData'',AVWVIEW);'</span>,<span class="keyword">...</span>
0285         <span class="string">'end; '</span>,<span class="keyword">...</span>
0286         <span class="string">'clear currentpoint imgvalue AxiSlice CorSlice SagSlice AVWVIEW;'</span>));
0287 
0288     
0289     <span class="keyword">if</span> xdim &gt; 1,
0290         slider_step(1) = 1/(xdim);