CN108335338B - Experimental animal multi-mode fusion imaging system and using method - Google Patents
Experimental animal multi-mode fusion imaging system and using method Download PDFInfo
- Publication number
- CN108335338B CN108335338B CN201810268753.9A CN201810268753A CN108335338B CN 108335338 B CN108335338 B CN 108335338B CN 201810268753 A CN201810268753 A CN 201810268753A CN 108335338 B CN108335338 B CN 108335338B
- Authority
- CN
- China
- Prior art keywords
- information
- experimental animal
- module
- fmt
- mri
- 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.)
- Active
Links
- 238000010171 animal model Methods 0.000 title claims abstract description 138
- 238000000034 method Methods 0.000 title claims abstract description 80
- 238000003384 imaging method Methods 0.000 title claims abstract description 75
- 230000004927 fusion Effects 0.000 title claims abstract description 58
- 238000012634 optical imaging Methods 0.000 claims abstract description 45
- 238000001914 filtration Methods 0.000 claims abstract description 24
- 230000003287 optical effect Effects 0.000 claims abstract description 23
- 241001465754 Metazoa Species 0.000 claims abstract description 22
- 238000012545 processing Methods 0.000 claims abstract description 15
- 238000003325 tomography Methods 0.000 claims abstract description 14
- 238000012360 testing method Methods 0.000 claims abstract description 13
- 230000005284 excitation Effects 0.000 claims description 59
- 239000011159 matrix material Substances 0.000 claims description 55
- 238000013519 translation Methods 0.000 claims description 19
- 239000007850 fluorescent dye Substances 0.000 claims description 11
- 239000013307 optical fiber Substances 0.000 claims description 10
- 239000000126 substance Substances 0.000 claims description 9
- 238000002939 conjugate gradient method Methods 0.000 claims description 7
- 238000010276 construction Methods 0.000 claims description 7
- 238000013507 mapping Methods 0.000 claims description 7
- 238000004364 calculation method Methods 0.000 claims description 4
- 238000001514 detection method Methods 0.000 claims description 3
- 239000011521 glass Substances 0.000 claims description 3
- 238000005457 optimization Methods 0.000 claims description 3
- 230000000149 penetrating effect Effects 0.000 claims description 3
- 238000011109 contamination Methods 0.000 claims 1
- 230000008569 process Effects 0.000 abstract description 20
- 238000002595 magnetic resonance imaging Methods 0.000 description 48
- 210000001519 tissue Anatomy 0.000 description 15
- 238000012937 correction Methods 0.000 description 10
- 238000005516 engineering process Methods 0.000 description 9
- 230000000877 morphologic effect Effects 0.000 description 8
- 241000699666 Mus <mouse, genus> Species 0.000 description 7
- 210000001015 abdomen Anatomy 0.000 description 6
- 210000000038 chest Anatomy 0.000 description 6
- XEBWQGVWTUSTLN-UHFFFAOYSA-M phenylmercury acetate Chemical compound CC(=O)O[Hg]C1=CC=CC=C1 XEBWQGVWTUSTLN-UHFFFAOYSA-M 0.000 description 6
- 238000010586 diagram Methods 0.000 description 5
- 238000009792 diffusion process Methods 0.000 description 5
- 230000005540 biological transmission Effects 0.000 description 4
- 238000009543 diffuse optical tomography Methods 0.000 description 4
- 239000000835 fiber Substances 0.000 description 4
- 210000000056 organ Anatomy 0.000 description 4
- 238000005481 NMR spectroscopy Methods 0.000 description 3
- 239000002390 adhesive tape Substances 0.000 description 3
- 238000002474 experimental method Methods 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 230000003595 spectral effect Effects 0.000 description 3
- IGXWBGJHJZYPQS-SSDOTTSWSA-N D-Luciferin Chemical compound OC(=O)[C@H]1CSC(C=2SC3=CC=C(O)C=C3N=2)=N1 IGXWBGJHJZYPQS-SSDOTTSWSA-N 0.000 description 2
- CYCGRDQQIOGCKX-UHFFFAOYSA-N Dehydro-luciferin Natural products OC(=O)C1=CSC(C=2SC3=CC(O)=CC=C3N=2)=N1 CYCGRDQQIOGCKX-UHFFFAOYSA-N 0.000 description 2
- BJGNCJDXODQBOB-UHFFFAOYSA-N Fivefly Luciferin Natural products OC(=O)C1CSC(C=2SC3=CC(O)=CC=C3N=2)=N1 BJGNCJDXODQBOB-UHFFFAOYSA-N 0.000 description 2
- 108060001084 Luciferase Proteins 0.000 description 2
- 239000005089 Luciferase Substances 0.000 description 2
- DDWFXDSYGUXRAY-UHFFFAOYSA-N Luciferin Natural products CCc1c(C)c(CC2NC(=O)C(=C2C=C)C)[nH]c1Cc3[nH]c4C(=C5/NC(CC(=O)O)C(C)C5CC(=O)O)CC(=O)c4c3C DDWFXDSYGUXRAY-UHFFFAOYSA-N 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 2
- 210000004556 brain Anatomy 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 238000002591 computed tomography Methods 0.000 description 2
- 230000005281 excited state Effects 0.000 description 2
- 210000002216 heart Anatomy 0.000 description 2
- 210000003734 kidney Anatomy 0.000 description 2
- 210000004185 liver Anatomy 0.000 description 2
- 210000004072 lung Anatomy 0.000 description 2
- 230000011218 segmentation Effects 0.000 description 2
- 230000000007 visual effect Effects 0.000 description 2
- 241000699670 Mus sp. Species 0.000 description 1
- 206010028980 Neoplasm Diseases 0.000 description 1
- 230000003187 abdominal effect Effects 0.000 description 1
- 238000010521 absorption reaction Methods 0.000 description 1
- 210000003484 anatomy Anatomy 0.000 description 1
- 210000000988 bone and bone Anatomy 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000013170 computed tomography imaging Methods 0.000 description 1
- 238000002596 diffuse optical imaging Methods 0.000 description 1
- 201000010099 disease Diseases 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- 230000036267 drug metabolism Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000012632 fluorescent imaging Methods 0.000 description 1
- 230000005283 ground state Effects 0.000 description 1
- 238000001727 in vivo Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000013421 nuclear magnetic resonance imaging Methods 0.000 description 1
- 230000001737 promoting effect Effects 0.000 description 1
- 230000001902 propagating effect Effects 0.000 description 1
- 239000000523 sample Substances 0.000 description 1
- 210000002356 skeleton Anatomy 0.000 description 1
- 239000000758 substrate Substances 0.000 description 1
- 230000008685 targeting Effects 0.000 description 1
- 210000000115 thoracic cavity Anatomy 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
- G06T11/003—Reconstruction from projections, e.g. tomography
- G06T11/005—Specific pre-processing for tomographic reconstruction, e.g. calibration, source positioning, rebinning, scatter correction, retrospective gating
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/0033—Features or image-related aspects of imaging apparatus classified in A61B5/00, e.g. for MRI, optical tomography or impedance tomography apparatus; arrangements of imaging apparatus in a room
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/0033—Features or image-related aspects of imaging apparatus classified in A61B5/00, e.g. for MRI, optical tomography or impedance tomography apparatus; arrangements of imaging apparatus in a room
- A61B5/0035—Features or image-related aspects of imaging apparatus classified in A61B5/00, e.g. for MRI, optical tomography or impedance tomography apparatus; arrangements of imaging apparatus in a room adapted for acquisition of images from more than one imaging mode, e.g. combining MRI and optical tomography
-
- 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/0071—Measuring for diagnostic purposes; Identification of persons using light, e.g. diagnosis by transillumination, diascopy, fluorescence by measuring fluorescence emission
-
- 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/0073—Measuring for diagnostic purposes; Identification of persons using light, e.g. diagnosis by transillumination, diascopy, fluorescence by tomography, i.e. reconstruction of 3D images from 2D projections
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/05—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves
- A61B5/055—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves involving electronic [EMR] or nuclear [NMR] magnetic resonance, e.g. magnetic resonance imaging
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
- A61B6/02—Devices for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
- A61B6/03—Computerised tomographs
- A61B6/032—Transmission computed tomography [CT]
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
- A61B6/44—Constructional features of apparatus for radiation diagnosis
- A61B6/4417—Constructional features of apparatus for radiation diagnosis related to combined acquisition of different diagnostic modalities
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
- A61B6/50—Clinical applications
- A61B6/508—Clinical applications for non-human patients
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61D—VETERINARY INSTRUMENTS, IMPLEMENTS, TOOLS, OR METHODS
- A61D3/00—Appliances for supporting or fettering animals for operative purposes
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/50—Image enhancement or restoration by the use of more than one image, e.g. averaging, subtraction
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/30—Determination of transform parameters for the alignment of images, i.e. image registration
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B2503/00—Evaluating a particular growth phase or type of persons or animals
- A61B2503/40—Animals
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B2503/00—Evaluating a particular growth phase or type of persons or animals
- A61B2503/42—Evaluating a particular growth phase or type of persons or animals for laboratory research
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61D—VETERINARY INSTRUMENTS, IMPLEMENTS, TOOLS, OR METHODS
- A61D3/00—Appliances for supporting or fettering animals for operative purposes
- A61D2003/003—Appliances for supporting or fettering animals for operative purposes with head or neck restraining means
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10072—Tomographic images
- G06T2207/10081—Computed x-ray tomography [CT]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10072—Tomographic images
- G06T2207/10084—Hybrid tomography; Concurrent acquisition with multiple different tomographic modalities
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10072—Tomographic images
- G06T2207/10088—Magnetic resonance imaging [MRI]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20212—Image combination
- G06T2207/20221—Image fusion; Image merging
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Medical Informatics (AREA)
- Veterinary Medicine (AREA)
- Public Health (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- General Health & Medical Sciences (AREA)
- Animal Behavior & Ethology (AREA)
- Surgery (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Molecular Biology (AREA)
- Pathology (AREA)
- Biophysics (AREA)
- Radiology & Medical Imaging (AREA)
- High Energy & Nuclear Physics (AREA)
- Theoretical Computer Science (AREA)
- Optics & Photonics (AREA)
- General Physics & Mathematics (AREA)
- Dentistry (AREA)
- Oral & Maxillofacial Surgery (AREA)
- Animal Husbandry (AREA)
- Wood Science & Technology (AREA)
- Zoology (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Pulmonology (AREA)
- Investigating, Analyzing Materials By Fluorescence Or Luminescence (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
The utility model provides a multimode fusion imaging system of experimental animal and a using method, comprising: an optical imaging unit for acquiring optical tomography information of an experimental animal body, comprising: the FMT module is used for collecting FMT information of an experimental animal body, and the FMT module reduces pollution of exciting light reflected light to collected signals by using a mirror filtering method; the MRI unit is arranged at a distance from the optical imaging unit and is used for acquiring MRI information of the experimental animal body; the fixed bed is used for fixing a test animal body and is respectively matched with the optical imaging unit and the MRI unit so that the test animal body can move rigidly between the optical imaging unit and the MRI unit; and the data processing unit is used for carrying out pixel-level fusion on the optical tomography information and the MRI information of the experimental animal body. According to the experimental animal multi-mode fusion imaging system and the using method, exciting light reflected light in emitted light is removed through a mirror filtering method, and then the problem of pollution of the reflected light of the exciting light to the emitted light in the FMT collection process is solved.
Description
Technical Field
The disclosure relates to the technical field of medical molecular imaging, in particular to an experimental animal multi-mode fusion imaging system and a using method thereof.
Background
The medical tomography technology for small animals is a tomography technology which takes experimental animals (such as mice, or other animals with the volume of imaging organs below 5cm by 5 cm) as imaging objects. The three-dimensional distribution and the shape of the focus are obtained by shooting and reconstructing a disease model established on an experimental animal body. The corresponding biomedical experiments are assisted, thereby being beneficial to the research of relevant aspects. The technology has wide application in tumor research, drug metabolism research and other aspects.
Magnetic Resonance Imaging (MRI) and diffuse optical tomography are common techniques for medical tomography of experimental animals. The nuclear magnetic resonance imaging technology utilizes the characteristic that protons in tissues generate nuclear magnetic resonance phenomenon after receiving radio frequency excitation pulse excitation under a uniform magnetic field to generate nuclear magnetic resonance signals, and high-resolution anatomical tomography imaging is further carried out on the tissues by acquiring the nuclear magnetic resonance signals. Among the diffuse optical tomography techniques, bio-autofluorescence computed tomography (BLT) and excitation fluorescence computed tomography (FMT) are more commonly used. The BLT imaging technology utilizes the chemical reaction between luciferase in an organism and a Luciferin substrate (Luciferin) injected into the organism to generate fluorescence, generates light spots on the surface of the organism and is received by a detector, and further obtains the spatial distribution information of the luciferase in the organism. FMT imaging techniques utilize targeted fluorescent probes to bring fluorophores into targeted biological tissue. Fluorophores in biological tissue are irradiated by an exogenous laser (i.e., excitation light) and their electrons transition to an excited state. Then, in the process of returning the electrons from the excited state to the ground state, fluorescence is released (i.e., light is emitted). After the fluorescence propagates in the organism, a light spot is generated on the surface of the organism and is received by the detector, so that the spatial distribution information of the fluorescence targeting probe in the organism is obtained.
However, in the process of implementing the present disclosure, the inventors found that since fluorescence is absorbed and scattered by various tissues in a living body while passing through the living body, a light spot propagating to the surface of the living body has been severely distorted, which makes it difficult to perform reconstruction imaging of morphological information of biological tissues in the living body by diffuse optical tomography. Moreover, although the FMT imaging technology selects the excitation light and the emission light in different spectral bands, the emission light is collected and filtered by using an optical filter. But the excitation light still pollutes the emitted light and affects the accuracy of the biological tissue morphological reconstruction. In addition, in the magnetic resonance imaging technology, an imaging object needs to be in a uniform magnetic field during imaging, and the magnetic field can affect electronic equipment such as a camera used for diffuse optical tomography imaging, so that other modalities cannot be acquired.
BRIEF SUMMARY OF THE PRESENT DISCLOSURE
Technical problem to be solved
Based on the technical problems, the present disclosure provides an experimental animal multimode fusion imaging system and a use method thereof, so as to alleviate the technical problem that excitation light pollutes emitted light and affects the accuracy of biological tissue morphological reconstruction in the FMT imaging technology in the prior art.
(II) technical scheme
According to one aspect of the present disclosure, there is provided a multimode fusion imaging system for laboratory animals, comprising: an optical imaging unit for acquiring optical tomography information of an experimental animal body, comprising: the FMT module is used for collecting FMT information of an experimental animal body, and the FMT module reduces pollution of reflected light of exciting light to collected signals by using a mirror filtering method; the MRI unit is arranged at a distance from the optical imaging unit and is used for acquiring MRI information of the experimental animal body; the fixed bed is used for fixing a test animal body and is respectively matched with the optical imaging unit and the MRI unit so that the test animal body can move rigidly between the optical imaging unit and the MRI unit; and the data processing unit is used for carrying out pixel-level fusion on the optical tomography information and the MRI information of the experimental animal body.
In some embodiments of the present disclosure, the mirror filtering method includes: calculating mirror filter wavelength lambdaf,λf=2λ1-λ2(ii) a Respectively collecting the wavelengths of the optical filters as lambda2And λfTime of excitation of fluorescent surface light spot result b2And bf(ii) a And calculating the final surface fluorescence spot bp=b2-α*bf(ii) a Wherein λ is1Is the wavelength of the excitation light, λ2To send outThe wavelength of the emitted light is determined,b2nto use a wavelength of lambda without using a fluorescent probe2The light spot result on the surface of the fluorescence excited by the light filter bfnTo use a wavelength of lambda without using a fluorescent probefThe obtained result of exciting the fluorescent surface light spot is collected when the optical filter is used.
In some embodiments of the present disclosure, the optical imaging unit further comprises: the BLT module is used for collecting BLT information of the experimental animal; and the CT module is used for acquiring CT information of the experimental animal.
In some embodiments of the present disclosure, the data processing unit includes: the system matrix construction module is used for constructing a system matrix by utilizing FMT information and BLT information; the reconstruction module is used for constructing a regularization matrix by adopting a Gaussian weight Laplace matrix method, and calculating the light intensity distribution of the fluorescent light source in the experimental animal body by using a conjugate gradient method in combination with a system matrix; the registration module is used for registering the reconstruction result of the FMT information and the BLT information and the MRI information by the reconstruction module by taking the CT information as a reference; and the fusion module is used for performing pixel-level fusion on the reconstruction results of the FMT information and the BLT information and the MRI information by utilizing the mapping relation of the FMT information, the BLT information and the MRI information obtained by the registration module.
In some embodiments of the present disclosure, the system matrix construction module comprises: the correction submodule corrects the FMT system matrix by adopting an excitation light threshold correction method, and the excitation light threshold correction method comprises the following steps: acquisition of phix(j) Maximum value of (1)xmax(j) (ii) a Let all phix(j)<0.4φxmax(j) Phi ofx(j) Are all set to 0.4 phixmax(j) (ii) a Correcting the FMT system matrix according to the following formula;
wherein phi isx(j) Is exciting light in the body of the experimental animalIntensity distribution of fluorescent light, [ phi ]m(i) Eta mu to measure the intensity distribution of the obtained surface fluorescence spotaf(j) Is the distribution of fluorescence light source in the body of the experimental animal.
In some embodiments of the present disclosure, the regularization matrix L is constructed by the reconstruction module by using a laplacian matrix with gaussian weights as follows:
LG=(li,j)N×N
wherein lijRepresents an attenuating element, dijRepresenting the spatial distance between two voxels, ρ SkFor the radial softening function, R is the gaussian kernel radius, and in the gaussian weighted laplace matrix method, R is 0.4.
In some embodiments of the present disclosure, the calculating the light intensity distribution of the fluorescent light source inside the experimental animal body by using a conjugate gradient method in combination with the system matrix in the reconstruction module includes: step A: establishing a target optimization function;
wherein the content of the first and second substances,is a least-squares equation of the equation,the regularization term is a Gihono Voff regularization term, lambda is a regularization parameter, L is the regularization matrix, and A is the system matrix; and B: taking an initial value x(0)C, calculating an initial residual vector r(0)=b-Ax(0)And let p stand for(0)=r(0)Wherein p is(0)Is the initial iteration direction;
and C: if r(k)||2C is less than or equal to c, let x be x(k)The calculation terminates, otherwise the following is calculated:
x(k+1)=x(k)+tk*p(k)
wherein r is(k)Is the residual vector for the k-th iteration,step size, p, for the k-th iteration(k)Is the k iteration direction;
and correcting the obtained elements in x by the following formula:
wherein the content of the first and second substances,representing the result of the kth iteration of the ith element of vector x,is the corrected result vector;
step D: calculating r(k+1)=b-Ax(k+1) And let p stand for(k+1)=r(k+1)+αkp(k)
step E: and C, returning the step C to judge if k is k + 1.
In some embodiments of the present disclosure, the fixed bed comprises: the two sides of the bed board are provided with baffles, and the bed board and the baffles are used for placing experimental animals; the optical imaging window is correspondingly arranged in the middle of the bed plate and on the baffle and is used for collecting imaging information of the experimental animal from the bottom and the side surface; the fixing interface is arranged on the bed board and/or the baffle and is used for assisting in fixing the experimental animal; the positioning point is arranged on the bed plate and/or the baffle and is used for positioning the optical imaging unit or the MRI unit; wherein, the bed board with the baffle contains organic glass.
In some embodiments of the present disclosure, the optical imaging unit includes: the BLT module, the FMT module and the CT module are all fixedly arranged on the plane turntable; the translation table reciprocates along the normal vector direction of the plane turntable and is detachably connected with the bed plate and/or the baffle; the translation table moves along the normal vector direction of the plane turntable, so that the optical imaging window on the fixed bed is exposed in the imaging range of the BLT module, the FMT module and the CT module, the BLT module, the FMT module and the CT module are driven to rotate through the plane turntable, and the optical tomography information of the experimental animal is acquired.
In some embodiments of the disclosure, wherein: the BLT module is a fluorescence camera and is used for collecting BLT information of the experimental animal body; the FMT module comprises: excitation light emitting means for emitting excitation light, comprising: the first rotating device is connected with the plane turntable and drives the rocker arm to rotate in a plane parallel to the plane turntable; the second rotating device is connected with the first rotating device through the rocker arm, and the rotating shaft of the second rotating device is vertical to the normal vector of the plane turntable; the optical fiber translation clamp is connected with the second rotating device, stretches along the rotating shaft of the second rotating device and forms a W-degree included angle with the rotating shaft of the second rotating device, wherein W is more than or equal to 0 degree and less than or equal to 15 degrees; the laser emergent port is connected with the optical fiber translation clamp, and is provided with a speed expanding head wheel for emitting a point light source or a surface light source; and a fluorescence camera for collecting FMT information of the experimental animal; the CT module includes: an X-ray bulb for emitting X-rays; the X-ray detection plate and the X-ray bulb tube are respectively arranged on two sides of the translation table and used for receiving X-rays penetrating through the experimental animal body; wherein the BLT module and the FMT module share the same fluorescence camera.
In some embodiments of the present disclosure, the MRI unit includes: a nuclear magnetic cavity for collecting MRI information; and the experimental animal bracket is detachably connected with the bed board and/or the baffle and drives the fixed bed to stretch into the nuclear magnetic cavity for nuclear magnetic imaging.
According to another aspect of the present disclosure, there is also provided a method for using a multimode fusion imaging system for a laboratory animal, the multimode fusion imaging system for a laboratory animal provided by the present disclosure includes: step 100: fixing a test animal to the optical imaging unit through the fixed bed; step 200: respectively collecting BLT information, FMT information and CT information of the experimental animal body, and correcting the FMT information by using a mirror filtering method; step 300: fixing a test animal body to the MRI unit through the fixed bed; step 400: collecting MRI information of the experimental animal; and step 500: performing pixel-level fusion on the BLT information, the FMT information, and the MRI information using the data processing unit.
In some embodiments of the present disclosure, collecting FMT information of the experimental animal body comprises: collecting FMT information by using a transmission model and selecting point laser; and/or collecting FMT information by using a reflection model and selecting planar laser; wherein the mirror filtering method is used when FMT information is collected using a reflection model.
In some embodiments of the present disclosure, selecting a spot laser to acquire FMT information using a transmission model comprises: the experimental animal body is a mouse, and the four limbs of the mouse are fixed by a wire through the fixing interface, so that the chest and abdomen of the mouse are exposed at the optical imaging window; adjusting the first rotating device to enable the laser emergent port and the fluorescence camera to be respectively positioned at two sides of the mouse; and adjusting the telescopic distance and angle of the optical fiber translation clamp, and acquiring the thoracic and abdominal FMT information under different excitation conditions.
In some embodiments of the present disclosure, selecting a facet laser to collect FMT information using a reflection model includes: the experimental animal body is a mouse, and the neck of the experimental animal body is fixed by using an adhesive tape, so that the head of the experimental animal body is exposed at the optical imaging window; selecting a filter with a wave band of 50nm wavelength less than the exciting light; collecting reflected light spot information of the exciting light under the condition of not using a fluorescent probe; using a fluorescent probe to collect FMT information of the head; and correcting the FMT information of the head by utilizing the reflected light spot information of the exciting light in combination with a mirror filtering method.
(III) advantageous effects
According to the technical scheme, the experimental animal multimode fusion imaging system and the using method have one or part of the following beneficial effects:
(1) fitting the reflected light distribution information of the excitation light under the emission light wavelength by a mirror filtering method, namely selecting the reflected light distribution information of the excitation light collected under the spectral band far away from the emission light wavelength, wherein the reflected light of the excitation light does not enter the body of the experimental animal, so that the reflected light of the excitation light under different wavelengths belongs to a linear relation, and fitting the reflected light distribution of the excitation light polluting the emission light by the reflected light distribution information of the excitation light collected under the band far away from the emission light wavelength, so that the reflected light of the excitation light in the emission light is removed, and the pollution problem of the reflected light of the excitation light on the emission light in the FMT collection process;
(2) the ill-posed property of a system matrix under a planar excitation light reflection model is reduced by using an excitation light threshold correction method;
(3) the method is based on the assumption that the correlation between two voxels is reduced along with the increase of the space distance between the two voxels, a Gaussian weight Laplacian matrix is designed as a regular matrix, the problem of the smoothness of diffusion optical reconstruction is successfully solved, the morphological information of a reconstruction object is provided, and the problems that the light spots transmitted to the surface of the organism are seriously distorted due to the absorption and scattering of various tissues in the organism when fluorescence passes through the organism, and the morphological information of the biological tissues in the organism is difficult to reconstruct and image are solved;
(4) the fixed bed is arranged to provide convenience for carrying the experimental animal body between the optical imaging unit and the MRI unit, so that the experimental animal body is guaranteed to be in rigid change under the imaging visual fields of the two units, registration and fusion among multi-mode imaging results are facilitated, and the influence of a uniform strong magnetic field of the MRI unit on electronic equipment in the optical imaging unit can be relieved.
Drawings
Fig. 1 is a schematic structural diagram of a multimode fusion imaging system of an experimental animal body according to an embodiment of the present disclosure.
Fig. 2 is a schematic structural diagram of a fixed bed in the experimental animal multi-mode fusion imaging system shown in fig. 1.
Fig. 3 is a schematic front view of an excitation light emitting device in the experimental animal multi-mode fusion imaging system shown in fig. 1.
Fig. 4 is a left side view of the excitation light emitting device in the experimental animal multi-mode fusion imaging system shown in fig. 1.
Fig. 5 is a schematic top view of an excitation light emitting device in the experimental animal multi-mode fusion imaging system shown in fig. 1.
Fig. 6 is a schematic flow chart of information acquisition in a method for using the experimental animal multi-mode fusion imaging system according to the embodiment of the disclosure.
Fig. 7 is a flowchart of data processing in the experimental animal multimodal fusion imaging system according to the embodiment of the disclosure.
[ description of main reference numerals in the drawings ] of the embodiments of the present disclosure
100-an optical imaging unit;
110-a planar turntable;
120-a translation stage;
130-a fluorescent camera;
140-excitation light emitting means;
141-first rotating means; 142-a rocker arm;
143-second rotating means; 144-translating the fiber clamp;
145-laser exit port; 146-speed expanding head wheel;
150-X-ray tube;
160-X-ray detection plate;
200-an MRI unit;
210-a nuclear magnetic cavity; 220-laboratory animal body holder;
300-fixed bed;
310-bed board; 320-a baffle plate;
330-an optical imaging window; 340-a fixed interface;
350-locating points;
400-data processing unit.
Detailed Description
In the experimental animal multimode fusion imaging system and the using method provided by the embodiment of the disclosure, the reflected light of the excitation light in the emission light is removed by a mirror filtering method, so that the pollution problem of the reflected light of the excitation light to the emission light in the FMT collection process is further alleviated, and the accuracy of biological tissue morphological reconstruction is improved.
For the purpose of promoting a better understanding of the objects, aspects and advantages of the present disclosure, reference is made to the following detailed description taken in conjunction with the accompanying drawings.
Fig. 1 is a schematic structural diagram of a multimode fusion imaging system of an experimental animal body according to an embodiment of the present disclosure.
According to an aspect of the present disclosure, as shown in fig. 1, there is provided a multimode fusion imaging system for laboratory animals, comprising: an optical imaging unit 100 for collecting optical tomographic imaging information of the experimental animal; an MRI (magnetic Resonance imaging) unit 200 spaced apart from the optical imaging unit 100 for collecting MRI information of the experimental animal; a fixed bed 300 for fixing the experimental animal body, which is respectively matched with the optical imaging unit 100 and the MRI unit 200 to make the experimental animal body perform rigid movement between the two; and a data processing unit 400 for performing pixel-level fusion of the optical tomographic imaging information and the MRI information of the experimental animal.
In some embodiments of the present disclosure, as shown in fig. 1, the optical imaging unit 100 includes: the FMT module is used for collecting FMT information of an experimental animal body, and the FMT module reduces pollution of reflected light of exciting light to collected signals by using a mirror filtering method; the BLT module is used for collecting BLT information of the experimental animal; and the CT module is used for acquiring CT information of the experimental animal.
In the present disclosureIn some embodiments, a mirror filtering method includes: calculating mirror filter wavelength lambdaf,λf=2λ1-λ2(ii) a Respectively collecting the wavelengths of the optical filters as lambda2And λfTime of excitation of fluorescent surface light spot result b2And bf(ii) a And calculating the final surface fluorescence spot bp=b2-α*bf。
Wherein λ is1Is the wavelength of the excitation light, λ2In order to emit the wavelength of the light,b2nto use a wavelength of lambda without using a fluorescent probe2The light spot result on the surface of the fluorescence excited by the light filter bfnTo use a wavelength of lambda without using a fluorescent probefThe filter is used for collecting the obtained result of exciting the fluorescent surface light spot, the reflected light distribution information of the exciting light under the emission wavelength is fitted by selecting the reflected light distribution information of the exciting light collected under the spectral waveband far away from the emission wavelength through a mirror filtering method, the reflected light distribution information of the exciting light under the emission wavelength is fitted, the reflected light of the exciting light polluting the emission light is fitted through the reflected light distribution information of the exciting light collected under the waveband far away from the emission wavelength because the reflected light of the exciting light does not enter the body of the experimental animal, and therefore the problem of pollution of the reflected light of the exciting light in the emission light to the emission light in the FMT collection process is solved.
In some embodiments of the present disclosure, a data processing unit includes: the system matrix construction module is used for constructing a system matrix by utilizing FMT information and BLT information; the reconstruction module is used for constructing a regularization matrix by adopting a Gaussian weight Laplace matrix method, and calculating the light intensity distribution of the fluorescent light source in the experimental animal body by using a conjugate gradient method in combination with a system matrix; the registration module is used for registering the reconstruction result of the FMT information and the BLT information and the MRI information by the reconstruction module by taking the CT information as a reference; and the fusion module is used for performing pixel-level fusion on the reconstruction results of the FMT information and the BLT information and the MRI information by utilizing the mapping relation of the FMT information, the BLT information and the MRI information obtained by the registration module.
The specific registration steps of the registration module are as follows: obtaining a three-dimensional anatomical tomographic structure diagram of organs (brain, lung, liver, skeleton, kidney, heart and the like) by segmentation from the reconstructed CT tomographic image, wherein each segmented image of the organs obtained by segmentation is marked as a subgraph; registering the subgraph with the FMT information and the BLT information to obtain a mapping relation between the FMT information and the BLT information; and then, the CT image and the MRI information are registered, so that the mapping relation between the FMT information and the BLT information and the MRI information is obtained.
The FMT and BLT imaging methods include a process with two cores, a forward problem and a reverse problem. The forward problem is to describe a complex physical process of transmitting fluorescence from a light source to the surface of a living body by a diffusion filter equation, establish a propagation model of the light in a biological tissue, solve the model and finally establish a linear relation between a fluorescence acquisition signal on the surface of the living body and the three-dimensional distribution of the fluorescence in the living body. The linear relationship can be expressed by a linear matrix equation:
Ax=b
where a is a system matrix describing the forward problem, x represents the light intensity distribution of the fluorescent light source inside the imaged object, and b represents the light intensity distribution of the fluorescent light spot on the surface of the imaged object.
The reconstruction process of BLT and FMT can be summarized as the following steps:
respectively carrying out BLT and FMT data acquisition experiments, and acquiring the light intensity distribution b of the fluorescence light spots on the surface of the imaging object under different conditions;
constructing a system matrix A of a BLT and an FMT forward process, wherein FMT imaging in principle comprises two associated processes, excitation and emission (BLT only comprises emission), can be described by two coupled equations:
two equations are formed byx(r) coupled together, derived to give the following equation:
wherein phim(i) Eta mu to measure the intensity distribution of the obtained surface fluorescence spotaf(j) I.e. the distribution of the fluorescent light source inside the object to be imaged, phix(j) Is the intensity distribution of the fluorescence of the excitation light in vivo, and phix(j) Decays continuously as the depth of the imaged object increases (in BLT reconstruction, this value is constant). When phi isx(j) If it is too small, the matrix of the system will be ill-conditioned.
Thus, in some embodiments of the present disclosure, the system matrix construction module comprises: the correction submodule corrects the FMT system matrix by adopting an excitation light threshold correction method, and the excitation light threshold correction method comprises the following steps:
acquisition of phix(j) Maximum value of (1)xmax(j);
Let all phix(j)<0.4φxmax(j) Phi ofx(j) Are all set to 0.4 phixmax(j);
Correcting the FMT system matrix according to the following formula;
by using the excitation light threshold correction method, the ill-posed property of the system matrix under the planar excitation light reflection model is reduced.
In some embodiments of the present disclosure, the regularization matrix L is constructed in the reconstruction module by using a laplacian matrix with gaussian weights as follows:
LG=(li,j)N×N
wherein li,jRepresents an attenuating element, di,jRepresenting the spatial distance between two voxels, ρ skFor the radial softening function, R is the gaussian kernel radius, and in the gaussian weighted laplace matrix method, R is 0.4.
In some embodiments of the present disclosure, calculating the light intensity distribution of the fluorescent light source inside the experimental animal (i.e. solving the inverse problem in the FMT and BLT imaging methods using the conjugate gradient method) using the conjugate gradient method in combination with the system matrix in the reconstruction module comprises:
step A: establishing a target optimization function;
wherein the content of the first and second substances,is a least-squares equation of the equation,the regularization term is a Gihono Voff regularization term, lambda is a regularization parameter, L is the regularization matrix, and A is the system matrix;
and B: taking an initial value x(0)C, calculating an initial residual vector r(0)=b-Ax(0)And let p stand for(0)=r(0)Wherein p is(0)Is the initial iteration direction
And C: if r(k)||2C is less than or equal to c, let x be x(k)The calculation is terminated, otherwise the calculation is
x(k+1)=x(k)+tk*p(k)
Wherein r is(k)Is the residual vector for the k-th iteration,step size, p, for the k-th iteration(k)Is the k iteration direction;
according to the characteristic that the fluorescence light source value has nonnegativity, correcting the obtained x elements by the following formula:
wherein the content of the first and second substances,representing the result of the kth iteration of the ith element of vector x,is the corrected result vector;
step D: calculating r(k+1)=b-Ax(k+1)And let p stand for(k+1)=r(k+1)+αkp(k)
step E: c, enabling k to be k +1, and returning to the step C for judgment; the Gaussian weight Laplace reconstruction method is used, based on the assumption that the correlation between two voxels is reduced along with the increase of the space distance between the two voxels, a Gaussian weight Laplace matrix is designed to serve as a regular matrix, the problem of over-smooth diffusion optical reconstruction is successfully solved, morphological information of a reconstruction object is provided, and the problem that due to the fact that fluorescence is absorbed and scattered by various tissues in a living body when passing through the living body, light spots transmitted to the surface of the living body are seriously distorted, and then the morphological information of the biological tissues in the living body is difficult to reconstruct and image is solved.
Fig. 2 is a schematic structural diagram of a fixed bed in the experimental animal multi-mode fusion imaging system shown in fig. 1.
In some embodiments of the present disclosure, as shown in fig. 2, the fixed bed 300 comprises: the two sides of the bed plate 310 are provided with baffle plates 320, and the bed plate 310 and the baffle plates 320 are used for placing experimental animals; the optical imaging window 330 is correspondingly arranged in the middle of the bed plate 310 and on the baffle 320 and is used for collecting imaging information of the experimental animal from the bottom and the side; the fixing interface 340 is arranged on the bed plate 310 and/or the baffle 320 and is used for assisting in fixing the experimental animal; the positioning point 350 is arranged on the bed plate 310 and/or the baffle 320 and is used for positioning the optical imaging unit 100 or the MRI unit 200; the bed plate 310 and the baffle plate 320 contain organic glass, the connection mode of the fixed bed 300 and the optical imaging unit 100 or the MRI unit 200 is clamping connection, bolt connection or buckle connection and the like, the connection mode of collecting optical tomography information or MRI information is not affected at will, convenience is brought to carrying of the experimental animal body between the optical imaging unit 100 and the MRI unit 200 through the fixed bed 300, accordingly, the experimental animal body is guaranteed to be in rigid change under the imaging visual fields of the two units, registration and fusion between multi-mode imaging results are facilitated, and the influence of a uniform magnetic field of the MRI unit 200 on electronic equipment in the optical imaging unit 100 can be relieved.
In some embodiments of the present disclosure, as shown in fig. 1, the optical imaging unit 100 includes: the flat rotary table 110 is vertically arranged, and the BLT module, the FMT module and the CT module are all fixedly arranged on the flat rotary table 110; and a translation stage 120 which reciprocates in the normal vector direction of the plane turntable 110 and is detachably connected to the bed plate 310 and/or the baffle 320.
The translation table 120 moves along the normal vector direction of the plane turntable 110, so that the optical imaging window 330 on the fixed bed 300 is exposed in the imaging range of the BLT module, the FMT module and the CT module, the BLT module, the FMT module and the CT module are driven to rotate by the plane turntable 110, and the optical tomography information of the experimental animal is acquired.
In some embodiments of the present disclosure, as shown in fig. 1, wherein: the BLT module is a fluorescent camera 130 and is used for collecting BLT information of experimental animals; the FMT module comprises: an excitation light emitting device 140 for emitting excitation light; and a fluorescence camera 130 for collecting FMT information of the experimental animal; the CT module includes: an X-ray tube 150 for emitting X-rays; and X-ray detecting plates 160 respectively disposed at both sides of the fixed bed 300 together with the X-ray bulb 150 for receiving X-rays penetrating through the test animal; wherein the BLT module and the FMT module share the same fluorescence camera 130.
Fig. 3 is a schematic front view of an excitation light emitting device in the experimental animal multi-mode fusion imaging system shown in fig. 1. Fig. 4 is a left side view of the excitation light emitting device in the experimental animal multi-mode fusion imaging system shown in fig. 1. Fig. 5 is a schematic top view of an excitation light emitting device in the experimental animal multi-mode fusion imaging system shown in fig. 1.
In some embodiments of the present disclosure, as shown in fig. 3 to 5, the excitation light emission device 140 includes: the first rotating device 141 is connected with the plane turntable 110 and drives the rocker arm 142 to rotate in a plane parallel to the plane turntable 110; a second rotating device 143 connected to the first rotating device 141 through a swing arm 142, wherein a rotation axis of the second rotating device 143 is perpendicular to a normal vector of the plane turntable 110; the optical fiber translation clamp 144 is connected with the second rotating device 143, extends and retracts along the rotating shaft of the second rotating device 143, and forms a W-degree included angle with the rotating shaft of the second rotating device 143, wherein W is more than or equal to 0 degree and less than or equal to 15 degrees; and a laser exit 145 connected to the fiber translation clamp 144, on which a speed expanding head wheel 146 is provided for emitting a point light source or a surface light source.
In some embodiments of the present disclosure, as shown in fig. 1, the MRI unit 200 includes: a nuclear magnetic cavity 210 for collecting MRI information; and the experimental animal bracket 220 is detachably connected with the bed plate 310 and/or the baffle 320 and drives the fixed bed to extend into the nuclear magnetic cavity 210 for nuclear magnetic imaging.
According to another aspect of the present disclosure, there is also provided a method for using a multimode fusion imaging system for a laboratory animal, where the multimode fusion imaging system for a laboratory animal provided by an embodiment of the present disclosure includes: step 100: fixing the experimental animal body to the optical imaging unit 100 through the fixed bed 300; step 200: respectively collecting BLT information, FMT information and CT information of the experimental animal body, and correcting the FMT information by using a mirror filtering method; step 300: fixing the experimental animal body to the MRI unit 200 through the fixed bed 300; step 400: collecting MRI information of the experimental animal; and step 500: the BLT information, FMT information, and MRI information are pixel-level fused using the data processing unit 400.
In some embodiments of the present disclosure, collecting FMT information of the experimental animal body comprises: collecting FMT information by using a transmission model and selecting point laser; and/or collecting FMT information by using a reflection model and selecting planar laser; wherein, when the FMT information is collected by using the reflection model, a mirror filtering method is used.
In some embodiments of the present disclosure, selecting a spot laser to acquire FMT information using a transmission model comprises: the experimental animal body is a mouse, and the four limbs of the experimental animal body are fixed by a wire through the fixing interfaces 340, so that the chest and abdomen of the experimental animal body are exposed at the optical imaging window 330; adjusting the first rotating device 141 to enable the laser emitting port 145 and the fluorescence camera 130 to be respectively located at two sides of the experimental animal body; the telescopic distance (plus or minus 1 cm) and the angle (plus or minus 10 degrees) of the optical fiber translation clamp 144 are adjusted, and the chest and abdomen FMT information under different excitation conditions is collected.
In some embodiments of the present disclosure, selecting a facet laser to collect head FMT information using a reflectance model comprises: the experimental animal body is a mouse, and the neck of the experimental animal body is fixed by using an adhesive tape, so that the head of the experimental animal body is exposed at the optical imaging window 330; selecting a filter (selecting a wave band with the wavelength less than 50nm of the exciting light) far away from the wavelength of the exciting light and the wavelength of the emitted light; collecting reflected light spot information of the exciting light under the condition of not using a fluorescent probe; using a fluorescent probe to collect FMT information of the head; and correcting the FMT information of the head by utilizing the reflected light spot information of the exciting light in combination with a mirror filtering method.
Fig. 6 is a schematic view of an acquisition process in a using method of the experimental animal multi-mode fusion imaging system according to the embodiment of the disclosure.
As shown in FIG. 6, the information acquisition process of the multimode fusion imaging system of the experimental animal body according to the embodiment of the disclosure is as follows
Step 1, starting the optical imaging unit 100 and starting the MRI unit 200;
step 2, fixing the experimental animal on the fixed bed 300, and using different fixing modes according to different shooting positions, which can be specifically divided into:
step 21: if the chest and the abdomen are shot, the limbs of the subject are bound by using a medical suture, and the thread passes through the fixing interface 340 to be bound and fixed;
step 23: if the head is shot, fixing the neck of the shot object by using a medical adhesive tape;
and step 3: after the experimental animal body is fixed, the fixed bed 300 is placed on the translation table 120 and moved to the position below the fluorescence camera 130, the camera shoots a two-dimensional picture of the experimental animal body under natural light, the proper placing position is determined, and whether all positioning points 350 on the fixed bed 300 are clear or not is determined;
and 4, step 4: closing natural light, and shooting an autofluorescence spot of the experimental animal;
and 5: securing the laser fiber in the translating fiber clamp 144 of the excitation light emitting apparatus 140;
step 6: if the chest and abdomen part is shot, entering the step 7; if the head is shot, the first rotating device 141 is adjusted to adjust the optical fiber to the same side as the fluorescence camera 130, and the reflection model shooting is performed. Entering a step 8;
and 7: if the chest and abdomen part is shot, the speed expanding head wheel 146 is adjusted, a proper speed expanding head is selected, laser irradiates a shooting target in a point laser mode, and the first rotating device 141 is adjusted to adjust the optical fiber to the opposite side of the fluorescence camera 130;
and 8: translating the optical fiber translation clamp 144 (plus or minus 1 cm) and the second rotating device 143 (plus or minus 10 degrees), and shooting the result of the fluorescent surface light spot excited by the experimental animal under different excitation conditions;
and step 9: if the head is photographed, that is, photographed using a reflection model, information on a reflection spot of excitation light necessary for mirror filtering is photographed. The method comprises the following specific steps:
step 91: selecting an optical filter (namely, a wave band with the wavelength less than 50nm of the excitation light) far away from the wavelength of the excitation light and the emission light, fixing the optical filter at the front end of the fluorescence camera 130, and filtering the emission light;
and step 92: shooting reflected light spot information of the exciting light, wherein the reflected light spot information of the exciting light is used for filtering the reflected light spot information of the exciting light mixed in the emitting light in the reflection model, and the reflected light spot information is specifically seen in a mirror filtering method;
step 10: after the fluorescent imaging acquisition is finished, the plane turntable 110 starts to rotate, and simultaneously the CT module starts to continuously acquire (1 degree and 1 piece) X-ray images of the experimental animal body, and the acquired X-ray image data is stored until the acquisition of 360 degrees is finished;
step 11: processing the acquired result, and reconstructing a three-dimensional tomographic anatomical structure (namely CT) of the shooting object by using a 360-degree X-ray image;
step 12: moving the fixed bed 300, transferring the fixed bed 300 into the experimental animal holder 220 in the MRI unit 200, and after fixing the fixed bed 300 in the experimental animal holder 220, sleeving a radio frequency coil;
step 13: the test object holder 220 is moved into the nuclear magnetic cavity 210 and MRI information of the test object is captured.
Fig. 7 is a flowchart of data processing in the experimental animal multimodal fusion imaging system according to the embodiment of the disclosure.
As shown in fig. 7, the data processing flow in the multimode fusion imaging system of the experimental animal in the embodiment of the disclosure is as follows:
step 1: selecting different acquisition modes according to the diffuse optical imaging modality to be acquired: if the FMT is the FMT, entering the step 2; otherwise, entering step 3;
step 2: in the FMT data acquisition process, if the FMT data acquisition process is head acquisition, a mirror filtering method is used, and pollution of reflected light of exciting light to light spot results on the surface of emitted light is reduced. Then entering step 4;
and step 3: in the BLT data acquisition process, an optical filter with the wavelength of 620nm is adopted to carry out a data acquisition experiment to obtain the autofluorescence intensity distribution of the surface of the organism under the wavelength;
and 4, step 4: acquiring a CT image, segmenting the acquired CT image according to organs (heart, lung, brain, bone, liver and kidney respectively), and registering the CT image with an FMT or BLT image;
and 5: after CT collection is completed, moving the shot object to an MRI module, and shooting an MRI image of the experimental animal;
step 6: after all the images are acquired, synchronously performing FMT reconstruction and BLT reconstruction, wherein if the FMT reconstruction is performed, the step 7 is performed, and otherwise, the step 9 is performed;
and 7: solving diffusion equations of an excitation process and an emission process of the FMT, and solving a forward problem to obtain a system matrix of the FMT;
and 8: correcting the fluorescence intensity distribution of the exciting light in the organism by adopting an exciting light threshold correction method so as to correct the system matrix of the FMT;
and step 9: solving a diffusion equation in the BLT transmitting process, and solving a forward problem to obtain a system matrix of the BLT;
step 10: after the system matrix construction is completed, entering a reconstruction module, namely solving the BLT and FMT reconstruction problems by adopting a Gaussian weight Laplace matrix method;
step 11: registering CT and MRI by using a mutual information registration method, and obtaining the mapping relation of BLT, FMT and MRI by registration;
step 12: and according to the mapping relation obtained by the registration of the CT and the MRI, carrying out pixel-level fusion on the BLT reconstruction result, the FMT reconstruction result and the MRI result to form the multi-modal imaging result of the experimental animal.
So far, the embodiments of the present disclosure have been described in detail with reference to the accompanying drawings. It is to be noted that, in the attached drawings or in the description, the implementation modes not shown or described are all the modes known by the ordinary skilled person in the field of technology, and are not described in detail. Further, the above definitions of the various elements and methods are not limited to the various specific structures, shapes or arrangements of parts mentioned in the examples, which may be easily modified or substituted by those of ordinary skill in the art.
From the above description, those skilled in the art should clearly recognize that the experimental animal multimodal fusion imaging system and the method of use provided by the present disclosure.
In summary, the experimental animal multimode fusion imaging system and the using method provided by the disclosure remove the excitation light reflected light in the emission light through a mirror filtering method, correct the fluorescence intensity distribution of the excitation light in the living body by adopting an excitation light threshold correction method, and solve the BLT and FMT reconstruction problems by adopting a gaussian weight laplacian matrix method, thereby alleviating the problem of pollution of the reflected light of the excitation light to the emission light in the FMT acquisition process and improving the accuracy of biological tissue morphology reconstruction.
It should also be noted that directional terms, such as "upper", "lower", "front", "rear", "left", "right", and the like, used in the embodiments are only directions referring to the drawings, and are not intended to limit the scope of the present disclosure. Throughout the drawings, like elements are represented by like or similar reference numerals. Conventional structures or constructions will be omitted when they may obscure the understanding of the present disclosure.
And the shapes and sizes of the respective components in the drawings do not reflect actual sizes and proportions, but merely illustrate the contents of the embodiments of the present disclosure. Furthermore, in the claims, any reference signs placed between parentheses shall not be construed as limiting the claim.
Similarly, it should be appreciated that in the foregoing description of exemplary embodiments of the disclosure, various features of the disclosure are sometimes grouped together in a single embodiment, figure, or description thereof for the purpose of streamlining the disclosure and aiding in the understanding of one or more of the various disclosed aspects. However, the disclosed method should not be interpreted as reflecting an intention that: that is, the claimed disclosure requires more features than are expressly recited in each claim. Rather, as the following claims reflect, disclosed aspects lie in less than all features of a single foregoing disclosed embodiment. Thus, the claims following the detailed description are hereby expressly incorporated into this detailed description, with each claim standing on its own as a separate embodiment of this disclosure.
The above-mentioned embodiments are intended to illustrate the objects, aspects and advantages of the present disclosure in further detail, and it should be understood that the above-mentioned embodiments are only illustrative of the present disclosure and are not intended to limit the present disclosure, and any modifications, equivalents, improvements and the like made within the spirit and principle of the present disclosure should be included in the scope of the present disclosure.
Claims (10)
1. A laboratory animal multimodal fusion imaging system comprising:
an optical imaging unit for acquiring optical tomography information of an experimental animal body, comprising:
an FMT module for collecting FMT information of an experimental animal, the FMT module reducing contamination of collected signals by reflected light of excitation light by using a mirror filtering method, wherein the mirror filtering method comprises:
calculating mirror filter wavelength lambdaf,λf=2λ1-λ2;
Respectively collecting the wavelengths of the optical filters as lambda2And λfTime of excitation of fluorescent surface light spot result b2And bf(ii) a And
calculating the final surface fluorescence spot bp=b2-α*bf;
Wherein λ is1Is the wavelength of the excitation light, λ2In order to emit the wavelength of the light,b2nto use a wavelength of lambda without using a fluorescent probe2The light spot result on the surface of the fluorescence excited by the light filter bfnTo use a wavelength of lambda without using a fluorescent probefThe obtained result of exciting the fluorescent surface light spot is collected when the optical filter is used;
the MRI unit is arranged at a distance from the optical imaging unit and is used for acquiring MRI information of the experimental animal body;
the fixed bed is used for fixing a test animal body and is respectively matched with the optical imaging unit and the MRI unit so that the test animal body can move rigidly between the optical imaging unit and the MRI unit; and
and the data processing unit is used for carrying out pixel-level fusion on the optical tomography information and the MRI information of the experimental animal body.
2. The experimental animal multimode fusion imaging system of claim 1, said optical imaging unit further comprising:
the BLT module is used for collecting BLT information of the experimental animal; and
and the CT module is used for acquiring CT information of the experimental animal.
3. The experimental animal multimodal fusion imaging system according to claim 2, the data processing unit comprising:
the system matrix construction module is used for constructing a system matrix by utilizing FMT information and BLT information;
the reconstruction module is used for constructing a regularization matrix by adopting a Gaussian weight Laplace matrix method, and calculating the light intensity distribution of the fluorescent light source in the experimental animal body by using a conjugate gradient method in combination with a system matrix;
the registration module is used for registering the reconstruction result of the FMT information and the BLT information and the MRI information by the reconstruction module by taking the CT information as a reference; and
and the fusion module is used for performing pixel-level fusion on the reconstruction results of the FMT information and the BLT information and the MRI information by utilizing the mapping relation of the FMT information, the BLT information and the MRI information obtained by the registration module.
4. The experimental animal multimode fusion imaging system according to claim 3, wherein the regularization matrix L constructed by the reconstruction module by using the Laplace matrix method with Gaussian weight is as follows:
LG=(li,j)N×N
wherein lijRepresents an attenuating element, dijRepresenting the spatial distance between two voxels, ρ skFor radial softening functions, R is the radius of the Gaussian kernel, in the Gaussian-weighted Laplace matrix methodAnd taking R as 0.4.
5. The multimode fusion imaging system of experimental animal body as claimed in claim 4, wherein the calculating the light intensity distribution of the fluorescent light source inside the experimental animal body by using the conjugate gradient method in combination with the system matrix in the reconstruction module comprises:
step A: establishing a target optimization function;
wherein the content of the first and second substances,is a least-squares equation of the equation,the regularization term is a Gihono Voff regularization term, lambda is a regularization parameter, L is the regularization matrix, and A is the system matrix;
and B: taking an initial value x(0)C, calculating an initial residual vector r(0)=b-Ax(0)And let p stand for(0)=r(0)Wherein p is(0)Is the initial iteration direction;
and C: if r(k)||2C is less than or equal to c, let x be x(k)The calculation terminates, otherwise the following is calculated:
x(k+1)=x(k)+tk*p(k)
wherein r is(k)Is the residual vector for the k-th iteration,step size, p, for the k-th iteration(k)Is the k iteration direction;
and correcting the obtained elements in x by the following formula:
wherein the content of the first and second substances,representing the result of the kth iteration of the ith element of vector x,is the corrected result vector;
step D: calculating r(k+1)=b-Ax(k+1)And let p stand for(k+1)=r(k+1)+αkp(k)
step E: and C, returning the step C to judge if k is k + 1.
6. The experimental animal multimode fusion imaging system of claim 2, said fixed bed comprising:
the two sides of the bed board are provided with baffles, and the bed board and the baffles are used for placing experimental animals;
the optical imaging window is correspondingly arranged in the middle of the bed plate and on the baffle and is used for collecting imaging information of the experimental animal from the bottom and the side surface;
the fixing interface is arranged on the bed board and/or the baffle and is used for assisting in fixing the experimental animal; and
the positioning point is arranged on the bed plate and/or the baffle and is used for positioning the optical imaging unit or the MRI unit;
wherein, the bed board with the baffle contains organic glass.
7. The experimental animal multimode fusion imaging system of claim 6, said optical imaging unit comprising:
the BLT module, the FMT module and the CT module are all fixedly arranged on the plane turntable; and
the translation table reciprocates along the normal vector direction of the plane turntable and is detachably connected with the bed plate and/or the baffle;
the translation table moves along the normal vector direction of the plane turntable, so that the optical imaging window on the fixed bed is exposed in the imaging range of the BLT module, the FMT module and the CT module, the BLT module, the FMT module and the CT module are driven to rotate through the plane turntable, and the optical tomography information of the experimental animal is acquired.
8. The experimental animal multimodal fusion imaging system of claim 7 wherein:
the BLT module is a fluorescence camera and is used for collecting BLT information of the experimental animal body;
the FMT module comprises:
excitation light emitting means for emitting excitation light, comprising:
the first rotating device is connected with the plane turntable and drives the rocker arm to rotate in a plane parallel to the plane turntable;
the second rotating device is connected with the first rotating device through the rocker arm, and the rotating shaft of the second rotating device is vertical to the normal vector of the plane turntable;
the optical fiber translation clamp is connected with the second rotating device, stretches along the rotating shaft of the second rotating device and forms a W-degree included angle with the rotating shaft of the second rotating device, wherein W is more than or equal to 0 degree and less than or equal to 15 degrees; and
the laser emergent port is connected with the optical fiber translation clamp, and is provided with a speed expanding head wheel for emitting a point light source or a surface light source; and
the fluorescence camera is used for collecting FMT information of the experimental animal body;
the CT module includes:
an X-ray bulb for emitting X-rays; and
the X-ray detection plates and the X-ray bulb tubes are respectively arranged on two sides of the translation table and used for receiving X-rays penetrating through the experimental animal body;
wherein the BLT module and the FMT module share the same fluorescence camera.
9. The experimental animal multimodal fusion imaging system according to any one of claims 6 to 8, the MRI unit comprising:
a nuclear magnetic cavity for collecting MRI information; and
the experimental animal bracket is detachably connected with the bed plate and/or the baffle plate and drives the fixed bed to stretch into the nuclear magnetic cavity for nuclear magnetic imaging.
10. A method for using a multimodal fusion imaging system of a laboratory animal, using the multimodal fusion imaging system of a laboratory animal according to any one of claims 1 to 9, comprising:
step 100: fixing a test animal to the optical imaging unit through the fixed bed;
step 200: respectively collecting BLT information, FMT information and CT information of the experimental animal body, and correcting the FMT information by using a mirror filtering method;
step 300: fixing a test animal body to the MRI unit through the fixed bed;
step 400: collecting MRI information of the experimental animal; and
step 500: performing pixel-level fusion on the BLT information, the FMT information, and the MRI information using the data processing unit.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810268753.9A CN108335338B (en) | 2018-03-28 | 2018-03-28 | Experimental animal multi-mode fusion imaging system and using method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810268753.9A CN108335338B (en) | 2018-03-28 | 2018-03-28 | Experimental animal multi-mode fusion imaging system and using method |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108335338A CN108335338A (en) | 2018-07-27 |
CN108335338B true CN108335338B (en) | 2021-01-05 |
Family
ID=62931609
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810268753.9A Active CN108335338B (en) | 2018-03-28 | 2018-03-28 | Experimental animal multi-mode fusion imaging system and using method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108335338B (en) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109044277B (en) * | 2018-08-10 | 2021-02-02 | 中国科学院自动化研究所 | Near-infrared two-region fluorescence tomography system |
CN109615673B (en) * | 2018-11-22 | 2021-01-05 | 中国科学院自动化研究所 | FMT reconstruction method and device based on self-adaptive Gaussian Laplace regularization |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103735252B (en) * | 2014-01-14 | 2016-01-20 | 中国科学院自动化研究所 | A kind of Optical multi-mode state imaging system and method |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9008397B2 (en) * | 2010-12-23 | 2015-04-14 | Institute of Automation, Chinese Academy of Science | Tomography system based on Cerenkov luminescence |
-
2018
- 2018-03-28 CN CN201810268753.9A patent/CN108335338B/en active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103735252B (en) * | 2014-01-14 | 2016-01-20 | 中国科学院自动化研究所 | A kind of Optical multi-mode state imaging system and method |
Non-Patent Citations (3)
Title |
---|
《Bioluminescence Tomography Based on》;Yuan Gao等;《IEEE TRANSACTIONS ON MEDICAL IMAGING》;20171130;第36卷(第11期);第2343-2354页 * |
《Dynamic Measurement of Tumor Vascular》;Wuwei Ren等;《MOLECULAR IMAGING AND BIOLOGY》;20150917;第191-200页 * |
《In vivo pentamodal tomographic imaging for》;Muhan Liu等;《Biomedical Optics EXPRESS 1356》;20170301;第8卷(第3期);第1356-1371页 * |
Also Published As
Publication number | Publication date |
---|---|
CN108335338A (en) | 2018-07-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US11659997B2 (en) | Instrument for acquiring co-registered orthogonal fluorescence and photoacoustic volumetric projections of tissue and methods of its use | |
Guo et al. | A combined fluorescence and microcomputed tomography system for small animal imaging | |
US20120302880A1 (en) | System and method for specificity-based multimodality three- dimensional optical tomography imaging | |
Cao et al. | An adaptive Tikhonov regularization method for fluorescence molecular tomography | |
Liu et al. | A parallel excitation based fluorescence molecular tomography system for whole-body simultaneous imaging of small animals | |
WO2013134949A1 (en) | Device and method for endoscopic x ray luminescence tomography imaging | |
CN110327020B (en) | Near-infrared two-zone/one-zone bimodal fluorescence tomography system | |
JP2013505468A (en) | System and method for virtual refractive index matching of diffusing media | |
JP2013505468A5 (en) | ||
Allard et al. | Combined magnetic resonance and bioluminescence imaging of live mice | |
CN108335338B (en) | Experimental animal multi-mode fusion imaging system and using method | |
CN109044277A (en) | Two area's fluorescence computed tomography (SPECT) system of near-infrared | |
CN105662354B (en) | A kind of wide viewing angle optical molecular tomographic navigation system and method | |
Zhang et al. | Three-dimensional reconstruction in free-space whole-body fluorescence tomography of mice using optically reconstructed surface and atlas anatomy | |
Lo et al. | Three-dimensional fluorescence diffuse optical tomography using the adaptive spatial prior approach | |
CN201641971U (en) | Multi-angle imaging device | |
CN114581553B (en) | Fluorescent molecular tomography reconstruction method based on magnetic particle imaging prior guidance | |
CN109615673B (en) | FMT reconstruction method and device based on self-adaptive Gaussian Laplace regularization | |
CN103750824B (en) | A kind of information extracting method for toy fluorescent chromatographic imaging system | |
Lo et al. | Small-animal 360-deg fluorescence diffuse optical tomography using structural prior information from ultrasound imaging | |
CN205795645U (en) | A kind of wide viewing angle optical molecular tomographic navigation system | |
CN114468998B (en) | Single-view angle reflection type near infrared two-region fluorescence dynamic tomography system and method | |
TWI433051B (en) | Image reconstruction method | |
CN206151439U (en) | Fluorescence molecule imaging system | |
CN117653933A (en) | Multi-mode image-guided irradiation method and system |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |