CN111524167B - 一种移动放射源定位校正与优化方法 - Google Patents
一种移动放射源定位校正与优化方法 Download PDFInfo
- Publication number
- CN111524167B CN111524167B CN202010213464.6A CN202010213464A CN111524167B CN 111524167 B CN111524167 B CN 111524167B CN 202010213464 A CN202010213464 A CN 202010213464A CN 111524167 B CN111524167 B CN 111524167B
- Authority
- CN
- China
- Prior art keywords
- time
- radiation
- radioactive source
- covariance
- source
- 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.)
- Active
Links
- 230000002285 radioactive effect Effects 0.000 title claims abstract description 81
- 238000000034 method Methods 0.000 title claims abstract description 68
- 238000005457 optimization Methods 0.000 title claims abstract description 11
- 230000005855 radiation Effects 0.000 claims abstract description 84
- 238000001914 filtration Methods 0.000 claims abstract description 40
- 238000009499 grossing Methods 0.000 claims abstract description 21
- 230000003287 optical effect Effects 0.000 claims abstract description 14
- 238000001514 detection method Methods 0.000 claims abstract description 11
- 238000005259 measurement Methods 0.000 claims description 30
- 239000011159 matrix material Substances 0.000 claims description 9
- 230000007704 transition Effects 0.000 claims description 7
- 230000000875 corresponding effect Effects 0.000 claims description 4
- 230000002596 correlated effect Effects 0.000 claims description 3
- 238000010586 diagram Methods 0.000 description 10
- 238000003384 imaging method Methods 0.000 description 6
- 230000000694 effects Effects 0.000 description 5
- 230000004927 fusion Effects 0.000 description 4
- 230000000052 comparative effect Effects 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 2
- 230000001133 acceleration Effects 0.000 description 1
- 238000009825 accumulation Methods 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 230000008092 positive effect Effects 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000017105 transposition Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/20—Analysis of motion
- G06T7/277—Analysis of motion involving stochastic approaches, e.g. using Kalman filters
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/22—Matching criteria, e.g. proximity measures
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/70—Determining position or orientation of objects or cameras
- G06T7/77—Determining position or orientation of objects or cameras using statistical methods
-
- 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/20—Special algorithmic details
- G06T2207/20076—Probabilistic image processing
-
- 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/20—Special algorithmic details
- G06T2207/20212—Image combination
- G06T2207/20221—Image fusion; Image merging
-
- 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/30—Subject of image; Context of image processing
- G06T2207/30241—Trajectory
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- General Physics & Mathematics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Life Sciences & Earth Sciences (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Evolutionary Biology (AREA)
- Data Mining & Analysis (AREA)
- Multimedia (AREA)
- Probability & Statistics with Applications (AREA)
- Artificial Intelligence (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Nuclear Medicine (AREA)
Abstract
本发明公开了一种移动放射源定位校正与优化方法,其步骤包括:1)获取放射源移动过程中的辐射图像,将每一所述辐射图像中像素最大值位置(x,y)作为放射源在辐射图像中的中心位置;2)采用卡尔曼滤波算法,从x、y方向分别对各辐射图像中的放射源位置进行校正和优化;3)将步骤2)校正和优化后的t时刻采集的辐射图像中的放射源位置与t时刻辐射图像对应的光学图像进行匹配融合。本发明使用卡尔曼滤波算法和平滑方法对放射源运动位置进行校正与优化,获得高准确,无跳动且运动轨迹连续的放射源定位;最后将放射源位置与光学图像进行匹配融合,可应用于移动放射源定位检测的多种场合。
Description
技术领域
本发明应用于编码孔径的伽马相机,涉及核安全检测领域,具体针对移动放射源检测的应用场景,是一种新的对移动放射源定位进行校正和优化的方法。
背景技术
基于编码孔径的伽马相机成像需要一定的统计计数才能得到信噪比高、真实可信的核辐射分布图像。传统的伽马相机,采用静止成像的方式,在探测器和放射源保持相对位置固定不变的情况下,经过长时间探测累积充分的计数进行成像,从而获得高质量的成像结果。
在运动成像过程中,由于探测器与放射源之间的相对位置不断变化,放射源在探测器视野中位置也不断发生变化。放射源的运动一方面将运动伪影引入数据中,另一方面数据的信号与噪声比值也发生变化。通常认为,核辐射分布图像中像素最大值所在位置即代表放射源位置。由于伽马相机采集数据运动伪影与噪声的干扰,辐射分布图像中放射源位置存在定位不准确的问题。在连续采集时放射源位置往往在某几帧发生跳动,会大幅度偏离正确位置。放射源被行人携带,伽马相机采集所得放射源位置与行人位置偏差较大,定位不准确。此外,由于辐射图像像素远小于光学图像像素,在进行热点与移动目标的匹配与融合时,热点定位的变化往往不够连续,不能紧跟移动的目标物体。
为了解决伽马相机在移动放射源成像中放射源定位跳动问题与定位不连续问题,本发明提供了基于卡尔曼滤波算法的放射源定位优化方法。针对移动放射源定位,一方面提高定位的连续性,另一方面可校正放射源定位偏离移动目标正确位置的数据点,使得伽马相机对移动放射源定位更连续与准确。
使用伽马相机对移动放射源进行检测与定位,传统方法通常使用逆互相关解码算法对采集数据进行解码重建,并叠加一段时间内的重建图像作为放射源在某点的辐射图像,其中选取辐射图像像素的最大值所在位置作为放射源位置。
传统方法未对辐射图像的放射源定位做优化,常出现放射源定位跳动与定位不连续问题。若通过提高探测器面积,在同等采集时长内可获取更多计数减少定位跳动,但势必会极大增加成本。
发明内容
针对伽马相机测量移动放射源的应用场景,本发明提供一种新的对移动放射源定位进行校正和优化的方法。对伽马相机采用逆互相关解码后的重建图像的放射源定位进行优化,从而提高放射源定位的准确性。同时该算法也可应用于移动放射源定位检测的多种场合。
本发明首先使用伽马相机或其他探测设备对移动反射源数据进行采集;其次通过逆相关解码获得放射源重建辐射图像;接着使用卡尔曼滤波算法和平滑方法对放射源运动位置进行校正与优化,获得高准确,无跳动且运动轨迹连续的放射源定位;最后将放射源位置与光学图像进行匹配融合。
本发明对技术方案为:
一种移动放射源定位校正与优化方法,其步骤包括:
1)获取放射源移动过程中的辐射图像,将每一所述辐射图像中像素最大值位置(x,y)作为放射源在辐射图像中的中心位置;
2)采用卡尔曼滤波算法,从x、y方向分别对各辐射图像中的放射源位置进行校正和优化;
3)将步骤2)校正和优化后的t时刻采集的辐射图像中的放射源位置与t时刻辐射图像对应的光学图像进行匹配融合。
进一步的,步骤2)中,从x方向对各辐射图像中的放射源位置进行校正和优化的方法为:
21)根据放射源的移动特性设置卡尔曼滤波算法中初始时刻t0的状态转移矩阵均方误差P0、过程噪声的协方差Q0和测量噪声的协方差R0;
22)采用卡尔曼滤波算法根据放射源x方向tk-1时刻的最优估计值xk-1和tk时刻的观测值zk,计算x方向tk时刻的状态最优估计值xk;
23)计算tk时刻的状态预测值和观测值zk之间的正则化残差ek,然后利用正则化残差ek调整tk时刻的过程噪声的协方差Qk和测量噪声的协方差Rk。
进一步的,当前时刻状态预测值和状态观测值的正则化差值ek表达式为: 其中,/>Hk-1为卡尔曼滤波算法tk-1时刻的滤波增益参数,Pk-1为tk-1时刻的均方误差,Rk-1为tk-1时刻的测量噪声的协方差;设定阈值为Thresh,当ek>Thresh时,调整当前时刻tk的过程噪声的协方差Qk和测量噪声的协方差Rk的值,使得Qk增大,Rk减小。
进一步的,步骤2)中,从y方向分别对各辐射图像中的放射源位置进行校正和优化的方法为:
21)根据放射源的移动特性设置卡尔曼滤波算法中初始时刻t0的状态转移矩阵均方误差P0、过程噪声的协方差Q0和测量噪声的协方差R0;
22)采用卡尔曼滤波算法根据放射源y方向tk-1时刻的最优估计值yk-1和tk时刻的观测值Lk,计算y方向tk时刻的状态最优估计值yk;
23)计算tk时刻的状态预测值和观测值Lk之间的正则化残差ek,然后利用正则化残差ek调整tk时刻的过程噪声的协方差Qk和测量噪声的协方差Rk。
进一步的,当前时刻状态预测值和状态观测值的正则化差值ek表达式为: 其中,/>Hk-1为卡尔曼滤波算法tk-1时刻的滤波增益参数,Pk-1为tk-1时刻的均方误差,Rk-1为tk-1时刻的测量噪声的协方差;设定阈值为Thresh,当ek>Thresh时,调整当前时刻tk的过程噪声的协方差Qk和测量噪声的协方差Rk的值,使得Qk增大,Rk减小。
进一步的,根据校正和优化后各辐射图像中的放射源位置生成该放射源的运动轨迹。
进一步的,采用平滑方法对放射源的运动轨迹进行平滑处理。
进一步的,所述平滑方法为RTS平滑方法。
进一步的,采集放射源移动过程中的光学图像,并对采集的光学图像进行逆相关解码获得该放射源移动过程中的辐射图像。
进一步的,使用伽马相机或探测设备采集放射源的光学图像。
与现有技术相比,本发明的积极效果为:
(1)在伽马相机对移动放射源检测和定位中,提高放射源定位的准确性;
(2)在伽马相机对移动放射源检测和定位中,提高放射源定位的连续性;
(3)在伽马相机对移动放射源检测和定位中,减少放射源定位跳动;
附图说明
图1为伽马相机对移动放射源采集并定位流程图;
图2为伽马相机逆相关解码重建辐射图像;
(a)t=0s,(b)t=8s,(c)t=9s,(d)t=10s,(e)t=11s,(f)t=12s;
图3伽马相机逆相关解码重建辐射图像中放射源位置;
(a)t=0s、x=7、y=6,(b)t=8s、x=6、y=7,(c)t=9s、x=5、y=7,(d)t=10s、x=7、y=1,(e)t=11s、x=3、y=8,(f)t=12s、x=2、y=9;
图4为对放射源x方向位置进行卡尔曼滤波效果图;
图5为对放射源y方向位置进行卡尔曼滤波效果图;
图6为对放射源x方向位置进行RTS滤波效果图;
图7为对放射源y方向位置进行RTS滤波效果图;
图8为t=7s时刻对比示意图;
(a)为滤波处理前匹配示意图,(b)为滤波处理后匹配示意图;
图9为t=8s时刻对比示意图;
(a)为滤波处理前匹配示意图,(b)为滤波处理后匹配示意图;
图10为t=10s时刻对比示意图;
(a)为滤波处理前匹配示意图,(b)为滤波处理后匹配示意图;。
具体实施方式
下面结合附图及实施例对本申请作进一步详细说明,但不仅限于此。
本发明提出了一种对移动放射源定位进行校正和优化的方法,可将该方法应用于移动放射源检测定位中。伽马相机对移动放射源采集、定位、优化、融合的流程如图1所示,主要用于对逆互相关解码算法重建出的辐射图像的定位进行校正和优化。
1伽马相机数据采集与解码重建
令小车或行人携带放射源移动,使用伽马相机对移动放射源进行检测和定位,获取数据并进行逆相关解码重建,所得的辐射图像如图2所示,分别展示t=0s,8s,9s,10s,11s,12s时刻的辐射图像。
2放射源定位
辐射图像长宽分别为19个像素大小,计算辐射图像中像素最大值所在位置,如图3所示,其中矩形框中的点为计算的放射源位置,以图片左上角为坐标原点,x与y即辐射图像中矩形框的中心坐标。
3定位优化
本方法主要采用卡尔曼滤波算法,对放射源位置的x,y方向分别进行定位的校正和优化。
卡尔曼滤波算法是一种利用线性系统状态方程,通过系统输入输出观测数据,对系统状态进行最优估计的算法。以伽马相机采集的辐射图像位置的x坐标为例。
已知前一时刻tk-1放射源位置的最优估计值为xk-1。对于当前时刻tk,观测得放射源位置为zk。计算tk时刻的放射源位置最优估计值xk,可利用前一时刻tk-1的最优估计值xk-1和当前时刻tk的观测值zk来更新对现时刻的最优估计值xk。
本系统的状态方程为:
其中为由前一时刻的放射源位置的最优估计值xk-1预测当前时刻状态的状态转移矩阵。wk-1为系统的过程噪声,表征了伽马相机对放射源定位系统中的系统固有噪声特点,Γk-1为本系统噪声驱动矩阵。该方程即意味当前时刻放射源的位置由前一时刻放射源的最优估计与系统固有噪声确定。
本系统的测量方程为:
zk=Ck*xk+vk
其中Ck为对当前时刻系统状态xk进行测量的测量矩阵,vk为测量过程的噪声。系统固有噪声与测量过程噪声满足
E[wk]=0 Cov[wk,wj]=Qkδkj
E[vk]=0 Cov[vk,vj]=Rkδkj
Cov[wk,vj]=0
其中Qk为当前时刻tk的过程噪声的协方差,Rk为当前时刻tk的测量噪声的协方差,参数Qk、Rk可根据物理模型的特点自己设定,可以是固定值,也可以进行调整。Qk与Rk可取固定值或根据系统物理特性进行调整。在本发明根据放射源运动加速度的变化Rk与Qk进行调整。
在已知伽马相机所测量的辐射热点前一时刻位置的最优估计xk-1和当前时刻位置的测量值zk,可以求得当前时刻的最优估计值xk。计算方法如下:
(1)根据前一时刻tk-1的位置最优估计值xk-1计算当前时刻tk的位置的状态预测
(2)计算当前时刻tk的状态预测均方误差
其中Pk-1为时刻tk-1的状态最优均方误差,Qk-1为时刻tk-1的过程噪声的协方差。Γk-1含义如上文所述,/>分别为/>Γk-1转置。
(3)计算当前时刻tk的滤波增益参数Hk:
其中为(2)式所得当前时刻tk的状态预测均方误差。Ck、Rk含义如上文所述,/>为Ck的转置。
(4)计算当前时刻tk位置的状态最优估计值xk:
其中Hk为(3)中所得当前时刻tk的滤波增益参数,zk、Ck含义如上文所述。
(5)更新当前时刻tk的状态最优估计的均方误差Pk:
其中Hk为(3)中所得当前时刻tk的滤波增益参数,为(2)式所得当前时刻tk的状态预测均方误差,Ck含义如上文所述。
根据以上原理,可分别对伽马相机采集所得的放射源x,y方向的位置坐标进行校正和优化,其中状态转移矩阵过程噪声的协方差Qk,测量噪声的协方差Rk均根据实际应用场景选择相应的值。
具体步骤如下:
(1)根据放射源携带者运动的物理特性(即放射源的移动特性)设置初始时刻t0的状态转移矩阵均方误差P0、过程噪声的协方差Q0和测量噪声的协方差R0。
(2)对伽马相机采集所得的放射源x,y方向的坐标分别使用上述原理进行优化。即根据放射源x方向tk-1时刻的最优估计值xk-1和tk时刻的观测值zk,计算x方向tk时刻的状态最优估计值xk。同样根据放射源y方向tk-1时刻的最优估计值yk-1和tk时刻的观测值Lk,计算y方向tk时刻的状态最优估计值yk。
(3)与传统卡尔曼滤波算法不同的是,在本发明中,通过计算当前时刻tk的状态预测值和状态观测值zk之间的正则化残差ek,来调整当前时刻的过程噪声的协方差Qk和测量噪声的协方差Rk。当前时刻状态预测值和状态观测值的正则化差值ek表达式如下:
设定阈值为Thresh,当ek>Thresh时,调整当前时刻的过程噪声的协方差Qk和测量噪声的协方差Rk的值,使得Qk增大,Rk减小,从而提升卡尔曼滤波器对当前时刻状态的计算。对伽马相机采集所得的放射源x,y方向位置坐标优化效果如图4、图5所示。图中黑色圆圈表示伽马相机测量获取的放射源位置,该放射源位置是离散的,且存在偏差较大的错误数据点。黑色线条表示通过卡尔曼滤波并且平滑后的放射源位置,可看出该曲线是连续变化且不存在偏差较大的点。卡尔曼滤波算法校正了偏离大错误的数据点,且使放射源轨迹随时间连续,提高放射源定位的准确性。
可进一步对卡尔曼滤波后的曲线进行优化,即根据卡尔曼滤波曲线的前后一段时间内的值来优化当前时刻的值。平滑的方法有固定点平滑(Fixed-Point Smoothing)、固定时延平滑(Fixed-Lag Smoothing)、固定时间间隔平滑(Fixed-Interval Smoothing)等方式。本发明以固定时间间隔平滑中的RTS平滑方法为例进行说明。
RTS平滑原理是,记录卡尔曼滤波器各个时刻输出的最优状态估计和对应的协方差,并遍历当前时刻之后的数据,把对未来的预测信息结合到当前时刻,对当前时刻的最优状态估计进行平滑。对以上数据进行RTS滤波所得结果如图6、图7所示,虚线所示为RTS滤波的结果。通过RTS滤波,放射源位置信息得到了进一步的平滑和优化,较卡尔曼滤波输出结果更为准确。
4匹配融合
使用伽马相机获取放射源位置(x,y),使用本专利中卡尔曼滤波RTS平滑方法对放射源位置(x,y)进行处理获得的优化后的放射源位置(x′,y′)。分别将放射源位置(x,y)和优化后的放射源位置(x′,y′)与光学图像中目标行人位置进行匹配融合,所得融合示意图对比如图8~10所示。根据实验对比结果可以发现,通过卡尔曼滤波器平滑后的放射源位置与光学图像中放射源位置更符合,说明基于卡尔曼滤波算法的移动放射源定位校正与优化方法可校正偏离放射源真实位置的测量点,能有效减少放射源定位的跳动并且使放射源运动轨迹更连续与平滑。
以上内容是结合具体的实施方式对本申请所作的进一步详细说明,不能认定本申请的具体实施只局限于这些说明。对于本申请所属技术领域的普通技术人员来说,在不脱离本申请发明构思的前提下,还可以做出若干简单推演或替换。
Claims (6)
1.一种移动放射源定位校正与优化方法,其步骤包括:
1)获取放射源移动过程中的辐射图像,将每一所述辐射图像中像素最大值位置(x,y)作为放射源在辐射图像中的中心位置;
2)采用卡尔曼滤波算法,从x、y方向分别对各辐射图像中的放射源位置进行校正和优化;
a)从x方向对各辐射图像中的放射源位置进行校正和优化的方法为:21a)根据放射源的移动特性设置卡尔曼滤波算法中初始时刻t0的状态转移矩阵均方误差P0、过程噪声的协方差Q0和测量噪声的协方差R0;22a)采用卡尔曼滤波算法根据放射源x方向tk-1时刻的最优估计值xk-1和tk时刻的观测值zk,计算x方向tk时刻的状态最优估计值xk;23a)计算tk时刻的状态预测值/>和观测值zk之间的正则化残差ek,然后利用正则化残差ek调整tk时刻的过程噪声的协方差Qk和测量噪声的协方差Rk;/>其中,Hk-1为卡尔曼滤波算法tk-1时刻的滤波增益参数,Pk-1为tk-1时刻的均方误差,Rk-1为tk-1时刻的测量噪声的协方差;设定阈值为Thresh,当ek>Thresh时,调整当前时刻tk的过程噪声的协方差Qk和测量噪声的协方差Rk的值,使得Qk增大,Rk减小;
b)从y方向分别对各辐射图像中的放射源位置进行校正和优化的方法为:21b)根据放射源的移动特性设置卡尔曼滤波算法中初始时刻t0的状态转移矩阵均方误差P0、过程噪声的协方差Q0和测量噪声的协方差R0;22b)采用卡尔曼滤波算法根据放射源y方向tk-1时刻的最优估计值yk-1和tk时刻的观测值Lk,计算y方向tk时刻的状态最优估计值yk;23b)计算tk时刻的状态预测值/>和观测值Lk之间的正则化残差ek,然后利用正则化残差ek调整tk时刻的过程噪声的协方差Qk和测量噪声的协方差Rk;/>其中,Hk-1为卡尔曼滤波算法tk-1时刻的滤波增益参数,Pk-1为tk-1时刻的均方误差,Rk-1为tk-1时刻的测量噪声的协方差;设定阈值为Thresh,当ek>Thresh时,调整当前时刻tk的过程噪声的协方差Qk和测量噪声的协方差Rk的值,使得Qk增大,Rk减小;
3)将步骤2)校正和优化后的t时刻采集的辐射图像中的放射源位置与t时刻辐射图像对应的光学图像进行匹配融合。
2.如权利要求1所述的方法,其特征在于,根据校正和优化后各辐射图像中的放射源位置生成该放射源的运动轨迹。
3.如权利要求2所述的方法,其特征在于,采用平滑方法对放射源的运动轨迹进行平滑处理。
4.如权利要求3所述的方法,其特征在于,所述平滑方法为RTS平滑方法。
5.如权利要求1所述的方法,其特征在于,采集放射源移动过程中的光学图像,并对采集的光学图像进行逆相关解码获得该放射源移动过程中的辐射图像。
6.如权利要求5所述的方法,其特征在于,使用伽马相机或探测设备采集放射源的光学图像。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010213464.6A CN111524167B (zh) | 2020-03-24 | 2020-03-24 | 一种移动放射源定位校正与优化方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010213464.6A CN111524167B (zh) | 2020-03-24 | 2020-03-24 | 一种移动放射源定位校正与优化方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111524167A CN111524167A (zh) | 2020-08-11 |
CN111524167B true CN111524167B (zh) | 2023-11-24 |
Family
ID=71901041
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010213464.6A Active CN111524167B (zh) | 2020-03-24 | 2020-03-24 | 一种移动放射源定位校正与优化方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111524167B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112365459B (zh) * | 2020-11-05 | 2024-07-23 | 中国科学院高能物理研究所 | 一种放射源位置与光学图像行人位置关联匹配的方法 |
CN112799116A (zh) * | 2020-12-22 | 2021-05-14 | 武汉第二船舶设计研究所(中国船舶重工集团公司第七一九研究所) | 一种基于互相关技术提高寻源距离的方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102854522A (zh) * | 2012-08-23 | 2013-01-02 | 成都理工大学 | 基于双重遗忘卡尔曼滤波的核辐射脉冲基线估计方法 |
EP3034984A1 (de) * | 2014-12-19 | 2016-06-22 | MBDA Deutschland GmbH | Verfahren und vorrichtung zur lokalen stabilisierung eines strahlungsflecks auf einem entfernten zielobjekt |
CN108931806A (zh) * | 2017-05-25 | 2018-12-04 | 中国辐射防护研究院 | 一种用于无线定位个人辐射剂量仪的定位坐标算法 |
US10365219B1 (en) * | 2018-02-06 | 2019-07-30 | Xidian University | Two-mode Raman optical projection tomography system |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7755057B2 (en) * | 2007-03-07 | 2010-07-13 | General Electric Company | CT gantry mounted radioactive source loader for PET calibration |
-
2020
- 2020-03-24 CN CN202010213464.6A patent/CN111524167B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102854522A (zh) * | 2012-08-23 | 2013-01-02 | 成都理工大学 | 基于双重遗忘卡尔曼滤波的核辐射脉冲基线估计方法 |
EP3034984A1 (de) * | 2014-12-19 | 2016-06-22 | MBDA Deutschland GmbH | Verfahren und vorrichtung zur lokalen stabilisierung eines strahlungsflecks auf einem entfernten zielobjekt |
CN108931806A (zh) * | 2017-05-25 | 2018-12-04 | 中国辐射防护研究院 | 一种用于无线定位个人辐射剂量仪的定位坐标算法 |
US10365219B1 (en) * | 2018-02-06 | 2019-07-30 | Xidian University | Two-mode Raman optical projection tomography system |
Non-Patent Citations (2)
Title |
---|
任波 ; 冯国强 ; 周焘 ; 于雷 ; .基于双波段红外辐射信息的机载无源定位算法.系统仿真学报.(第06期),第1427-1433页. * |
董广云 ; 张晓莉 ; 李长星 ; .放射源在线自动监控技术的研究.石油仪器.(第06期),第10-12页、18页. * |
Also Published As
Publication number | Publication date |
---|---|
CN111524167A (zh) | 2020-08-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110796010B (zh) | 一种结合光流法和卡尔曼滤波的视频稳像方法 | |
US8508605B2 (en) | Method and apparatus for image stabilization | |
JP3934151B2 (ja) | 画像生成装置および画像生成方法 | |
US8422795B2 (en) | Quality evaluation of sequences of images | |
US8368766B2 (en) | Video stabilizing method and system using dual-camera system | |
US9094606B2 (en) | Motion compensation in range imaging | |
US6628715B1 (en) | Method and apparatus for estimating optical flow | |
CN111524167B (zh) | 一种移动放射源定位校正与优化方法 | |
US20150254872A1 (en) | Method for Aligning and Tracking Point Regions in Images with Radial Distortion that Outputs Motion Model Parameters, Distortion Calibration, and Variation in Zoom | |
CN109272520B (zh) | 一种联合运动指导与边缘检测的自适应红外焦平面非均匀校正方法 | |
JP5107409B2 (ja) | 動き領域の非線形スムージングを用いた動き検出方法およびフィルタリング方法 | |
US10614554B2 (en) | Contrast adaptive video denoising system | |
CN102968765B (zh) | 一种基于sigma滤波器的红外焦平面非均匀性校正方法 | |
CN106017695B (zh) | 基于运动状态估计的自适应红外非均匀性校正方法 | |
WO2014069103A1 (ja) | 画像処理装置 | |
CN115170619B (zh) | 一种基于密集光流法的云遮挡预测方法 | |
CN110782403A (zh) | 一种红外图像非均匀性校正方法 | |
CN109978908B (zh) | 一种适应大尺度形变的单目标快速跟踪定位方法 | |
CN103868601B (zh) | Irfpa探测器非均匀响应的双边全变分正则化校正方法 | |
CN104980759A (zh) | 侦测遮挡区域的方法以及装置 | |
US10339640B2 (en) | Image processing apparatus, image processing method, and storage medium | |
Zhang et al. | Dehazing with improved heterogeneous atmosphere light estimation and a nonlinear color attenuation prior model | |
CN108830884B (zh) | 一种多视觉传感器协同目标跟踪方法 | |
US20110043649A1 (en) | Image processing apparatus, image processing method,and image processing program | |
JP2590705B2 (ja) | 動き補償予測装置 |
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 |