CN104156908B - 高时间分辨率的频域瞬态成像方法 - Google Patents
高时间分辨率的频域瞬态成像方法 Download PDFInfo
- Publication number
- CN104156908B CN104156908B CN201410428227.6A CN201410428227A CN104156908B CN 104156908 B CN104156908 B CN 104156908B CN 201410428227 A CN201410428227 A CN 201410428227A CN 104156908 B CN104156908 B CN 104156908B
- Authority
- CN
- China
- Prior art keywords
- image
- frequency
- imaging system
- fundamental wave
- harmonic
- 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
Landscapes
- Image Processing (AREA)
Abstract
本发明公开了一种高时间分辨率的频域瞬态成像方法,包括以下步骤:1)通过调制成像***对被测对象采集一组图像集Ix,y(f,φ);2)从频率为f的所述图像集Ix,y(f,φ)中提取出基波图像和谐波图像In x,y(nf);3)对所述基波图像和所述谐波图In x,y(nf)像进行幅值和相位的校正;4)将校正后的所述基波图像和所述谐波图像In x,y(nf)组成频谱图像集Ip x,y(f),然后沿着所述频谱图像集的频率f做傅里叶变换得到瞬态图像ix,y(t)。该方法从采集得到的有限工作频率的图像数据中提取出高次谐波分量,获得更高频率的信息,从而提高瞬态图像的时间分辨率。
Description
技术领域
本发明涉及光学和计算机视觉领域,特别涉及一种高时间分辨率的频域瞬态成像方法。
背景技术
对于一般成像***来说光速可以认为是无穷大的,因此传统图像记录的是在光照稳定的条件下的对象的光响应。但是对象的光瞬态响应包含了更多的信息。采集对象的光瞬态响应的信息的技术称为瞬态成像。瞬态成像的一种方法是采用超高速的相机,并且需要超短脉冲光源。这种方法属于时域采集方案,其时间分辨率取决于相机的成像速度。缺点为所采用的超高速成像******和光源非常昂贵,而且采集过程复杂非常耗时。另一种方法是采用调制成像***,所需的***与TOF(飞行时间)深度相机原理相同,区别在于深度相机采用固定频率的脉冲进行成像,而瞬态成像所用的调制成像***需要用多个频率分别进行成像。这种方法属于频域采集方案,成像***在每个工作频率上采集的数据就是瞬态图像的在该频率的频谱,通过傅里叶变换及相关去噪算法或采用优化方法可以求解出瞬态图像。频域瞬态成像的时间分辨率主要取决于调制成像***的工作频率。由于硬件***的限制,调制成像***只能工作在有限的频率范围,这限制了得到的瞬态图像的时间分辨率。
发明内容
本发明的目的是为了克服现有技术中调制成像***工作的频率范围有限导致瞬态图像的时间分辨率较低的技术缺陷,发明了一种高时间分辨率的频域瞬态成像方法。
为解决上述技术问题,本发明提供了一种高时间分辨率的频域瞬态成像方法,该方法包括以下步骤:1)通过调制成像***对被测对象采集一组图像集Ix,y(f,φ);2)从频率为f的所述图像集Ix,y(f,φ)中提取出基波图像和谐波图像In x,y(nf);3)对所述基波图像和所述谐波图In x,y(nf)像进行幅值和相位的校正;4)将校正后的所述基波图像和所述谐波图像In x,y(nf)组成频谱图像集Ip x,y(f),然后沿着所述频谱图像集的频率f做傅里叶变换得到瞬态图像ix,y(t)。
上述技术方案中,步骤1)中所述调制成像***至少包括一个信号发生器、一个光源及一个成像***。
上述技术方案中,步骤1)中所述通过调制成像***对被测对象采集一组频率为f
的图像集Ix,y(f,φ)的步骤为:a)所述信号发生器控制所述光源以频率f发射方波信号,并
控制所述成像***的快门以同频率成像;b)采集所述光源的控制信号和所述成像***的快
门的控制信号的相位差为φ时得到的图像及φ-90°得到的图像组成所
述图像集其中j为虚数符号。
上述技术方案中,步骤1)中所述调制成像***选择从低到高的频率范围选择等间隔工作频率以及从0到2π的相位选择等间隔相位对被测对象进行采集图像。
上述技术方案中,步骤2)中从所述频率为f的图像集Ix,y(f,φ)中提取出基波图像和谐波图像In x,y(nf)的步骤为:
上述技术方案中,步骤3)中在对所述基波图像和所述谐波图In x,y(nf)像进行幅值和相位的校正之前,先对所述调制成像***进行标定,其步骤为:
A)选取一个表面均匀的漫反射板,放置在所述调制成像***之前;
B)对所述漫反射板进行数据采集,得到一组图像Cx,y(f,φ);
C)采用参数化模型拟合得到参数An
(nf)和
上述技术方案中,步骤4)将校正后的所述基波图像和所述谐波图像In x,y(nf)组成频谱图像集Ip x,y(f)的步骤中,如果一个频率点有多幅图像,则保留基波图像或最低次谐波的图像;如果一个频率点没有图像,则用幅值线性插值和相位线性插值补足该频率点的图像。
本发明具有如下有益效果:本发明的高时间分辨率的频域瞬态成像方法,从采集得到的有限工作频率的图像数据中提取出高次谐波分量,获得更高频率的信息,从而提高瞬态图像的时间分辨率。
下面通过附图和实施例,对本发明的技术方案做进一步的详细描述。
附图说明
附图用来提供对本发明的进一步理解,并且构成说明书的一部分,与本发明的实施例一起用于解释本发明,并不构成对本发明的限制。在附图中:
图1是本发明调制成像***示意图;
图2是本发明高时间分辨率的频域瞬态成像方法流程图。
具体实施方式
下面结合附图,对本发明的一个具体实施方式进行详细描述,但应当理解本发明的保护范围并不受具体实施方式的限制。
除非另有其它明确表示,否则在整个说明书和权利要求书中,术语“包括”或其变换如“包含”或“包括有”等等将被理解为包括所陈述的元件或组成部分,而并未排除其它元件或其它组成部分。
如图1所示,本发明提供了一种调制成像***,该***至少包括:一个信号发生器、一个光源和一个成像***;
如图2所示,为利用上述调制成像***实施的高时间分辨率的频域瞬态成像方法流程图,该方法包括以下步骤:
步骤S100:通过调制成像***对被测对象采集一组图像集Ix,y(f,φ),其中,f为工作频率,φ是可调相位;
更为具体地,首先信号发生器控制光源以频率f发射方波信号,并控制成像***的快门以同频率成像;若采集光源的控制信号和成像***的快门的控制信号的相位差为φ时得到的图像若采集光源的控制信号和成像***的快门的控制信号的相位差为φ-90°时得到的图像因此,采集得到的图像集为j为虚数符号;其中,该调制成像***选择从低到高的频率范围选择等间隔工作频率(即工作频率f=nfs,f∈[fL,fH],其中n为正整数,fs为频率间隔,fL,fH分别为调制成像***的最小和最大调制频率)以及从0到2π的相位选择等间隔相位(即φ=mφs,φs为相位间隔,m=0,1,...,K,K=2π/φs)对被测对象进行采集图像。
步骤S102:从频率为f的所述图像集Ix,y(f,φ)中提取出基波图像和谐波图像In x,y(nf),其中n为正整数;
其中,从频率为f的图像集Ix,y(f,φ)中提取出基波图像和谐波图像In x,y(nf)的具体方法步骤为:
步骤S104:对基波图像和谐波图In x,y(nf)像进行幅值和相位的校正;
在进行校正之前,首先对调制成像***进行标定(对于一个调制成像***,标定仅需执行一次),获得在每个工作频率下光源信号的基波和谐波的幅值An(nf)和相移然后将In x,y(nf)除以标定的步骤为:
1)选取一个表面均匀的漫反射板,放置在调制成像***之前;
2)对漫反射板进行数据采集,得到一组图像Cx,y(f,φ);
3)采用以下参数化模型拟合得到参数An(nf)和
步骤S106:将校正后的基波图像和谐波图像In x,y(nf)组成频谱图像集Ip x,y(f),然后沿着频谱图像集的频率f做傅里叶变换得到瞬态图像ix,y(t);
在将校正后的基波图像和谐波图像In x,y(nf)组成频谱图像集Ip x,y(f)的进程中,如果一个频率点有多幅图像,则保留基波图像或最低次谐波的图像;如果一个频率点没有图像,则用幅值线性插值和相位线性插值补足该频率点的图像,最终得到频率间隔为fs的频谱图像集Ip x,y(f)。
本发明的高时间分辨率的频域瞬态成像方法,从采集得到的有限工作频率的图像数据中提取出高次谐波分量,获得更高频率的信息,从而提高瞬态图像的时间分辨率。
下面对本发明实施例高时间分辨率的频域瞬态成像方法进行实例分析,具体为对一个对象进行瞬态成像的实例,步骤如下:
步骤200:设置调制成像***的工作频率5MHz~165MHz,间隔fs=1MHz,相位0~π,间隔φs=π/200,采集得到一组图像Ix,y(f,φ)。
步骤202:对于频率为f的图像集Ix,y(f,φ),提取出基波图像和谐波图像;
由于高次谐波幅值较小,只选择3次谐波,即n的最大值为3。
步骤204:用标定参数对基波图像和谐波图像进行校正;
步骤206:用校正后的基波图像和谐波图像组成频谱图像集Ip x,y(f),频率间隔为fs=1MHz。其中5~165MHz频段采用基波图像,166~330MHz频段采用二次谐波图像,331~495MHz频段采用三次谐波图像。0~5MHz频段和166~495MHz频段中不存在的图像,用幅值线性插值和相位线性插值补足。
步骤208:重构瞬态图像ix,y(t)=FFT[Ip x,y(f)]。
最后应说明的是:以上实施例仅用以说明本发明而并非限制本发明所描述的技术方案;因此,尽管本说明书参照上述的各个实施例对本发明已进行了详细的说明,但是,本领域的普通技术人员应当理解,仍然可以对本发明进行修改或等同替换;而一切不脱离本发明的精神和范围的技术方案及其改进,其均应涵盖在本发明的权利要求范围中。
Claims (1)
1.一种高时间分辨率的频域瞬态成像方法,其特征在于,包括以下步骤:
1)通过调制成像***对被测对象采集一组图像集Ix,y(f,φ);
2)从频率为f的所述图像集Ix,y(f,φ)中提取出基波图像和谐波图像In x,y(nf);
3)对所述基波图像和所述谐波图像In x,y(nf)进行幅值和相位的校正;
4)将校正后的所述基波图像和所述谐波图像In x,y(nf)组成频谱图像集Ip x,y(f),然后沿着所述频谱图像集的频率f做傅里叶变换得到瞬态图像ix,y(t);其中,如果一个频率点有多幅图像,则保留基波图像或最低次谐波的图像;如果一个频率点没有图像,则用幅值线性插值和相位线性插值补足该频率点的图像;
步骤1)中所述调制成像***至少包括一个信号发生器、一个光源及一个成像***;
步骤1)中所述通过调制成像***对被测对象采集一组图像集Ix,y(f,φ)的步骤为:
a)所述信号发生器控制所述光源以频率f发射方波信号,并控制所述成像***的快门以同频率成像;
b)采集所述光源的控制信号和所述成像***的快门的控制信号的相位差为φ时得到的图像及φ-90°得到的图像组成所述图像集其中j为虚数符号;
步骤1)中所述调制成像***选择从低到高的频率范围选择等间隔工作频率以及从0到2π的相位选择等间隔相位对被测对象进行采集图像;
步骤2)中从所述频率为f的图像集Ix,y(f,φ)中提取出基波图像和n次谐波图像In x,y(nf)的步骤为:
步骤3)中在对所述基波图像和所述谐波图像In x,y(nf)进行幅值和相位的校正之前,先对所述调制成像***进行标定,其步骤为:
A)选取一个表面均匀的漫反射板,放置在所述调制成像***之前;
B)对所述漫反射板进行数据采集,得到一组图像Cx,y(f,φ);
C)采用参数化模型拟合得到参数An(nf)和
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410428227.6A CN104156908B (zh) | 2014-08-27 | 2014-08-27 | 高时间分辨率的频域瞬态成像方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410428227.6A CN104156908B (zh) | 2014-08-27 | 2014-08-27 | 高时间分辨率的频域瞬态成像方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104156908A CN104156908A (zh) | 2014-11-19 |
CN104156908B true CN104156908B (zh) | 2018-11-02 |
Family
ID=51882399
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410428227.6A Active CN104156908B (zh) | 2014-08-27 | 2014-08-27 | 高时间分辨率的频域瞬态成像方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104156908B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2017054147A1 (zh) * | 2015-09-30 | 2017-04-06 | 温州医科大学 | 单次快照多频解调方法 |
CN105245761B (zh) * | 2015-09-30 | 2018-08-24 | 温州医科大学 | 单次快照多频解调方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1769925A (zh) * | 2004-11-05 | 2006-05-10 | 清华大学 | 一种合成孔径雷达运动目标成像方法 |
CN101623203A (zh) * | 2009-08-07 | 2010-01-13 | 西安交通大学 | 一种瞬态物理过程中多模式和多参量同步检测成像监控***及监控方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP5823175B2 (ja) * | 2011-06-06 | 2015-11-25 | 株式会社東芝 | 超音波診断装置、医用画像処理装置および医用画像処理プログラム |
-
2014
- 2014-08-27 CN CN201410428227.6A patent/CN104156908B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1769925A (zh) * | 2004-11-05 | 2006-05-10 | 清华大学 | 一种合成孔径雷达运动目标成像方法 |
CN101623203A (zh) * | 2009-08-07 | 2010-01-13 | 西安交通大学 | 一种瞬态物理过程中多模式和多参量同步检测成像监控***及监控方法 |
Non-Patent Citations (1)
Title |
---|
Fourier Analysis on Transient Imaging with a Multifrequency Time-of-Flight Camera;Jingyu Lin,et al;《2014 IEEE Conference on Computer Vision and Pattern Recognition》;20140628;参见摘要,正文第三部分 * |
Also Published As
Publication number | Publication date |
---|---|
CN104156908A (zh) | 2014-11-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US10362984B2 (en) | Measuring and monitoring skin feature colors, form and size | |
RU2020111209A (ru) | Устройство обработки информации, способ генерации изображения, способ управления и носитель данных | |
IL274426B2 (en) | Discovering and describing a fully complex point of interest using homographic processing | |
CN104156908B (zh) | 高时间分辨率的频域瞬态成像方法 | |
EP2860694A3 (en) | Image processing apparatus, image pickup apparatus, image processing method, image processing program, and non-transitory computer-readable storage medium | |
JP2014068330A5 (zh) | ||
EP3067862A3 (en) | Image processing method and apparatus | |
EP2518719A3 (en) | Image range expansion control methods and apparatus | |
CN103634544A (zh) | 一种投影方法及电子设备 | |
MY196475A (en) | Imaging Device, Electronic Device, and Method for Obtaining Image by The Same | |
SG11201909071UA (en) | Image processing methods and apparatuses, computer readable storage media and eletronic devices | |
US20220084225A1 (en) | Depth Map Processing Method, Electronic Device and Readable Storage Medium | |
EP3879494A3 (en) | Method, apparatus, electronic device, computer readable medium and program for calibrating external parameter of camera | |
CN103829966A (zh) | 用于自动确定侦测图像中的定位线的方法和*** | |
JP2016048904A5 (zh) | ||
SG11201907985VA (en) | Risk assessment device, risk assessment method, and risk assessment program | |
JP2017116420A5 (zh) | ||
CN104299231A (zh) | 多传感器实时图像配准方法和*** | |
CN105763785A (zh) | 动态调整快门速度的拍照***及方法 | |
Sych et al. | Web-based data processing system for automated detection of oscillations with applications to the solar atmosphere | |
CN106651962A (zh) | 一种基于投影仪‑相机***的自适应光补偿方法 | |
JP2012083233A5 (zh) | ||
CN105472245A (zh) | 一种拍照方法、电子设备 | |
AU2017260525A1 (en) | Method and system for body scanning and display of biometric data | |
WO2014203192A3 (en) | Cortical bone segmentation from mr dixon data |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |