CN106940895A - 一种应用于无线层析成像系统的降质函数的估计方法 - Google Patents
一种应用于无线层析成像系统的降质函数的估计方法 Download PDFInfo
- Publication number
- CN106940895A CN106940895A CN201710021570.2A CN201710021570A CN106940895A CN 106940895 A CN106940895 A CN 106940895A CN 201710021570 A CN201710021570 A CN 201710021570A CN 106940895 A CN106940895 A CN 106940895A
- Authority
- CN
- China
- Prior art keywords
- function
- matrix
- degrades
- convolution
- rti
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 59
- 238000003384 imaging method Methods 0.000 title claims abstract description 44
- 230000006870 function Effects 0.000 claims abstract description 96
- 239000011159 matrix material Substances 0.000 claims abstract description 92
- 230000000694 effects Effects 0.000 claims abstract description 26
- 230000008569 process Effects 0.000 claims abstract description 24
- 238000011084 recovery Methods 0.000 claims abstract description 11
- 230000009466 transformation Effects 0.000 claims abstract description 10
- 230000006866 deterioration Effects 0.000 claims abstract description 7
- 238000013507 mapping Methods 0.000 claims abstract description 4
- 230000000593 degrading effect Effects 0.000 claims description 10
- 238000004458 analytical method Methods 0.000 claims description 6
- 238000009826 distribution Methods 0.000 claims description 5
- 238000003325 tomography Methods 0.000 claims description 5
- 230000007704 transition Effects 0.000 claims description 5
- 238000006243 chemical reaction Methods 0.000 claims description 4
- 238000005192 partition Methods 0.000 claims description 3
- 238000012546 transfer Methods 0.000 claims description 3
- 238000013519 translation Methods 0.000 claims description 3
- 238000005562 fading Methods 0.000 claims 4
- 238000012545 processing Methods 0.000 abstract description 3
- 238000004891 communication Methods 0.000 description 7
- 238000004088 simulation Methods 0.000 description 7
- 238000005516 engineering process Methods 0.000 description 6
- 238000005259 measurement Methods 0.000 description 6
- 238000005457 optimization Methods 0.000 description 6
- 230000008859 change Effects 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 2
- 230000007423 decrease Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 230000004888 barrier function Effects 0.000 description 1
- 238000004587 chromatography analysis Methods 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 230000002452 interceptive effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
- G06T11/003—Reconstruction from projections, e.g. tomography
- G06T11/008—Specific post-processing after tomographic reconstruction, e.g. voxelisation, metal artifact correction
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S11/00—Systems for determining distance or velocity not using reflection or reradiation
- G01S11/02—Systems for determining distance or velocity not using reflection or reradiation using radio waves
- G01S11/06—Systems for determining distance or velocity not using reflection or reradiation using radio waves using intensity measurements
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10072—Tomographic images
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Image Processing (AREA)
- Image Analysis (AREA)
Abstract
本发明公开了一种应用于无线层析成像系统的降质函数的估计方法,属于无线网络中图像处理领域。令无线层析成像过程为线性的图像降质过程,结合约束最小二乘法准则,使用图像复原方法可以解决无线层析成像系统目标成像模糊、具有膨胀效应的问题,获得较清晰的“原始”图像。本发明提出将函数卷积过程转换为矩阵相乘的形式,推导出降质函数的卷积变换矩阵与无线层析成像系统线性解的关系;然后由矩阵变换理论得到降质函数与卷积变换矩阵之间的映射关系,并通过建立混合高斯模型估计出降质函数。本发明应用于无线层析成像系统后,有效地减弱目标成像的膨胀效应,提高无线层析成像系统的成像精度,可应用于实际环境中。
Description
技术领域
本发明涉及一种基于无线层析成像系统的降质函数的估计方法,应用于无线层析成像系统的图像复原过程中,属于无线网络中的图像处理领域。
背景技术
无线层析成像(RTI)是一种用来对无线传感器网络区域内的目标进行定位的无源定位技术,其利用无线节点间的通信链路被目标遮挡引起的信号强度衰减来反演出区域内各个位置处的信号衰落强度,得到区域内的信号衰减图,进而得到目标所有的位置信息。近年来,随着各个领域对位置服务需求的增加,无线定位技术相关的研究越来越活跃。RTI技术作为一种新兴的无源定位技术,引起了众多领域的兴趣,并已经取得了大量研究成果。最早是由美国犹他大学的Wilson和Patwari提出的基于阴影衰减的RTI方法(SRTI),SRTI是利用接收信号强度(RSS)变化在无线网络中成像的方法。SRTI方法假设无线链路被目标遮挡会引起很大的阴影衰减,没被障碍物遮挡的链路则会保持稳定。
这种经典的基于阴影的RTI成像技术过程为:
将S个无线节点等距离部署在监测区域周围,这些节点具有相同的物理结构,工作在相同频段并支持相同的通信协议。每个节点坐标已知,为 (xs,ys),s∈{1,2,...,S},这些无线节点共构成L=S(S-1)/2条无线链路,每个节点按照预先设定的协议和时序发送数据,并且接收及测量其他节点发的无线信号的RSS 值。
测量时先测量监控区域内没有目标时,每条链路的RSS值rl,其中l是链路编号;然后测量目标出现在监测区域内时每条链路在离散时刻t的RSS值rl,t,从而得到第l条链路在t时刻的RSS变化量为:△rl,t=rl-rl,t,l=1,2,...,L.。将监测区域分割成N个小的区域,每个小的区域称为像素,每个像素用△xj,j∈{1,2,...,N}表示,则每条链路的RSS变化量看成每个像素的RSS变化量的加权和,那么第l条链路在t时刻的RSS变化量就可以用公式表示为:其中△xj,t是在t时刻发生在像素j内的RSS衰减值,nl,t是时刻t时链路l上的测量噪声,wlj是链路l里的像素j的权重值。
所有链路的RSS变化量的表达式写成矩阵的形式:
r=Wx+n, (1)
r和n都是L×1向量,分别表示所由链路的RSS变化量和噪声;x是N×1向量,表示衰减图像;W是L×N维的权重矩阵;式(1)为RTI系统对无线网络中目标成像的线性方程表达式。
利用正则化方法求解线性方程,得到衰减图像向量的线性解为:
Π是线性转移矩阵,DX和DY分别为水平方向和垂直方向的差分算子,式(2) 即为RTI系统的线性方程解,也是基于阴影衰减的RTI系统对目标的成像结果。
但是在上述RTI方法中,由于无线通信链路中存在相互干扰和测量链路数目有限的问题,导致得到的目标衰落图像的精度不高,成像后的图像是目标的膨胀的图像区域,这样的成像结果不能提供足够的细节来识别想要的目标。这个膨胀效应是由无线通信链路有限导致的,然而当增加链路时,无线传感器节点和扫描所有通信链路的时间将会增加,传感器节点之间的干扰也会增强,这些也会降低成像的精度。所以考虑采用保持链路数目不变的方法提高成像效果,图像复原方法就是这样一种方法。由于成像出的图像是“膨胀”的,也可以说是被“降质”的,图像复原方法就是对此图像进行优化,减弱它的“膨胀效应”,从而获得更清晰的图像。
图像复原技术是试图利用图像降质过程中的先验知识将“降质”图像复原为“原始”清晰图像的过程。具体过程为:将图像降质过程建模为一个线性过程,即降质图像是由“原始”图像与降质函数卷积再叠加噪声得到的,而图像复原就是图像降质的逆过程。所以问题的关键就是如何得到降质函数。
发明内容
为了解决上述问题,本发明提出了一种应用矩阵理论和混合高斯模型来估计无线层析成像系统的降质函数的方法。估计降质函数的具体方法为:首先推导获得降质函数的卷积变换矩阵与RTI系统线性解的关系,计算出卷积变换矩阵;然后由矩阵变换理论得到降质函数与其卷积变换矩阵之间的数量关系,再通过建立混合高斯模型估计出降质函数。获得降质函数后,用降质函数结合约束最小二乘法准则进行图像目标优化处理,获得“原始”成像图形。实验证明,使用本发明求出降质函数再进行图像复原,复原后的图像精度比不优化的RTI 成像精度提高很多。
本发明所述的应用于无线层析成像(RTI)系统的降质函数的估计方法,具体包含如下步骤:
步骤一:建模RTI系统降质函数;
由于基于阴影衰减的RTI方法得到的RTI图像具有膨胀效应,所以使用图像复原方法减弱膨胀效应;在图像复原中,将图像降质过程看成一个线性模型,降质后的图像表示为:
g(x,y)=h(x,y)*f(x,y)+n(x,y), (4)
h(x,y)是空间域的降质函数,*表示卷积;f(x,y)表示“原始”的没有膨胀效应的图像;g(x,y)是降质后的实际成像出来的RTI图像,n(x,y)是噪声函数;
步骤二:获得降质函数和RTI系统线性解的关系,得到卷积变换矩阵;
根据卷积定义,将(4)式写成向量矩阵相乘的形式:
g=Hf+η, (5)
f是N×1的矩阵,表示“原始”输入图像;g和η也是N×1矩阵;H是大小为 N×N的卷积变换矩阵;
将(1)式代入(2)式,得到:
对比式(6)和式(5),由于g和都表示的是成像出来的RTI图像,所以假设噪声相同,得到如下方程:
Hf=ΠWx, (7)
且f和x表示的是相同的原始图像向量,故获得卷积变换矩阵和RTI系统线性解的关系如下:
H=ΠW, (8)
即卷积变换矩阵由RTI系统的线性变换矩阵和权重矩阵相乘得到;
步骤三:根据矩阵理论由降质函数获得卷积变换矩阵;
卷积转移矩阵H的元素由二维圆周卷积定义为:
f(m,n)表示原始图像在(m,n)处的原始像素值,表示在降质后图像在坐标(x,y)处的像素值;假设降质函数h大小为(2K+1)×(2K+1),它的元素即为:
由于二维卷积运算的实质是将卷积模板翻转180度,然后将该卷积模板依次从上到下、从左到右滑动,计算在模板与原始图像产生交集元素的乘积和,作为卷积以后的数值;因此可以得到翻转后的点扩散函数h-为:
再计算h-与图像产生交集元素的乘积和就可得到卷积的结果;此时的运算方式仍然为两个二维矩阵的平移,计算乘积和;下面推导如何由h-构成H;
由卷积定义可知,图像卷积结果的第一个像素点为矩阵h-中的区域与图像区域f(m,n)对应元素的乘积和;区域为:
其中为的列向量,因此H的第一行数据为:
同理可以求得H的第二行数据为:
其中它与相比多了一个元素h1,i,由此可以分析出矩阵H的前M行前M列的数据为:
同理分析矩阵H的前M行,(M+1):2M列的数据为:
由此类推可以定义由矩阵h生成的矩阵Hi(-K≤i≤K)为:
结合上述分析,可以将卷积转换矩阵H简化为:
到此就实现了二维卷积运算与矩阵相乘运算之间的转换,得到卷积变换矩阵H;且知H中的元素都是降质函数h(x,y)的特定元素分布在特定位置;
步骤四:用混合高斯模型估计出降质函数的元素,得到降质函数h(x,y);
由于矩阵H中的元素来自降质函数h(x,y)中的特定元素,将矩阵H中的特定元素的值看作降质函数h(x,y)的特定元素的样本过程;然后将降质函数中的每个元素建模成一个混合Q高斯模型:
其中hi,j是降质函数h(x,y)在(x,y)处的元素值,P(hi,j)是hi,j的概率密度函数,Q是高斯分布的个数,wq是第q个高斯分布的权重并都为1;g(hi,j,μq,σq)表示高斯概率密度函数,μq和σq为均值和标准差;
由于降质函数的每个元素的很多样本都处在矩阵H中的很多不同位置上,故引入一个学习算法来不断更新模型参数;
首先根据上式(18)和(19)的映射关系将H中的元素分配到相应的数据组hi,j,m(1≤m≤M)中;然后检验第一个数据组中的每个元素值,如果其中的元素值和已有的Q高斯分布不匹配,那么就将新加入的样本数据和原有的数据一起作为新的数据集合,并用新数据集合的均值作为新的高斯分布的均值,用新数据集合的无偏样本方差作为新的高斯分布的方差;若元素值和其中一个Q高斯分布相匹配,则将第m个元素值的Q分布的先验权重更新为:
wq,m=(1-β)wq,m-1+β(Mq,m), (21)
β是一个学习速度参数,决定分布参数收敛的快慢;Mq,m匹配时为1,不匹配为0;
若不匹配,高斯分布的均值和标准差保持不变;若匹配,更新为:
μm=(1-ρ)μm-1+ρhi,j,m, (22)
其中ρ=βg(hi,j,m|μq,σq)为第二个学习参数,描述数据与估计的模型间匹配的程度;处理完一个数据组中的元素值之后,计算出相应的降质函数的元素hi,j:
以此类推估计出h(x,y)中的全部元素,得到RTI系统的降质函数;
对比现有技术,本发明有益效果在于能通过矩阵理论和混合高斯模型估计出降质函数,降质函数结合约束最小二乘法优化RTI图像,能有效降低成像中目标的膨胀效应,提高RTI系统的定位精度,且能广泛应用于很多RTI系统中。
附图说明
图1:基于无线层析成像的降质函数算法流程图;
图2:实验过程中节点布置;
图3:实验过程复原前的RTI成像效果图;
图4:实验过程复原后的RTI成像效果图;
图5:模拟过程中节点布置;
图6:模拟过程中模拟的三种目标形状;
图7:模拟过程图像复原前的RTI成像效果图;
图8:模拟过程使用降质函数复原后的RTI成像效果图;
具体实施方式
下面将结合附图和实施例对本发明加以详细说明,同时也叙述了本发明技术方案解决的技术问题及有益效果。需要指出的是,所描述的实施例仅旨在便于对本发明的理解,而对其不起任何限定作用。
本发明所述的应用于无线层析成像(RTI)系统的降质函数的估计方法流程图如图1所示,具体包含如下实施步骤:
步骤一:配置节点,由基于阴影的RTI方法得到RTI图像;
本实施例中节点放置如图2所示,监测区域位于二维xoy平面;将20个具有相同的物理结构,工作在相同频段并支持相同的通信协议的无线节点等距离部署在9.5×9.5m2的监测区域周围,相邻节点之间相距1.9m;
每个节点按照预先设定的协议和时序发送数据,并且接收及测量其他节点发的无线信号的RSS值,即:在t时刻,编号为s的节点发送数据,其他节点接收数据并测量接收信号强度;在下个时刻,编号为s+1的节点发送数据,其他节点接收数据并测量接收信号强度;
测量时先测量监控区域内没有目标时,每条链路的RSS值rl,其中l是链路编号;然后测量目标出现在监测区域内时每条链路在离散时刻t的RSS值rl,t,从而得到第l条链路在t时刻的RSS变化量为:△rl,t=rl-rl,t,l=1,2,...,L.;将监测区域分割成N个小的区域,每个小的区域称为像素,每个像素用△xj,j∈{1,2,...,N}表示,每条链路的RSS变化量看成每个像素的RSS变化量的加权和,那么第l条链路在t时刻的RSS变化量就可以用公式表示为:其中△xj,t是在t时刻发生在像素j内的RSS衰减值,nl,t是时刻t时链路l上的测量噪声,wlj是链路l里的像素j的权重值;
所有链路的RSS变化量的表达式为:
r=Wx+n, (1)
r和n都是L×1向量,分别表示所由链路的RSS变化量和噪声;x是N×1向量,表示衰减图像;W是L×N维的权重矩阵;
利用正则化方法求解线性方程,得到衰减图像向量的线性解为:
Π是线性转移矩阵,DX和DY分别为水平方向和垂直方向的差分算子;
分别令单个目标和两个目标在监测区域中运动,运用上述方法对运动目标成像,得到单目标和双目标的RTI图像分别如图3的(a)图和(b)图所示;
步骤二:建模RTI系统降质函数;
由于基于阴影衰减的RTI方法得到的RTI图像具有膨胀效应,所以使用图像复原方法减弱膨胀效应;在图像复原中,将图像降质过程看成一个线性模型,降质后的图像表示为:
g(x,y)=h(x,y)*f(x,y)+n(x,y), (4)
h(x,y)是空间域的降质函数,*表示卷积;f(x,y)表示“原始”的没有膨胀效应的图像;g(x,y)是降质后的实际成像出来的RTI图像,n(x,y)是噪声函数;
步骤三:获得降质函数和RTI系统线性解的关系,得到卷积变换矩阵;
根据卷积定义,将(4)式写成向量矩阵相乘的形式:
g=Hf+η, (5)
f是N×1的矩阵,表示“原始”输入图像;g和η也是N×1矩阵;H是大小为N×N的卷积变换矩阵;
将(1)式代入(2)式,得到:
对比式(6)和式(5),由于g和都表示的是成像出来的RTI图像,所以假设噪声相同,得到如下方程:
Hf=ΠWx, (7)
且f和x表示的是相同的原始图像向量,故获得卷积变换矩阵和RTI系统线性解的关系如下:
H=ΠW, (8)
即卷积变换矩阵由RTI系统的线性变换矩阵和权重矩阵相乘得到;
步骤四:根据矩阵理论由降质函数获得卷积变换矩阵;
卷积转移矩阵H的元素由二维圆周卷积定义为:
f(m,n)表示原始图像在(m,n)处的原始像素值,表示在降质后图像在坐标(x,y)处的像素值;假设降质函数h大小为(2K+1)×(2K+1),它的元素即为:
由于二维卷积运算的实质是将卷积模板翻转180度,然后将该卷积模板依次从上到下、从左到右滑动,计算在模板与原始图像产生交集元素的乘积和,作为卷积以后的数值;因此可以得到翻转后的点扩散函数h-为:
再计算h-与图像产生交集元素的乘积和就可得到卷积的结果;此时的运算方式仍然为两个二维矩阵的平移,计算乘积和;下面推导如何由h-构成H;
由卷积定义可知,图像卷积结果的第一个像素点为矩阵h-中的区域与图像区域f(m,n)对应元素的乘积和;区域为:
其中为的列向量,因此H的第一行数据为:
同理可以求得H的第二行数据为:
其中它与相比多了一个元素h1,i,由此可以分析出矩阵H的前M行前M列的数据为:
同理分析矩阵H的前M行,(M+1):2M列的数据为:
由此类推可以定义由矩阵h生成的矩阵Hi(-K≤i≤K)为:
结合上述分析,可以将卷积转换矩阵H简化为:
到此就实现了二维卷积运算与矩阵相乘运算之间的转换,得到卷积变换矩阵H;且知H中的元素都是降质函数h(x,y)的特定元素分布在特定位置;
步骤五:用混合高斯模型估计出降质函数的元素,得到降质函数h(x,y);
由于矩阵H中的元素来自降质函数h(x,y)中的特定元素,将矩阵H中的特定元素的值看作降质函数h(x,y)的特定元素的样本过程;然后将降质函数中的每个元素建模成一个混合Q高斯模型:
其中hi,j是降质函数h(x,y)在(x,y)处的元素值,P(hi,j)是hi,j的概率密度函数,Q是高斯分布的个数,wq是第q个高斯分布的权重并都为1;g(hi,j,μq,σq)表示高斯概率密度函数,μq和σq为均值和标准差;
由于降质函数的每个元素的很多样本都处在矩阵H中的很多不同位置上,故引入一个学习算法来不断更新模型参数;
首先根据上式(18)和(19)的映射关系将H中的元素分配到相应的数据组hi,j,m(1≤m≤M)中;然后检验第一个数据组中的每个元素值,如果其中的元素值和已有的Q高斯分布不匹配,那么就将新加入的样本数据和原有的数据一起作为新的数据集合,并用新数据集合的均值作为新的高斯分布的均值,用新数据集合的无偏样本方差作为新的高斯分布的方差;若元素值和其中一个Q高斯分布相匹配,则将第m个元素值的Q分布的先验权重更新为:
wq,m=(1-β)wq,m-1+β(Mq,m), (21)
β是一个学习速度参数,决定分布参数收敛的快慢;Mq,m匹配时为1,不匹配为0;
若不匹配,高斯分布的均值和标准差保持不变;若匹配,更新为:
μm=(1-ρ)μm-1+ρhi,j,m, (22)
其中ρ=βg(hi,j,m|μq,σq)为第二个学习参数,描述数据与估计的模型间匹配的程度;处理完一个数据组中的元素值之后,计算出相应的降质函数的元素hi,j:
以此类推估计出h(x,y)中的全部元素,得到RTI系统的降质函数;
步骤六:使用步骤五中获得的降质函数结合约束最小二乘准则进行图像复原,得到原始RTI图像;
由步骤二可知降质后的图像函数可表示为(4),且空间域的卷积相当于频域相乘,则频域表达式为:
G(u,v)=H(u,v)F(u,v)+N(u,v), (25)
以上函数均为(4)中对应函数的傅里叶变换;用约束最小二乘法准则强制原始输入图像f(x,y)变得平滑,约束最小二乘准则表达式为:
其中第二项l(x,y)*f(x,y)为估计图像进行滤波后的图像,常使用的是拉普拉斯高通滤波器,可使得到的图像更为平滑且局部区域具有高相关性;γ是一个可调参数,用来控制恢复图像的平滑程度,值越大图像越平滑;将上述时域表达式转换到频域得到:
以上函数均为对应函数的傅里叶变换;对上式求导并令导数等于零,得到恢复估计图像的频域解为:
再通过傅里叶逆变换得到估计出的原始图像为:
按上述步骤复原后的RTI图像如图4所示,(a)和(b)分别是单目标和双目标的复原后的RTI图像,分别与图3中的对应图像对比,可以看出复原后的 RTI图像减小了目标成像的膨胀效应,成像图像更清晰。
利用均方误差MSE指数来评价两种图像,结果显示当人类半径为0.4m时,计算不同位置的平均均方误差MSE值,RTI成像时值为0.0135,用降质函数复原后的RTI图像的平均MSE值为0.0101。说明复原后的图像的定位精度更高,证明了本发明提出的应用于RTI图像复原的降质函数是有效的。
本实施例中还包含一个模拟过程。模拟过程中节点放置如图5所示,监测区域位于二维xoy平面;将20个具有相同的物理结构,工作在相同频段并支持相同的通信协议的无线节点等距离部署在8×8m2的监测区域周围,相邻节点之间相距1.6m。
模拟过程包含三种类型的模拟目标形状,如图6所示:(a)是条形,(b) 空心矩形,(c)L型。对这三种形状的目标进行RTI成像,成像效果图如图7 所示,目标具有膨胀效应,效果较差,尤其目标边缘处细节都很模糊。然后应用本发明所述的方法求出降质函数,再将降质函数应用于图像复原中,得到的“原始”的RTI图像如图8所示。我们可以看出优化后的“原始”的RTI图像在定位精度上表现更好,目标更清晰,有效地削减了目标的膨胀效应,目标的成像区域更接近于目标的中心线。
引入边缘保持指数EPI和结构相似指数SSIM来评估优化前后的RTI成像效果,发现与实际目标相比,优化的RTI成像后的边缘保持指数EPI和结构相似指数SSIM都由大幅度提高,说明优化后的成像更接近目标的真实形态。证明了本发明提出的应用于RTI图像复原的降质函数是有效的。
以上所述,仅为本发明的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉该技术的人在本发明所揭露的技术范围内,可理解想到的变换和替换,都应涵盖在本发明的包含范围之内,因此,本发明的保护范围应该以权利要求书的保护范围为准。
Claims (7)
1.应用于无线层析成像(RTI)系统的降质函数的估计方法,其特征在于,包含如下步骤:
步骤一:建模RTI系统降质函数;
由于基于阴影衰减的RTI方法得到的RTI图像具有膨胀效应,所以使用图像复原方法减弱膨胀效应;在图像复原中,将图像降质过程看成一个线性模型,降质后的图像表示为:
g(x,y)=h(x,y)*f(x,y)+n(x,y), (4)
h(x,y)是空间域的降质函数,*表示卷积;f(x,y)表示“原始”的没有膨胀效应的图像;g(x,y)是降质后的实际成像出来的RTI图像,n(x,y)是噪声函数;
步骤二:获得降质函数和RTI系统线性解的关系,得到卷积变换矩阵;
根据卷积定义,将(4)式写成向量矩阵相乘的形式:
g=Hf+η, (5)
f是N×1的矩阵,表示“原始”输入图像;g和η也是N×1矩阵;H是大小为N×N的卷积变换矩阵;
将(1)式代入(2)式,得到:
对比式(6)和式(5),由于g和都表示的是成像出来的RTI图像,所以假设噪声相同,得到如下方程:
Hf=ΠWx, (7)
且f和x表示的是相同的原始图像向量,故获得卷积变换矩阵和RTI系统线性解的关系如下:
H=ΠW, (8)
即卷积变换矩阵由RTI系统的线性变换矩阵和权重矩阵相乘得到;
步骤三:根据矩阵理论由降质函数获得卷积变换矩阵;
卷积转移矩阵H的元素由二维圆周卷积定义为:
f(m,n)表示原始图像在(m,n)处的原始像素值,表示在降质后图像在坐标(x,y)处的像素值;假设降质函数h大小为(2K+1)×(2K+1),它的元素即为:
由于二维卷积运算的实质是将卷积模板翻转180度,然后将该卷积模板依次从上到下、从左到右滑动,计算在模板与原始图像产生交集元素的乘积和,作为卷积以后的数值;因此可以得到翻转后的点扩散函数h-为:
再计算h-与图像产生交集元素的乘积和就可得到卷积的结果;此时的运算方式仍然为两个二维矩阵的平移,计算乘积和;下面推导如何由h-构成H;
由卷积定义可知,图像卷积结果的第一个像素点为矩阵h-中的区域与图像区域f(m,n)对应元素的乘积和;区域为:
其中为的列向量,因此H的第一行数据为:
同理可以求得H的第二行数据为:
其中它与相比多了一个元素h1,i,由此可以分析出矩阵H的前M行前M列的数据为:
同理分析矩阵H的前M行,(M+1):2M列的数据为:
由此类推可以定义由矩阵h生成的矩阵Hi(-K≤i≤K)为:
结合上述分析,可以将卷积转换矩阵H简化为:
到此就实现了二维卷积运算与矩阵相乘运算之间的转换,得到卷积变换矩阵H;且知H中的元素都是降质函数h(x,y)的特定元素分布在特定位置;
步骤四:用混合高斯模型估计出降质函数的元素,得到降质函数h(x,y);
由于矩阵H中的元素来自降质函数h(x,y)中的特定元素,将矩阵H中的特定元素的值看作降质函数h(x,y)的特定元素的样本过程;然后将降质函数中的每个元素建模成一个混合Q高斯模型:
其中hi,j是降质函数h(x,y)在(x,y)处的元素值,P(hi,j)是hi,j的概率密度函数,Q是高斯分布的个数,wq是第q个高斯分布的权重并都为1;g(hi,j,μq,σq)表示高斯概率密度函数,μq和σq为均值和标准差;
由于降质函数的每个元素的很多样本都处在矩阵H中的很多不同位置上,故引入一个学习算法来不断更新模型参数;
首先根据上式(18)和(19)的映射关系将H中的元素分配到相应的数据组hi,j,m(1≤m≤M)中;然后检验第一个数据组中的每个元素值,如果其中的元素值和已有的Q高斯分布不匹配,那么就将新加入的样本数据和原有的数据一起作为新的数据集合,并用新数据集合的均值作为新的高斯分布的均值,用新数据集合的无偏样本方差作为新的高斯分布的方差;若元素值和其中一个Q高斯分布相匹配,则将第m个元素值的Q分布的先验权重更新为:
β是一个学习速度参数,决定分布参数收敛的快慢;Mq,m匹配时为1,不匹配为0;
若不匹配,高斯分布的均值和标准差保持不变;若匹配,更新为:
μm=(1-ρ)μm-1+ρhi,j,m, (22)
其中ρ=βg(hi,j,m|μq,σq)为第二个学习参数,描述数据与估计的模型间匹配的程度;处理完一个数据组中的元素值之后,计算出相应的降质函数的元素hi,j:
以此类推估计出h(x,y)中的全部元素,得到RTI系统的降质函数。
2.根据权利要求1所述一种无线层析成像系统的降质函数的估计方法,其特征在于,步骤二中,将RTI成像系统表达式变换成矩阵向量相乘的形式。
3.根据权利要求1所述一种无线层析成像系统的降质函数的估计方法,其特征在于,步骤二中,降质函数的卷积变换矩阵和RTI线性方程解的关系为卷积变换函数由RTI线性方程解矩阵和权重矩阵相乘得出。
4.根据权利要求1所述一种基于阴影衰落的无线探测网络节点定位方法,其特征在于,步骤三中,应用矩阵变换理论中的二维卷积运算与矩阵相乘运算之间的转换关系,推导得出降质函数和卷积变换矩阵关系的运算过程。
5.根据权利要求1所述一种基于阴影衰落的无线探测网络节点定位方法,其特征在于,步骤三中,卷积变换矩阵由很多分块子矩阵构成,详见公式(23)。
6.根据权利要求1所述一种基于阴影衰落的无线探测网络节点定位方法,其特征在于,步骤三中,卷积变换矩阵的子矩阵中的元素是降质函数的元素在特定位置的分布,详见公式(22)。
7.根据权利要求1所述一种基于阴影衰落的无线探测网络节点定位方法,其特征在于,步骤四中,应用混合高斯模型求解降质函数矩阵中元素的估计方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710021570.2A CN106940895B (zh) | 2017-01-12 | 2017-01-12 | 一种应用于无线层析成像系统的降质函数的估计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710021570.2A CN106940895B (zh) | 2017-01-12 | 2017-01-12 | 一种应用于无线层析成像系统的降质函数的估计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106940895A true CN106940895A (zh) | 2017-07-11 |
CN106940895B CN106940895B (zh) | 2021-06-29 |
Family
ID=59469446
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710021570.2A Active CN106940895B (zh) | 2017-01-12 | 2017-01-12 | 一种应用于无线层析成像系统的降质函数的估计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106940895B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108761391A (zh) * | 2018-05-29 | 2018-11-06 | 南京信息工程大学 | 一种模型类无设备目标定位方法 |
CN110678899A (zh) * | 2017-06-02 | 2020-01-10 | 三菱电机株式会社 | 成像系统和图像处理方法 |
CN112818514A (zh) * | 2021-01-05 | 2021-05-18 | 河南工业大学 | 无线层析成像方法及结合水平距离衰减的椭圆权重模型 |
CN113610221A (zh) * | 2021-06-29 | 2021-11-05 | 西安电子科技大学 | 一种基于fpga的可变膨胀卷积运算硬件系统 |
US11540798B2 (en) | 2019-08-30 | 2023-01-03 | The Research Foundation For The State University Of New York | Dilated convolutional neural network system and method for positron emission tomography (PET) image denoising |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20130272490A1 (en) * | 2010-12-15 | 2013-10-17 | Fujifilm Corporation | Radio tomographic image generation method and device |
CN103559699A (zh) * | 2013-11-18 | 2014-02-05 | 首都师范大学 | 一种基于投影估计的多能谱ct图像重建方法 |
CN104008539A (zh) * | 2014-05-29 | 2014-08-27 | 西安理工大学 | 基于多尺度几何分析的图像超分辨率重建方法 |
US20160041949A1 (en) * | 2014-08-06 | 2016-02-11 | International Business Machines Corporation | Dynamic highlighting of repetitions in electronic documents |
CN105719251A (zh) * | 2016-01-19 | 2016-06-29 | 浙江大学 | 一种用于大像移线性模糊的压缩降质图像复原方法 |
CN105913403A (zh) * | 2016-06-27 | 2016-08-31 | 乐视控股(北京)有限公司 | 图像降噪方法及系统 |
-
2017
- 2017-01-12 CN CN201710021570.2A patent/CN106940895B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20130272490A1 (en) * | 2010-12-15 | 2013-10-17 | Fujifilm Corporation | Radio tomographic image generation method and device |
CN103559699A (zh) * | 2013-11-18 | 2014-02-05 | 首都师范大学 | 一种基于投影估计的多能谱ct图像重建方法 |
CN104008539A (zh) * | 2014-05-29 | 2014-08-27 | 西安理工大学 | 基于多尺度几何分析的图像超分辨率重建方法 |
US20160041949A1 (en) * | 2014-08-06 | 2016-02-11 | International Business Machines Corporation | Dynamic highlighting of repetitions in electronic documents |
CN105719251A (zh) * | 2016-01-19 | 2016-06-29 | 浙江大学 | 一种用于大像移线性模糊的压缩降质图像复原方法 |
CN105913403A (zh) * | 2016-06-27 | 2016-08-31 | 乐视控股(北京)有限公司 | 图像降噪方法及系统 |
Non-Patent Citations (1)
Title |
---|
CHENG SUN: "Super-Resolution Reconstruction of Radio Tomographic Image", 《2016 IEEE 83RD VEHICULAR TECHNOLOGY CONFERENCE》 * |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110678899A (zh) * | 2017-06-02 | 2020-01-10 | 三菱电机株式会社 | 成像系统和图像处理方法 |
CN110678899B (zh) * | 2017-06-02 | 2023-04-28 | 三菱电机株式会社 | 成像系统和图像处理方法 |
CN108761391A (zh) * | 2018-05-29 | 2018-11-06 | 南京信息工程大学 | 一种模型类无设备目标定位方法 |
CN108761391B (zh) * | 2018-05-29 | 2022-04-01 | 南京信息工程大学 | 一种模型类无设备目标定位方法 |
US11540798B2 (en) | 2019-08-30 | 2023-01-03 | The Research Foundation For The State University Of New York | Dilated convolutional neural network system and method for positron emission tomography (PET) image denoising |
CN112818514A (zh) * | 2021-01-05 | 2021-05-18 | 河南工业大学 | 无线层析成像方法及结合水平距离衰减的椭圆权重模型 |
CN113610221A (zh) * | 2021-06-29 | 2021-11-05 | 西安电子科技大学 | 一种基于fpga的可变膨胀卷积运算硬件系统 |
CN113610221B (zh) * | 2021-06-29 | 2024-02-13 | 西安电子科技大学 | 一种基于fpga的可变膨胀卷积运算硬件系统 |
Also Published As
Publication number | Publication date |
---|---|
CN106940895B (zh) | 2021-06-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106940895A (zh) | 一种应用于无线层析成像系统的降质函数的估计方法 | |
Molina et al. | Bootstrap for estimating the MSE of the Spatial EBLUP | |
CN108520197A (zh) | 一种遥感图像目标检测方法及装置 | |
US20170299681A1 (en) | Sparse reconstruction strategy for multi-level sampled mri | |
Pigoli et al. | Wavelets in functional data analysis: estimation of multidimensional curves and their derivatives | |
Vallejos et al. | Effective sample size of spatial process models | |
CN108510013A (zh) | 基于低秩核心矩阵的改进稳健张量主成分分析方法 | |
CN105120517A (zh) | 基于多维尺度分析的室内wlan信号平面图构建与定位方法 | |
Matsuo et al. | Nonstationary covariance modeling for incomplete data: Monte Carlo EM approach | |
Hassani et al. | Deflation strategies for multi-block principal component analysis revisited | |
Wang et al. | Efficient dimension reduction for high-dimensional matrix-valued data | |
Bolin et al. | Latent Gaussian random field mixture models | |
Dimeglio et al. | A robust algorithm for template curve estimation based on manifold embedding | |
CN109581282B (zh) | 基于贝叶斯半监督深度学习的室内定位方法 | |
CN112991394B (zh) | 基于三次样条插值和马尔科夫链的kcf目标跟踪方法 | |
Arbia | Pairwise likelihood inference for spatial regressions estimated on very large datasets | |
CN110120066A (zh) | 面向监视系统的鲁棒群目标跟踪方法及跟踪系统 | |
Moran et al. | Two dimensional inverse imaging (2DII) of current sources in magnetoencephalography | |
CN114070262B (zh) | 附加扰动的集成混合集合Kalman滤波天气预报同化方法及装置 | |
Aghajani et al. | Intensity based image registration by minimizing the complexity of weighted subtraction under illumination changes | |
Guhaniyogi et al. | Large multi-scale spatial modeling using tree shrinkage priors | |
CN110780604B (zh) | 一种基于空时平滑性和时间相关性的空时信号恢复方法 | |
DeSarbo et al. | Constrained canonical correlation | |
Aleotti et al. | Fractional graph Laplacian for image reconstruction | |
Mutanen et al. | Tree height estimates in boreal forest using Gaussian process regression |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |