CN101816552B - Multi-spectral time-domain fluorescence molecular tomography measurement system imitating CT scanning mode - Google Patents

Multi-spectral time-domain fluorescence molecular tomography measurement system imitating CT scanning mode Download PDF

Info

Publication number
CN101816552B
CN101816552B CN 201010167623 CN201010167623A CN101816552B CN 101816552 B CN101816552 B CN 101816552B CN 201010167623 CN201010167623 CN 201010167623 CN 201010167623 A CN201010167623 A CN 201010167623A CN 101816552 B CN101816552 B CN 101816552B
Authority
CN
China
Prior art keywords
wavelength
optical fiber
fluorescence
laser
scanning mode
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Expired - Fee Related
Application number
CN 201010167623
Other languages
Chinese (zh)
Other versions
CN101816552A (en
Inventor
高峰
李娇
赵会娟
张丽敏
周仲兴
张伟
易茜
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Tianjin University
Original Assignee
Tianjin University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Tianjin University filed Critical Tianjin University
Priority to CN 201010167623 priority Critical patent/CN101816552B/en
Publication of CN101816552A publication Critical patent/CN101816552A/en
Application granted granted Critical
Publication of CN101816552B publication Critical patent/CN101816552B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/0059Measuring for diagnostic purposes; Identification of persons using light, e.g. diagnosis by transillumination, diascopy, fluorescence
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B2503/00Evaluating a particular growth phase or type of persons or animals
    • A61B2503/40Animals

Landscapes

  • Life Sciences & Earth Sciences (AREA)
  • Health & Medical Sciences (AREA)
  • Medical Informatics (AREA)
  • Biophysics (AREA)
  • Pathology (AREA)
  • Engineering & Computer Science (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Physics & Mathematics (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • General Health & Medical Sciences (AREA)
  • Public Health (AREA)
  • Veterinary Medicine (AREA)
  • Investigating, Analyzing Materials By Fluorescence Or Luminescence (AREA)

Abstract

The invention relates to the field of small animal molecular imaging, and provides an imitating CT scanning mode multispectral time domain fluorescence molecular tomography measuring system with ultrasensitive and high time-space resolution measurement performances so as to realize the high-quality multi-component multi-parameter imaging function. The imitating CT scanning mode multispectral time domain fluorescence molecular tomography measuring system comprises two picosecond semiconductor lasers, an optical fiber coupler, an incident optical fiber, an imaging cavity, a rotating platform, a receiving optical fiber, a motor driving filter wheel and a detection module, wherein the two picosecond semiconductor lasers and the optical fiber coupler are used for providing wavelength ultrashort laser required, the optical fiber coupler is used for coupling two same-frequency ultrashort laser pulses of different wavelengths into a bundle of laser according to fixed time interval, the incident optical fiber is used for projecting emergent laser of two or one of both picosecond semiconductor lasers onto an object, the imaging cavity and the rotating platform are used for realizing the imitating CT scanning mode, the receiving optical fiber is used for receiving reflective or transmitted laser from the object, and the detection module is used for receiving emergent light passing through corresponding filter pieces. The invention is mainly applied to the molecular imaging of small animals.

Description

The imitating CT scanning mode multispectral time domain fluorescence molecular tomography measuring system
Technical field
The present invention relates to the small animal molecule image-forming field, especially, be specifically related to the imitating CT scanning mode multispectral time domain fluorescence molecular tomography measuring system towards the time-domain fluorescent diffusion chromatography system and the time resolution optical chromatography system of small animal molecule image-forming.
Background technology
The molecular imaging qualitative and quantitative research that to be image application method carry out at cell and molecular level the bioprocess under the condition of living organism, be with the generation mechanism of body internal specific molecule, reflect intravital physiology and pathological change process at cell or molecular level by image as image contrast.The advantage of molecular imaging can mainly be summarised as 3 points: one, the molecular image technology can become intuitive image with the process of complexity such as gene expression, bio signal transmission, makes people can understand the mechanism and the feature of disease better on the molecular cell level; Its two, can find variation of the early stage molecular cell of disease and pathological change process; Its three, mechanism and the effect that can on live body, observe medicine or gene therapy continuously.The molecular image technology has continuous, quick, remote, undamaged advantage as a kind of detection method at live body (in vivo), and the 3-D view of molecule in the body or cell more may be provided.It can disclose the early molecule biological property of pathological changes, has promoted the early diagnosis and therapy of disease, has also introduced new notion for clinical diagnosis.
Optical molecular imaging is important molecular image mode, has multiple advantages such as ultra-sensitivity, high specific, real-time and labelling motility concurrently, is particularly suitable for the biomedical basic research based on the toy pathological model.。Especially (Near-Infrared Spectroscopy, NIRS) development is very rapid on biological tissue's imaging in vivo is used for near-infrared spectrum technique.The effect of near infrared light and biological tissue mainly shows as scattering and absorption, wherein in the 600-900nm wavelength band, has tangible irrelevant absorption spectra feature, formed a biological tissue " optical diagnostics window ", (DiffuseOptical Tomography DOT) can obtain the closely-related optical parametric 3-D quantitative information of organizer deep and physiological pathology index to the diffuse optical chromatographic technique of Jian Liing on this basis.Meanwhile, along with of the deepen continuously application of optical molecular iconography in biological study, be based upon combine with specific molecular fluorescent labeling technology optical molecular chromatographic technique on the basis of DOT method and then become the research field of optical molecular imaging forefront, comprise fluorescent chromatographic (Fluorescence Molecular Tomography, FMT) and the archebiosis photosphere analyse (Bioluminescence Tomography, BLT) two kinds of basic mode.FMT is by simultaneous excitation light-fluorescence measurement and fluorescent diffusion chromatography process on exciting light and two wave bands of emitting fluorescence, produce the spatial distribution of the fluorescence parameter of inner special molecular biochemical process intensity of reflection organizer and microenvironment feature thereof, comprise emissivity (fluorescent agent quantum efficiency and absorptance are long-pending) and life-span etc.Difference according to mode of excitation, FMT has three kinds of measurement patterns: time domain (Time Domain, TD), frequency domain (FrequencyDomain, FD) and continuous wave (Continuous Wave, CW), wherein simple and direct in CW mode principle and the realization, and FD and TD mode are expected to realize higher function and index.Though can realize that for the FD mode multiparameter rebuilds simultaneously, but when using this mode in small animal imaging, need the above modulating frequency of 1GHz to reach suitable signal to noise ratio, it is almost nil to exchange amplitude under the so high frequency, the limitation of visible FD mode.Therefore the time domain FMT technology that possesses multiparameter imaging capability and function admirable has obtained to pay attention to widely.Time domain FMT technology is intended to develop high sensitive time resolved detection technique, excite the measurement of lower surface fluorescence transient state " projection " and based on the image inversion algorithm of accurate photon transport model by many " angles ", rebuild the spatial distribution of the inner many fluorescence parameters of specific molecular target of complex organization's body (fluorescent yield and life-span), and the finiteness that remedies spatial sampling quantity by time resolution information realizes the application requirements of FMT with effective raising image quality and detection sensitivity.
(Time-Correlated Single PhotonCounting, TCSPC) technology are a kind of discrete channel measuring techniques with ultra-high sensitive, reasonable time resolution to the existing employing of time domain FMT implementation time correlation single photon meter; Adopt time gated image to strengthen CCD camera (Time-gated Intensifier CCD camera) in addition, it is a kind of measuring technique with high spatial sample rate, but its linearity and dynamic range are not high, and compare the TCSPC technology and do not have advantage yet on sensitivity and temporal resolution.
Summary of the invention
For overcoming the deficiencies in the prior art, remedy the defective of prior art on spatial sampling density, make system with ultra-sensitivity and high time-space resolution measurement performance, realize high-quality multicomponent, multiparameter imaging function then.For achieving the above object, the technical solution used in the present invention is: imitating CT scanning mode multispectral time domain fluorescence molecular tomography measuring system system comprises:
Two the psec semiconductor lasers and the fluorescent probe of required wavelength ultrashort laser are provided, are used to send laser pulse, the wavelength of one of them laser instrument is at the exciting light wave band of fluorescent probe, and another wavelength is at the fluorescence wave band of fluorescent probe;
Be used for the same frequency ultrashort laser pulse of two kinds of different wave lengths was coupled at interval by regular time the fiber coupler of beam of laser;
Be used for two psec semiconductor lasers or the shoot laser of one of them are projected the incident optical of objective body;
Be used to realize the imaging chamber and the rotation platform of imitating CT scanning mode;
Be used to receive from the reflection of objective body or the reception optical fiber of transmission laser;
Be used for the optical filter that filtering or decling phase stress be luminous and realize that the motor that the different performance optical filter switches mutually drives filter wheel;
Be used to receive detection module, comprise: photoelectron counting PMT detector set, 8 passage routing modules, multidimensional time correlation single photon meter TCSPC module by emergent light after the respective filter.
It is respectively corresponding relatively with exciting with emission spectra of different fluorescent dyes that two psec semiconductor laser operation wavelengths and filter set wavelength are provided with, as follows according to the combination that different fluorescent probes carry out:
1) Cy5.5 near infrared fluorescent dye, its peak excitation wavelength is about 670nm, the peak fluorescence wavelength is about 700nm, the wavelength of two psec semiconductor lasers then selecting for use is respectively 670nm and 700nm, there is halfwidth to be spaced apart 10nm in the filter set for 10nm centre wavelength, wave-length coverage is 8 bandpass filters and long pass filter sheet that cutoff wavelength is 690nm of 680-760nm;
2) indocyanine green Indocyanine Green, the ICG near infrared fluorescent dye, its peak excitation wavelength is about 780nm, the peak fluorescence wavelength is about 830nm, the wavelength of two psec semiconductor lasers then selecting for use is respectively 780nm and 830nm, have halfwidth to be spaced apart 10nm for 10nm centre wavelength in the filter set, wave-length coverage is 8 bandpass filters of 790-870nm, and long pass filter sheet that cutoff wavelength is 800nm;
3) Alexa Fluor 750 near infrared fluorescent dyes, its peak excitation wavelength is about 750nm, the peak fluorescence wavelength is about 785nm, the wavelength of two psec semiconductor lasers then selecting for use is respectively 750nm and 785nm, there is halfwidth to be spaced apart 10nm's in the filter set for 10nm centre wavelength, wave-length coverage is 8 bandpass filters of 760-840nm, and the long pass filter sheet that cutoff wavelength is 770m.
Two psec semiconductor lasers are connected to incident optical by a bonder after the optical fiber of different length produces the regular time interval, receive optical fiber and be 8 detection optical fibers or fibre bundle in the isoplanar on incident optical opposite location arrangements, the imaging chamber places on the rotating-/ lifting platform, by the stepped intervals control of the rotating-/ lifting platform anglec of rotation and adjustable height, realize the time discrimination measurement of different spaces sampling density.
8 detection optical fibers or fibre bundle are corresponding to respectively one in 8 the 12 hole filter wheel, a corresponding passage that inserts in 8 PMT sense channels, a hole is a full impregnated in the 12 hole filter wheel, remaining hole is provided with the fluorescent belt pass filter group of long pass filter group of fluorescence and different centre wavelengths respectively, and be used for the neutral density ND attenuator that exciting light detects, be used to carry out the measurement of exciting light, full gloss spectrum fluorescence or multispectral fluorescence signal; 8 PMT detectors insert multidimensional time correlation single photon meter TCSPC module by one 8 passage router.
Two psec semiconductor lasers are correlated with control to realize exporting synchronously the light pulse of same frequency by computer, produce regular time through the optical fiber of different length again and form beam of laser by bonder in the back at interval, to realize measuring the background tissues optical parameter value that obtains under this two wavelength by the DOT before the fluorescence measurement.
The boundary face of incident optical and reception optical fiber and objective body Ω
Figure GDA0000021224100000031
Certain distance is arranged, be positioned at r dLuminous flux Γ (the r that the fibre-optical probe of place's reception optical fiber records d, t) with imaging chamber emissivity
Figure GDA0000021224100000032
Satisfy following ray theory relation:
Figure GDA0000021224100000033
Wherein r is Position vector arbitrarily on the surface,
Figure GDA0000021224100000035
With
Figure GDA0000021224100000036
Be respectively the unit vector on the normal direction unit vector and fibre-optical probe direction outside the surface, imaging chamber, V (r d, r) be r point and r dVisibility function between the point, μ A0The aerial attenuation quotient of light, R ARadius for cylindrical target.
System of the present invention adopts optical fiber coupling noncontact mode to realize the measurement of time domain DOT/FMT prototype, and its characteristics are:
1, imaging chamber of the present invention and structure of fiber_optic are inside and outside concentric design, and postanesthetic tumor-bearing mice is vertically put into the imaging chamber with slight extruding during measurement, and fill the gap with optical match liquid.Two are respectively the picosecond pulse laser device that excites with wavelength of fluorescence and are connected to incident optical by a bonder, in the isoplanar on this opposite location arrangements 8 detection optical fibers or fibre bundle are arranged.The imaging chamber places on the rotating-/ lifting platform, stepped intervals control by the anglec of rotation and adjustable height, can realize the time discrimination measurement of different spaces sampling density, the data that obtain can be carried out the monoplane two-dimension image rebuild in groups by Different Plane, or the 3 D stereo reconstruction is carried out in comprehensive utilization.
2, the present invention can realize the multispectral survey technology, outgoing optical fiber inserts 8 PMT sense channels respectively, one 12 hole filter wheel of coupling before each PMT, one of them hole is a full impregnated, remaining hole is provided with the fluorescent belt pass filter group of long pass filter group of fluorescence and different centre wavelengths respectively, and be used for neutral density (ND) attenuator that exciting light detects, can carry out the measurement of exciting light, full gloss spectrum fluorescence or multispectral fluorescence signal.
3, the present invention has adopted single channel multidimensional TCSPC measurement pattern, and 8 PMT detectors insert multidimensional TCSPC module by one 8 passage router, carry out the high speed time discrimination measurement.
4, incident and outgoing optical fiber are the non-cpntact measurement pattern among the present invention, as shown in Figure 6 the boundary face of luminous source optical fiber and detection fiber and objective body Ω
Figure GDA0000021224100000037
Certain distance is arranged, promptly use optical fiber noncontact mode incentives target body and sensed light signal.Be positioned at r dLuminous flux Γ (the r that the fibre-optical probe at place records d, t) with imaging chamber emissivity
Figure GDA0000021224100000038
Satisfy following ray theory relation:
Figure GDA0000021224100000041
Wherein r is
Figure GDA0000021224100000042
Position vector arbitrarily on the surface,
Figure GDA0000021224100000043
With
Figure GDA0000021224100000044
Be respectively the unit vector on the normal direction unit vector and fibre-optical probe direction outside the surface, imaging chamber, V (r d, r) be r point and r dVisibility function between the point, μ A0The aerial attenuation quotient of light, R ARadius for cylindrical target.
5, two of light-source system employings of the present invention are respectively the picosecond pulse laser device that excites with wavelength of fluorescence and are connected to incident optical by a bonder.Two wavelength have regular time at interval during common the use, measure the background tissues optical parameter value that can obtain under this two wavelength by the DOT before the fluorescence measurement.The priori value of background during as the fluorescent chromatographic image reconstruction with this optical parametric can improve the resolution and the quantified precision of fluorescent chromatographic image significantly.
6, non-cpntact measurement pattern of the present invention, not only effective simplified measurement process is got rid of unnecessary systematic error and loaded down with trivial details coupling calibration process, and is made data acquisition become versatile and flexible.Can realize that high-density acquisition and key area emphasize by discrete angle intervals of gated sweep and scope, thereby effectively remedy the deficiency of TCSPC technology on spatial sampling density, make system with ultra-sensitivity and high time-space resolution measurement performance, realize high-quality multicomponent, multiparameter imaging function then.
Description of drawings
Fig. 1 is an imitating CT scanning mode multispectral time domain fluorescence molecular tomography measuring system system block diagram.
The basic structure sketch map of Fig. 2 light-source system.
Fig. 3 motor drives filter wheel and photomultiplier tube group structural representation.
Fig. 4 motor drives optical filter branch sketch map in the filter wheel.
Fig. 5 imaging chamber and scan control scheme sketch map.
Fig. 6 optical fiber non-cpntact measurement pattern chromatography floor map.
The specific embodiment
The contactless pattern of development imitating CT scanning realizes having the meaning of particular importance in the time domain FMT measuring system based on the TCSPC technology: on the one hand, the same with all non-cpntact measurement patterns, it is the simplified measurement process effectively, gets rid of unnecessary systematic error and loaded down with trivial details coupling calibration process; On the other hand, non-contact measurement makes data acquisition become versatile and flexible, can realize that high-density acquisition and key area emphasize by discrete angle intervals of gated sweep and scope, thereby effectively remedy the deficiency of TCSPC technology on spatial sampling density, make system with ultra-sensitivity and high time-space resolution measurement performance, realize high-quality multicomponent, multiparameter imaging function then.
The present invention is an imitating CT scanning mode, i.e. FMT/DOT measuring system under the optical fiber noncontact mode.Fig. 1 is the basic structure of imitating CT scanning mode multispectral time domain fluorescence molecular tomography measuring system system.It is by light-source system 1, incident optical 2, and imaging chamber 3, outgoing optical fiber 4, detection system, computer 11 etc. are partly formed.FMT/DOT measuring process under the concrete optical fiber noncontact mode is described below.
A. when Cy5.5 is fluorescent dye, the measuring process of this system (two optical maser wavelengths that are system are respectively 670nm and 700nm).
1, by before the fluorescence measurement, utilizes two wavelength to carry out the DOT measurement earlier, thereby obtain the background tissues optical parameter value under this two wavelength.The wavelength of excitation wavelength picopulse semiconductor laser 13 is 670nm as shown in Figure 2, and the wavelength of wavelength of fluorescence picopulse semiconductor laser 14 is 700nm.Two laser instrument switches are controlled by related software, and its frequency is identical and receive the synchronizing signal end of multidimensional TCSPC module 10 synchronously by power combiner 12.Two synchronous laser that laser instrument sent are connected to fiber coupler 15 with the optical fiber by different length, and (synchronous signal frequency is elected 40MHZ as, and two wavelength laser pulse intervals are elected 12.5ns as.Because optical fiber internal diameter material is quartzy (refractive index is 1.54), therefore two laser instrument should differ to the fiber lengths of fiber coupler and be about 2.4 meters), thus make two kinds of different wave length laser blend together at interval with a branch of laser emitting to incident optical 2 (core diameter 62.5 μ m) with regular time.
2, incident optical 2 is beaten laser in imaging chamber 3 and (is made by translucent lucite, the about 2-3cm of diameter) on, imaging chamber 3 is inside and outside concentric design with structure of fiber_optic 45 (diameter 4cm) as shown in Figure 5, postanesthetic tumor-bearing mice is vertically put into imaging chamber 3 with slight extruding during measurement, and with optical match liquid filling gap, 101.25 °~258.75 ° the position, isoplanar on this opposite (establishing incident illumination in position, 0 ° of angle) is furnished with 8 outgoing optical fiber 4 (detection optical fibers, core diameter 500 μ m), the angle, space is 22.5 °.The imaging chamber places on the rotating-/ lifting platform 46, and the stepped intervals control by the anglec of rotation and adjustable height realizes imitating CT scanning mode.To survey optical signal by outgoing optical fiber 4 and import detection module.
3,8 outgoing optical fiber 4 insert motor and drive optical filtering wheels 16-23, and motor drives the filter wheel structure as shown in Figure 4, and its mesopore 33 is full impregnated via holes; Hole 34 is the neutral-density filters that decay to 30dB; Hole 35 is the neutral-density filters that decay to 50dB; Before measuring fluorescence, carry out DOT and measure, then drive filter wheel, transfer to the optical signal after suitably being decayed in hole 33/34/35 by filter wheel controller 24 governor motors.
4, the optical signal that as shown in Figure 3 motor is driven the outgoing of optical filtering wheels is sent into the PMT detector 25-32 that has refrigerating function, and the low light level conversion of signals during PMT will test is an electronic pulse signal.8 PMT detectors are linked into multidimensional TCSPC module 10 by 8 passage routing modules 9, electronic impulse to input is counted, and deposit corresponding memory block in by the channel address signal, obtain the temporal extension curve of 8 dual wavelength exciting lights of the test surface outgoing that records from 8 passages.
5, after dual wavelength DOT measured, carrying out the FMT measurement was fluorescence measurement, only opens excitation wavelength picopulse semiconductor laser 13 (wavelength is 670nm).To imaging chamber 3, the stepped intervals control by rotating-/ lifting platform 46 anglecs of rotation and adjustable height realizes imitating CT scanning mode to laser through incident optical 2.Import detection system by outgoing optical fiber 4 again.Drive filter wheel by filter wheel controller 24 governor motors, its mesopore 36 is the long pass filter of 690nm for cut-off frequency; Hole 37-44 is halfwidth is spaced apart 10nm for 10nm centre wavelength 8 bandpass filters (wave-length coverage is 680-760nm); Obtain full gloss spectrum fluorescence or multispectral fluorescence signal by different choice to hole 36-44.Low light level conversion of signals in will testing by PMT again is an electronic pulse signal.8 PMT detectors are linked into multidimensional TCSPC module 10 by 8 passage routing modules 9, electronic impulse to input is counted and is deposited corresponding memory block in by the channel address signal, thereby obtains 8 fluorescence temporal extension curves of the test surface outgoing that records from 8 passages.
6, the each rotation in an imaging chamber angle (generally electing 1 ° as) repeats above-mentioned measuring process, after 360 ° of rotations, by repeating the process of wheel measuring after the step-length of automatically controlled lifting platform stepping (being generally 1mm) again.The generally long 5cm of imaging cylinder, so the stepping scope of lifting platform is greater than 5cm.Be that general available detection data amount is 360 * 50 * 8=144000.Owing to can select more accurate automatically controlled rotation and lifting platform for use, it is not difficult further to improve the spatial discrimination of surveying.
B. when ICG is fluorescent dye, the measuring process of this system (two optical maser wavelengths that are system are respectively 780nm and 830nm).
1, by before the fluorescence measurement, utilizes two wavelength to carry out the DOT measurement earlier, thereby obtain the background tissues optical parameter value under this two wavelength.The wavelength of excitation wavelength picopulse semiconductor laser 13 is 780nm as shown in Figure 2, and the wavelength of wavelength of fluorescence picopulse semiconductor laser 14 is 830nm.Two laser instrument switches are controlled by related software, and its frequency is identical and receive the synchronizing signal end of multidimensional TCSPC module 10 synchronously by power combiner 12.Two synchronous laser that laser instrument sent are connected to fiber coupler 15 with the optical fiber by different length, and (synchronous signal frequency is elected 40MHZ as, and two wavelength laser pulse intervals are elected 12.5ns as.Because optical fiber internal diameter material is quartzy (refractive index is 1.54), therefore two laser instrument should differ to the fiber lengths of fiber coupler and be about 2.4 meters), thus make two kinds of different wave length laser blend together at interval with a branch of laser emitting to incident optical 2 (core diameter 62.5 μ m) with regular time.
2, incident optical 2 is beaten laser in imaging chamber 3 and (is made by translucent lucite, the about 2-3cm of diameter) on, imaging chamber 3 is inside and outside concentric design with structure of fiber_optic 45 (diameter 4cm) as shown in Figure 5, postanesthetic tumor-bearing mice is vertically put into imaging chamber 3 with slight extruding during measurement, and with optical match liquid filling gap, 101.25 °~258.75 ° the position, isoplanar on this opposite (establishing incident illumination in position, 0 ° of angle) is furnished with 8 outgoing optical fiber 4 (detection optical fibers, core diameter 500 μ m), the angle, space is 22.5 °.The imaging chamber places on the rotating-/ lifting platform 46, and the stepped intervals control by the anglec of rotation and adjustable height realizes imitating CT scanning mode.To survey optical signal by outgoing optical fiber 4 and import detection module.
3,8 outgoing optical fiber 4 insert motor and drive optical filtering wheels 16-23, and motor drives the filter wheel structure as shown in Figure 4, and its mesopore 33 is full impregnated via holes; Hole 34 is the neutral-density filters that decay to 30dB; Hole 35 is the neutral-density filters that decay to 50dB; Before measuring fluorescence, carry out DOT and measure, then drive filter wheel, transfer to the optical signal after suitably being decayed in hole 33/34/35 by filter wheel controller 24 governor motors.
4, the optical signal that as shown in Figure 3 motor is driven the outgoing of optical filtering wheels is sent into the PMT detector 25-32 that has refrigerating function, and the low light level conversion of signals during PMT will test is an electronic pulse signal.8 PMT detectors are linked into multidimensional TCSPC module 10 by 8 passage routing modules 9, electronic impulse to input is counted, and deposit corresponding memory block in by the channel address signal, obtain the temporal extension curve of 8 dual wavelength exciting lights of the test surface outgoing that records from 8 passages.
5, after dual wavelength DOT measured, carrying out the FMT measurement was fluorescence measurement, only opens excitation wavelength picopulse semiconductor laser 13 (wavelength is 780nm).To imaging chamber 3, the stepped intervals control by rotating-/ lifting platform 46 anglecs of rotation and adjustable height realizes imitating CT scanning mode to laser through incident optical 2.Import detection system by outgoing optical fiber 4 again.Drive filter wheel by filter wheel controller 24 governor motors, its mesopore 36 is the long pass filter of 800nm for cut-off frequency; Hole 37-44 is halfwidth is spaced apart 10nm for 10nm centre wavelength 8 bandpass filters (wave-length coverage is 790-870nm); Obtain full gloss spectrum fluorescence or multispectral fluorescence signal by different choice to hole 36-44.Low light level conversion of signals in will testing by PMT again is an electronic pulse signal.8 PMT detectors are linked into multidimensional TCSPC module 10 by 8 passage routing modules 9, electronic impulse to input is counted and is deposited corresponding memory block in by the channel address signal, thereby obtains 8 fluorescence temporal extension curves of the test surface outgoing that records from 8 passages.
6, the each rotation in an imaging chamber angle (generally electing 1 ° as) repeats above-mentioned measuring process, after 360 ° of rotations, by repeating the process of wheel measuring after the step-length of automatically controlled lifting platform stepping (being generally 1mm) again.The generally long 5cm of imaging cylinder, so the stepping scope of lifting platform is greater than 5cm.Be that general available detection data amount is 360 * 50 * 8=144000.Owing to can select more accurate automatically controlled rotation and lifting platform for use, it is not difficult further to improve the spatial discrimination of surveying.
C. when Alexa Fluor 750 is fluorescent dye, the measuring process of this system (two optical maser wavelengths that are system are respectively 750nm and 785nm).
1, by before the fluorescence measurement, utilizes two wavelength to carry out the DOT measurement earlier, thereby obtain the background tissues optical parameter value under this two wavelength.The wavelength of excitation wavelength picopulse semiconductor laser 13 is 750nm as shown in Figure 2, and the wavelength of wavelength of fluorescence picopulse semiconductor laser 14 is 785nm.Two laser instrument switches are controlled by related software, and its frequency is identical and receive the synchronizing signal end of multidimensional TCSPC module 10 synchronously by power combiner 12.Two synchronous laser that laser instrument sent are connected to fiber coupler 15 with the optical fiber by different length, and (synchronous signal frequency is elected 40MHZ as, and two wavelength laser pulse intervals are elected 12.5ns as.Because optical fiber internal diameter material is quartzy (refractive index is 1.54), therefore two laser instrument should differ to the fiber lengths of fiber coupler and be about 2.4 meters), thus make two kinds of different wave length laser blend together at interval with a branch of laser emitting to incident optical 2 (core diameter 62.5 μ m) with regular time.
2, incident optical 2 is beaten laser in imaging chamber 3 and (is made by translucent lucite, the about 2-3cm of diameter) on, imaging chamber 3 is inside and outside concentric design with structure of fiber_optic 45 (diameter 4cm) as shown in Figure 5, postanesthetic tumor-bearing mice is vertically put into imaging chamber 3 with slight extruding during measurement, and with optical match liquid filling gap, 101.25 °~258.75 ° the position, isoplanar on this opposite (establishing incident illumination in position, 0 ° of angle) is furnished with 8 outgoing optical fiber 4 (detection optical fibers, core diameter 500 μ m), the angle, space is 22.5 °.The imaging chamber places on the rotating-/ lifting platform 46, and the stepped intervals control by the anglec of rotation and adjustable height realizes imitating CT scanning mode.To survey optical signal by outgoing optical fiber 4 and import detection module.
3,8 outgoing optical fiber 4 insert motor and drive optical filtering wheels 16-23, and motor drives the filter wheel structure as shown in Figure 4, and its mesopore 33 is full impregnated via holes; Hole 34 is the neutral-density filters that decay to 30dB; Hole 35 is the neutral-density filters that decay to 50dB; Before measuring fluorescence, carry out DOT and measure, then drive filter wheel, transfer to the optical signal after suitably being decayed in hole 33/34/35 by filter wheel controller 24 governor motors.
4, the optical signal that as shown in Figure 3 motor is driven the outgoing of optical filtering wheels is sent into the PMT detector 25-32 that has refrigerating function, and the low light level conversion of signals during PMT will test is an electronic pulse signal.8 PMT detectors are linked into multidimensional TCSPC module 10 by 8 passage routing modules 9, electronic impulse to input is counted, and deposit corresponding memory block in by the channel address signal, obtain the temporal extension curve of 8 dual wavelength exciting lights of the test surface outgoing that records from 8 passages.
5, after dual wavelength DOT measured, carrying out the FMT measurement was fluorescence measurement, only opens excitation wavelength picopulse semiconductor laser 13 (wavelength is 750nm).To imaging chamber 3, the stepped intervals control by rotating-/ lifting platform 46 anglecs of rotation and adjustable height realizes imitating CT scanning mode to laser through incident optical 2.Import detection system by outgoing optical fiber 4 again.Drive filter wheel by filter wheel controller 24 governor motors, its mesopore 36 is the long pass filter of 770nm for cut-off frequency; Hole 37-44 is halfwidth is spaced apart 10nm for 10nm centre wavelength 8 bandpass filters (wave-length coverage is 760-840nm); Obtain full gloss spectrum fluorescence or multispectral fluorescence signal by different choice to hole 36-44.Low light level conversion of signals in will testing by PMT again is an electronic pulse signal.8 PMT detectors are linked into multidimensional TCSPC module 10 by 8 passage routing modules 9, electronic impulse to input is counted and is deposited corresponding memory block in by the channel address signal, thereby obtains 8 fluorescence temporal extension curves of the test surface outgoing that records from 8 passages.
6, the each rotation in an imaging chamber angle (generally electing 1 ° as) repeats above-mentioned measuring process, after 360 ° of rotations, by repeating the process of wheel measuring after the step-length of automatically controlled lifting platform stepping (being generally 1mm) again.The generally long 5cm of imaging cylinder, so the stepping scope of lifting platform is greater than 5cm.Be that general available detection data amount is 360 * 50 * 8=144000.Owing to can select more accurate automatically controlled rotation and lifting platform for use, it is not difficult further to improve the spatial discrimination of surveying.

Claims (6)

1.一种仿CT扫描模式多光谱时域荧光分子层析测量系统,其特征是,包括:1. A multi-spectral time-domain fluorescent molecular tomography measurement system imitating CT scan mode, characterized in that it comprises: 提供所需波长超短激光的两个皮秒半导体激光器及荧光探针,用于发出激光脉冲,其中一个激光器的波长在荧光探针的激发光波段,另一个波长在荧光探针的荧光波段;Provide two picosecond semiconductor lasers and fluorescent probes with ultra-short wavelength lasers for emitting laser pulses, one of the lasers has a wavelength in the excitation light band of the fluorescent probe, and the other has a wavelength in the fluorescence band of the fluorescent probe; 用于将两种不同波长的同频率超短激光脉冲按固定的时间间隔耦合成一束激光的光纤耦合器;A fiber coupler used to couple two ultrashort laser pulses of the same frequency with different wavelengths into one laser beam at a fixed time interval; 用于将两个皮秒半导体激光器或其中一个的出射激光投射到目标体的入射光纤;An input optical fiber for projecting the output laser light of two picosecond semiconductor lasers or one of them to the target; 用于实现仿CT扫描模式的成像腔及旋转平台;Imaging cavity and rotating platform for imitating CT scanning mode; 用于接收来自目标体的反射或透射激光的接收光纤;A receiving fiber for receiving reflected or transmitted laser light from the target; 用于滤除或衰减相应激发光的滤光片及实现不同性能滤光片相互切换的马达驱动滤光轮;Filters for filtering or attenuating the corresponding excitation light and motor-driven filter wheels for switching between filters with different performances; 用于接收通过相应滤光片后出射光的检测模块,包括:光电子计数(PMT)检测器组、8通道路由模块、多维时间相关单光子计(TCSPC)模块。The detection module used to receive the emitted light after passing through the corresponding filter includes: a photoelectron counting (PMT) detector group, an 8-channel routing module, and a multidimensional time-correlated single photon meter (TCSPC) module. 2.根据权利要求1所述的一种仿CT扫描模式多光谱时域荧光分子层析测量系统,其特征是,两个皮秒半导体激光器工作波长以及滤光片组波长设置分别是与不同荧光染料的激发和发射光谱相对应,根据不同荧光探针进行的组合如下:2. A kind of imitative CT scanning mode multi-spectral time-domain fluorescence molecular tomography measurement system according to claim 1, it is characterized in that, two picosecond semiconductor laser operating wavelengths and filter set wavelengths are respectively set with different fluorescence The excitation and emission spectra of the dyes correspond to each other, and the combination of different fluorescent probes is as follows: 1)Cy5.5近红外荧光染料,其峰值激发波长约为670nm,峰值荧光波长约为700nm,则选用的两个皮秒半导体激光器的波长分别为670nm和700nm,滤光片组中有半高宽为10nm,中心波长间隔为10nm,波长范围是680-760nm的8个带通滤光片,和一个截止波长为690nm的长通滤波片;1) Cy5.5 near-infrared fluorescent dye has a peak excitation wavelength of about 670nm and a peak fluorescence wavelength of about 700nm. The wavelengths of the two picosecond semiconductor lasers selected are 670nm and 700nm respectively, and there are half-height 8 band-pass filters with a width of 10nm, a center wavelength interval of 10nm, a wavelength range of 680-760nm, and a long-pass filter with a cut-off wavelength of 690nm; 2)吲哚氰绿(Indocyanine Green,ICG)近红外荧光染料,其峰值激发波长约为780nm,峰值荧光波长约为830nm,则选用的两个皮秒半导体激光器的波长分别为780nm和830nm,滤光片组中有半高宽为10nm,中心波长间隔为10nm,波长范围是790-870nm的8个带通滤光片,及一个截止波长为800nm的长通滤波片;2) Indocyanine Green (Indocyanine Green, ICG) near-infrared fluorescent dye has a peak excitation wavelength of about 780nm and a peak fluorescence wavelength of about 830nm. In the light sheet group, there are 8 band-pass filters with a half-height width of 10nm, a center wavelength interval of 10nm, and a wavelength range of 790-870nm, and a long-pass filter with a cut-off wavelength of 800nm; 3)Alexa Fluor 750近红外荧光染料,其峰值激发波长约为750nm,峰值荧光波长约为785nm,则选用的两个皮秒半导体激光器的波长分别为750nm和785nm,滤光片组中有半高宽为10nm,中心波长间隔为10nm,波长范围是760-840nm的8个带通滤光片,及一个截止波长为770m的长通滤波片。3) Alexa Fluor 750 near-infrared fluorescent dye has a peak excitation wavelength of about 750nm and a peak fluorescence wavelength of about 785nm. The wavelengths of the two picosecond semiconductor lasers selected are 750nm and 785nm respectively, and there are half-height The width is 10nm, the center wavelength interval is 10nm, 8 bandpass filters with a wavelength range of 760-840nm, and a long-pass filter with a cutoff wavelength of 770m. 3.根据权利要求1所述的一种仿CT扫描模式多光谱时域荧光分子层析测量系统,其特征是,两个皮秒半导体激光器经不同长度的光纤产生固定的时间间隔后,通过一个耦合器连接到入射光纤,接收光纤为在入射光纤对面的同平面位置布置的8个探测光纤或光纤束,成像腔置于一个旋转/升降台上,通过旋转/升降台旋转角度和升降高度的步进间隔控制,实现不同空间采样密度的时间分辨测量。3. A kind of imitation CT scanning mode multi-spectral time-domain fluorescence molecular tomography measurement system according to claim 1, it is characterized in that, after two picosecond semiconductor lasers produce fixed time intervals through optical fibers of different lengths, pass through a The coupler is connected to the incident fiber, and the receiving fiber is 8 detection fibers or fiber bundles arranged on the same plane opposite to the incident fiber. The imaging cavity is placed on a rotating/elevating platform, and the rotation angle and lifting height of the rotating/elevating platform are adjusted. Step interval control to realize time-resolved measurement of different spatial sampling densities. 4.根据权利要求3所述的一种仿CT扫描模式多光谱时域荧光分子层析测量系统,其特征是,8个探测光纤或光纤束分别对应通过8个12孔滤波轮中的一个,对应接入8个光电子计数(PMT)检测通道中的一个通道,12孔滤波轮中一个孔为全透,其余的孔分别设置荧光长通滤光片及不同中心波长的荧光带通滤光片组,以及用于激发光检测的中性密度(ND)衰减片组,用于进行激发光、全光谱荧光或多光谱荧光信号的测量;8个光电子计数(PMT)检测器通过一个8通道路由器接入多维时间相关单光子计(TCSPC)模块。4. A kind of imitative CT scanning mode multi-spectral time-domain fluorescence molecular tomography measurement system according to claim 3, is characterized in that, 8 detection optical fibers or optical fiber bundles correspond to one of the 8 12-hole filter wheels respectively, Corresponding to one of the 8 photoelectron counting (PMT) detection channels, one hole in the 12-hole filter wheel is fully transparent, and the remaining holes are respectively equipped with fluorescent long-pass filters and fluorescent band-pass filters of different central wavelengths group, and a neutral density (ND) attenuation film set for excitation light detection, used to measure excitation light, full-spectrum fluorescence, or multi-spectrum fluorescence signals; 8 photoelectron counting (PMT) detectors through an 8-channel router Access to the multidimensional time-correlated single photon meter (TCSPC) module. 5.根据权利要求1所述的一种仿CT扫描模式多光谱时域荧光分子层析测量系统,其特征是,两个皮秒半导体激光器由计算机进行相关控制以实现同步输出相同频率的光脉冲,再经过不同长度的光纤产生固定的时间间隔后通过耦合器形成一束激光,以实现通过荧光测量前的扩散光学层析技术(Diffuse Optical Tomography,DOT)测量获得该两波长下的背景组织光学参数值。5. A multi-spectral time-domain fluorescence molecular tomography measurement system imitating CT scanning mode according to claim 1, characterized in that two picosecond semiconductor lasers are controlled by a computer to achieve synchronous output of light pulses of the same frequency , and then through the optical fiber of different lengths to generate a fixed time interval, a beam of laser light is formed through the coupler, so as to achieve the measurement of the background tissue optics at the two wavelengths by the Diffuse Optical Tomography (DOT) measurement before the fluorescence measurement. parameter value. 6.根据权利要求1所述的一种仿CT扫描模式多光谱时域荧光分子层析测量系统,其特征是,入射光纤和接收光纤与目标体Ω的边界面
Figure FDA0000046967340000021
有一定距离,位于rd处接收光纤的光纤探头测得的光流量Γ(rd,t)与成像腔表面辐射率
Figure FDA0000046967340000022
满足下列射线理论关系:
6. A kind of CT-like scanning mode multi-spectral time-domain fluorescence molecular tomography measurement system according to claim 1, characterized in that, the incident optical fiber and receiving optical fiber and the boundary surface of the target body Ω
Figure FDA0000046967340000021
With a certain distance, the optical flow Γ(r d , t) measured by the optical fiber probe located at the receiving optical fiber at r d is related to the surface radiance of the imaging cavity
Figure FDA0000046967340000022
The following ray theory relations are satisfied:
Figure FDA0000046967340000023
其中r为
Figure FDA0000046967340000024
表面上任意的位置矢量,
Figure FDA0000046967340000025
Figure FDA0000046967340000026
分别为成像腔表面外法向单位矢量和光纤探头方向上的单位矢量,V(rd,r)为r点与rd点之间的视见函数,μa0为光在空气中的衰减系数,辐射率
Figure FDA0000046967340000027
中,
Figure FDA0000046967340000028
表示指向某个方向单位向量;t表示时间参量。
Figure FDA0000046967340000023
where r is
Figure FDA0000046967340000024
an arbitrary position vector on the surface,
Figure FDA0000046967340000025
and
Figure FDA0000046967340000026
are the normal unit vector outside the surface of the imaging cavity and the unit vector in the direction of the fiber optic probe, V(r d , r) is the visibility function between point r and point r d , μ a0 is the attenuation coefficient of light in air , emissivity
Figure FDA0000046967340000027
middle,
Figure FDA0000046967340000028
Represents a unit vector pointing to a certain direction; t represents a time parameter.
CN 201010167623 2010-05-11 2010-05-11 Multi-spectral time-domain fluorescence molecular tomography measurement system imitating CT scanning mode Expired - Fee Related CN101816552B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN 201010167623 CN101816552B (en) 2010-05-11 2010-05-11 Multi-spectral time-domain fluorescence molecular tomography measurement system imitating CT scanning mode

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN 201010167623 CN101816552B (en) 2010-05-11 2010-05-11 Multi-spectral time-domain fluorescence molecular tomography measurement system imitating CT scanning mode

Publications (2)

Publication Number Publication Date
CN101816552A CN101816552A (en) 2010-09-01
CN101816552B true CN101816552B (en) 2011-06-22

Family

ID=42651948

Family Applications (1)

Application Number Title Priority Date Filing Date
CN 201010167623 Expired - Fee Related CN101816552B (en) 2010-05-11 2010-05-11 Multi-spectral time-domain fluorescence molecular tomography measurement system imitating CT scanning mode

Country Status (1)

Country Link
CN (1) CN101816552B (en)

Families Citing this family (19)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102327111A (en) * 2011-09-23 2012-01-25 天津析像光电科技有限公司 Multichannel near-infrared brain functional imaging parallel detection system
CN102525420B (en) * 2011-12-16 2013-07-10 天津大学 Calibration method for multi-passage time domain fluorescence chromatography imaging system
CN102579012A (en) * 2012-03-15 2012-07-18 天津大学 Dynamic fluorescence-optical breast tomography device on basis of phase locking photon counting measurement
CN103356165A (en) * 2012-03-29 2013-10-23 锐珂(上海)医疗器材有限公司 Multispectral fluorescence three-dimensional tomography method and device
CN102764138B (en) * 2012-08-02 2014-07-02 北京大学 Multi-mode little animal molecular image imaging device and imaging method
CA2929677C (en) * 2013-11-13 2023-10-31 Danmarks Tekniske Universitet Method for surface scanning in medical imaging and related apparatus
CN105286797B (en) * 2015-11-03 2018-06-08 北京大学 A kind of vertical rotating fluorescent molecular tomography system and its imaging method
WO2017173315A1 (en) * 2016-04-01 2017-10-05 Black Light Surgical, Inc. Systems, devices, and methods for time-resolved fluorescent spectroscopy
US10578463B2 (en) * 2016-12-19 2020-03-03 Microchip Technology Incorporated Detection of defects in motor position decoder system
DE102018202239A1 (en) * 2018-02-14 2019-08-14 Dr. Johannes Heidenhain Gmbh Bearing angle measuring device
CN108433704B (en) * 2018-04-10 2024-05-14 西安维塑智能科技有限公司 Three-dimensional human body scanning equipment
TWI666445B (en) * 2018-07-26 2019-07-21 緯創資通股份有限公司 Detection device and detection cage module thereof
CN109040698B (en) * 2018-09-07 2020-12-22 上海联影医疗科技股份有限公司 A monitoring system and method for use in medical equipment
CN109580439B (en) * 2018-10-10 2023-12-22 金华职业技术学院 An ultrafast spectroscopy experimental method for microparticles
CN111449623B (en) * 2020-03-26 2021-11-02 天津大学 Subdiffuse tissue domain spatially resolved optical measurement system for rapid diagnosis of cervical cancer
CN112229827B (en) * 2020-09-07 2022-02-08 南京大学 Real-time multispectral tomography method and device
CN112130316B (en) * 2020-10-27 2022-02-01 华北电力大学 Multi-channel multi-spectral-band optical filter structure and application and method thereof
CN115962874B (en) * 2022-11-21 2025-01-28 兰州空间技术物理研究所 A rolling torque measurement device for electric thrusters based on laser induced fluorescence
CN116849624B (en) * 2023-08-31 2023-11-10 南京诺源医疗器械有限公司 Image sensor fluorescence imaging method and system based on 4CMOS

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8357351B2 (en) * 2004-04-21 2013-01-22 Ananth Annapragada Nano-scale contrast agents and methods of use
CN100469324C (en) * 2006-09-29 2009-03-18 清华大学 Data acquisition system for small animal molecular imaging with nuclide and fluorescence dual-mode integration
JP2011506325A (en) * 2007-12-05 2011-03-03 マーバル バイオサイエンシーズ インコーポレイテッド Nanoscale contrast agent and method of use thereof
CN101236160B (en) * 2008-03-04 2011-06-01 天津大学 Time Domain Fluorescence Diffusion Tomography System for Molecular Imaging of Small Animals

Also Published As

Publication number Publication date
CN101816552A (en) 2010-09-01

Similar Documents

Publication Publication Date Title
CN101816552B (en) Multi-spectral time-domain fluorescence molecular tomography measurement system imitating CT scanning mode
US7242997B2 (en) Diffuse optical tomography system and method of use
US9134243B2 (en) System and method for sub-surface fluorescence imaging
JP6035006B2 (en) Method and apparatus for diffusion excitation in imaging
US6825928B2 (en) Depth-resolved fluorescence instrument
CN105358947B (en) Integration spectrum for Raman, reflection and fluorescence spectral measuring is popped one's head in
US9635349B2 (en) Second generation hand held optical imager
Sharma et al. Design and characterization of a novel multimodal fiber-optic probe and spectroscopy system for skin cancer applications
US8330962B2 (en) System and method for optical coherence tomography with modulated detector sensitivity
CN102551671B (en) Photon counting-type dynamic diffusion fluorescence tomography method and device
US20050143663A1 (en) Depth-resolved fluorescence instrument with angled excitation
Lapointe et al. A multi-view time-domain non-contact diffuse optical tomography scanner with dual wavelength detection for intrinsic and fluorescence small animal imaging
He et al. Novel endoscopic optical diagnostic technologies in medical trial research: recent advancements and future prospects
CN101612034B (en) Time-resolved measurement system and method for reconstructing optical parameters of turbid media
JP2007528500A (en) Methods and systems for tomographic imaging using fluorescent proteins
CN102892348A (en) Method and device for multi-spectral photonic imaging
CN102499645A (en) Photoacoustic and fluorescence dual-mode integrated tomography imaging system and imaging method
CN108956564B (en) Photosensitizer concentration detection device, system and method
US8964183B2 (en) Systems and methods for screening of biological samples
CN201042433Y (en) Optical detection device for living tissue
CN101303315B (en) Optical fiber common focusing micro spectrum and imaging apparatus of cell analysis
CN101940463B (en) Fluorescence endoscopic imaging system for living body
KR100968611B1 (en) Hybrid Imaging Device Using Double Cladding Fiber Optic Device
CN1953701A (en) Method and system for tomographic imaging using fluorescent proteins
CN207163913U (en) A kind of optical microphotograph imaging system

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20110622

Termination date: 20210511