CN104541153A - 涉及利用光纤拉曼光谱的内窥镜进行实时癌症诊断的方法 - Google Patents
涉及利用光纤拉曼光谱的内窥镜进行实时癌症诊断的方法 Download PDFInfo
- Publication number
- CN104541153A CN104541153A CN201380035488.5A CN201380035488A CN104541153A CN 104541153 A CN104541153 A CN 104541153A CN 201380035488 A CN201380035488 A CN 201380035488A CN 104541153 A CN104541153 A CN 104541153A
- Authority
- CN
- China
- Prior art keywords
- spectrum
- raman
- probe
- spectrometer
- optical fiber
- 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
- 238000000034 method Methods 0.000 title claims abstract description 148
- 238000001069 Raman spectroscopy Methods 0.000 title claims abstract description 123
- 206010028980 Neoplasm Diseases 0.000 title claims description 39
- 201000011510 cancer Diseases 0.000 title claims description 29
- 238000001839 endoscopy Methods 0.000 title description 17
- 238000001228 spectrum Methods 0.000 claims abstract description 245
- 239000000523 sample Substances 0.000 claims abstract description 190
- 239000000835 fiber Substances 0.000 claims abstract description 46
- 238000005259 measurement Methods 0.000 claims abstract description 19
- 238000001237 Raman spectrum Methods 0.000 claims description 138
- 239000013307 optical fiber Substances 0.000 claims description 84
- 230000003595 spectral effect Effects 0.000 claims description 49
- 238000004458 analytical method Methods 0.000 claims description 45
- 238000000513 principal component analysis Methods 0.000 claims description 40
- VYPSYNLAJGMNEJ-UHFFFAOYSA-N Silicium dioxide Chemical compound O=[Si]=O VYPSYNLAJGMNEJ-UHFFFAOYSA-N 0.000 claims description 28
- 238000003860 storage Methods 0.000 claims description 26
- 238000012360 testing method Methods 0.000 claims description 26
- 239000010980 sapphire Substances 0.000 claims description 16
- 229910052594 sapphire Inorganic materials 0.000 claims description 16
- 230000004044 response Effects 0.000 claims description 14
- 239000000377 silicon dioxide Substances 0.000 claims description 11
- 235000012239 silicon dioxide Nutrition 0.000 claims description 11
- 230000005540 biological transmission Effects 0.000 claims description 10
- 230000005236 sound signal Effects 0.000 claims description 10
- 230000002159 abnormal effect Effects 0.000 claims description 9
- 238000012937 correction Methods 0.000 claims description 8
- 238000012706 support-vector machine Methods 0.000 claims description 8
- 238000005457 optimization Methods 0.000 claims description 7
- 230000001419 dependent effect Effects 0.000 claims description 5
- 238000010202 multivariate logistic regression analysis Methods 0.000 claims description 5
- 230000001613 neoplastic effect Effects 0.000 claims description 5
- 230000009826 neoplastic cell growth Effects 0.000 claims description 4
- 238000010239 partial least squares discriminant analysis Methods 0.000 claims description 2
- 238000012546 transfer Methods 0.000 abstract description 9
- 238000004445 quantitative analysis Methods 0.000 abstract description 2
- 210000001519 tissue Anatomy 0.000 description 117
- 230000006870 function Effects 0.000 description 83
- 238000003745 diagnosis Methods 0.000 description 39
- 238000001727 in vivo Methods 0.000 description 29
- 230000005284 excitation Effects 0.000 description 28
- 230000008569 process Effects 0.000 description 27
- 238000010586 diagram Methods 0.000 description 22
- 230000008859 change Effects 0.000 description 20
- 210000001989 nasopharynx Anatomy 0.000 description 17
- 210000002784 stomach Anatomy 0.000 description 17
- 108090000623 proteins and genes Proteins 0.000 description 16
- 102000004169 proteins and genes Human genes 0.000 description 16
- 230000000875 corresponding effect Effects 0.000 description 15
- 238000003384 imaging method Methods 0.000 description 15
- 238000004422 calculation algorithm Methods 0.000 description 12
- 210000000981 epithelium Anatomy 0.000 description 12
- 210000000867 larynx Anatomy 0.000 description 12
- 239000000203 mixture Substances 0.000 description 12
- 230000008520 organization Effects 0.000 description 12
- 230000007170 pathology Effects 0.000 description 12
- 238000010238 partial least squares regression Methods 0.000 description 11
- 108010010803 Gelatin Proteins 0.000 description 10
- 238000005516 engineering process Methods 0.000 description 10
- 229920000159 gelatin Polymers 0.000 description 10
- 239000008273 gelatin Substances 0.000 description 10
- 235000019322 gelatine Nutrition 0.000 description 10
- 235000011852 gelatine desserts Nutrition 0.000 description 10
- 150000002632 lipids Chemical class 0.000 description 10
- 230000003287 optical effect Effects 0.000 description 10
- 230000015654 memory Effects 0.000 description 9
- XKRFYHLGVUSROY-UHFFFAOYSA-N Argon Chemical compound [Ar] XKRFYHLGVUSROY-UHFFFAOYSA-N 0.000 description 8
- IJGRMHOSHXDMSA-UHFFFAOYSA-N Atomic nitrogen Chemical compound N#N IJGRMHOSHXDMSA-UHFFFAOYSA-N 0.000 description 8
- 208000005718 Stomach Neoplasms Diseases 0.000 description 8
- 230000014509 gene expression Effects 0.000 description 8
- 210000003128 head Anatomy 0.000 description 8
- 230000003902 lesion Effects 0.000 description 8
- 210000000056 organ Anatomy 0.000 description 8
- 238000012545 processing Methods 0.000 description 8
- RNAMYOYQYRYFQY-UHFFFAOYSA-N 2-(4,4-difluoropiperidin-1-yl)-6-methoxy-n-(1-propan-2-ylpiperidin-4-yl)-7-(3-pyrrolidin-1-ylpropoxy)quinazolin-4-amine Chemical compound N1=C(N2CCC(F)(F)CC2)N=C2C=C(OCCCN3CCCC3)C(OC)=CC2=C1NC1CCN(C(C)C)CC1 RNAMYOYQYRYFQY-UHFFFAOYSA-N 0.000 description 7
- 206010058314 Dysplasia Diseases 0.000 description 7
- 238000002790 cross-validation Methods 0.000 description 7
- 238000001514 detection method Methods 0.000 description 7
- 230000002496 gastric effect Effects 0.000 description 7
- 238000009499 grossing Methods 0.000 description 7
- 238000012544 monitoring process Methods 0.000 description 7
- 238000011160 research Methods 0.000 description 7
- 229920006395 saturated elastomer Polymers 0.000 description 7
- 150000001408 amides Chemical class 0.000 description 6
- JJWKPURADFRFRB-UHFFFAOYSA-N carbonyl sulfide Chemical compound O=C=S JJWKPURADFRFRB-UHFFFAOYSA-N 0.000 description 6
- 238000002189 fluorescence spectrum Methods 0.000 description 6
- 239000011159 matrix material Substances 0.000 description 6
- 210000000214 mouth Anatomy 0.000 description 6
- 238000012216 screening Methods 0.000 description 6
- 238000004611 spectroscopical analysis Methods 0.000 description 6
- 238000003841 Raman measurement Methods 0.000 description 5
- 238000009825 accumulation Methods 0.000 description 5
- 230000008878 coupling Effects 0.000 description 5
- 238000010168 coupling process Methods 0.000 description 5
- 238000005859 coupling reaction Methods 0.000 description 5
- 239000005350 fused silica glass Substances 0.000 description 5
- 206010017758 gastric cancer Diseases 0.000 description 5
- 238000000491 multivariate analysis Methods 0.000 description 5
- 230000035945 sensitivity Effects 0.000 description 5
- 201000011549 stomach cancer Diseases 0.000 description 5
- 239000000126 substance Substances 0.000 description 5
- 239000013598 vector Substances 0.000 description 5
- 210000001260 vocal cord Anatomy 0.000 description 5
- 229910052786 argon Inorganic materials 0.000 description 4
- 238000004590 computer program Methods 0.000 description 4
- 238000009826 distribution Methods 0.000 description 4
- 210000005081 epithelial layer Anatomy 0.000 description 4
- 238000001914 filtration Methods 0.000 description 4
- 239000007789 gas Substances 0.000 description 4
- 210000002175 goblet cell Anatomy 0.000 description 4
- 238000007689 inspection Methods 0.000 description 4
- 239000007788 liquid Substances 0.000 description 4
- QSHDDOUJBYECFT-UHFFFAOYSA-N mercury Chemical compound [Hg] QSHDDOUJBYECFT-UHFFFAOYSA-N 0.000 description 4
- 229910052753 mercury Inorganic materials 0.000 description 4
- 229910052757 nitrogen Inorganic materials 0.000 description 4
- 210000003800 pharynx Anatomy 0.000 description 4
- 238000007430 reference method Methods 0.000 description 4
- 208000023514 Barrett esophagus Diseases 0.000 description 3
- 208000023665 Barrett oesophagus Diseases 0.000 description 3
- 241000282414 Homo sapiens Species 0.000 description 3
- 230000005856 abnormality Effects 0.000 description 3
- 238000000540 analysis of variance Methods 0.000 description 3
- 238000001574 biopsy Methods 0.000 description 3
- 239000008280 blood Substances 0.000 description 3
- 210000004369 blood Anatomy 0.000 description 3
- 238000012512 characterization method Methods 0.000 description 3
- 238000007405 data analysis Methods 0.000 description 3
- 238000013480 data collection Methods 0.000 description 3
- 238000011156 evaluation Methods 0.000 description 3
- 238000002474 experimental method Methods 0.000 description 3
- 230000004907 flux Effects 0.000 description 3
- 230000003760 hair shine Effects 0.000 description 3
- 210000004072 lung Anatomy 0.000 description 3
- 239000000463 material Substances 0.000 description 3
- 230000005055 memory storage Effects 0.000 description 3
- 210000003097 mucus Anatomy 0.000 description 3
- 238000010606 normalization Methods 0.000 description 3
- 150000004032 porphyrins Chemical class 0.000 description 3
- 238000002310 reflectometry Methods 0.000 description 3
- 210000003296 saliva Anatomy 0.000 description 3
- 238000007619 statistical method Methods 0.000 description 3
- 201000000498 stomach carcinoma Diseases 0.000 description 3
- 230000001360 synchronised effect Effects 0.000 description 3
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Chemical compound O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 3
- COLNVLDHVKWLRT-QMMMGPOBSA-N L-phenylalanine Chemical compound OC(=O)[C@@H](N)CC1=CC=CC=C1 COLNVLDHVKWLRT-QMMMGPOBSA-N 0.000 description 2
- 238000000692 Student's t-test Methods 0.000 description 2
- 210000002659 acromion Anatomy 0.000 description 2
- 210000001124 body fluid Anatomy 0.000 description 2
- 239000010839 body fluid Substances 0.000 description 2
- 210000003679 cervix uteri Anatomy 0.000 description 2
- 210000001072 colon Anatomy 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 238000002405 diagnostic procedure Methods 0.000 description 2
- 230000004069 differentiation Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000002708 enhancing effect Effects 0.000 description 2
- 210000003238 esophagus Anatomy 0.000 description 2
- 239000000284 extract Substances 0.000 description 2
- 206010020718 hyperplasia Diseases 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 210000000936 intestine Anatomy 0.000 description 2
- 238000002095 near-infrared Raman spectroscopy Methods 0.000 description 2
- COLNVLDHVKWLRT-UHFFFAOYSA-N phenylalanine Natural products OC(=O)C(N)CC1=CC=CC=C1 COLNVLDHVKWLRT-UHFFFAOYSA-N 0.000 description 2
- 230000000750 progressive effect Effects 0.000 description 2
- 239000007787 solid Substances 0.000 description 2
- 238000012306 spectroscopic technique Methods 0.000 description 2
- 238000010183 spectrum analysis Methods 0.000 description 2
- 210000003932 urinary bladder Anatomy 0.000 description 2
- 230000000007 visual effect Effects 0.000 description 2
- RZVAJINKPMORJF-UHFFFAOYSA-N Acetaminophen Chemical compound CC(=O)NC1=CC=C(O)C=C1 RZVAJINKPMORJF-UHFFFAOYSA-N 0.000 description 1
- 241000894006 Bacteria Species 0.000 description 1
- VYZAMTAEIAYCRO-UHFFFAOYSA-N Chromium Chemical compound [Cr] VYZAMTAEIAYCRO-UHFFFAOYSA-N 0.000 description 1
- XDTMQSROBMDMFD-UHFFFAOYSA-N Cyclohexane Chemical compound C1CCCCC1 XDTMQSROBMDMFD-UHFFFAOYSA-N 0.000 description 1
- 108020004414 DNA Proteins 0.000 description 1
- 238000012327 Endoscopic diagnosis Methods 0.000 description 1
- 208000032843 Hemorrhage Diseases 0.000 description 1
- PMMYEEVYMWASQN-DMTCNVIQSA-N Hydroxyproline Chemical compound O[C@H]1CN[C@H](C(O)=O)C1 PMMYEEVYMWASQN-DMTCNVIQSA-N 0.000 description 1
- 206010054949 Metaplasia Diseases 0.000 description 1
- 108010063954 Mucins Proteins 0.000 description 1
- 102000015728 Mucins Human genes 0.000 description 1
- 241000158526 Nasalis Species 0.000 description 1
- 238000012952 Resampling Methods 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 230000004913 activation Effects 0.000 description 1
- 230000004075 alteration Effects 0.000 description 1
- 238000004164 analytical calibration Methods 0.000 description 1
- 230000000454 anti-cipatory effect Effects 0.000 description 1
- 238000005452 bending Methods 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 210000000621 bronchi Anatomy 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 210000004027 cell Anatomy 0.000 description 1
- 210000003467 cheek Anatomy 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 229910052804 chromium Inorganic materials 0.000 description 1
- 239000011651 chromium Substances 0.000 description 1
- 238000003759 clinical diagnosis Methods 0.000 description 1
- 238000000576 coating method Methods 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 230000001276 controlling effect Effects 0.000 description 1
- 230000002596 correlated effect Effects 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 229910003460 diamond Inorganic materials 0.000 description 1
- 239000010432 diamond Substances 0.000 description 1
- 201000010099 disease Diseases 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 239000012153 distilled water Substances 0.000 description 1
- PMMYEEVYMWASQN-UHFFFAOYSA-N dl-hydroxyproline Natural products OC1C[NH2+]C(C([O-])=O)C1 PMMYEEVYMWASQN-UHFFFAOYSA-N 0.000 description 1
- 229940079593 drug Drugs 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 201000006549 dyspepsia Diseases 0.000 description 1
- 238000013399 early diagnosis Methods 0.000 description 1
- 238000011846 endoscopic investigation Methods 0.000 description 1
- YQGOJNYOYNNSMM-UHFFFAOYSA-N eosin Chemical compound [Na+].OC(=O)C1=CC=CC=C1C1=C2C=C(Br)C(=O)C(Br)=C2OC2=C(Br)C(O)=C(Br)C=C21 YQGOJNYOYNNSMM-UHFFFAOYSA-N 0.000 description 1
- 210000002919 epithelial cell Anatomy 0.000 description 1
- 210000000416 exudates and transudate Anatomy 0.000 description 1
- 238000000799 fluorescence microscopy Methods 0.000 description 1
- 235000013305 food Nutrition 0.000 description 1
- 210000004051 gastric juice Anatomy 0.000 description 1
- 210000001156 gastric mucosa Anatomy 0.000 description 1
- 239000011521 glass Substances 0.000 description 1
- 229910052736 halogen Inorganic materials 0.000 description 1
- 210000001983 hard palate Anatomy 0.000 description 1
- 201000000615 hard palate cancer Diseases 0.000 description 1
- 230000036541 health Effects 0.000 description 1
- HLUCICHZHWJHLL-UHFFFAOYSA-N hematein Chemical compound C12=CC=C(O)C(O)=C2OCC2(O)C1=C1C=C(O)C(=O)C=C1C2 HLUCICHZHWJHLL-UHFFFAOYSA-N 0.000 description 1
- 230000001744 histochemical effect Effects 0.000 description 1
- 238000010562 histological examination Methods 0.000 description 1
- 230000003118 histopathologic effect Effects 0.000 description 1
- 229960002591 hydroxyproline Drugs 0.000 description 1
- 238000012623 in vivo measurement Methods 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 230000001678 irradiating effect Effects 0.000 description 1
- 238000012804 iterative process Methods 0.000 description 1
- 210000004185 liver Anatomy 0.000 description 1
- 230000004807 localization Effects 0.000 description 1
- 238000007726 management method Methods 0.000 description 1
- 239000003550 marker Substances 0.000 description 1
- 230000011278 mitosis Effects 0.000 description 1
- 230000000877 morphologic effect Effects 0.000 description 1
- 238000012627 multivariate algorithm Methods 0.000 description 1
- VFBILHPIHUPBPZ-UHFFFAOYSA-N n-[[2-[4-(difluoromethoxy)-3-propan-2-yloxyphenyl]-1,3-oxazol-4-yl]methyl]-2-ethoxybenzamide Chemical class CCOC1=CC=CC=C1C(=O)NCC1=COC(C=2C=C(OC(C)C)C(OC(F)F)=CC=2)=N1 VFBILHPIHUPBPZ-UHFFFAOYSA-N 0.000 description 1
- 238000013450 outlier detection Methods 0.000 description 1
- 125000004430 oxygen atom Chemical group O* 0.000 description 1
- 229960005489 paracetamol Drugs 0.000 description 1
- 230000035515 penetration Effects 0.000 description 1
- 229910052698 phosphorus Inorganic materials 0.000 description 1
- 230000000258 photobiological effect Effects 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 238000004321 preservation Methods 0.000 description 1
- 230000001737 promoting effect Effects 0.000 description 1
- 238000012113 quantitative test Methods 0.000 description 1
- 238000011897 real-time detection Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000000985 reflectance spectrum Methods 0.000 description 1
- 238000002271 resection Methods 0.000 description 1
- 230000029058 respiratory gaseous exchange Effects 0.000 description 1
- 238000005316 response function Methods 0.000 description 1
- 230000000717 retained effect Effects 0.000 description 1
- 210000003491 skin Anatomy 0.000 description 1
- 210000001584 soft palate Anatomy 0.000 description 1
- 230000000638 stimulation Effects 0.000 description 1
- 238000003756 stirring Methods 0.000 description 1
- 238000005728 strengthening Methods 0.000 description 1
- 210000004003 subcutaneous fat Anatomy 0.000 description 1
- 238000001356 surgical procedure Methods 0.000 description 1
- 238000012549 training Methods 0.000 description 1
- FGMPLJWBKKVCDB-UHFFFAOYSA-N trans-L-hydroxy-proline Natural products ON1CCCC1C(O)=O FGMPLJWBKKVCDB-UHFFFAOYSA-N 0.000 description 1
- 230000001052 transient effect Effects 0.000 description 1
- UFTFJSFQGQCHQW-UHFFFAOYSA-N triformin Chemical compound O=COCC(OC=O)COC=O UFTFJSFQGQCHQW-UHFFFAOYSA-N 0.000 description 1
- 230000002792 vascular Effects 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
- 210000001835 viscera Anatomy 0.000 description 1
- 238000011179 visual inspection Methods 0.000 description 1
- 230000001755 vocal effect Effects 0.000 description 1
- 229910052724 xenon Inorganic materials 0.000 description 1
- FHNFHKCVQCLJFQ-UHFFFAOYSA-N xenon atom Chemical compound [Xe] FHNFHKCVQCLJFQ-UHFFFAOYSA-N 0.000 description 1
Classifications
-
- 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
-
- 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/62—Systems in which the material investigated is excited whereby it emits light or causes a change in wavelength of the incident light
- G01N21/63—Systems in which the material investigated is excited whereby it emits light or causes a change in wavelength of the incident light optically excited
- G01N21/65—Raman scattering
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B1/00—Instruments for performing medical examinations of the interior of cavities or tubes of the body by visual or photographical inspection, e.g. endoscopes; Illuminating arrangements therefor
- A61B1/00163—Optical arrangements
- A61B1/00165—Optical arrangements with light-conductive means, e.g. fibre optics
-
- 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
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
- G01J3/00—Spectrometry; Spectrophotometry; Monochromators; Measuring colours
- G01J3/28—Investigating the spectrum
- G01J3/44—Raman spectrometry; Scattering spectrometry ; Fluorescence 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/27—Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands using photo-electric detection ; circuits for computing concentration
- G01N21/274—Calibration, base line adjustment, drift correction
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B2560/00—Constructional details of operational features of apparatus; Accessories for medical measuring apparatus
- A61B2560/02—Operational features
- A61B2560/0223—Operational features of calibration, e.g. protocols for calibrating sensors
- A61B2560/0228—Operational features of calibration, e.g. protocols for calibrating sensors using calibration standards
- A61B2560/0233—Optical standards
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
- G01J3/00—Spectrometry; Spectrophotometry; Monochromators; Measuring colours
- G01J3/28—Investigating the spectrum
- G01J3/44—Raman spectrometry; Scattering spectrometry ; Fluorescence spectrometry
- G01J2003/4424—Fluorescence correction for Raman 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/62—Systems in which the material investigated is excited whereby it emits light or causes a change in wavelength of the incident light
- G01N21/63—Systems in which the material investigated is excited whereby it emits light or causes a change in wavelength of the incident light optically excited
- G01N21/65—Raman scattering
- G01N2021/653—Coherent methods [CARS]
- G01N2021/656—Raman microprobe
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2201/00—Features of devices classified in G01N21/00
- G01N2201/08—Optical fibres; light guides
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2201/00—Features of devices classified in G01N21/00
- G01N2201/12—Circuits of general importance; Signal processing
- G01N2201/127—Calibration; base line adjustment; drift compensation
- G01N2201/12746—Calibration values determination
- G01N2201/12753—Calibration values determination and storage
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2201/00—Features of devices classified in G01N21/00
- G01N2201/12—Circuits of general importance; Signal processing
- G01N2201/129—Using chemometrical methods
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- General Health & Medical Sciences (AREA)
- Pathology (AREA)
- Surgery (AREA)
- Engineering & Computer Science (AREA)
- Medical Informatics (AREA)
- Heart & Thoracic Surgery (AREA)
- Biomedical Technology (AREA)
- Molecular Biology (AREA)
- Biophysics (AREA)
- Animal Behavior & Ethology (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- General Physics & Mathematics (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- Immunology (AREA)
- Optics & Photonics (AREA)
- Radiology & Medical Imaging (AREA)
- Mathematical Physics (AREA)
- Theoretical Computer Science (AREA)
- Investigating, Analyzing Materials By Fluorescence Or Luminescence (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
- Endoscopes (AREA)
- High Energy & Nuclear Physics (AREA)
- Spectrometry And Color Measurement (AREA)
- Measuring And Recording Apparatus For Diagnosis (AREA)
Abstract
本发明公开一种实现用于光纤拉曼光谱系统的定量分析的仪器独立测量的方法,该系统包括激光源、光谱仪和光纤探针,该光纤探针将光从激光源发送光到目标,并且将散射光返回到光谱仪,该方法包括从激光源发送光到具有已知光谱的标准目标,记录来自标准目标的散射光的校准光谱,比较已知的光谱和校准系统,并且生成探针和/或探针系统传递函数,以及存储传递函数。本发明还提供了结合其它公开方法可选地执行实时诊断拉曼光谱的方法。
Description
本申请要求2012年7月2日提交的美国临时申请61/667,384和在2013年4月23日提交的GB申请No.1307338.2的权益,其两者通过引用整体在此并入此案以供用于所有目的之参考。
技术领域
本公开涉及用于在内窥镜实时癌症诊断的在线生物光谱软件平台,和用于在光纤拉曼光谱中用于定量分析的与仪器无关的测量方法。
背景技术
拉曼光谱术是使用单色光的非弹性或拉曼散射的技术。传统地,单色光源是在可见或近红外(“NIR”)范围中的激光。散射光子的能量响应于与在照射材料中的振动模式或激发的相互作用,向上或向下移动,改变散射光子的波长。因此,来自散射光的光谱可提供关于散射材料的信息。
NIR拉曼光谱术已知为用于在若干器官中体内癌前病变和癌性细胞和组织的表征和诊断。该技术是理想的,因为它可以无创或微创,不需要活组织检查或组织的其它去除。已知的是使用在两个波长范围中的NIR拉曼光谱。第一个是所谓的指纹(“FP”)范围,其中由于例如来自在用于组织表征和诊断的该光谱区域中包含的蛋白质、DNA和脂质内容的高度特异性生物分子信息的丰富度,波长数从800到1800cm-1。该波长范围的缺点在于当与通常使用的785nm激光源一起使用时,可生成强烈的组织自发荧光背景信号。此外,在探针使用光纤时,拉曼信号从在光纤中的熔融二氧化硅散射。特别地,在电荷耦合器件(“CCD”)用于测量散射光谱时,自发荧光信号可使CCD饱和并且与在该波长范围中固有的非常弱的拉曼信号的检测干扰。
采用光纤拉曼光谱术作为一种技术的另一个问题在于仪器的标准化。光纤拉曼光谱技术已经主要限制为单个系统,并且尚未有转用于多中心临床试验或常规医疗诊断中的相关尝试。这主要是因为拉曼光谱仪仪器通常存在差异(即光学、响应函数、对准、通量等)并且一般产生非常不同的拉曼光谱。此外,光纤拉曼探针的使用寿命有限并且必须定期更换或互换。不幸地,使用不同的光纤探针采集的拉曼数据无法进行比较,因为每个光纤探针具有其本身的独特背景,以及与不同的透射光谱特性相关联。不同的透射特性显著地使光谱强度失真,使采用不同光纤探针获得的组织拉曼光谱不可比较。因此,在主临床平台上开发的多变量诊断算法不能施加到次级临床平台。特别地,组织拉曼强度的定量测量是在光纤生物医学拉曼应用中最具挑战性的问题之一。仪器/光纤探针独立强度校准和标准化对于在生物医学中光纤拉曼光谱的全球使用的实现至关重要。出于这个原因,使用“总”探针(‘master’probe)构建的多变量统计诊断模型不能施加到采用“从属”探针(‘slave’probe)测量的光谱。为了使得拉曼技术成为用于在全球范围的癌症筛查的广泛工具,有必要将拉曼光谱仪和光纤探针这两者进行标准化,特别是在生物医学应用中。大多数的报道研究已经聚焦在用于不具有光纤探针的简单化学混合物的测量的跨拉曼光谱仪标准化。一般简单化学混合物的拉曼光谱不能与异种生物组织样本的光纤拉曼光谱比较。
采用仪器间的标准化结果的进一步问题在于与激光激发功率相关联的光谱变化。传统地,拉曼光谱的标准化,其保留了通常光谱的形状,但这消除了绝对定量光谱特性。已知的是尝试监测在光纤拉曼探针中传送的激光功率,例如通过在光纤尖端中嵌入钻石或在激光路径中定位聚合物盖作为参考来进行监测。然而,这些解决方案不令人满意,并且可能导致在所需光谱区域中的干涉。
在使用用于体内癌症和癌前病变条件诊断的光学光谱技术(包括反射荧光和拉曼)中的其它问题在于数据分析大多限于后处理和离线算法开发。对于内窥镜分析正是如此,因为在内窥镜期间收集的大量光谱是异常值。具有允许用于内窥镜的实时诊断的系统将是有用的。
发明内容
根据第一方面,本发明提供了校准光纤拉曼光谱系统的方法,该系统包括激光源、光谱仪和光纤探针,该光纤探针将光从激光源发送到目标,并将散射光返回到光谱仪;该方法包括从激光源发送光到具有已知光谱的标准目标,记录来自标准目标的散射光的校准光谱,比较已知的光谱和校准系统,并且生成传递函数,以及存储传递函数。
本方法还可包括随后照射测试对象、记录光谱并且根据存储的传递函数校准光谱的步骤。
该方法可包括记录用于多个光纤探针的每个的记录校准光谱,计算用于包括每个所述探针的系统的传递函数,并且将传递函数与对应探针相关联。
所述光谱仪具有相关联的光谱仪传递函数,所述探针可具有相关联的探针传递函数,传递函数可以是光谱仪传递函数和探针传递函数中的一种函数。
该方法可包括在主光谱仪系统上计算主光纤探针的第一传递函数,和次级光纤探针的第二传递函数,并且基于第一传递函数和第二传递函数计算(探针间)校准函数。
该方法可包括将校准函数与次级光纤探针相关联。
该方法可包括在次级光谱仪系统上使用主光纤探针,生成次级系统传递函数,以及存储所述次级系统传递函数。
该方法可包括使用具有次级光谱仪系统的次级光纤探针,以及根据校准函数修正存储的次级系统传递函数。
该方法可包括根据主光谱仪系统执行次级光谱仪系统的波长轴校准的初始步骤。
根据第二方面,本发明提供了操作拉曼光谱仪系统的方法,该系统包括激光源、光谱仪和光纤探针,所述光纤探针将光从激光源发送到目标,并且将散射光返回到光谱仪;该方法包括从激光源发送光到具有已知光谱的目标,记录来自目标的散射光的光谱,以及根据存储的传递函数修正记录光谱。
所述存储的传递函数可将光谱仪和光纤探针相关联。
所述存储的传递函数可将光谱仪和主光纤探针相关联,并且该方法还可包括根据与光纤探针相关联的存储校准函数来修正存储的传递函数。
根据第三方面,本发明提供了一种拉曼光谱系统,其包括激光源、光谱仪和光纤探针,以及存储的传递函数;所述光纤探针将光从激光源发送到目标并且将散射光返回到光谱仪;,该系统可操作地用于从激光源发送光到具有已知光谱的目标,记录来自目标的散射光的光谱,以及根据存储的传递函数修正记录光谱。
所述存储的传递函数可将光谱仪和光纤探针相关联。
所述存储的传递函数可将光谱仪和主光纤探针相关联,并且该方法还可包括根据与光纤探针相关联的存储校准函数来修正存储的传递函数。
根据第四方面,本发明提供了估计拉曼光谱仪系统中发送的激光功率的方法,该系统包括激光源、光谱仪和光纤探针,该光纤探针将光从激光源发送到目标,并且将散射光返回到光谱仪;该方法包括从激光源发送光到多个目标,对于每个测试目标,测量来自激光源的光的发送功率和在光谱仪处散射光的光谱,执行捕捉光谱的多变量分析,其中测量的发送功率用作为因变量,,以及存储所得的模型。
该方法可包括如下步骤:发送激光到测试目标,向模型供应捕捉的光谱,以及计算发送功率的估值。
根据第五方面,本发明提供了从具有光纤探针的光纤拉曼光谱系统减去背景信号的方法,该方法包括如下步骤:
a)存储背景光谱,
b)接收测试光谱,
c)使用一个或多个参考峰来估计背景影响,
d)基于估计的背景影响将背景光谱与校准因子相乘,并将其从测试光谱减去,
e)检查测试光谱余留的背景影响,以及
f)如果背景影响可忽略不计,则输出测试光谱,否则重复步骤(c)到(e)。
所述一个或多个参考峰可包括与光纤探针中的二氧化硅或蓝宝石对应的一个或多个峰。
根据第六方面,本发明提供了在内窥镜期间使用拉曼光谱用于实时诊断的计算机实施方法。该方法包括接收与组织相关联的至少一个光谱;在使用光谱的模型中分析至少一个光谱以确定得分,其中所述得分指示组织是癌性的可能性;以及输出所述得分。
在一些实施例中,模型使用解读函数来生成,所述解读函数选自由偏最小二乘判别分析、主成分分析线性判别分析、蚁群优化线性判别分析、分类和回归树、支持向量机和自适应增强组成的群组。
在一些实施例中,所述至少一个光谱由拉曼光谱生成。分析模型中的至少一个光谱可包括分析在第一模型和第二模型中的至少一个光谱。在一些实施例中所述模型基于分析的组织来选择。在一些实施例中,得分指示组织是否为正常、异常增生或瘤形成。
在某些方面中,分析至少一个光谱包括:执行异常分析;响应于确定至少一个光谱是异常的异常分析,抑制光谱。执行异常分析可包括主成分分析。
在某些方面中,音频发射装置响应于确定至少一个光谱是异常的异常分析,来发射音频信号。响应于确定至少一个光谱是异常的异常分析来获得至少一个附加光谱,该获得的至少一个附加光谱由系统接收以用于分析。
在一些实施例中,音频发射装置发射识别组织为正常的、异常增生或瘤形成的音频信号。在一些实施例中,与正常的、异常增生或瘤形成的每个相关联的音频信号是不同的。
在一些实施例中,诊断在内窥镜程序期间发生。
本发明还提供了用于执行计算机实施方法的系统以及具有在其上用于执行计算机实施方法的指令的非暂时性计算机可读介质。
附图说明
本发明所公开的系统和方法的实施例通过仅参考附图的示例来描述。
图1是根据一个实施例的拉曼光谱系统的示意图。
图1a是在更大比例上图1的内窥镜端部视图。
图1b是更详细的图1a的内窥镜拉曼探针的视图。
图2是示出测量的荧光光谱与参考标准的比较的曲线图。该校准函数也被示出。
图3是第一校准方法的示意图。
图4a是示出作为第一校准方法使用的第一过程的流程图。
图4b是示出作为第一校准方法使用的第二过程的第一部分的流程图。
图4c是示出作为第一校准方法使用的第二过程的第二部分的流程图。
图5是示出在主光谱仪和次级光谱仪之间氩气/汞灯的波长对准的曲线图。
图6是示出使用第二校准方法的主光谱仪和次级光谱仪的光谱校准的图。
图7a是示出作为第一校准方法使用的第一过程的流程图。
图7b是示出作为第二校准方法使用的第二过程的第一部分的流程图。
图7c是示出作为第二校准方法使用的第二过程的第二部分的流程图。
图8是示出采用主探针和从属探针和探针校准传递函数测量的荧光标准的曲线图。
图9a是分别比较使用主探针和从属探针的未校准的主光谱仪和次级光谱仪的组织拉曼光谱的曲线图。
图9b是在使用第一校准方法的重新校准之后,分别来自使用主探针和从属探针的主光谱仪和次级光谱仪的组织拉曼光谱的曲线图。
图9c是在使用第二校准方法的重新校准之后,分别来自使用主探针和从属探针的主光谱仪和次级光谱仪的光谱的曲线图。
图10是在校准之前和之后在来自胃部的体内组织拉曼光谱上的主成分分析计分散射点图。
图11是示出由于在拉曼光谱中的光纤探针导致的背景光谱峰的曲线图。
图12是示出采用激发激光功率的拉曼光谱变化的曲线图。
图13a是示出生成用于预测激光功率的模型的方法的流程图。
图13b是示出预测激光功率的方法的流程图,
图14a是示出用于任何数量的包括潜在变量的均方根误差的曲线图。
图14b示出用于图10的方法的潜在变量的负载因子和回归因子,
图15是示出在体内测试对象中测量激光功率对预测激光功率的的曲线图。
图16是示出减去探针背景信号的方法的流程图。
图17是示出从手掌和光纤二氧化硅和蓝宝石背景接收的光谱的曲线图。
图18是比较图16的拉曼光谱和在背景消除之后的光谱的曲线图。
图19是示出方法组合的流程图。
图20是用于根据一个实施例的实时癌症诊断的光谱采集和处理流程的系统的架构图。
图21是示出用于根据一个实施例的光谱采集和处理流程的示意图的流程图。
图22A和22B是用于使用根据两个实施例的实时癌症诊断的图形用户界面(GUI)。
图23是从305个胃部患者采集的正常(n=2465)和癌症(n=283)胃部组织的体内平均拉曼光谱。
图24示出从光谱训练数据库计算的主成分(PC)负荷。
图25是两个诊断显著的PC计分(PC1VS PC2)的散射点图。
图26是展示用于从10个预期癌症样本采集的105个拉曼光谱(45个正常,30个癌症,30个异常)的Hotelling’s T2与Q-残差。
图27是基于与留一交叉验证一起的PLS-DA建模,属于预期正常(n=45)和癌症(n=30)胃部组织的后验概率值的散射点图。
图28示出从用于预期预测的光谱数据库计算的受试者工作特征曲线(ROC)以及用于正常和癌症胃部组织预期预测的曲线。
图29示出通过使用不同的785-nm激光激发功率(即10mW、30mW和60mW)在内唇的体内组织拉曼光谱±1SD的减去自体荧光和强度校准平均值。
图30a示出基于留一交叉验证以及对数据的线性拟合,在实际的和使用PLS回归模型的预测激光激发功率之间的关系。
图30b示出基于独立验证,在实际的和使用PLS回归的预测激光激发功率之间的关系。
图31示出用在60mW激光激发功率处测量的制备的不同浓度的(即以重量计20%、25%、30%、35%、40%、45%和50%)明胶组织模型(phantom)的拉曼光谱。
图32示出在采用预测激光功率的校准之后在实际的和在组织模型中预测的明胶浓度之间的相关关系。
图33示出在临床内窥镜检查期间从咽隐窝采用0.1s采集的代表性体内原始拉曼光谱。图33的插图是在去除强度自发荧光背景之后的处理组织拉曼光谱。
图34示出后鼻咽部(PN)(n=521)、咽隐窝(FOR)(n=157)和喉头声带(LVC)(n=196)的体内(主体间)平均拉曼光谱±1标准偏差(SD)。应注意的是,平均拉曼光谱垂直移动以便更好地观察。同样示出了在白光反射率(WLR)和窄带(NB)成像引导下来自后鼻咽部(上)、咽隐窝(中)和喉头声带(下)的体内光纤拉曼内窥镜采集。
图35示出PN(n=18)、FOR(n=18)和LVC(n=17)的体内(主体内)平均拉曼光谱±1SD。应注意的是,平均拉曼光谱垂直移动以便更好地观察。
图36示出不同解剖组织类型(主体间)的差异光谱±1SD的比较:[后鼻咽部(PN)-喉声带(LVC)];[后鼻咽部(PN)-咽隐窝(FOR)]和[喉声带(LVC)-咽隐窝(FOR)]。
图37示出来自人体体液的(鼻涕、唾液和血液)的可能混杂因素的体内拉曼光谱。
图38示出分辨在头部和颈部中不同组织之间生物分子变异的PC负荷,表示光谱方差57.41%(PC1:22.86%;PC2:16.16%;PC3:8.13%;PC4:6.22%;PC5:4.05%)的总计。
图39提供用于不同组织类型(即PN、FOR和LVC)的5个PCA计分的方框图。在每个槽口箱内的线表示中位数,但箱子的下边界和上边界分别指示第一(25.0%百分位数)和第三(75.0%百分位数)四分位数。误差条(须状物)表示1.5倍四分范围。p值同样在不同的组织类型之间给出。
图40A示出在临床内窥镜检查期间采集的鳞状上皮层(n=165)、柱状上皮层(n=907)、巴雷特食管(n=318)、高度异型增生(n=77)的平均体内共焦拉曼光谱。
图40B-E示出B)与测量组织部位对应的代表性组织学切片(苏木精和曙红(H&E)染色)。鳞状上皮层;(C)在缺乏杯状细胞情况下的柱状层食管(n=907),×200;(D)Barrett’s食管,其中正常的复层鳞状上皮由包含杯状细胞的肠化生上皮取代,×200;(E)示出架构和细胞学异型两者以及具有分支和乳突状形成、细胞学pleomophism和极性损失的高度异型增生,×100。
图41A示出使用共焦拉曼内窥镜技术、属于“正常”柱状上皮层(CLE),(ii)“低风险”肠化生(IM)(iii)‘高风险’高度异型增生(HGD)的预测后验概率的二维三角图。
图41B示出‘正常’CLE,(ii)‘低风险’IM(iii)‘高风险’HGD的二分法判别(dichotomous discriminations)的受试者工作特征曲线(ROC)曲线。在ROC曲线下方的面积(AUC)分别是0.88、0.84和0.90。
具体实施方式
本文所提供的是用于实现在内窥镜处在不同的器官(例如,胃肠道(胃、食道、结肠)、膀胱、肺、口腔、鼻咽、喉、子宫颈、肝、皮肤等)中肿瘤性病变的实时检测的生物医学光谱(即,反射率、荧光和拉曼光谱)的在线系统和方法。诊断方法整合了激发源的同步、积分时间调整、数据采集、预处理、异常分析和概率多变量诊断(即,偏最小二乘-判别分析(PLS-DA)、主成分分析(PCA)-线性判别分析(LDA)、蚁群优化(ACO)-LDA、分类和回归树(CART)、支持向量机(SVM)、自适应增强(AdaBoost)等),包括基于不同器官的全面光谱数据库(即,拉曼、荧光、反射率)的多级诊断。
在一个实施例中,所公开的系统和方法结合了在线诊断框架与最近开发的用于在内窥镜处在上消化道中癌前病变和癌症的早期诊断和检测的多模态图像引导(WLR/NBI/AFI)的拉曼光谱平台。组织拉曼光谱的积累和具有0.5秒预定义上限的积分时间的自动比例缩放允许具有改进的SNR的体内组织光谱的即时采集,同时防止CCD信号饱和。这是对于内窥镜诊断特别重要,其中自体荧光强度在由组织中不同内生荧光可能导致的不同解剖区域(例如,在胃部中的胃窦和胃体,肺中的支气管)之间显著改变。
现在详细参考具体的附图,需要强调的是所示的细节是通过举例的方式,并且为了对优选实施例进行说明性讨论的目的,并且呈现出为了提供被认为是所公开的系统和方法的最有用和最容易理解的描述。在这点上,没有尝试来比对于所公开系统和方法的基本理解所必要的公开所更为详细的结构细节,采用附图的描述使得对于本领域技术人员而言,所公开的系统和方法的几种形式如何在时间中进行实施,是显而易见的。
在详细解释所公开的系统和方法的至少一个实施例之前,应当理解的是,本公开并不将其应用限于在下面的描述中或在附图中示出的组件构造和布置的细节。所公开的系统和方法适用于其它实施例或以各种方式实践或执行。此外,可以理解的是本文使用的措辞和术语是为了描述的目的,而不应被视为限制。
现在参考图1,诊断仪器包括根据一个实施例的内窥镜系统,在10处示出。内窥镜本身在11处示出,并且内窥镜11的测量头通常在图1a中示出。为提供引导和被测试区域的可见查看,内窥镜11装有合适的视频系统。来自氙灯光源的光被发送到在内窥镜12末端中的照明窗口15。响应于白光反射成像、窄带成像或自发荧光成像的CCD 16和17接收反射光并发送视频数据,以允许所测试组织的可视检查和引导内窥镜至期望的位置。共焦拉曼探头在18处显示,并且在图1b中更详细。
拉曼光谱系统通常在20处示出。在具有大约785nm的输出波长的二极管激光器的本示例中,单色激光源在21处示出。来自激光二极管21的光通过近端带通滤波器22,其包括具有±2.5nm的半宽度而以在785nm处为中心的窄带通滤波器。该光通过穿过耦合件23而进入作为光纤束一部分的激发光纤25中。激发光纤25具有200μm的直径和0.22的数值孔径(numerical aperture,‘NA’)。由激发光纤25传输的光进入在内窥镜11端部的球透镜26,在本实施例中,包括具有约1.0mm直径和折射率n=1.77的蓝宝石球透镜。如在图1b中所示,从激发光纤25传输的光在球透镜26内被内部反射。在球透镜与被测试的组织接触时,如这里在27处示出,来自激发光纤25的透射光的至少一部分经过在组织27内的拉曼散射到约140μm的深度。散射光再次在球透镜26中被内部反射,并在也提供为光纤束一部分的多个收集光纤28中接收。在本示例中26个100μm的收集光纤被使用,其NA值为0.22。收集光纤28可以任何适当的结构进行安置,例如在包围激发光纤25的圆形布置中。
由收集光纤28返回的收集散射光穿过长波通内联收集滤波器29,其类似地具有在约800nm处的截断值(cutoff)。所述蓝宝石球透镜26,激发和收集光纤25、28,带通滤波器22、长波通滤波器29的配置提供了用于从组织27选择性地收集反向散射的拉曼光子的良好系统。
散射返回光然后在光谱仪30处分开,并且所得光谱在感光阵列34处成像,在本示例中是电荷耦合器件(“CCD”)。在35处示出的计算机控制系统的操作、处理并存储光谱和控制数据,并且提供结果和数据给用户。
在一个实施例中,计算机35包括至少一个耦合到芯片组的处理器。也耦合到芯片组的是内存、存储装置、键盘、图形适配器、指向装置、声音发射装置和网络适配器。显示器耦合到图形适配器。在一个实施例中,芯片组的功能由内存控制器集线器和I/O控制器集线器提供。在另一个实施例中,内存直接耦合到处理器,而不是芯片组。
存储装置是能够保持数据的任何装置,例如硬盘驱动器、压缩光盘只读内存(CD-ROM)、DVD或固态存储器装置。所述存储器保持由处理器使用的指令和数据。指向装置可以是鼠标、轨迹球,或其它类型的指示装置,并且与键盘组合使用以输入数据到计算机系统中。图形适配器显示图像和在显示器上的其它信息。网络适配器耦合计算机系统到广域网或局域网。
如本领域已知的,计算机35可具有除前面描述的不同和/或其它组件。此外,计算机可以缺少某些组件。另外,该存储装置可以是本地和/或远程计算机(诸如在存储区域网(SAN)内实施)。
如本领域已知的,该计算机适于执行用于提供本文所描述的功能性的计算机程序模块。如本文所使用的,术语“模块”是指用于提供特定功能性的计算机程序逻辑。因此,模块可以在硬件、固件和/或软件中来实现。在一个实施例中,程序模块被存储在存储装置上,加载到内存中,并且由处理器执行。
本文这里所描述的实体的实施例可以包括除了在此描述的模块的其它和/或不同的模块。此外,归属于模块的功能性可以通过在其它实施例中其它或不同的模块来执行。此外,为了清楚和方便性,本说明书中偶尔省略了术语“模块”。
计算机35还对光谱数据执行预处理。因为测量组织的拉曼光谱基本上被组织的自体荧光背景遮蔽,所以体内组织的拉曼光谱的预处理有必要提取弱拉曼信号。从体内组织测量的原始拉曼光谱代表了弱拉曼信号、强烈的自体荧光背景和噪声的组合。光谱被首先归一化成积分时间和激光功率。光谱然后使用一阶的Savitzky-Golay平滑滤波器(3个像素的窗口宽度)平滑化以减少噪声。五阶多项式被发现对于拟合在噪声平滑光谱中的自体荧光背景是最优的,并且随后这个多项式从原始光谱中减去,以单独产生组织的拉曼光谱。计算机35还可以包括用于癌前病变和癌症检测的诊断算法。
光谱仪和光纤探针校准
已知的是,不同的光谱仪将具有不同的传递函数,即:即使在使用相同的源进行照射时,也会在光谱内呈现出不同的强度变化。如在图2中所示,其示出的是来自标准源的光谱。在本示例中的标准源是指:当被激光,如激光源2,激发时,能够发射出已知荧光光谱的荧光标准目标。所述荧光标准目标必须是一致和稳定的,并在激光激发下(例如785nm)发射宽荧光光谱。随着时间的推移,荧光光谱必须是稳定的,并且有效地表征在感兴趣的整个光谱区域上(例如,400cm-1-1800cm-1,2000cm-1-3800cm-1)的光谱透射特性。一个示例是铬掺杂的玻璃。从两个光谱仪所得的光谱可以看出,这是明显不同的。为了补偿光谱仪响应或传递函数,已知的是施加校准函数,所述校准函数能够校准从光谱仪接收的光谱。示例和校准函数在图2中示出,当施加到光谱仪的对应光谱时,这将使光谱带到具有已知标准光谱的线中。
使用荧光标准源,传递函数,即:光谱仪的波长相关响应可由如下给出:
(方程式1)其中F(λ)是校正的荧光标准光谱,S(λ)是荧光标准源的测量光谱,并且T(λ)是光谱仪的传递函数。因此,如T(λ)是已知的,新样本R(λ)的正确校准的拉曼光谱由如下计算:
(方程式2),其中S(λ)是测量样本光谱。
传递函数T(λ)是光谱仪传递函数TS(λ)和探针传递函数TP(λ)两者的函数。方程式2可因此写为
(方程式3)。因为光纤探针是可更换的并且可以被消耗,所以将显而易见的是当具有新探针传递函数TP的新探针插入时,系统的总体传递函数将改变。
现在参考图3,主要的或主光谱仪在50处示出,并且次级的或从属光谱仪在51处示出。光谱仪50、51每个具有与在图1中示出的配置类似的配置,但可以具有不同的光纤探针和光谱特性。理想情况下,控制每个光谱仪的个人计算机35使用程序中的通用库以提供对该系统及数据处理的控制,并且主要和次级光谱仪50、51的特性一致的特性是令人满意的。在本示例中,主光谱仪50与主要或主探针52相关联,并且次级光谱仪51与在53a、53a、53b处示出的多个次级或从属探针相关联。在每一种情况下,参考在54处图解示出的标准荧光源来执行校准。
校准的第一种方法在图4a中示出。在步骤60处,根据主光谱仪对次级光谱仪进行波长校准。在这种情况下,次级光谱仪51的波长轴校准是通过,例如使用氩气汞光谱灯或具有确定光谱线的化学样本来执行,并且接下来执行使用线性内插的像素分辨率匹配以确保次级光谱仪的轴的尺寸与主光谱仪的轴的尺寸匹配。该校准的结果在图5中示出,其中来自主要和次级光谱仪50、51的光谱示出了精确对齐的灯的光谱线。在步骤61处,使用荧光源54对于第二光谱仪和探针53a执行校准。以与图2的图标类似的方式,光谱将从荧光源记录,然后计算传递函数以将测量的光谱带到具有已知光谱的线中,并且通过例如个人计算机35来存储光谱。在步骤62处,光谱仪51后续可用于体内拉曼测试或其它方式,并且测量的拉曼光谱可以使用在步骤61处记录的校准函数来校正。
当探针53a被丢弃并期望对新对象执行测试时,更换的探针53b可被取代,在这种情况下图4a的方法是重复的。
如在图4b和4c中示出的替代过程中,多个校准函数可以首先记录用于次级光谱仪和多个次级探针。在图4b中的步骤60处,如在图4a中,次级光谱仪51被校准,以使其与主光谱仪60一致。在步骤61处,测量用于次级探针53a的校准函数,并在步骤63处该校准函数以某种方式(例如,通过将校准函数保存为用与次级探针53a对应的参考标记来标签的计算机文件56a)来进行存储并且与探针53a相关联。如箭头64所示,这个过程然后被重复用于任何数目的探针53b,...,53n以提供探针的储存或储备。如在图4c中所示,当需要使用光谱仪51执行测试时,在步骤60处光谱仪根据如上所述的主光谱仪50来校准。在步骤65处,探针53n安装在系统上,并且相应的存储传递函数56n被检索。在步骤66处,使用次级光谱仪51的测试可被执行,并使用检索的校准函数56n来进行校准。
替代方法参考图6示出,其中从属或次级探针53a,…,53n在主要或主探针50上校准。根据方程式2,在主光谱仪采用具有传递函数TPP(λ)的主要或主探针和具有传递函数TSP(λ)的次级或从属探针测试时,来自荧光源F(λ)的光谱将产生用于主探针的光谱SPP(λ),其中
(方程式4)和使用次级探针的光谱SSP(λ),其中
(方程式5)。方程式4和5可相除,以通过探针校准函数TCF来关联两个探针传递值,其中
(方程式6)。因此,从方程式2和6,当次级光谱仪与次级探针一起使用时,测量的光谱S(λ)和拉曼光谱R(λ)由如下来相关
(方程式7),其中T(λ)=TS(λ)TPP(λ)是使用主探针用于次级光谱仪测量的存储系统传递函数。
如在图6至图7c中所示,这允许任何数量的次级或从属探针53a、53b、53n匹配到任何数量的次级光谱仪51a、51b、51n。如在图7a中所示,在第一步骤70处,采用与步骤60类似的方式,使用主探针52,按照主光谱仪50对次级光谱仪51a进行校准。以与图3至4c的方法类似的方式,通过测试相对荧光标准源54的次级光谱仪,在步骤72处发现系统传递函数71a。系统传递函数71a以任何适当的方式与相应的光谱仪51a相关联,例如在步骤73处在控制软件中或以其它方式。如由箭头74所示,这可以重复用于任何数量的次级光谱仪系统51b,...51n,以生成合适的系统传递函数71b,...71n。
如在图7b中所示,所述次级或从属探针53a、53b,...,53n相对主探针52来校准。在步骤75处,主光谱仪系统50采用相对荧光源54的主探针来适当地校准,尽管如果这一步已经执行并且与主探针相关联的传递函数已经存储,则该步骤可以省略。在步骤76处,主探针被探针53a替换,并且接下来,主光谱仪系统和相应的从属探针的组合相对荧光标准54进行检测。在步骤77处,校准函数TCF从主探针光谱和次级探针光谱的比率来计算,并且在78处这个被记录,并且与如在79a处示出的次级探针进行相关联的存储。如箭头80所示,这可以重复用于任何数量的次级探针53b,...53n,并且相应的校准函数TCF存储,如在79b,…79b处所示。
如在图7c中所示,次级光谱仪系统51n中的一个可与在次级探针53n中的任何一个一起使用,因为使用主探针52的系统传递函数71a是已知的,并且将次级探针53n与主探针52相关联的校准函数TCF是已知的。如在步骤81中所示的,次级光谱仪系统51n根据主光谱仪系统50来校准。在步骤82处,次级探针校准函数TCF被检索,并且存储系统的传递函数71n根据存储的校准函数TCF来修改。在步骤83处,体内拉曼测试或以其它方式可然后被执行,并且捕捉的拉曼光谱被校正。
因此在任何方法中,通过将次级光谱特性匹配回主光谱仪特性,存储用于光谱仪和探针组合的传递函数或存储用于合并了主探针的系统的传递函数,以及保存与刺激探针联用的校准函数,使用不同光谱仪和探针组合所捕捉的光谱将仍然一致并且可比较。
从图8和图9a至9c这是显然的。图8示出在主探针52和次级或从属探针53n之间不同的响应。强度响应随着光谱而改变,并且如图所示校准函数将次级探针的光谱映射到主要的或主探针的光谱。来自主光谱仪50和次级光谱仪51的未校准组织光谱在图9a中示出,并且在它们之间的差异是明显的。图9b和9c示出使用上述所述的方法中的每一种方法的校准结果,并且来自主要的和次级光谱仪的光谱基本上一致。
拉曼光谱从具有两种不同的探针(n=902光谱)的胃部测量。主成分分析(PCA)在次级探针的校准之前和之后进行。图10示出在光纤探针的校准之前和之后的PCA分析,以及对不同计分的95%置信区间。显而易见的是,在光纤探针的校准之后,光谱落在指示在主要和从属光纤拉曼探针之间的成功传递的相同置信区间内。
监测激光功率
图11是示出来自光纤探针的背景光谱(即在不存在组织信号的情况下)的曲线图。对应于拉曼散射或光纤二氧化硅内的荧光的峰和对应于远端球透镜的蓝宝石的峰是明显的。图12示出采用不同水平的发送功率从体内组织接收的拉曼光谱的曲线图。图11的峰在图12的不同线中是明显的,但将明显的是,峰的相对高度和连续背景随着发射功率而改变。
有利的是,已经发现在捕捉的拉曼光谱中的所述光纤探头和蓝宝石球透镜的光谱特性可以用作内部参考以导出发射的激光功率,而无需在光学系统中任何其它组件的提供。如在图13A的方法中所示,在步骤90处,在本示例352中适当大数量的光谱被收集,并且发射的激光功率被测量。在步骤91处,在本示例的偏最小二乘(“PLS”)回归中,执行适当的多变量统计分析。PLS回归将光谱数据的维数减到若干的潜在变量(“LV”)。在此情况下,光谱的变化和因变量之间的差异,所述激光功率被最大化,使得潜在变量将更高的权重给到与激光功率充分相关的光谱峰。通过选择适当数量的潜在变量,作为光谱特性函数的激光功率模型可以得出,并在步骤92处如图所示被存储。因此,在步骤93处如在附图13b中示出的操作中,测试光谱被捕捉,例如从体内对象或以其它方式被捕捉,并且在步骤94处,光谱值被提供给存储模型。在步骤95处,激光功率被导出并展示给操作者,例如在个人计算机35上显示给操作者。
在本示例中,所述包含的若干潜在变量相对均方误差的曲线图在图14a中示出,并且四个变量被选作给出误差和复杂度之间的最佳平衡。四个潜在变量和回归向量的相对负荷在图14b中示出。在图15中,示出了来自在五个对象中的166个光谱的实时测量的示例数据,其中测量的激光功率相对由模型预测功率进行绘制。将显而易见的是,大致的线性拟合表明预测所述的功率是实际传送功率的良好指标。
迭代背景减法
参考图16至18,其示出了减去从荧光、以及探针的二氧化硅和透镜的蓝宝石中的拉曼散射所产生的背景的方法。该背景信号对于每个特定光纤探针是唯一的。理想的是,从组织的拉曼光谱去除背景,而不过度或不足地减去背景。
如在图16中步骤110处所示,例如在不存在目标的情况下通过探针从激光源发送光,背景光谱被捕捉和存储。在步骤111处,从测试对象(例如从组织)接收拉曼光谱。在步骤112处,通过使用一个或多个不同的参考峰的强度来预测在测试对象拉曼光谱中的光纤背景信号量。在本示例中,峰可能是由于二氧化硅和/或蓝宝石(例如,417或490cm-1)。使用背景信号的预测量,所存储的背景信号可以乘以合适的,且可能与波长相关的校正因子,并且从测试光谱中减去(步骤113)。
在步骤114中,检查光谱剩余背景的存在。如果背景已经被完全去除(即,当二氧化硅和蓝宝石信号对于组织拉曼光谱的贡献可以忽略不计时),则光谱被传递以进行输出或进一步的分析,如在步骤115处所示。如果背景信号仍然存在,则步骤112至114如由箭头116所示的进行重复。
该方法不需要被限于单个二氧化硅/蓝宝石峰。多变量分析(例如,偏最小二乘和曲线解析方法等)也可以用于此目的。
举例来说,图17是示出从手掌组织接收的拉曼光谱和来自探针的背景光谱的曲线图。来自荧光、在探针的二氧化硅和球透镜的蓝宝石中的拉曼散射的峰是显而易见的,叠加在来自手掌的拉曼光谱上。该背景信号对于每个特定的光纤探针是唯一的。
如在图18中所示,在已经执行图16的迭代过程之后,平滑的拉曼光谱被示出,没有背景信号的独特的峰,但保持所需的基本拉曼光谱信息。
组合的系统
各种公开的方法可以一起使用。这种组合的一个实施例在图19中示出。在步骤100处,所述校准方法可以被执行,使得根据主系统或主要系统50的系统传递函数是已知的,并且后续的光谱可以被适当地校正。在步骤101处,可以执行信号的预处理,包括平滑和组织背景减除。在步骤103处,可对如上讨论的光谱进行功率监测,并且并行地,在步骤102处,可对探针背景进行减除。如在步骤104处所示,步骤102和103中的信息提供给在个人计算机35上的合适程序,以执行其它诊断或输出步骤。
在将所公开的校准方法与诊断方法进行组合的过程中,仪器独立的光纤拉曼光谱对于定量组织分析和表征是可能的。这允许由不同仪器采用的光谱和同样由相同仪器但是不同探针采用的光谱的比较。这对于诊断是重要的,因为它允许使用在不同机器上或者通过使用不同探针得到的光谱来用于比较。这对于提高诊断的精度是重要的。
实时癌症诊断
在线的生物医学光谱(即,反射率、荧光和拉曼光谱)系统和方法实现了临床内窥镜的实时癌症诊断,并且可以使用听觉反馈以及采用预测病理的概率诊断算法的结果的图形显示界面与临床医生相互交流。采用在胃部中的拉曼内窥镜作为示例(图22A);该方法能够预测几种病理:正常、肠化生、异常增生和肿瘤形成。这种在线诊断方法实时将信息提供给组织病理学的临床医生,其可用于做出决策,如活组织检查指导或肿瘤切除。包括GUI的该系统优化用于允许实时诊断(<0.1s),例如用于临床内窥镜的快速数据处理。
为了解决跨解剖和脏器间的光谱方差,在线框架实现了器官特异性诊断模型并且在不同的器官(例如,食道、胃、结肠、子宫颈、膀胱、肺、鼻咽、喉,和口腔(硬腭、软腭、颊、内唇、腹侧和舌))的光谱数据库之间切换。因此,所公开的拉曼平台是用于在内窥镜进行癌症诊断的通用诊断工具。
图20是用于根据一个实施例的实时癌症诊断的光谱采集和处理流程的诊断系统115的架构图。所述诊断系统115可以在个人计算机35上实施。诊断系统115包括光谱采集模块120、光谱预处理模块125、异常值分析模块130、多变量分析模块135、病理学模块140和数据库142。为简单起见,仅一个光谱采集模块120、光谱预处理模块125、异常分析模块130、多变量分析模块135、病理学模块140和数据库142示出,但在实践中对于其中的每种模块在操作过程中都可以设置为多个。
参考图20-22,在步骤145处,光谱采集模块120将激光激发源与CCD电子同步,并将合并读出(binned read-out)从CCD存储到数据库142中以用于进一步处理。通过缩放到基于前述组织测量的总光子计数的约85%内,所述光谱采集模块进一步自动调节曝光时间和光谱积累,而0.5秒的上限被设定为实现临床上可接受的条件。多光谱积累和曝光时间的自动调节提供了快速和简单的方法,以防止信号饱和,并获得用于内窥镜应用的高信噪比。如果光谱信号饱和,则该方法以减少的积分时间启动新的数据采集以防止饱和。在光谱采集之后,该方法识别并消除宇宙射线(例如,使用设定为最大阈值的整个光谱范围内具有95%置信区间(CI)光谱的一阶导数)。所识别的宇宙射线通过线性内插去除。短光谱采集时间帧对于内窥镜应用是特别有用的。在图22A处示出的GUI示出在180处的光谱采集。
对于其它应用,所述光谱采集框架也可以用于外部或内部外科手术介入或评估外科手术期间的组织类型。该实时功能允许在现场诊断,并且因此可用于指导肿瘤切除术的切除范围。所述诊断信息可以在线给出(即<0.5秒)以辅助医疗决策,这是很关键的。对于皮肤测量,有较宽松的测量时间的要求,因为皮肤光谱在更可控的实验设置下以更长的曝光时间的概率来采集。在线软件架构还可以适用于需要快速光谱测量(包括荧光、反射光谱或在如过程分析技术、食品科学、取证等不同的领域)的其它领域,从而不间断的实时筛选是必要的。在步骤150中,光谱采集模块120确定信号是否已达到饱和。如果是这样,它以减少的积分时间启动新数据采集以防止饱和。在步骤155中,所述光谱采集模块120识别并消除宇宙射线(即,使用设定为最大阈值的整个光谱范围内具有95%置信区间(CI)光谱的一阶导数)。在一个实施例中,识别的宇宙射线通过线性内插去除。宇宙射线可以通过其它方法,包括多变量分析、平滑、均值滤波、中值滤波、傅立叶变换、小波分析等去除。
在步骤160中,所述光谱预处理模块125以积分时间和激光功率缩放采集的光谱。一阶Savitzky-Golay平滑滤波器进一步用于去除在强度校正光谱中的噪声。然后,限制平滑光谱的下边界的5阶修改的多项式被减去,以解析单独组织的拉曼光谱。拉曼光谱在从800cm-1到1800cm-1的曲线下最后归一化到积分区域以解析光谱的线形状和相对强度,减少在临床内窥镜处的探针处理变化。GUI(图22A)示出在185处的归一化光谱。在一些实施例中,光谱预处理模块125利用用于预处理的附加方法,包括但不限于,多重散射校正(MSC)、FIR滤波、加权基线减法、降噪、平均中心化(mean centering)、差分(differentiation)等。
在步骤165中,异常分析模块130使用与Hotelling’s T2和Q-残差统计耦合的主成分分析(principal component,PCA)来检测异常光谱。GUI(图22A)示出在190处的异常分析。异常分析的实现用作在使用与Hotelling’s T2和Q-残差统计耦合的主成分分析(PCA)的在线框架中的高度模型特异性反馈工具。Hotelling’s T2和Q-残差是提供在拟合模型内部和外部的信息的两个独立参数。使用这些参数作为光谱质量(即,探针接触模式、混杂因素、白光干涉等)的指标,听觉反馈集成到促进实时的光谱筛选和用于临床医生的探针处理建议的在线拉曼诊断系统中。软件系统提供用于不同诊断结果的不同声音反馈。例如,如果光谱是异常值,则一定的声音将出现。如果光谱被诊断分类为“正常”,则第二种不同的声音将出现。如果光谱被分类为“癌前病变”或“癌症”,则第三或第四声音将出现。声音频率可以与“后验概率”成比例。这是非常有用的,因为它在接收诊断信息时,提供实时指导给内窥镜医生。因此,内窥镜医师并不需要关注拉曼平台监测器,而是将重点放在具有声音指导的内窥镜操作步骤上。如果异常分析模块确定采集的光谱是异常值,则诊断系统115在步骤145处重新开始。
如果该光谱被验证用于进一步分析,则它们被馈送到用于体内癌症诊断的概率模型。在步骤170中,多变量分析模块135应用用于体内癌症诊断的概率模型。基于大量患者的光谱数据库,多变量分析模块135在不同的预定义模型之间切换,包括偏最小二乘-判别分析(PLS-DA)、PCA-线性判别分析(LDA)、蚁群算法(ACO)-LDA、分类和回归树(CART)、支持向量机(SVM)、自适应增强(AdaBoost)等。
在步骤175中,病理模块140实现可以在概率癌症诊断的不同器官的光谱数据库之间切换的器官特异性诊断模型。除了音频反馈,GUI(图22A)提供由在195处的病理模块140的输出给临床医生。
图22B提供根据第二实施例的GUI。
所述数据库142存储采集的光谱以及用于诊断的存储光谱。
在一些实施例中,多个光谱被采集并进行分析。例如在5-15之间采集。每个被分析,并且如果超过阈值百分比则提供相同的输出:癌症VS正常,即所确定的诊断。例如,如果10个光谱被采集,并且7个或更多提供相同的答案,则这是诊断。如果只有5个或6个提供相同的答案,则该过程被重复。
实施例
用于实时诊断的集成式拉曼光谱和三峰宽视场成像系统包括与USB 6501数字I/O(美国国家仪器公司,美国德克萨斯州奥斯汀)电子同步的光谱稳定的785nm的二极管激光器(最大输出:300mW,B&W TEK公司,美国特拉华州纽瓦克市),配有液氮冷却、NIR优化、背照和深耗尽电荷耦合器件(CCD)摄像机(在每个像素20μm×20μm处1340×400像素;规格-10:400BR/LN,美国新泽西州Trenton的Princetion仪器公司)的透射成像光谱仪(Holospec,f/1.8,美国密歇根州Ann Arbor的Karser光学系统公司),以及用于激光递送和体内组织拉曼信号收集的专门设计的拉曼内窥镜探针。1.8mm的拉曼内窥镜探针由32个收集光纤组成,所述收集光纤围绕中央光递送光纤,且该中央光递送光纤在探针的近端和远端具有用于最大化收集组织拉曼信号的光学滤波的两级,同时减少瑞利散射光、光纤荧光和二氧化硅拉曼信号的干扰。拉曼探针可以很容易地向下传递给医疗内窥镜的仪器通道,并且在宽视场内窥镜成像(WLR/AFI/NBI)模态的指导下,引导到可疑的组织部位。使用具有约9cm-1的光谱分辨率的1.5W/cm2(“光点尺寸”(spot size)为200μm)的785nm激发功率,系统采集在0.5s内来自体内上消化道组织的800cm-1-1800cm-1波数范围中的拉曼光谱。
拉曼系统的硬件组件(例如,激光功率控制装置、光谱仪、CCD快门和摄像头读取同步)通过用于不同光谱仪/摄像机(例如,PVCAM库(美国新泽西州Trenton的Princeton仪器,Roper科学有限公司)和OmniDriver(美国佛罗里达州Dunedin的海洋光学公司)等)接入到Matlab软件。激光与CCD快门电子同步。激光功率、曝光时间和光谱积累的自动调节是通过缩放到基于前述组织拉曼测量的总光子数的85%(例如65000个中55250的光子)来实现的,而0.5秒的上限被设定为实现临床可接受的条件。多光谱的积累和曝光时间的自动调节提供了一种快速和简单的方法,以防止CCD饱和并获得用于内窥镜应用的高信噪比(SNR)。拉曼位移轴(波长)使用水银/氩气校准灯(美国佛罗里达州Dunedin的海洋光学公司)来校准。用于系统的波长依存性的光谱响应校正使用标准灯(RS-10,美国加利福尼亚州圣地亚哥的EG&G伽玛科学公司)进行。该平台的再现性采用作为波数标准的环己烷和对乙酰氨基酚的激光频率和拉曼光谱进行连续地监测。包括CCD温度、积分时间、激光功率、CCD对准的所有系统性能测量因此经由SQL服务器登陆到中央数据库。
拉曼信号的实时预处理通过使用设定为最大阈值的整个光谱范围内具有95%置信区间(CI)的一阶导数,以宇宙射线的快速检测来实现。位于阈值外侧的数据点内插到2阶。光谱进一步用积分时间和激光功率缩放。一阶5点的Savitzky-Golay平滑滤波器用于去除在强度校正的光谱中的噪声,而限制到平滑光谱下界的5阶修改后多项式被减去以解析单独组织的拉曼光谱。拉曼光谱被归一化到从800cm-1到1800cm-1曲线下的积分区域,使得能够在不同组织病理学之间的光谱形状和相对拉曼基带强度进行更好的比较。光谱随后根据具体的数据库局部地平均值中心化以去除在数据库中的常见变化。预处理之后,拉曼光谱馈送到模型特定的异常分析。
通过使用耦合到Hotelling’s T2与Q-残差统计的PCA,异常检测方案被结合到生物医学光谱中作为在线框架中的高度模型特定反馈工具。通过将它们分解到正交分量(主成分(PCs))的线性组合中,PCA减少了拉曼光谱的维度,使得在数据集中的光谱变化被最大化。数据矩阵X的PCA模型被定义为:
X=TPT+E
其中T和P表示计分和负荷,并且E包含残差。负荷对应于新的旋转轴,而计分表示数据投影值。因此,Hotelling’s T2统计值是由PCA模型(样本到模型的距离)捕捉的方差测量,并且由如下定义:
其中tik是用于使用分量k的第i个样本光谱的PC计分,并且是用于分量k的协方差矩阵的归一化特征值的对角线矩阵。因此,Hotelling’s T2给出在PCA模型内的极值指示。另一方面,Q残差是没有由PCA模型捕捉(缺少模型拟合统计)的方差测量,并且由如下定义:
其中,xi是样本光谱,Qik是用于使用分量k的第i个样本光谱的平方重建误差的总和,并且Pk是PC负荷。对于Hotelling’s T2与Q-残差两者,标准化的99%CI被用作上阈值以截断异常的拉曼光谱。因此,Hotelling’s T2与Q-残差是两种提供关于模型拟合的定量信息的独立参数。使用这些参数作为光谱质量(即,探针接触模式、混杂因素、白光干扰等)的指标,听觉反馈已结合到在线拉曼诊断系统中,有利于实时探针处理建议和在临床内窥镜程序期间用于临床医生的内窥镜筛选。
在组织的拉曼光谱质量的验证之后,那些合格的拉曼光谱立即馈入用于在线体内诊断和病理预测的概率模型。该GUI可以在不同的模式之间快速切换,包括偏最小二乘-判别分析(PLS-DA)、PCA-线性判别分析(LDA)、蚁群优化(ACO)-LDA、分类和回归树(CART)、支持向量机(SVM)、自适应增强(AdaBoost)等,用于在临床内窥镜程序处的预测分类。作为例子,概率的PLS-DA被用于胃癌诊断。PLS-DA采用PCA的基本原理,但是进一步通过将光谱变化和组关联性之间的协方差最大化来进一步旋转成分,以获得与诊断有关的变化,而不是在光谱数据集中的最显著的变化。该系统支持库分类,一对全部(one-against-all)和一对一多类(one-against-one multiclass)(即良性、异型增生和癌症)的概率PLS-DA判别分析以预测特定的组织病理。
实施例1
总计2748个体内胃组织光谱(2465个正常,和283个癌症)从被招募以构建用于开发用于胃癌诊断的诊断算法的光谱数据库的305个患者中采集。组织病理被作为体内组织诊断和表征的拉曼技术的性能评估的金标准。
胃呈现有用于光谱诊断的许多混杂因素(即胃液、食物残渣、出血、渗出液等)的最有挑战性的器官中的一个。从305个胃部病人(正常(n=2465)和癌症(n=283))中采集的用于算法开发的体内平均拉曼光谱在图23中示出。胃组织的拉曼光谱示出在875cm-1(羟基脯氨酸的v(C-C))、936cm-1(蛋白质的v(C-C))、1004cm-1(苯丙氨酸的环呼吸Vs(C-C))、1078cm-1(脂类的v(C-C))、1265cm-1(蛋白质的酰胺III v(C-N)和δ(N-H))、1302cm-1和1335cm-1(蛋白质和脂类的δ(CH2)变形)、1445cm-1(蛋白质和脂类的δ(CH2))、1618cm-1(卟啉的v(C=C))、1652cm-1(蛋白质的酰胺I v(C=O))和1745cm-1(脂类的v(C=O))。胃组织拉曼光谱包含来自可能反映在胃壁中皮下脂肪的甘油三酯的巨大贡献(即在1078cm-1、1302cm-1、1445cm-1、1652cm-1和1745cm-1处主要的峰)。胃癌的拉曼光谱揭示在上述拉曼光谱性质(例如,强度、光谱形状、带宽和峰值位置)中的显著变化,再次确认了在体内拉曼研究中我们前述的内容。
自动异常检测使用具有Hotelling’s T2与Q-残差(99%CI)的PCA实现了可预测在线分析。为了使在线诊断有效,包括最大组织的光谱变化的双组分PCA模型被定义。在数据集(n=2748个拉曼光谱)中总变异性的38.71%(PC1:30.33%、PC2:8.38%)最大方差所占的这些选择的显著PC(P<0.0001)和相应的PC负荷在图24中示出。
图25示出了用于例示用于将癌症光谱与正常光谱分开的PC计分功能的正常(n=2465)和癌症组织光谱(n=283)的计分散点图(即PC1与PC2)。Hotelling’s T2与Q-残差99%的CI因此从训练数据集计算,并且固定为用于预期的在线光谱验证的阈值。我们然后定义用于胃癌预测的概率PLS-DA模型。训练数据库随机地重新采样多次(n=10)成学习集(80%)和测试集(20%)。生成的PLS-DA模型分别提供用于胃癌诊断的85.6%(95%CI:82.9%-88.2%)(80.5%(95%CI:71.4%-89.6%)的灵敏度和86.2%(95%CI:83.6%-88.7%)的特异性)。然后,我们进一步测试异常检测以及在10个准胃癌患者中的概率PLS-DA。用于预期正常(n=45)和癌症(n=30)组织光谱的PC计分散点图(即,PC1VS PC2)同样在图25中示出。
图26示出用于从10个预期胃部样本采集的105个(45个正常,30个癌症,30个异常)光谱、具有99%CI边界的Hotelling’s T2(38.71%)与Q-残差(61.29%)的散点图。虚线表示验证预期拉曼光谱是否处于主成分分析(PCA)模型的常见组织变化内的99%的置信区间(CI)。可以观察到,大量的非接触光谱位于99%的CI外侧,并因此被实时丢弃,无需用于组织诊断。被验证的组织拉曼光谱大量落入T2与Q-残差的99%CI内侧,表明该在线数据分析提供快速并且高效的生物医学组织光谱的实时验证手段。
由在线异常分析验证的预期采集光谱进一步馈入概率PLS-DA用于即时疾病预测,实现用于胃癌检测(图27)的80.0%(60/75)的诊断精度,如由组织病理学检查所证实的。单独的虚线给出了用于从正常胃组织体内分离癌瘤的90.0%(27/30)的诊断灵敏度和73.3%(33/45)的特异性。
受试者工作特征曲线(ROC)被进一步生成以评估组分离。图28示出从用于回顾性预测的光谱数据库的每个随机分割计算的ROC曲线以及用于预期数据集预测所计算的ROC。在用于回顾性和预测数据集所生成的ROC曲线下的积分区域分别是0.90和0.92,示出用于体内胃癌诊断的PLS-DA算法的稳定性。
用于对组织病理预测的全部前述在线数据采集的总处理时间是0.13秒。用于在图21中流程表的每个步骤的处理时间在表1中给出。可以实现自由运行光学诊断和处理时间<0.5s,这对于在内窥镜处实现实时的体内组织诊断是至关重要的。
表1.用于在具有64位I7四核4GB内存的个人计算机上在线生物医学拉曼光谱框架的平均处理时间
实施例2
拉曼光谱系统包括光谱稳定的785nm二极管激光器(最大输出:300mW,美国特拉华州Newark的B&W TEK公司),配有液氮冷却、背照和深耗尽CCD摄像机(在每个像素20μm×20μm处1340×400像素;规格-10:400BR/LN,美国新泽西州Trenton的Princetion仪器公司)的透射成像光谱仪(Holospec,f/1.8,美国密歇根州Ann Arbor的Karser光学系统公司)。该系统还由特别设计的熔融二氧化硅光纤拉曼内窥镜探针(外径1.8mm,长度1.30米)组成,该内窥镜探针包括围绕中央光传送光纤(直径200μm,N.A.=0.22)的9×200μm的收集光纤(N.A.=0.22)。1.0mm的蓝宝石球透镜(折射率1.76)被耦合到用于增强上皮组织拉曼测量的拉曼探针的光纤探头。所述系统用9cm-1的光谱分辨率采集在800cm-1至1800cm-1范围内的拉曼光谱。在这项研究中每个拉曼光谱在785nm的激光激发下用0.5s的积分时间测量。快速拉曼光谱技术使用氩气/水银光谱灯(AR-1和HG-1,美国佛罗里达州Dunedin的海洋光学公司)来进行波长校准。所有波长校准光谱使用钨-卤校准灯(RS-10,美国加利福尼亚州圣地亚哥的EG&G伽玛科学公司)对于系统的强度响应进行校正。
使用在图20-22中描述的系统和方法用于控制用于实时数据采集和分析的拉曼光谱系统。从体内组织测量的原始拉曼光谱表示弱拉曼信号、强烈的自体荧光背景和噪声的组合。原始光谱由一阶Savitzky-Goly平滑滤波器进行预处理(选择具有3个像素的窗口宽度来匹配光谱分辨率)以减少光谱噪声。在指纹区域(800cm-1-1800cm-1)中,五阶多项式被发现对于拟合在噪音平衡光谱中的自体荧光背景最优,并且这个多项式随后从原始光谱中减去以产生单独组织的拉曼光谱。所有的上述预处理在100ms内完成,并且处理结果可在计算机屏幕上实时显示。
PLS回归被用作多变量方法以从光纤拉曼探针提取特征内部参考的背景信号。简单地说,PLS采用PCA的基本原理,但通过将光谱变化和因变量(例如,激光激发功率)之间的协方差最大化,来进一步旋转成分LV,使得LV负荷解释相关变化,而不是在光谱数据集中的最突出的变化。与激光激发功率相关联的重要光谱参考信号被保留在最初的几个LV中。在这项研究中,在建模之前执行平均中心化以减小PLS回归模型的复杂性。PLS回归模型的最优复杂性通过留一交叉验证确定,并且PLS回归模型的性能通过计算确定(R2)的系数、校正均方根误差(RMSEC)、交叉验证均方根误差(RMSECV)和预测均方根误差(RMSEP)来检查。需要注意的是,最优PLS模型具有高的R2,但具有低的RMSEC、RMSECV和RMSEP。在这项研究中用于解析参考信号开发的PLS回归模型也实现为在我们的实时临床拉曼软件中的在线激光激发功率预测器,并且以无偏差的方式(unbiased manner)预期测试。多变量统计分析在Matlab(马萨诸塞州纳提克的MathWorks公司)程序环境中进行。
总计30名正常健康的对象(16名女性和14名男性)被招募用于口腔中的体内组织的拉曼光谱测量。在体内组织的拉曼光谱测量之前,所有对象进行进行充分漱口,以减少混杂因素(如食物残渣、细菌包覆等)。体内组织的拉曼光谱(n=783)收集有25个对象的内唇。对于25个对象,体内口腔组织的拉曼光谱(n=~5)在5mW-65mW(约10mW的间隔)的范围中在六个功率水平处收集。在每一个组织的拉曼测量之前,使用功率计以±0.5%的线性度和±3%(0.1mW至100mW的范围)的精度,激光的激发功率水平在光纤探针的远端头处测量。其它混杂因素(例如在组织表面上的探针压力、光漂白、组织的光学性能和光纤探针的弯曲)没有特意监测,但结合到用于在原位的参考信号的鲁棒提取的PLS建模中。在在线拉曼采集框架中部署开发的PLS模型之后,用于激光激发功率监测的内部参考信号的预期和独立验证在5个新对象(n=166个光谱)上实时执行。
为了进一步验证在这项工作中开发的内部参考方法的定量值,我们还进行了组织模型实验。各种明胶浓度的组织模型从牛的皮肤(B型明胶(美国Sigma公司G9391))制备。明胶以预定浓度(以重量计20、25、30、35、40、45和50%)溶解在蒸馏水中。溶解的明胶在水浴中采用连续搅拌加热至50℃1小时。接着,将熔融明胶倒入一预冷的模具(4℃)中,并存储2-3小时以产生固体明胶模型。然后,执行组织模型的定量光纤拉曼光谱分析。使用光纤拉曼探针,以不同的激光功率从各个组织模型测量得到总数n=133的拉曼光谱。激光激发功率在10-60mW的范围中改变,并且测量得到的光谱归一化到如由内部参考方法预测的激光泵浦功率。
图11示出当由785nm的二极管激光器激发时,所使用的球透镜光纤拉曼探针的背景光谱。源自远端球透镜的不同蓝宝石(AL2O3)拉曼峰可以在417cm-1和646cm-1(以A1g对称的声子模式)以及380cm-1(和751cm-1(Eg声子模式)处被发现。从熔融二氧化硅光纤以及相对较弱的光纤荧光背景存在两个主导拉曼分量。在490cm-1和606cm-1处表示为D1和D2的熔融二氧化硅的尖锐“缺陷峰”已经分别分配在四元环和三元环中的氧原子的呼吸振动。与无定形二氧化硅物质的一般特征相关联的强烈玻色子拉曼频带的肩峰(shoulder)(约130cm-1)同样从光纤拉曼探针的背景光谱观察。二氧化硅玻色子频带在约60cm-1附近达到峰值,但由于我们的拉曼探针设计的光滤波,仅肩峰是明显的。来自光纤拉曼探针本身的这些特征背景拉曼峰(比指纹区更短(800cm-1至1800cm-1))可以用作用于体内组织拉曼测量的内部参考信号。
为了开发PLS回归模型并解析内部参考信号,我们采用作为独立参数的激光激发功率测量在口腔中25个对象体内的拉曼光谱。对于每一个对象,体内组织拉曼光谱(n=~5)以在5mW-65mW(约10mW的间隔)范围中的不同功率水平来采集。图12示出使用不同的激光激发功率(例如10mW、30mW和60mW)从內唇测量的平均体内原始拉曼光谱±1标准偏差(SD)的示例。可观察到在方差的广阔自体荧光背景上叠加的弱组织拉曼信号。图29示出校准的无背景平均拉曼光谱±1SD。内唇缘的体内拉曼光谱示出在约853cm-1(v(C-C))、1004cm-1(Vs(C-C))、1245cm-1(蛋白质的酰胺III v(C-N)和δ(N-H)、1302cm-1(CH3CH2扭曲和摇摆)、1443cm-1(δ(CH2)变形)、1655cm-1(蛋白质的酰胺I v(C=O))和1745cm-1(v(C=O))。另一方面,原始的体内组织的拉曼光谱(图12)也包含来自光纤拉曼探针的突出熔融二氧化硅和蓝宝石拉曼峰,即:380cm-1、417cm-1、490cm-1、606cm-1、646cm-1和751cm-1。
一种从口腔组织的拉曼光谱中提取广泛范围特性的内部参考峰的PLS回归模型。瑞利散射光从PLS分析排除。测量的体内原始组织拉曼光谱被设置在具有逐行光谱和逐列波数的矩阵中。参考激光功率水平被设置在表示因变量的列向量中。在平均中心化之后,使用留一交叉验证开发PLS回归模型,以便建立用于定义用于激光泵浦功率预测的鲁棒参考信号的最优算法。图14a示出作为保留LV函数的激光功率预测的RMSEC和RMSECV。PLS回归分析示出,最优模型(RMSECV=2.5mW)可以使用4个LV获得。图14b显示占最大拉曼光谱方差(即LV1:94.8%、LV2:3.0%、LV3:0.9%和LV4:0.2%)和激光激发功率方差(LV1:80.1%、LV2:16.8%、LV3:0.8%、LV4:0.7%)的前四个LV负荷。还示出了所计算出的PLS回归向量。图30a使用留一交叉验证示出体内激光功率的监测结果(即测量的激光功率vs预测的激光功率)。数据可以通过指示相当的线性关系(R=0.98)的方程(y=0.551+0.984x)来拟合。4个LV的PLS模型的复杂性提供准确的内部参考用于采用2.5mW的RMSECV和0.981的R2的激光激发功率。相同的PLS回归模型随后在拉曼软件中实时在线实施,用于5个新对象(n=166个光谱)的独立验证。图30b示出使用开发的PLS回归模型在测量的实际激光激发功率和预测的激光激发功率之间的关系。2.4mW的RMSEP和线性关系(y=0.342+1.011x;R2=0.985)可获得,再次确认在体内组织拉曼测量期间作为内部参考方法的PLS回归应用。
内部参考方法的定量值被开发用于组织模型定量光谱分析。具有不同浓度(即以重量计20%、25%、30%、35%、40%、45%和50%)的明胶的七种组织模型被构造和测试。来自明胶模型的拉曼光谱(n=133个光谱)被测量并归一化到实时预测的激光功率。图31示出在60mW激发激光功率处从具有不同浓度的明胶组织模型测量的拉曼光谱。正如预期的那样,这些拉曼光谱示出拉曼峰强度和明胶浓度之间的线性关系(R=0.992)。图32示出随着变化的激发激光功率(从10mW到60mW)在实际的明胶浓度和预测浓度之间的相关性。显而易见的是,通过原位实时的激光激发功率来校正激光功率,明胶组织模型的准确定量分析可实现(RMSEP=1.9%和R2=0.985)。上述结果表明,基于多变量的内部参考信号,开发的实时功率监测方法可以实现在光纤组织拉曼光谱中的鲁棒定量组成分析。
实施例3-体内实时经鼻的图像引导拉曼内窥镜:定义在鼻咽部和喉部的光谱特性
这项研究表明在经鼻内窥镜应用中的拉曼光谱中的可行性,提供在头部和颈部的大规模临床研究的基础。采用开发的小型化光纤拉曼探针集成的图像引导的拉曼内窥镜平台提供了在临床内窥镜检查期间在分子水平处的头部和颈部的内源性组织成分的迅速和微创评估。这显著促进了临床医生来获得在头部和颈部的组织的详细生物分子指纹,反映真正的组成和形态特征,而不会引入由血管穿刺或组织脱水、形态和解剖效应等导致的人为现象。
拉曼光谱系统由光谱稳定的785nm的二极管激光器(最大输出:300mW,美国特拉华州Newark的B&W TEK公司),配有液氮冷却、NIR优化、背照和深耗尽电荷耦合器件(CCD)摄像机(在每个像素20μm×20μm处1340×400像素;规格-10:400BR/LN,美国新泽西州Trenton的Princetion仪器公司)的透射成像光谱仪(Holospec,f/1.8,Karser光学系统公司)组成。新颖的光谱仪光纤输入耦合件由58个光纤(100μm)的抛物线对齐阵列构成,以校正用于改进拉曼信号的光谱分辨率和信噪比两者的光谱仪像差。利用1.8mm光纤拉曼探针进行经鼻内窥镜应用,可以使得组织激发和体内组织拉曼收集两者均最大化。拉曼光纤探针拟合到柔性经鼻内窥镜的仪器通道,并且可以在宽视场成像(即,白光反射率(WLR)和窄带成像(NBI))引导下安全地引导到在鼻咽和喉部中的不同位置。临床拉曼内窥镜平台已经集成有我们最近开发的在线数据处理软件,以方便实时(处理时间<0.1s)探针处理提醒和声音反馈给临床医生。简单地说,在线拉曼内窥镜框架将光谱采集(即,激光曝光、积分时间、CCD快门和读出等)同步,并使用包括平滑、五阶多项式基线减法等建立的预处理方法,从体内的拉曼光谱(包括强的自体荧光背景和弱拉曼信号)自动地提取原始组织光谱。体内的拉曼信号和多变量算法的结果(例如,主成分分析)可以在临床经鼻拉曼内窥镜期间实时显示在易于理解的图形用户界面(GUI)中。
总共有不同的种族(22个亚洲人和一个白人)的23例健康男性对象被招募用于以经鼻内窥镜的体内组织拉曼光谱测量。在这些对象招募中,在WLR和NB成像检查下没有可疑病灶被识别。共有假定正常(或良性)的组织的三个主要测量位点被预定义以用于体内拉曼采集,包括真正的喉声带(LVC)、后鼻咽(PN)和此外其中NPC通常启动的咽隐窝(即咽隐窝(FOR))。光纤拉曼探头可以放置与实时询问内源生物分子组合物的内部组织轻柔接触。对活组织检查部位的准确定位由主管的内窥镜医师对WLR/NBI监视器验证。探针允许拉曼光谱从具有约1mm3和约800μm的穿透深度的区域(直径200μm)收集。使用具有在组织表面上约50mW功率的785nm激光,每个光谱在0.5s内采集。
拉曼光谱在线显示并且储存用于后程序检查。这种快速拉曼内窥镜技术是非破坏性的,并且现在可以常规地在用于临床评价的经鼻内窥镜检查中使用。为了评估组织内部位的方差,几个拉曼光谱(约18)也从每个组织部位获得。因此,来自47个部位的总计874个体内拉曼光谱在经鼻内窥镜处测量并且用于来自23个对象的光谱分析[PN(n=521),FOR(n=157)和LVC(n=196)。
在数据分析之前,原始拉曼光谱首先使用线性Savitzky Golay滤波器平滑,并且使用5阶多项式拟合将组织自体荧光背景从平滑光谱中减去。背景减去的拉曼光谱在曲线下归一化到积分区域,以相对于不同的对象和组织部位,最大限度地减少对临床拉曼测量的拉曼探针的处理变化的效果。所有处理的拉曼光谱组成矩阵,并且然后执行整个拉曼数据集的平均中心化。为了减少光谱数据的维度,主成分分析(PCA)被用来提取说明在用于组织表征的拉曼光谱数据集中的最大方差的一组正交主成分(PC)。因此,在PC上的负荷表示在占逐渐减小的方差的数据集中最突出光谱变化的正交基础光谱,而在PC上的计分表示在相应负荷上的组织拉曼光谱的投影值。因此,PCA可以有效地用于解析光谱方差,同时降低数据集的维度为最小。基于方差(ANOVA)和在0.05水平处t分布检验(Student’s t-test)的分析来选择保留PC的数量。我们采用事后(post-hoc)Fisher’s的最小二乘差异(LSD)的测试来评估平均值中的差。使用在Matlab(马萨诸塞州纳提克的MathWorks公司)编程环境中的PLS工具箱(华盛顿州Wenatchee的特征向量研究公司)执行多变量统计分析。
可常规地在经鼻图像引导(即,WLR和NBI)内窥镜检查期间在鼻咽和喉部中采集高质量的体内拉曼光谱。图1示出在内窥镜处以0.1s的采集时间从后鼻咽采集的体内原始拉曼光谱(叠加在大组织自体荧光背景上的弱拉曼信号)的示例。具有>10(图33的插图)的信噪比(SNR)的背景减去的组织拉曼光谱可获得并且在临床内窥镜测量期间在线显示。图34示出当拉曼探针在WLR/NB成像指导下与组织轻柔接触时正常的鼻咽[PN(n=521)和FOR(n=157)]和喉组织[LVC(n=196)]的主体间体内平均拉曼光谱±1个标准差(SD)。收集的鼻咽和喉组织拉曼光谱的比较(图34)展示出:在体液中的那些生化物质没有显著贡献于在经鼻内窥镜处的体内组织拉曼光谱。如所示出的是从相应的解剖位置采集的WLR图像。与蛋白质和脂质相关联的突出拉曼频带以暂时的生物分子分配被识别,如在表2中列出的。
表2 参与鼻咽和喉组织的拉曼散射的分子振动和生化物质的暂定分配(Tentativeassignment)(其中v是拉伸模式;Vs是对称的拉伸模式;δ是弯曲模式)。
图35示出随机选择的对象的主体内平均光谱±1SD。体内组织的拉曼光谱被认为可以以在鼻咽和喉中的微小的主体间和主体内的方差(<10%)再现。进一步拉曼内窥镜测试表明,后鼻咽内的不同组织部位之间的变异性是微小的(<5%)(数据未示出)。我们还计算如在图36中示出的不同组织类型(即,PN-LVC、LV-FOR和PN-FOR)之间的计算差光谱±1SD,解析在生物分子层面不同的解剖组织部位的独特成分和形态分布。ANOVA揭示以812cm-1、875cm-1、948cm-1、986cm-1、1026cm-1、1112cm-1、1254cm-1、1340cm-1、1450cm-1、1558cm-1、1655cm-1和1745cm-1为中心的三个解剖组织部位之间显著变异的十二个显著和广泛的拉曼光谱次级区域,再次确认表征朝向准确的体内组织诊断的鼻咽和喉的拉曼光谱特性的重要性。
在从健康志愿者获得的血液、唾液和鼻涕的体内拉曼光谱进行测量,如图37所示。在唾液和鼻涕中最突出的拉曼频带是在1638cm-1、(v2是水的弯曲模式),而血液显示在1560cm-1、和1620cm-1附近的卟啉拉曼频带。为进一步评价在头部和颈部中不同组织之间的光谱差,基于占总方差的57.41%的ANOVA和t分布检验(Student’s t-test)(P<0.05)(PC1:22.86%;PC2:16.16%;PC3:8.13%;PC4:6.22;PC5:4.05%)的五组分PCA模型已开发来解析不同解剖位置的显著峰值方差。图38示出揭示与蛋白质(即853cm-1、940cm-1、1004cm-1、1265cm-1、1450cm-1和1660cm-1)和脂质(即,1078cm-1、1302cm-1、1440cm-1、1655cm-1和1745cm-1)相关联的拉曼频带。图39(A至E)显示不同组织类型(即PN、FOR和LVC)的PCA计分的图表。在每个槽口箱内的线表示中位数,但箱子的下边界和上边界分别指示第一(25.0%百分位数)和第三(75.0%百分位数)四分位数。误差条(须状物)表示1.5倍四分范围。p值同样在不同的组织类型之间给出。集成线性判别分析(LDA)的二分法PCA算法分化使用留一交叉验证提供了分别用于在PNvs FOR与LVCvsPN之间差的77.0%(401/521)、67.3%(132/192)的灵敏度和89.2%(140/157)和76.0%(396/521)的特异性。总的来说,这些结果表明,在头部和颈部中的鼻咽和喉的拉曼光谱可以在经鼻内窥镜处体内测量,并且诊断算法的发展应是组织部位特异性的,以确保最小的算法复杂性。
实施例4–用于在Barrett食管中异型增生的实时体内诊断的光纤共焦拉曼光谱仪
光纤共焦拉曼诊断可以实时实现(<0.5秒),并揭示原位Barrett癌变中上皮细胞和组织的渐进生物分子和功能变化。组织病理学表征:预测的已测量的组织部位中的152个为柱状上皮层(n=597个光谱),48个为肠上皮化生(n=123个光谱),9个高度异常增生(n=77个光谱)。使用受试者工作特征曲线(ROC)分析,可以成功实现对高度异常增生的识别,产生在光谱基础上的87.0%的灵敏度和84.7%的特异性。ROC曲线下的面积被认为是0.90。这种新的具有实时功能的生物分子特异性内窥镜方法为肠胃病学家们在正在进行的内窥镜过程中,客观定位Barrett病人体内的高风险组织提供了一种可靠的工具。
共焦拉曼光谱系统由近红外(NIR)激光二极管(λex=785nm),高通量透射成像光谱仪和专门设计的1.8mm光纤共焦拉曼探针组成;所述高通量投射成像光谱仪配备有液氮冷却的以及NIR优化电荷耦合器件(CCD)摄像机。该系统以约9cm-1的分辨率在范围800cm-1-1800cm-1中采集拉曼光谱。所开发的光纤共焦拉曼内窥镜探针用于激光递送以及体内组织拉曼信号收集。
1.8mm(外直径)的共焦拉曼内窥镜探针包括围绕中央光传送光纤(直径200μm,NA=0.22)的9×200μm滤波器涂覆的收集光纤(NA=0.22)。微型的1.0mm蓝宝石球透镜(NA=1.78)被耦合到共焦探针的光纤探头以将激发光紧密地聚焦到组织上,使来自上皮层(<200μm)的有效拉曼光谱收集成为可行的。所述光纤共焦拉曼探针可以插入到常规内窥镜的仪器通道中,并放置与用于体内组织表征和诊断的上皮轻柔接触。这种共焦拉曼探针的深度选择性提供了令人信服的实验优点,包括:(i)光纤共焦拉曼光谱可以选择性地定位与早发型Barrett’s癌变相关联的上皮层,优于询问较大的组织体积的常规体积类型光纤拉曼探针;(ⅱ)由于来自更深组织层(例如,基质)的大幅减少的组织自体荧光的贡献,共焦拉曼技术的浅组织询问能力提供了更高的组织拉曼相对自体荧光背景比值,和(iii)结合这种新光纤共焦拉曼光谱平台与良好记载的多变量分析使上皮分子信息能够被提取并实时体内分析。整个共焦显微拉曼内窥镜系统在直观的软件框架中控制,其允许采用到内窥镜检查者的听觉概率反馈来允许在内窥镜筛选设置中的快速调查,推动拉曼光谱的前沿进入常规临床诊断。
总共有450名患者在用于各种迹象的监视或筛选的拉曼内窥镜检查中登记,包括消化不良和上消化道瘤形成。在可疑病变的典型检查期间,每个组织的拉曼光谱测量可在0.5秒内采集,其允许较大组织区域的快速检查。从具有在上消化道中不同组织学亚型的373个患者采集的拉曼光谱数据已被用来构造综合拉曼库(>12000拉曼光谱)。对于招募用于BE筛选和监视的患者,拉曼光谱分为以下三种病理风险的类别:(i)“普通”-柱状上皮层(CLE),(ii)“低风险”BE-定义为杯状细胞的存在,(iii)“高风险”-低度异常增生(LGD)和高度异型增生(HGD)。例如,图40A示出从以不同的组织类型(例如,鳞状层(n=165),CLE(n=907),肠化生(IM)(n=318)和HGD(n=77))呈现的我们数据库中的患者所测量的平均体内共焦拉曼光谱,如由组织病理学特征所证实的。每个拉曼光谱在0.5秒内被采集。为了比较的目的,光谱已经在1445cm-1处被归一化为拉曼峰。显著的组织拉曼峰可以在如下周围观察到:936cm-1(v(C-C)蛋白质)、1004cm-1(苯丙氨酸的Vs(C-C)环呼吸)、1078cm-1(脂质的v(C-C))、1265cm-1(蛋白质的酰胺III v(C-N)和δ(N-H)、1302cm-1(蛋白质的CH2扭曲和摇摆)、1445cm-1(蛋白质和脂类的δ(CH2)变形)、1618cm-1(卟啉的v(C=C))、1655cm-1(蛋白质的酰胺I v(C=O))和1745cm-1(脂质的v(C=O))。显著的拉曼光谱差异(例如,峰值强度、移位和频带变宽)可以在不同的组织类型之间进行辨别。这些丰富的光谱特征描绘伴随Barrett’s癌变的上皮中发生的生物分子和功能变化。虽然组织学检查识别杯状细胞和渐进的架构和细胞学异型的存在(图40(B、C、D、E),但是光纤共焦拉曼光谱揭示,上皮在整个Barrett’s癌变序列中经历主要的功能和生物分子的变化。有趣的是,BE的拉曼生物分子标志高度类似于异型增生,证实到肠化生表型的转化是在Barrett’s癌变中的关键事件。这些高度特异性上皮分子标志可能反映了大量的内源光生物标志物(即,癌蛋白、DNA、粘蛋白表达、有丝分裂等)。因此,上皮拉曼光谱标志与组织病理学或组织化学的相关性可以加深在生物分子水平的原位Barrett’s发生和进展的理解。目前,没有其它竞争光学光谱技术(如荧光、弹性散射光谱)可在内窥镜处提供在体内的这种详尽的分子特征。
预测的(即独立)测量的组织部位的组织病理学表征的152个为CLE(n=597个光谱),48个为IM(n=123个光谱),并且9个为HGD(N=77个光谱)。图41A示出在属于正常、低风险和高风险病变的共焦拉曼光谱的77个患者中的预测测量的风险计分的二维三元散点图。相应的二分法受试者工作特征曲线(ROC)曲线也分别以在0.88、0.84和0.90曲线(AUC)下的区域从图3A中产生的,用于在正常的、低风险和高风险病变之间的辨别。不仅共焦拉曼技术区分具有BE低风险病变(图41A),还能够客观定位包含异常增生上皮的特异性组织区域。上述ROC分析示出,以高风险组织的检测为目标可以实时成功实现,产生在光谱基础上的87.0%(67/77)的诊断灵敏度和84.7%(610/720)的特异性。
在上述的说明中,实施例是所公开的系统和方法的例子或实施方式。“一个实施例”、“实施例”或“一些实施例”的各种出现不一定都指代相同的实施例。
虽然所公开的系统和方法的各种特征可在单个实施例的上下文中进行描述,但是该特征也可以单独地或以任何适当的组合提供。相反地,虽然为清楚起见,所公开的系统和方法可以在本文中在单独实施例的上下文中描述,但是所公开的系统和方法也可以在单个实施例中实现。
此外,还应当理解,所公开的系统和方法可以各种方式执行或实施,并且可以在除了上述说明中概述的之外,在实施例中实现。
在此使用的技术和科学术语含义应通常为如由属于本领域的普通技术人员所理解,除非另有定义。
本发明的某些方面包括以一个方法形式在此描述的过程步骤和指令。应当指出的是,本公开的过程步骤和指令可体现在软件、固件或硬件中,并且当以软件体现时,可以被下载以驻留在并且由实时网络操作系统所使用的不同平台所操作。
本发明还涉及用于执行这里操作的设备。该设备可以被特别地构造用于所需目的,或者它可以包括通用计算机,该通用计算机选择性激活或由存储在能够由计算机访问的计算机可读介质上的计算机程序来重新配置。这种计算机程序可以存储在有形的非临时性计算机可读存储介质中,诸如但不限于任何类型的盘,包括软盘、光盘、CD-ROM、磁光盘、只读内存(ROM)、随机存取内存(RAM)、EPROM、EEPROM、磁卡或光卡、专用集成电路(ASIC),或适于存储电子指令的任何类型介质,并且每一个都耦合到计算机系统总线。此外,在本说明书中提到的计算机可以包括单个处理器或者可以是采用多处理器设计以提高计算能力的体系结构。
本文中所呈现的方法和操作并不固有地涉及任何特定计算机或其它装置。各种通用系统也可根据本文的教导与程序一起使用,或者它可以证明构造更加专业的设备以执行所需方法步骤的便利性。用于多种这些系统的所需结构对于本领域的那些技术人员连同等效变体将是显而易见。另外,本公开不参考任何特定的编程语言来描述。可以理解的是,各种编程语言可以用于实现如本文所述的本公开的教导,并且对于特定语言的任何参考被提供用于启用的公开和本公开的最优模式。
本申请非常适合在众多的拓扑上各种各样的计算机网络系统。在这个领域中,大型网络的配置和管理包括存储设备和计算机,其通过网络,例如因特网、公共网络、私有网络或使计算系统之间的通信可行的其它网络,通信地耦合到不同的计算机和存储设备。最后,应该指出的是,在本说明书中使用的语言原则上选择用于可读性和指导性的目的,并且可能没有被选择以描绘或限制发明的主题。因此,本发明的公开意在是说明性,而不是限制在所附权利要求中阐述的本公开范围。
Claims (33)
1.一种用于癌症的光谱鉴定的计算机实施的方法,其特征在于,包括:
接收与组织相关联的至少一个光谱;
在使用所述光谱的模型中分析所述至少一个光谱以确定得分,其中所述得分指示所述组织是癌性的可能性;以及
输出所述得分。
2.根据权利要求1所述的方法,其特征在于,所述模型使用解释函数来生成,所述解释函数选自由偏最小二乘判别分析、主成分分析线性判别分析、蚁群优化线性判别分析、分类和回归树、支持向量机和自适应增强组成的群组。
3.根据权利要求1所述的方法,其特征在于,所述至少一个光谱由拉曼光谱生成。
4.根据权利要求1所述的方法,其特征在于,所述分析模型中的至少一个光谱包括分析第一模型和第二模型中的至少一个光谱。
5.根据权利要求1所述的方法,其特征在于,所述模型基于组织来选择。
6.根据权利要求1所述的方法,其特征在于,所述得分指示所述组织是否正常、异常增生或瘤形成。
7.根据权利要求1所述的方法,其特征在于,分析所述至少一个光谱包括:
执行异常分析;以及
响应于确定所述至少一个光谱是异常的所述异常分析,抑制所述光谱。
8.根据权利要求7所述的方法,其特征在于,执行异常分析包括主成分分析。
9.根据权利要求7所述的方法,其特征在于,还包括指示音频发射装置发射音频信号,所述音频信号响应于所述异常分析,所述异常分析确定所述至少一个光谱为异常。
10.根据权利要求7所述的方法,其特征在于,还包括指示光谱仪获得至少一个附加光谱,所述至少一个附加光谱响应于所述异常分析,所述异常分析确定所述至少一个光谱为异常并接收所述至少一个附加光谱。
11.根据权利要求6所述的方法,其特征在于,还包括指示音频发射装置发射音频信号,所述音频信号识别所述组织为正常的、异常增生或瘤形成。
12.根据权利要求6所述的方法,其特征在于,与正常的、异常增生或瘤形成的每个相关联的所述音频信号是不同的。
13.根据权利要求12所述的方法,其特征在于,与正常的、异常增生或瘤形成的每个相关联的所述音频信号还与异常光谱相关联的音频信号不同。
14.一种校准光纤拉曼光谱系统的方法,其特征在于,
所述系统包括激光源、光谱仪和光纤探针,所述光纤探针将光从所述激光源发送到目标,并且将散射光返回到所述光谱仪,
所述方法包括将光从所述激光源发送到具有已知光谱的标准目标,
记录来自所述标准目标的散射光的校准光谱,
比较所述已知光谱和校准系统,并且生成传递函数,以及
存储所述传递函数。
15.根据权利要求14所述的方法,其特征在于,包括记录多个光纤探针的每个的校准光谱,计算包括每个所述探针的系统的传递函数,并且将所述传递函数与对应探针相关联。
16.根据权利要求14或15所述的方法,其特征在于,包括随后照射测试对象、记录光谱并且根据存储的传递函数校准所述光谱的步骤。
17.根据权利要求14~16中任一项所述的方法,其特征在于,所述光谱仪具有相关联的光谱仪传递函数,所述探针具有相关联的探针传递函数,所述传递函数是光谱仪传递函数和探针传递函数中的一种函数。
18.根据权利要求14~17中任一项所述的方法,其特征在于,包括:在主光谱仪系统上计算主光纤探针的第一传递函数、次级光纤探针的第二传递函数、以及计算基于所述第一传递函数和第二传递函数的校准函数。
19.根据权利要求18所述的方法,其特征在于,包括将所述校准函数与所述次级光纤探针相关联。
20.根据权利要求18或19所述的方法,其特征在于,包括:在次级光谱仪系统上使用所述主光纤探针、、生成次级系统传递函数、以及存储所述次级系统传递函数。
21.根据权利要求20所述的方法,其特征在于,包括使用具有所述次级光谱仪系统的次级光纤探针,以及根据所述校准函数修正存储的次级系统传递函数。
22.根据权利要求18~21中任一项所述的方法,其特征在于,包括根据所述主光谱仪系统执行所述次级光谱仪系统的波长轴校准的初始步骤。
23.一种操作拉曼光谱系统的方法,其特征在于,
所述系统包括激光源、光谱仪和光纤探针,所述光纤探针将光从所述激光源发送到目标,并且将散射光返回到所述光谱仪,
所述方法包括从所述激光源发送光到具有已知光谱的目标,
记录来自所述目标的散射光的光谱,以及
根据存储的传递函数修正记录的光谱。
24.根据权利要求23所述的方法,其特征在于,所述存储的传递函数与所述光谱仪和所述光纤探针对应。
25.根据权利要求23所述的方法,其特征在于,所述存储的传递函数与所述光谱仪和主光纤探针对应,所述方法还包括根据与所述光纤探针相关联的存储校准函数来修正存储的传递函数。
26.一种拉曼光谱系统,其特征在于,包括激光源、光谱仪、光纤探针和存储的传递函数,所述光纤探针将光从所述激光源发送到目标并且将散射光返回到所述光谱仪,
所述系统可操作将光从所述激光源发送到具有已知光谱的目标,
记录来自所述目标的散射光的光谱,以及
根据所述存储的传递函数修正记录的光谱。
27.根据权利要求26所述的系统,其特征在于,所述存储的传递函数与所述光谱仪和所述光纤探针对应。
28.根据权利要求26所述的系统,其特征在于,所述存储的传递函数与所述光谱仪和主光纤探针对应,并且所述方法还包括根据与所述光纤探针相关联的存储的校准函数来修正存储的传递函数。
29.一种预测在光纤拉曼光谱系统中发送的激光功率的方法,其特征在于,
所述系统包括激光源、光谱仪和光纤探针,所述光纤探针将光从所述激光源发送到目标,并且将散射光返回到所述光谱仪,
所述方法包括从所述激光源发送光到多个目标,
对于每个测试目标,测量来自所述激光源的光的发送功率、和所述光谱仪处散射光的光谱,
执行所述捕捉光谱的多变量分析,测量的发送功率作为因变量,
以及存储所得的模型。
30.根据权利要求29所述的方法,其特征在于,包括如下步骤,
发送激光到测试目标,
向所述模型供应捕捉光谱,以及
计算所述发送功率的估量。
31.一种从具有光纤探针的光纤拉曼光谱系统减去背景信号的方法,其特征在于,所述方法包括如下步骤:
a)存储背景光谱,
b)接收测试光谱,
c)使用一个或多个参考峰来估计背景影响,
d)基于估计的背景影响,将所述背景光谱与校准因子相乘,并将其从所述测试光谱减去,
e)检查所述测试光谱余留的背景影响,以及
f)如果所述背景贡献可忽略不计,则输出所述测试光谱,否则重复步骤(c)到(e)。
32.根据权利要求31所述的方法,其特征在于,所述一个或多个参考峰包括与所述光纤探针中二氧化硅或蓝宝石对应的一个或多个峰。
33.根据紧随权利要求1至13中任一项方法的权利要求14至18中任一项所述的校准方法。
Applications Claiming Priority (5)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US201261667384P | 2012-07-02 | 2012-07-02 | |
US61/667,384 | 2012-07-02 | ||
GB1307338.2A GB2513343A (en) | 2013-04-23 | 2013-04-23 | Methods related to instrument-independent measurements for quantitative analysis of fiber-optic Raman spectroscopy |
GB1307338.2 | 2013-04-23 | ||
PCT/SG2013/000273 WO2014007759A1 (en) | 2012-07-02 | 2013-07-02 | Methods related to real-time cancer diagnostics at endoscopy utilizing fiber-optic raman spectroscopy |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104541153A true CN104541153A (zh) | 2015-04-22 |
CN104541153B CN104541153B (zh) | 2017-12-05 |
Family
ID=48537691
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201380035488.5A Active CN104541153B (zh) | 2012-07-02 | 2013-07-02 | 涉及利用光纤拉曼光谱的内窥镜进行实时癌症诊断的方法 |
Country Status (10)
Country | Link |
---|---|
US (1) | US9931039B2 (zh) |
EP (1) | EP2870462B1 (zh) |
JP (1) | JP6230077B2 (zh) |
KR (1) | KR101947859B1 (zh) |
CN (1) | CN104541153B (zh) |
CA (1) | CA2917308C (zh) |
ES (1) | ES2930428T3 (zh) |
GB (1) | GB2513343A (zh) |
SG (1) | SG11201408793XA (zh) |
WO (1) | WO2014007759A1 (zh) |
Cited By (24)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105534470A (zh) * | 2015-12-22 | 2016-05-04 | 精微视达医疗科技(武汉)有限公司 | 一种共焦显微内窥镜系统及其调节方法 |
CN105997000A (zh) * | 2016-05-30 | 2016-10-12 | 福建师范大学 | 一种基于纤维内窥镜的拉曼光谱检测装置及其实现方法 |
CN106092996A (zh) * | 2016-08-03 | 2016-11-09 | 复旦大学 | 一种基于自体荧光寿命的癌症诊断系统 |
CN106420041A (zh) * | 2015-08-13 | 2017-02-22 | 复旦大学附属肿瘤医院 | 一种带有光纤探针的高频手术电刀 |
CN106691549A (zh) * | 2015-08-13 | 2017-05-24 | 复旦大学附属肿瘤医院 | 一种带有光纤探针的穿刺活检针 |
CN107664532A (zh) * | 2016-07-29 | 2018-02-06 | 唯亚威解决方案股份有限公司 | 使用稀疏转移集合转移校准模型 |
CN108007910A (zh) * | 2017-12-28 | 2018-05-08 | 中国科学院合肥物质科学研究院 | 一种用于皮肤组织荧光光谱和后向散射光光谱原位测量的深度敏感光纤探针 |
CN108169200A (zh) * | 2016-12-07 | 2018-06-15 | 同方威视技术股份有限公司 | 用于校准拉曼光谱检测装置的方法 |
CN108287129A (zh) * | 2018-03-22 | 2018-07-17 | 中国计量大学 | 多通道荧光谱生物气溶胶粒子的检测装置 |
CN108692813A (zh) * | 2017-04-05 | 2018-10-23 | 凯塞光学系统股份有限公司 | 用于近红外检测器的辐射护罩 |
CN109068969A (zh) * | 2016-03-10 | 2018-12-21 | 比奥普-医疗有限公司 | 诊断组织的装置 |
CN109073536A (zh) * | 2016-05-06 | 2018-12-21 | 索尼公司 | 信息处理装置、信息处理方法、程序以及信息处理系统 |
CN109259718A (zh) * | 2018-11-06 | 2019-01-25 | 王瑞贞 | 一种可拆分组合的内窥镜检测系统 |
CN109730694A (zh) * | 2019-01-08 | 2019-05-10 | 北京信息科技大学 | 一种诊断人体血管内斑块类型的拉曼光谱测量探针 |
CN109781701A (zh) * | 2019-01-18 | 2019-05-21 | 拉曼兄弟(深圳)科技发展有限公司 | 一种基于拉曼光谱技术的甲状旁腺术中实时检测方法 |
CN110337258A (zh) * | 2017-02-28 | 2019-10-15 | 威里利生命科学有限责任公司 | 使用可编程光源对图像进行多类分类的系统和方法 |
CN110957034A (zh) * | 2018-09-26 | 2020-04-03 | 金敏 | 疾病预测系统 |
CN111281313A (zh) * | 2020-03-05 | 2020-06-16 | 中国科学院长春光学精密机械与物理研究所 | 基于拉曼光谱的病灶在线检测系统 |
CN111770720A (zh) * | 2018-01-17 | 2020-10-13 | 奥迪思医疗公司 | 用于实时拉曼光谱法以进行癌症检测的系统和方法 |
CN112089404A (zh) * | 2020-10-10 | 2020-12-18 | 南京航空航天大学 | 一种内窥OCT-Raman双模成像装置及成像方法 |
CN113624734A (zh) * | 2021-06-25 | 2021-11-09 | 港湾之星健康生物(深圳)有限公司 | 兼具灵敏度和响应速度的拉曼光谱检测方法与系统 |
CN114449968A (zh) * | 2019-08-05 | 2022-05-06 | 捷锐士阿希迈公司(以奥林巴斯美国外科技术名义) | 激光光纤至目标的距离控制 |
CN114615946A (zh) * | 2019-08-05 | 2022-06-10 | 捷锐士阿希迈公司(以奥林巴斯美国外科技术名义) | 激光治疗的信号协调传送 |
CN117152362A (zh) * | 2023-10-27 | 2023-12-01 | 深圳市中安视达科技有限公司 | 内窥镜多光谱的多路成像方法、装置、设备及存储介质 |
Families Citing this family (51)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP6560485B2 (ja) * | 2014-04-17 | 2019-08-14 | ローム株式会社 | 診断システム |
WO2014175223A1 (ja) | 2013-04-22 | 2014-10-30 | ローム株式会社 | 癌診断装置、診断システム、診断装置 |
GB201315195D0 (en) | 2013-08-23 | 2013-10-09 | Perkinelmer Uk Ltd | Identifying presence of substances |
DE102014003386B4 (de) | 2014-03-07 | 2016-06-23 | Celltool Gmbh | Verfahren und Vorrichtung zur Qualitätskontrolle eines Blutprodukts |
EP3128917A4 (en) * | 2014-04-08 | 2018-01-03 | Polyvalor, Limited Partnership | System and method for assessing a cancer status of biological tissue |
AU2015284810B2 (en) * | 2014-07-02 | 2020-03-12 | National University Of Singapore | Raman spectroscopy system, apparatus, and method for analyzing, characterizing, and/or diagnosing a type or nature of a sample or a tissue such as an abnormal growth |
CN104116497B (zh) * | 2014-07-22 | 2016-10-05 | 中国科学院自动化研究所 | 内窥式光学分子影像导航系统和多光谱成像方法 |
GB2547350B (en) * | 2014-08-29 | 2022-05-04 | Synaptive Medical Inc | Molecular cell imaging using optical spectroscopy |
EP3258834A4 (en) * | 2015-02-19 | 2018-09-26 | Ramot at Tel-Aviv University Ltd. | Chest radiograph (cxr) image analysis |
WO2017024145A1 (en) | 2015-08-05 | 2017-02-09 | Canon U.S.A., Inc. | Forward and angle view endoscope |
CN117434627A (zh) * | 2015-08-26 | 2024-01-23 | 新加坡国立大学 | 用于保持微球体的膜 |
EP3377865A4 (en) * | 2015-11-16 | 2019-06-12 | ChemImage Corporation | SYSTEMS AND METHODS FOR RAMAN-BASED IMMUNOASSAYS |
WO2017117203A1 (en) | 2015-12-28 | 2017-07-06 | Canon U.S.A., Inc. | Optical probe, light intensity detection, imaging method and system |
US10365229B2 (en) * | 2015-12-31 | 2019-07-30 | Kaiser Optical Systems, Inc. | Real-time characterization of pharmaceutical tablet coatings using Raman spectroscopy |
US10551245B2 (en) | 2016-02-12 | 2020-02-04 | Canon U.S.A., Inc. | Simple monolithic optical element for forward-viewing spectrally encoded endoscopy |
KR20170114229A (ko) * | 2016-04-04 | 2017-10-13 | 스페클립스 주식회사 | 질병 진단과 피부 나이 측정이 가능한 시스템과 이에 사용되는 핸드피스 |
WO2017212307A1 (en) | 2016-06-06 | 2017-12-14 | Hajdu Imre | Diagnostic medical device |
US10321810B2 (en) | 2016-06-13 | 2019-06-18 | Canon U.S.A., Inc. | Spectrally encoded endoscopic probe having a fixed fiber |
US10234694B2 (en) | 2016-07-15 | 2019-03-19 | Canon U.S.A., Inc. | Spectrally encoded probes |
US10984519B2 (en) | 2016-08-03 | 2021-04-20 | Panduit Corp. | Method for evaluating the optical insertion loss of a mechanical splice joint of two optical fibers |
US10646111B2 (en) | 2016-09-23 | 2020-05-12 | Canon U.S.A., Inc. | Spectrally encoded endoscopy apparatus and methods |
KR20180042700A (ko) | 2016-10-18 | 2018-04-26 | 삼성전자주식회사 | 스펙트럼 안정성 모니터링 장치 및 방법 |
JP2018094395A (ja) | 2016-11-03 | 2018-06-21 | キヤノン ユーエスエイ, インコーポレイテッドCanon U.S.A., Inc | 診断用スペクトル符号化内視鏡検査装置およびシステム、ならびにこれらと共に使用するための方法 |
US10371642B2 (en) * | 2016-12-08 | 2019-08-06 | The Board Of Trustees Of The Leland Stanford Junior University | Raman topography system and methods of imaging |
CN106923781A (zh) * | 2017-03-28 | 2017-07-07 | 戎创前沿科技(北京)有限公司 | 一种拉曼胃镜诊断仪 |
FR3065285B1 (fr) * | 2017-04-14 | 2019-07-05 | Commissariat A L'energie Atomique Et Aux Energies Alternatives | Procede d’acquisition d’un spectre |
US10801963B2 (en) | 2018-08-22 | 2020-10-13 | Paul Bartholomew | Raman spectroscopy for minerals identification |
US11583177B2 (en) | 2017-08-06 | 2023-02-21 | Biop—Medical Ltd. | Optical probe for cervical examination |
CN107320062B (zh) * | 2017-08-11 | 2023-07-25 | 山东大学齐鲁医院 | 一种带有活体实时拉曼光谱检测的膀胱镜及方法 |
CA3066429A1 (en) * | 2017-08-22 | 2019-02-28 | Leica Biosystems Melbourne Pty Ltd | A system and method of monitoring tissue samples to be processed by a tissue processor |
US10825152B2 (en) | 2017-09-14 | 2020-11-03 | Canon U.S.A., Inc. | Distortion measurement and correction for spectrally encoded endoscopy |
CN108155936B (zh) * | 2018-01-03 | 2023-05-26 | 成都航天通信设备有限责任公司 | 一种多通道综合数据链光纤测试方法 |
KR102167946B1 (ko) * | 2018-02-15 | 2020-10-20 | 스페클립스 주식회사 | 인-비보 병변 조직 검출을 위한 레이저 분광 기반의 독립 장치 및 방법 |
KR102258181B1 (ko) * | 2018-02-15 | 2021-05-28 | 스페클립스 주식회사 | 레이저 분광 기반의 인-비보 병변 조직 검출을 위한 넌-디스크리트 스펙트럼 분석 알고리즘 및 이를 위한 방법 |
US11672425B2 (en) | 2018-02-15 | 2023-06-13 | Speclipse, Inc. | Stand-alone apparatus and methods for in vivo detection of tissue malignancy using laser spectroscopy |
US11754567B2 (en) * | 2018-04-30 | 2023-09-12 | City Of Hope | Cancer detection and ablation system and method |
US10928309B2 (en) * | 2018-06-29 | 2021-02-23 | Viavi Solutions Inc. | Cross-validation based calibration of a spectroscopic model |
WO2020077445A1 (en) * | 2018-10-16 | 2020-04-23 | Polyvalor, Limited Partnership | Methods for performing a raman spectroscopy measurement on a sample and raman spectroscopy systems |
RS63454B1 (sr) * | 2018-12-17 | 2022-08-31 | Evonik Operations Gmbh | Postupak za identifikaciju neispravno kalibrisanog ili nekalibrisanog infracrvenog spektrometra i za korekciju spektra snimljenog na neispravno kalibrisanom ili nekalibrisanom infracrvenom spektrometru |
US11726034B2 (en) * | 2019-03-07 | 2023-08-15 | Missouri State University | IR spectra matching methods |
AU2019437773B2 (en) | 2019-03-22 | 2021-10-28 | Speclipse, Inc. | Diagnostic method using laser induced breakdown spectroscopy, and diagnostic device for performing same |
US20220238187A1 (en) * | 2019-06-03 | 2022-07-28 | Frederic Leblond | Methods and Systems for Assessing a Phenotype of a Biological Tissue of a Patient Using Raman Spectroscopy |
CA3148141A1 (en) * | 2019-08-16 | 2021-02-25 | Paul Bartholomew | Raman spectroscopy for minerals identification |
WO2021067677A1 (en) * | 2019-10-02 | 2021-04-08 | Chemimage Corporation | Fusion of molecular chemical imaging with rgb imaging |
EP4078508A4 (en) * | 2019-12-18 | 2023-11-22 | ChemImage Corporation | SYSTEMS AND METHODS FOR COMBINING IMAGING MODALITIES FOR ENHANCED TISSUE DETECTION |
US20240023791A1 (en) * | 2020-08-26 | 2024-01-25 | Baylor College Of Medicine | End expandable optical fiber bundle |
KR102514246B1 (ko) * | 2021-03-26 | 2023-03-27 | 전주대학교 산학협력단 | 머신러닝에 기반한 x-선 상 척추 회전각도 측정시스템 및 이를 이용한 측정방법 |
WO2023052608A1 (en) * | 2021-10-01 | 2023-04-06 | Trinamix Gmbh | Multi single detector application specific spectrometer |
GB202203993D0 (en) * | 2022-03-22 | 2022-05-04 | Keit Ltd | Methods of analysing a sample and calibrating a spectrometer |
EP4309568A1 (en) * | 2022-07-19 | 2024-01-24 | Canostix GmbH | Method and apparatus for spectroscopic cancer screening |
CN116849612B (zh) * | 2023-07-07 | 2024-01-09 | 广东省新黄埔中医药联合创新研究院 | 一种多光谱舌象图像采集分析系统 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5976885A (en) * | 1995-11-13 | 1999-11-02 | Bio-Rad Laboratories, Inc. | Method for the detection of cellular abnormalities using infrared spectroscopic imaging |
CN1264466A (zh) * | 1997-03-14 | 2000-08-23 | 伊斯特曼化学公司 | 使喇曼光谱仪标准化,以获得稳定、可转移的校准的方法 |
CN1105296C (zh) * | 1996-08-22 | 2003-04-09 | 伊斯曼化学公司 | 定量监测化学组合物的组份的方法 |
US6621574B1 (en) * | 2000-05-25 | 2003-09-16 | Inphotonics, Inc. | Dual function safety and calibration accessory for raman and other spectroscopic sampling |
US20030191398A1 (en) * | 2002-04-05 | 2003-10-09 | Massachusetts Institute Of Technology | Systems and methods for spectroscopy of biological tissue |
CN102072894A (zh) * | 2009-11-25 | 2011-05-25 | 欧普图斯(苏州)光学纳米科技有限公司 | 用基于纳米结构的光谱检测方法检测化学物和生化物杂质 |
WO2011072380A1 (en) * | 2009-12-17 | 2011-06-23 | British Columbia Cancer Agency Branch | Apparatus and methods for in vivo tissue characterization by raman spectroscopy |
Family Cites Families (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5357336A (en) * | 1991-08-08 | 1994-10-18 | The Dow Chemical Company | Method and apparatus for multivariate characterization of optical instrument response |
US5638172A (en) * | 1994-05-27 | 1997-06-10 | Eastman Chemical Company | On-line quantitative analysis of chemical compositions by raman spectrometry |
JPH09145619A (ja) * | 1995-09-20 | 1997-06-06 | Kdk Corp | 散乱光等の分光測定方法及び装置 |
US6351306B1 (en) * | 1999-10-19 | 2002-02-26 | Kaiser Optical Systems, Inc. | Optical measurement probe calibration configurations |
JP4588324B2 (ja) * | 2002-04-05 | 2010-12-01 | マサチユセツツ・インスチチユート・オブ・テクノロジイ | 組織測定用プローブ |
WO2006064446A1 (en) * | 2004-12-15 | 2006-06-22 | Koninklijke Philips Electronics N.V. | Calibration of optical analysis making use of multivariate optical elements |
US7583378B2 (en) * | 2006-02-23 | 2009-09-01 | Itt Manufacturing Enterprises, Inc. | Spectrograph calibration using known light source and Raman scattering |
CN102297856B (zh) * | 2010-06-25 | 2015-09-30 | 清华大学 | 用于自动校准拉曼光谱检测系统的方法及拉曼光谱检测系统 |
DE102011077564B4 (de) | 2011-06-15 | 2016-08-25 | Sirona Dental Systems Gmbh | Verfahren zur optischen dreidimensionalen Vermessung eines dentalen Objekts |
-
2013
- 2013-04-23 GB GB1307338.2A patent/GB2513343A/en not_active Withdrawn
- 2013-07-02 US US14/412,053 patent/US9931039B2/en active Active
- 2013-07-02 CN CN201380035488.5A patent/CN104541153B/zh active Active
- 2013-07-02 CA CA2917308A patent/CA2917308C/en active Active
- 2013-07-02 EP EP13813251.9A patent/EP2870462B1/en active Active
- 2013-07-02 WO PCT/SG2013/000273 patent/WO2014007759A1/en active Application Filing
- 2013-07-02 JP JP2015520131A patent/JP6230077B2/ja active Active
- 2013-07-02 SG SG11201408793XA patent/SG11201408793XA/en unknown
- 2013-07-02 ES ES13813251T patent/ES2930428T3/es active Active
- 2013-07-02 KR KR1020157002866A patent/KR101947859B1/ko active IP Right Grant
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5976885A (en) * | 1995-11-13 | 1999-11-02 | Bio-Rad Laboratories, Inc. | Method for the detection of cellular abnormalities using infrared spectroscopic imaging |
CN1105296C (zh) * | 1996-08-22 | 2003-04-09 | 伊斯曼化学公司 | 定量监测化学组合物的组份的方法 |
CN1264466A (zh) * | 1997-03-14 | 2000-08-23 | 伊斯特曼化学公司 | 使喇曼光谱仪标准化,以获得稳定、可转移的校准的方法 |
US6621574B1 (en) * | 2000-05-25 | 2003-09-16 | Inphotonics, Inc. | Dual function safety and calibration accessory for raman and other spectroscopic sampling |
US20030191398A1 (en) * | 2002-04-05 | 2003-10-09 | Massachusetts Institute Of Technology | Systems and methods for spectroscopy of biological tissue |
CN102072894A (zh) * | 2009-11-25 | 2011-05-25 | 欧普图斯(苏州)光学纳米科技有限公司 | 用基于纳米结构的光谱检测方法检测化学物和生化物杂质 |
WO2011072380A1 (en) * | 2009-12-17 | 2011-06-23 | British Columbia Cancer Agency Branch | Apparatus and methods for in vivo tissue characterization by raman spectroscopy |
Cited By (34)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106420041A (zh) * | 2015-08-13 | 2017-02-22 | 复旦大学附属肿瘤医院 | 一种带有光纤探针的高频手术电刀 |
CN106691549A (zh) * | 2015-08-13 | 2017-05-24 | 复旦大学附属肿瘤医院 | 一种带有光纤探针的穿刺活检针 |
CN105534470A (zh) * | 2015-12-22 | 2016-05-04 | 精微视达医疗科技(武汉)有限公司 | 一种共焦显微内窥镜系统及其调节方法 |
CN109068969A (zh) * | 2016-03-10 | 2018-12-21 | 比奥普-医疗有限公司 | 诊断组织的装置 |
CN109073536A (zh) * | 2016-05-06 | 2018-12-21 | 索尼公司 | 信息处理装置、信息处理方法、程序以及信息处理系统 |
CN105997000A (zh) * | 2016-05-30 | 2016-10-12 | 福建师范大学 | 一种基于纤维内窥镜的拉曼光谱检测装置及其实现方法 |
CN107664532A (zh) * | 2016-07-29 | 2018-02-06 | 唯亚威解决方案股份有限公司 | 使用稀疏转移集合转移校准模型 |
US11378452B2 (en) | 2016-07-29 | 2022-07-05 | Viavi Solutions Inc. | Transfer of a calibration model using a sparse transfer set |
CN106092996A (zh) * | 2016-08-03 | 2016-11-09 | 复旦大学 | 一种基于自体荧光寿命的癌症诊断系统 |
CN106092996B (zh) * | 2016-08-03 | 2024-03-15 | 复旦大学 | 一种基于自体荧光寿命的癌症诊断系统 |
CN108169200A (zh) * | 2016-12-07 | 2018-06-15 | 同方威视技术股份有限公司 | 用于校准拉曼光谱检测装置的方法 |
CN108169200B (zh) * | 2016-12-07 | 2020-12-22 | 同方威视技术股份有限公司 | 用于校准拉曼光谱检测装置的方法 |
US11699102B2 (en) | 2017-02-28 | 2023-07-11 | Verily Life Sciences Llc | System and method for multiclass classification of images using a programmable light source |
CN110337258A (zh) * | 2017-02-28 | 2019-10-15 | 威里利生命科学有限责任公司 | 使用可编程光源对图像进行多类分类的系统和方法 |
US10753805B2 (en) | 2017-04-05 | 2020-08-25 | Kaiser Optical Systems Inc. | Radiation shield for near-infrared detectors |
CN108692813A (zh) * | 2017-04-05 | 2018-10-23 | 凯塞光学系统股份有限公司 | 用于近红外检测器的辐射护罩 |
CN108007910A (zh) * | 2017-12-28 | 2018-05-08 | 中国科学院合肥物质科学研究院 | 一种用于皮肤组织荧光光谱和后向散射光光谱原位测量的深度敏感光纤探针 |
CN108007910B (zh) * | 2017-12-28 | 2023-12-29 | 中国科学院合肥物质科学研究院 | 一种用于皮肤组织荧光光谱和后向散射光光谱原位测量的深度敏感光纤探针 |
CN111770720A (zh) * | 2018-01-17 | 2020-10-13 | 奥迪思医疗公司 | 用于实时拉曼光谱法以进行癌症检测的系统和方法 |
CN108287129A (zh) * | 2018-03-22 | 2018-07-17 | 中国计量大学 | 多通道荧光谱生物气溶胶粒子的检测装置 |
CN110957034A (zh) * | 2018-09-26 | 2020-04-03 | 金敏 | 疾病预测系统 |
CN110957033A (zh) * | 2018-09-26 | 2020-04-03 | 金敏 | 疾病预测系统 |
CN109259718A (zh) * | 2018-11-06 | 2019-01-25 | 王瑞贞 | 一种可拆分组合的内窥镜检测系统 |
CN109730694A (zh) * | 2019-01-08 | 2019-05-10 | 北京信息科技大学 | 一种诊断人体血管内斑块类型的拉曼光谱测量探针 |
CN109781701A (zh) * | 2019-01-18 | 2019-05-21 | 拉曼兄弟(深圳)科技发展有限公司 | 一种基于拉曼光谱技术的甲状旁腺术中实时检测方法 |
CN114615946A (zh) * | 2019-08-05 | 2022-06-10 | 捷锐士阿希迈公司(以奥林巴斯美国外科技术名义) | 激光治疗的信号协调传送 |
CN114449968A (zh) * | 2019-08-05 | 2022-05-06 | 捷锐士阿希迈公司(以奥林巴斯美国外科技术名义) | 激光光纤至目标的距离控制 |
US11957410B2 (en) | 2019-08-05 | 2024-04-16 | Gyrus Acmi, Inc. | Laser fiber-to-target distance control |
CN114615946B (zh) * | 2019-08-05 | 2024-05-17 | 捷锐士阿希迈公司(以奥林巴斯美国外科技术名义) | 激光治疗的信号协调传送 |
CN111281313A (zh) * | 2020-03-05 | 2020-06-16 | 中国科学院长春光学精密机械与物理研究所 | 基于拉曼光谱的病灶在线检测系统 |
CN112089404A (zh) * | 2020-10-10 | 2020-12-18 | 南京航空航天大学 | 一种内窥OCT-Raman双模成像装置及成像方法 |
CN113624734A (zh) * | 2021-06-25 | 2021-11-09 | 港湾之星健康生物(深圳)有限公司 | 兼具灵敏度和响应速度的拉曼光谱检测方法与系统 |
CN117152362A (zh) * | 2023-10-27 | 2023-12-01 | 深圳市中安视达科技有限公司 | 内窥镜多光谱的多路成像方法、装置、设备及存储介质 |
CN117152362B (zh) * | 2023-10-27 | 2024-05-28 | 深圳市中安视达科技有限公司 | 内窥镜多光谱的多路成像方法、装置、设备及存储介质 |
Also Published As
Publication number | Publication date |
---|---|
KR101947859B1 (ko) | 2019-02-13 |
EP2870462A1 (en) | 2015-05-13 |
EP2870462B1 (en) | 2022-08-17 |
JP6230077B2 (ja) | 2017-11-15 |
SG11201408793XA (en) | 2015-01-29 |
CN104541153B (zh) | 2017-12-05 |
US20150335248A1 (en) | 2015-11-26 |
US9931039B2 (en) | 2018-04-03 |
KR20150036345A (ko) | 2015-04-07 |
WO2014007759A1 (en) | 2014-01-09 |
GB2513343A (en) | 2014-10-29 |
CA2917308A1 (en) | 2014-01-09 |
GB201307338D0 (en) | 2013-05-29 |
ES2930428T3 (es) | 2022-12-12 |
EP2870462A4 (en) | 2016-03-23 |
JP2015526135A (ja) | 2015-09-10 |
CA2917308C (en) | 2021-07-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104541153A (zh) | 涉及利用光纤拉曼光谱的内窥镜进行实时癌症诊断的方法 | |
US11419499B2 (en) | Optical coherence tomography for cancer screening and triage | |
Fabelo et al. | Deep learning-based framework for in vivo identification of glioblastoma tumor using hyperspectral images of human brain | |
Terry et al. | Detection of dysplasia in Barrett's esophagus with in vivo depth-resolved nuclear morphology measurements | |
US8649849B2 (en) | Optical methods to intraoperatively detect positive prostate and kidney cancer margins | |
Heng et al. | Advances in real‐time fiber‐optic Raman spectroscopy for early cancer diagnosis: Pushing the frontier into clinical endoscopic applications | |
Barman et al. | Selective sampling using confocal Raman spectroscopy provides enhanced specificity for urinary bladder cancer diagnosis | |
US20170138860A1 (en) | Raman spectroscopy system, apparatus, and method for analyzing, characterizing, and/or diagnosing a type or nature of a sample or a tissue such as an abnormal growth | |
US20120022367A1 (en) | Chemically-selective, label free, microendoscopic system based on coherent anti-stokes raman scattering and microelectromechanical fiber optic probe | |
Yan et al. | Tongue squamous cell carcinoma discrimination with Raman spectroscopy and convolutional neural networks | |
CN105190262A (zh) | 诊断仪器和方法 | |
CN102740762A (zh) | 由拉曼光谱法表征肺组织的设备和方法 | |
CN104603601A (zh) | 与拉曼光谱学有关的诊断仪器和方法 | |
US11445915B2 (en) | Compact briefcase OCT system for point-of-care imaging | |
Heidari et al. | The use of optical coherence tomography and convolutional neural networks to distinguish normal and abnormal oral mucosa | |
US20150062320A1 (en) | Diffuse reflectance hyperspectral imaging system | |
Chang et al. | Deep learning methods for oral cancer detection using Raman spectroscopy | |
Mangotra et al. | Hyperspectral imaging for early diagnosis of diseases: A review | |
Yu et al. | Deep multi-feature fusion residual network for oral squamous cell carcinoma classification and its intelligent system using Raman spectroscopy | |
Bailey et al. | Detection of precancerous lesions in the oral cavity using oblique polarized reflectance spectroscopy: a clinical feasibility study | |
Valor et al. | Study of methionine choline deficient diet-induced steatosis in mice using endogenous fluorescence spectroscopy | |
Goodacre et al. | Clinical Spectroscopy: general discussion | |
Geldof et al. | Toward the use of diffuse reflection spectroscopy for intra-operative tissue discrimination during sarcoma surgery | |
van Doorn et al. | Raman spectroscopy for guidance of vulvar cancer surgery: a pilot study | |
Cordero | Methods and instrumentation for raman characterization of bladder cancer tumor |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
EE01 | Entry into force of recordation of patent licensing contract | ||
EE01 | Entry into force of recordation of patent licensing contract |
Application publication date: 20150422 Assignee: Ente Fu Pte Ltd Assignor: NATIONAL University OF SINGAPORE Contract record no.: X2020990000498 Denomination of invention: A method for real-time cancer diagnosis using an endoscope using fiber Raman spectroscopy is described Granted publication date: 20171205 License type: Exclusive License Record date: 20200918 |