在医院实际环境中,经常遇到有问题的患者,对于一些特殊的场景,比如骨折,肺结节,心脑血管问题 需要图像对比增强来更为清晰的显示病灶助于医生确诊,先看效果: 肺纹理增强: 肺结节增强: 血管对比增强: 骨骼对比增强: 根据参考资料: MATLAB版本: https://ww2.mathworks.cn/matlabcentral/fileexchange/24409-hessian-based-frangi-vesselness-filter 算法原理: https://baike.baidu.co…
根据本系列教程文章上一篇说到,在完成C++和Opencv对Hessian矩阵滤波算法的实现和封装后, 再由C#调用C++ 的DLL,(参考:C#处理医学图像(一):基于Hessian矩阵的血管肺纹理骨骼增强对比) 功能虽然已经实现,但在实际应用中要考虑到性能和耦合,本篇将介绍性能方面的注意点以及和其他功能的联动. 我们将Demo里面的功能集成到正式工程中: 1.新建一个新窗体,用来显示结果和调整滤波参数: 其中滑块控件选择工具箱中的Slider,定义好控件样式,变化事件选择PreviewMous…
matlab中repmat函数的用法 B = repmat(A,m,n) B = repmat(A,[m n]) B = repmat(A,[m n p...]) 这是一个处理大矩阵且内容有重复时使用,其功能是以A的内容堆叠在(MxN)的矩阵B中,B矩阵的大小由MxN及A矩阵的内容决定,如果A是一个3x4x5的矩阵,有B = repmat(A,2,3)则最后的矩阵是6x12x5 例如: >>B=repmat( [1 2;3 4],2,3) B = 1      2      1     2  …
所用的函数非常简单,只需要用到mesh函数,示例代码如下: Ima=imread('F:\pathto\test.jpg'); surf_ima = surf(rgb2gray(Ima)); %黑色的3D图 title('3D') mesh_ima = mesh(rgb2gray(Ima)); %有色彩的3D图 参考链接:https://cn.mathworks.com/matlabcentral/answers/17998-image-processing-how-can-i-create-a…
Hessian矩阵与多元函数极值 海塞矩阵(Hessian Matrix),又译作海森矩阵,是一个多元函数的二阶偏导数构成的方阵.虽然它是一个具有悠久历史的数学成果.可是在机器学习和图像处理(比如SIFT和SURF特征检測)中,我们也经常遇到它.所以本文就来向读者道一道Hessian Matrix的来龙去脉.本文的主要内容包括: 多元函数极值问题 泰勒展开式与Hessian矩阵 多元函数极值问题 回忆一下我们是怎样处理一元函数求极值问题的. 比如.f(x)=x2,我们会先求一阶导数,即f′(x)…
1. 前言 什么是特殊矩阵? C++,一般使用二维数组存储矩阵数据. 在实际存储时,会发现矩阵中有许多值相同的数据或有许多零数据,且分布呈现出一定的规律,称这类型的矩阵为特殊矩阵. 为了节省存储空间,可以设计算法,对这类特殊矩阵进行压缩存储,让多个相同的非零数据只分配一个存储空间:对零数据不分配空间. 本文将讲解如何压缩这类特殊矩阵,以及压缩后如何保证矩阵的常规操作不受影响. 2. 压缩对称矩阵 什么是对称矩阵? 在一个n阶矩阵A中,若所有数据满足如下述特性,则可称A为对称矩阵. a[i][j]…
最近在做关于数据采集方面的东西,这就不免涉及到了滤波的算法,在网上找到了关于几种算法的比较. 数字滤波方法有很多种,每种方法有其不同的特点和使用范围.从大的范围可分为3类. 1.克服大脉冲干扰的数字滤波法 ㈠.限幅滤波法 ㈡.中值滤波法 2.抑制小幅度高频噪声的平均滤波法 ㈠.算数平均 ㈡.滑动平均 ㈢.加权滑动平均 ㈣一阶滞后滤波法 3.复合滤波法 四.介绍 在这我选用了常用的8种滤波方法予以介绍 (一)克服大脉冲干扰的数字滤波法: 克服由仪器外部环境偶然因素引起的突变性扰动或仪器内部不稳定引…
原文地址为:matlab矩阵的表示和简单操作 一.矩阵的表示在MATLAB中创建矩阵有以下规则: a.矩阵元素必须在”[ ]”内: b.矩阵的同行元素之间用空格(或”,”)隔开: c.矩阵的行与行之间用”;”(或回车符)隔开: d.矩阵的元素可以是数值.变量.表达式或函数: e.矩阵的尺寸不必预先定义. 二,矩阵的创建: 1.直接输入法 最简单的建立矩阵的方法是从键盘直接输入矩阵的元素,输入的方法按照上面的规则.建立向量的时候可以利用冒号表达式,冒号表达式可以产生一个行向量,一般格式是: e1:…
目前比较经典的图像去噪算法主要有以下三种: 均值滤波:也称线性滤波,主要思想为邻域平均法,即用几个像素灰度  的平均值来代替每个像素的灰度.有效抑制加性噪声,但容易引起图像模糊,  可以对其进行改进,主要避开对景物边缘的平滑处理. %x是需要滤波的图像,n是模板大小(即n×n) function d=avg_filter(x,n)    a(1:n,1:n)=1;   %a即n×n模板,元素全是1 [height, width]=size(x);   %输入图像是hightxwidth的,且hi…
 基于Vivado HLS在zedboard中的Sobel滤波算法实现 平台:zedboard  + Webcam 工具:g++4.6  + VIVADO HLS  + XILINX EDK + XILINX SDK 系统:ubuntu12.04 总体设计思路 sobel 算法理论基础       索贝尔算子(Sobel operator)主要用作边缘检测,在技术上,它是一离散性差分算子,用来运算图像亮度函数的灰度之近似值.在图像的任何一点使用此算子,将会产生对应的灰度矢量或是其法矢量. 该…