CN106037657A - Space-time-frequency multi-coupling high-density functional near-infrared spectroscopy brain function imaging method - Google Patents
Space-time-frequency multi-coupling high-density functional near-infrared spectroscopy brain function imaging method Download PDFInfo
- Publication number
- CN106037657A CN106037657A CN201610490279.5A CN201610490279A CN106037657A CN 106037657 A CN106037657 A CN 106037657A CN 201610490279 A CN201610490279 A CN 201610490279A CN 106037657 A CN106037657 A CN 106037657A
- Authority
- CN
- China
- Prior art keywords
- light source
- time
- space
- detection channels
- detector
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000010168 coupling process Methods 0.000 title claims abstract description 19
- 238000005859 coupling reaction Methods 0.000 title claims abstract description 19
- 238000003384 imaging method Methods 0.000 title claims abstract description 19
- 230000003925 brain function Effects 0.000 title claims abstract description 7
- 238000004497 NIR spectroscopy Methods 0.000 title abstract description 3
- 238000001514 detection method Methods 0.000 claims abstract description 62
- 238000000034 method Methods 0.000 claims abstract description 27
- 230000002490 cerebral effect Effects 0.000 claims description 18
- 230000015572 biosynthetic process Effects 0.000 claims description 17
- 230000008878 coupling Effects 0.000 claims description 17
- 230000006870 function Effects 0.000 claims description 17
- 238000002329 infrared spectrum Methods 0.000 claims description 15
- 230000005540 biological transmission Effects 0.000 claims description 11
- 238000012545 processing Methods 0.000 claims description 7
- 230000008569 process Effects 0.000 claims description 5
- 230000017531 blood circulation Effects 0.000 claims description 3
- 238000001914 filtration Methods 0.000 claims description 3
- 238000003475 lamination Methods 0.000 claims description 3
- 238000006243 chemical reaction Methods 0.000 claims description 2
- 241000208340 Araliaceae Species 0.000 claims 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 claims 1
- 235000003140 Panax quinquefolius Nutrition 0.000 claims 1
- 239000004744 fabric Substances 0.000 claims 1
- 235000008434 ginseng Nutrition 0.000 claims 1
- 230000008901 benefit Effects 0.000 abstract description 6
- 238000002059 diagnostic imaging Methods 0.000 abstract description 2
- 230000009286 beneficial effect Effects 0.000 abstract 1
- 238000013467 fragmentation Methods 0.000 abstract 1
- 238000006062 fragmentation reaction Methods 0.000 abstract 1
- 238000010586 diagram Methods 0.000 description 8
- 239000000523 sample Substances 0.000 description 8
- 230000002123 temporal effect Effects 0.000 description 8
- 238000003491 array Methods 0.000 description 5
- 230000009977 dual effect Effects 0.000 description 5
- 210000003128 head Anatomy 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 3
- 230000003287 optical effect Effects 0.000 description 3
- 238000005481 NMR spectroscopy Methods 0.000 description 2
- 238000010521 absorption reaction Methods 0.000 description 2
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 description 2
- 239000008280 blood Substances 0.000 description 2
- 210000004369 blood Anatomy 0.000 description 2
- 210000004556 brain Anatomy 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 230000007423 decrease Effects 0.000 description 2
- 238000009543 diffuse optical tomography Methods 0.000 description 2
- 239000000203 mixture Substances 0.000 description 2
- 229910052760 oxygen Inorganic materials 0.000 description 2
- 239000001301 oxygen Substances 0.000 description 2
- 238000012546 transfer Methods 0.000 description 2
- 102000001554 Hemoglobins Human genes 0.000 description 1
- 108010054147 Hemoglobins Proteins 0.000 description 1
- 239000011324 bead Substances 0.000 description 1
- 238000003759 clinical diagnosis Methods 0.000 description 1
- 238000002591 computed tomography Methods 0.000 description 1
- 230000001054 cortical effect Effects 0.000 description 1
- 238000009792 diffusion process Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 230000005284 excitation Effects 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 230000004060 metabolic process Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 210000004761 scalp Anatomy 0.000 description 1
- 238000002603 single-photon emission computed tomography Methods 0.000 description 1
- 210000003625 skull Anatomy 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
- 210000001519 tissue Anatomy 0.000 description 1
- 238000005303 weighing Methods 0.000 description 1
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/0059—Measuring for diagnostic purposes; Identification of persons using light, e.g. diagnosis by transillumination, diascopy, fluorescence
- A61B5/0075—Measuring for diagnostic purposes; Identification of persons using light, e.g. diagnosis by transillumination, diascopy, fluorescence by spectroscopy, i.e. measuring spectra, e.g. Raman spectroscopy, infrared absorption spectroscopy
-
- 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/004—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 image acquisition of a particular organ or body part
- A61B5/0042—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 image acquisition of a particular organ or body part for the brain
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/02—Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
- A61B5/026—Measuring blood flow
- A61B5/0261—Measuring blood flow using optical means, e.g. infrared light
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7203—Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal
Abstract
The invention discloses a space-time-frequency multi-coupling high-density functional near-infrared spectroscopy brain function imaging method, and relates to the technical field of medical imaging. The imaging method comprises steps of detection channel division, light source time and space fragmentation, intelligent allocation of light source lightening time, power and frequency and detector coding. The imaging method disclosed by the invention can relieve cross-talk among detection channels and can improve time resolution, spatial resolution and a signal-to-noise ratio; by conducting the intelligent allocation of light source lightening time, power and frequency, the method is beneficial for improving spatial resolution and signal-to-noise ratio of a system; and by coding a detector, a parallel detection frequency is improved; therefore, the method has incomparable application advantages in the high-density fNIRS system.
Description
Technical field
The present invention relates to medical imaging technology field, be specifically related to the high density near infrared light of a kind of time space frequency multiple coupling
Spectrum cerebral function formation method, is the implementation method of a kind of parallel detection data.
Background technology
In recent years, near infrared spectrum (Functional near-infrared spectroscopy, fNIRS) develops
For image mode of new generation indispensable in brain function research and clinical diagnosis.FNIRS utilize wavelength 650-1000nm with
Interior near infrared light, can detect the change of the optical absorption characteristic that cerebral metabolism causes, thus it is blood red to calculate this region deoxidation
Albumen (HbR), oxygen carrying hemoglobin (HbO2) concentration relative variation.Compared to tradition cerebral function imaging method, such as brain electricity
Figure, function NMR (Nuclear Magnetic Resonance)-imaging etc., fNIRS has cost performance height, temporal resolution is high, spatial resolution is moderate, functional parameter
Abundant, to advantages such as motion-insensitives.
In fNIRS imaging system, place light source probe at head corresponding site, and putting at a certain distance from light source probe
Put detector.The light (usually dual wavelength or three wavelength) that system source is launched, after wavelength-division multiplex, is arrived by fiber-optic transfer and puts
In the light source probe of head, penetrate scalp layer, skull etc. subsequently and incide cerebral tissue, after going through a series of absorption, scattering, still
Some photon and arrive head cortical surface, these optical signals can be received by a detector, then according to Beer-
Lambert law, calculates HbO2, the concentration relative variation of HbR.Light source probe lay down location is light source incidence point, but one
Light incident at light source incidence point comprises different wave length (usually dual wavelength or three wavelength), therefore real at a light source incidence point
There is the light source that 2 (dual wavelengths) or 3 (three wavelength) overlap on border, and for making a distinction, the multi wave length illuminating source closed of weighing is light sources,
And single wavelength light source is i.e. referred to as light source, detector lay down location is sensing point, just can reflect correspondence from the information of sensing point detection
The information of detection channels.
At present, relatively simple equidistant probe topology arrangement is widely used in multichannel fNIRS system.?
In equidistant topological imaging, adjacent light source and detector one detection channels of composition, the interval of source-detector about 3cm, it is possible to
Effectively detect the change of cortex blood oxygen parameter.Different ray detector (source-detector) arrangements, can form not
Same detection channels, multiple light sources and the arrangement of multiple detector, multi-channel system can be formed, thus obtain more detection and lead to
The brain signal in road.But, the spatial resolution of fNIRS system based on simple equidistant topological structure is the most relatively low (about
3cm), its spatial resolution still has bigger room for promotion.
Improving spatial resolution strategy is to improve sampling density, i.e. high density fNIRS formation method.With simple topology imaging
Difference, high density imaging is a kind of fault imaging mode, and for given light source, its light launched can not only be by adjacent detection
Device detects, moreover it is possible to detected by detector farther out, and therefore, a light source can different spies multiple with multiple detectors composition
Find range from detection channels.The detection channels that system uses is the most, and the coincidence between optical pickocff is the most, to different depth information
The ability of layering is the highest, and the picture quality reconstructed is the highest.In high density fNIRS, by densely arranged light source and detection
Device, is simultaneously introduced cross sectional reconstruction technology, it is possible to imaging depth rises to 2～3cm, and spatial resolution rises to about 1cm.So
And, light source, detector and detection channels needed for highdensity light source and detector arrangement make fNIRS sharply increase, but visit
Survey number of channels, signal to noise ratio and temporal resolution these three parameter mutually to restrict so that improve temporal resolution and spatial discrimination
Rate is the most challenging.
Eggebrecht et al. (Mapping distributed brain function and networks with
Diffuse optical tomography, Nature Photonics 8,448 454 (2014)) develop high density diffusion
Computed tomography (SPECT) system (high-density diffuse optical tomography, HD-DOT), have employed timesharing subregion
The method of light source excitation.This system uses dual wavelength, containing 96 light source probes, according to the spatial arrangement of ray detector, and will
96 light sources coincidence light source of two different wave lengths (each light source is actual be) are divided into 6 rectangular areas, numbered region 1～
Region 6,16, each region light source, respectively numbered light source 1～light source 16,16 light source timesharing point successively in each region
Bright, the light source of the identical numbering of zones of different is lighted simultaneously, and the light source in odd zone dual numbers region has between certain frequency
Every.But the temporal resolution of this method can only obtain limited raising, and comprises the signal cross-talk between detection channels, makes
Signal to noise ratio reduces.
Summary of the invention
Owing to detection channels quantity, signal to noise ratio and temporal resolution can not improve simultaneously, it is multiple that the present invention proposes time space frequency
The high density near infrared spectrum cerebral function formation method of coupling captures this difficult problem.The present invention according to ray detector arrange,
The parameters such as detection channels connection and connecting length, classify the detection channels of different connecting lengths, and when light source is carried out
Between coding and space encoding, the lighting time of smart allocation light source, power and frequency etc., can be effectively improved the time of imaging simultaneously
Resolution and spatial resolution, and have the highest signal to noise ratio, parallel detecting efficiency can be significantly improved, in high density fNIRS system
System has unrivaled application advantage.
The high density near infrared spectrum cerebral function formation method of a kind of time space frequency multiple coupling that the present invention provides, mainly wraps
Include following steps:
The first step, detection channels divides.According to parameters such as ray detector arrangement, detection channels connection and connecting lengths,
The detection channels of different connecting lengths is classified, generally the detection channels of equal length can be divided into a class；
Second step, light source time and spatial lamination.According to the result of detection channels classification, by light source and detection channels time
Carry out burst and sequence etc. between, and according to parameters such as space lengths, spatially carry out light source space array code, make difference
The light source of array was lighted in the different time, marked off finer timeslice, so can greatly reduce the interference between light source, and
Temporal resolution, spatial resolution and signal to noise ratio can be improved；
3rd step, light source igniting time, power and frequency smart allocation.In different connecting lengths, light source is not according to
With power drive, and the time of light source igniting also differ, connecting length is the longest, and the light source igniting time is the longest, and power is more
Greatly, be conducive to improving spatial resolution, and in same timeslice, according to the spatial distribution of light source, carry out frequency intelligent modulation, make
Must be the biggest apart from the nearest light source, frequency shift (FS), thus improve signal to noise ratio；
4th step, detector encodes.Divide according to detection channels and the result of light source time and space burst, detector is entered
Row coding, improves parallel detecting efficiency.
5th step, light source is according to the routine lights encoded, and detector also according to the programmed acquisition signal encoded, is carried out
Imaging.
The invention have the advantage that
(1) the high density near infrared spectrum cerebral function formation method of the time space frequency multiple coupling that the present invention proposes, according to light
The parameters such as source detector-arrangement, detection channels connection and connecting length carry out detection channels division and the time of light source and space
Burst, decreases the crosstalk between detection channels, can improve temporal resolution, spatial resolution and signal to noise ratio；
(2) in the time space frequency multiple coupling formation method that the present invention proposes, to light source igniting time, power and frequency intelligence
Distribution, the spatial resolution of beneficially raising system and signal to noise ratio；
(3) the time space frequency multiple coupling formation method that the present invention proposes, owing to encoding detector, can improve parallel
Detection efficient, has unrivaled application advantage in high density fNIRS system.
Accompanying drawing explanation
Figure 1A is the schematic diagram that light source probe and detector are positioned over head.
Figure 1B is ray detector arrangement schematic diagram.
The short type of attachment schematic diagram of Fig. 2 A detection channels.
Fig. 2 B is the long type of attachment schematic diagram of detection channels.
Fig. 3 A and two array of source coding schematic diagrams that Fig. 3 B is short connection.
Fig. 4 A～Fig. 4 C is long four the exemplary light sources array code schematic diagrams connected.
Fig. 5 A is detector coded system in short connection array.
Fig. 5 B connects detector coded system in array for long.
Fig. 6 is fNIRS image system hardware control figure.
Detailed description of the invention
High density near infrared spectrum cerebral function imaging to the time space frequency multiple coupling that the present invention proposes below in conjunction with the accompanying drawings
Method is described further.
The present embodiment, as a example by the ray detector of 192 light sources is arranged, has 3 different wave lengths at a light source incidence point
The light source of (785nm, 808nm, 850nm) overlaps, and therefore has 64 light sources to be distributed on 64 different light source incidence points, as
Shown in Figure 1A and Figure 1B, Figure 1B hollow core tetragon is detector, and filled black circle represents light source 1～light source 64, light source and spy
The line surveyed between device represents detection channels.Figure 1A is the schematic diagram that light source probe and detector are positioned over head, and Figure 1B is light
Source detector arrangement schematic diagram.Generally, the detection channels of high density fNIRS has four kinds of different connecting lengths, is respectively
1.3cm, 3.0cm, 3.9cm and 4.7cm, for the sake of simplicity, the present embodiment is as a example by two kinds of connecting lengths (1.3cm and 3.0cm)
The high density near infrared spectrum cerebral function formation method of the time space frequency multiple coupling of the present invention is described.
The high density near infrared spectrum cerebral function formation method of described time space frequency multiple coupling, implements step such as
Under:
The first step, detection channels divides.According to parameters such as ray detector arrangement, detection channels and connecting lengths, will not
Classify with the detection channels of connecting length, generally the detection channels of equal length can be divided into a class.
In the arrangement of this example ray detector, detection channels connecting length 1.3cm and 3.0cm respectively, by all connecting lengths
Detection channels for 1.3cm is divided into short connection, as shown in Figure 2 A.All connecting lengths are that the detection channels of 3.0cm is divided into
Long connection, division result is as shown in Figure 2 B.
Second step, light source space arrangement array code and time slicing.
According to the result of detection channels classification, light source and detection channels are carried out burst and sequence etc., and root in time
According to parameters such as space lengths, spatially carry out light source-detection channels space array coding, make the light source of different array in difference
Time light, mark off finer timeslice, so can greatly reduce the interference between light source, and time resolution can be improved
Rate, spatial resolution and signal to noise ratio.
In this example, the detection channels of short connection encodes 2 light source space arrays, as shown in Figure 3 A and Figure 3 B, each short
Connecting in light source space array, coding rule is:
Light source on same diagonal is lighted in same timeslice, it is meant that the light source of first light source space array is at certain
One timeslice is lighted simultaneously, and the light source of second light source space array is lighted in another timeslice simultaneously.This light source space array
It is distant that coded system makes between light source, and interference from each other is less, can improve signal to noise ratio.
The long detection channels connected encodes nine light source space arrays, as Fig. 4 A～4C gives three principal light source skies
Between array, the length at light source interval that these three light source space array level respectively is moved right, can be formed three new
Light source space array；On this basis, then level moves right the length at light source interval, forms last three light source skies
Between array.Move one new light source space array of formation each time, thus can form other six light source space arrays, one
Forming nine light source space arrays altogether, the most corresponding nine timeslices, different time sheet sequentially starts.Each long connection detection is logical
In the light source space array in road, coding rule is: ensure to be spaced between the light source all directions lighted two non-lit up light sources.
According to the light source space array of coding, dividing timeslice, each light source space array occupies a timeslice.This
The light source space array of 2 short connections in example and the light source space array of 9 long connections need altogether 11 timeslices, make this
11 light source space arrays are corresponding with timeslice, and in the shortest connection, light source space array 1 is lighted at time slice 1, in short connection
Light source space array 2 is lighted at time slice 2, and in long connection, light source space array 1 is lighted at time slice 3 ..., long connection
Middle light source space array 9 is lighted at time slice 11.As shown in table 1.
The high density near infrared spectrum cerebral function formation method step of table 1 time space frequency multiple coupling
3rd step, light source igniting time, power and frequency smart allocation.In different connecting lengths, light source is not according to
With power drive, and the time of light source igniting also differ, the connecting length of detection channels is the longest, and the light source igniting time is more
Long, power is the biggest, is conducive to improving spatial resolution, and in same timeslice, according to the spatial distribution of light source, carries out frequency intelligence
Can modulate so that the light source that distance is the nearest, frequency shift (FS) is the biggest, thus improves signal to noise ratio.
In short connection, as shown in figs.3 a and 3b, the power of light source is relatively low, and the light source igniting time is short；In long connection, such as Fig. 4 A
～shown in 4C, the power of light source is higher, the light source igniting time is long.The lighting time of light source and power all with connecting length approximation in
Exponential relationship, i.e. lighting time t=BeAd, power P=CeAd, wherein d is connecting length, and A, B, C are constant.
Further, in same timeslice, according to parameters such as light source space distance and detection channels connections, the light that distance is the nearest is made
Between source, frequency shift (FS) is the biggest, and the frequency shift (FS) between the most remote light source is the least, and distance about becomes to bear finger with frequency shift (FS)
Number relation.
4th step, effective detection channels and detector coding.According to detection channels classification and the result of light source burst, to having
Effect detection channels and detector encode, and improve parallel detecting efficiency.
As shown in Figure 5A, in each light source space array of short connection, light source is lighted simultaneously, and lighting time is shorter, power
Less, signal that each detector detects is mainly from the surrounding four light sources away from its nearest neighbours, from other light sources
Signal is due to distant, and decay is serious, and signal is capped, thus is easily achieved and compiles detector and target acquisition passage
Code, confirms target acquisition passage.
As shown in Figure 5 B, in the long each light source space array connected, light source is lighted simultaneously, and lighting time is longer, power
Relatively big, the 3cm passage that illuminating source is formed is destination channel.Owing to having carried out light source space array code in advance and luminous
Light source arrangement is sparse, it is easy to accomplish encode detector and target acquisition passage (effective detection channels), confirms that target is visited
Survey passage.
5th step, light source is according to the routine lights encoded, and detector also according to the programmed acquisition signal encoded, is carried out
Imaging.
After signal that detector collects is filtered and A/D changes, use FPGA (field programmable gate array) or DSP
(Digital Signal Processing) chip etc. carries out bottom layer signal process, and uses high-speed transmission interface to carry out data transmission.Then
Utilize computer PC or directly utilize the bottom layer signal processing module such as bottom FPGA, bottom DSP data are demodulated rebuild, from
And obtain effective brain function blood flow parameter.The operational order of user also can be sent to high-speed transmission interface, warp by computer
The bottom layer signal such as FPGA, dsp chip processes rear drive light source, it is achieved the real-time control to light source.
The high density near infrared spectrum cerebral function formation method of the time space frequency multiple coupling that the present invention provides, relies on fNIRS
Imaging system realizes.This fNIRS imaging system is mainly formed by with lower part: light source, detector, signal filtering and modulus of conversion
Block, bottom layer signal processing module, high-speed transmission interface and computer, as shown in Figure 6.
Light source can use laser diode or LED (light emitting diode).Detector can use PMT (photomultiplier transit
Pipe), APD (avalanche diode) or PD (photodiode) etc..After signal that detector collects is filtered and A/D changes, make
Bottom layer signal process is carried out, including number with FPGA (field programmable gate array) or DSP (Digital Signal Processing) chip etc.
According to demodulation etc., and high-speed transmission interface is used to carry out data transmission.Described high-speed transmission interface can use kilomega network, USB3.0
Deng.Then utilize computer PC or directly utilize the bottom layer signal processing module such as bottom FPGA, bottom DSP data are demodulated
Rebuild, thus obtain effective brain function blood flow parameter.The operational order of user also can be connect to high-speed transfer by computer transmission
Mouthful, process rear drive light source through the bottom layer signal such as FPGA, dsp chip, it is achieved the real-time control to light source.
The high density near infrared spectrum cerebral function formation method of the time space frequency multiple coupling of the present invention, according to ray detector
Arrangement, passage connect and the parameter such as connecting length carries out passage division and the time of light source and spatial lamination, decrease passage it
Between crosstalk, improve the temporal resolution of system, spatial resolution and signal to noise ratio.And to the light source igniting time, power and
Frequency smart allocation, the spatial resolution of beneficially raising system and signal to noise ratio.Owing to detector is encoded, can improve also
Row detection efficient, has unrivaled application advantage in high density fNIRS system.
Claims (4)
1. the high density near infrared spectrum cerebral function formation method of a time space frequency multiple coupling, it is characterised in that: include as follows
Step,
The first step, detection channels divides；
Connect and connecting length according to ray detector arrangement, detection channels, the detection channels of different connecting lengths is carried out point
Class, is divided into a class by the detection channels of equal length；
Second step, light source time and spatial lamination；
According to the result of detection channels classification, light source and detection channels are carried out burst and sequence in time, and according to space
Distance, spatially carries out light source space array code, makes the light source of different array light in the different time；
3rd step, light source igniting time, power and frequency smart allocation；
In different connecting lengths, light source is according to different power drive, and the time of light source igniting also differs, detection
The connecting length of passage is the longest, and the light source igniting time is the longest, and power is the biggest, and in same timeslice, divides according to the space of light source
Cloth, carries out frequency intelligent modulation so that the light source that distance is the nearest, frequency shift (FS) is the biggest；
4th step, detector encodes；
Divide according to detection channels and the result of light source time and space burst, detector is encoded；
5th step, light source is according to the routine lights encoded, and detector also according to the programmed acquisition signal encoded, becomes
Picture.
The high density near infrared spectrum cerebral function formation method of a kind of time space frequency multiple coupling the most according to claim 1,
It is characterized in that: described detection channels is divided into long connection and short connection two class, the volume of the light source space array of each short connection
Code rule is:
Light source on same diagonal is lighted in same timeslice, it is meant that the light source of first light source space array is at certain for the moment
Between sheet light simultaneously, the light source of second light source space array is lighted in another timeslice simultaneously；
The coding rule of each long light source space array connected is: ensure to be spaced two between the light source all directions lighted not
The light source lighted.
The high density near infrared spectrum cerebral function formation method of a kind of time space frequency multiple coupling the most according to claim 1,
It is characterized in that: the lighting time of light source and power all approximate exponentially relation, i.e. lighting time t=Be with connecting lengthAd, merit
Rate P=CeAd, wherein d is connecting length, and A, B, C are constant；
Further, in same timeslice, make frequency shift (FS) between the nearest light source the biggest, and between the most remote light source
Frequency shift (FS) is the least, and distance about becomes negative exponent relation with frequency shift (FS).
4. the high density near infrared spectrum cerebral function imaging system of a time space frequency multiple coupling, it is characterised in that: include light source,
Detector, signal filtering and modular converter, bottom layer signal processing module, high-speed transmission interface and computer；
The laser signal that detector collects, after signal filtering and modular converter carry out signal conversion, uses FPGA or DSP core
Sheet carries out bottom layer signal process, and uses high-speed transmission interface to carry out data transmission；Then utilize computer PC or directly utilize
Data are demodulated rebuilding by bottom FPGA, bottom DSP bottom layer signal processing module, thus obtain effective brain function blood flow ginseng
Number.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610490279.5A CN106037657B (en) | 2016-06-28 | 2016-06-28 | A kind of high density near infrared spectrum cerebral function imaging method of time space frequency multiple coupling |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610490279.5A CN106037657B (en) | 2016-06-28 | 2016-06-28 | A kind of high density near infrared spectrum cerebral function imaging method of time space frequency multiple coupling |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106037657A true CN106037657A (en) | 2016-10-26 |
| CN106037657B CN106037657B (en) | 2017-11-21 |
Family
ID=57166830
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610490279.5A Active CN106037657B (en) | 2016-06-28 | 2016-06-28 | A kind of high density near infrared spectrum cerebral function imaging method of time space frequency multiple coupling |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106037657B (en) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107252305A (en) * | 2017-06-18 | 2017-10-17 | 天津大学 | The NIRS cerebral function imaging systems of phase photon counting detection pattern are locked based on full parellel |
| CN114723926A (en) * | 2021-01-07 | 2022-07-08 | 株式会社日本显示器 | Biometric authentication device |
| CN114841899A (en) * | 2022-07-01 | 2022-08-02 | 南京智谱科技有限公司 | Method for removing infrared image horizontal stripes by space-time-frequency combined compact coding and infrared equipment |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2006071891A2 (en) * | 2004-12-23 | 2006-07-06 | The General Hospital Corporation | Evaluating central nervous system |
| US20090292210A1 (en) * | 2008-04-17 | 2009-11-26 | Washington University In St. Louis | Task-less optical mapping of brain function using resting state functional connectivity |
| CN103735274A (en) * | 2013-12-25 | 2014-04-23 | 电子科技大学 | Device and method for detecting absolute amount of blood oxygen and blood volume of local brain tissue |
| CN103968862A (en) * | 2013-02-06 | 2014-08-06 | 约翰内斯﹒海德汉博士有限公司 | Optical position-measuring device |
| EP2888994A1 (en) * | 2013-12-26 | 2015-07-01 | Fundació Institut de Ciències Fotòniques | Speckle contrast optical tomography |
| US20150208964A1 (en) * | 2014-01-27 | 2015-07-30 | Covidien Lp | Systems and methods for determining respiration information |
| US9498134B1 (en) * | 2016-03-28 | 2016-11-22 | Cephalogics, LLC | Diffuse optical tomography methods and system for determining optical properties |
-
2016
- 2016-06-28 CN CN201610490279.5A patent/CN106037657B/en active Active
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2006071891A2 (en) * | 2004-12-23 | 2006-07-06 | The General Hospital Corporation | Evaluating central nervous system |
| US20090292210A1 (en) * | 2008-04-17 | 2009-11-26 | Washington University In St. Louis | Task-less optical mapping of brain function using resting state functional connectivity |
| CN103968862A (en) * | 2013-02-06 | 2014-08-06 | 约翰内斯﹒海德汉博士有限公司 | Optical position-measuring device |
| CN103735274A (en) * | 2013-12-25 | 2014-04-23 | 电子科技大学 | Device and method for detecting absolute amount of blood oxygen and blood volume of local brain tissue |
| EP2888994A1 (en) * | 2013-12-26 | 2015-07-01 | Fundació Institut de Ciències Fotòniques | Speckle contrast optical tomography |
| US20150208964A1 (en) * | 2014-01-27 | 2015-07-30 | Covidien Lp | Systems and methods for determining respiration information |
| US9498134B1 (en) * | 2016-03-28 | 2016-11-22 | Cephalogics, LLC | Diffuse optical tomography methods and system for determining optical properties |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107252305A (en) * | 2017-06-18 | 2017-10-17 | 天津大学 | The NIRS cerebral function imaging systems of phase photon counting detection pattern are locked based on full parellel |
| CN107252305B (en) * | 2017-06-18 | 2020-11-10 | 天津大学 | NIRS brain function imaging system based on full parallel phase-locked photon counting detection mode |
| CN114723926A (en) * | 2021-01-07 | 2022-07-08 | 株式会社日本显示器 | Biometric authentication device |
| CN114841899A (en) * | 2022-07-01 | 2022-08-02 | 南京智谱科技有限公司 | Method for removing infrared image horizontal stripes by space-time-frequency combined compact coding and infrared equipment |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106037657B (en) | 2017-11-21 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106037657B (en) | A kind of high density near infrared spectrum cerebral function imaging method of time space frequency multiple coupling | |
| CN101816552B (en) | Imitating CT scanning mode multispectral time domain fluorescence molecular tomography measuring system | |
| CN104107051B (en) | Portable multiple brain region blood oxygen detection and real-time display system | |
| RU2603240C2 (en) | Radiation detector | |
| CN103156620A (en) | Multi-channel parallel near infrared spectrum imaging system | |
| CN107252305B (en) | NIRS brain function imaging system based on full parallel phase-locked photon counting detection mode | |
| CN1380821A (en) | Method and device for localizing deviant region in turbid medium | |
| CN100506150C (en) | Optical detection device and working method for tissue of living body | |
| CN100370278C (en) | Detector element for spatially resolved detection of gamma radiation | |
| RU2014133352A (en) | X-ray examination using fiber scintillation sensors with shift of wavelengths | |
| CN102920434A (en) | Fluorescence-optical combined tomography system and measuring method | |
| CN1870929A (en) | Automated endoscopy device, diagnostic method, and uses | |
| CN102946794A (en) | Device and method for detecting and monitoring ingredients or properties of a measurement medium, in particular of physiological blood values | |
| US20120232402A1 (en) | Functional Near Infrared Spectroscopy Imaging System and Method | |
| CN102551671B (en) | Photon counting-type dynamic diffusion fluorescence tomography method and device | |
| CN110123281A (en) | Actuated in parallel based on locking phase photon counting technique spreads optical tomography imaging device | |
| CN103679152B (en) | The true and false fingerprint verification method of optically-based chart system | |
| CN109222865A (en) | A kind of multi-modality imaging endoscopic system | |
| CN107167832A (en) | A kind of nuclear detector | |
| CN1550211A (en) | Methods and apparatus for radiation detecting and imaging using monolithic detectors | |
| CN109069076A (en) | Tumor hypoxia is continuously monitored using near infrared spectroscopy and the tomographic imaging with photonic mixer equipment | |
| CN106999117A (en) | Pulse oximetry device and the method for operating pulse oximetry device | |
| CN104055489A (en) | Blood vessel imaging device | |
| CN102599888A (en) | Optical topology imaging system and method on basis of digital phase locking detection technology | |
| CN209172253U (en) | A kind of multi-modality imaging endoscopic 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 | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |