CN1395683A - 用可见光线和/或近红外线诊断有无乳房炎的方法与装置 - Google Patents
用可见光线和/或近红外线诊断有无乳房炎的方法与装置 Download PDFInfo
- Publication number
- CN1395683A CN1395683A CN01800748A CN01800748A CN1395683A CN 1395683 A CN1395683 A CN 1395683A CN 01800748 A CN01800748 A CN 01800748A CN 01800748 A CN01800748 A CN 01800748A CN 1395683 A CN1395683 A CN 1395683A
- Authority
- CN
- China
- Prior art keywords
- ray
- milk
- near infrared
- raw milk
- light
- 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.)
- Pending
Links
- 208000004396 mastitis Diseases 0.000 title claims abstract description 41
- 238000000034 method Methods 0.000 title claims abstract description 37
- 235000020185 raw untreated milk Nutrition 0.000 claims abstract description 64
- 210000002700 urine Anatomy 0.000 claims abstract description 33
- 238000000611 regression analysis Methods 0.000 claims abstract description 16
- 230000005855 radiation Effects 0.000 claims abstract 4
- 235000013336 milk Nutrition 0.000 claims description 44
- 239000008267 milk Substances 0.000 claims description 44
- 210000004080 milk Anatomy 0.000 claims description 44
- 241000283690 Bos taurus Species 0.000 claims description 40
- 210000004907 gland Anatomy 0.000 claims description 37
- 230000006651 lactation Effects 0.000 claims description 37
- 238000003745 diagnosis Methods 0.000 claims description 23
- 238000005070 sampling Methods 0.000 claims description 22
- 230000005540 biological transmission Effects 0.000 claims description 17
- 238000002329 infrared spectrum Methods 0.000 claims description 16
- 230000003287 optical effect Effects 0.000 claims description 14
- 239000013307 optical fiber Substances 0.000 claims description 8
- 238000001514 detection method Methods 0.000 claims description 7
- 238000012417 linear regression Methods 0.000 claims description 4
- 230000008676 import Effects 0.000 claims description 2
- 238000005259 measurement Methods 0.000 abstract description 17
- 238000001228 spectrum Methods 0.000 abstract description 17
- 210000000481 breast Anatomy 0.000 abstract 3
- 210000001082 somatic cell Anatomy 0.000 abstract 1
- 230000000392 somatic effect Effects 0.000 description 41
- 239000011159 matrix material Substances 0.000 description 22
- 210000004027 cell Anatomy 0.000 description 16
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 5
- 238000000862 absorption spectrum Methods 0.000 description 4
- 238000009792 diffusion process Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 4
- 238000000491 multivariate analysis Methods 0.000 description 4
- 230000003595 spectral effect Effects 0.000 description 4
- 238000006243 chemical reaction Methods 0.000 description 3
- 239000012530 fluid Substances 0.000 description 3
- 239000000126 substance Substances 0.000 description 3
- 239000006059 cover glass Substances 0.000 description 2
- 238000013016 damping Methods 0.000 description 2
- 238000007405 data analysis Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000004043 dyeing Methods 0.000 description 2
- 230000004446 light reflex Effects 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 102000004169 proteins and genes Human genes 0.000 description 2
- 108090000623 proteins and genes Proteins 0.000 description 2
- 238000004611 spectroscopical analysis Methods 0.000 description 2
- YBNMDCCMCLUHBL-UHFFFAOYSA-N (2,5-dioxopyrrolidin-1-yl) 4-pyren-1-ylbutanoate Chemical compound C=1C=C(C2=C34)C=CC3=CC=CC4=CC=C2C=1CCCC(=O)ON1C(=O)CCC1=O YBNMDCCMCLUHBL-UHFFFAOYSA-N 0.000 description 1
- 229910001218 Gallium arsenide Inorganic materials 0.000 description 1
- 229910000530 Gallium indium arsenide Inorganic materials 0.000 description 1
- XUIMIQQOPSSXEZ-UHFFFAOYSA-N Silicon Chemical compound [Si] XUIMIQQOPSSXEZ-UHFFFAOYSA-N 0.000 description 1
- 238000010521 absorption reaction Methods 0.000 description 1
- 239000004411 aluminium Substances 0.000 description 1
- XAGFODPZIPBFFR-UHFFFAOYSA-N aluminium Chemical compound [Al] XAGFODPZIPBFFR-UHFFFAOYSA-N 0.000 description 1
- 229910052782 aluminium Inorganic materials 0.000 description 1
- 230000033228 biological regulation Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 230000031709 bromination Effects 0.000 description 1
- 238000005893 bromination reaction Methods 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 150000001720 carbohydrates Chemical class 0.000 description 1
- 235000014633 carbohydrates Nutrition 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 239000004020 conductor Substances 0.000 description 1
- 235000020247 cow milk Nutrition 0.000 description 1
- 239000007850 fluorescent dye Substances 0.000 description 1
- 238000005558 fluorometry Methods 0.000 description 1
- 230000004907 flux Effects 0.000 description 1
- 238000009499 grossing Methods 0.000 description 1
- 229910052736 halogen Inorganic materials 0.000 description 1
- 238000010438 heat treatment Methods 0.000 description 1
- 229910052500 inorganic mineral Inorganic materials 0.000 description 1
- 210000004962 mammalian cell Anatomy 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000000386 microscopy Methods 0.000 description 1
- 239000011707 mineral Substances 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 230000032696 parturition Effects 0.000 description 1
- 238000004321 preservation Methods 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 229910052710 silicon Inorganic materials 0.000 description 1
- 239000010703 silicon Substances 0.000 description 1
- 230000006641 stabilisation Effects 0.000 description 1
- 238000011105 stabilization Methods 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 239000012780 transparent material Substances 0.000 description 1
- 239000010937 tungsten Substances 0.000 description 1
- 229910052721 tungsten Inorganic materials 0.000 description 1
- -1 tungsten halogen Chemical class 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/25—Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
- G01N21/31—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry
- G01N21/35—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light
- G01N21/359—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light using near infrared light
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/25—Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
- G01N21/31—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/25—Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
- G01N21/31—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry
- G01N21/35—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light
- G01N21/3577—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light for analysing liquids, e.g. polluted water
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
- G01N33/02—Food
- G01N33/04—Dairy products
Landscapes
- Physics & Mathematics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Engineering & Computer Science (AREA)
- Food Science & Technology (AREA)
- Medicinal Chemistry (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
- Investigating Or Analysing Biological Materials (AREA)
Abstract
一种以光学手段通过对牛的尿液、生奶或乳房中体细胞计数而在短时间内以高准确度诊断奶牛乳房炎的方法和装置。装置包括:发射测量光的光源(1)、通过使光色散仅从光源的光中提取所需波长的分光仪(9)设置在透镜(11)与反光镜(12)之间的光斩波器(14)、由光漫射元件形成的积分球(13)、和固定诸如生奶的样品的样品固定器(40)。将可见光或近红外辐射射在奶牛的尿液、生奶或乳房上;从来自生奶或乳房的光谱测量预定波长的透射光、反射光、或透射与反射光的强度;作多元回归分析;由此诊断乳房炎。
Description
发明的技术领域
本发明涉及根据来自奶牛的尿液、生奶或哺乳腺的可见光和/或红外光谱来诊断乳房炎的方法与装置。
现有技术
生奶中体细胞的数量是乳房炎诊断的主要标志。迄今为止,测量体细胞数量一直是应用直接显微法、CMT修正法和荧光测定法。
目前,一般使用荧光型体细胞计数器(Fossomatic)测量生奶的体细胞数。该装置能计算和显示每毫升的体细胞数,方法是将缓冲液与染色液(溴化3,8-二氨基-5-乙基-6-苯基菲啶鎓溶液)混入生奶,对体细胞的细胞核作荧光染色,将得到的混合奶用微注射器分散加到连续旋转盘的外围部分,再用荧光显微镜自动测量体细胞数。
日本规定,在用荧光型体细胞订数器(Fossomatic)测量生奶时,若每毫升的体细胞数达300000或更多,就断定该奶牛已患乳房炎,禁止挤奶。
本发明要解决的问题
然而,基于用荧光型体细胞计数器测量体细胞的常规乳房炎诊断法有各种问题要解决:(1)挤取的生奶要加缓冲液与染色液等化学品作预处理;(2)不能以非破坏方式测量生奶样品;(3)生奶易被其它物质干扰;(4)化学品价格高,从成本性能的观点来看是不利的;和(5)操纵装置与样品要用熟练的技术方法。
因此,本发明的目的是提出一种测量方法与装置以及判断方法,它通过对奶牛的尿液、生奶或哺乳腺发出的可见光和/或近红外光谱作光学测量,能在短时间内高度精密地诊断乳房炎。
解决问题的措施
本发明涉及诊断奶牛乳房炎的方法,包括以下步骤:将400~2500nm波长范围内的可见光线和/或近红外线辐照入奶牛的尿液、生奶或哺乳腺里,检测来自所述尿液、生奶或哺乳腺的透射光线、反射光线或者透射与反射光线的强度,作多元回归分析,诊断奶牛存在乳房炎。相对于用于检测的可见光线和/或近红外线,选用能有效判断乳房炎的波长的光线。
根据本发明,检测来自奶牛的尿液、生奶或哺乳腺的透射光线、反射光线或者透射与反射光线的强度,可以确定随尿液、生奶或哺乳腺中体细胞数变化的吸收比,这样通过多元回归分析和诊断存在奶牛乳房炎,可以判定奶牛的乳房炎。因此,无须作常规而麻烦的预处理、使用贵重化学品和熟练操作样品。从哺乳腺反射的光线等的强度表示由活体组织(包括哺乳动物细胞)反射的光线等的强度。在哺乳腺中,哺乳腺细胞(包括生奶)和活组织分别认为是牛乳和透明小容器。至于透射光线的测量,例如,将入射光线经光纤加到哺乳腺的一个右侧,同时把该光纤加到左侧,通过加到哺乳腺另一左侧的另一根光纤测量透射光线(在检测器一侧)。近红外范围的光线在哺乳腺中通过甚至更厚的路径,视波长范围而定。
本发明中用于诊断乳房炎的可见光线与近红外线的波长范围为400~2500nm。若应用400~1000nm范围的可见光线与红外线,就使用硅光探测器,若应用700~2500nm波长范围的近红外线,PbSe、InGaAs或GaAs等光探测器。
在上述可见光线与近红外线波长范围内,由于400~700nm范围的可见光线与近红外线有噪声,所以最好使用700~2500nm范围的近红外线。再者,由于生奶含各种成分,如水、蛋白质、脂肪、碳水化合物、矿物质等,在各种波长区内,光线主要被作为主成分的水吸收,担心会妨碍测量近红外光谱。然而,在700~2500nm波长区内,与水相关的干扰比其它波长区要小。在1100~2500纳以长区内,尿液、生奶或哺乳腺中体细胞的吸收比变化作为分子振动的第一谐音或组合谐音。因此,最好用1100~2500nm波长范围的近红外线作测量,它能在短时间里测量尿液、生奶或哺乳腺中的体细胞。
再者,在红外线波长区中,由于在尿液、生奶或哺乳腺中吸收比较小,在测量透射光线或者透射与反射光线时,可将样品厚度保持为几毫米,因而易于搬动与设定样品容器。
通过对尿液、生奶或哺乳腺作光学测量和基于多元回归的数据处理,很容易高度精密地诊断奶牛的乳房炎。下面将描述这种诊断奶牛乳房炎的方法。
本发明的特征在于,在光学测量中,用分光仪在波长上扫描来自尿液、生奶或哺乳腺的入射光线、透射光线、反射光线或透射与反射光线,对得到的分光光谱作多元回归分析。
根据本发明,由于用分光仪扫描波长可得到波长分辨率高的基本上连续的分光光谱,故能得到数据分析必需的大量数据。例如,若在1100~2500nm波长区内以2nm波长分辨率作扫描,每次扫描可得到701个数据点,提高了数据分析的精密度。
本发明还涉及诊断奶牛乳房炎的装置,包括:(1)近红外线发生器,产生波长范围为400~2500nm的可见光线和/或近红外线;(2)将可见光线和/或近红外线引入奶牛的尿液、生奶或哺乳腺的光学系统;(3)检测器,检测来自所述尿液、生奶或哺乳腺的透射光线、反射光线或透射与反射光线的强度;和(4)数据处理器,用于接收来自所述检测器的信号,作多元回归分析,以诊断奶牛存在乳房炎。
本发明的乳房炎诊断装置最好还包括一光纤,用于将来自奶牛的所述尿液、生奶或哺乳腺的可见光线和/或近红外线引入光探测器,从而用所述光探测器通过该光纤检测与来自所述尿液、生奶或哺乳腺的透射光线、反射光线或透射与反射光线的强度。
利用光纤可提供便携小型的乳房炎诊断装置。
奶牛乳房炎诊断装置最好还包括一供料器,将所述生奶经在线或线上引入生奶样品容器。
设置供料器将所述生奶经在线或线上引入样品容器,可按消逝的时间连续地测量可见光线和/或近红外线。
奶牛乳房炎诊断装置还包括保存生奶的样品容器和使样品容器内的生奶稳定到指定温度的温控器。
样品容器中生奶温度稳定可防止生奶的吸收比因温度而变化,可提高诊断乳房炎的精度。在测量哺乳腺以诊断乳房炎时,将哺乳腺固定于挤奶机,必要时,按上述方法控制温度。
附图简介
图1是本发明一实施例的生奶光谱测量装置的结构图。
图2是一例生奶光谱测量装置的电气结构的框图。
图3是样品架40的剖视图。
图4是一曲线图,示出一例若干生奶样品在400~1100nm波长范围内的近红外线光谱。
图5是一曲线图,同一例若干生奶样品在1100~2500nm波长范围内的近红外线光谱。
图6示出执行PLS法的步骤,该方法是本发明多元回归分析法的一种形式。
图7给出的曲线表示用于各计算的若干PLS因素和分析误差。
图8是一曲线图,表示PLS法推定的体细胞数与使用取自奶牛的生奶样品时实际体细胞数之间的相关性。
图9是一曲线图,表示PLS法推定的体细胞数与使用取自奶牛的尿样的实际体细胞数之间的相关性。
实施本发明的实施例
现参照图1描述本发明的乳房炎诊断装置。
图1是本发明一实施例的生奶光谱测量装置的结构图。以光行进方向看,该装置包括产生测量光线的光源、使光源1发出的光线相互平行的透镜2,通过分离光源1的光线而取得所需光线的分光仪9、截断分光仪发射的高光部分光线的滤光器10、汇集分离光线的透镜11、插在透镜11与反射镜12之间的光斩波器14、由光漫射材料形成的积分球13、保存样品的样品架40等。
光源1由钨卤素灯等组成,产生波长范围很宽的光线,包括近红外线。分光仪9包括汇集入射光线的透镜3、调节光线通量与大小的缝隙4、对通过缝隙4的光线作反射的反射镜5、带曲线表面的衍射光栅6、控制衍射光栅6的衍射角的电机6a、在衍射光栅6上衍射的光线中只通过所需光部分的缝隙7、相互平行地发射衍射光线的光发射透镜8等。通过用电机6a作角度控制,能有选择地只取出所需波长的光线。
光斩波器14设计成转盘形式,其中交替安置了光反射部分与光透过部分,驱动电机14a而转动光学斩波器14,周期性反射或通过来自透镜11的光线。斩波传感器14b检测光斩波器14的旋到相位,根据来自斩波传感器14b的信号,同步电路14c输出同步信号Sa与Sb,分别表示来自透镜11的光线的反射与通过状态。
积分球13包括向上打开的入射光窗口13a、向下打开的光发射窗口13b和多个把接收的光线量转换为电信号的光探测器20。积分球13的作用是使进入该球的光线发生漫射而减少测量误差。光探测器20由PbS等构成,在近红外线区具有灵敏度。样品架40位于光发射窗口13b附近。
若用光斩波器14反射经分光仪9分离的光线,则进入样品架40的光线就像经入射光窗口13a通过积分球13一样,结果返回的光线在积分球13中漫射,一部分光线被光探测器20接收。另一方面,若分光仪9分离的光线通过光斩波器14,它被反射镜12反射而通过入射光窗口13a斜向进入积分球13,因而该光线漫射而不到达样品,部分光线被光探测器20接收。斩波器的上述操作取出了受样品影响和不受样品影响的信号。
图2是一例框图,表示生奶光谱测量装置中数据处理器的电气结构。来自光探测器20的检测信号被放大器21放大后,输入采样保持电路22以同步信号Sa采样,并输入采样保持电路23以同步信号Sb采样。采样保持电路22在光线由分光仪9进入样品时只在采样周期内保持信号电压,而采样保持电路23在光线未从分光仪9进入样品时只在采样周期内保持电信号。然后,从采样保持电路22与23输出的信号分别由对数转换电路24、25作对数转换,在减法电路26中将它们相减。通过与光斩波器14同步的检测,可以消除扰动分量。
减法电路的输出信号经AD(模/数)转换器27量化,导入个人计算机(PC)30,机内安装的各种程序可按多元回归法作数据处理。与PC30相接的有输入数据的键盘38、显示数据的显示器29等。
图3是一剖视图,表示样品架40的结构。该样品架40适合积分球13的光发射窗口13b的结构。由铝等热传导材料制作的样品架包括:样品容器41,可保持生奶等液体样品SP;透明盖玻璃板42,用于遮盖样品容器41的开口;对样品容器41加热或冷动的珀耳帖元件;温度传感器45,可控制样品容器41的温度;温度电路44,通过根据来自温度传感器45等的温度信号驱动珀耳帖元件,稳定样品SP的温度。
当光斩波器14反射的光线经盖玻璃板42进入样品SP时,它们在经过衰减与散射后(取决于样品SP的吸收光谱)至新返回积分球13,因此光探测器210接收一部分返回的光线并把它们转换成电信号。
由于生奶的吸收比对温度变化很敏感,而且必须得到生奶中很少受影响的脂肪,若测量环境温度在每次测量时都发生变化,就失去了测量的意义。因此,根据该实施例,用温度传感器45、温控电路44与珀耳帖元件43构成的温度反馈系统来稳定样品SP的温度,从而提高测量精密度。
图4是表示一例生奶近红外光谱的曲线图,其中纵坐标表示由指数代表的吸收比,指数是对光反射比的倒数作对数转换得到的,横坐标表示波长(nm)。曲线对应的吸收光谱,通过用图1的分光仪9在400~1100nm波长内扫描而得到。图4中,以重迭状态显示通过测量多个生奶样品而得到的结果。图5也是表示一例生奶近红外光谱的曲线图,其中纵坐标给出由指数表示的吸收比,指数通过对光反射比的例数作对数到换而得到,横坐标表示波长(nm)。曲线对应的吸收光谱是用分光仪9在1100~2500nm波长上扫描得到的。图5中,显示并重迭了测量多个生奶样品所得到的结果。
所有的曲线都可归因于水的吸收光谱,尤其是接近1400~1500nm和1850~2050的大的峰值可归因于水的分子振动。
上述说明是对透射与反射型结构作的,被测光线通过样品SP,在样品容器41的内表面上反射,再次通过样品SP。此外,也可用透射型测量,样品容器41由透明材料制作,检测通过样品SP的透射光,或采用反射型,测量从样品SP表面反射的光线。
上述说明对这样的结构例子作的,其中将分光仪9安置在光源1与样品SP之间,分离进入样品SP的光线。此外,可以应用结构例子,其中将分光仪9安置在样品SP与光探测器20之间,并分开从样品SP透射的光线或透射与反射光线。
下面说明多元回归分析。
根据多元回归法,建立某种样品的定量特性(如对应于本申请中体细胞浓度或数量)与一个或多个示例变量之间的关系。在本申请中,示例变量表示测得的近红外光谱。然而,近红外光谱实际上在波长下不是由单一吸收比而是由若干吸收比构成的。因此,由于多元数据称为多元回归分析,在本发明方法中,用近红外光谱从完全不同的测量值里推定已知样品的定量特性(本申请中为体细胞数)。在本申请中,基于将近红外光谱作为多元数据,对要推定的特性的体细胞数,用多元回归分析法找出相互关系。
作为多元回归分析法之一,下面将说明PLS法(局部最小平方回归)(参照:Tetsuro Aijima,“Chemometrics”1991年Marugen出版,和H.Martens &T.Naes,“多元校正法”,John Willy出版,纽约,1991)。如在图4中,若在700~1100nm波长范围内以2nm分辨率扫描,一次扫描可得201数据。
作为光谱预处理,必要时可作Kubelka-Munk转换(把视在反射指数转换成吸收率K(λ)与散射指数S(λ)之间的比率)、移动平均等平滑处理、一维或二维差分处理、基线校正处理等。
图6表示本发明PLS法的执行步骤。一个波长点的吸收比是一维数据,由多个分量构成(牛奶中的体细胞、蛋白质、类脂体等或哺乳腺的构成组分),一个光谱包括在M个波长点对波长扫描一次的多维数据,因而在体细胞的多元分析方法中,一个光谱数据构成了一维数据×M个直角坐标空间。
接着,运用正交型PLS法,以相互正交的系数记分t1,t2,t3,……ta,……,th代表光谱数据X中的变量和体细胞数里的变量。图6示出形成第a系数的系数参量的步骤。E是光谱残余矩阵(波长数据点×采样的样品数据阵),t为系数记分,w为加载权重(Loading weight),p与q分别为x与y加载。
a-1系数的谱残余矩陈Ea-1表示为系数记分ta与第a系数的加载权重wa的乘积、系数记分ta与X加载的乘积和第a系的谱残余矩阵Ea之和。第(a-1)系数的体细胞误差矩阵fa-1表示为系数记分ta与Y加载qa的乘积和第a系数的体细胞数误差矩阵fa之和。
因此,作为分析的一个步骤,可根据谱残余矩阵Ea-1和第(a-1)系数的体细胞数误差矩阵fa-1,来计算第(a-1)系数的系数记分、加载权重wa、X加载pa与Y加载qa,再计算谱残余矩阵Ea与体细胞数误差矩阵fa。
更具体地说,计算出该谱数据的第一残余矩阵E0=X-Xm(下标m表示“平均值”)和体细胞数的第一误差矩阵f0=Y-Ym,并算出第一系数的参量t1、w1、p1与q1。然后去除第一系数的影响,形成新的谱残余矩阵E1和体细胞误差矩阵f1。接着,运用谱残余矩阵E1和体细胞数误差矩阵f1,计算第二系数的参量t2、w2、p2与q2。然后去除第二系数的影响,形成新的谱残余矩阵E2与体细胞数误差矩阵f2,于是应用谱残余矩阵E2与体细胞数误差矩阵f2,计算第三系数的参量t3、w3、p3与q3。接着去除第三系数的影响,形成新的谱残余矩阵E3与体细胞数误差矩阵f3。以后以同样方法重复计算各系数的参量并形成每个残余矩阵与误差矩阵,直到提取完有效的系数。
图7的曲线示出了若干PLS系数与分析误差,该曲线图与体细胞数相关,纵坐标表示体细胞的对数。图7(a)给出的数据(SEC:标准校正误差)是在样品具有已知体细胞数时通过测量体细胞数而得到的,而图7(b)给出的数据(SEP:标准预测误差)是在使用体细胞数未知的样品时通过测量体细胞数而得到的。
从这些曲线可看出,当PLS系数为5时,即当误差计算重复到第五系数时,SEC与SEP两值都收敛在约0.25的误差范围内。
当使用图7中达到第五的PLS系数时,可构建一模型,它可以根据近红外谱逆向地计算体细胞数。因此,对该模型应用与未知的体细胞数相关的近红外谱,可算出某样品推定的体细胞数。
图8是一曲线图,表示按PLS法推定的体细胞数与实际体细胞数的相关性。横坐标表示实际的体细胞数,纵坐标表示推定的体细胞数。从该图可看出,在近乎线性分布的两种情况之间,图中有很强的相关性,说明按PLS法作多元分析极为有效。
下面描述MLR法与PCR法。根据MLR法,对一已知样品相对于测量得到的体细胞数编制一线性回归公式,在由吸收比组成的新红外谱中,只把多个特定波长的吸收比作为多个变量。将要预测的未知体细胞数的近红外谱应用于这样编制的线性回归公式,可预测出来知的体细胞数。在PCR法中,作基本上与PLS法一样的计算(参照:Yoshikatu Miyashita&Sinichi Sasaki,“Kemometrix-化学模式的识别与多元分析”,1995年由Kyoritsu出版者出版)。
图9是一曲线图,表示按PLS法推定的体细胞数与使用奶牛尿液时实际的体细胞数的相关性。尿液光谱应用上述生奶光谱一样的测量装置与同样的测量方法确定。由图9可知,横坐标的实际体细胞数与纵坐标的推定体细胞数的强相关性,呈几乎线性的方式分布,表明PS法的多元分析极其有效。
发明效果
如上所述,根据本发明,通过检测来自尿液、生奶或哺乳腺的透射光线、反射光线或透射与反射光线的强度,对得到的吸收比作多元回归分析,可以测量体细胞数。
通过对尿液作光学测量并对生奶或哺乳腺作数据处理,可高度精密地测量出体细胞数。
Claims (11)
1、一种诊断奶牛乳房炎的方法,其特征在于包括以下步骤:将波长范围为400~2500nm的可见光线和/或近红外光线辐射入奶牛的尿液、生奶或哺乳腺里,检测来自所述尿液、生奶或哺乳腺的透射光线、反射光线或透射与反射光线的强度,作多元回归,并诊断奶牛存在乳房炎。
2、如权利要求1的诊断奶牛乳房炎的方法,其特征在于包括以下步骤:将波长范围为700~2500nm的近红外线辐射入奶牛的尿液、生奶或哺乳腺里,检测来自所述尿液、生奶或哺乳腺的透射光线、反射光线或透射与反射光线的强度,作多元回归,并诊断奶牛存在乳房炎。
3、一种诊断奶牛乳房炎的方法,其特征在于包括以下步骤:将波长范围为700~1100nm或1100~2500nm的近红外线辐射入奶牛的尿液、生奶或哺乳腺里,检测来自所述尿液、生奶或哺乳腺透射光线、反射光线或透射与反射光线的强度,作多元回归分析,并诊断奶牛存在乳房炎。
4、如权利要求1~3之一的奶牛乳房炎诊断法,其特征在于包括,对来自所述尿液、生奶或哺乳腺的入射光线、透射光线、反射光线或透射与反射光线的波长作扫描,并对如此获得的可见光线和/或近红外线的光谱作多元回归分析。
5、如权利要求1~4之一的奶牛乳房炎诊断法,其中多元回归分析由多重线性回归(MLR)、主分量回归(PCR)或局部最小平方回归(PLS)执行。
6、一种诊断奶牛乳房炎的装置,包括:
近红外线发生器,产生波长范围为400~2500nm的可见光线和/或近红外线;
光学系统,将可见光线和/或近红外线引入奶牛的尿液、生奶或哺乳腺;
检测器,检测来自所述尿液、生奶或哺乳腺的透射光线、反射光线或透射与反射光线的强度;和
数据处理器,接收来自所述检测器的信号,作多元回归分析,以诊断奶牛存在乳房炎。
7、如权利要求6的奶牛乳房炎诊断装置,其中所述近红外线发生器是产生波长范围为700~2500nm的近红外线的红外线发生器、产生波长范围为700~1100nm的近红外线的红外线发生器,或产生波长范围为1100~2500nm的近红外线的红外线发生器。
8、如权利要求6或7的奶牛乳房炎诊断装置,其特征在于还包括一根接至所述近红外线发生器的光纤,用于将可见光线和/或近红外线引入所述的奶牛的尿液、生奶或哺乳腺,并将来自所述的尿液、生奶或哺乳腺的透射光线、反射光线或透射与反射光线通过该光纤导入所述检测器。
9、如权利要求6~8之一的奶牛乳房炎诊断装置,其特征在于还包括将所述生奶经在线或线上引入样品容器的供料器。
10、如权利要求6~9之一的奶牛乳房炎诊断装置,其特征在于还包括保存生奶的样品容器和使样品容器里的生奶稳定至指定温度的温控器。
11、如权利要求6~10之一的奶牛乳房炎诊断装置,其中多元回归分析由多重线性回归(MLR)、主分量回归(PCR)或局部最小平方回归(PLS)执行。
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2000099142 | 2000-03-31 | ||
JP99142/2000 | 2000-03-31 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN1395683A true CN1395683A (zh) | 2003-02-05 |
Family
ID=18613532
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN01800748A Pending CN1395683A (zh) | 2000-03-31 | 2001-03-14 | 用可见光线和/或近红外线诊断有无乳房炎的方法与装置 |
Country Status (14)
Country | Link |
---|---|
US (1) | US6748251B2 (zh) |
EP (1) | EP1191326A4 (zh) |
JP (1) | JP3574851B2 (zh) |
CN (1) | CN1395683A (zh) |
AU (1) | AU769362B2 (zh) |
CA (1) | CA2376173A1 (zh) |
CZ (1) | CZ20014305A3 (zh) |
HU (1) | HUP0201283A2 (zh) |
IL (1) | IL146786A (zh) |
NO (1) | NO20015820L (zh) |
NZ (1) | NZ515725A (zh) |
PL (1) | PL351030A1 (zh) |
RU (1) | RU2248554C2 (zh) |
WO (1) | WO2001075421A1 (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101918830A (zh) * | 2007-05-31 | 2010-12-15 | S.A.E.阿菲金公司 | 用于分析流体的系统和方法 |
CN102384889A (zh) * | 2011-08-03 | 2012-03-21 | 长春迪瑞医疗科技股份有限公司 | 一种全自动便潜血分析仪的积分球式光学系统 |
CN108027848A (zh) * | 2015-07-16 | 2018-05-11 | 皇家飞利浦有限公司 | 用于管理对人的炎症性自身免疫性疾病的处置的设备、系统和方法 |
Families Citing this family (26)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
FR2823565B1 (fr) * | 2001-04-12 | 2003-10-24 | Agronomique Inst Nat Rech | Procede de tracabilite d'une alimentation avec de l'herbe fraiche chez un animal |
WO2004017066A1 (en) * | 2002-08-16 | 2004-02-26 | Lattec I/S | A system and a method for observing and predicting physiological state of an animal |
US7302349B2 (en) | 2002-08-16 | 2007-11-27 | Lattec I/S | System and a method for observing and predicting a physiological state of an animal |
WO2006085897A2 (en) | 2004-05-13 | 2006-08-17 | Advanced Animal Diagnostics | Microfluidic device and leucocyte antigen mediated microfluidic assay |
US20070020357A1 (en) | 2005-07-20 | 2007-01-25 | Yissum Research Development Company | Method and device for determining milk quality |
JP5306651B2 (ja) | 2005-09-02 | 2013-10-02 | ポーラ化成工業株式会社 | 皮膚の状態の判定方法及び皮膚の厚みの予測方法 |
IL173604A (en) | 2006-02-08 | 2013-01-31 | E Afikim Milking Systems Agricultural Cooperative Ltd Sa | A device and method for recording animal poses, especially for live animals |
US20090233329A1 (en) | 2006-03-24 | 2009-09-17 | Rodriguez Rodolfo R | Microfluidic chamber assembly for mastitis assay |
WO2009111194A1 (en) * | 2008-02-29 | 2009-09-11 | Arkema Inc. | High efficiency photovoltaic modules |
EP2416644B1 (en) * | 2009-04-09 | 2015-06-03 | DeLaval Holding AB | Milking system |
US20120115184A1 (en) * | 2009-06-09 | 2012-05-10 | University of Tarty | Method for the detection of mastitis and milk quality and mastitis sensor |
EP2579705B2 (en) * | 2010-06-14 | 2022-08-10 | Milfos International Limited | Improved milking apparatus and system |
RU2485914C2 (ru) * | 2011-05-10 | 2013-06-27 | Государственное образовательное учреждение высшего профессионального образования "Московский государственный университет технологий и управления" | Способ ранней диагностики мастита у коров |
KR20130094083A (ko) * | 2012-02-15 | 2013-08-23 | 삼성전기주식회사 | 레이저 스캔 장치 |
US9297749B2 (en) * | 2012-03-27 | 2016-03-29 | Innovative Science Tools, Inc. | Optical analyzer for identification of materials using transmission spectroscopy |
JP2015057591A (ja) * | 2013-08-09 | 2015-03-26 | 株式会社島津製作所 | 懸濁液中の懸濁物濃度の分析方法および分析装置 |
ES2431242B1 (es) * | 2013-09-02 | 2014-09-09 | Centro De Investigación Y Tecnología Agroalimentaria De Aragón (Cita) | Procedimiento de identificación de animales gestantes |
CZ2013773A3 (cs) | 2013-10-03 | 2015-04-15 | Labmediaservis S.R.O. | Testovací sada a způsob roztřídění stáda krav |
US9574993B2 (en) * | 2014-06-24 | 2017-02-21 | Shimadzu Corporation | Method and apparatus for analyzing the concentration of materials in suspension |
CN104568794B (zh) * | 2015-01-01 | 2017-06-23 | 浙江工商大学 | 带鱼肉储存时间检测方法 |
CN104568796B (zh) * | 2015-01-01 | 2017-06-23 | 浙江工商大学 | 猪肉保存时间检测方法 |
CN105342574A (zh) * | 2015-12-11 | 2016-02-24 | 无限极(中国)有限公司 | 一种自动寻肌筋膜扳机点光学仪 |
US10371628B2 (en) * | 2017-08-07 | 2019-08-06 | The Boeing Company | Apparatus for measuring spectral hemispherical reflectance of samples at grazing angles |
JP7405555B2 (ja) * | 2019-10-08 | 2023-12-26 | 株式会社アドバンテスト | 体細胞計、体細胞測定方法、プログラムおよび記録媒体 |
CN111220565B (zh) * | 2020-01-16 | 2022-07-29 | 东北大学秦皇岛分校 | 一种基于cpls的红外光谱测量仪器标定迁移方法 |
WO2024023770A1 (en) * | 2022-07-28 | 2024-02-01 | «Uray Technologies» Llc | Systems, devices and methods for analyzing urine |
Family Cites Families (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
HU172951B (hu) | 1976-05-25 | 1979-01-28 | Gyogyszerkutato Intezet | Sposob poluchenija farmacevticheskojj kompozicii aktivnojj protiv mastitisa |
US4385590A (en) * | 1981-12-11 | 1983-05-31 | Bruce Mortensen | Apparatus for on-site detection of mastitis in milk animals |
CA2025330C (en) * | 1989-09-18 | 2002-01-22 | David W. Osten | Characterizing biological matter in a dynamic condition using near infrared spectroscopy |
JPH0447254A (ja) * | 1990-06-15 | 1992-02-17 | Snow Brand Milk Prod Co Ltd | 近赤外線を用いて脱脂乳、牛乳、クリーム及びチーズの成分含量を測定する方法及び装置 |
US5596992A (en) * | 1993-06-30 | 1997-01-28 | Sandia Corporation | Multivariate classification of infrared spectra of cell and tissue samples |
JPH07151677A (ja) * | 1993-12-01 | 1995-06-16 | Hitachi Ltd | 濃度計 |
WO1995022051A1 (en) * | 1994-02-09 | 1995-08-17 | Abbott Laboratories | Diagnostic flow cell device |
JP3462573B2 (ja) * | 1994-05-11 | 2003-11-05 | アークレイ株式会社 | 液体試料の成分濃度等を測定する方法及び装置 |
JPH11503236A (ja) | 1995-04-06 | 1999-03-23 | アルファ・ラヴァル・アグリ・アクチボラゲット | 流体中の粒子の定量決定方法及びその装置 |
JP3639721B2 (ja) * | 1998-06-26 | 2005-04-20 | 東亜ディーケーケー株式会社 | 測定システム |
NL1010540C2 (nl) * | 1998-11-12 | 2000-05-15 | Maasland Nv | Werkwijze voor het vaststellen van de aanwezigheid van bepaalde stoffen in melk en inrichting voor het toepassen van deze werkwijze. |
-
2001
- 2001-03-14 CN CN01800748A patent/CN1395683A/zh active Pending
- 2001-03-14 CZ CZ20014305A patent/CZ20014305A3/cs unknown
- 2001-03-14 RU RU2001135823A patent/RU2248554C2/ru not_active IP Right Cessation
- 2001-03-14 PL PL35103001A patent/PL351030A1/xx unknown
- 2001-03-14 EP EP01912358A patent/EP1191326A4/en not_active Ceased
- 2001-03-14 NZ NZ51572501A patent/NZ515725A/xx not_active IP Right Cessation
- 2001-03-14 WO PCT/JP2001/002018 patent/WO2001075421A1/ja not_active Application Discontinuation
- 2001-03-14 US US09/980,213 patent/US6748251B2/en not_active Expired - Lifetime
- 2001-03-14 HU HU0201283A patent/HUP0201283A2/hu unknown
- 2001-03-14 IL IL14678601A patent/IL146786A/xx not_active IP Right Cessation
- 2001-03-14 AU AU41138/01A patent/AU769362B2/en not_active Ceased
- 2001-03-14 JP JP2001572852A patent/JP3574851B2/ja not_active Expired - Fee Related
- 2001-03-14 CA CA 2376173 patent/CA2376173A1/en not_active Abandoned
- 2001-11-29 NO NO20015820A patent/NO20015820L/no not_active Application Discontinuation
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101918830A (zh) * | 2007-05-31 | 2010-12-15 | S.A.E.阿菲金公司 | 用于分析流体的系统和方法 |
CN101918830B (zh) * | 2007-05-31 | 2014-10-08 | S.A.E.阿菲金公司 | 用于分析流体的系统和方法 |
CN102384889A (zh) * | 2011-08-03 | 2012-03-21 | 长春迪瑞医疗科技股份有限公司 | 一种全自动便潜血分析仪的积分球式光学系统 |
CN102384889B (zh) * | 2011-08-03 | 2014-06-25 | 长春迪瑞医疗科技股份有限公司 | 一种全自动便潜血分析仪的积分球式光学系统 |
CN108027848A (zh) * | 2015-07-16 | 2018-05-11 | 皇家飞利浦有限公司 | 用于管理对人的炎症性自身免疫性疾病的处置的设备、系统和方法 |
Also Published As
Publication number | Publication date |
---|---|
EP1191326A4 (en) | 2005-04-06 |
JP3574851B2 (ja) | 2004-10-06 |
RU2248554C2 (ru) | 2005-03-20 |
WO2001075421A1 (fr) | 2001-10-11 |
US6748251B2 (en) | 2004-06-08 |
AU769362B2 (en) | 2004-01-22 |
PL351030A1 (en) | 2003-02-24 |
NO20015820L (no) | 2002-01-29 |
IL146786A (en) | 2005-03-20 |
NZ515725A (en) | 2003-01-31 |
IL146786A0 (en) | 2002-07-25 |
AU4113801A (en) | 2001-10-15 |
US20020183600A1 (en) | 2002-12-05 |
EP1191326A1 (en) | 2002-03-27 |
CZ20014305A3 (cs) | 2002-05-15 |
NO20015820D0 (no) | 2001-11-29 |
HUP0201283A2 (en) | 2002-08-28 |
CA2376173A1 (en) | 2001-10-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN1395683A (zh) | 用可见光线和/或近红外线诊断有无乳房炎的方法与装置 | |
CN1184467C (zh) | 谷物质量监测仪 | |
CN1366604A (zh) | 用可见光线和/或近红外线诊断有无乳房炎的方法与装置 | |
CN101339198B (zh) | 自动分析装置以及自动分析方法 | |
CN101482572A (zh) | 自动分析装置以及自动分析方法 | |
JP2013108855A (ja) | 判別フィルタ設計方法、判別方法、判別フィルタセット、判別装置、および、プログラム | |
WO2020000712A1 (zh) | 固定波长的拉曼散射快速采集和成像设备 | |
US20070030482A1 (en) | Spectrophotometer with adjustable light pathlength | |
CN1295494C (zh) | 集成化微型光学分析仪 | |
RU92191U1 (ru) | Спектрометрический экспресс-анализатор показателей качества молока и молочного напитка | |
JPH11173982A (ja) | 血清中蛋白質濃度の測定方法および装置 | |
JP3462573B2 (ja) | 液体試料の成分濃度等を測定する方法及び装置 | |
KR101893219B1 (ko) | 체액분석 방법 및 이를 이용한 체액분석 시스템 | |
CN2711726Y (zh) | 集成化微型光学分析仪 | |
CN218726656U (zh) | 能够提供全波长光束的pcr光路及检测装置 | |
RU2410671C1 (ru) | Спектрометрический экспресс-анализатор показателей качества молока и молочного напитка | |
RU103924U1 (ru) | Датчик со смещенными каналами | |
Misto et al. | Using multisample refractometer to determine the sugar content of sugarcane juice in sugar factory Besuki | |
CN117470801B (zh) | 基于多光谱同时检测的母乳成分快速分析装置 | |
KR100406838B1 (ko) | 다채널 분광분석을 위한 고속 주사형 겹살 분광광도계 | |
CN214408685U (zh) | 一种基于全反射折光法的食品分析仪 | |
CN211206249U (zh) | 一种便携式甲状腺功能快速检测仪 | |
RU2713022C1 (ru) | Способ контроля качества меда | |
Dortu et al. | A compact multichannel spectrometer for label-free monitoring of biochips for point-of-care testing | |
RU109564U1 (ru) | Спектрометрический экспресс-анализатор показателей качества молока и молочного напитка |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C02 | Deemed withdrawal of patent application after publication (patent law 2001) | ||
WD01 | Invention patent application deemed withdrawn after publication |