Dear SPM experts,
I want to overlay my segmented images to the original unsegmented image to
get an idea about how well the segmentation worked. I just recently started
with SPM, so the only way I found to do this was running this file (and
slightly different versions if i just wanted to overlay one segmented image):
function ov = overlay_grey(P1,P2,P3,P4)
P1 = spm_select(1,'image'); % normalised unsegmented image
P2 = spm_select(1,'image'); % normalised modulated grey matter
P3 = spm_select(1,'image'); % normalised modulated white matter
P4 = spm_select(1,'image'); % normalised modulated CSF
% Clear graphics window..
spm_clf
% Display background image..
h = spm_orthviews('Image', P1,[0.05 0.05 0.9 0.9]);
% Display grey matter in red. Use [0 1 0] for green, [0 0 1] for blue
% [0.6 0 0.8] for purple etc..
spm_orthviews('AddColouredImage', h, P2,[1 0 0]);
% Display white matter in green.
spm_orthviews('AddColouredImage', h, P3,[0 1 0]);
% Display CSF in blue.
spm_orthviews('AddColouredImage', h, P4,[0 0 1]);
% Update the display..
spm_orthviews('Redraw');
return;
This doesn´t allow me to switch between the overlay(s) and the underlay. Can
anybody tell me how to solve this problem and is it also possible to display
the voxel coordinates?
Thank You in advance!
Beate
|