精华内容
下载资源
问答
  • Matlab神经网络自适应算法的Matla盲源分离-神经网络自适应算法的Matla盲源分离.rar 神经网络自适应算法的Matla盲源分离
  • Matlab超强经典教程2014

    2019-03-07 23:37:45
    Matlab超强经典教程2014
  • Matla用于导航的simple滤波和PCRB!Matla用于导航的simple滤波和PCRB
  • 图像灰度三维matla显示,图像长宽和灰度值组成三维图像
  • 数控伺服系统的ADAMS MATLA联合仿真研究zip,数控伺服系统的ADAMS MATLA联合仿真研究
  • Matla下data结构变量与XML互换 specor_m1.m Mike Mayer's
  • 人脸表情识别 matla程序代码 非常有用
  • matlab超声波源码
  • MATLAB超前校正在自动控制原理课程中的应用
  • matla绘图教程和实例

    2011-07-05 15:52:46
    matla绘图教程和实例,很经典使用,包括二维和三维图......
  • MATLAB超前校正在自动控制原理课程中的应用.pdf
  • MATLAB深度学习教程,Matlab超强经典教程,matlab源码
  • OFDM MATLA

    2014-03-08 17:59:37
    基于MATLAB的OFDM系统仿真,完成对OFDM系统性能的分析
  • BCD_GroupLasso_grouplasso_BCD算法matlab_matlabgrouplasso_bcd_matla
  • 基于matla的FHN相平面仿真设计(FHN神经元)
  • matlab超级应用指南

    2010-02-23 18:30:27
    matlab超级应用指南 十分的实用,浅显易懂,用以理解
  • matlab超强实用教程

    2011-09-16 19:55:28
    matlab超强实用教程,新手必备哦
  • MATLA 2009B 中安装Sedumi

    2021-02-09 14:49:39
    MATLA 2009B 中安装Sedumi为了安装Sedumi安装过程遇到的问题及解决 为了安装Sedumi 历时两天,终于把Sedumi装好了,主要是安装C/C++编译器试了各种办法。我这里主要讲成功安装的方法,试错的过程就省略了,如果遇到...

    为了安装Sedumi

    历时两天,终于把Sedumi装好了,主要是安装C/C++编译器试了各种办法。我这里主要讲成功安装的方法,试错的过程就省略了,如果遇到同样问题的道友可参考。

    安装过程遇到的问题及解决

    1 安装C/ C++编译器
    MATLAB没有为MATLAB 2015以前的版本提供编译器,打开C:\Program Files\MATLAB\R2009b\bin\win64\mexopts可以看到MATLAB 2019B支持的编译器版本,msvc90就是VS2008,刚好我电脑里有VS2008的安装包就选择了它。
    MATLAB 2009B支持的编译器版本

    1. 安装VS2008
      MATLAB 命令窗口运行mex -setup,"Would you like mex to locate installed compilers [y]/n? "选n也可以看到MATLAB 2019B支持的编译器版本。
    >> mex -setup
    Please choose your compiler for building external interface (MEX) files: 
     
    Would you like mex to locate installed compilers [y]/n? n
     
    Select a compiler: 
    [1] Intel C++ 9.1 (with Microsoft Visual C++ 2005 SP1 linker) 
    [2] Intel Visual Fortran 10.1 (with Microsoft Visual C++ 2005 SP1 linker) 
    [3] Microsoft Visual C++ 2005 SP1 
    [4] Microsoft Visual C++ 2008 Express 
    [5] Microsoft Visual C++ 2008 SP1 
    Compiler: 
    

    可以看到[5]就是 Microsoft Visual C++ 2008 SP1 。因此,安装完VS2008后我又安装了VS2008 SP1补丁(https://www.onlinedown.net/soft/1153213.htm)。

    1. 选择编译器
      接着我们选择了[5]号编译器,即Compiler: 5,回车后:
    The default location for Microsoft Visual C++ 2008 SP1 compilers is C:\Program Files (x86)\Microsoft Visual Studio 9.0, 
    but that directory does not exist on this machine.  
     
    Use C:\Program Files (x86)\Microsoft Visual Studio 9.0 anyway [y]/n? 
    The default location for Microsoft Visual C++ 2008 SP1 compilers is C:\Program Files (x86)\Microsoft Visual Studio 9.0, 
    but that directory does not exist on this machine.  
    

    第一句说VS2008的默认安装路径是“C:\Program Files (x86)\Microsoft Visual Studio 9.0”,但是这个路径在电脑中不存在。
    这个地方折腾了好几遍才搞定,虽然我把VS2008安装在默认路径“C:\Program Files (x86)\Microsoft Visual Studio 9.0”,但是MATLAB找不到它……

    Use C:\Program Files (x86)\Microsoft Visual Studio 9.0 anyway [y]/n? n
    Please enter the location of your compiler: [C:\Program Files (x86)\Microsoft Visual Studio 9.0] 
    

    解决办法是第二句仍然使用默认路径选择n,回车后把VS2008的安装路径粘贴过去(其实跟默认路径是一样的呀……),回车

     Please verify your choices: 
     
    Compiler: Microsoft Visual C++ 2008 SP1  
    Location: C:\Program Files (x86)\Microsoft Visual Studio 9.0 
     
    Are these correct [y]/n? y
    Use C:\Program Files (x86)\Microsoft Visual Studio 9.0 anyway [y]/n? 
    

    ,然后Are these correct [y]/n?回答y确认正确就可以了

    *************************************************************************** 
      Warning: MEX-files generated using Microsoft Visual C++ 2008 require 
               that Microsoft Visual Studio 2008 run-time libraries be  
               available on the computer they are run on. 
               If you plan to redistribute your MEX-files to other MATLAB 
               users, be sure that they have the run-time libraries. 
    *************************************************************************** 
     
    Trying to update options file: C:\Users\shuji\AppData\Roaming\MathWorks\MATLAB\R2009b\mexopts.bat 
    From template:              C:\PROGRA~1\MATLAB\R2009b\bin\win64\mexopts\msvc90opts.bat 
     
    Done . . . 
    

    3.安装Sedumi
    安装Sedumi的具体过程就不多说了,网上的资料很多。主要讲一下安装过程中遇到的3个问题。
    (1)运行install_sedumi,报错:

    Error: fatal error LNK1113: invalid machine type
    搜到此答案https://zhidao.baidu.com/question/337992282.html,觉得很有道理,运行VS2008的安装文件,在安装界面选择“更改或移除Visual Studio 2008”,“添加或删除功能”,果然X64编译器和工具没有勾选安装,勾选后“更新”,成功!

    在这里插入图片描述
    在这里插入图片描述
    (2)报错2:
    LINK : fatal error LNK1104: cannot open file ‘libsvmread.mexw64’

    解决办法:以管理员身份运行MATLAB

    (3)报错3:
    C:\PROGRA~1\MATLAB\R2009B\BIN\MEX.PL: Error: ‘ddot.c’ not found.
    WHY? WHY? WHY?
    解决办法:在Sedumi路径下运行instal_sedumi.m

    总算安装成功!!菜鸟用个软件好难……

    展开全文
  • 语音端点检测算法的研究与matla程序仿真与实现设计说明.doc
  • matlab超声成像代码骨重建 从超声图像重建 3D 骨骼结构 主要的 MATLAB 代码是RunDemo.m 。 结果位于/Result文件夹中。 使用 ImageVis3D 查看.uvf文件。
  • matlab超声成像代码超声跟踪 概述:该项目使用 MATLAB 的计算机视觉工具箱在动态超声成像过程中跟踪肌肉结构。 注意:我会尝试评论此代码并发布任何改进。 如果您有兴趣使用此代码,请随时使用和参考我们的工作。 ...
  • matlab超声波源码
  • matla简介

    2014-10-31 21:47:11
    matlab简介 介绍简单方便 有例子和解释
  • 基于MATLA的图像复原系统二、图像复原的基本概念图像在形成、传输和记录过程中,由于受到多种原因的影响,图像的质量就会有所下降,典型的表现为图像模糊、失真、有噪声等,这一过程称为图像的退化[5]。图像复原是...

    基于MATLA的图像复原系统



    二、图像复原的基本概念

    图像在形成、传输和记录过程中,由于受到多种原因的影响,图像的质量就会有所下降,典型的表现为图像模糊、失真、有噪声等,这一过程称为图像的退化[5]。

    图像复原是试图利用退化过程的先验知识使已退化的图像恢复本来面目,即根据退化的原因,分析引起退化的环境因素,建立相应的数学模型,并沿着使图像降质的逆过程恢复图像[5]。目的在于消除或减轻在图像获取以及传输的过程中造成的图像品质下降,恢复图像的本来面目。因此,复原技术就是把退化模型化,并采用相反的过程进行处理,以便尽可能复原被退化图像的本来面目。

    广义上讲,图像复原是一个求逆问题,逆问题经常存在非唯一解,甚至无解。要想恢复全真的景物图像比较困难。为了得到逆问题的有用解,图像复原本身往往需要一个质量标准,即衡量接近全真景物图像的程度,或者说,对图像的估计是否达到最佳的程度。需要有先验知识以及对解的附加约束条件。

    典型的图像复原是根据图像退化的先验知识建立一个退化模型,以此模型为基础,采用各种逆退化处理的方法进行恢复,使图像质量得到改善。


    2.2图像退化的数学模型

    一般来说,图像的生成可以简单地被描述为如下数学模型:


    v2-9d69335795ee2cc61240505a8ba24eef_b.jpg


    f(x,y)是成像景物,H是综合退化因子,g(x,y)是退化图像。

    图像f(x,y)可以表示为:


    v2-a6252b6f8d140cfa2d1ec7802816e060_b.jpg


    用卷积符号表示为:

    v2-b6a06f0542d13f162ccb0e84fa6f7a95_b.jpg


    因此还有:

    v2-5892e1fef7e44d1610da1d661adbf5e5_b.jpg


    式中,

    v2-646615125453d9dbc5f7d30dea3c3c78_b.jpg

    是像素点的特性函数,

    v2-fed18522d151443137c40941d3088628_b.jpg

    为冲击响应。

    假定成像系统是线性移不变系统:

    退化模型如图所示

    v2-3464051f776533ae713a65c3cac59404_b.jpg


    +



    v2-6d5193941103531144a829b2803e5c8f_b.jpg

    v2-87073d11cd34de6aad926056a613d7d4_b.jpg


    不考虑加性噪声:

    v2-d50d0e1322624f36e7b9d805165c8674_b.jpg


    考虑加性噪声:

    v2-72855ea6528a61a8b7983c68f9e096b7_b.jpg


    卷积等同于频域内乘积:

    v2-3c6ed6f34116f0d75bcbc810a90f8d56_b.jpg


    2.3逆滤波复原

    逆滤波复原法也叫做反向滤波法,其主要过程是首先将要处理的数字图像从空间域转换到傅里叶频域中,进行反向滤波后再由频率域转回到空间域,从而得到复原的图像信号[5]。

    1.在不考虑噪声的情况下:

    v2-51086af7bb4478ee2424197fd5dc8637_b.jpg


    上式两边进行傅里叶变换得

    v2-ba177111947e99851d3548c4a35df6ba_b.jpg


    则原始图像

    v2-ae8c89f99852384dc6b04628a4c3555e_b.jpg

    =

    v2-ce89e92055bbce3fa3b13c977a3ac529_b.jpg

    然后进行傅里叶逆变换,就可以得到原始图像。由此可看出,如果已知退化图像的傅里叶变换和“滤波”传递函数,则可以求得原始图像的傅里叶变换,经反傅里叶变换就可以求得原始图像f(x,y),这就是逆滤波法的基本原理。但在实际中用逆滤波法存在病态的情况:当H(u,v)=0时,或非常小的数值点上,F(u,v)将变成无穷大或非常大的数。

    2.在有噪声的情况下:

    逆滤波原理可以写成:G(u,v)=F(u,v)H(u,v)+N(u,v)

    写成逆滤波的方式:

    (u,v)=F(u,v)+

    v2-d3004cc508f6e1bda19bb0a02d009fa6_b.jpg


    但实际用逆滤波存在病态的情况:

    噪声存在,当H(u,v)很小或为零时,则噪声被放大。这意味着退化图像中小噪声的干扰在H(u,v)较小时,会对逆滤波恢复的图像产生很大的影响,有可能使恢复的图像和f(x,y)相差很大,甚至面目全非。

    实验证明,当退化图像的噪声较小,即轻度降质时,采用逆滤波复原的方法可以获得较好的结果。通常,在离频率平面原点较远的地方数值较小或为零,因此图像复原在原点周围的有限区域内进行,即将退化图像的傅里叶频谱限制在没出零点而且数值又不是太小的有限范围内。

    2.4 维纳滤波复原

    逆滤波比较简单,但没有清楚地说明如何处理噪声,而维纳滤波综合了退化函数和噪声统计特性两个方面进行复原处理。维纳滤波是维纳在1949年提出的,并应用于一维平稳时间序列,获得了满意的结果。这是最早也是最著名的线性滤波技术。

    采用维纳滤波是假设图像信号可以近似看成平稳随机过程的前提下,按照使f(x,y)和之间的均方误差达到最小的准则函数来实现图像复原的,即

    v2-0c7eb8e95358347bf62c1ae78b9951f7_b.jpg


    式中,E()代表求期望值。因此维纳滤波又称为最小均方误差滤波器。

    维纳滤波需要假定下述条件成立:

    1. 系统为线性空间移不变系统。
    2. 退化图像、原始图像、噪声都是均匀随机场,噪声的均值为零,且与图像不相关。

    维纳滤波的复原滤波函数,即滤波器的传递函数为:

    v2-ad26c61383903986ab0313a85b233589_b.jpg

    没有噪声时,维纳滤波退化为逆滤波。有噪声时,维纳滤波利用信噪功率比对恢复过程进行修正,在信噪功率比很小的区域内,P(u,v)的值也很小,这使恢复图像较小地依赖于退化图像。在H(u,v)很小或等于零时,P(u,v)的分母不为零,维纳滤波没有病态问题[7]。

    在实际系统中,维纳滤波经常用下式近似:


    v2-94ed7a4b14ef5337c78205d2e0dc32e3_b.jpg

    ,K为特殊常数。 ②

    ①的维纳滤波要求未退化图像和噪声的功率必须是已知的。虽然用②近似的方法能得到好的结果,但功率谱比常数K的估计一般没有合适的解[8]。


    MATLAB语言程序

    clear;

    I=imread('C:\ok\原始图.jpg');

    imshow(I);

    I=rgb2gray(I); %将原图像转化为黑白图

    figure;

    subplot(2,2,1);

    imshow(I);

    title('转成黑白图像');

    [m,n]=size(I);

    F=fftshift(fft2(I));

    k=0.0025;

    for u=1:m

    for v=1:n

    H(u,v)=exp((-k)*(((u-m/2)^2+(v-n/2)^2)^(5/6)));

    end

    end

    G=F.*H;

    I0=real(ifft2(fftshift(G)));

    I1=imnoise(uint8(I0),'gaussian',0,0.001)

    subplot(2,2,2);

    imshow(uint8(I1));

    title('模糊退化且添加高斯噪声的图像');

    F0=fftshift(fft2(I1));

    F1=F0./H;

    I2=ifft2(fftshift(F1));

    subplot(2,2,3);

    imshow(uint8(I2));

    title('全逆滤波复原图');

    K=0.1;

    for u=1:m

    for v=1:n

    H(u,v)=exp(-k*(((u-m/2)^2+(v-n/2)^2)^(5/6)));

    H0(u,v)=(abs(H(u,v)))^2;

    H1(u,v)=H0(u,v)/(H(u,v)*(H0(u,v)+K));

    end

    end

    F2=H1.*F0;

    I3=ifft2(fftshift(F2));

    subplot(2,2,4);

    imshow(uint8(I3));

    title(维纳滤波复原图);

    运行结果如下:

    原始图:


    v2-8fd8293c2c4e56d36c5aecbf3be45aef_b.jpg


    复原后图像:


    v2-90049722dd3bee3f67223883b43c893e_b.jpg


    经过仿真,如上图所示,可以看出逆滤波复原与维纳滤波复原的区别和联系。维纳滤波后虽然仍有一些噪声存在,但已经和原图很接近了。因为原图像和噪声函数都是已知的,可以正确的估算参量。

    2.5使用Lucy-Richardson算法的迭代非线性复原

    L-R算法是一种迭代非线性复原算法,它是从最大似然公式印出来的,图像用泊松分布加以模型化的。当下面这个迭代收敛时模型的最大似然函数就可以得到一个令人满意的方程:


    v2-e087fa7f868abff27d93c94f39cd208f_b.jpg

    *代表卷积,

    v2-096dc38884e6d32999b6c9446733fecd_b.jpg

    代表未退化图像的估计,g和h和以前定义一样。这个算法的本质是显而易见的。它的非线性本质是在方程右边用

    v2-d1eb373a068aa5ddf22a09b1e2acdc9f_b.jpg

    来除产生的[4]。

    在IPT中,L-R算法是由名为deconvlucy的函数完成的,此函数的语法为

    fr= deconvlucy(g,PSF,NUMIT,DAMPAR,WEIGHT)

    其中,fr代表复原的图像,g代表退化的图像,PSF是点扩散函数,NUMIT为迭代次数(默认为10次),DAMPAR是一个标量,它指定了结果图像与原图像g之间的偏离阈值。WEIGHT是一个与g同样大小的数组,它为每一个像素分配一个权重来反映其重量。

    L-R算法程序:

    I=imread('C:\ok\苹果.jpg');

    PSF=fspecial('gaussian',5,5) ;

    BN=imnoise(Blurred,'gaussian',0,V);

    luc=deconvlucy(BN,PSF,5);

    figure

    subplot(2,2,1);

    imshow(I);

    title('原始图像');

    subplot(2,2,2);

    imshow (Blurred);

    title('模糊后的图像');

    subplot(2,2,3);

    imshow (BN);

    title('加噪后的图像');

    subplot(2,2,4);

    imshow (luc);

    title('恢复后的图像');

    模拟实验结果如下:


    v2-1f502a7cfd1a824b8ff846c89a21ef48_b.jpg


    用Lucy-Richardson算法可以较好的恢复图像[1]。但由于迭代产生的噪声痕迹是最大化可能性数据逼近法的常见问题,在低信噪比条件下,恢复图像可能会出现一些斑点,这些斑点并不代表图像的真实结构,只不过是恢复图像过于逼近噪声所产生的结果。另外此法存在一些较严重的缺陷,一是噪声放大问题,二是对于恢复图像中的不同部分,分别执行多少迭代才合适的问题。因为图像中信噪比高的部分可能需要数百次迭代才能获得满意的结果;而另一些光滑的对象可能只需很少次数即可达到满意的结果,所以适当选择迭代次数对图像恢复也很重要。这两个问题若得不到解决,将会对最终结果产生不利影响。

    2.6 盲去卷积

    通常图像恢复方法均在成像系统的点扩展函数PSF已知下进行, 实际上它通常是未知的. 在 PSF未知的情况下, 盲去卷积是实现图像恢复的有效方法。因此,盲去卷积算法就是那些不以PSF知识为基础的图像复原的方法。

    在过去的20年里,一种盲去卷积的方法已经受到了人们的极大重视,它是以最大似然估计(MLE)为基础,即一种用被随机噪声所干扰的量进行估计的最优化策略。简要的说,关于MLE方法的一种解释就是将图像数据看成随机量,它们与另外一族可能的随机量之间有着某种似然性。似然函数用

    v2-689ce1c8698be4cff32f38e9e5ecc513_b.jpg

    v2-99ec3b12c2b69383f597e42cb84587e2_b.jpg

    v2-ee7d05d064a7847e276039816d0ada4b_b.jpg

    来加以表达,然后,问题就变成了寻求最大似然函数。在盲去卷积中,最优化问题规定的约束条件并假定收敛时通过迭代来求解,得到的最大

    v2-355a343a958724d9a7dfa2d67a3fbf8a_b.jpg

    v2-ee7d05d064a7847e276039816d0ada4b_b.jpg

    就是还原的图像和PSF[3]。

    工具箱通过函数deconvblind来执行盲去卷积,它有如下语法:

    [f,PSFe]=deconvblind(g,INITPSF)

    其中,g代表退化函数,INITPSF是点扩散函数的出事估计。PSFe是这个函数最终计算到的估计值,fr是利用估计的PSF复原的图像。 用来去的复原图像的算法是L-R迭代复原算法[13]。PSF估计受其初始推测尺寸的巨大影响,而很少受其值的影响。

    盲去卷积程序:

    I=checkerboard(8);

    PSF=fspecial('gaussian',7,10);

    V=.0001;

    BlurredNoisy=imnoise(imfilter(I,PSF),'gaussian',0,V);

    WT=zeros(size(I));

    WT(5:end-4,5:end-4)=1;

    INITPSF=ones(size(PSF));

    FUN=inline('PSF+P1','PSF','P1');

    subplot(221);

    imshow(BlurredNoisy);

    title('A=Blurred and Noisy');

    subplot(222);

    imshow(PSF,[]);

    title('True PSF');

    subplot(223);

    imshow(J);

    title('Deblured Image');

    subplot(224);

    imshow(P,[]);

    模拟实验结果如下:


    v2-3839ab91663f5e37a0cccb222997ba52_b.jpg


    该算法优点是,同时恢复了图像和点扩张函数,在对失真情况毫无先验知识的情况下, 仍能实现对模糊图像的恢复操作。利用 MATLAB实现的图像恢复, 并对恢复图像的失真情况做了改善。在进行图像恢复时,重建 PSF,对图像进行重建, 得到恢复的图像。

    三、图像重建

    3.1 图像重建的概述

    图像处理中的一个重要研究分支是物体图像的重建,它被广泛的应用于检测和观察中,而这种重建方法一般是根据物体的一些横截面部分的投影而进行的。在一些应用中,某个物体的内部结构图像的检测只能通过这种重建才不会有任何物理上的损伤。由于这种无损检测技术的显著优点,因此,它的适用面非常广泛,它在各个不同的领域都显示出独特的重要性。例如:医疗放射学、核医学、电子显微、无线和雷达天文学、光显微和全息成像学及理论视觉等等领域都有应用。在医学影像处理中是医学图像获取的重要方法。如医疗放射学、核医学,电子显微等领域是必不可少的技术,在工业生产中的无损检测技术图像重建也扮演重要的角色。

    图像重建经过多年研究已取得巨大进展,产生了许多有效的算法,如:傅里叶反投影法、卷积反投影法、代数法、迭代法等,其中以卷积反投影法运用最为广泛。近年来,由于与计算机图形学相结合,把多个二维图像合成三维图像,并加以光照模型和各种渲染技术,已能生成各种具有强烈真实感的高质量三维人工合成图像。

    3.2傅里叶反投影重建

    傅里叶反投影重建方法可以说是最简单的一种变换重建方法,一个三维(或二维)物体,它的二维(或一维)投影的傅里叶变换恰好与此物体的傅里叶变换的主题部分相等,傅里叶变换的重建方法正是以此为基础的。该重建方法最早于1974年由Shepp和Logan提出,该方法是建立在“投影切片定理”这一理论基础之上的。根据此结论,可以先对投影进行旋转和傅里叶变换以构造整个傅里叶变换域中的各个方向的切片数据,然后再对其进行傅里叶反变换即可得到重建后的目标(原空间域中的图像)。

    傅里叶变换重建的原理如下:

    令f(x,y)代表一图像函数,则此二维函数的傅里叶变换为:

    v2-1fc660ae520eaa9a645fea9f6d222ef0_b.jpg


    而图像在x轴上的投影为:

    v2-2e3d9d0e12f6a7e08d33f01e65792159_b.jpg


    投影的一维傅氏变换为:

    v2-1aa1d7b48158eca640ab87a760ac996a_b.jpg

    ,它恰与二维傅氏变换的表达式一致。

    即:

    v2-5a2b757e30c25dac400b6cd3b6cb96e4_b.jpg


    现在假设将函数投影到一条经过旋转的直线上,该直线的旋转角度为。


    v2-0d5ab07221ee1a8a70a8690dc708c415_b.jpg


    t


    s


    v2-0910d3f6b0f5d31b82adcb00bf4d3fc4_b.jpg



    v2-047f57a5054bc3198f7cce1d23941f68_b.jpg



    v2-b849ab9adfc9bf35f81256a86410c07c_b.jpg













    定义旋转坐标为:

    v2-db99cb464418a21fdd25da01317ccc9e_b.jpg


    而将函数投影的直线选为x轴。投影点通过对距离t轴为处的一平行线进行函数积分,因此,该投影可如下表示:

    v2-acf11c58543e16a85eaac7cb2c838c43_b.jpg

    这里,积分路径是沿着

    v2-536cdea79bace87e2214916b6e8001fb_b.jpg

    直线进行的。此投影的一维傅氏变换为:

    v2-8ca221080d9880a8f8ed7809e0d09e16_b.jpg

    展开后为:

    v2-c6a3b44499ebf40c0f3bfdaad96de8d2_b.jpg


    为使展开式与投影的二维傅里叶变换相等,把指数项做某种代换得到:

    v2-5e33d8ac893a0fc9312d52bcbbbf3505_b.jpg


    因而,若点(u,v)在一条角一定而距原点距离为r的直线上,投影变换将于二维变换中的一直线有相同的傅氏变换,即:

    v2-907657a0cfdd4adaabcb5081446fea29_b.jpg


    若投影变换

    v2-b7c8d5f018c031fff552df4e931839e2_b.jpg

    中的所有

    v2-59cd8c66d273630540b93a60053d3dc5_b.jpg

    v2-ba51b8bbcf5fe2ab95d3428f2ed00a22_b.jpg

    值都是已知的,则图像的二维变换也是可以确定的。为得到图像函数,我们必须进行反变换运算,即:


    v2-184a0ff3257abda8cacba4a7b51cf019_b.jpg


    这些结论很容易推广到三维情形中。令:

    v2-23cd5c26a36c77c0ecd6a018e5ade87e_b.jpg

    表示一物体,这里f可以为实数或负数。它的三维傅氏变换由下式给出


    v2-41c27a8bae84c9fa12d6416fbe4474a5_b.jpg


    而变换的核心部分是:

    v2-741a322b10c557838c3923341372c293_b.jpg


    通过定义,纵剖面或在

    v2-b41a8fdf93ab84d8d723c22d82312571_b.jpg

    面上的投影是:

    v2-d29b16a4704ec044b0c2585128f84741_b.jpg


    注意到

    v2-58663006c37d4a40f01ad4423f4b65fa_b.jpg

    的二维傅里叶变换正好等于上述三维变换的核心部分。这也说明了如果投影在

    v2-999f338a3d72534350443ef15090c535_b.jpg

    平面上旋转了

    v2-d6793025889efad959721c09ffddcc6f_b.jpg

    角度,相应的傅里叶变换部分也将在变换域内的

    v2-761012ae59744987550b78fca1ffccc7_b.jpg

    平面内转过

    v2-89f76d021f4940f364596505044fc347_b.jpg

    角。这样,投影可以采用不同的方向角

    v2-7a89f3f3db632be38f7dc55008b66e6c_b.jpg

    插入到三维变换域中。建立一个傅里叶变换空间需要很多的投影。最后通过傅里叶反变换重建图像

    v2-9f508a550a8fce7b0c67b9cb5f3963f5_b.jpg

    。既然在三维空间中的任意平面都可以被重建,那么,一个二维图像

    v2-1b893337b1f49e6c3fa7e8ad1f781920_b.jpg

    的重建也不失一般性。

    我们可以重写二维投影方程,定出

    v2-4282f15edc204707fd570dfaffac6a65_b.jpg

    即投影平面

    v2-d8a914a6ac73d53c26abe3ecc47fdb29_b.jpg

    v2-83697be12cf1212374147666b31ffd64_b.jpg

    ,这里ds是光线几何路径中的微分长度。

    傅里叶变换的结论由下面给出:


    v2-da94df915f7db6a51c719661a708dc0d_b.jpg



    v2-d8c2d0f910c26dfd5bbbf3199453cacc_b.jpg


    若已知无数的投影,从极坐标

    v2-2f83efd90c5bebbb2c8b78140c5feab1_b.jpg

    中计算得到的投影变换推出在矩形平面

    v2-49c195cf9c2be4fed6a87823acf4535c_b.jpg

    中的傅里叶变换并不困难。但是,若只有有限个投影是有效的,则可能需要在变换中插入一些数据。另外需要注意的是,虽然只需一维傅里叶变换的投影数据就可构成变换空间,但图像重建则需要二维反变换。由此,我们得出一个结论,即:三维图像不能在得到部分投影数据的过程中局部的重建,而必须延迟到所有投影数据都获得之后才能重建。

    3.3 卷积法重建

    卷积法重建的基本思路和方法:

    逆投影原理:从各个方向得到的投影逆向返回到该方向的各个位置,如果对多个投影方向中的每个方向都进行这样的逆投影,就可能建立平面上的一个部分。典型的方法是卷积逆投影重建。

    卷积重建法是一种变换重建法,可以根据傅里叶变换投影定理推出。

    按照二维傅里叶反变换标准定义,有

    v2-4f7d2ecbcae5078f54592bf0a1a99fb7_b.jpg

    作代换:

    v2-cae19eefbbec490f4c8b40b8926cef53_b.jpg

    v2-646f9f1172dab8090a4f834e1e527f03_b.jpg


    写成极坐标(R,)的形式:

    v2-2259d8867861eb211eaaee0ee0c520ec_b.jpg

    利用傅里叶变换共轭对称性,有:

    v2-db92d3e59d2620a0928088019e825fbb_b.jpg

    令:

    v2-185c27208591d6868d1c44ca12746e4e_b.jpg

    则③可以表示为:

    v2-0e19ecf7eb35e76a1b0d414723f1b104_b.jpg

    , ④

    在②中,当用FFT计算投影数据的傅里叶变换F(R,)时,投影数据g()总被有限截断。当的采样间隔为d时,在变换域R的变化范围为

    v2-d71482d48a3d6d60d7beef40304ce06a_b.jpg

    v2-b91547b8965379ac1ec13ffda5038e75_b.jpg

    ,于是投影反变换重建公式可以近似写成:


    v2-38cb414f5fa0b233f425756f9062fd68_b.jpg


    采用标记:

    v2-6804a6784e3e20240db93889ac34dedd_b.jpg

    根据前式③,结合傅里叶投影定理可知:


    v2-2d092332456c80a226c7cffd36255dca_b.jpg


    由上式可以得出,要实现对已经得到的投影数据实现图像重建,则可以采取两步:首先将投影数据

    v2-2ba1fa7deb6922d70a0f3c5ed38b270e_b.jpg

    和相应脉冲滤波器⑤进行卷积,然后由式④对不同旋转角求和,就能实现图像重建。

    卷积可以看成一种滤波手段,卷积投影相当于先对数据滤波再将结果逆投影回来,这样可以使模糊得到校正。


    MATLAB程序:

    P=phantom(256);

    subplot(2,2,1);

    imshow(P);

    title('大脑的幻影图');

    theta1=0:10:170;

    [R1,xp]=radon(P,theta1);

    theta2=0:5:175;

    [R2,xp]=radon(P,theta2);

    theta3=0:2:178;

    [R3,xp]=radon(P,theta3);

    figure;

    imagesc(theta3,xp,R3);

    title('大脑幻影图的90条投影光速的Radon变换');

    colormap(hot);

    colorbar;

    xlabel('\theta');

    ylabel('x\prime');

    I1=iradon(R1,10);%R1有18条投影光速

    I2=iradon(R2,5);% R2有36条投影光速

    I3=iradon(R3,2);%R3有90条投影光速

    figure(1);

    subplot(2,2,2);

    imshow(I1);

    title('用R1重建图像');

    subplot(2,2,3);

    imshow(I2);

    title('用R2重建图像');

    subplot(2,2,4);

    imshow(I3);

    title('用R3重建图像');


    v2-987f6e9ecafd2d533039f7021cf66d20_b.jpg




    v2-82ffd210c1801a19759575ea652f2264_b.jpg


    通过MATLAB演示我们可以看出:用卷积逆投影来实现,数据质量高的情况下可重建出准确清晰的图像,而用傅里叶反变换重建法实现,不是很容易,且重建的图像质量很差。


    3.4代数重建方法

    代数重建技术就是事先对未知图像的各像素给予一个初始估值,然后利用这些假设数据去计算各射线穿过对象时可能得到的投影值(射影和),再用它们和实测投影值进行比较,根据差异获得一个修正值,利用这些修正值,修正各对应射线穿过的诸像素值。如此反复迭代,直到计算值和实测值接近到要求的精确度为止[5]。

    具体实施步骤:

    (1)、对于未知图像各像素均给予一个假定的初始值,从而得到一组初始计算图像;

    (2)、根据假设图像,计算对应各射线穿过时,应得到的各个相应投影值,,……,;

    (3)、将计算值和对应的实测值进行比较,然后取对应差值作为修正值;

    (4)、用每条射线的修正值修正和该射线相交的诸像素值;

    (5)、用修正后的像素值重复1~4各步,直到计算值和实测值之差,即修正值小到所期望的值为止。

    只要所测的射线投影值组成一个独立的集合,那么代数重建便将收敛于唯一解。


    参考文献

    [1] 龚声蓉,刘纯平,王强.数字图像处理[M].北京:清华大学出版社,2006.46-84,123-144[2] 张德丰,雷晓平.MATLAB基础与工程应用.北京:清华大学出版社,2012.

    [3] 冈萨雷斯等.数字图像处理(MATLAB版)[M].北京:电子工业出版社,2005.5:123-133

    [4] 张铮,倪红霞,苑春苗,杨立红.精通MATLAB数字图像处理与识别[M].北京:人民邮电出版社,2013.8:194-196

    [5] 陈传波,金先级.数字图像处理[M].北京:机械工业出版社,2004.119-124,161-176

    [6] 林旭梅,葛广英.MATLAB实用教程[M].山东:中国石油大学出版社,2010.3:49-75

    [7] 张德丰,张葡青.维纳滤波图像恢复的理论分析与实现[J].中山大学学报,2006.45(6):44-47

    [8] 何东健.数字图像处理[M].西安:西安电子科技大学出版社,2003.261-279

    [9] 徐昕.MATLAB工具箱应用指南:控制工程篇[M].北京:电子工业出版社,2000.1-12

    [10] 罗军辉,冯平.MATLAB7.0在图像处理中的应用[M].北京:机械工业出版社,2005.257-277

    [11] 赵怡红,张常年编著.数字信号处理及其MATLAB实现.北京:化学工业出版社,2002.1:38

    -125

    [12]徐晰.MATLAB工具箱应用指南:控制工程篇[M].北京:电子工业出版社,2000.1-12

    [13]刘维一,于德月等.用迭代法消除数字图像放大后的模糊[J].光电子.激光,2002,13(4):389-400

    [14]陶洪.数字共焦显微技术及其图像复原算法研究四川大学硕士学位论文CNKI::CDMD:10610.2.2003.6632

    [15]苏开亮.运动模糊图像的恢复及恢复质量评价[D].西安:西安电子科技大学出版社,2010,21-25

    展开全文
  • matlab超声波原始码超声发声在mTBI / PTSD和神经性疼痛动物模型中的情感状态分析 阿比盖尔·辛德勒(Abigail Schindler)(主管),瓦伦蒂娜(Valentina Staneva)(科学主管),迈克尔·帕克(Michael Park),本·...
  • MATLA技巧

    2021-01-14 20:15:45
    如何将MATLAB矩阵每列都归一化为单位向量? function matrix_normalized = norm_matrix_by_column(matrix) column_norms = sqrt(sum(abs(matrix).^2, 1)); % or vecnorm(matrix) matrix_normalized = matrix ./ ...
    如何将MATLAB矩阵每列都归一化为单位向量?
    function matrix_normalized = norm_matrix_by_column(matrix)
    
        column_norms = sqrt(sum(abs(matrix).^2, 1)); % or vecnorm(matrix)
        matrix_normalized = matrix ./ column_norms;
        
    end
    
    展开全文
  • 基于粒子群算法的PID参数寻优(MATLA程序),利用混合粒子去算法寻优参数。
  • matla错误

    2020-06-17 22:50:46
    matlab出现:(Undefined function ‘函数名’ forinput arguments of type ‘double’) 查资料的时候看到过两种解决办法 我是遇到了第一种情况 一 检查函数名和你引用这个函数的时候的名字一定要一样 ...

    matlab出现:(Undefined function ‘函数名’ forinput arguments of type ‘double’)

    查资料的时候看到过两种解决办法
    我是遇到了第一种情况

    一 检查函数名和你引用这个函数的时候的名字一定要一样

    二 在file->set path->add floder在这里选择你M文件所在的目录

    原文链接

    展开全文
  • matla绪论

    2020-03-30 13:20:31
    绪论 1.1简介 MATLAB是matrix&laboratory两个词的组合,意为矩阵工厂(矩阵实验室)。是由美国Mathworks公司开发出品的商用数学软件。可用于数据分析处理,算法开发,数据可视化以及进行图形科学处理的高级...
  • matlab超声成像代码HDR-美国 在这项工作中,我们将 HDR 技术应用于超声成像,我们将在不同功率级别获取的超声图像结合起来,以提高最终图像中可见的细节水平。 我们的结果强烈表明 HDR-US 成像可以提高超声在基于...

空空如也

空空如也

1 2 3 4 5 ... 20
收藏数 3,862
精华内容 1,544
关键字:

MATLA