How to remove optic disc from retinal images?

Question

Sidra Aleem 2018-4-3

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/392277-how-to-remove-optic-disc-from-retinal-images

回答： Image Analyst 2018-8-15

I am trying to do diabteic retinopathy grading. Before feature extraction, I have to remove optic disc. I am trying to do so using the paper titled "Automatic Optic Disk Segmentation in Presence of Disk Blurring" https://pdfs.semanticscholar.org/44cc/5980ecadaca668d2158d6bde6b6dd01efa63.pdf . They have used inverted circular LoG filter to enhance OD region. After the candidate circular regions have been enhanced using the LoG filter, they are binarized. The threshold that is used for this binarization is given by (T = 0.6 ∗ mLoG). mLoG is the maximum value in LoG ﬁltered image. Below is the code that I have written.

retina = imread('1.png');
retinaRband = retina(:,:,1); %Red Chanel extraction
retinaRband = imresize(retinaRband,[512 512]);
H= imcomplement(fspecial('log',[120 120],5));
H=(H-min(H(:)))/(max(H(:))-min(H(:)));
log_filtered = imfilter(retinaRband,H);
log_filtered=(log_filtered-min(log_filtered(:)))/(max(log_filtered(:))-min(log_filtered(:)));
max_log_filtered = max(log_filtered (:));
T = max_log_filtered * 0.6; %Threshold
Thresholded_image = imbinarize(log_filtered,0.6);
imshow(Thresholded_image);