Hi Randall, List,
Try the isosurface function in matlab, it is great to extract polygon
renderings from volumatric data (e.g. T-maps, gray matter volumes) etc.
Type 'help isosurface' or see:
http://www.mathworks.com/access/helpdesk/help/techdoc/ref/isosurface.html
Together with reducepatch (reduces too large number of polygon meshes)
this gives realy cool results (see attachment for an example).
I.e. something like gives a skin rendering (set alpha blending to a
value between 0-1 for opaque surfaces):
img='/home/me/mydata/fmri/someT1volume.img';
V=spm_vol(img);
[Data,XYZ]=spm_read_vols(V);
X=zeros(size(Data));
Y=X;
Z=X;
X(:)=XYZ(1,:);
Y(:)=XYZ(2,:);
Z(:)=XYZ(3,:);
FV=isosurface(X,Y,Z,Data,100,'verbose');
NFV=reducepatch(FV,0.1);
myp = patch(NFV);
set(myp, 'FaceColor', [1 0.9 0.8], 'EdgeColor', 'none');
daspect([1 1 1])
view(3)
camlight; lighting phong
Then do something similar with gray matter volume, T-maps, and you are
in business.
When you have matlab you don't need nothing else...
Cheers,
Bas
]Op ma, 15-08-2005 te 16:35 +0100, schreef Randall Benson:
> SPMers,
>
> I am trying in vain to use Medx and SPM to create surface renderings of 3
> different contrasts (phonetics, FM, object processing) displayed
> simultaneously using 3 different colors. Contrasts and analyses were
> created with SPM99 or SPM2 (either). Maybe there is an extension to SPM
> that does it but don't know what it is. It seems like such a common need I
> would not think I would have to look far for this. I know Analyze does it
> but I don't have the program any more. In Medx it's easy to display slices
> and orthogonals with multiple overlays but no easy way to turn this into a
> rendering that I can see.
>
> Appreciate any help.
>
> Randy Benson, M.D.
>
> Wayne State University & Detroit Medical Center
|