please help me for retina optic disc segmentation code

[getsanjeev]

@moorthysmart How do you expect us to help you? You mean by some paper or algorithm?

[moorthysmart]

no.. i have the code for optic disc cup detection and location.. so help me for it

On Tue 13 Mar, 2018, 2:40 PM Sanjeev Dubey, notifications@github.com wrote:

@moorthysmart https://github.com/moorthysmart How do you expect us to help you? You mean by some paper or algorithm?

— You are receiving this because you were mentioned. Reply to this email directly, view it on GitHub https://github.com/getsanjeev/retina-features/issues/6#issuecomment-372595232, or mute the thread https://github.com/notifications/unsubscribe-auth/APhUHZ1YUFRq1fN7FBp-cmbqTaCKM_u-ks5td4ysgaJpZM4SKQg4 .

[moorthysmart]

i have the code for only detection but crop and detection i needed/...

On Tue, Mar 13, 2018 at 9:47 PM, Moorthysmart M moorthy1778@gmail.com wrote:

no.. i have the code for optic disc cup detection and location.. so help me for it

On Tue 13 Mar, 2018, 2:40 PM Sanjeev Dubey, notifications@github.com wrote:

@moorthysmart https://github.com/moorthysmart How do you expect us to help you? You mean by some paper or algorithm?

— You are receiving this because you were mentioned. Reply to this email directly, view it on GitHub https://github.com/getsanjeev/retina-features/issues/6#issuecomment-372595232, or mute the thread https://github.com/notifications/unsubscribe-auth/APhUHZ1YUFRq1fN7FBp-cmbqTaCKM_u-ks5td4ysgaJpZM4SKQg4 .

[sidraaleem12]

Please help me with optic disc removal. Below is the code i tried following a paper tilted "Feature extraction and classification of retinal images for automated detection of Diabetic Retinopathy" link "https://ieeexplore.ieee.org/document/7802862/". The code is working fine, but not with all the images. By following mentioned paper I am creating meshgrid and then finding the x-y cordinates that corresponds to maximum intensity values (as optic disc has the highest intenstiy in retinal images). and then creating a circle of raius 45 around the brightest point. I used "diaretdb0 database" for testing. Out of 126 images, the result for 8 image is wrong. retina = imread('Input_Image.png'); retinaGband = retina(:,:,2); %Green Chanel extraction retinaGband = imresize(retinaGband,[512 512]); mask = imread('fmask.tif'); %mask mask = imresize(mask,size(retinaGband)); %% Histogram Equalization heImage = adapthisteq(retinaGband,'clipLimit',0.02,'Distribution','rayleigh'); % Mask the image. heImage = heImage .* uint8(mask); med2 = medfilt2(heImage); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%; [imageSizeX, imageSizeY] = size(med2); [x, y] = meshgrid(1:imageSizeX, 1:imageSizeY); %meshgrid creation [r, c] = find(med2 == max(med2(:))); %finding max values row and column wise centerX = max(r); %x-cordniate with max intensity value centerY = max(c); % y-cordinate with max intensity value radius = 45; circlePixels = (x - centerY).^2 ...

• (y - centerX).^2 <= radius.^2; circlePixels = imcomplement(circlePixels); %% masking mask1 = med2 .* uint8(circlePixels); figure,imshow(mask1);