CN103877677A - 一种新生儿黄疸治疗仪用led光源的光谱构造方法 - Google Patents
一种新生儿黄疸治疗仪用led光源的光谱构造方法 Download PDFInfo
- Publication number
- CN103877677A CN103877677A CN201410135960.9A CN201410135960A CN103877677A CN 103877677 A CN103877677 A CN 103877677A CN 201410135960 A CN201410135960 A CN 201410135960A CN 103877677 A CN103877677 A CN 103877677A
- Authority
- CN
- China
- Prior art keywords
- absorption spectrum
- spectrum
- monochromatic led
- light source
- several monochromatic
- 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.)
- Granted
Links
Images
Landscapes
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
本发明公开了一种新生儿黄疸治疗仪用LED光源的光谱构造方法,包括:通过实验测得标准胆红素溶液的吸收光谱,并构造该吸收光谱的数据组;选用若干种单色LED,并构造所述若干种单色LED的光谱数据矩阵;以所述若干种单色LED按一定的比例组合后能够模拟产生所述吸收光谱为目标,设定所述若干种单色LED的未知比例系数矩阵,并根据所述若干种单色LED的光谱数据矩阵、未知比例系数矩阵及所述吸收光谱的数据组构造超定方程组;采用简单遗传算法求解所述超定方程组的非负最小二乘解,得到若干种单色LED的最优比例组合及其在最优比例组合下的拟合光谱。本发明的算法具有简单易操作、效率高、匹配度高等优点。
Description
技术领域
本发明涉及光谱构造方法,特别涉及一种新生儿黄疸治疗仪用LED光源的光谱构造方法。
背景技术
随着光电技术的发展,发光二极管(Lighting-emitting diode,LED)发光效率不断提高,其单色性好、体积小、寿命长、能量消耗低、使用直流电、发热量低和控制灵活等特点,使LED光源在医疗、美容、植物栽培等领域的研究受到广泛关注,LED光源技术在生物医学上的研究和应用越来越广。自从1958年光照疗法用于治疗新生儿高胆红素血症(黄疸病)以后,应用于特殊部位或特殊病症治疗的特种波长、强度的光疗设备开始广泛被研制,其应用大幅降低了换血的需要。
而传统研究认为,游离胆红素对440nm波长的光或紫外线敏感;中国药典上描述:在400~500nm波长处,测定胆红素吸收曲线,其最大吸收为453nm;美国儿科医学会(AAP)定义的蓝光波长为430~490nm;也有研究表明,胆红素的吸收光谱在蓝色到绿色区域(460~490nm)。目前现有的LED黄疸治疗仪波长也很不一致,如BabyBlue(450~470nm)、NeoBlue(450~470nm)、Billitron(400~550nm)、BiliLED(453~487nm)等,因此,提出一种新型的新生儿黄疸治疗仪用LED光源光谱的形状具有重要的现实意义。
近年来,利用LED来模拟所需光源的目标光谱已成为研究的热点。应用LED的窄波段特性,根据人们的需要采用不同波段的LED模拟产生各种形状光谱分布的光谱匹配技术的研究及应用也显得尤为重要。目前进行光谱匹配的方法主要是利用光谱匹配算法和通过改变电流进行动态可调的光谱匹配技术,主要应用在太阳模拟技术、视觉检测和光谱图像处理等方面,在生物医学领域的应用研究还比较少,尤其是医学治疗领域,如新生儿黄疸治疗仪用LED光源光谱匹配方面的应用。
发明内容
本发明的目的在于克服现有技术之不足,提供一种新生儿黄疸治疗仪用LED光源的光谱构造方法,以实验测试的标准胆红素溶液的吸收光谱为目标光谱,基于智能算法中的简单遗传算法对LED光源光谱匹配进行优化组合,获得一种新型新生儿黄疸治疗仪用LED光源光谱的形状,大大简化了通过改变电流来调节光源光谱匹配问题的流程。
本发明解决其技术问题所采用的技术方案为:一种新生儿黄疸治疗仪用LED光源的光谱构造方法,包括如下步骤:
A1)通过实验测得标准胆红素溶液的吸收光谱,并构造该吸收光谱的数据组;
A2)根据所述吸收光谱的波长范围,选用不同峰值波长、峰值波长间隔、半高宽的若干种单色LED,并构造所述若干种单色LED的光谱数据矩阵;
A3)以所述若干种单色LED按一定的比例组合后能够模拟产生所述吸收光谱为目标,设定所述若干种单色LED的未知比例系数矩阵,并根据所述若干种单色LED的光谱数据矩阵、未知比例系数矩阵及所述吸收光谱的数据组构造超定方程组;
A4)采用简单遗传算法求解所述超定方程组的非负最小二乘解,得到所述若干种单色LED的最优比例组合;
A5)根据所述若干种单色LED的光谱数据矩阵及所述若干种单色LED的最优比例组合,构造所述若干种单色LED在该最优比例组合下所产生的拟合光谱。
进一步的,步骤A2)中,所述若干种单色LED的波长范围在350nm到800nm之间。
进一步的,步骤A4)中,所述简单遗传算法的求解步骤包括:
B1)随机产生十进制编码规模为100的初始群体;
B2)以最小二乘误差平方和为目标,建立用于评价群体优劣的适应度函数;
B3)根据所述适应度函数,反复对所述初始群体进行遗传运算,最后得到满意解或最优解。
进一步的,步骤A4)中,采用Matlab软件实现所述简单遗传算法。
进一步的,还包括:步骤A6)通过计算残差平方和与相关指数来评价所述若干种单色LED采用最优比例组合所产生的拟合光谱与所述吸收光谱之间的匹配结果。
由上述对本发明的描述可知,与现有技术相比,本发明具有如下有益效果:
本发明的一种新生儿黄疸治疗仪用LED光源的光谱构造方法,以实验测试的标准胆红素溶液的吸收光谱为目标光谱,基于智能算法中的简单遗传算法对LED光源光谱匹配进行优化组合,获得一种新型新生儿黄疸治疗仪用LED光源光谱的形状,大大简化了通过改变电流来调节光源光谱匹配问题的流程。该算法应用于光源光谱匹配技术简单易操作、效率高、误差小,拟合光谱和目标光谱十分近似,具有很高的光谱匹配度。
附图说明
图1为采用本发明的方法模拟标准胆红素溶液吸收光谱的匹配图。
具体实施方式
实施例:
本发明提供了一种新生儿黄疸治疗仪用LED光源的光谱构造方法,包括如下步骤:
A1)通过实验测得标准胆红素溶液的吸收光谱,并构造该吸收光谱的数据组,记为b=(y1,y2,…,ym)T;其中,测量标准胆红素溶液的吸收光谱的具体操作为:取纯度为99.28%的标准胆红素(Bilibubin)试剂(本实施例采用大连美仑生物技术有限公司的试剂)约20mg,精密称定;用少量三氯甲烷研磨后转移至100ml棕色量瓶中,超声处理使溶解;取出,迅速放冷,再加三氯甲烷稀释至刻度,摇匀;精密量取2.5ml,置100ml棕色瓶中,加三氯甲烷稀释至刻度,摇匀;用紫外可见分光光度计(U-3010,HITACHI)在常温18℃、避光条件下测定350-800nm波长吸收曲线,也就是所述标准胆红素溶液的吸收光谱;
A2)根据所述吸收光谱的波长范围,选用不同峰值波长、峰值波长间隔、半高宽的若干种单色LED,并构造所述若干种单色LED的光谱数据矩阵;假设Si(λ)是其中一种单色LED的光谱分布,记Si=(Si(λ1),Si(λ2),…,Si(λm))T,则所述若干种单色LED的光谱数据矩阵为A=(S1,S2,…,Sn);其中,单色LED的光谱分布模型为
A3)以所述若干种单色LED按一定的比例组合后能够模拟产生所述吸收光谱为目标,设定所述若干种单色LED的未知比例系数矩阵,记为X=(k1,k2,…,kn)T,ki为其中一种LED的未知比例系数,并根据所述若干种单色LED的光谱数据矩阵、未知比例系数矩阵及所述吸收光谱的数据组构造超定方程组;由于所述若干种单色LED的光谱数据矩阵为A=(S1,S2,…,Sn),其未知比例系数矩阵为X=(k1,k2,…,kn)T,所述吸收光谱的数据组为b=(y1,y2,…,ym)T,且根据光谱的叠加原理可知,所述若干种单色LED的光谱合成的基本数学模型为因此,当m>n时,可构造超定方程组为
AX=b (2)
其矩阵形式表示为
A4)采用简单遗传算法求解所述超定方程组的非负最小二乘解,得到所述若干种单色LED的最优比例组合;在实际应用中,因为LED的比例系数ki只能取非负值,因此只能取所述超定线性方程组的非负最小二乘解;
A5)根据所述若干种单色LED的光谱数据矩阵及所述若干种单色LED的最优比例组合,构造所述若干种单色LED在该最优比例组合下所产生的拟合光谱。
通常情况下,超定方程组是无解的,但它的近似解是可以求得的。超定方程组的最小二乘解是一种广义解,是指使残差r=b-AX的2-范数达到极小值时的解,因此它可以使得总体误差达较小。即:
取非负最小二乘解X*,则所述若干种单色LED在最优比例组合的情况下所产生的拟合光谱为
本实施例还包括:步骤A6)通过计算残差平方和与相关指数来评价所述若干种单色LED采用最优比例组合所产生的拟合光谱与所述吸收光谱之间的匹配结果。
残差平方和的定义为
相关指数的定义为
将和b代入公式(6)和(7),即可获取所述若干种单色LED采用最优比例组合所产生的拟合光谱(即拟合光谱)与所述标准胆红素溶液的吸收光谱(即目标光谱)的光谱匹配的残差平方和与相关指数。对于光谱匹配而言,相关指数越大,目标光谱和拟合光谱相似度就越高。
进一步的,步骤A2)中,所述若干种单色LED的波长范围在350nm到800nm之间。
进一步的,步骤A4)中,所述简单遗传算法的求解步骤包括:
B1)随机产生十进制编码规模为100的初始群体;
B2)以最小二乘误差平方和为目标,建立用于评价群体优劣的适应度函数;
B3)根据所述适应度函数,反复对所述初始群体进行遗传运算,最后得到满意解或最优解。
进一步的,步骤A4)中,采用Matlab软件实现所述简单遗传算法。
本实施例中所选取的若干种单色LED的峰值波长分别为400nm、425nm、450nm、470nm、490nm,它们的半高宽均为25nm。则根据本发明的方法对目标光谱(所述标准胆红素溶液的吸收光谱)进行匹配模拟,可计算得到所述若干种单色LED的最优比例组合为
0.2071:0.3763:0.4336:0.3299:0.0587,
取整数,可得所需五种波长的LED的数量分别约为
21、38、43、33、6;
所述五种波长的LED采用该最优比例组合即可模拟产生一种拟合光谱,用于新生儿黄疸治疗仪,该拟合光谱与目标光谱的匹配图如图1所示。该拟合光谱与目标光谱的残差平方和为RSS=0.0365、相关指数为R2=0.9839,可见,目标光谱和拟合光谱匹配度较高。
上述实施例仅用来进一步说明本发明的一种新生儿黄疸治疗仪用LED光源的光谱构造方法,但本发明并不局限于实施例,凡是依据本发明的技术实质对以上实施例所作的任何简单修改、等同变化与修饰,均落入本发明技术方案的保护范围内。
Claims (5)
1.一种新生儿黄疸治疗仪用LED光源的光谱构造方法,其特征在于,包括如下步骤:
A1)通过实验测得标准胆红素溶液的吸收光谱,并构造该吸收光谱的数据组;
A2)根据所述吸收光谱的波长范围,选用不同峰值波长、峰值波长间隔、半高宽的若干种单色LED,并构造所述若干种单色LED的光谱数据矩阵;
A3)以所述若干种单色LED按一定的比例组合后能够模拟产生所述吸收光谱为目标,设定所述若干种单色LED的未知比例系数矩阵,并根据所述若干种单色LED的光谱数据矩阵、未知比例系数矩阵及所述吸收光谱的数据组构造超定方程组;
A4)采用简单遗传算法求解所述超定方程组的非负最小二乘解,得到所述若干种单色LED的最优比例组合;
A5)根据所述若干种单色LED的光谱数据矩阵及所述若干种单色LED的最优比例组合,构造所述若干种单色LED在该最优比例组合下所产生的拟合光谱。
2.如权利要求1所述的一种新生儿黄疸治疗仪用LED光源的光谱构造方法,其特征在于,步骤A2)中,所述若干种单色LED的波长范围在350nm到800nm之间。
3.如权利要求1所述的一种新生儿黄疸治疗仪用LED光源的光谱构造方法,其特征在于,步骤A4)中,所述简单遗传算法的求解步骤包括:
B1)随机产生十进制编码规模为100的初始群体;
B2)以最小二乘误差平方和为目标,建立用于评价群体优劣的适应度函数;
B3)根据所述适应度函数,反复对所述初始群体进行遗传运算,最后得到满意解或最优解。
4.如权利要求1所述的一种新生儿黄疸治疗仪用LED光源的光谱构造方法,其特征在于,步骤A4)中,采用Matlab软件实现所述简单遗传算法。
5.如权利要求1所述的一种新生儿黄疸治疗仪用LED光源的光谱构造方法,其特征在于,还包括:步骤A6)通过计算残差平方和与相关指数来评价所述若干种单色LED采用最优比例组合所产生的拟合光谱与所述吸收光谱之间的匹配结果。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410135960.9A CN103877677B (zh) | 2014-04-04 | 2014-04-04 | 一种新生儿黄疸治疗仪用led光源的光谱构造方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410135960.9A CN103877677B (zh) | 2014-04-04 | 2014-04-04 | 一种新生儿黄疸治疗仪用led光源的光谱构造方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103877677A true CN103877677A (zh) | 2014-06-25 |
CN103877677B CN103877677B (zh) | 2017-01-04 |
Family
ID=50947005
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410135960.9A Active CN103877677B (zh) | 2014-04-04 | 2014-04-04 | 一种新生儿黄疸治疗仪用led光源的光谱构造方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103877677B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104991988A (zh) * | 2015-05-21 | 2015-10-21 | 大连工业大学 | 基于多颗单色大功率led实现类日光光源的方法 |
CN108644661A (zh) * | 2018-05-18 | 2018-10-12 | 杭州电子科技大学 | 利用多led的功率输出可调整的太阳光光谱合成方法 |
CN109442240A (zh) * | 2018-10-19 | 2019-03-08 | 南京理工大学 | 一种光谱可调模拟夜天光光谱的面光源 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20080089480A1 (en) * | 2006-10-16 | 2008-04-17 | Oraya Therapeutics, Inc. | Portable orthovoltage radiotherapy |
US20110117025A1 (en) * | 2008-05-20 | 2011-05-19 | Ralph Sebastian Dacosta | Device and method for fluorescence-based imaging and monitoring |
CN103251505A (zh) * | 2012-11-28 | 2013-08-21 | 苏州科医世凯半导体技术有限责任公司 | 一种led针灸治疗仪 |
-
2014
- 2014-04-04 CN CN201410135960.9A patent/CN103877677B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20080089480A1 (en) * | 2006-10-16 | 2008-04-17 | Oraya Therapeutics, Inc. | Portable orthovoltage radiotherapy |
US20110117025A1 (en) * | 2008-05-20 | 2011-05-19 | Ralph Sebastian Dacosta | Device and method for fluorescence-based imaging and monitoring |
CN103251505A (zh) * | 2012-11-28 | 2013-08-21 | 苏州科医世凯半导体技术有限责任公司 | 一种led针灸治疗仪 |
Non-Patent Citations (2)
Title |
---|
朱继亦等: "基于LED的光谱可调光源的光谱分布合成", 《发光学报》 * |
王福昌: "遗传算法在求解超定方程组中的应用", 《韶关学院学报》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104991988A (zh) * | 2015-05-21 | 2015-10-21 | 大连工业大学 | 基于多颗单色大功率led实现类日光光源的方法 |
CN104991988B (zh) * | 2015-05-21 | 2019-01-18 | 大连工业大学 | 基于多颗单色大功率led实现类日光光源的方法 |
CN108644661A (zh) * | 2018-05-18 | 2018-10-12 | 杭州电子科技大学 | 利用多led的功率输出可调整的太阳光光谱合成方法 |
CN109442240A (zh) * | 2018-10-19 | 2019-03-08 | 南京理工大学 | 一种光谱可调模拟夜天光光谱的面光源 |
Also Published As
Publication number | Publication date |
---|---|
CN103877677B (zh) | 2017-01-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Baranoski et al. | Light and skin interactions: simulations for computer graphics applications | |
Nishidate et al. | Estimation of melanin and hemoglobin using spectral reflectance images reconstructed from a digital RGB image by the Wiener estimation method | |
Archambault et al. | A mathematical formalism for hyperspectral, multipoint plastic scintillation detectors | |
Grossweiner et al. | The science of phototherapy: an introduction | |
US7835786B2 (en) | Methods, systems, and computer program products for optimization of probes for spectroscopic measurement in turbid media | |
CN106290244A (zh) | 实现炉膛内温度和气体浓度二维分布重建的系统及方法 | |
CN103877677A (zh) | 一种新生儿黄疸治疗仪用led光源的光谱构造方法 | |
CN104778349B (zh) | 一种用于水稻表土氮肥施用等级评定方法 | |
Ni et al. | Discrimination of Radix Isatidis and Rhizoma et Radix Baphicacanthis Cusia samples by near infrared spectroscopy with the aid of chemometrics | |
CN110414729A (zh) | 基于特征波长的植物潜在最大光合能力预测方法 | |
Olvera-Gonzalez et al. | A LED-based smart illumination system for studying plant growth | |
Gu et al. | Portable TDLAS sensor for online monitoring of CO2 and H2O using a miniaturized multi-pass cell | |
Hsu et al. | Developing low-cost mobile device and Apps for accurate skin spectrum measurement via low-cost spectrum sensors and deep neural network technology | |
Quistián‐Vázquez et al. | Retrieval of Absorption or Scattering Coefficient Spectrum (RASCS) program: a tool to monitor optical properties in real time | |
US20110178408A1 (en) | Accurate Low-Cost Non-Invasive Body Fat Measurement | |
Hegemann et al. | Using spectral sensors to determine photosynthetic photon flux density in daylight–A theoretical approach | |
Lisenko et al. | Method for estimating optimal spectral and energy parameters of laser irradiation in photodynamic therapy of biological tissue | |
CN105842203A (zh) | 基于玉米冠层特定光谱构建氮营养诊断植被指数 | |
Tarasov et al. | Combined non-invasive optical oximeter and flowmeter with basic metrological equipment | |
CN110487403A (zh) | 一种led光谱功率分布的预测方法 | |
Nedrebø | Use of market ready light dosimeters for patients with erythropoietic protoporphyria disorder | |
Dupouy et al. | Microbial absorption and backscattering coefficients from in situ and POLDER satellite data during an El Nino–Southern Oscillation cold phase in the equatorial Pacific (180) | |
Yang et al. | Method of matching different background spectra based on spectrally tunable light source | |
Dillon et al. | Determination of the rate coefficients for the reactions IO+ NO2+ M (air)→ IONO2+ M and O (3P)+ NO2→ O2+ NO using laser-induced fluorescence spectroscopy | |
Zhao et al. | High-precision fitting of plant growth spectra by multi-LED lighting systems based on segment fitting method |
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 |