清空筛选条件
清空筛选条件

Error while evaluationg UIcontrolm Callback.

3 次查看（过去 30 天）
khushboo chhikara
khushboo chhikara 2018-4-20
评论： khushboo chhikara 2018-4-21
function varargout = Glaucoma_Detection(varargin)
% GLAUCOMA_DETECTION MATLAB code for Glaucoma_Detection.fig
% GLAUCOMA_DETECTION, by itself, creates a new GLAUCOMA_DETECTION or raises the existing
% singleton*.
%
% H = GLAUCOMA_DETECTION returns the handle to a new GLAUCOMA_DETECTION or the handle to
% the existing singleton*.
%
% GLAUCOMA_DETECTION('CALLBACK',hObject,eventData,handles,...) calls the local
% function named CALLBACK in GLAUCOMA_DETECTION.M with the given input arguments.
%
% GLAUCOMA_DETECTION('Property','Value',...) creates a new GLAUCOMA_DETECTION or raises the
% existing singleton*. Starting from the left, property value pairs are
% applied to the GUI before Glaucoma_Detection_OpeningFcn gets called. An
% unrecognized property name or invalid value makes property application
% stop. All inputs are passed to Glaucoma_Detection_OpeningFcn via varargin.
%
% *See GUI Options on GUIDE's Tools menu. Choose "GUI allows only one
% instance to run (singleton)".
%
% See also: GUIDE, GUIDATA, GUIHANDLES
% Edit the above text to modify the response to help Glaucoma_Detection
% Last Modified by GUIDE v2.5 19-Apr-2018 10:21:01
% Begin initialization code - DO NOT EDIT
gui_Singleton = 1;
gui_State = struct('gui_Name', mfilename, ...
'gui_Singleton', gui_Singleton, ...
'gui_OpeningFcn', @Glaucoma_Detection_OpeningFcn, ...
'gui_OutputFcn', @Glaucoma_Detection_OutputFcn, ...
'gui_LayoutFcn', [] , ...
'gui_Callback', []);
if nargin && ischar(varargin{1})
gui_State.gui_Callback = str2func(varargin{1});
end
if nargout
[varargout{1:nargout}] = gui_mainfcn(gui_State, varargin{:});
else
gui_mainfcn(gui_State, varargin{:});
end
% End initialization code - DO NOT EDIT
% --- Executes just before Glaucoma_Detection is made visible.
function Glaucoma_Detection_OpeningFcn(hObject, eventdata, handles, varargin)
% This function has no output args, see OutputFcn.
% hObject handle to figure
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
% varargin command line arguments to Glaucoma_Detection (see VARARGIN)
% Choose default command line output for Glaucoma_Detection
handles.output = hObject;
% Update handles structure
guidata(hObject, handles);
% UIWAIT makes Glaucoma_Detection wait for user response (see UIRESUME)
% uiwait(handles.figure1);
% Choose default command line output for testfinal
handles.output = hObject;
a=ones(256,256);
b=imread('images.jpeg');
axes(handles.axes3);
imshow(b);
axes(handles.axes2);
imshow(a);
axes(handles.axes1);
imshow(a);
%Update handles structure
guidata(hObject, handles);
% --- Outputs from this function are returned to the command line.
function varargout = Glaucoma_Detection_OutputFcn(hObject, eventdata, handles)
% varargout cell array for returning output args (see VARARGOUT);
% hObject handle to figure
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
% Get default command line output from handles structure
varargout{1} = handles.output;
function ANGLE_Callback(hObject, eventdata, handles)
% hObject handle to ANGLE (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,'String') returns contents of ANGLE as text
% str2double(get(hObject,'String')) returns contents of ANGLE as a double
% --- Executes during object creation, after setting all properties.
function ANGLE_CreateFcn(hObject, eventdata, handles)
% hObject handle to ANGLE (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles empty - handles not created until after all CreateFcns called
% Hint: edit controls usually have a white background on Windows.
% See ISPC and COMPUTER.
if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor'))
set(hObject,'BackgroundColor','white');
end
function STATUS1_Callback(hObject, eventdata, handles)
% hObject handle to STATUS1 (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,'String') returns contents of STATUS1 as text
% str2double(get(hObject,'String')) returns contents of STATUS1 as a double
% --- Executes during object creation, after setting all properties.
function STATUS1_CreateFcn(hObject, eventdata, handles)
% hObject handle to STATUS1 (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles empty - handles not created until after all CreateFcns called
% Hint: edit controls usually have a white background on Windows.
% See ISPC and COMPUTER.
if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor'))
set(hObject,'BackgroundColor','white');
end
% --- Executes on button press in delete.
function delete_Callback(hObject, eventdata, handles)
% hObject handle to delete (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
[file, pathname] = uigetfile('*.bmp', 'Pick an Image');
if isequal(file,0) | isequal(pathname,0)
warndlg('User pressed cancel')
else
a=file;
delete(a);
uiwait(msgbox('File Deleted Succesfully','Message','modal'));
end
% --- Executes on button press in clear.
function clear_Callback(hObject, eventdata, handles)
% hObject handle to clear (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
a=ones(256,256);
axes(handles.axes1);
imshow(a);
axes(handles.axes2);
imshow(a);
bb='0'
set(handles.ANGLE,'String',bb);
bb='0'
set(handles.STATUS1,'String',bb);
% --- Executes on button press in help.
function help_Callback(hObject, eventdata, handles)
% hObject handle to help (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
helpwin ch1;
% --- Executes on button press in exist.
function exist_Callback(hObject, eventdata, handles)
% hObject handle to exist (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
exit;
% --- Executes on button press in load.
function load_Callback(hObject, eventdata, handles)
% hObject handle to load (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
[file, pathname] = uigetfile('*.jpeg', 'Pick an Image');
if isequal(file,0) | isequal(pathname,0)
warndlg('User pressed cancel')
else
a=imread(file);
axes(handles.axes1);
imshow(a);
handles.file=file;
%Update handles structure
guidata(hObject, handles);
end
% --- Executes on button press in enhacement.
function enhacement_Callback(hObject, eventdata, handles)
% hObject handle to enhacement (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
file = handles.file;
a =imread(file);
[r c p]=size(a); % If it is color, convert to gray....
if p==3
a=rgb2gray(a);
end
b=histeq(a); % Applying the Histogram
imwrite(b,'enhanced.bmp');
axes(handles.axes2);
imshow(b);
% --- Executes on button press in thresh.
function thresh_Callback(hObject, eventdata, handles)
% hObject handle to thresh (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
file = handles.file;
a=imread(file);
[r c p]=size(a);
b=a;
if p==3
b=rgb2gray(a);
end
Out =zeros(r,c);
for i=1:r
for j=1:c
if b(i,j) >150
Out(i,j)=1;
else
Out(i,j)=0;
end
end
end
axes(handles.axes2);
imshow(Out);
imwrite(Out,'thresholded.bmp');
% --- Executes on button press in erosion.
function erosion_Callback(hObject, eventdata, handles)
% hObject handle to erosion (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
file = handles.file;
Out=imread(file);
B1 = medfilt2(Out,[5 5]);
axes(handles.axes2);
imshow(B1);
imwrite(B1,'eroded.bmp');
% --- Executes on button press in holefill.
function holefill_Callback(hObject, eventdata, handles)
% hObject handle to holefill (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
file = handles.file;
bw1=imread(file);
BW15 = imfill(bw1,8,'holes');
axes(handles.axes2);
imshow(BW15);
BW15=uint8(BW15);
imwrite(BW15,'FILLED.bmp');
uiwait(msgbox('Process completed','Hole filling Status','modal'));
% --- Executes on button press in segment.
function segment_Callback(hObject, eventdata, handles)
% hObject handle to segment (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
file = handles.file;
file=imread(file);
J = imcomplement(file);
figure(1);
imshow(J);
J=bwmorph(J,'FILL');
figure(2);
imshow(J);
[r1 c1]=size(J); %%Here B&W image
[L2 NUM2] = bwlabel(J);
disp(NUM2);
removed=0;
STATS = regionprops(L2,'area'); %Finds the number of pixels on each region
aA=[STATS.Area]; %Conversion from structure to array
for j=1:1:NUM2
bw= aA(j);
if bw <1500
L2(L2==j)=0;
removed = removed + 1;
end
end
NUM2=NUM2-removed;
[L3 NUM3] = bwlabel(L2);
STATS = regionprops(L3,'area');
disp(NUM3);
stats1 = regionprops(L3, 'PixelList'); %getting the pixel co-ordinates
CC=[];
for j=1:1:NUM3 %NUM3
aA1=stats1(j);
CC{j} = [aA1.PixelList]; %Conversion from structure to array
end
removed1=0;
tmp = [];
for j=1:NUM3 % THIS LOOP FORMS THE TRI_ANGLE
tmp = CC{j}(:);
B=find(tmp(:,1)==1);%%%%%%%%%
disp(B);
if isempty(B)
imshow(L3==j);
else
L3(L3==j)=0;
removed1=removed1+1;
end
end
NUM3=NUM3-removed1;
disp(NUM3);
[L4 NUM4] = bwlabel(L3);
for i=1:NUM4
[r,c] = find(L4 == i);
BW2 = bwselect(L3,c,r,4);
imwrite(BW2,'SEGMENTED.BMP');
end
% figure(7);
% imshow(L4);
if NUM4~=0
Image=BW2;
th1=1;
th2=0;
[rows columns]=size(Image);
len=rows;
Rimage = Image;
OutputImage=ones(rows,columns);
Boundary_Pixel=[];
%detecting ones
h = waitbar(0,'Please wait...');
index=1;
for y= 1:columns
for x=1:rows
%for left top corner
if (x==1 & y==1)
if (Rimage(x, y+1)==th1)
if (Rimage(x+1, y+1)==th1)
if (Rimage(x+1,y)==th1)
Boundary_Pixel{index}=[x,y];
OutputImage(x,y)=0;
index=index+1;
end
end
end
%for right top corner
elseif (x==1 & y==len)
if (Rimage(x+1, y)==th1)
if (Rimage(x+1, y-1)==th1)
if (Rimage(x,y-1)==th1)
Boundary_Pixel{index}=[x,y];
OutputImage(x,y)=0;
index=index+1;
end
end
end
%for bottom left corner[2,3,4]
elseif (x==len & y==1)
if (Rimage(x-1, y)==th1)
if (Rimage(x-1, y+1) ==th1)
if (Rimage(x,y+1) ==th1)
Boundary_Pixel{index}=[x,y];
OutputImage(x,y)=0;
index=index+1;
end
end
end
%for bottom right corner[8,1,2]
elseif (x==len & y==len)
if (Rimage(x, y-1)==th1)
if (Rimage(x-1, y-1)==th1)
if (Rimage(x-1,y)==th1)
Boundary_Pixel{index}=[x,y];
OutputImage(x,y)=0;
index=index+1;
end
end
end
%for top edge[8,7,6,5,4]
elseif (x==1)
if (Rimage(x, y-1)==th1)
if (Rimage(x+1, y-1)==th1)
if (Rimage(x+1, y)==th1)
if (Rimage(x+1, y+1)==th1)
if (Rimage(x,y+1)==th1)
Boundary_Pixel{index}=[x,y];
OutputImage(x,y)=0;
index=index+1;
end
end
end
end
end
%for right edge[2,1,8,7,6]
elseif (y==len)
if (Rimage(x-1, y)==th1)
if (Rimage(x-1, y-1)==th1)
if (Rimage(x, y-1)==th1)
if (Rimage(x+1, y-1)==th1)
if (Rimage(x+1,y)==th1)
Boundary_Pixel{index}=[x,y];
OutputImage(x,y)=0;
index=index+1;
end
end
end
end
end
%for bottom edge[8,1,2,3,4]
elseif (x==len)
if (Rimage(x, y-1)==th1)
if (Rimage(x-1, y-1)==th1)
if (Rimage(x-1, y)==th1)
if (Rimage(x-1, y+1)==th1)
if (Rimage(x,y+1)==th1)
Boundary_Pixel{index}=[x,y];
OutputImage(x,y)=0;
index=index+1;
end
end
end
end
end
%for left edge[2,3,4,5,6]
elseif (y==1)
if (Rimage(x-1, y)==th1)
if (Rimage(x-1, y+1)==th1)
if (Rimage(x, y+1)==th1)
if (Rimage(x+1, y+1)==th1)
if (Rimage(x+1,y)==th1)
Boundary_Pixel{index}=[x,y];
OutputImage(x,y)=0;
index=index+1;
end
end
end
end
end
%for remaining pixels
else
if (Rimage(x-1,y-1)==th1)
if (Rimage(x-1,y)==th1)
if (Rimage(x-1,y+1)==th1)
if (Rimage(x,y+1)==th1)
if (Rimage(x+1,y+1)==th1)
if (Rimage(x+1,y)==th1)
if (Rimage(x+1,y-1)==th1)
if (Rimage(x,y-1)==th1)
Boundary_Pixel{index}=[x,y];
OutputImage(x,y)=0;
index=index+1;
end
end
end
end
end
end
end
end
end
end
end
% figure(221);
% imshow(OutputImage);
%detecting zeros
index=1;
for y=1:columns
for x=1:rows
%for left top corner[4,5,6]
if (x==1 & y==1)
if (Rimage(x, y+1)==th2)
if (Rimage(x+1, y+1)==th2)
if (Rimage(x+1,y)==th2)
Boundary_Pixel{index}=[x,y];
OutputImage(x,y)=0;
index=index+1;
end
end
end
%for right top corner[8,7,6]
elseif (x==1 & y==len)
if (Rimage(x+1, y)==th2)
if (Rimage(x+1, y-1)==th2)
if (Rimage(x,y-1)==th2)
Boundary_Pixel{index}=[x,y];
OutputImage(x,y)=0;
index=index+1;
end
end
end
%for bottom left corner[2,3,4]
elseif (x==len & y==1)
if (Rimage(x-1, y)==th2)
if (Rimage(x-1, y+1) ==th2)
if (Rimage(x,y+1) ==th2)
Boundary_Pixel{index}=[x,y];
OutputImage(x,y)=0;
index=index+1;
end
end
end
%for bottom right corner[8,1,2]
elseif (x==len & y==len)
if (Rimage(x, y-1)==th2)
if (Rimage(x-1, y-1)==th2)
if (Rimage(x-1,y)==th2)
Boundary_Pixel{index}=[x,y];
OutputImage(x,y)=0;
index=index+1;
end
end
end
%for top edge[8,7,6,5,4]
elseif (x==1)
if (Rimage(x, y-1)==th2)
if (Rimage(x+1, y-1)==th2)
if (Rimage(x+1, y)==th2)
if (Rimage(x+1, y+1)==th2)
if (Rimage(x,y+1)==th2)
Boundary_Pixel{index}=[x,y];
OutputImage(x,y)=0;
index=index+1;
end
end
end
end
end
%for right edge[2,1,8,7,6]
elseif (y==len)
if (Rimage(x-1, y)==th2)
if (Rimage(x-1, y-1)==th2)
if (Rimage(x, y-1)==th2)
if (Rimage(x+1, y-1)==th2)
if (Rimage(x+1,y)==th2)
Boundary_Pixel{index}=[x,y];
OutputImage(x,y)=0;
index=index+1;
end
end
end
end
end
%for bottom edge[8,1,2,3,4]
elseif (x==len)
if (Rimage(x, y-1)==th2)
if (Rimage(x-1, y-1)==th2)
if (Rimage(x-1, y)==th2)
if (Rimage(x-1, y+1)==th2)
if (Rimage(x,y+1)==th2)
Boundary_Pixel{index}=[x,y];
OutputImage(x,y)=0;
index=index+1;
end
end
end
end
end
%for left edge[2,3,4,5,6]
elseif (y==1)
if (Rimage(x-1, y)==th2)
if (Rimage(x-1, y+1)==th2)
if (Rimage(x, y+1)==th2)
if (Rimage(x+1, y+1)==th2)
if (Rimage(x+1,y)==th2)
Boundary_Pixel{index}=[x,y];
OutputImage(x,y)=0;
index=index+1;
end
end
end
end
end
%for remaining pixels
else
if (Rimage(x-1,y-1)==th2)
if (Rimage(x-1,y)==th2)
if (Rimage(x-1,y+1)==th2)
if (Rimage(x,y+1)==th2)
if (Rimage(x+1,y+1)==th2)
if (Rimage(x+1,y)==th2)
if (Rimage(x+1,y-1)==th2)
if (Rimage(x,y-1)==th2)
Boundary_Pixel{index}=[x,y];
OutputImage(x,y)=0;
index=index+1;
end
end
end
end
end
end
end
end
end
waitbar(y/columns,h)
end
end
close(h);
% figure(224);
% imshow(OutputImage);
% DETECTION OF BOUNDARY DETECTED IMAGE
% imshow(OutputImage)
imwrite(OutputImage,'Bound1.bmp');
% save a
% load a
Out1=OutputImage;
[m n]=size(Out1);
Out2=zeros(m,n);
for i=1:m
for j=1:n
if Out1(i,j)>=1
Out2(i,j)=255;
end
end
end
index=1;
for x = 1:256
for y=1:256
bb=Out2(x,y);
if bb>=10
Boundary_Pixel{index}=[x,y];
X1(index)=x;
Y1(index)=y;
index=index+1;
end
end
end
% save a1
% load a1
X = X1; %X co-ordinates
Y =Y1; %Y co-ordinates
% XAA = X1;
% YAA =Y1;
save X X
save Y Y
X1=X;
Y1 =Y;
XB=X;
YB =Y;
Xmin = min(XB)
[Ymin Ind] = min(YB) %FIND OUT THE CORNER FIRST POINTS OF TRI_ANGLE
Dec=X(Ind)
Xmax = max(XB);
Ymax = max(YB);
[Ymax Ind] = max(YB);
Dec1=X(Ind);
axes(handles.axes2);
% X = [Xmin Xmax];
X = [Dec Dec1];
Y = [Ymin Ymax];
line(Y,X) ;
X = [Dec Dec];
Y = [Ymin Ymax];
line(Y,X) ;
X = [Dec Dec1];
Y = [Ymax Ymax];
line(Y,X) ;
X = [Dec Xmax];
Y = [Ymax Ymax];
line(Y,X) ;
X = [Dec Xmax];
Y = [Ymin Ymax];
b=ones(256,256);
line(Y,X) ;
handles.Xmin=Xmin;
handles.Xmax=Xmax;
handles.Dec=Dec;
handles.Dec1=Dec1;
handles.Ymin=Ymin;
handles.Ymax=Ymax;
else
bb='0'
set(handles.ANGLE,'String',bb);
warndlg('anterior chamber is almost closed ... glaucoma detected..no further processing');
bb='GLAUCOMA DETECTED'
set(handles.STATUS1,'String',bb);
disp('no anterior chamber');
end
% Update handles structure
guidata(hObject, handles);
% --- Executes on button press in mask.
function mask_Callback(hObject, eventdata, handles)
% hObject handle to mask (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
Xmin=get(handles.Xmin);
Xmax=get(handles.Xmax);
Dec =handles.Dec;
Dec1 =handles.Dec1;
Ymin =handles.Ymin;
Ymax =handles.Ymax;
file = handles.file;
aa=imread(file);
%%%%%%%%%%%%%%%%%%%CALCUATING BOUNDNG BOX
X1=Ymin;
Y1=Dec;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
LTX1=X1-25;
LTY1=Y1-25;
RTX1=X1-25;
RTY1=Y1+25;
LBX1=X1+25;
LBY1=Y1-25;
RBX1=X1+25;
RBY1=Y1+25;
ROIIMAGE11=aa(LTY1:RBY1,LTX1:LBX1);
axes(handles.axes2);
imshow(ROIIMAGE11);
imwrite(ROIIMAGE11,'masked.bmp');
% --- Executes on button press in aod.
function aod_Callback(hObject, eventdata, handles)
% hObject handle to aod (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
file=handles.file;
a=imread(file);
a=uint8(a);
b=histeq(a);
figure(2);
imshow(b);
%pixval on;
[r c]=size(a);
for i=1:r
for j=1:c
if b(i,j) >50
aOut2(i,j)=1;
aOut21(i,j)=0;
else
aOut2(i,j)=0;
aOut21(i,j)=1;
end
end
end
aOut22 = imcomplement(imfill(imcomplement(aOut2),'holes'));
% imshow(aOut22);
aOut22 = imfill(aOut22,'holes');
aJ2 = imcomplement(aOut22);
% INVERTED IMAGE
[L2 NUM2] = bwlabel(aJ2);
% disp(NUM2);
removed=0;
removed1=0;
STATS = regionprops(L2,'area');
%STATS1 = imfeature(L2,'all');
aA=[STATS.Area];
for j=1:1:NUM2
bw= aA(j);
if bw <150
L2(L2==j)=0;
removed = removed + 1;
end
end
NUM2=NUM2-removed;
[L3 NUM3] = bwlabel(L2);
imshow(L3==1);
[r,c] = find(L3 == 1);
% figure;
BW2 = bwselect(L3,c,r,4);
% a=imread(filename);
Image=BW2;
imshow(Image);
% pixval on;
th1=1;
th2=0;
[rows columns]=size(Image);
len=rows;
Rimage = Image;
j=columns;
Num=[];
index1=1;
for i=1:r
if Image(i,j)==1
Boundary_Pixel_1{index1}=[i,j];
Num(index1)=i;
index1=index1+1;
end
end
Mini=min(Num);
Maxi=max(Num);
OutputImage=ones(rows,columns);
Boundary_Pixel=[];
% Boundary_y=[];
% detecting ones
h = waitbar(0,'Please wait...');
index=1;
for y= 1:columns
for x=1:rows
% for left top corner[4,5,6]
if (x==1 & y==1)
if (Rimage(x, y+1)==th1)
if (Rimage(x+1, y+1)==th1)
if (Rimage(x+1,y)==th1)
Boundary_Pixel{index}=[x,y];
OutputImage(x,y)=0;
index=index+1;
end
end
end
% for right top corner[8,7,6]
elseif (x==1 & y==len)
if (Rimage(x+1, y)==th1)
if (Rimage(x+1, y-1)==th1)
if (Rimage(x,y-1)==th1)
Boundary_Pixel{index}=[x,y];
OutputImage(x,y)=0;
index=index+1;
end
end
end
% for bottom left corner[2,3,4]
elseif (x==len & y==1)
if (Rimage(x-1, y)==th1)
if (Rimage(x-1, y+1) ==th1)
if (Rimage(x,y+1) ==th1)
Boundary_Pixel{index}=[x,y];
OutputImage(x,y)=0;
index=index+1;
end
end
end
% for bottom right corner[8,1,2]
elseif (x==len & y==len)
if (Rimage(x, y-1)==th1)
if (Rimage(x-1, y-1)==th1)
if (Rimage(x-1,y)==th1)
Boundary_Pixel{index}=[x,y];
OutputImage(x,y)=0;
index=index+1;
end
end
end
% for top edge[8,7,6,5,4]
elseif (x==1)
if (Rimage(x, y-1)==th1)
if (Rimage(x+1, y-1)==th1)
if (Rimage(x+1, y)==th1)
if (Rimage(x+1, y+1)==th1)
if (Rimage(x,y+1)==th1)
Boundary_Pixel{index}=[x,y];
OutputImage(x,y)=0;
index=index+1;
end
end
end
end
end
% for right edge[2,1,8,7,6]
elseif (y==len)
if (Rimage(x-1, y)==th1)
if (Rimage(x-1, y-1)==th1)
if (Rimage(x, y-1)==th1)
if (Rimage(x+1, y-1)==th1)
if (Rimage(x+1,y)==th1)
Boundary_Pixel{index}=[x,y];
OutputImage(x,y)=0;
index=index+1;
end
end
end
end
end
% for bottom edge[8,1,2,3,4]
elseif (x==len)
if (Rimage(x, y-1)==th1)
if (Rimage(x-1, y-1)==th1)
if (Rimage(x-1, y)==th1)
if (Rimage(x-1, y+1)==th1)
if (Rimage(x,y+1)==th1)
Boundary_Pixel{index}=[x,y];
OutputImage(x,y)=0;
index=index+1;
end
end
end
end
end
% for left edge[2,3,4,5,6]
elseif (y==1)
if (Rimage(x-1, y)==th1)
if (Rimage(x-1, y+1)==th1)
if (Rimage(x, y+1)==th1)
if (Rimage(x+1, y+1)==th1)
if (Rimage(x+1,y)==th1)
Boundary_Pixel{index}=[x,y];
OutputImage(x,y)=0;
index=index+1;
end
end
end
end
end
else
% for remaining pixels
if (Rimage(x-1,y-1)==th1)
if (Rimage(x-1,y)==th1)
if (Rimage(x-1,y+1)==th1)
if (Rimage(x,y+1)==th1)
if (Rimage(x+1,y+1)==th1)
if (Rimage(x+1,y)==th1)
if (Rimage(x+1,y-1)==th1)
if (Rimage(x,y-1)==th1)
Boundary_Pixel{index}=[x,y];
OutputImage(x,y)=0;
index=index+1;
end
end
end
end
end
end
end
end
end
end
end
%detecting zeros
index=1;
for y=1:columns
for x=1:rows
% for left top corner[4,5,6]
if (x==1 & y==1)
if (Rimage(x, y+1)==th2)
if (Rimage(x+1, y+1)==th2)
if (Rimage(x+1,y)==th2)
Boundary_Pixel{index}=[x,y];
OutputImage(x,y)=0;
index=index+1;
end
end
end
% for right top corner[8,7,6]
elseif (x==1 & y==len)
if (Rimage(x+1, y)==th2)
if (Rimage(x+1, y-1)==th2)
if (Rimage(x,y-1)==th2)
Boundary_Pixel{index}=[x,y];
OutputImage(x,y)=0;
index=index+1;
end
end
end
% for bottom left corner[2,3,4]
elseif (x==len & y==1)
if (Rimage(x-1, y)==th2)
if (Rimage(x-1, y+1) ==th2)
if (Rimage(x,y+1) ==th2)
Boundary_Pixel{index}=[x,y];
OutputImage(x,y)=0;
index=index+1;
end
end
end
% for bottom right corner[8,1,2]
elseif (x==len & y==len)
if (Rimage(x, y-1)==th2)
if (Rimage(x-1, y-1)==th2)
if (Rimage(x-1,y)==th2)
Boundary_Pixel{index}=[x,y];
OutputImage(x,y)=0;
index=index+1;
end
end
end
% for top edge[8,7,6,5,4]
elseif (x==1)
if (Rimage(x, y-1)==th2)
if (Rimage(x+1, y-1)==th2)
if (Rimage(x+1, y)==th2)
if (Rimage(x+1, y+1)==th2)
if (Rimage(x,y+1)==th2)
Boundary_Pixel{index}=[x,y];
OutputImage(x,y)=0;
index=index+1;
end
end
end
end
end
% for right edge[2,1,8,7,6]
elseif (y==len)
if (Rimage(x-1, y)==th2)
if (Rimage(x-1, y-1)==th2)
if (Rimage(x, y-1)==th2)
if (Rimage(x+1, y-1)==th2)
if (Rimage(x+1,y)==th2)
Boundary_Pixel{index}=[x,y];
OutputImage(x,y)=0;
index=index+1;
end
end
end
end
end
% for bottom edge[8,1,2,3,4]
elseif (x==len)
if (Rimage(x, y-1)==th2)
if (Rimage(x-1, y-1)==th2)
if (Rimage(x-1, y)==th2)
if (Rimage(x-1, y+1)==th2)
if (Rimage(x,y+1)==th2)
Boundary_Pixel{index}=[x,y];
OutputImage(x,y)=0;
index=index+1;
end
end
end
end
end
% for left edge[2,3,4,5,6]
elseif (y==1)
if (Rimage(x-1, y)==th2)
if (Rimage(x-1, y+1)==th2)
if (Rimage(x, y+1)==th2)
if (Rimage(x+1, y+1)==th2)
if (Rimage(x+1,y)==th2)
Boundary_Pixel{index}=[x,y];
OutputImage(x,y)=0;
index=index+1;
end
end
end
end
end
else
% for remaining pixels
if (Rimage(x-1,y-1)==th2)
if (Rimage(x-1,y)==th2)
if (Rimage(x-1,y+1)==th2)
if (Rimage(x,y+1)==th2)
if (Rimage(x+1,y+1)==th2)
if (Rimage(x+1,y)==th2)
if (Rimage(x+1,y-1)==th2)
if (Rimage(x,y-1)==th2)
Boundary_Pixel{index}=[x,y];
OutputImage(x,y)=0;
index=index+1;
end
end
end
end
end
end
end
end
end
waitbar(y/columns,h)
end
end
close(h);
% DETECTION OF BOUNDARY DETECTED IMAGE
% imshow(OutputImage)
% imwrite(OutputImage,'Bound1.bmp')
Out1=OutputImage;
[m n]=size(Out1);
% pixval on
Out2=zeros(m,n);
for i=1:m
for j=1:n
if Out1(i,j)>=1
Out2(i,j)=255;
end
end
end
% imshow(Out2,[]);
% pixval on
index=1;
threshold=255;
for x = 1:rows
for y=1:columns
bb=Out2(x,y);
if bb>=10
Boundary_Pixel{index}=[x,y];
X1(index)=x;
Y1(index)=y;
index=index+1;
end
end
end
save a11
load a11
%
% X = X1;
% Y =Y1;
% XAA = X1;
% YAA =Y1;
%
% save XAA XAA
% save YAA YAA
%
X1=X;
Y1 =Y;
XB=X;
YB =Y;
Xmin = min(XB);
[Ymin Ind] = min(YB);
Dec=X(Ind);
Xmax = max(XB);
Ymax = max(YB);
[Ymax Ind] = max(YB);
Dec1=X(Ind);
aa=imread(file);
axes(handles.axes2);
imshow(aa);
% X = [Xmin Xmax];
X = [Dec Mini];
Y = [Ymin Ymax];
line(Y,X) ;
X = [Dec Dec];
Y = [Ymin Ymax];
line(Y,X) ;
X = [Dec Dec1];
Y = [Ymax Ymax];
line(Y,X) ;
X = [Dec Xmax];
Y = [Ymax Ymax];
line(Y,X) ;
X = [Dec Maxi];
Y = [Ymin Ymax];
b=ones(256,256);
line(Y,X) ;
% pixval on
Distance1=sqrt((Dec-Mini)^2+(Ymin-Ymax)^2);%opposite of traingle 1hyp
disp(Distance1);
Distance2=sqrt((Dec-Maxi)^2+(Ymin-Ymax)^2);%opposite of traingle 1hyp2
disp(Distance2);
Distance3=sqrt((Dec-Dec)^2+(Ymin-Ymax)^2);%adjacet of traingle 1
disp(Distance3);
%tan=opp/adj
Thete1=acos(Distance3/Distance1);
Thete2=acos(Distance3/Distance2);
AA=Thete1+Thete2;
BB=(180/pi)*AA;
% Thete=atan(Distance1/Distance3);
% Theta=acos(adjacent/hyp);
Theta1=(180/pi)*Thete1;
Theta2=(180/pi)*Thete2;
disp(Theta1);
disp(Theta2);
% Theta=Theta1-Theta2;
Theta3=Theta1-Theta2;
Theta13=num2str(Theta3);
set(handles.ANGLE,'String',Theta13);
% disp(Theta);
Theta3=round(Theta3);
if Theta3 > 19
bb='NORMAL EYE'
set(handles.STATUS1,'String',bb);
else
bb='GLAUCOMA DETECTED'
set(handles.STATUS1,'String',bb);
end
% --- Executes on button press in radian.
function radian_Callback(hObject, eventdata, handles)
% hObject handle to radian (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
% Hint: get(hObject,'Value') returns toggle state of radian
% --- Executes on button press in togglebutton2.
function togglebutton2_Callback(hObject, eventdata, handles)
% hObject handle to togglebutton2 (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
% Hint: get(hObject,'Value') returns toggle state of togglebutton2
button_state = get(hObject,'Value');
if button_state == get(hObject,'Max')
display('clicked')
elseif button_state == get(hObject,'Min')
display('rleased');
end
After running this code .. i am getting some errors. Please help. The errors are:
  • Error in gui_mainfcn (line 95)
  • feval(varargin{:});
  • Error in Glaucoma_Detection (line 42)
  • gui_mainfcn(gui_State, varargin{:});
  • Error in Glaucoma_Detection>mask_Callback (line 779)
  • Xmin=get(handles.Xmin);

请先登录，再回答此问题。

采纳的回答

Walter Roberson
Walter Roberson 2018-4-20

You have

function erosion_Callback(hObject, eventdata, handles)
  % hObject    handle to erosion (see GCBO)
  % eventdata  reserved - to be defined in a future version of MATLAB
  % handles    structure with handles and user data (see GUIDATA)
  file = handles.file;
  Out=imread(file);
  B1 = medfilt2(Out,[5 5]);

that code relies upon the result of the imread() being 2 dimensional. That is not impossible, but it does rely upon the user having selected a grayscale file. If the user happened to select a jpeg image, then that would rarely be grayscale (not impossible, but I have only ever seen one grayscale jpeg file "in real life").

Your image enhancement callback has code to adjust for this, so you already know how to do this, just copy the technique from there.

segment_Callback has the same problem.

Your Masking routine will fail if the segmentation routine has not been completed, because it relies upon handles.Xmin having been set. This is a more general problem in your code: most of your routines will fail if you have not Load, for example.

There is a truism in writing user interfaces: if an action is not valid because the user has not done something yet, then the user should not be able to select the action.

In other words, you should be configuring the state of all of the pushbuttons along the left side to 'Enable', 'off', and you should only set them to 'Enable', 'on' once they become meaningful by the user having completed all previous steps that are relevant. (And remember to disable them again when the user loads a new image.)

  6 个评论
显示 4更早的评论
Walter Roberson
Walter Roberson 2018-4-20
Updated / cleaned up file attached.
You need to repair the two lines marked %FIXME as the bug in your code prevents AOD from working.
I have doubts that your code works properly. At each stage you are going back to the original image. It is much more likely that you want to process the result of the previous stage.

请先登录，再进行评论。

更多回答（0 个）

类别

在 Help Center 和 File Exchange 中查找有关 Graphics Object Properties 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by