技术领域
本发明涉及成像测量系统领域，尤其涉及一种应用于乳腺的方波频率编码的高光谱图 成像测量系统。
背景技术
现有技术中通过光对物体内部进行成像，特别是对人体内部进行成像，具有无损无创 无辐射的突出优点，但至今尚未有能够进入家庭使用的乳腺透射成像光测量系统用于经常 性的肿瘤自检，其原因在于现有乳腺成像测量系统的成本高、运算复杂、精度有限，无法 满足实际应用中的乳腺自检需要。
发明内容
本发明提供了一种应用于乳腺的方波频率编码的高光谱图成像测量系统，本发明实现 了高速度、大信息的乳腺透射光成像的高精度测量，满足了实际应用中的需要，详见下文 描述：
一种应用于乳腺的方波频率编码的高光谱图成像测量系统，所述成像测量系统包括： 一组单色光源、一个摄像头，以及与摄像头外接的计算机；一组单色光源分布在乳腺组织 样本的一侧，摄像头分布在乳腺组织样本的另一侧；
其中，一组单色光源中的各个单色光源密集排列在一个所设半球面上，采用透镜汇聚 成一束光，构成光源；摄像头构成光源接收器件；
采用不同频率且成2倍比率关系的方波分别驱动一组单色光源中的各个单色光源，摄 像头接收到的图像中每个像素点是每个单色光源透过乳腺的单色光组合；
计算机对单色光组合进行分离得到单色光组合中的各个单色光源的贡献，据此实现对 乳腺的透射高光谱图成像。
其中，所述单色光源、所述摄像头在所述乳腺组织样本两面对称设置。
所述计算机对单色光组合进行分离得到单色光组合中的各个单色光源的贡献，据此实 现对乳腺的透射高光谱图成像的操作具体为：
以4种波长的发光二极管2进行说明，λ1、λ2、λ3和λ4波长的发光二极管的驱动方 波频率分别为8f、4f、2f和f；
摄像头采样频率为fS，且fS＝2f，在λ1驱动信号高、低电平中间采样；
采样频率fS远高于调制方波信号和背景光的变化频率，在最低驱动信号频率的一个周 期各路方波信号的幅值和背景光信号的幅值不变；
以顺序每16个数字信号为一组进行运算，得到8倍的波长λ1、λ2、λ3和λ4的光信号。 其中，所述单色光源为激光二极管。
另一实施例，所述单色光源为单色二极管。
另一实施例，所述摄像头为手机摄像头。
本发明提供的技术方案的有益效果是：本发明采用不同频率且成2倍比率关系的方波 驱动单色光源，对摄像头检测到的光电信号进行分离可以得到单色光组合中的各个单色光 源的贡献，进而实现对乳腺组织样本的成像，本发明实现了高速度、大信息的线阵成像光 的高精度测量，且本发明具有成本低、应用方便等优点，适宜于家庭自检。
附图说明
图1为一种应用于乳腺的方波频率编码的高光谱图成像测量系统的结构示意图；
图2为本发明提供的单色光源、乳腺组织样本和摄像头相对位置示意图；
图3为方波的示意图。
附图中，各标号所代表的部件列表如下：
1：一组单色光源； 2：透镜；
3：乳腺组织样本； 4：摄像头。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚，下面对本发明实施方式作进一步地详 细描述。
实施例1
一种应用于乳腺的方波频率编码的高光谱图成像测量系统，参见图1和图2，成像测 量系统包括：一组n个单色光源1(用LD表示，单色光源的波长处于600～1200nm之间， 此波段为“光学窗口”，穿透深度比较深)和一个摄像头4。
其中，n的具体取值与乳腺组织样本3对不同波长的敏感程度相关，本发明实施例对 此不做限制。一组单色光源1分布在乳腺组织样本3的一侧，一个摄像头4分布在乳腺组 织样本3的另一侧；
其中，一组单色光源1中的每个单色光源LD1…LDn成线性排列，采用透镜2汇聚成 一束光，构成线阵光源；摄像头4构成多波长光源接收器件，摄像头4为CMOS或CCD， 现有摄像头4已能采集近红外波段的光。该成像测量系统还包括与摄像头4外接的计算机 (图中未示出)。
参见图3，采用不同频率且成2倍比率关系的方波分别驱动一组单色光源1中的各个 单色光源LDi，摄像头4进行采集，采集到的图像中，每个像素点为每个单色光源LDi透 过乳腺组织样本3的单色光组合Iij；计算机对单色光组合Iij进行分离可以得到单色光组合 Iij中的各个单色光源LDi的贡献，构成各波长下的透射图像。根据各波长下乳腺组织中正 常组织、肿瘤组织等的光学特性不同，进行乳腺的透射光成像测量，实现肿瘤的早期检测。
其中，计算机对单色光组合Iij进行分离可以得到单色光组合Iij中的各个单色光源LDi的贡献的步骤具体为：
为简便说明起见以4种波长的发光二极管2为例进行说明，假定λ1、λ2、λ3和λ4波 长的发光二极管的驱动方波频率分别为8f、4f、2f和f。
假定摄像头4采样频率为fS，且fS＝2f，并保证在λ1驱动信号高、低电平中间采样。
数字信号序列可以表示为：
D i t = D i λ 1 + D i λ 2 + D i λ 3 + D i λ 4 + D i B - - - ( 1 ) ]]>
其中，和分别为波长λ1、λ2、λ3和λ4的电压信号，为低频 信号，包括：背景光、摄像头4的本底噪声。
假定采样频率fS远高于调制方波信号和背景光的变化频率，在最低驱动信号频率的一 个周期可以近似认为各路方波信号的幅值和背景光信号的幅值不变。以最前16个采样数 据为例：
D 1 λ 1 = D 3 λ 1 = D 5 λ 1 = D 7 λ 1 = D 9 λ 1 = D 11 λ 1 = D 13 λ 1 = D 15 λ 1 = D A λ 1 D 2 λ 1 = D 4 λ 1 = D 6 λ 1 = D 8 λ 1 = D 10 λ 1 = D 12 λ 1 = D 14 λ 1 = D 16 λ 1 = 0 D 1 λ 2 = D 2 λ 2 = D 5 λ 2 = D 6 λ 2 = D 9 λ 2 = D 10 λ 2 = D 13 λ 2 = D 14 λ 2 = D A λ 2 D 3 λ 2 = D 4 λ 2 = D 7 λ 2 = D 8 λ 2 = D 11 λ 2 = D 12 λ 2 = D 15 λ 2 = D 16 λ 2 = 0 D 1 λ 3 = D 2 λ 3 = D 3 λ 3 = D 4 λ 3 = D 9 λ 3 = D 10 λ 3 = D 11 λ 3 = D 12 λ 3 = D A λ 3 D 5 λ 3 = D 6 λ 3 = D 7 λ 3 = D 8 λ 3 = D 13 λ 3 = D 14 λ 3 = D 15 λ 3 = D 16 λ 3 = 0 D 1 λ 4 = D 2 λ 4 = D 3 λ 4 = D 4 λ 4 = D 5 λ 4 = D 6 λ 4 = D 7 λ 4 = D 8 λ 4 = D A λ 4 D 9 λ 4 = D 10 λ 4 = D 11 λ 4 = D 12 λ 4 = D 13 λ 4 = D 14 λ 4 = D 15 λ 4 = D 16 λ 4 = 0 D 1 B = D 2 B = D 3 B = D 4 B = D 5 B = D 6 B = D 7 B = D 8 B = D 9 B = D 10 B = D 11 B = D 12 B = D 13 B = D 14 B = D 15 B = D 16 B = D A B - - - ( 2 ) ]]>
其中，和分别为波长λ1、λ2、λ3和λ4的光信号和背景信号 的幅值。
换言之，以顺序每16个数字信号为一组进行运算：
D 16 n + 1 - D 16 n + 2 + D 16 n + 3 - D 16 n + 4 + D 16 n + 5 - D 16 n + 6 + D 16 n + 7 - D 16 n + 8 + D 16 n + 9 - D 16 n + 10 + D 16 n + 11 - D 16 n + 12 + D 16 n + 13 - D 16 n + 14 + D 16 n + 15 - D 16 n + 16 = 8 D An λ 1 , n = 0,1,2 . . . . . . - - - ( 3 ) ]]>
即得到8倍的波长λ1的光信号而且完全消除了背景信号的影响。
D 16 n + 1 + D 16 n + 2 - D 16 n + 3 - D 16 n + 4 + D 16 n + 5 + D 16 n + 6 - D 16 n + 7 - D 16 n + 8 + D 16 n + 9 + D 16 n + 10 - D 16 n + 11 - D 16 n + 12 + D 16 n + 13 + D 16 n + 14 - D 16 n + 15 - D 16 n + 16 = 8 D An λ 2 , n = 0,1,2 . . . . . . - - - ( 4 ) ]]>
即得到8倍的波长λ2的光信号而且完全消除了背景信号的影响。
D 16 n + 1 + D 16 n + 2 + D 16 n + 3 + D 16 n + 4 - D 16 n + 5 - D 16 n + 6 - D 16 n + 7 - D 16 n + 8 + D 16 n + 9 + D 16 n + 10 + D 16 n + 11 + D 16 n + 12 - D 16 n + 13 - D 16 n + 14 - D 16 n + 15 - D 16 n + 16 = 8 D An λ 3 , n = 0,1,2 . . . . . . - - - ( 5 ) ]]>
即得到8倍的波长λ3的光信号而且完全消除了背景信号的影响。
D 16 n + 1 + D 16 n + 2 + D 16 n + 3 + D 16 n + 4 + D 16 n + 5 + D 16 n + 6 + D 16 n + 7 + D 16 n + 8 - D 16 n + 9 - D 16 n + 10 - D 16 n + 11 - D 16 n + 12 - D 16 n + 13 - D 16 n + 14 - D 16 n + 15 - D 16 n + 16 = 8 D An λ 4 , n = 0,1,2 . . . . . . - - - ( 6 ) ]]>
即得到8倍的波长λ4的光信号而且完全消除了背景信号的影响。
本发明实施例对各器件的型号除做特殊说明的以外，其他器件的型号不做限制，只要 能完成上述功能的器件均可。
实施例2
一种应用于乳腺的方波频率编码的高光谱图成像测量系统，参见图1和图2，该实施 例以激光二极管作为单色光源LD1…LDn为例进行说明。
采用不同频率且成2倍比率关系的方波分别驱动一组单色光源1中的各个激光二极管 LDi，摄像头4中每个像素点接收到每个激光二极管LDi透过乳腺组织样本3的单色光组合 Iij；计算机对单色光组合Iij进行分离可以得到单色光组合Iij中的各个激光二极管LDi的贡 献，据此可以对乳腺组织样本3进行透射高光谱图成像。
本实施例中的计算机处理步骤均与实施例1相同，本实施例在此不再赘述。
本发明实施例对各器件的型号除做特殊说明的以外，其他器件的型号不做限制，只要 能完成上述功能的器件均可。
实施例3
一种应用于乳腺的方波频率编码的高光谱图成像测量系统，参见图1和图2，该实施 例以单色二极管作为单色光源LD1…LDn为例进行说明。
采用不同频率且成2倍比率关系的方波分别驱动一组单色光源1中的各个单色二极管 LDi，摄像头4中每个像素点接收到每个单色二极管LDi透过乳腺组织样本3的单色光组合 Iij；计算机对单色光组合Iij进行分离可以得到单色光组合Iij中的各个单色二极管LDi的贡 献，据此可以对乳腺组织样本3进行透射高光谱图成像。
本实施例中的计算机处理步骤均与实施例1相同，本实施例在此不再赘述。
本发明实施例对各器件的型号除做特殊说明的以外，其他器件的型号不做限制，只要 能完成上述功能的器件均可。
实施例4
一种应用于乳腺的方波频率编码的高光谱图成像测量系统，参见图1和图2，该实施 例以激光二极管作为单色光源LD1…LDn，手机摄像头作为摄像头为例进行说明。
采用不同频率且成2倍比率关系的方波分别驱动一组单色光源1中的各个激光二极管 LDi，手机摄像头4中每个像素点接收到每个激光二极管LDi透过乳腺组织样本3的单色光 组合Iij；计算机对单色光组合Iij进行分离可以得到单色光组合Iij中的各个激光二极管LDi的贡献，据此可以对乳腺组织样本3进行透射高光谱图成像。
本实施例中的计算机处理步骤均与实施例1相同，本实施例在此不再赘述。
本发明实施例对各器件的型号除做特殊说明的以外，其他器件的型号不做限制，只要 能完成上述功能的器件均可。
实施例5
一种应用于乳腺的方波频率编码的高光谱图成像测量系统，参见图1和图2，该实施 例以单色二极管作为单色光源LD1…LDn，手机摄像头作为摄像头为例进行说明。
采用不同频率且成2倍比率关系的方波分别驱动一组单色光源1中的各个单色二极管 LDi，手机摄像头4中每个像素点接收到每个单色二极管LDi透过乳腺组织样本3的单色光 组合Iij；计算机对单色光组合Iij进行分离可以得到单色光组合Iij中的各个激光二极管LDi的贡献，据此可以对乳腺组织样本3进行透射高光谱图成像。
本实施例中的计算机处理步骤均与实施例1相同，本实施例在此不再赘述。
本发明实施例对各器件的型号除做特殊说明的以外，其他器件的型号不做限制，只要 能完成上述功能的器件均可。
本领域技术人员可以理解附图只是一个优选实施例的示意图，上述本发明实施例序号 仅仅为了描述，不代表实施例的优劣。
以上所述仅为本发明的较佳实施例，并不用以限制本发明，凡在本发明的精神和原则 之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。