《基于连接突触计算网络的医学图像融合算法-高媛.pdf》由会员分享，可在线阅读，更多相关《基于连接突触计算网络的医学图像融合算法-高媛.pdf（5页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、Joumal of Computer AppIications计算机应用，2017，37(12)：35543557，3568ISSN 1001908lCODEN JYUDU201712一10http：wwwjocacn文章编号：1001908l(2017)12355404 DOI：1011772jissn100l-90812017123554基于连接突触计算网络的医学图像融合算法高 媛，贾紫婷，秦品乐，王丽芳(中北大学大数据学院，太原030051)(通信作者电子邮箱949763497qqcom)摘要：针对传统的脉冲耦合神经网络(PcNN)融合方法中参数过多，以及参数和网络迭代次数难以准确设置、
2、融合效果差等缺点，提出了一种用连接突触计算网络(LscN)模型的连接项(L项)进行图像融合的算法。首先，把两幅待融合图像分别输入到【scN模型中；其次，使用L项代替传统PcNN中的点火频率作为输出；然后，使用多通工作方式终止迭代；最后，通过比较L项的值得到融合后图像的像素。理论分析与实验结果表明，与改进的PcNN模型和在PCNN模型的基础上提出的新模型进行图像融合的算法进行比较，所提算法得到的融合图像更有利于人眼观察；特别是与点火频率作为输出的LscN方法相比，所提算法在边缘信息评价因子、信息熵、标准差、空间频率、平均梯度上均较优。该算法简单易行，不仅减少了待定参数数目，降低了计算复杂度，而且
3、解决了传统模型中迭代次数难以确定的问题。关键词：连接突触计算网络；点火频率；图像融合；计算机断层扫描；磁共振成像；脉冲耦合神经网络中图分类号：TP3914l 文献标志码：AMedicm image fhsion algOritll】m b嬲ed On li州I【ing synaptic computation networkGAO YuaIl。，JIA Ziting，QIN Pinle，WANG Lifang(&00f 0，眈比sc如，lce，舶衲，耐wm妙矿吼i，泖nn Skwi 03005l，傩im)A瞻tnct：I。he traditional lusion methods based
4、on Pulse Coupled Neural Network(PCN N) haVe the shortcomings of toom锄y p眦meters，tle par锄eters and number of network itemtions dimcult to accurately set，and p0哪fusion eectIn orderto solVe the pmblems， a new image fusion algoritlm using the connection item(L item) of Linking Synaptic ComputingNetwork(
5、LSCN)model was pmposed Firsdy， the two images to be fused was input into the LsCN model respectivelySecondly， tJle L te册w鹊used to replace t11e i伊ition frequency in the叽ditional PcNN as the outputThen，the itemtion wastem“nated by the multi-pass opemtion Finally， the pixels of the fused image were obt
6、ained by compdng the values of Lte册sThe tIleoretical analysis and experimental results show that， compared with the image fusion algorithms using theimpmVed PCNN model and the new model pmposed on the b鹊is of PCNN model，the fusion images generated by the pmposedalgorithm have better visual effectsIn
7、 addiIion，compared诵th the fusion a190rithm of LsCN using ignition frequence as theoutput， the pmposed algothm is all supeor in edge info丌nation evaluation factor， inf0邢ation entmpy， standard deviation，space fhquency， average grads The proposed algorithm is simple and convenient， which not only reduc
8、es the number ofp锄硼eters to be detenllined， reduces the computational complexity， but also solves the pmblem that the number of iterationsin the traditional model is di仿cult to be detenllinedKey woIcIs：Linking Synaptic Computation Network(LSCN)； ignition fIequency； image fusion； Computed Tomography(
9、CT)；Magnetic Resonance Ima舀ng(MRI)；Pulse Coupled Neural Network(PCNN)0 引言目前国内图像融合领域研究的方法较多，如基于神经网络、基于多尺度分解2j、基于统计3 o的融合算法等。其中神经网络具有较强的生物学背景，能够考虑到图像的全局特性。一种特殊类型的神经网络来自猫的视觉皮层，称为脉冲耦合神经网络(Pulse coupled Neural Network，PcNN)Hj已经成功地被应用于图像融合，取得了不错的成果o。尽管如此，PCNN模型本身存在一些弊端，由于其复杂的结构与实际生物脑神经网络的运行机制类似，所以如何使大量参数、
10、迭代次数得到准确的设置以及降低计算复杂度仍然是迫切需要解决的问题哺。目前基于PCNN融合方法的改进主要包括两个方面：一种是改进PcNN模型自身，例如参数设置的改进。Qu等在非下采样轮廓波变换(Nonsubsampled contourletTransfo丌rI，NSCT)域内采用空间频率激励PCNN模型来融合图像，但该算法需要对图像进行NscT，提高了算法的计算复杂度。王俭等o采用粒子群算法优化PcNN中的参数，在很大程度上解决了参数设置问题，但粒子群算法的加入也增加了算法的复杂度。薛寺中等9。将小波变换与PcNN相结合融合图像得到了更好的视觉效果，但迭代次数设置为200次，浪收稿日期：201
11、706-0l；修回日期：201708-30。 基金项目：山西省自然科学基金资助项目(2叭5011045)。作者简介：高嫒(1972一)，女，山西太原人，副教授，硕士，主要研究方向：图像处理、人工智能；贾紫婷(1992一)，女，山西吕梁人，硕士研究生，主要研究方向：医学图像融合； 秦品乐(1978一)，男，山西长治人，副教授，博士，主要研究方向：机器视觉、大数据处理； 王丽芳(1977一)，女，山西长治人，副教授，博士，主要研究方向：机器视觉、大数据处理。万方数据第12期 高媛等：基于连接突触计算网络的医学图像融合算法 3555费了大量时间，而迭代次数过小又不能使PCNN得到充分的执行。另一种改
12、进是在PCNN理论的基础上提出新模型，如相交皮质模型(Intersecting Cortical Model，IcM)、脉冲发放皮层模型(Spil(ing Conical M0del，scM)1“，这些模型也已经应用于图像融合中2“J。这两种模型均简化了PcNN，减少了参数的数目，但文献10一11都没有解决迭代次数的问题。基于以上分析，为了降低传统算法的计算复杂度，减少待定参数的数目，解决网络迭代次数难确定的问题，本文提出了一种用连接突触计算网络(Linking synaptic computationNetwork，LscN)的连接项进行图像融合的算法。该算法具有更简洁的模型结构，相比传统算
13、法减少了待定参数的数目，降低了计算复杂度，并通过多通工作方式解决了迭代次数难以确定的问题。通过三组计算机断层扫描(computedTomog阻phy，cT)与磁共振成像(Maglletic Reson蚰ce Imaging，MRI)图像进行实验，结果表明无论从是主观还是客观上本文算法均有更优良的图像融合效果。1 LSCN基本结构zh蛐等驯于2017年1月提出了LScN模型。在设计该模型时，受Ecknom等1 6。思想的启发，为了更恰当地描述生物电传输的指数衰减和视觉皮层受到一定刺激时相邻神经元同步发放脉冲等现象，用反馈输入作为外部激励，并且用内部活动项本身表示神经元状态的衰减项。不同于前几代神
14、经网络，LscN与PcNN模型同属第三代神经网络的范畴，该代神经网络主要利用脉冲同步发放现象实现对相邻神经元的激发和抑制。LsCN模型和PcNN有两个主要的不同点：第一个不同点是在该模型中膜电位由漏电积分器表示；第二个不同点是反馈输入仅仅由激励表示。基本的PCNN模型涉及到馈送输入、连接输入、内部活动项、点火阈值以及脉冲输出5个方面的内容，需要设定的参数数目高达10个。所以，众多参数设置的准确性将对融合结果产生严重的影响，同时运行过程中会浪费大量的计算资源。相比之下，作为PcNN模型的改进模型，LscN在输入信号采集机制上不断简化，待定参数的总量也大大减少。传统的PCNN模型中有三个漏电积分器
15、，需要进行两次卷积运算。而在LScN模型中，也有三个漏电积分器，但只需进行一次卷积运算。这一点决定了LScN模型的时间复杂度低于传统模型，同时可以看出内部活动项与外部激励的关系更直接。不仅如此，与传统PcNN不同，LscN模型的迭代过程是自动停止而不是手动设置的，在多次迭代中，操作更加方便。LscN的基本神经元结构如图1所示。图1 LSCN基本神经元结构Fig l B聃ic neumnal st兀lcture in LSCNLsCN模型的数学表达式如下：F。(n)=S。岛(凡)=儿i(n一1)+形咖k(n一1)(凡)=i(，l1)+s#(1+肚。(n)pi(凡)=Di(n1)一6+(n一1)驰
16、)=矗搿朋“曲(2)(3)(4)(5)式中：，。，(n)为神经元的反馈输入；(i，J)表示每个神经元的索引；s。携带刺激信息；匕(n)为连接输入；(p，q)表示其相邻神经元；f为连接常数；形抽为应用于连接突触的权重；K(n1)是突触后动作电位。在本文中，膜电位由漏电积分器表示，(n)为内部活动项，其中厂是膜电位衰减常数，卢是连接强度；神经元的阈值由漏电积分器表示，9i(n)是点火阈值；占是一个小的正常数；是幅度调整；匕(n)为脉冲输出。2 本文算法21算法原理在LscN网络中，图像里每个像素点相当于一个神经元。当某个神经元的(n)Oi(凡)时，该神经元点火，发出一个脉冲。此时阈值9。(凡)通过
17、反馈输入迅速提高，致使该神经元停止发射脉冲。阈值随着时间开始衰减，当其再次小于内部活动项时，神经元再次点火，如此反复。神经元点火频率越大说明对应像素点的亮度越大。在图像融合中融合指标的好坏直接影响到融合结果。以往的基于PCNN及其改进模型的融合算法中均是根据若干次迭代后神经元的点火频率来决定融合后图像的像素。选择点火频数较大的那幅图像所对应的像素值作为融合后图像的像素值。点火频率如式(6)：肘F(n)=肘i(，l1)+K(，l1) (6)但在本文中，使用L项代替传统的点火频率作为输出。由式(2)与(6)可知，当f取l时，L项和点火频率相似。但L项在每次迭代中均对y作了卷积处理，实验证明L项代替
18、点火频率输出能够取得更好的效果。22算法流程假定cT图像为A。MRI为B，并且已经过几何配准处理。融合过程如下：1)分别将待融合图像A和B作为LsCN模型的反馈输入；2)初始值设置为k(o)=o，(o)=o，K(o)=o，(0)=0，pi(0)=l；3)分别按照式(1)一(5)，计算k(n)、(，1)、90(n)、匕(凡)；4)使用L项代替传统的点火频率作为输出；5)迭代使用式(1)一(5)，使用多通工作方式作为迭代终止条件；6)利用比较选择算子得到融合后图像的像素，如式(7)所示：眠j)：f抛儿驯，J)， (7)L占(；，j)， 厶(i，J)k(i，)式(7)中LscN的L项数值越大，表示该
19、像素点的亮度越大，所以选择该点为融合后图像中的像素点。在使用LscN模型处理图像时，迭代次数一直是一个棘手的问题，传统的算法中迭代次数总是通过大量的实验或者个人经验来设置。但迭代效果的好坏与迭代次数有很大的关系，如果迭代次数过大，将会浪费计算资源，迭代次数过少，将万方数据3556 计算机应用 第37卷会导致算法运行不充分。LScN模型中使用多通工作方式作为迭代终止条件，即当所有神经元都点火迭代终止。设，是神经元的总数，在每一次迭代中神经元点火的次数被统计，l，】，矩阵的初始值是0，其大小与y相同，yy和y之间是或运算关系，当l，y矩阵中所有元素均为1时迭代终止，也就是yy中l的个数总和为、，。
20、本文算法中LScN的其他参数为：神经元的连接范围为33，=02，厣=3，6=l255，f=l，=3，r=004 004 0041lo04 100 o04 I。本文算法中的参数除卢采用了经验值【o04 o04 o04 J外其他均是参考文献153 实验结果与分析本文分别使用正常脑部、脑膜瘤、脑弓形虫病的CT和MRI图像进行融合仿真实验：实验环境为：64位wjndows 7操作系统，Matlab R2叭5b，Intel core i5-2430，24 GHz处理器，8 GB运行内存。实验图片来自哈佛大学医学院。实验结果如图24所示。(e)文献【l 7】算 (f)scM (g)点火频率 (h)本文算法
21、融合结果 融合结果 输出融合结果 法融合结果图3脑膜瘤患者的cT与MRI图像及融合结果Fig3 CT and MRIimages and fusion results of menin百oma从图24的子图(a)和(b)中可以看出：cT图像中骨骼成像突出，但软组织对比度较低；MRI图像则不同，骨结构显示较差，但软组织分辨力高，更有利于对病灶位置的定位。将文献7算法、在NscT下LscN模型的融合算法、文献17算法、基于scM的融合算法、使用点火频率作为输出的LScN融合方法和本文算法进行比较。其中文献7算法为在NscT下基于PcNN模型的融合算法，该模型的输出为点火频率；NscT下LscN模型
22、的融合算法与文献7的融合框架相同，其输出为L项；文献17为没有分频的PcNN融合算法，该模型的输出为点火频率；基于scM的融合算法中，scM模型是在PCNN上改进而来的，其输出仍为点火频率；使用点火频率作为输出的LscN算法和本文中使用L项作输出的算法的区别仅在于输出项，目的是为了比较LSCN模型分别以点火频率和L项来指导融合的结果。传统的算法均是通过点火频率的大小来选取融合后图像的像素值。本文算法和分频的LScN算法是用L项的大小来决定融合后图像的像素值。将本文算法与四种传统算法进行对比，目的是比较分别使用L项和点火频率进行图像融合的结果。(a)cT图像 (b)MRI图像 (c)文献【7】算
23、 (d)分频LscN法融合结果 融合结果(e)文献【17】算 (f)scM (g)点火频率 (h)本文算法融合结果 融合结果 输出融合结果 法融合结果图4脑弓形虫病患者的CT与MRI图像及融合结果Fig 4 CT and MRJ images and fusion results ofbrajn t似opl鼬m09is分频的LScN算法、点火频率输出的LScN算法中参数和本文算法相同。参考文献18，scM的参数为：，=O2，g=0109 l 01409 0109 11o9，=20，w=l o1409 0 o1409 l。【o109 1 o1409 o109 1 J由图24可见：6种算法都成功地
24、融合了两幅图像，既突出了脑部的骨骼部分又包含了软组织的信息，取得了不错的效果。仔细观察可以发现：图2(c)和2(d)中视觉效果差不多，但是均丢失了源MRI中的部分信息。图2(e)中软组织和骨骼信息较暗，图像质量较差；图2(f)(h)的融合结果较好，但图2(h)的视觉效果最好，可以很好地保留图像的空间细节信息，骨骼信息突出，软组织对比度高，成像质量最好。图3(c)、(d)中在对MRI源图像信息的提取和表达上表现不佳，存在一些伪影，图3(e)视觉效果相对好一些，但是存在同样的问题。相比之下，图3(f)、(g)、(h)的清晰度、对比度较高，但本文算法中软组织部分与源图像较接近，融合效果最好。由图4中
25、脑弓形虫病患者MRI图像可知，图4(c)、(d)中细节信息丰富，但是图像灰暗难以观察；图4(e)中源MRI中的很多细节被弱化；图4(f)、(g)、(h)在视觉效果上优于其他对比算法，但本文算法中信息丰富成像质量最好。以上是从人眼角度的主观评价，为了更加全面地对融合结果进行评估，本文选取了客观评价指标互信息(MutualInfo丌11ation，MI)、边缘信息评价因子(Edge Info丌IlationEvaluation Factor，Q)、信息熵(Info咖ation Entropy，IE)、标准差(standard Deviation，SD)、空间频率(Space Frequency，S
26、F)、平均梯度(Average Grads，AG)(AG可衡量图像的清晰度)对融合图像进行评价。以上6个指标的值越大融合效果越好，评价结果如表l一3所示。万方数据第12期 高媛等：基于连接突触计算网络的医学图像融合算法 3557表l 正常脑部的客观指标Tab1 0bjective indicators of no丌rIal bmin由表1可以看出，本文算法中选用L项作为输出在Q“、IE、sD、sF、AG均高于点火频率作为输出时的值。本文算法有4项指标高于scM算法，有5项指标高于文献17算法；本文算法有4项指标高于分频后的LscN，其中IE、AG低于该算法；本文算法有4项指标高于文献7算法，I
27、E、AG低于该算法。表2脑膜瘤的客观指标Tab2 ObjectiVe indicators of menin西oma由表2可得本文算法在MI、Q、IE、sD、SF、Ac上均高于点火频率作为输出的值。本文算法有4项指标高于scM，1项与其相同，有5项指标的值高于文献17算法；本文算法有4项指标高于分频后的LScN，其中IE和AG均低于该算法；本文算法有4项指标高于文献7算法，其中在IE和AG上低于该算法。表3脑弓形虫病的客观指标Tab3 Objective indicator of brain tox叩l如mosis由表3可得本文算法在Q“、IE、sD、sF、AG均高于点火频率作为输出的值。本文
28、算法中5项指标均高于scM，有6项指标均高于文献17算法；本文算法有4项指标高于分频后的LSCN，其中IE和AG均低于该算法；本文算法有4项高于文献7算法，其中文献7算法中IE和AG高于本文算法。由表l3可知，本文算法和文献7算法、文献17算法、scM算法、以点火频率输出的LScN算法这些均是以点火频率作为指标来融合图像的算法相比，6项客观评价指标中均有大于等于4项优于对比算法。本文算法和分频LScN算法均是使用L项作为指标来融合图像，但本文算法也有4项指标明显高于对比算法。综上所述，本文算法取得了不错的结果，不仅能够较好地融合图像信息而且保留了细节部分，客观评价与直观观察结果保持一致。通过观
29、察文献7算法、分频LsCN算法实验结果可知，这两种算法在IE、AG上均明显高于其他算法。分频的LscN算法在IE、AG上均高于本文算法，这两种算法的区别仅在于分频与否，因此可得到一个结论：在利用LscN模型进行图像融合时，如果首先对图像进行NscT分频处理，分频后的LScN算法在IE和AG指标上较高；如果不分频，则该算法在MI、Q删7、sD、sF上较高。因此在使用该模型进行图像融合时可以针对融合图像的需求选择合适的算法。如果注重融合图像的IE和AG则可对图像进行分频处理；如果更注重MI、Q、SD、SF则不需分频。4 结语本文首次使用LScN模型进行图像融合实验，并使用多通工作方式作为迭代终止条
30、件，使用该模型中L项代替点火频率作为输出，解决了传统PcNN模型参数过多、参数和迭代次数难设置及融合效果差等问题，同时得出使用LscN模型进行图像融合时是否需要分频处理的结论。通过与五种算法进行对比得出，本文算法不论是从人眼还是客观评价指标均优于对比算法。目前基于PcNN模型进行图像融合的文献较多，而基于两种经典的简化模型IcM和scM的融合算法研究较少，特别是，截止目前尚未发现基于LScN模型进行图像融合的文献和相关研究成果。因此，接下来的工作就是对该模型进一步地深入研究，比如对激励和连接强度的改进等。参考文献(References)【l】 wANG z B，wANc s，zHu Y，et
31、a1ReView of image fusionbased on pulsecoupled neural network【J】 ArchiVes of computational Metlods in Engineering，2016，23(4)：659671【2】 UU Y，UU S P，WANG Z FA generalfhmeworkforimagefusio“b鹊ed on mul6一scaIe咖sfonn and sparse representationj】 Inf0卜mation Fusion201524：147一l“【3】 余二永，王润生，徐学文基于统计模型的遥感图像多分辨融合
32、方法【J】电子学报，2008，36(3)：490493(SHE E Y，wANG RS，XU X WMultiResolution fusion of remote sensing image basedon 8tadstical modelJ】Acta Electronica sinica，2008，36(3)：490一4931【4】 ZHAN K，SHI J H，WANG H B，et a1Computational mechallismsof puIsecoupled neuraI networks： a comprehensive review【J】 A卜chives of Compu
33、tational Methods in En矛neering，201 7，24(3)：573588【5】 薛寺中，周爱平，梁久祯基于小波变换的自适应脉冲耦合神经网络图像融合【J】计算机应用，2010，30(12)：32253228(XUE S Z，ZHOU A P，UANG J ZImage fusion based on wavelettmIlsf0册aIld adaptive PCNN【J】 Joumal of C0mputer AppIica-tions，2010，30(12)：32253228)【6】 邓翔字，马义德PCNN参数自适应设定及其模型的改进【J】电子学报，2012，40(5
34、)：955964(DENG x Y，MA Y DPcNNmodel automatic par蛐eters dete丌TIinatjon and its modined modelJ】Acta Electmnica siIIica，2012，40(5)：955964)【7】 Qu x B，YAN J w，xIA0 H z，et a1Image fusion algorithmbased on spa“al frequency-motivated pulSe coupled neural networksin nonsubs啪pled contourlet transf0珊domain【J】Ac
35、ta AutomaticaSinica，2008，34(12)：1508一1514(下转第3568页)万方数据3568 计算机应用 第37卷【3】 MLAO L D，QI H REndmember extr8ction hm higIly mixed datausing minimum volume constmined nonnegative ma硒x f如toriz毗ion【J】IEEE Tr帅明ctions on Geoscience蚰d Remote Se璐ing2007，45(3)：765777【4】 ZYMNls A，KIM s J，sl(AF J，et a1Hypenpec叫im
36、age u啪ix-ing vi8 a】temBling p叫ecled su蛔mdi朗ts f C】ACSSC 2007：P肛ceedin8s of the 2007 Confe陀nce Record of tlle For【yFirst Asilom盯Conference on Signals， SystemsCompute玛 PiscBtaway， NJ：IEEE2007：1164一1168【5】 JLA S，QLAN Y Tconstrained nonne鹳tive ma试x factorization forhype措pectral u叫dxing【J】 IEEE Trans舵tio
37、n8 on Geoscience明dRemote Sensing，2009，47(1)：16l173【6】 杜世强，石玉清，王维兰，等基于图正则化的半监督非负矩阵分解【J】计算机工程与应用，2012，48(36)：194200(Du s Q，SHI Y Q，WANG w L，et a1Gmph regIll耐舱dbased selllisuper_“sed nonnq弘tive ma喇x factori龃cion【J】 Compucer En垂neering鞠d Application8，2012，48(36)：194200)【7】 刘建军，吴泽彬，韦志辉，等基于空间相关性约束稀疏表示的高光谱
38、图像分类【J】电子与信息学报，2012，34(11)：26662671(UU J J，WU Z B，wEI Z H，et a1Spatial conl砒i叩cons乜ained印ar鸵repre鸵ntabon for hype玛pectml image cla黯ification【J】J咖止ofElectIDnicsI山n11ation khlogy，2012，34(11)：266626711【8】 CHEN B L，U M，WANG J X，et a1Di骶：ase gene identificaon byusing graph kemelsd Markov mdom fields【J】sc
39、ience chinaL讧e sciences，2014，57(11)：1054一1063【9】 姜小燕，孙福明，李豪杰基于图正则化和稀疏约束的半监督非负矩阵分解【J】计算机科学，2016，43(7)：7782(JIANG x Y，SUN F M，U H JSemisupervisedflrIegative瑚硒x船tofi越donb塔ed on gmph陀gul耐勰tion明d spar靶ness constraints【J】 Computer Science，2016，43(7)：7782)(上接第3557页)【8】王俭，聂仁灿，周冬明，等多目标粒子群优化PcNN参数的图像融合算法【J】中国
40、图象图形学报，2016，2l(10)：1298一1306(wANG Q，NIE R C，zH0u D M，et a1Image fIlsion a190rithmusing PCNN model p出锄ete璐of multiobjective p丑nicle gwa珊opti一面ljon【J】J叫n瑚0f Image and G强phjcs，2016，2l(10)：1298一13061【9】 薛寺中，周爱平，梁久祯基于小波变换的PcNN多传感器图像融合【J】计算机工程与应用，201l，47(35)：2lO一213(xuE sZ，ZHOU A P，UANG J ZFusion of mIllt
41、i一鸵n$0r imag姻based onwavelet tr蚰8fb唧柚d PcNN【J】computer EIlgineering帅d Ap出cations，20ll，47(35)：2lO一213)【lO】 EKBLAD U，KINSER J M，ATMER J，et a1The intersecting cortical model in image proee牺ing【J】Nucle盯In8trumen协Meth-ods in P11ysics Rese8rch，2004，525(12)：392396IlI zHAN K，ZHANG H J，MA Y D卜kw spiking cort
42、ical m伽lel forinv丽肌t textu托re晡eval卸d image processingIEEE 7rr帅锄tionsNeu糟l Networks，2009，20(12)：19801986f12】 KONc w WMuln-sensor j删lge fusj蚰based on NSsT d伽“n12cM【J】Electmnics kters，2013，49(13)：802803【13】 wANG N Y，MA Y D，ZHAN KSpiking co而cal model fbr mIll“-fhus image fusion【J】Neurocomputing，2014，130
43、(3)：4451【14】 wANG R，Wu Y，DING M Y，et a1MedicaI image fusion b啪d0n spiking eomcal model【C】Proceedin伊of tlle SPIE 8676，Medical I肌酒ng 2013：D酒t丑l PatlologyBellingIlm：SPIE，【10】DOBRuSHIN P LThe descption of a mndom field by means olconditional pmbabilities锄d conditions of its”gIIl耐ty【J】卟eoryof Prob&bility
44、蛐d its Appljcations，1968，13(2)：】97224【l l】 sPrrZER F Markov mdom fields and gibbs ensembles【J】American Mathematical Monthly，197l，78(2)：142一154f 12】 DENG H w， CLAuSI D A UnsuperVised image segmentationusing a simple MRF model witIl a new implementation scheme【J】Pattem Reco印ition，2004，37(12)：23232335【
45、13】 B10UCASDIAS J M，NASCIMEN7r0 J M PHype鸺pectral sub-spe identmcati【J1IEEE TraIl船ctions on&：0science蛐d Re-mote sensing，2008，46(8)：24352445【14】 袁博高光谱像元解混盲处理方法研究面向盲分解的相关性分析与应用【D1北京：中国科学院大学，2015：4548(YUAN BResearch佣blind proce鸽ing a190ritlms 0f hype聘pec-tral unIIIi】【ingcoHlation蛐alysis蛐d印plication 0f
46、blid deco-position(D】B叫ing：universify of Chine辩Ac8d绷y 0f sciences，2015：45481【15】 LEE D D，SEuNG H SA190rithms for nonnegative眦咄factori-zalion【J】Advces in Neuml Info丌nation Proce她System，200l，13(6)：556562【16】 吴金彪DN交替迭代法及其收敛性分析【J1数值计算与计算机应用，2002，23(2)：12l130(WU J BDN c0蚴utalivemethod柚d its c仰ve学nce【J】J叫
47、fnal鲫Numerical M如0dsand C鲫puter Applications，2002，23(2)： 12l130)17】 朱子明，祁新华基于Mor”I的闽南三角洲空间发展研究【J】经济地理，2009，29(12)：1977一1980(zHu z M，QI x HRe-鸵arch on s邮al development of Ininnan delta bl峙ing on MoranI【J】EcoIIIic G盱lphy，2009。29(12)：1977一1980)YUAN跏。b唧in 1982，PhD，Lect哦rHi8 reaIh interestsinclude hype玛p
48、ectml data processing，Intemet of tlings engineeTing2013：8676lO【15】 zHAN K，SHI J H，TENG J C，ea1“nking syn卵fjc compulatifor im唔e enhcement【J】NeurDco“puting，2017，238：l一12【16】【17】【18】ECKHORN R。REnBOECK H J，ARNIM，et a1A neural networl【for f朗tu陀linking via synchron叫s activiIy【J】C锄adi蛐J叫mal of Microbiology，1989，46(8)：759763WANG Z B， MA Y D， GU J Multifbcus image fusion usingPcNN【J】Pattem RecogIlition，20lO，43(6)：20032016王念一脉冲发放皮层模型图像融合技术研究【D】兰州：兰州大学，2014：2847(wANG N YS