登陆注册
7554300000016

第16章 基于小波变换的PET图像处理算法(4)

为明确,我们首先以一维的信号为例说明小波算法。对于一维信号的小波变换就是用高通和低通小波滤波器分别对信号进行卷积,尺度函数对应的滤波器为低通滤波器,对于db2小波其滤波器的非零参数即为上面我们计算出的值,一般记为h(n)或hn,小波函数对应的滤波器为高通滤波,其参数可以通过低通滤波器的参数计算出来,记为g(n)或gn。而相应的重建滤波器即在前面加上I表示Ihn和Ign。

定理3db2小波滤波器参数之间的关系定理

定理中对每一个参数分母乘上2,是为了对参数进行归一化。对于滤波器与信号的卷积操作我们可以通过将滤波器参数反向而转化为内积操作。

设原始信号为{s(n)},它被低通滤波器的滤波操作可用下式表示:

这里L为滤波器的长度,对于db2小波L=4;{s1(j)}为滤波以后的信号。根据Nyquist采样定理:

Fs(=1/Ts)≥2Fmax或Fs≥2Bs

式中:Ts为采样周期;

Fs为采样频率;

Fmax为原始信号的最高频率;

Bs(=Fmax-Fmin)为原始信号的带宽。

滤波后的信号带宽减少了一倍,因而我们对信号的采样频率也可减少一倍,所以对于滤波后得到的信号{s1(j)}我们可以对其进行下采样(Subsample),即只保留信号中的奇数位置上的点或偶数位置上的点,这样可以将信号的采样点减少一半,从而降低计算的复杂性。

下面我们用线性代数的方法来进行小波变换,虽然实际的信号采样点是有限的,但我们首先对理想的无限长度的信号进行讨论,采用的是db2小波。

在进行反变换时,我们用的滤波器对可从定理3中得到,其他不同类型的小波变换所用的方法与上面的方法基本类似,只是采用的滤波器参数的值和长度不同罢了。

对于有限的信号,以我们采用的db2小波对一共有八个采样点的一维信号进行正向变换为例,由于我们进行变换的信号采样点是有限的点,而滤波器长度为4所以在进行内积操作时在n-2的地方会超出信号的长度,这种效应在Haar小波中是不存在的,因为Haar小波的滤波器参数长度为2。

有限长度信号的小波反变换也存在这种边界问题:

在下面一节我们将对如何处理有限长度信号的边界问题进行一些讨论。

对于二维PET图像信号的处理与一维信号的处理相似,只是要分别对行像素和列像素用滤波器进行滤波。由于小波滤波器是成对使用的,如果行像素和列像素均采用低通滤波器滤波,得到的子图为概貌图像;如果行像素采用低通滤波器滤波而列像素采用高通滤波器滤波,得到的子图为水平细节(水平边缘)图像;如果行像素采用高通滤波器滤波而列像素采用低通滤波器滤波,得到的子图为垂直细节(垂直边缘)图像;如果行像素和列像素均采用高通滤波器滤波,得到的子图为对角细节(对角边缘)图像。

实际在进行PET图像的小波变换时,我们可以有以下两种分解方法:

第一种分解方法在滤波后未进行下采样的操作,因此每一次分解后的四个子图大小与原图像大小一样,四个子图像像素个数总和被扩大为原来图像的四倍,计算复杂度增加。根据采样定理滤波后各子图的带宽减小的原理,不需要采用与原图像相同的采样频率,因此采用这种方法给出的结果是含有冗余信息的。

第二种方法在对行和列进行滤波操作时均进行了下采样的操作,因此得到的各个子图像的大小都为原图像大小的1/4,而四个子图像像素的总和等于原图像的大小,这在信息表示是无冗余的,并且这种表示方法也不会丢失原图像的任何信息。我们在进行医学图像处理时,由于图像一般都比较大,计算速度是一个重要的算法评价指标,因此在本书中我们采用的是第二种分解方法。

8.6.3PET图像处理中小波变换的边界问题

仔细分析小波变换的算法可以发现,我们所要处理的原始真实离散信号S总是只有有限个非零点,而对实际信号实施小波分解时要进行下采样,在对分解后的信号进行合成时又要先实施插入零操作。本书所处理的PET图像信号也是一个长度有限的信号,它的边界效应也很严重,所以在进行小波变换时需要解决PET图像的有限边界问题。

目前,解决边界问题的方案大致有两类,一是假设非零样本点以外的点为零点;二是根据已知实际信号向外延拓。在图像处理领域,由于样本点的边界条件是固定不能改变的,为了能够重构原来的信号,即要使重构后的信号与原信号具有同样多样本点数,就必须考虑边界延拓问题。一般采用对称延拓技术,即假设信号关于信号两边缘的点对称。将延拓的信号用表示为:

对于非周期信号不论采用哪种延拓方案,都无法从采样截断得到数据延拓与被截断的实际信号相同的信号,因而经过小波分解或重构得到的信号的边界附近的点都会受到影响。受影响的样本点数与滤波器的选择直接相关。对于基于Daubechies所给出的滤波器系数,小波变换结果受影响的点数与滤波器的长度成正比关系,而且分别对离散逼近信号的右边界和离散细节信号的左边界接受边界效应的影响。通过上述分析可以看出,由于边界效应的影响,应用小波分析的信号不宜过短。另外因为受边界效应影响的点数随滤波器长度的增大而增多,所以边界效应的影响也应该作为选择小波基函数的一个考虑因素。在本书中由于图像信号的长度不会太长,为了尽量避免边界效应我们采用了滤波器长度较短的db2小波,这在前面已经介绍了。

8.6.4模拟PET图像的激活区提取结果与讨论

1.激活度为5%时的模拟计算结果

两个结果均是在检验水平为0.01的情况下得到的。比较这两幅图可以看出在这种条件下直接在空间域中检验得到的激活区域由于受噪声的影响,只得到一些孤立、杂乱的激活点,无法得到类似的激活区域。而采用小波域方法检验得到的结果已明显地形成了六个不同的激活区域,即使是半径为1的激活区域也被准确地检测到了。并且可以看到小波域上检验得到的激活区域大致是成片的,不像是一些孤立的点,这一优点是由小波变换的多尺度(Multi-scale)特性决定的,小波变换后的图像相当于是从大尺度上去观察一幅图,可以避免细节和高频噪声的影响,变换层数越多,尺度越大,这就像我们从远处看物体可以看到全貌,而看不见细节一样。

为比较两种方法对激活区检测的灵敏度,列出了在激活度为5%、噪声方差为0.001时不同检验水平下两种方法得到的结果中犯Ⅱ类错误的像素个数,即未被检测到的激活区域的像素个数。从表中可以看出在不同的检验水平下小波域方法犯Ⅱ类错误的像素个数均比空间域上的方法少,并且随着检验水平的上升小波域方法犯Ⅱ类错误的像素个数迅速下降。这说明小波域方法是一种灵敏的检测图像序列激活区的方法。

同类推荐
  • SEO智慧

    SEO智慧

    本书不仅是系统讲述搜索引擎优化(SEO)的原理、技法、思维、经验的畅销作品,而且论述了市场学、营销学、心理学怎样与搜索营销形成综合战略的经典著作。从说明网站的搭建该怎样围绕搜索引擎优化开始,到网站优化的具体策划和技术细节,再到搜索营销的持续开展,说明打造投资回报率高的网站应该怎样具备必需的搜索营销思维和战略做法,是作者12年网站营销经验,特别是搜索引擎优化(SEO)实现营销方面的技术和实践总结。
  • 我看电商

    我看电商

    本书是作者近30年从事零售及电子商务管理的总结和分享。近年来电商行业在中国迅猛发展,2012年网络零售市场规模达到13,000亿人民币,诞生了淘宝、天猫、京东、当当、凡客、唯品会……等一大批全新的网络公司,电子商务正在日益深入的影响着越来越多人的生活。如今这位自称跨界老兵的电商操盘手第一次全方位剖析电商风云,细说中美电子商务发展的不同路径,评点国内各大电商企业的是是非非,详解其中的融资、模式、运营、效率……个中经验,毫无保留的娓娓道来,笔触至深,让人不禁掩卷沉思。
  • 如何处理电脑故障

    如何处理电脑故障

    本书以问答的方式介绍了电脑会出现的各种故障,内容包括了音箱声音失真,如何处理?如何做好电脑的日常维护等等问题。
  • 计算机与人工智能(科学新知丛书)

    计算机与人工智能(科学新知丛书)

    本书是专门为中小学生编写的科普读物。全书共分八个部分,分别就计算机诞生的历程,计算机的种类、构造、作用及其发展趋势等方面进行了详尽的论述,全书论述深入浅出,行文活泼流畅。为本书的主要读者——广大的中小学生朋友认识计算机、了解计算机和掌握计算机提供了一把钥匙,是新时期中小学生开展素质教育不可多得的科普读物,中小学生家长及教师阅读此书也能从中获得有益的启示和感悟。
  • 体感交互技术

    体感交互技术

    本书提出了体感计算的研究,探讨了基于体验性认知的个性化电子学习方式,并且分析了面向儿童的艺术创造体验,在体感虚拟化身的帮助下,用户可以通过身体运动,促进情感心智的全面发展,
热门推荐
  • 复兴在灭亡之后

    复兴在灭亡之后

    我们都生活在这个平凡的世界里,有的时候感觉一切已经注定,从生下来那一瞬间开始。如果有一天,我们的生活不再平凡。如果有一天,一切的秩序重新的建立。你,期待吗?
  • 传奇附身

    传奇附身

    他,黑道中人闻风丧胆!谈之色变!失去“亲人”之后,一心想过平凡生活,继续寻找时光隧道,回古代报家门之仇的他,却在迫不得已的情况下,改变了身份和姓名。重出江湖,带着曾经的兄弟一路拼杀……站在了Z国黑道金字塔顶端时,只有他知道这条路是痛?是血?还是泪?
  • 好想爱这个世界啊致郁

    好想爱这个世界啊致郁

    对自我的认知,对生活的态度,对世界的感知……感谢有你。我不想爱这个世界。但我想陪你一起路过这个世界。——邢雨——悼念我,忘掉我,逃离我;寻找我,放掉我,靠近我。人生至幸就是确信自己被爱。彼此救赎,感知被爱。——华果——无论我,在这里,在那里。那些旁白,那些姿态,那些伤害。
  • 都市热血保镖

    都市热血保镖

    李延庭受何雪梅所托,刚到城里担负保护华莱尔公司安全的职责,就遇到了仗势欺人的保安马文将,以及高冷自傲的黄琛,被他们阻挠并且奚落。何雪梅却出现与李延庭相认姐妹,并且在饭店吃饭时,遇到酒鬼调戏服务员,李延庭出手相助。
  • 当二次元总裁穿越到三次元世界

    当二次元总裁穿越到三次元世界

    【脑洞短篇,非太监文】龚俱仁,是一部沙雕漫里的霸道总裁,有一天,他“不小心”打破了次元壁,然后……
  • 光下影子

    光下影子

    你有没有注意过你的影子呢?黑黑的一团,形影不离地追随着你,可是如果我爱你,你会接受我吗?
  • 天行

    天行

    号称“北辰骑神”的天才玩家以自创的“牧马冲锋流”战术击败了国服第一弓手北冥雪,被誉为天纵战榜第一骑士的他,却受到小人排挤,最终离开了效力已久的银狐俱乐部。是沉沦,还是再次崛起?恰逢其时,月恒集团第四款游戏“天行”正式上线,虚拟世界再起风云!
  • 英雄联盟一曾经的故事

    英雄联盟一曾经的故事

    符文之地曾经的回忆:每个英雄都有一个属于自己故事…从“埃卡西亚”到“恕瑞玛”再到如今“瓦洛兰”大陆的每个一个国家。战争的硝烟四处弥漫所到之处除了死亡或许就只剩下仇恨了……然而“英雄联盟”的成立将会为这个世界带来怎样的改变…那就让我带你们一起去见证那段不为人知的故事吧!
  • 三殇缘

    三殇缘

    一世:她说情深缘浅,来世再相见,他提笔写半阕相思。二世:他说若得佳人心,江山为聘,后来香消玉陨,帝王何心?三世:不管天上人间,只愿生死毋忘,两情长。
  • 新大众菜谱:微波炉美味600款

    新大众菜谱:微波炉美味600款

    本套菜谱共30本,汇集了南北方以及各种风味的菜系,每本600余种做法。本书介绍了微波炉菜的做法,简单好学易做,是符合大众口味的家居生活常备书籍。