- 博客(16)
- 资源 (9)
- 问答 (1)
- 收藏
- 关注
原创 空间三维散点数据的线性拟合
clc;clear all;close all;num = 50; % num个随机点Rand1 = randi([-1,1],num,3); %噪声范围Rand2 = randi([-1,1],num,3);Point1 = [1:0.5:0.5*(num+1); 1:0.5:0.5*(num+1); 1:0.5:0.5*(num+1)]'+ Rand1;plot3(Point1(
2016-03-31 22:19:40 26717 3
原创 深度理解Powell优化算法
1.鲍威尔优化法综述鲍威尔法又称方向加速法,它由Powell于1964年提出,是利用共轭方向可以加快收敛速度的性质形成的一种搜索方法。该方法不需要对目标函数进行求导,当目标函数的导数不连续的时候也能应用,因此,鲍威尔算法是一种十分有效的直接搜索法。Powell法可用于求解一般无约束优化问题,对于维数n不同于其他的直接法,Powell法有一套完整的理论体系,故其计算效率高于其他直接法。该
2016-03-31 22:04:30 55802 12
原创 互信息的数学解释以及matlab编程
1.互信息互信息是信息论中的一个基本概念,通常用于描述两个系统间的统计相关性,或者是一个系统中所包含另一个系统中信息的多少。引入互信息对图像进行配准是因为:不需要假定两幅图像灰度的对应关系,不需要对图像进行分割。但是这对Ct-MR配准、或者CT-PET配准是可行的,对于CT-US配准是不可行的。在概率论中,两个随机变量A和B,他们的边缘概率分布为p_A (a)和p_
2016-03-30 20:46:50 27879 9
原创 matlab 3D绘图详解
1、三维曲线: plot3( ) →plot3(X1,Y1,Z1,...) →plot3(X1,Y1,Z1,LineSpec,...) →plot3(...,'PropertyName',PropertyValue,...) →h = plot3(...) plot3( )的用法与plot( )类似,只是多了一个 Z 数组。举例:
2016-03-30 11:44:06 40118 1
原创 鼠标取点——ginput()和getpts()详解
1 ginput() 函数ginput提供了一个十字光标使我们能更精确的选择我们所需要的位置,并返回坐标值。函数调用形式为:[x,y] = ginput(n),能使你从当前的坐标系中读取n个点,返回这n个点的x,y坐标,均为nX1的向量。可以按回车提前结束读数。[x,y] = ginput 可以无限的读取坐标直到按下回车键。[x,y,button] = ginput
2016-03-30 10:35:42 48013 1
原创 特征点提取之Harris角点提取法
1. 特征点提取的意义2.角点3. Harris角点检测的基本原理4.Harris角点检测算法的步骤5.Harris角点提取算法设计function [ptx,pty] = HarrisPoints(ImgIn,threshold)% Harri
2016-03-29 21:39:53 12808 2
原创 深度理解高斯滤波器
1.高斯滤波器综述 高斯滤波器是一类根据高斯函数的形状来选择权值的线性平滑滤波器。高斯平滑滤波器对于抑制服从正态分布的噪声非常有效。一维零均值高斯函数为: g(x)=exp( -x^2/(2 sigma^2)其中,高斯分布参数Sigma决定了高斯函数的宽度。对于图像处理来说,常用二维零均值离散高斯函数作平滑滤波器。 高斯函数具有五个重要的性质,
2016-03-28 16:21:16 16967
原创 uint8、double、mat2uint、imhist绘直方图w恩替藐视
1.问题描述ImgUS = imread('050.bmp');Igray = rgb2gray(ImgUS);Igray = double(Igray);figure,imshow(Igray ,[]);figure;imhist(Igray);图像可以正常显示,但是灰度直方图却无法显示。2.解答:2.1 网友的说法: 默认情况下,matlab将图像中的数据存
2016-03-25 17:41:53 2630
原创 矩阵的特征值和特征向量
1.目的:1.1理解矩阵的特征值和矩阵向量的概念及性质1.2会求矩阵的特征值以及特征向量2. 特征矩阵、特征向量的定义及求解3.经典例题4. 讨论由以上讨论可知,对于方阵的每一个特征值,我们都可以求出其全部的特征向量.但对于属于不同特征值的特征向量,它们之间存在什么关系呢?这一问题的讨论在对角化理论中有很重要的作用.对此有结论如下:
2016-03-22 11:45:53 6847
原创 如何获得文件路径名、带后缀的文件名、不带后缀的文件名。
global pathname imname;%fileneme 指的就是带文件后缀的文件名,如‘050.bmp’%pathname指的是文件路径名[filename, pathname] = uigetfile('*.bmp','Input filename');%寻找后缀名前面的标志‘.’i = find('.'==filename);%去除文件后缀,提取单纯的文件名
2016-03-19 22:08:20 12489
原创 基于光学导航系统,矩阵变换和3D-2D配准研究
1. 手术导航任务中,3D-2D配准的数学框架2. 框架分析三维实时影响导航的理解,完全考虑刚体的情况下:1. 利用标准体模对导航系统进行校准。此处的校准是指得到上文中提到的校准矩阵,即把2维超声数据从图像坐标系变换到传感器坐标系下。2. 术前采集得到CT数据,经三维重建可以得到3DCT图像坐标系下的三维立体结构。3. 术中实时采集2D超声图像,利用
2016-03-19 21:08:46 3477 2
转载 我可以蹲下来,陪你做一只蘑菇
有一个精神病人,以为自己是一只蘑菇。于是他每天都撑着一把伞蹲在房间的墙角里,不吃也不喝,像一只真正的蘑菇一样。 心理医生想了一个办法。有一天,心理医生也撑了一把伞,蹲坐在了病人的旁边。病人很奇怪的问:你是谁呀?医生回答:我也是一只蘑菇呀。病人点点头,继续做他的蘑菇。过了一会儿,医生站了起来,在房间里走来走去。病人就问他:你不是蘑菇么,怎么可以走来走去?医生回答说:蘑菇当然可以走来走去啦!病人觉
2016-03-19 10:32:40 1185 1
原创 自适应中值滤波用于超声图像降噪
自适应中值滤波原理:RAMF主要通过以下两步来处理图像。1.首先确定最大的滤波半径,然后用一个合适的半径r对图像进行滤波。计算当前滤波半径像素灰度的Imin,Imax,Imed,然后判断Imed是否在[Imin,Imax]中间,如果在则向下进行,否则扩大当前半径r继续滤波直到r等于最大滤波半径。2.如果当前处理的像素img(i,j)在[Imin,Imax]之间,则输出当前像素,否则输出
2016-03-16 11:46:09 4801
原创 区域生长算法原理及MATLAB实现
1. 基于区域生长算法的图像分割原理数字图像分割算法一般是基于灰度值的两个基本特性之一:不连续性和相似性。前一种性质的应用途径是基于图像灰度的不连续变化分割图像,比如图像的边缘。第二种性质的主要应用途径是依据实现指定的准则将图像分割为相似的区域。区域生长算法就是基于图像的第二种性质,即图像灰度值的相似性。1.1 基本公式令R表示整幅图像区域,那么分割可以看成将区域R划分为n个子区域R1
2016-03-13 21:20:57 104831 50
转载 双边滤波Matlab实现<The Bilateral Filter>
原文:http://blog.csdn.net/lifeitengup/article/details/8902326#comments双边滤波与一般的高斯滤波的不同就是:双边滤波既利用了位置信息又利用了像素信息来定义滤波窗口的权重。像素值越接近,权重越大。双边滤波会去除图像的细节信息,又能保持边界。对于彩色图像,像素值的接近与否不能使用RG
2016-03-10 17:29:01 9312
转载 双边滤波算法原理
一、引言 双边滤波在图像处理领域中有着广泛的应用,比如去噪、去马赛克、光流估计等等,最近,比较流行的Non-Local算法也可以看成是双边滤波的一种扩展。自从Tomasi et al等人提出该算法那一天起,如何快速的实现他,一直是人们讨论和研究的焦点之一,在2011年及2012年Kunal N. Chaudhury等人发表的相关论文中,提出了基于三角函数关系的值域核算法,能
2016-03-09 21:51:18 30689 1
Python+OpenCV实现Selective Search算法
2018-10-29
英伟达深度神经网络加速库 与 matconvnet配套使用 (v5.1)
2018-07-11
21个项目玩转深度学习
2018-06-15
TA创建的收藏夹 TA关注的收藏夹
TA关注的人