CN104799817B - 应用于乳腺的正交三角波频率编码高光谱图成像测量系统 - Google Patents

应用于乳腺的正交三角波频率编码高光谱图成像测量系统 Download PDF

Info

Publication number
CN104799817B
CN104799817B CN201510196490.1A CN201510196490A CN104799817B CN 104799817 B CN104799817 B CN 104799817B CN 201510196490 A CN201510196490 A CN 201510196490A CN 104799817 B CN104799817 B CN 104799817B
Authority
CN
China
Prior art keywords
monochromatic light
light source
camera
monochromatic
mammary gland
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
Application number
CN201510196490.1A
Other languages
English (en)
Other versions
CN104799817A (zh
Inventor
李刚
杨雪
张启蕊
林凌
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Tianjin University
Original Assignee
Tianjin University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Tianjin University filed Critical Tianjin University
Priority to CN201510196490.1A priority Critical patent/CN104799817B/zh
Publication of CN104799817A publication Critical patent/CN104799817A/zh
Application granted granted Critical
Publication of CN104799817B publication Critical patent/CN104799817B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/43—Detecting, measuring or recording for evaluating the reproductive systems
    • A61B5/4306—Detecting, measuring or recording for evaluating the reproductive systems for evaluating the female reproductive systems, e.g. gynaecological evaluations
    • A61B5/4312—Breast evaluation or disorder diagnosis
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/0059—Measuring for diagnostic purposes; Identification of persons using light, e.g. diagnosis by transillumination, diascopy, fluorescence
    • A61B5/0075—Measuring for diagnostic purposes; Identification of persons using light, e.g. diagnosis by transillumination, diascopy, fluorescence by spectroscopy, i.e. measuring spectra, e.g. Raman spectroscopy, infrared absorption spectroscopy
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/0059—Measuring for diagnostic purposes; Identification of persons using light, e.g. diagnosis by transillumination, diascopy, fluorescence
    • A61B5/0082—Measuring for diagnostic purposes; Identification of persons using light, e.g. diagnosis by transillumination, diascopy, fluorescence adapted for particular medical purposes
    • A61B5/0084—Measuring for diagnostic purposes; Identification of persons using light, e.g. diagnosis by transillumination, diascopy, fluorescence adapted for particular medical purposes for introduction into the body, e.g. by catheters
    • A61B5/0086—Measuring for diagnostic purposes; Identification of persons using light, e.g. diagnosis by transillumination, diascopy, fluorescence adapted for particular medical purposes for introduction into the body, e.g. by catheters using infrared radiation

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Biomedical Technology (AREA)
  • Medical Informatics (AREA)
  • Biophysics (AREA)
  • Pathology (AREA)
  • Engineering & Computer Science (AREA)
  • Veterinary Medicine (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Public Health (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • General Health & Medical Sciences (AREA)
  • Reproductive Health (AREA)
  • Gynecology & Obstetrics (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Investigating Or Analysing Materials By Optical Means (AREA)

Abstract

本发明公开了一种应用于乳腺的正交三角波频率编码高光谱图成像测量系统，一组单色光源分布在乳腺组织的一侧，摄像头分布在乳腺组织的另一侧；一组单色光源中的各个单色光源密集排列在一个所设半球面上，采用透镜汇聚成一束光，构成光源；一个摄像头构成光源接收器件；采用不同频率的正交三角波分别驱动一组单色光源中的各个单色光源，摄像头接收到的图像中每个像素点为每个单色光源透过乳腺组织的单色光组合；计算机对单色光组合进行分离得到单色光组合中的各个单色光源的贡献，据此实现对乳腺的透射高光谱图成像。本发明实现了高速度、大信息的乳腺透射高光谱图成像的高精度测量，且本发明装置具有成本低、应用方便等优点，适宜于经常性家庭自检。

Description

应用于乳腺的正交三角波频率编码高光谱图成像测量系统
技术领域
本发明涉及成像测量系统领域，尤其涉及一种应用于乳腺的正交三角波频率编码高光谱图成像测量系统。
背景技术
现有技术中通过光对物体内部进行成像，特别是对人体内部进行成像，具有无损无创无辐射的突出优点，但至今尚未有能够进入家庭使用的乳腺透射成像光测量系统用于经常性的肿瘤自检，其原因在于现有乳腺成像测量系统的成本高、运算复杂、精度有限，无法满足实际应用中的乳腺自检需要。
发明内容
本发明提供了一种应用于乳腺的正交三角波频率编码高光谱图成像测量系统，本发明实现了高速度、大信息的乳腺透射光成像的高精度测量，满足了实际应用中的需要，详见下文描述：
一种应用于乳腺的正交三角波频率编码高光谱图成像测量系统，所述成像测量系统包括：一组单色光源、一个摄像头，以及与摄像头外接的计算机，一组单色光源分布在乳腺组织的一侧，摄像头分布在乳腺组织的另一侧；
其中，一组单色光源中的各个单色光源密集排列在一个所设半球面上，采用透镜汇聚成一束光，构成光源；一个摄像头构成光源接收器件；
采用不同频率的正交三角波分别驱动一组单色光源中的各个单色光源，摄像头接收到的图像中每个像素点是每个单色光源透过乳腺的单色光组合；
计算机对单色光组合进行分离得到单色光组合中的各个单色光源的贡献，据此实现对乳腺的透射高光谱图成像。
其中，所述单色光源、所述摄像头在乳腺组织样本两面对称设置。
所述计算机对单色光组合进行分离得到单色光组合中的各个单色光源的贡献，据此实现对乳腺的透射高光谱图成像的操作具体为：
以4种波长的发光二极管进行说明，λ1和λ2波长的发光二极管驱动正交三角波频率为2倍f，λ3和λ4波长的发光二极管驱动正交三角波频率为1倍f，相同驱动频率正交三角波的相位相差90；
摄像头采样频率为fS，且fS＝2f，在λ1驱动信号由低变高的过程中采样；
以顺序每8个数字信号为一组进行运算，分别得到4/3倍的波长λ1、λ2的光信号，2倍的波长λ3、λ4的光信号。
所述单色光源为激光二极管。
另一实施例，所述单色光源为单色二极管。
进一步地，所述摄像头为手机摄像头。
本发明提供的技术方案的有益效果是：本发明采用不同频率的正交三角波驱动单色光源，对摄像头检测到的光电信号进行分离可以得到单色光组合中的各个单色光源的贡献，进而实现对乳腺组织样本的成像，本发明实现了高速度、大信息的线阵成像光的高精度测量，且本系统具有成本低、应用方便等优点，适宜于家庭自检。
附图说明
图1为一种应用于乳腺的正交三角波频率编码高光谱图成像测量系统的结构示意图；
图2为本发明提供的单色光源、乳腺组织样本和摄像头相对位置示意图；
图3为采用不同频率的正交三角波分别驱动各个单色光源的示意图。
附图中，各标号所代表的部件列表如下：
1：一组单色光源； 2：透镜；
3：乳腺组织样本； 4：摄像头。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚，下面对本发明实施方式作进一步地详细描述。
实施例1
一种应用于乳腺的正交三角波频率编码高光谱图成像测量系统，参见图1和图2，成像测量系统包括：一组n个单色光源1(用LD表示)和一个摄像头4(CMOS或CCD，已能采集近红外波段的光)，单色光源1的波长处于600～1200nm之间，此波段为“光学窗口”，穿透深度比较深。一组单色光源1分布在乳腺组织样本3的一侧，摄像头4分布在乳腺组织样本3的另一侧。n的具体取值与乳腺样本3对不同波长的敏感程度相关，本发明实施例对此不做限制；
其中，一组单色光源1中的各个单色光源LD1…LDn密集排列在一个所设半球面上，采用透镜2汇聚成一束光，构成线阵光源；摄像头4构成多波长光源接收器件。该成像测量系统还包括与摄像头4外接的计算机(图中未示出)。
参见图3，采用不同频率的正交三角波分别驱动一组单色光源1中的各个单色光源LDi，摄像头4进行采集，采集到的图像中，每个像素点为每个单色光源LDi透过乳腺组织样本3的单色光组合Iij；计算机对单色光组合Iij进行分离可以得到单色光组合Iij中的各个单色光源LDi的贡献，构成各波长下的透射图像。根据各波长下乳腺组织中正常组织、肿瘤组织等的光学特性不同，进行乳腺的透射光成像测量，实现肿瘤的早期检测。
其中，计算机对单色光组合Iij进行分离可以得到单色光组合Iij中的各个单色光源LDi的贡献的步骤具体为：
为简便说明起见，以4种波长的发光二极管2为例进行说明，假定λ1(D1发光二极管)和λ2(D2发光二极管)波长的发光二极管2驱动正交三角波频率为2倍f，λ3(D3发光二极管)和λ4(D4发光二极管)波长的发光二极管2驱动正交三角波频率分别为1倍f，而且λ1、λ2波长的发光二极管2的驱动正交三角波频率相同但相位相差90°，λ3、λ4波长的发光二极管2的驱动正交三角波频率相同但相位相差90°。
假定摄像头4采样频率为fS，且fS＝2f，并保证在λ1驱动信号由低变高的过程中采样。
数字信号序列
Figure BDA0000704677310000031
可以表示为：
Figure BDA0000704677310000032
其中，
Figure BDA0000704677310000033
和
Figure BDA0000704677310000034
分别为波长λ1、λ2、λ3和λ4的电压信号，
Figure BDA0000704677310000035
为低频信号，包括：背景光、摄像头4的本底噪声。
假定采样频率fS远高于调制正交三角波信号和背景光的变化频率，在最低驱动信号频率的一个周期可以近似认为各路正交三角波信号的幅值和背景光信号的幅值不变。以最前8个采样值为例：
Figure BDA0000704677310000041
其中，
Figure BDA0000704677310000042
和
Figure BDA0000704677310000043
分别为波长λ1、λ2、λ3和λ4的光信号和背景信号的幅值。
换言之，以顺序每8个数字信号为一组进行运算：
Figure BDA0000704677310000044
即得到4/3倍的波长λ1的光信号
Figure BDA0000704677310000045
而且完全消除了背景信号
Figure BDA0000704677310000046
的影响。
Figure BDA0000704677310000047
即得到4/3倍的波长λ2的光信号
Figure BDA0000704677310000048
而且完全消除了背景信号
Figure BDA0000704677310000049
的影响。
Figure BDA00007046773100000410
即得到2倍的波长λ3的光信号
Figure BDA00007046773100000411
而且完全消除了背景信号
Figure BDA00007046773100000412
的影响。
Figure BDA00007046773100000413
即得到2倍的波长λ4的光信号
Figure BDA00007046773100000414
而且完全消除了背景信号
Figure BDA00007046773100000415
的影响。
对幅值为x的被采样值，如果在一定的时间内均匀采样N(>>1)点并进行平均，得到的平均值是
Figure BDA00007046773100000416
其中，[x]是模数转换器对x进行量化，也即按四舍五入圆整得到的正整数。xi是第i点的幅值，[xi]是模数转换器对xi进行量化，也即按四舍五入圆整得到的正整数。
(7)式表明，对一个比较“干净”的信号采样多次进行平均，并不能提高其精度，所得到的平均值的误差与单次采样的误差相同，为Δxi。
如果对幅值为x的被采样锯齿波，同样在一定的时间内均匀采样N(>>1)点并进行平均，得到的平均值是
Figure BDA0000704677310000051
其中，xi＝mi+Δxi，mi＝[xi]。也即mi是圆整得到正整数，而Δxi是被四舍五入后丢去的“随机”误差。
(8)式可以进一步利用等差级数求和公式得到：
Figure BDA0000704677310000052
(9)式中的前一项是量化后的值，虽然比(7)式的结果小了一半，但按照误差理论，一个数据的精度并不因乘以一个固定非零常数而改变。但后面一项中是零均值的随机数，相比(7)式中的要降低
Figure BDA0000704677310000053
倍，因此，对锯齿波或三角波激励信号进行过采样后同样可以得到提高精度的效果，且不需要另外加高频扰动信号。
本发明实施例对各器件的型号除做特殊说明的以外，其他器件的型号不做限制，只要能完成上述功能的器件均可。
实施例2
一种应用于乳腺的正交三角波频率编码高光谱图成像测量系统，参见图1和图2，该实施例以激光二极管作为单色光源LD1…LDn为例进行说明。
采用不同频率的正交三角波分别驱动一组单色光源1中的各个激光二极管LDi，摄像头4中每个像素点接收到每个激光二极管LDi透过乳腺组织样本3的单色光组合Iij；计算机对单色光组合Iij进行分离可以得到单色光组合Iij中的各个激光二极管LDi的贡献，据此可以对乳腺组织样本3进行透射高光谱图成像。
本实施例中的计算机处理步骤均与实施例1相同，本实施例在此不再赘述。
本发明实施例对各器件的型号除做特殊说明的以外，其他器件的型号不做限制，只要能完成上述功能的器件均可。
实施例3
一种应用于乳腺的正交三角波频率编码高光谱图成像测量系统，参见图1和图2，该实施例以单色二极管作为单色光源LD1…LDn为例进行说明。
采用不同频率的正交三角波分别驱动一组单色光源1中的各个单色二极管LDi，摄像头4中每个像素点接收到每个单色二极管LDi透过乳腺组织样本3的单色光组合Iij；计算机对单色光组合Iij进行分离可以得到单色光组合Iij中的各个单色二极管LDi的贡献，据此可以对乳腺组织样本3进行透射高光谱图成像。
本实施例中的计算机处理步骤均与实施例1相同，本实施例在此不再赘述。
本发明实施例对各器件的型号除做特殊说明的以外，其他器件的型号不做限制，只要能完成上述功能的器件均可。
实施例4
一种应用于乳腺的正交三角波频率编码高光谱图成像测量系统，参见图1和图2，该实施例以激光二极管作为单色光源LD1…LDn，手机摄像头作为摄像头为例进行说明。
采用不同频率的正交三角波分别驱动一组单色光源1中的各个单色二极管LDi，手机摄像头4中每个像素点接收到每个激光二极管LDi透过乳腺组织样本3的单色光组合Iij；计算机对单色光组合Iij进行分离可以得到单色光组合Iij中的各个激光二极管LDi的贡献，据此可以对乳腺组织样本3进行透射高光谱图成像。
本实施例中的计算机处理步骤均与实施例1相同，本实施例在此不再赘述。
本发明实施例对各器件的型号除做特殊说明的以外，其他器件的型号不做限制，只要能完成上述功能的器件均可。
实施例5
一种应用于乳腺的正交三角波频率编码高光谱图成像测量系统，参见图1和图2，该实施例以单色二极管作为单色光源LD1…LDn，手机摄像头作为摄像头为例进行说明。
采用不同频率的正交三角波分别驱动一组单色光源1中的各个单色二极管LDi，手机摄像头4中每个像素点接收到每个单色二极管LDi透过乳腺组织样本3的单色光组合Iij；计算机对单色光组合Iij进行分离可以得到单色光组合Iij中的各个激光二极管LDi的贡献，据此可以对乳腺组织样本3进行透射高光谱图成像。
本实施例中的计算机处理步骤均与实施例1相同，本实施例在此不再赘述。
本发明实施例对各器件的型号除做特殊说明的以外，其他器件的型号不做限制，只要能完成上述功能的器件均可。
本领域技术人员可以理解附图只是一个优选实施例的示意图，上述本发明实施例序号仅仅为了描述，不代表实施例的优劣。
以上所述仅为本发明的较佳实施例，并不用以限制本发明，凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

Claims (1)

1.一种应用于乳腺的正交三角波频率编码高光谱图成像测量系统，所述成像测量系统包括：一组单色光源、一个摄像头，以及与摄像头外接的计算机，其特征在于，一组单色光源分布在乳腺组织的一侧，摄像头分布在乳腺组织的另一侧；单色光源的波长处于600至1200nm；所述高光谱图成像测量系统可进入家庭使用，用于乳腺自检；
其中，一组单色光源中的各个单色光源密集排列在一个所设半球面上，采用透镜汇聚成一束光，构成光源；一个摄像头构成光源接收器件；
采用不同频率的正交三角波分别驱动一组单色光源中的各个单色光源，摄像头接收到的图像中每个像素点是每个单色光源透过乳腺的单色光组合；
计算机对单色光组合进行分离得到单色光组合中的各个单色光源的贡献，据此实现对乳腺的透射高光谱图成像；
所述计算机对单色光组合进行分离得到单色光组合中的各个单色光源的贡献，据此实现对乳腺的透射高光谱图成像的操作具体为：
包括4种不同波长的发光二极管，λ1和λ2波长的发光二极管驱动正交三角波频率为2倍f，λ3和λ4波长的发光二极管驱动正交三角波频率为1倍f，相同驱动频率正交三角波的相位相差90；
摄像头采样频率为fS，且fS＝2f，在λ1驱动信号由低变高的过程中采样；
在最低驱动信号频率的一个周期各路正交三角波信号的幅值和背景光信号的幅值不变，以顺序每8个数字信号为一组进行运算，分别得到4/3倍的波长λ1、λ2的光信号，2倍的波长λ3、λ4的光信号；
所述单色光源为激光二极管；或，所述单色光源为单色二极管；
所述单色光源、所述摄像头在乳腺组织样本两面对称设置；
所述摄像头为手机摄像头。
CN201510196490.1A 2015-04-23 2015-04-23 应用于乳腺的正交三角波频率编码高光谱图成像测量系统 Expired - Fee Related CN104799817B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510196490.1A CN104799817B (zh) 2015-04-23 2015-04-23 应用于乳腺的正交三角波频率编码高光谱图成像测量系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510196490.1A CN104799817B (zh) 2015-04-23 2015-04-23 应用于乳腺的正交三角波频率编码高光谱图成像测量系统

Publications (2)

Publication Number Publication Date
CN104799817A CN104799817A (zh) 2015-07-29
CN104799817B true CN104799817B (zh) 2021-08-27

Family

ID=53685343

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510196490.1A Expired - Fee Related CN104799817B (zh) 2015-04-23 2015-04-23 应用于乳腺的正交三角波频率编码高光谱图成像测量系统

Country Status (1)

Country Link
CN (1) CN104799817B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106404666A (zh) * 2016-10-20 2017-02-15 天津大学 抬高预设电平的正交三角波频率编码的高光谱成像系统

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102389313A (zh) * 2011-08-17 2012-03-28 天津大学 一种方波调制光电容积脉搏波测量装置和测量方法
CN102429646A (zh) * 2011-08-17 2012-05-02 天津大学 一种正交正弦波光电容积脉搏波测量装置和测量方法
CN102579012A (zh) * 2012-03-15 2012-07-18 天津大学 基于锁相光子计数测量的动态荧光-光学乳腺断层成像装置
CN203763052U (zh) * 2013-12-13 2014-08-13 深圳市瀚科创智医疗电子有限公司 一种乳腺组织检测仪及其检测探头

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102389313A (zh) * 2011-08-17 2012-03-28 天津大学 一种方波调制光电容积脉搏波测量装置和测量方法
CN102429646A (zh) * 2011-08-17 2012-05-02 天津大学 一种正交正弦波光电容积脉搏波测量装置和测量方法
CN102579012A (zh) * 2012-03-15 2012-07-18 天津大学 基于锁相光子计数测量的动态荧光-光学乳腺断层成像装置
CN203763052U (zh) * 2013-12-13 2014-08-13 深圳市瀚科创智医疗电子有限公司 一种乳腺组织检测仪及其检测探头

Also Published As

Publication number Publication date
CN104799817A (zh) 2015-07-29

Similar Documents

Publication Publication Date Title
US20150319347A1 (en) Methods and apparatus for imaging in scattering environments
US8809764B2 (en) Absolute position measuring encoder
CN107595297B (zh) 生理信息检测装置及方法
US10682081B2 (en) Acquiring and processing non-contact functional near-infrared spectroscopy data
KR102280072B1 (ko) Lidar 시스템들에 대한 간섭 검출 및 완화
TW201527784A (zh) 利用光脈衝壓縮的以飛行時間為基礎之系統
O’Connor et al. Underwater modulated pulse laser imaging system
Yu et al. Single-photon compressive imaging with some performance benefits over raster scanning
CN103860152A (zh) 一种对脉搏波信号进行处理的方法
WO2021092250A1 (en) Spatiotemporal antialiasing in photoacoustic computed tomography
CN104799817B (zh) 应用于乳腺的正交三角波频率编码高光谱图成像测量系统
CN104814721B (zh) 应用于乳腺的三角波频率编码的高光谱图成像测量系统
CN104783761B (zh) 应用于乳腺的正交方波频率编码的高光谱图成像测量系统
CN104921703B (zh) 应用于乳腺的方波频率编码的高光谱图成像测量系统
CN111045030B (zh) 一种深度测量装置和方法
KR101211050B1 (ko) 전역 광단층 영상기기의 이미지 해상도 향상방법
CN104799821B (zh) 应用于乳腺的正交正弦波频率编码高光谱图成像测量系统
CN105167741B (zh) 应用于乳腺的正弦波频率编码的高光谱图成像测量系统
Piao et al. Instrumentation for video-rate near-infrared diffuse optical tomography
EP3311138B1 (fr) Determination de coefficient d'absorption et de diffusion a l'aide d'un signal optique de relfexion calibre
US9503624B2 (en) Optical imaging system using multiple light sources and driving control method thereof
KR101322833B1 (ko) 이차원 실시간 공초점 영상에 기반한 티―레이 토모그래피 장치 및 방법
US8487274B2 (en) Stroboscopic optical image mapping system
Tsagkatakis et al. Active range imaging via random gating
Chiarelli et al. Flexible CW-fNIRS system based on Silicon Photomultipliers: In-Vivo characterization of sensorimotor response

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
EXSB Decision made by sipo to initiate substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20210827