一、简介
学相干层析技术(OCT)近年来快速发展,其中,眼前节组织 OCT 图像中中央角膜的边缘检测是急需解决的问题。本文根据传统形态学的边缘检测算法,先对所用的结构元素进行多尺度膨胀,再计算膨胀后各个图像的边缘检测结果的信息熵,最后用各个边缘检测结果的信息熵与相应图像的总信息熵之比作为图像融合的权值系数进行图像融合,得到最终的眼前节组织 OCT 图像的上下角膜边缘。实验结果表明,本文提出的算法能很好的适应眼前节组织OCT 图像的特征,边缘检测效果清晰,抗噪声能力强。
二、部分源码
% 清理工作空间
clc; clear all; close all;
% 检查预处理文件
if exist(fullfile(pwd, 'dt.mat'), 'file')
load(fullfile(pwd, 'dt.mat'));
figure; imshow(result, []);
return;
end
% 读取图片
Img = imread('images\image.bmp');
% 计算1~5的算子结果
Gf1 = Main_Process(Img, 1);
Gf2 = Main_Process(Img, 2);
Gf3 = Main_Process(Img, 3);
Gf4 = Main_Process(Img, 4);
Gf5 = Main_Process(Img, 5);
% 整合到cell
G{1} = Gf1;
G{2} = Gf2;
G{3} = Gf3;
G{4} = Gf4;
G{5} = Gf5;
% 计算系数
ua1 = Coef(Gf1, G);
ua2 = Coef(Gf2, G);
ua3 = Coef(Gf3, G);
ua4 = Coef(Gf4, G);
ua5 = Coef(Gf5, G);
% 组合
u = [ua1, ua2, ua3, ua4, ua5];
% 权重
u = u/sum(u);
% 加权融合
Gf = Edge_One(G, u);
% 显示结果
result = Gf5;
figure; imshow(Img, []);
figure; imshow(result, []);
三、运行结果
四、matlab版本
MATLAB R2019a
五、学习与交流
文中不足之处请大家多多指教,学习与交流可以联系企鹅:3752243968
文中部分源码仅供参考,若需要全部代码可以私信(有偿)
全部代码可以私信(有偿)