function mri_display_phasemap(anat, phasemap, ampmap, anat_thr, transp, radial, correction);

if nargin<5, transp = 1; end;
if nargin<6, radial = 0; end;
if nargin<7, correction = 0; end;

anat = conv3dToMosaic(anat, transp);
ampmap = conv3dToMosaic(ampmap, transp);
phasemap = conv3dToMosaic(phasemap, transp);

if size(anat,1)==2*size(phasemap,1),
    ampmap = ampmap(floor(1:0.5:size(phasemap,1)+0.5),:);
    phasemap = phasemap(floor(1:0.5:size(phasemap,1)+0.5),:);
end;
if size(anat,2)==2*size(phasemap,2),
    ampmap = ampmap(:,floor(1:0.5:size(phasemap,2)+0.5));
    phasemap = phasemap(:,floor(1:0.5:size(phasemap,2)+0.5));
end;

f = find(anat<anat_thr);
ampmap(f)=0;
% anat = anat-min(anat(:));
f = find(anat<0);
ima = 128*anat./max(1+anat(:));
ima(f) = 127;

f = find(ampmap>0);
if ~radial,
    ima(f) = 129+(phasemap(f)+pi)/pi*64;
    cmap = [gray(128); hsv(128)];
else
    cmap = [gray(128); jet(128)];
    ima(f) = 129+mod(phasemap(f)-correction, 2*pi)/pi*64;
end;

clf;
SUBPLOT('position',[0.11 0.11 0.8 0.815]);
f = find(ima<1); ima(f)=1;
image(ima); colormap(cmap); axis off;
if nargin>=2,
	SUBPLOT('position',[0.95 0.5 0.02 0.4]);
	image([1],[-64:64],[129:256]'); 
	set(gca, 'XTickLabel',[]);
	set(gca, 'XTick',[]);
    if ~radial,
		set(gca, 'YTick', [-64 0 64]);
		set(gca, 'YTickLabel', [90 0 -90]);
        title('a');
    else;
		set(gca, 'YTick', [-64 0 64]);
		set(gca, 'YTickLabel', [1 0.5 0]);
        title('e');
    end;
end;