Matlab 高斯_拉普拉斯滤波器处理医学图像
前言:本程序是我去年实现论文算法时所做。主要功能为标记切割肝脏区域。时间有点久,很多细节已经模糊加上代码做了很多注释,因此在博客中不再详述。
NOTE: 程序分几大段功能模块,仔细阅读,对解决医学图像还是有一定的借鉴意义
想借鉴本文的一定要仔细阅读代码和注释,中间有人机交互部分,空跑会抛异常
.dcm数据,我放到了我的百度云盘,有兴趣的可以下载,实测一下代码。dcm数据连接
clc,clear
img_1=dicomread('10011.dcm');%读取dcm文件 (所谓的灰度值)
metadata=dicominfo('10011.dcm');%获取dcm文件的信息
% figure
% imagesc(img_1);
% imshow(uint8(img_1));
Hu0=(int16(img_1)*1+(-1024));
%文档中的第二步,转为CT值的那个。
%% 窗宽窗位设置:c 窗位 w 窗宽 这个是按照你给的公式写的
Hu=double(Hu0);%此时有正负。转双精度的CT值。
gm=255;
c=60;
w=100;
Gv=0.*(Hu<c-w/2)+(gm/w)*(Hu+(w/2)-c).*(c-w/2<=Hu&Hu<=c+w/2)+gm.*(Hu>c+w/2);%Gv为显示灰度。
figure
imshow(uint8(Gv));
title('加窗')
%% 一次CTRL+T
%% 拉普拉斯高斯滤波 (有库函数) %这个是整体区域。
img_gray=uint8(Gv);
hsize=10;%滤波模板大小。自己可以修改
sigma=0.4; %滤波系数sigma 自己修改会得到不同的效果图
h = fspecial('log', hsize, sigma); %构造拉普拉斯_高斯滤波器,'log'是这个滤波器的标志
img_filter0=double(imfilter(Hu0, h)); %调用matlab中的imfilter函数 ,进行滤波。其中img_gray为要过滤的图像,h为滤波器。
figure
img_filter=uint8(img_filter0);
imshow((img_filter));%显示图像
title('拉普拉斯高斯滤波')
%%
%% 自作ROI区域 利用ginput函数标出点,然后包围所标区域。%%对滤波变换后进行抠图
figure
imshow(img_filter);%%此处为img_filter 就是对滤波后抠图,此处为img_gray就是对为滤波的灰度图进行抠图。
title('加窗图抠图')
hold on
x=[];y=[];
n=0;
while(1)
[xtemp,ytemp,button]=ginput(1);
plot(xtemp,ytemp,'r*');
x=[x xtemp];
y=[y ytemp];
n=n+1;
% text(xtemp+0.1,ytemp,int2str(n));
if(button==32) %button=32表示 当你点点完了。你就按空格,退出点点的状态。
break
end
end
line(x,y);%x y连线
% hold off
t=1:n;
tt=1:0.1:n;
xx=spline(t,x,tt); %因为手点 点数不够多,不够精确,需要插值 spline为插值函数
yy=spline(t,y,tt);
plot(xx,yy,'b:');
BW = roipoly(img_filter,xx,yy); %roipoly为抠图函数。以img_filiter为基础。xx yy为轮廓。包围的区域扣出来
figure, imshow(BW) %返回的是一个二值图像BW。就是一张黑白图。扣除的地方为1.其他地方均为0;
title('轮廓图')
%% 下面返回 被扣的地方的图像。 就是图上只有被扣除的东西
figure
img_last=BW.*(double(img_filter));%原理是这样:BW为 1 0的二值图,被扣的地方是兴趣区域都是1.这时
%与原图像进行.* ,除了为1的地方返回灰度值*1,其他地方都返回灰度值*0。
img_last1=uint8(img_last);%转为uint8格式。不然显示不了
img_last1(img_last1==0)=NaN;% 调整背景。
% figure
subplot(221)
imshow(img_last1)%显示最后的图像
title('轮廓返回图')
%%
%你需要的兴趣区域的数据如下: 。它的特点就是 兴趣区域为原图的CT值,非兴趣区域全是0;matlab图像是矩阵。所以必须是规则的。
%所以 要有0 来填充。
final_CT=BW.*(Hu);%%你要是要用 最后的抠图的CT值。这个就是。
%%论文中 有这么一句话:the mean gray-level intensity (m) and uniformity (u)
%%指的是平均灰度强度和均一性。实质求的是ROI区域的灰度值,和灰度值的均一性;
%% l是灰度值, p(l)是l灰度值出现的概率。
%% 先给图像加窗 转为灰度值Gv——>拉普拉斯高斯滤波——>抠图——>求抠图地方的灰度平均值m、均一性值n
%% 那个gui是演示用的。是在没滤波的基础上进行的抠图。没有关系。因为数据在这产生。
%% img_last1 的最后数据。就是 抠出得图。周围用 0填充的。所以周围为黑色。
[p,q]=size(img_last1);
RGB=zeros(p,q);
img_B=cat(3,RGB,RGB,img_last1);%实质就是讲RGB R=0,G=0,B=灰度值。下面类似。
subplot(222)
imshow(img_B);
title('B')
img_G=cat(3,RGB,img_last1,RGB);
subplot(223)
imshow(img_G);
title('G')
img_R=cat(3,img_last1,RGB,RGB);
subplot(224)
imshow(img_R);
title('R')
%% ROI 区域均值求取。 img_mean;
img_sum=sum(img_last1(:));%求出所有元素总和。填充区域为0 。不影响。
N=numel(find(BW==1));%BW中有1的地方 就是有灰度的地方。所以BW中1 的多少,就是灰度的总数。
img_mean=img_sum/N;
大致的结果图如下:
Matlab 高斯_拉普拉斯滤波器处理医学图像的更多相关文章
- MATLAB高斯混合数据的生成
MATLAB高斯混合数据的生成 作者:凯鲁嘎吉 - 博客园 http://www.cnblogs.com/kailugaji/ 高斯混合模型的基本原理:聚类——GMM,MATLAB中GMM聚类算法:M ...
- matlab函数_连通区域
1. matlab函数bwareaopen──删除小面积对象格式:BW2 = bwareaopen(BW,P,conn)作用:删除二值图像BW中面积小于P的对象,默认情况下使用8邻域.算法:(1)De ...
- matlab初学者_脚本文件调用函数文件
问题: matlab里面有两种文件,一种是脚本文件,一种是函数文件,为了模块化程序,我们需要把专门的功能写成一个函数封装到某个函数文件里面. 那么来看如何在脚本文件里调用函数文件中的函数. 注意点: ...
- 【转】matlab函数_连通区域
转载自einyboy的博文Matlab的regionprops详解 1. matlab函数bwareaopen──删除小面积对象格式:BW2 = bwareaopen(BW,P,conn)作用:删除二 ...
- matlab中实现Gabor滤波器
1.spatialgabor.m描述gabor函数 % SPATIALGABOR - applies single oriented gabor filter to an image%% Usage: ...
- matlab高斯模板生成,K是归一化系数,上面是一个半径200的高斯模板用来做MSR
R3=; F=zeros(*R3+,*R3+); sigma=R3/; r=; :*R3+ :*R3+ r=(x-R3)*(x-R3)+(y-R3)*(y-R3); ...
- SSE图像算法优化系列二十八:深度优化局部拉普拉斯金字塔滤波器。
基于局部拉普拉斯金字塔的Edge-aware滤波器是在2011年由Adobe 公司的研究员Sylvain Paris(大神级人物,写了很多文章)提出的,我在4年前曾经参考有关代码实现过这个算法,但是速 ...
- IIR滤波器软件实现(Matlab+C++)
使用C++来写一个IIR滤波器 我们首先要在MATLAB中设计一个IIR滤波器,并生成一个头文件,这个头文件中反映了IIR滤波器的频率响应特性 理论支持 IIR滤波叫做递归滤波器,它是一种具有反馈的滤 ...
- MATLAB用“fitgmdist”函数拟合高斯混合模型(一维数据)
MATLAB用“fitgmdist”函数拟合高斯混合模型(一维数据) 作者:凯鲁嘎吉 - 博客园 http://www.cnblogs.com/kailugaji/ 在MATLAB中“fitgmdis ...
随机推荐
- Tomcat一个BUG造成CLOSE_WAIT
之前应该提过,我们线上架构整体重新架设了,应用层面使用的是Spring Boot,前段日子因为一些第三方的原因,略有些匆忙的提前开始线上的内测了.然后运维发现了个问题,服务器的HTTPS端口有大量的C ...
- Html Agility Pack 解析Html
Hello 好久不见 哈哈,今天给大家分享一个解析Html的类库 Html Agility Pack.这个适用于想获取某网页里面的部分内容.今天就拿我的Csdn的博客列表来举例. 打开页面 用Fir ...
- HTTP协议系列(1)
一.为什么学习Http协议 首先明白我们为什么学习HTTP协议,也就是说明白HTTP协议的作用.HTTP协议是用于客户端与服务器之间的通讯.明白了HTTP协议的作用也就知道了为什么要学习H ...
- Android 判断一个 View 是否可见 getLocalVisibleRect(rect) 与 getGlobalVisibleRect(rect)
Android 判断一个 View 是否可见 getLocalVisibleRect(rect) 与 getGlobalVisibleRect(rect) [TOC] 这两个方法的区别 View.ge ...
- html与html5
HTML 是一种在 Web 上使用的通用标记语言.HTML 允许你格式化文本,添加图片,创建链接.输入表单.框架和表格等等,并可将之存为文本文件,浏览器即可读取和显示.HTML 的关键是标签,其作用是 ...
- 【解决方案】Myeclipse 10 安装 GIT 插件 集成 步骤 图解
工程开发中,往往要使用到集成GIT ,那么下面说说插件安装步骤 PS:以Myeclipse 10 为例,讲解集成安装步骤. ----------------------main------------ ...
- MVC Core 网站开发(Ninesky) 1、创建项目
又要开一个新项目了!说来惭愧,以前的东西每次都没写完,不是不想写完,主要是我每次看到新技术出来我都想尝试一下,看到.Net Core 手又痒了,开始学MVC Core. MVC Core最吸引我的有三 ...
- ASP.NET MVC5+EF6+EasyUI 后台管理系统(67)-MVC与ECharts
系列目录 ECharts 特性介绍 ECharts,一个纯 Javascript 的图表库,可以流畅的运行在 PC 和移动设备上,兼容当前绝大部分浏览器(IE8/9/10/11,Chrome,Fire ...
- Android消息传递之基于RxJava实现一个EventBus - RxBus
前言: 上篇文章学习了Android事件总线管理开源框架EventBus,EventBus的出现大大降低了开发成本以及开发难度,今天我们就利用目前大红大紫的RxJava来实现一下类似EventBus事 ...
- MySQL 系列(一) 生产标准线上环境安装配置案例及棘手问题解决
一.简介 MySQL是最流行的开放源码SQL数据库管理系统,它是由MySQL AB公司开发.发布并支持的.有以下特点: MySQL是一种数据库管理系统. MySQL是一种关联数据库管理系统. MySQ ...