基于最小误差阈值法的特定医学图像分割算法研究-含Matlab代码

news2024/10/7 8:22:21

⭕⭕ 目 录 ⭕⭕

  • ✳️ 一、引言
  • ✳️ 二、阈值分割原理
  • ✳️ 三、最小误差阈值法
  • ✳️ 四、基于最小误差阈值法图像分割的实验分析
    • ✳️ 4.1 图像分割结果
    • ✳️ 4.2 边缘曲线提取结果
  • ✳️ 五、参考文献
  • ✳️ 六、Matlab代码获取

✳️ 一、引言

图像是人们所熟悉的,但是图像里的内容并不全是人们所需要的,因此人们便想方设法从图像里提取自己所需要的部分,有学者在 20 世纪 60 年代就已经开始研究图像分割这一领域了,而且图像分割也是图像处理到分析的一个重要环节。因此,各界的研究者都对图像分割的研究很重视。

图像分割在数字图像处理领域中起着很大的作用,因此分割的最终结果决定着其他图像处理工作能否顺利完成。从上述图像分割的历史渊源可知,图像分割的研究意义明显,所以有大量学者从事这一领域的研究工作。而且图像分割应用广泛,各行各业均有相对应的应用范围这让它成为一门热门学科。所谓图像的分割,其具体是指在进行相关图像的处理时,将其中感兴趣、具有相对大价值对象加以分离,此项技术手段是针对图像进行处理时所应用到一项极为关键技术。在医学领域,图像分割技术对于医学图像的研究有着无法忽视的作用。在癌症研究领域,癌症的诊断目的是发现肿瘤。癌症研究机构的一份报告显示,在 70%至 80%的癌症死亡病例是由于早期肿瘤无法检测到造成的。

在学者们提出了图像分割的阈值处理技术之后,该项技术得以不断改进,被应用在各行各业。在阈值分割原理的基础上提出了最小误差法阈值法。该方法在实际应用中操作简单且分割效果较理想深受研究者的喜爱。

✳️ 二、阈值分割原理

阈值分割方法是图像分割领域中较为常见的算法。其主要原理:借助于图像所拥有的相应灰度分析得出不同区域对应的阈值,然后把所得阈值和分析图像对应像素具体灰度情况加以对比,最终将像素根据分割结果分到合适的分类中完成图像分割。在进行阈值的选取过程中,其中最为重要的是依据相应的准则函数对最佳阈值进行计算,只有确保了所设定阈值更为准确才可以将图像正确分割。通过原始图像f(𝑥, 𝑦)选取阈值 T 进行分割,分割后的图像为g(𝑥, 𝑦),可以定义为如下公式所示:

在这里插入图片描述

在图像f(𝑥, 𝑦)进行了分割处理之后,产生对应的二值图像,在这之中部分黑色区域属于分割的目标范围,而白色区域属于分割的背景范围,从而可以有效的将分割的目标和背景加以区别。通过上述分析能够了解到,采用阈值针对图像进行分割处理其原理为把图像对应的灰度数据加以分类处理,同时依照阈值选取的具体准则,找出能够确保目标以及背景更为精准与快速划分出来所对应的最佳阈值。

✳️ 三、最小误差阈值法

采用最小误差阈值技术过程中,要求应当设定所要分割的具体图像之中包含有目标对象以及背景,同时要符合高斯分布,在针对所要分割的目标以及图像中背景部分对应均值以及方差值加以计算,并根据最小误差思想得到最小误差目标函数,最后计算目标函数最小时的取值并将其作为最佳阈值,根据阈值将图像分割为二值图像。该算法能够很好地适应医学图像的灰度分布特点,取得有效的目标区域分割效果。

通过利用函数f(𝑥, 𝑦)代表尺寸值是 M × N 图像之中对应(𝑥, 𝑦)坐标位置像素点所拥有灰度值大小, 在这里插入图片描述 。通过利用直方图形h(g)进一步的对不同灰度所对应分布情况进行计算与统计。设定所构建的灰度模型其实际分布情况和混合正态分布相符和情况下,则 在这里插入图片描述,其中 Pi 代表的为子分布所对应先验概率值大小, p(g)所对应的不同子分布 在这里插入图片描述服从于均值大小是ui ,方差值大小是σi 正态分布,即 在这里插入图片描述

而在满足阈值t∈G基础上,学者 Kittler 以及学者 Illingworth 结合最小误差相关理论和计算方法,得出了如公式所示:

在这里插入图片描述

其中

在这里插入图片描述
最佳阈值选为使 J(t)最小的 t= t∗

在这里插入图片描述
上述方法称为一维最小误差阈值法。

✳️ 四、基于最小误差阈值法图像分割的实验分析

✳️ 4.1 图像分割结果

由于一维最小误差阈值法只依赖于一维直方图,所以容易受到图像中噪声的干扰。一维最小误差阈值法不关注像素平均灰度值,只关注像素的灰度值,这样对图像分割的性能和效果会有极大的影响。利用一维最小误差阈值法对胸片图像进行肺部分割,分割结果如图1所示。

在这里插入图片描述

图1 分割结果

✳️ 4.2 边缘曲线提取结果

经过最小误差法分割得到的二值图像包含一定的孔洞、噪声干扰,可通过形态学后处理进行孔洞填充、去噪操作,去噪效果如下图2(a)所示。形态学后处理过程得到了较好的分割效果,能对应到目标肺部区域。为了进一步定位肺部边缘,可继续执行边缘曲线提取操作,在原图中定位肺部的边缘曲线,结果图2(b)所示。

在这里插入图片描述

图2 边缘提取

✳️ 五、参考文献

[1] Kittler J, Illing; worth J. Minimum error thresholding. Pattern Recognition, 1986,19(1): 41-47.

✳️ 六、Matlab代码获取

上述Matlab代码,可私信博主获取。


博主简介:研究方向涉及智能图像处理、深度学习等领域,先后发表过多篇SCI论文,在科研方面经验丰富。任何与算法程序科研方面的问题,均可私信交流讨论


本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.coloradmin.cn/o/28008.html

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈,一经查实,立即删除!

相关文章

服务器稳定性测试-LTP压力测试方法及工具下载

简介 LTP(LinuxTest Project)是SGI、IBM、OSDL和Bull合作的项目,目的是为开源社区提供一个测试套件,用来验证Linux系统可靠性、健壮性和稳定性。LTP测试套件是测试Linux内核和内核相关特性的工具的集合。 该工具的目的是通过把测…

【HMS Core】构建SplitBill应用集成多个HMS Core服务,助力您更好的了解华为生态组成

一、介绍 Duration: 3:00 总览 通过构建本次的SplitBill应用,您可以更好地了解华为生态的组成部分,包括认证服务、云存储和云数据库等Serverless服务。此外您还可以了解如何使用近距离数据通信服务的Nearby Connection功能分享文件。无需使用现金&…

重回synchronized(源码解读与实战解析篇)

读前必知 文中锁,也称为对象锁,而锁对象就是指的承载这个锁的对象,如下面,用法中所指的Object o,在print3中就是锁对象。 以下源码分析仅适用于jdk8,版本未知,因为源码提供者的源码版本访问地址…

金仓数据库KingbaseES GIN 索引

目录 一、索引的逻辑结构 二、索引的物理结构 三、GIN索引使用例子 1、前后模糊查询 2、全文检索 四、gin 索引可用于超长的字段 GIN(Generalized Inverted Index, 通用倒排索引) 是一个存储对(key, posting list)集合的索引结构,其中key是一个键值&#xff0c…

【uniapp】利用Vuex实现购物车功能

实战项目名称:实现购物车功能 文章目录一、实战步骤1. 先编辑store.js文件2. 定义方法和基本的结构3. 编写SETSHPPING二、在项目中调用1. 触发相应的mutations2. 利用computed计算数量和总价的方法提示:本实战内容大部分为具体实现的思路,界面…

C语言 2 —— 常量

常量是什么? 常量就是在程序运行过程中,值不会发生改变,而且一眼可以辨识出值的量。 如: 20,‘a’ , 3.1415926 , "helloworld" 常量的分类: 整形,浮点型,字符型&#…

成电860考研专业课考前划重点-学长课程音频转文字-用科大讯飞花钱买的-三万五千字

成电860考研专业课考前划重点-学长课程音频转文字。 这个是我2021年把视频的音频扒下来后用科大讯飞音频转文字网站上花了几十块钱买的。 说话人1 03:04 对有回放,我这边开回放了,大家可以让大家下载,然后我怕这边回放有问题,大家…

Airtest poco 入门小结

目录 一、poco介绍 1、poco能做什么 2、三种定位方式 1)基本选择器 2)相对选择器 3)空间顺序选择器 4)正则表达式方式 3、poco支持平台 二、Airtest介绍 1、Airtest能做什么 2、Airtest不能做什么 3、Airtest的图像识别…

发挥数字化平台优势,电子元器件采购商城系统助力企业改变固有低效流程

我国是全球最大的电子元器件生产国,电子元器件在国民经济发展中占据着重要地位。近年来,随着数字经济的快速发展,电子元器件的需求量也在不断升高,但疫情与国际环境对电子元器件产业要素流通仍然造成了一定困扰。在此背景下&#…

win10系统下使用onnxruntime部署yolov5模型

文章目录前言一、环境1、硬件2、软件二、YOLO模型三、新建Qt项目1、pro文件2、mainwindow.h3、mainwindow.cpp四、YOLO 类封装1、yolov5.h2、yolov5.cpp3、class.names五、效果前言 上一篇介绍过使用opencv-dnn模块实现模型推理部署,但视频效果较差,本篇…

ps2023最新版免费滤镜插件Exposure安装下载教程

滤镜插件是ps的重要功能之一,它主要是用来制作不同的图片特效。那么,ps滤镜插件哪些好用,ps滤镜插件如何获取,下面我们一起来学习这些内容。 ps滤镜插件是比较多的,下面对几款常见的ps滤镜插件进行讲解,看…

TIA博途中的TRACE功能具体使用方法示例

TIA博途中的TRACE功能具体使用方法示例 我们可以利用TRACE曲线来监控程序、排查故障,那么具体怎样使用呢,可以参考以下内容。 如下图所示,打开TIA博途,新建项目后,在左侧项目树中可以看到TRACES, 如下图所示,双击添加新轨迹,然后在右侧窗口中,添加需要监视的信号,…

TOUGH2系列建模方法及在CO2地质封存、水文地球化学、地热、地下水污染等领域中的技术

TOUGH2系列软件是由美国劳伦斯伯克利实验室开发的,旨在解决非饱和带中地下水、热运移的通用模拟软件。和传统地下水模拟软件Feflow和Modflow不同,TOUGH2系列软件采用模块化设计和有限积分差网格剖分方法,通过配合不同EOS模块,软件…

【矩阵论】3. 矩阵函数——矩阵函数求导

3.6 矩阵函数求导 3.6.1 积分与求导定义 设 mnm\times nmn 阶矩阵 A(x)(aij(x))mnA(x)\left(a_{ij}(x)\right)_{m\times n}A(x)(aij​(x))mn​ 中的元素都是 x 的可导函数,则 A(x)A(x)A(x) 为关于 xxx 的求导为: A′(A)dA(x)dx(daij(x)dx)mnA(A)\frac{…

正点原子stm32F407学习笔记5——串口通信实验

一、串口通信实验1 上位机给开发板发送数据,开发板将收到的数据发回给上位机 串口设置的一般步骤可以总结为如下几个步骤: 串口时钟使能,GPIO 时钟使能。设置引脚复用器映射:调用 GPIO_PinAFConfig 函数。GPIO 初始化设置&#…

数据库性能翻3倍:Redis on Flash分层存储技术是如何做到的?

Redis on flash简介:Redis on Flash 涉及到的是Redis的分层存储技术,即将数据存放在不同地方。Redis自2016年以来支持Redis on Flash。从2019年开始, Redis企业版(Redis Enterprise)宣布支持英特尔Optane DC持久性内存&#xff0c…

基于NB-IoT的智能垃圾桶系统设计与实现

本设计是基于物联网的智能垃圾桶,主要实现以下功能: 1,压力传感器模块采集垃圾重量数据; 2,GPS定位模块采集垃圾桶所在的经纬度数据; 3,人体红外模块检测人体并返回是否有人通过的数据&#xf…

会议管理系统SSM记录(二)

目录: (1)整合Freemarker (2)用户登录 (3)提取头部 (4)提取菜单抽取 (1)整合Freemarker 在pom.xml中加入Freemark依赖: 创建free…

HTML+CSS大作业:基于HMTL校园学校网页设计题材【我的学校网站】

🎉精彩专栏推荐 💭文末获取联系 ✍️ 作者简介: 一个热爱把逻辑思维转变为代码的技术博主 💂 作者主页: 【主页——🚀获取更多优质源码】 🎓 web前端期末大作业: 【📚毕设项目精品实战案例 (10…

专精特新企业三个层级

专精特新企业也是分层级的。工信部2022年6月印发《优质中小企业梯度培育管理暂行办法》里面明确提出中小企业培育的3个梯度,分别是创新型中小企业、专精特新中小企业和专精特新小巨人企业,刚好构成中小企业发展层级金字塔。这就意味着企业想要发展崛起&a…