CN101785663B - 一种光声与x光检测双模态数字化成像***及成像方法 - Google Patents
一种光声与x光检测双模态数字化成像***及成像方法 Download PDFInfo
- Publication number
- CN101785663B CN101785663B CN2010101197184A CN201010119718A CN101785663B CN 101785663 B CN101785663 B CN 101785663B CN 2010101197184 A CN2010101197184 A CN 2010101197184A CN 201010119718 A CN201010119718 A CN 201010119718A CN 101785663 B CN101785663 B CN 101785663B
- Authority
- CN
- China
- Prior art keywords
- imaging
- optoacoustic
- ray
- subsystem
- ray check
- 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.)
- Expired - Fee Related
Links
Images
Landscapes
- Analysing Materials By The Use Of Radiation (AREA)
Abstract
本发明公开了一种光声与X光检测双模态数字化成像***及成像方法。本发明所述的光声与X光检测双模态数字化成像***包含光声成像子***、X光检测子***以及计算机;光声成像子***与X光检测子***通过共同的成像腔集成于一体;计算机处理采集到的光声成像和X光检测子***的数据,融合光声和X光检测信息,重建出带有光声与X光检测双重信息的图像。本发明利用高分辨率的X成像获取组织的结构信息;利用光声成像的高光谱选择性的获取组织的功能信息,实现双模态数字化成像,克服了现有单一成像模态的缺点与不足,结合了光声成像和X光检测成像的优点,得到原理上互补的、定位准确及分辨率高的成像***,造价较为低廉,易于推广。
Description
技术领域
本发明涉及一种多模态成像方法,特别涉及一种光声与X光检测双模态数字化成像***及成像方法,属于激光、X射线、电子以及图像重建在医学应用的综合技术领域。
背景技术
成像技术和手段的发展使得医学和药物研究可以在生物机体理各层面上观察结构和活动的过程。其目标瞄向了通过结构和功能成像来研究生物机制和疾病诊断及治疗方法。单一成像模态往往无法同时提供组织的结构与功能信息成像***,但由于每一种成像模态都有基于特定物质波和人体相互作用的规律,且这些规律不同;基于这些不同规律研制的成像技术和设备观测人体时得到的人体的信息也就不完全相同。因此,多模态成像可同时提供组织的结构和功能信息,是目前生物医学成像的发展趋势。
X光成像利用其在组织穿透过程中受到的吸收即衰减差异,能得到一般可见光不能穿透的各种不同密度物质的影像,可提供高分辨率的组织结构信息。数字化X光成像(Digital Radiography imaging)克服了传统X光成像使用胶片的缺点,成像分辨率高,对比度好。
光声成像(Photoacoustic imaging)是基于光声效应的新型无损医学成像方法。当脉冲激光照射到生物组织中时,组织吸收光能量引起温升,温升导致热膨胀而产生光声信号。生物组织的光声信号携带了组织的光吸收特征信息,通过测量光声信号能重建出组织中的光吸收分布图像。由于激发光的光谱选择性,可通过多波长实现反映组织生理和代谢活动的功能成像。
因此,利用高分辨的X光图像可以实现组织的精确定位;光声成像可提供光谱选择性的组织生理和病理改变的功能信息。因此,对于如何实现集成X光与光声双模态的成像方法于同一个***,有效利用两种成像模态优势互补的特性,对于同时研究组织结构、功能变化有着非常重要的现实意义。
发明内容
本发明的首要目的在于克服现有技术的单一成像模态的缺点与不足,提供一种光声与X光检测双模态数字化成像的***。所述的***通过共同的成像腔和特有的机械结构布局集成一体,可同时实现光声与X光检测双模态数字化成像。
本发明的另一目的在于提供所述***的成像方法。
本发明的目的通过下述技术方案实现:一种光声与X光检测双模态数字化成像***,包含光声成像子***、X光检测子***以及计算机;光声成像子***与X光检测子***通过共同的成像腔集成于一体;计算机处理采集到的光声成像子***和X光检测子***的数据,得到带有光声与X光检测双重信息的图像;
所述光声成像子***包含激光器、光声探测器和信号放大器;激光器、光声探测器和信号放大器依次连接;激光器作为光声激发光源产生脉冲激光,进入成像腔中聚焦到载物台上的样品区域;脉冲激光激发产生的光声信号经过声耦合液,通过光声探测器接收,送入信号放大器,通过高速数据采集卡采集,采集的信号储存到计算机缓存中等待图像重建软件进行图像重建,其中:
激光器位于成像腔的上方,输出脉冲激光的波长为400~2500nm;
光声探测器为均匀分布在同一圆周上的3个以上多元阵列探测器,每个阵列探测器的主频为20KHz~100MHz,阵元数为64,128或256;
所述的X光成像子***含有X射线源、光电转换器CMOS、信号放大器和A/D转换器;X射线源、光电转换器CMOS、信号放大器和A/D转换器依次连接;X射线源产生X射线进入成像腔,透过载物台上的样品,衰减后的X光被成像腔底部的光电转换器CMOS接收;接收的信号经信号放大器、A/D转换器转换后输入计算机中等待图像重建软件进行图像重建,其中:
X射线源位于成像腔的上方,工作电流为0.5mA~1mA,电压为5~45kVp可调;
光电转换器CMOS位于成像腔的底部,采用主动面阵CMOS光电转换器件,灵敏度高、像素量大(1兆~15兆),噪声低;
所述的成像腔含有载物台、光声探测器和光电转换器,其中:
载物台包含有机玻璃材质的透明水槽,透明塑料薄膜和小动物实验台;透明塑料薄膜将位于透明水槽底部的空洞密封,小动物实验台通过螺旋与透明水槽底部连接,通过螺旋可调节动物在成像腔中的高度;
光声探测器为3个以上的多元阵列探测器,固定在透明水槽壁,均匀分布在同一圆周上;探测的光声信号经放大器放大后,采集并储存到计算机缓存中;激光光源产生的脉冲激光进入成像腔并入射到载物台上的样品,激发产生的光声信号被光声探测器接收;每个阵列探测器的主频为20KHz~100MHz,阵元数为64,128或256;
光电转换器位于小动物实验台的下方;X射线源产生X射线透过载物台上的样品,衰减后的X射线被光电转换器CMOS接收;
所述的透明水槽装有声耦合液;脉冲激光激发产生的光声信号通过透明玻璃水槽,与其中的声耦合液耦合,进入光声阵列探测器;
所述的声耦合液为水;
所述的计算机含有高速数据采集卡、基于LABVIEW软件平台的控制软件和MATLAB平台的图像重建软件;
所述的采集采用基于LABVIEW控制平台的多通道并行实时采集***,其作用是将模拟信号转换为数字信号;
所述的处理是指通过MATLAB程序对光声采用小角度滤波反投影算法,X光采用直线投影算法,进行图像重建,融合两种模态的数据,从而得到待测部位的光声与X光检测双模态图像。
上述光声与X光检测双模态数字化成像***的成像方法包括以下步骤:
(1)激光器发出脉冲激光,进入成像腔中聚焦到载物台上的样品区域;脉冲激光激发产生的光声信号经过声耦合液,通过光声探测器接收,信号放大器放大;
(2)X射线源产生X射线进入成像腔,透过载物台上的样品,衰减后的X光被成像腔底部的光电转换器CMOS接收;接收的信号经信号放大器放大,再由A/D转换器转换;
(3)采集步骤(1)得到的放大的光声信号数据和步骤(2)得到的A/D转换后的信号,通过数据处理软件处理,得到光声与X光检测双模态图像。
所述的数据处理软件为MATLAB程序;
所述的采集采用基于LABVIEW控制平台的多通道并行实时采集***,其作用是将模拟信号转换为数字信号;
所述的处理是指通过MATLAB程序对光声采用小角度滤波反投影算法,X光采用直线投影算法,进行图像重建,融合两种模态的数据,从而得到待测部位的光声与X光检测双模态图像。
本发明并不直接用于疾病的诊断和治疗,只是为组织结构和功能活动参数改变提供提供检测数据;而且,本发明不仅限于此,本发明还可用于其他工业和生产中的无损检测等。
本发明利用高分辨率的X成像获取组织的结构信息;利用光声成像的高光谱选择性的获取组织的功能信息,实现双模态数字化成像。所述光声成像与数字化X光检测***通过共同的成像腔和特有的机械结构布局集成一体,可同时实现光声与X光检测双模态数字化成像。
本发明相对于现有技术具有如下的优点及效果:
(1)本发明中,X光成像可获取高分辨率的组织结构信息;光声成像的高光谱选择性可获取组织功能信息,因此可实现信息上互补的数字化双模态成像,有效克服了单一成像模态信息不足的缺点。
(2)本发明在结构上,有机的将光声成像和X光成像集成到一个***,可同步进行光声成像和荧光检测,使用方便。
(3)本发明集光声成像和X光检测于一体的关键部件是成像腔,共同空间定位结构,有利于双模态成像的信息融合和图像重建定位。
(4)本发明的各组件的造价较低,所以整体装置的造价亦相对较低,没有特殊限制,应用广泛。
附图说明
图1是本发明的框架图。
图2是本发明的结构图。
图3是本发明所述的成像腔的示意图。
图4是实施例2利用图1所示***得到了光声与X光检测图像:
(a)是单一的X光检测图像;
(b)是单一的光声图像;
(c)是光声与X光检测融合图像。
具体实施方式
下面结合实施例及附图对本发明作进一步详细的描述,但本发明的实施方式不限于此。
实施例1
图1为本发明所述光声与X光检测双模态数字化成像的***的框架图。
图2是本发明的结构图,包含光声成像子***、X光检测子***以及计算机;光声成像子***与X光检测子***通过共同的成像腔集成于一体;采用基于LABVIEW控制平台的多通道并行实时采集***,采集光声信号和X射线信号,计算机处理光声成像和X光检测子***的数据,利用基于MATLAB软件平台的光声小角度滤波反投影算法和X光直线投影算法,融合光声和X光检测信息,重建出带有光声与X光检测双重信息的图像。
光声成像的过程为(如图2和图3所示):位于成像腔上方的激光器1-1作为光声激发光源产生脉冲激光(激光器产生的激光脉冲触发采集卡同步采集光声信号),进入成像腔4中激发产生的光声信号,通过光声探测器1-2接收,送入信号放大器1-3,通过可位于计算机外部也可位于计算机内部的高速数据采集卡1-4,进入计算机3,利用基于Matlab软件平台的光声小角度滤波反投影算法图像重建软件进行图像重建。其中,激光器1-1选用Nd:YAG泵浦的OPO激光器(Vibrant 532I,Opotek,Carlsbad,Calif.),输出激光波长为532nm,690~960nm,脉宽为10ns,重复频率是10Hz。光声探测器1-2为多元阵列探测器EZU-PL21(Hitachi,Japan),由320个振元组成并被分割成64个振群,每5个振元组成一个振群。它的中心频率为3.5MHz,扫描宽度为102mm。高速数据采集卡1-4选用Gage Applied公司的Compuscope 12100型高速数据采集卡(DAS-CARD),采用外部触发方式,采样速率可达100MHz;采集控制程序用LABVIEW软件实现,数字图形信号处理用MATLAB软件实现。成像时利用小角度滤波反投影得到光声二维图像。
X光成像的过程为(如图2和图3所示):位于成像腔上方的X射线源2-1产生X射线进入到成像腔4,透过载物台上的小动物样品,衰减后的X光被成像腔底部的光电转换器(CMOS)2-2接收;接收的信号经信号放大器2-3、A/D转换器2-4转换后输入到带有图像重建和扫描控制软件的计算机3。其中,X射线源工作电流为0.5mA-1mA,电压为5-45kVp可调。光电转换器CMOS,采用主动面阵CMOS光电转换器件,它具有灵敏度高,分辨率高,像素量大(1024*1024),噪声低等优点。
图3是实施例1所述本发明光声与X光检测双模态成像***集成的关键部件——成像腔示意图。
成像腔4成像腔含有载物台、光声探测器及光电转换器。X射线源2-1产生X射线进入成像腔透过载物台上的样品,衰减后的X射线被光电转换器(CMOS)2-2接收。激光器1-1产生的脉冲激光进入成像腔并入射到载物台上的样品,激发产生的光声信号被光声探测器1-2接收。载物台,包含透明玻璃水槽4-1,透明塑料薄膜4-2和小动物实验台4-3;在透明玻璃水槽4-1底部留有空洞,空洞上密封有透明的塑料薄膜4-2,小动物实验台4-3通过螺旋与玻璃水槽底部连接,通过螺旋可调节动物在成像腔中的高度;透明玻璃水槽4-1装有声耦合液4-4,脉冲激光激发产生的光声信号通过透明玻璃水槽4-1,与其中的声耦合液4-4耦合,进入光声探测器;光声探测器为3个以上的多元阵列探测器,固定在玻璃水槽壁,均匀分布在同一圆周上;探测的光声信号经放大器放大后,采集并储存到计算机缓存中。
实施例2
应用实施例1所述的***进行模拟样品的光声与X光检测双成像。试验样品为一个带有移植肿瘤的Balb/c小鼠。实验中Balb/c小鼠被麻醉后,放置在小动物实验台4-3上。透明玻璃水槽内充满水4-4,作为光声信号的耦合液。激光器1-1输出激光波长为532nm,脉宽为10ns,重复频率是10Hz,进入成像腔4中并入射到载物台的样品上;脉冲激光激发产生的光声信号通过光声探测器1-2接收,送入信号放大器1-3,通过高速数据采集卡1-4采集,采集的信号在计算机3内部利用通过MATLAB程序对光声采用小角度滤波反投影算法进行图像重建,可得到如图4(b)所示的光声图像。图4(b)给出了小鼠肿瘤部位高对比度的光声横向层析图像。X光检测成像时,X射线源2-1产生X射线进入成像腔4,透过载物台上的样品,衰减后的X光被成像腔底部的光电转换器CMOS(2-2)接收;接收的信号经信号放大器2-3、A/D转换器2-4转换后输入到带有图像重建和扫描控制软件的计算机3,可得到如图4(a)所示X光图像,可得到高分辨率的小鼠骨架等结构信息。利用计算机3通过利用基于Matlab软件平台的光声小角度滤波反投影算法和X光直线投影算法分别重建出光声和X光检测图像,通过图像融合处理可在同一图像中同时得到光声与X光检测双模态的图像(如图4(c)所示),从而得到组织信息更丰富,更全面。
上述实施例为本发明较佳的实施方式,但本发明的实施方式并不受上述实施例的限制,其他的任何未背离本发明的精神实质与原理下所作的改变、修饰、替代、组合、简化,均应为等效的置换方式,都包含在本发明的保护范围之内。
Claims (9)
1.一种光声与X光检测双模态数字化成像***,其特征在于:所述的***包含光声成像子***、X光检测子***以及计算机;光声成像子***与X光检测子***通过共同的成像腔集成于一体;计算机处理采集到的光声成像和X光检测子***的数据,融合光声和X光检测信息,重建出带有光声与X光检测双重信息的图像;
所述光声成像子***包含激光器、光声探测器和信号放大器;激光器、光声探测器和信号放大器依次连接;
所述的X光检测子***含有X射线源、光电转换器CMOS、信号放大器和A/D转换器;X射线源、光电转换器CMOS、信号放大器和A/D转换器依次连接;
所述的成像腔含有载物台、光声探测器和光电转换器CMOS,光电转换器位于载物台的下方,光声探测器设置在载物台上;
所述光声成像子***中的激光器和所述X光检测子***中的X射线源位于所述成像腔的上方,与所述光电转换器分别居于所述载物台的两端。
2.根据权利要求1所述的光声与X光检测双模态数字化成像***,其特征在于:所述的激光器输出脉冲激光的波长为400~2500nm。
3.根据权利要求1所述的光声与X光检测双模态数字化成像***,其特征在于:所述的光声探测器为均匀分布在同一圆周上的3个以上多元阵列探测器,每个阵列探测器的主频为20KHz~100MHz,阵元数为64,128或256。
4.根据权利要求1所述的光声与X光检测双模态数字化成像***,其特征在于:所述X射线源的工作电流为0.5mA~1mA,电压为5~45kVp可调。
5.根据权利要求1所述的光声与X光检测双模态数字化成像***,其特征在于:所述的载物台包含有机玻璃材质的透明水槽,透明塑料薄膜和小动物实验台;透明塑料薄膜将位于透明水槽底部的空洞密封,小动物实验台通过螺旋与透明水槽底部连接;
所述的光声探测器为3个以上多元阵列探测器,固定在透明水槽壁上,均匀分布在同一圆周上。
6.根据权利要求1所述的光声与X光检测双模态数字化成像***,其特征在于:所述的光电转换器为主动面阵光电转换器件。
7.根据权利要求1所述的光声与X光检测双模态数字化成像***,其特征在于:所述的计算机含有高速数据采集卡、控制软件和图像重建软件。
8.权利要求1~7任一项所述光声与X光检测双模态数字化成像***的成像方法,其特征在于包括以下步骤:
(1)激光器发出脉冲激光,进入成像腔中聚焦到载物台上的样品区域;脉冲激光激发产生的光声信号经过装在载物台中的声耦合液,通过光声探测器接收,送入所述光声成像子***中的信号放大器;
(2)X射线源产生X射线进入成像腔,透过载物台上的样品,衰减后的X光被成像腔底部的光电转换器接收;接收的信号经所述X光检测子***中的信号放大器放大和A/D转换器转换;
(3)采集步骤(1)得到的放大的光声信号数据和步骤(2)得到的A/D转换后的信号,通过数据处理软件处理,得到光声与X光检测双模态图像;所述的处理是指通过MATLAB程序对光声采用小角度滤波反投影算法,X光采用直线投影算法,进行图像重建,融合两种模态的数据,从而得到待测部位的光声与X光检测双模态图像。
9.根据权利要求8所述的成像方法,其特征在于:所述的采集采用基于LABVIEW控制平台的多通道并行实时采集***。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2010101197184A CN101785663B (zh) | 2010-03-09 | 2010-03-09 | 一种光声与x光检测双模态数字化成像***及成像方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2010101197184A CN101785663B (zh) | 2010-03-09 | 2010-03-09 | 一种光声与x光检测双模态数字化成像***及成像方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101785663A CN101785663A (zh) | 2010-07-28 |
CN101785663B true CN101785663B (zh) | 2011-07-20 |
Family
ID=42529150
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2010101197184A Expired - Fee Related CN101785663B (zh) | 2010-03-09 | 2010-03-09 | 一种光声与x光检测双模态数字化成像***及成像方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101785663B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107485408A (zh) * | 2017-09-13 | 2017-12-19 | 北京大学 | 一种实现xct和光声成像的双模成像***及其成像方法 |
Families Citing this family (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102331430A (zh) * | 2011-06-23 | 2012-01-25 | 南京林业大学 | 不规则形状肉异物检测*** |
CN102621115B (zh) * | 2012-02-27 | 2014-08-27 | 华南师范大学 | 一种共焦的光声荧光同时成像方法与装置 |
TWM458203U (zh) * | 2012-12-17 | 2013-08-01 | Ind Tech Res Inst | 光聲檢測器、光聲板與使用此光聲板的檢測器 |
WO2014118674A2 (en) | 2013-01-30 | 2014-08-07 | Koninklijke Philips N.V. | Imaging system with hyperspectral camera guided probe |
CN103364107B (zh) * | 2013-07-31 | 2015-01-14 | 北京航空航天大学 | 一种衰减自补偿的光纤拉曼电缆温度监测与报警*** |
CN103389273A (zh) * | 2013-08-01 | 2013-11-13 | 中国科学院自动化研究所 | 一种光声和光学融合的多模态成像*** |
CN103519789A (zh) * | 2013-10-24 | 2014-01-22 | 中国科学院自动化研究所 | 一种集成光声与x射线断层成像的分离式成像*** |
CN103690244B (zh) * | 2013-12-24 | 2016-08-10 | 中国科学院自动化研究所 | 一种光声与x射线断层成像融合的成像*** |
CN103829961A (zh) * | 2014-03-21 | 2014-06-04 | 南京大学 | 一种结合有限角x射线成像、超声成像的多模式光声成像方法 |
CN111067523A (zh) * | 2019-12-25 | 2020-04-28 | 东软医疗***股份有限公司 | 门控装置及其控制方法、装置、医疗设备*** |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1470218A (zh) * | 2003-06-18 | 2004-01-28 | 华南师范大学 | 多元阵列电子扫描生物组织光声层析成像的方法及其装置 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20060262903A1 (en) * | 2005-05-20 | 2006-11-23 | Diebold Roger M | Method and apparatus for photothermal modification of x-ray images |
US20080123083A1 (en) * | 2006-11-29 | 2008-05-29 | The Regents Of The University Of Michigan | System and Method for Photoacoustic Guided Diffuse Optical Imaging |
-
2010
- 2010-03-09 CN CN2010101197184A patent/CN101785663B/zh not_active Expired - Fee Related
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1470218A (zh) * | 2003-06-18 | 2004-01-28 | 华南师范大学 | 多元阵列电子扫描生物组织光声层析成像的方法及其装置 |
Non-Patent Citations (1)
Title |
---|
邢达,向良忠.生物组织的光声成像技术及其在生物医学中的应用.《激光与光电子学进展》.2007,第44卷(第8期),第26-33页. * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107485408A (zh) * | 2017-09-13 | 2017-12-19 | 北京大学 | 一种实现xct和光声成像的双模成像***及其成像方法 |
Also Published As
Publication number | Publication date |
---|---|
CN101785663A (zh) | 2010-07-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101785663B (zh) | 一种光声与x光检测双模态数字化成像***及成像方法 | |
CN101785662A (zh) | 一种集成光声与荧光双模态的成像***及成像方法 | |
EP2579777B1 (en) | Second generation hand-held optical imager | |
EP2806803B1 (en) | Laser optoacoustic ultrasonic imaging system (louis) and methods of use | |
CN102621115B (zh) | 一种共焦的光声荧光同时成像方法与装置 | |
CN105011909B (zh) | 被检体信息获取设备和被检体信息获取方法 | |
US8712504B2 (en) | Hand-held optical probe based imaging system with 3D tracking facilities | |
JP6399753B2 (ja) | 被検体情報取得装置、表示方法、およびプログラム | |
JP5911196B2 (ja) | 光音響イメージング装置 | |
US20090018451A1 (en) | Dynamic Sampling System and Method for In Vivo Fluorescent Molecular Imaging | |
CN100546541C (zh) | 基于多环阵列光声传感的便携式血糖检测仪及检测方法 | |
CN106560160A (zh) | 被检体信息获取装置及其控制方法 | |
WO2009073979A1 (en) | Three-dimensional photoacoustic imager and methods for calibrating an imager | |
CN107960983A (zh) | 被检体信息获得装置、显示方法、程序和处理装置 | |
JP6859180B2 (ja) | 光音響装置および光音響装置の制御方法 | |
JP2015205136A (ja) | 光音響装置、光音響装置の制御方法、及びプログラム | |
CN109475345A (zh) | 用于显示超声图像和光声图像的装置、方法和程序 | |
JP6168802B2 (ja) | 処理装置、処理方法、及びプログラム | |
CN107184180A (zh) | 信息处理***和显示控制方法 | |
CN104027132A (zh) | 一种基于多谱光声成像的装置和方法 | |
CN201211188Y (zh) | 基于多环阵列光声传感的便携式血糖检测仪 | |
CN108403083A (zh) | 信息处理装置、信息处理方法和非暂时性存储介质 | |
CN105640496B (zh) | 光声装置和被检体信息取得方法 | |
JP2016067926A (ja) | 光音響装置、情報処理方法、およびプログラム | |
JP6740004B2 (ja) | 光音響装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20110720 |
|
CF01 | Termination of patent right due to non-payment of annual fee |