untitled.m

This code helps to fuse pet and ct scan images.The file untitled.m opens the interface to load and f

% --- Executes on button press in pushbutton4.
function pushbutton4_Callback(hObject, eventdata, handles)
% hObject    handle to pushbutton4 (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    structure with handles and user data (see GUIDATA)

 
% FUNCTION ZOOM2CURSOR
%
% zoom2cursor, without arguments, will activate the current axis, create a text box showing the
% current position of the mouse pointer (similar to pixval), and automatically zoom the image to the
% location of the cursor  as it is moved. The zoomed display dynamically scrolls with the motion of the cursor.
%
% By default, the function zooms to 50% of the image in the axis.
%
% BUTTON CLICKS:
% Left-clicking will zoom in further, and right-clicking will zoom out.
% Shift-clicking (or simultaneously clicking the left and right mouse buttons) at any point
% will display the original (un-zoomed) image, as will moving the
% cursor outside of the current axis. The zoom percentage is restored when the mouse is moved.
% Double-clicking zooms out to the original image, modifying the zoom percentage.
%
% Tested under R12.1 and R13.
%
% Written by Brett Shoelson, Ph.D. ([e-mail], [e-mail])
% 12/26/02
% 2/16/03; Rev 2: Program is more robust; fixes a bug when window is resized.
%                 Incremental increase/decrease in zoom percent (on mouseclick) has been reduced.
%                 Also: Now works with images, surfaces, lines (and thus plots), and patches (rather than just images)
global result;
j=result;
imshow(j);
currfig = findobj('type','figure');
if isempty(currfig)
    beep;
    h=warndlg('There are no figures open!');
    uiwait(h);
    return
end
currfig = currfig(1);
figure(currfig);
%zoomparams.currax = findobj(currfig,'type','axes');
zoomparams.currax = get(currfig,'currentaxes');
if isempty(zoomparams.currax)
    error('The current figure contains no axes!');
end

%Precedence: Images, surfaces, lines, patches
%Are there any images in the current axes?
zoomparams.currobj = findobj(currfig,'type','image');
if isempty(zoomparams.currobj)
    zoomparams.currobj = findobj(currfig,'type','surface');
end
if isempty(zoomparams.currobj)
    zoomparams.currobj = findobj(currfig,'type','line');
end
if isempty(zoomparams.currobj)
    zoomparams.currobj = findobj(currfig,'type','patch');
end
if isempty(zoomparams.currobj)
    error('The current axis doesn''t appear to contain any valid images, surfaces, lines, or patches.');
end

zoomparams.pct = 0.5; %Default value

zoomparams.currobj = zoomparams.currobj(1);
zoomparams.objtype = get(zoomparams.currobj,'type');
zoomparams.bdfcnold = get(zoomparams.currobj,'buttondownfcn');
zoomparams.baold = get(zoomparams.currobj,'busyaction');
zoomparams.oldpointer = get(currfig,'pointer');
set(currfig, 'Pointer', 'crosshair');
axes(zoomparams.currax);

currax = get(zoomparams.currax);
zoomparams.refax = copyobj(zoomparams.currax,gcf);

% For simplicity, I store the bdfunction in both the current axis AND the current object. For images, it makes sense to
% store it in the object (since the object covers the axis). For other objects (like lines), storing in the object forces
% a click directly on the line/point, but storing in the axis only means a click on the line does not trigger the callback.
warning off; %I turn this off because of the annoying (but erroneous)"Unrecognized OpenGL" message.
set(zoomparams.currax,'buttondownfcn','feval(getappdata(gcf,''bdfcnhandle''));','busyaction','queue');
set(zoomparams.currobj,'buttondownfcn','feval(getappdata(gcf,''bdfcnhandle''));','busyaction','queue');

set(findobj(zoomparams.refax,'type','children'),'handlevisibility','on');
set(zoomparams.refax,'visible','off');
axes(zoomparams.refax);
cla;

zoomparams.oldaxunits = get(zoomparams.currax,'units');
zoomparams.ydir = get(zoomparams.currax,'ydir');

zoomparams.oldxlim = get(zoomparams.currax,'xlim');
zoomparams.oldylim = get(zoomparams.currax,'ylim');
zoomparams.oldzlim = get(zoomparams.currax,'zlim');
zoomparams.dbold = get(currfig,'doublebuffer');
zoomparams.xrange = diff(zoomparams.oldxlim);
zoomparams.yrange = diff(zoomparams.oldylim);
zoomparams.zrange = diff(zoomparams.oldzlim);
zoomparams.xdist = zoomparams.pct*zoomparams.xrange;
zoomparams.ydist = zoomparams.pct*zoomparams.yrange;
zoomparams.zdist = zoomparams.pct*zoomparams.zrange;
zoomparams.oldwbmf = get(currfig,'windowbuttonmotionfcn');

setappdata(currfig, 'zoomfcnhandle', @zoomfcn);
setappdata(gcf, 'bdfcnhandle',@bdfcn);

set(currfig,'doublebuffer','on','windowbuttonmotionfcn','feval(getappdata(gcf,''zoomfcnhandle''));');
endbutton = uicontrol('style','pushbutton','string','X','units','normalized','position',[0 0 0.04 0.065],'foregroundcolor','r','backgroundcolor','k',...
    'fontweight','b','callback',...
    ['zoomparams = getappdata(gcf,''zoomparams'');set(gcf,''windowbuttonmotionfcn'',zoomparams.oldwbmf);',...
        'set(zoomparams.currax,''units'',zoomparams.oldaxunits,''xlim'',zoomparams.oldxlim,''ylim'',zoomparams.oldylim);',...
		'set(get(zoomparams.currax,''children''),''buttondownfcn'',zoomparams.bdfcnold,''busyaction'',zoomparams.baold);',...
        'set(gcf,''pointer'',zoomparams.oldpointer,''doublebuffer'',zoomparams.dbold);,',...
        'delete(zoomparams.dispbox1);delete(zoomparams.dispbox2);delete(gcbo);']);
zoomparams.dispbox1 = uicontrol('style','frame','backgroundcolor','k','units','normalized','position',[0.0475 0 0.35 0.065]);
if ~strcmp(zoomparams.objtype,'surface')
    msgstr = sprintf('x = %3.0f;  y = %3.0f',0,0);
else
    msgstr = sprintf('x = %3.0f;  y = %3.0f; z = %3.0f',0,0,0);
end
zoomparams.dispbox2 = uicontrol('style','text','backgroundcolor','k','foregroundcolor','y','units','normalized',...
	'position',[0.065 0.02 0.32 0.04],'string',msgstr,...
    'horizontalalignment','l');
setappdata(gcf,'zoomparams',zoomparams);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

function zoomfcn
zoomparams = getappdata(gcf,'zoomparams');
posn = get(zoomparams.refax,'currentpoint');
posn = posn(1,:);

x = posn(1,1);
y = posn(1,2);
z = posn(1,3);

switch zoomparams.objtype
    case 'image'
        % x and y are already in expressed in proper pixel coordinates
        x1 = min(max(1,x-0.5*zoomparams.xdist),zoomparams.xrange-zoomparams.xdist) + 0.5;
        y1 = min(max(1,y-0.5*zoomparams.ydist),zoomparams.yrange-zoomparams.ydist) + 0.5;
        z1 = min(max(1,z-0.5*zoomparams.zdist),zoomparams.zrange-zoomparams.zdist) + 0.5;
        x2 = x1 + zoomparams.xdist;
        y2 = y1 + zoomparams.ydist;
        z2 = z1 + zoomparams.zdist;
    case {'line','surface','patch'}
        % x, y and z are in normalized units; must be converted
        x = zoomparams.oldxlim(1) + x*zoomparams.xrange;
        y = zoomparams.oldylim(1) + y*zoomparams.yrange;
        z = zoomparams.oldzlim(1) + z*zoomparams.zrange;
        % now change the limits of currax, ensuring that the original limits are not exceeded
        x1 = max(x-zoomparams.xdist/2,zoomparams.oldxlim(1));
        y1 = max(y-zoomparams.ydist/2,zoomparams.oldylim(1));
        z1 = max(z-zoomparams.zdist/2,zoomparams.oldzlim(1));
        x2 = x1+zoomparams.xdist;
        y2 = y1+zoomparams.ydist;
        z2 = z1+zoomparams.zdist;
        % if new limits are out of range, adjust them:
        if x2 > zoomparams.oldxlim(2)
            x2 = zoomparams.oldxlim(2);
            x1 = x2 - zoomparams.xdist;
        end
        if y2 > zoomparams.oldylim(2)
            y2 = zoomparams.oldylim(2);
            y1 = y2 - zoomparams.ydist;
        end
        if z2 > zoomparams.oldzlim(2)
            z2 = zoomparams.oldzlim(2);
            z1 = z2 - zoomparams.zdist;
        end

        % now get the x,y positions in currax for display purposes
        posn = get(zoomparams.currax,'currentpoint');
        posn = posn(1,:);
        x = posn(1,1);
        y = posn(1,2);
        z = posn(1,3);
end

if x >= zoomparams.oldxlim(1) & x <= zoomparams.oldxlim(2) & ...
        y >= zoomparams.oldylim(1) & y <= zoomparams.oldylim(2) & ...
    z >= zoomparams.oldzlim(1) & z <= zoomparams.oldzlim(2)
    if strcmp(zoomparams.objtype,'surface')
        set(zoomparams.dispbox2,'string',sprintf('x = %3.2f;  y = %3.2f; z = %3.2f',x,y,z));
        set(zoomparams.currax,'xlim',[x1 x2],'ylim',[y1 y2],'zlim',[z1 z2]);
    else
        set(zoomparams.dispbox2,'string',sprintf('x = %3.2f;  y = %3.2f',x,y));
        set(zoomparams.currax,'xlim',[x1 x2],'ylim',[y1 y2]);
    end
else
    if strcmp(zoomparams.objtype,'surface')
        set(zoomparams.dispbox2,'string',sprintf('x = %3.f;  y = %3.0f; z = %3.0f',0,0,0));
        set(zoomparams.currax,'xlim',zoomparams.oldxlim,'ylim',zoomparams.oldylim,'zlim',zoomparams.oldzlim);
    else
        set(zoomparams.dispbox2,'string',sprintf('x = %3.0f;  y = %3.0f',0,0));
        set(zoomparams.currax,'xlim',zoomparams.oldxlim,'ylim',zoomparams.oldylim);
    end
end

%Note: up to this point, the only thing that has changed in refax is the currentpoint property
return

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function bdfcn
zoomparams = getappdata(gcf,'zoomparams');
% SelectionType
% normal: Click left mouse button
% extend: Shift - click left mouse button or click both left and right mouse buttons
% alt: Control - click left mouse button or click right mouse button
% open: Double click any mouse button

switch get(gcf,'selectiontype')
case 'normal'
	zoomparams.pct = max(0.01,zoomparams.pct*0.9);
case 'alt'
	zoomparams.pct = min(1,zoomparams.pct*1.1);
case 'extend'
	set(zoomparams.currax,'xlim',zoomparams.oldxlim,'ylim',zoomparams.oldylim,'zlim',zoomparams.oldzlim);
case 'open'
	zoomparams.pct = 1;
end

zoomparams.xdist = zoomparams.pct*zoomparams.xrange;
zoomparams.ydist = zoomparams.pct*zoomparams.yrange;
zoomparams.zdist = zoomparams.pct*zoomparams.zrange;

if ~strcmp(get(gcf,'selectiontype'),'extend')
    setappdata(gcf,'zoomparams',zoomparams);
    feval(getappdata(gcf,'zoomfcnhandle'));
end
return
 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%




% --- Executes on slider movement.
function slider6_Callback(hObject, eventdata, handles)
% hObject    handle to slider6 (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    structure with handles and user data (see GUIDATA)

% Hints: get(hObject,'Value') returns position of slider
%        get(hObject,'Min') and get(hObject,'Max') to determine range of slider


% --- Executes during object creation, after setting all properties.
function slider6_CreateFcn(hObject, eventdata, handles)
% hObject    handle to slider6 (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    empty - handles not created until after all CreateFcns called

% Hint: slider controls usually have a light gray background, change
%       'usewhitebg' to 0 to use default.  See ISPC and COMPUTER.
usewhitebg = 1;
if usewhitebg
    set(hObject,'BackgroundColor',[.9 .9 .9]);
else
    set(hObject,'BackgroundColor',get(0,'defaultUicontrolBackgroundColor'));
end




% --- Executes on slider movement.
function slider7_Callback(hObject, eventdata, handles)
% hObject    handle to slider7 (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    structure with handles and user data (see GUIDATA)
global pet;
input=pet;
 vals = get(handles.slider7,'Value');
 vals=vals/1000;
 set(handles.edit3,'String',vals)
 if(vals==0)
    
         vals=0.001;
 end      
     
  low_limit=vals;
high_limit=1-vals;
%CONTRAST_CONT Modulates the contrast of an image represented by its unit8
%matrix,'input' to give 'output'.
%
% NOTE:
% (1)The matrices 'output' and 'input' have class 'uint8'.
% (2)'low_limit' and 'high_limit' are the limiters.
% (3)'low_limit'<'high_limit' and both have values ranging between 0 and 1,except 0.

%--------------------------------------------------------------------------
img=input;
[m1 n1 r1]=size(img);
img2=double(img);
%--------------------calculation of vmin and vmax--------------------------
for(k=1:r1)
arr=sort(reshape(img2(:,:,k),m1*n1,1));
vmin(k)=arr(ceil(low_limit*m1*n1));
vmax(k)=arr(ceil(high_limit*m1*n1));
end
%--------------------------------------------------------------------------
if(r1==3)
    v_min=rgb2ntsc(vmin);
    v_max=rgb2ntsc(vmax);
end
if(r1==1)
    v_min=vmin;
    v_max=vmax;
end   
%--------------------------------------------------------------------------
for(i=1:m1)
     for(j=1:n1)
         for(k=1:r1)
                img2(i,j,k)=255*(img2(i,j,k)-v_min(1))/(v_max(1)-v_min(1));
         end
     end
end
%--------------------------------------------------------------------------
global output;
output=uint8(img2);
axes(handles.axes4);
%place image onto the axes
image(output);
%remove the axis tick marks
axis off
 global pet
 global result;
global temp
val=0.52789;
global outpu;
if exist('imlincomb')
    result = imlincomb(val, output, 1-val, outpu);
else
     result = val .* output + (1-val) .* outpu;
end
 
axes(handles.axes6);
%place image onto the axes
image(result);
%remove the axis tick marks
axis off
% Hints: get(hObject,'Value') returns position of slider
%        get(hObject,'Min') and get(hObject,'Max') to determine range of slider


% --- Executes during object creation, after setting all properties.
function slider7_CreateFcn(hObject, eventdata, handles)
% hObject    handle to slider7 (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    empty - handles not created until after all CreateFcns called

% Hint: slider controls usually have a light gray background, change
%       'usewhitebg' to 0 to use default.  See ISPC and COMPUTER.
usewhitebg = 1;
if usewhitebg
    set(hObject,'BackgroundColor',[.9 .9 .9]);
else
    set(hObject,'BackgroundColor',get(0,'defaultUicontrolBackgroundColor'));
end




% --- Executes on slider movement.
function slider8_Callback(hObject, eventdata, handles)
% hObject    handle to slider8 (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    structure with handles and user data (see GUIDATA)
 
 

 

global temp;
input=temp;
 vals = get(handles.slider8,'Value');

 vals=vals/1000;
  set(handles.edit2,'String',vals)
 if(vals==0)
    
         vals=0.001;
 end      
     
  low_limit=vals;
high_limit=1-vals;
%CONTRAST_CONT Modulates the contrast of an image represented by its unit8
%matrix,'input' to give 'output'.
%
% NOTE:
% (1)The matrices 'output' and 'input' have class 'uint8'.
% (2)'low_limit' and 'high_limit' are the limiters.
% (3)'low_limit'<'high_limit' and both have values ranging between 0 and 1,except 0.

%--------------------------------------------------------------------------
img=input;
[m1 n1 r1]=size(img);
img2=double(img);
%--------------------calculation of vmin and vmax--------------------------
for(k=1:r1)
arr=sort(reshape(img2(:,:,k),m1*n1,1));
vmin(k)=arr(ceil(low_limit*m1*n1));
vmax(k)=arr(ceil(high_limit*m1*n1));
end
%--------------------------------------------------------------------------
if(r1==3)
    v_min=rgb2ntsc(vmin);
    v_max=rgb2ntsc(vmax);
end
if(r1==1)
    v_min=vmin;
    v_max=vmax;
end   
%--------------------------------------------------------------------------
for(i=1:m1)
     for(j=1:n1)
         for(k=1:r1)
                img2(i,j,k)=255*(img2(i,j,k)-v_min(1))/(v_max(1)-v_min(1));
         end
     end
end
%--------------------------------------------------------------------------
global outpu;
outpu=uint8(img2);
axes(handles.axes5);
%place image onto the axes
image(outpu);
%remove the axis tick marks
axis off
 global pet