CN1792323A - Method and equipment for transcranial cerebral blood flow high-resolution imaging - Google Patents

Method and equipment for transcranial cerebral blood flow high-resolution imaging Download PDF

Info

Publication number
CN1792323A
CN1792323A CNA2005101205758A CN200510120575A CN1792323A CN 1792323 A CN1792323 A CN 1792323A CN A2005101205758 A CNA2005101205758 A CN A2005101205758A CN 200510120575 A CN200510120575 A CN 200510120575A CN 1792323 A CN1792323 A CN 1792323A
Authority
CN
China
Prior art keywords
blood flow
time
pixel
cerebral blood
field picture
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
Application number
CNA2005101205758A
Other languages
Chinese (zh)
Other versions
CN100502771C (en
Inventor
李鹏程
曾绍群
骆清铭
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Huazhong University of Science and Technology
Original Assignee
Huazhong University of Science and Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Huazhong University of Science and Technology filed Critical Huazhong University of Science and Technology
Priority to CNB2005101205758A priority Critical patent/CN100502771C/en
Publication of CN1792323A publication Critical patent/CN1792323A/en
Application granted granted Critical
Publication of CN100502771C publication Critical patent/CN100502771C/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)
  • Investigating Or Analysing Materials By Optical Means (AREA)

Abstract

A high-definition imaging method for the craniocerebral blood flow includes such steps as irradiating the near infrared collimated laser beam onto an object to be tested, continuously taking several frames of reflected laser speckle image at intervals, calculating the time contrast of light intensity variation with time for each pixel, calculating the speed of its cerebral blood flow, and using it as grayscale to configure the 2D speed distribution chart of cerebral blood flow. Its apparatus is composed of laser device, two polarizers, bench, photoelectric imaging system and computer.

Description

A kind of transcranial cerebral blood flow high-resolution imaging method and device thereof
Technical field
The invention belongs to the cerebral blood flow detection technique, be specially a kind of transcranial cerebral blood flow high-resolution imaging method and device thereof, it is particularly suitable for that cortex regional flow under studying physiological and the pathological state distributes and cerebrovascular form etc.
Background technology
Obtain the local two-dimentional blood flow distributed intelligence of high-resolution cortex to the nerve-blood vessel coupling of studying physiological and morbid state hypencephalon and adjusting, cerebral function imaging, medicine Effect Evaluation, and the diagnosis and the pathological study of great disease of brain such as cerebral ischemia, cerebral hemorrhage are significant to cerebrovascular and microcirculation.But traditional cerebral blood flow detection method scarcely possesses imaging capability, promptly there is not spatial resolving power, admittance is traced as volume, rheoencephalogram based on the brain resistance measurement, transcranial doppler, laser-Doppler etc., and the higher micro tv video recording of resolution, scan laser Doppler, the laser speckle imaging space contrasts technology such as analysis and then needs measurand is carried out operation of opening cranium, remove the skull (even cerebral dura mater) of detected part, expose cortex, perhaps with the skull wear down of detected part to pellucidity, thereby be the method for damage, can influence the normal physiological state of measurand, imaging causes adverse effect to cerebral blood flow, and does not possess the clinical expansion potentiality.
Summary of the invention
The object of the present invention is to provide a kind of transcranial cerebral blood flow high-resolution imaging method, this method has solved the problem that existing high-resolution cerebral blood flow formation method need be opened cranium or wear down skull operation to experimental subject, and the high-resolution of having realized the cerebral blood flow Two dimensional Distribution is through the imaging of cranium not damaged; The present invention also provides the device of realizing this method.
The invention provides a kind of transcranial cerebral blood flow high-resolution imaging method, the steps include:
(1) collimated laser beam of wavelength in 700～950nm scope shone on the measurand;
(2) with identical time of exposure and the some frame measurands of frame period time continuous acquisition laser light reflected speckle image;
(3) each two field picture to gathering, it is capable to take out in every two field picture i, and (i j), utilizes N the gray value that is taken out to calculate on this pixel the time dependent time of light intensity and contrasts K the gray value I of j column position place pixel t(i, j), its computing formula is:
K t ( i , j ) = &sigma; t ( i , j ) < I ( i , j ) > = 1 N - 1 &Sigma; n = 1 N [ I n ( i , j ) - < I ( i , j ) > ] 2 < I ( i , j ) >
I wherein n(i j) represents in the n two field picture i capable, and the gray value of j row place pixel, N are the number of image frames of being gathered, σ t(i, j) represent the N two field picture (i j) locates the standard deviation of pixel gray scale,＜I (i, j)〉be the N two field picture (i j) locates the time average of pixel gray scale, and computing formula is:
< I ( i , j ) &GreaterEqual; 1 N &Sigma; n = 1 N I n ( i , j )
(4) the cerebral blood flow velocity V that utilizes the gained laser speckle time to contrast to calculate this pixel place (i, j), formula is as follows:
V ( i , j ) = c K t 2 ( i , j ) , Wherein c is a correction coefficient
(5) set by step all pixels in the traversing graph picture of (2)～(4), the laser speckle time that obtains all pixel correspondences is contrasted value and blood flow rate value;
(6) be gray scale with the blood flow rate value of each pixel correspondence respectively, make up the cerebral blood flow velocity scattergram of two dimension.
Each two field picture time of exposure is 1ms～20ms in the above-mentioned steps (1), and the frame period time is less than 200ms, frame number N 〉=20 of collection.
Carry out dynamic cerebral blood flow distribution monitoring if desired, then comprise step (7), need repeat the step of (2)～(6) to the time point of cerebral blood flow distribution carrying out imaging, obtain difference two-dimensional brain blood flow distributed images constantly at each.
Realize the device of said method, its structure is: laser beam, first linear polarizer and workbench are positioned on the illumination path successively, and first linear polarizer is vertical with the incident laser light beam; Workbench, second polaroid and photo electric imaging system are positioned on the imaging optical path successively, and second linear polarizer is vertical with the photo electric imaging system optical axis direction, and be concentric with photo electric imaging system, and the polarization direction of its polarization direction and first linear polarizer is perpendicular; Computer links to each other with photo electric imaging system by image pick-up card, is used for data acquisition and processing (DAP).
The present invention is based on dynamic laser speckle imaging time domain statistical characteristic analysis, see through skull and directly cortex medium vessels under the skull and intracortical blood flow are carried out the high spatial resolution imaging.Compare with other existing high-resolution cerebral blood flow detection technique, the advantage of transcranial cerebral blood flow high-resolution imaging method provided by the present invention is: need not the operation that cranium or wear down skull are opened in enforcement to measurand, can directly see through blood flow distribution of high-resolution two dimension and cerebrovascular form that complete skull obtains cortex under the skull with no damage, the cortex hemodynamics be changed carry out the monitoring of real-time, dynamic, high time, spatial resolution.The present invention has the advantage of not damaged and high spatial resolution, and need not to scan.Its range of application is can be in order to laboratory animals such as the rat under studying physiological and the pathological state, mice, rabbit, cat, monkeys, and people's cortex regional flow distributes, and the cortex blood flow that neural activity, disease of brain cause changes.The present invention is applicable to the research of cerebral function imaging, neuro physiology, disease of brain pathology and evaluating drug effect.
Description of drawings
Fig. 1 is a transcranial cerebral blood flow high-resolution imaging apparatus structure sketch map.
Fig. 2 is an image acquisition control software flow pattern.
Fig. 3 a is applied to the experimental result of rat through cranium brain blood cortex stream high-resolution imaging with method disclosed by the invention.
Fig. 3 b adopts other existing method (the laser speckle imaging space contrasts analysis) gained experimental result down with the same experiment condition of Fig. 3 a.
The specific embodiment
The reconstruction of cerebral blood flow distributed image need utilize at each need measure the multiframe laser speckle image that the time point of cerebral blood flow is gathered, the laser speckle image sequence of being gathered is carried out laser speckle time domain statistical characteristic analysis, calculate in the laser speckle image the time dependent time statistic of light intensity (being gradation of image) on each pixel, with this statistic reflect this pixel the blood flow rate at corresponding cortex place; So all pixels in the traversing graph picture can obtain high-resolution two-dimensional brain blood flow distributed image.
Particularly, the step of the inventive method is:
(1) collimated laser beam of wavelength in 700～950nm scope shone on the measurand;
(2) with identical time of exposure and frame period time continuous acquisition N frame measurand laser light reflected speckle image.For reaching better technique effect, each two field picture time of exposure should be 1ms～20ms, and the frame period time should be less than 200ms, and the frame number N of collection generally should be more than or equal to 20;
(3) to collection N two field picture, it is capable to take out in every two field picture i, and (i j), utilizes N the gray value that is taken out to calculate on this pixel the time dependent time of light intensity and contrasts K the gray value I of j column position place pixel t(i, j), its computing formula is:
K t ( i , j ) = &sigma; t ( i , j ) < I ( i , j ) > = 1 N - 1 &Sigma; n = 1 N [ I n ( i , j ) - < I ( i , j ) > ] 2 < I ( i , j ) > ,
I wherein n(i j) represents in the n two field picture i capable, and the gray value of j row place pixel, N are the number of image frames of being gathered, σ t(i, j) represent the N two field picture (i j) locates the standard deviation of pixel gray scale,＜I (i, j)〉be the N two field picture (i j) locates the time average of pixel gray scale, and computing formula is:
< I ( i , j ) &GreaterEqual; 1 N &Sigma; n = 1 N I n ( i , j ) .
(4) the cerebral blood flow velocity V that utilizes the gained laser speckle time to contrast to calculate this pixel place (i, j), formula is as follows:
V ( i , j ) = c K t 2 ( i , j ) , Wherein c is a correction coefficient.
(5) set by step all pixels in the traversing graph picture of (2)～(4), the laser speckle time that obtains all pixel correspondences is contrasted value and blood flow rate value.As every frame laser speckle image size is that M is capable, the K row, then obtaining M * K laser speckle time altogether contrasts value, and M * K blood flow rate value.
(6) be gray scale with the blood flow rate value of each pixel correspondence respectively, can make up the cerebral blood flow velocity scattergram of two dimension.
(7) carry out dynamic cerebral blood flow distribution monitoring if desired, the variation that the different cerebral blood flow constantly of i.e. observation distribute, then need repeat the step of (2)～(6) to the time point of cerebral blood flow distribution carrying out imaging, obtain different high-resolution two-dimensional brain blood flow distributed images constantly at each.
The device of realization said method as shown in Figure 1.Wavelength shines the object to be measured that is placed on the workbench 3 in collimation near-infrared laser light beam 1 warp first polaroid 2 backs of 700～950nm scope, form laser speckle in object surface to be measured, this laser speckle through after second linear polarizer 4 in computer 7 under the control of respective image acquisition software by photo electric imaging system 5 imagings, and pass through image pick-up card 6 with gained laser speckle image input computer 7, wherein photo electric imaging system 5 can be the microscope of electrically charged coupled apparatus camera, the charge-coupled device camera of band photographic lens, analog video camera or DV.It should be noted that the direction that first linear polarizer 2 is placed should be vertical with incident laser light beam 1; Second linear polarizer 4 should vertically be placed along photo electric imaging system 5 optical axis directions, and maintenance and photo electric imaging system 5 is concentric, and should be perpendicular along the polarization direction that makes its polarization direction and first linear polarizer 1, the adverse effect that blood flow detection is caused with the direct reflection of removing region surface to be measured.
The workflow of image acquisition control software as shown in Figure 2 in the computer 7.
Implement one:
Experimental subject is the Wistar rat, is fixed on the workbench 3, and be light source with the semiconductor laser of 780nm, the laser beam 1 behind the collimator and extender is through retread rat head after inciding peeling of first linear polarizer 2.Photo electric imaging system 5 imagings that the laser speckle that rat head reflects to form is made of the charge-coupled device camera of being with the macrovision camera lens after second linear polarizer 4, the optical system amplification is made as 0.5 times.Time of exposure 5ms, frame period time 25ms, continuous acquisition 40 frame laser speckle images.Utilize the 40 frame laser speckle images of being gathered, contrast by the method for the invention calculating laser speckle time, and further rebuilding two-dimensional brain blood flow distributed image, gained result such as Fig. 3 a show that the position of gray value low more (dark more) means that this place's cerebral blood flow is fast more among the figure.The rat skull that sees through that this method can be successful as can be seen from Fig. 3 a obtains the Two dimensional Distribution of cortex blood flow, and can clearly tell the cerebrovascular form of cortex, and resolution is better than 50 microns.For comparing with other method, under identical experiment condition, utilize the laser speckle space to contrast and analyze and carry out blood flow detection, the gained result is shown in Fig. 3 b.From the comparison of Fig. 3 a and Fig. 3 b, the solution of finding out method success provided by the present invention that can be perfectly clear other method can not see through skull and obtain and carry out the problem that the high-resolution cerebral blood flow is scattered in picture under the skull.

Claims (4)

1, a kind of transcranial cerebral blood flow high-resolution imaging method the steps include:
(1) collimated laser beam of wavelength in 700～950nm scope shone on the measurand;
(2) with identical time of exposure and the some frame measurands of frame period time continuous acquisition laser light reflected speckle image;
(3) each two field picture to gathering, it is capable to take out in every two field picture i, and (i j), utilizes N the gray value that is taken out to calculate on this pixel the time dependent time of light intensity and contrasts K the gray value I of j column position place pixel t(i, j), its computing formula is:
K t ( i , j ) = &sigma; t ( i , j ) < I ( i , j ) > = 1 N - 1 &Sigma; n = 1 N [ I n ( i , j ) - < I ( i , j ) > ] 2 < I ( i , j ) > ,
I wherein n(i j) represents in the n two field picture i capable, and the gray value of j row place pixel, N are the number of image frames of being gathered, σ t(i, j) represent the N two field picture (i j) locates the standard deviation of pixel gray scale,＜I (i, j)〉be the N two field picture (i j) locates the time average of pixel gray scale, and computing formula is:
< I ( i , j ) > = 1 N &Sigma; n = 1 N I n ( i , j )
(4) the cerebral blood flow velocity V that utilizes the gained laser speckle time to contrast to calculate this pixel place (i, j), formula is as follows:
V ( i , j ) = c K t 2 ( i , j ) , Wherein c is a correction coefficient,
(5) set by step all pixels in the traversing graph picture of (2)～(4), the laser speckle time that obtains all pixel correspondences is contrasted value and blood flow rate value;
(6) be gray scale with the blood flow rate value of each pixel correspondence respectively, make up the cerebral blood flow velocity scattergram of two dimension.
2, method according to claim 1 is characterized in that: each two field picture time of exposure is 1ms～20ms in the step (1), and the frame period time is less than 200ms, frame number N 〉=20 of collection.
3, method according to claim 1 and 2 is characterized in that: it also comprises step (7), need repeat the step of (2)～(6) to the time point of cerebral blood flow distribution carrying out imaging at each, obtains different two-dimensional brain blood flow distributed images constantly.
4, a kind of transcranial cerebral blood flow high-resolution imaging device, its structure is: laser beam (1), first linear polarizer (2) are positioned on the illumination path successively with workbench (3), and first linear polarizer (2) is vertical with incident laser light beam (1); Workbench (3), second polaroid (4) and photo electric imaging system (5) are positioned on the imaging optical path successively, second linear polarizer (4) is vertical with photo electric imaging system (5) optical axis direction, concentric with photo electric imaging system (5), and its polarization direction is vertical with the polarization direction of first linear polarizer (2); Computer (7) links to each other with photo electric imaging system (5) by image pick-up card (6), is used for data acquisition and processing (DAP).
CNB2005101205758A 2005-12-30 2005-12-30 Method and equipment for transcranial cerebral blood flow high-resolution imaging Expired - Fee Related CN100502771C (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CNB2005101205758A CN100502771C (en) 2005-12-30 2005-12-30 Method and equipment for transcranial cerebral blood flow high-resolution imaging

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CNB2005101205758A CN100502771C (en) 2005-12-30 2005-12-30 Method and equipment for transcranial cerebral blood flow high-resolution imaging

Publications (2)

Publication Number Publication Date
CN1792323A true CN1792323A (en) 2006-06-28
CN100502771C CN100502771C (en) 2009-06-24

Family

ID=36804103

Family Applications (1)

Application Number Title Priority Date Filing Date
CNB2005101205758A Expired - Fee Related CN100502771C (en) 2005-12-30 2005-12-30 Method and equipment for transcranial cerebral blood flow high-resolution imaging

Country Status (1)

Country Link
CN (1) CN100502771C (en)

Cited By (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102357033A (en) * 2011-09-27 2012-02-22 华中科技大学 Laser speckle blood stream imaging processing system and method
CN103393404A (en) * 2013-07-26 2013-11-20 重庆蝶波科技有限公司 Multi-wavelength infrared polarized light real-time image guide system
CN104873192A (en) * 2015-06-10 2015-09-02 上海大学 Miniature brain function monitoring and imaging apparatus and method
CN105395184A (en) * 2015-12-04 2016-03-16 华中科技大学 Biological tissue blood flow, blood oxygen and blood volume multi-parameter detection method and device
CN105433906A (en) * 2015-12-14 2016-03-30 华中科技大学 Scanning dark field laser speckle blood flow imaging method and device
CN106419890A (en) * 2016-11-14 2017-02-22 佛山科学技术学院 Blood speed measuring device and method based on space-time modulation
CN107260155A (en) * 2016-04-07 2017-10-20 动心医电股份有限公司 Detect the method and its device of velocity of blood flow
CN107320112A (en) * 2017-07-07 2017-11-07 广州医软智能科技有限公司 The multi-parameter imaging detection method and device of microcirculation
CN110507305A (en) * 2019-08-27 2019-11-29 北京大学 Contrast the measuring blood flow rate method of waveform conduction time difference based on laser speckle
CN112617789A (en) * 2020-07-28 2021-04-09 上海大学 Laser speckle blood flow imaging method and system
CN113197564A (en) * 2021-04-27 2021-08-03 燕山大学 Portable neurovascular coupling detection device for conscious animals
CN116718514A (en) * 2023-04-26 2023-09-08 广州医科大学附属第一医院（广州呼吸中心） Simulation detection device for depth of water laser boosting atomization steam inlet nose lung

Cited By (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102357033A (en) * 2011-09-27 2012-02-22 华中科技大学 Laser speckle blood stream imaging processing system and method
CN103393404A (en) * 2013-07-26 2013-11-20 重庆蝶波科技有限公司 Multi-wavelength infrared polarized light real-time image guide system
CN104873192A (en) * 2015-06-10 2015-09-02 上海大学 Miniature brain function monitoring and imaging apparatus and method
CN105395184A (en) * 2015-12-04 2016-03-16 华中科技大学 Biological tissue blood flow, blood oxygen and blood volume multi-parameter detection method and device
CN105395184B (en) * 2015-12-04 2018-05-08 华中科技大学 The multi-parameter detecting method and device of biological tissue's blood flow, blood oxygen and blood volume
CN105433906A (en) * 2015-12-14 2016-03-30 华中科技大学 Scanning dark field laser speckle blood flow imaging method and device
CN107260155A (en) * 2016-04-07 2017-10-20 动心医电股份有限公司 Detect the method and its device of velocity of blood flow
CN106419890A (en) * 2016-11-14 2017-02-22 佛山科学技术学院 Blood speed measuring device and method based on space-time modulation
CN106419890B (en) * 2016-11-14 2024-04-30 佛山科学技术学院 Blood flow velocity measuring device and method based on space-time modulation
CN107320112A (en) * 2017-07-07 2017-11-07 广州医软智能科技有限公司 The multi-parameter imaging detection method and device of microcirculation
CN110507305A (en) * 2019-08-27 2019-11-29 北京大学 Contrast the measuring blood flow rate method of waveform conduction time difference based on laser speckle
CN110507305B (en) * 2019-08-27 2024-04-16 北京大学 Blood flow velocity measurement method based on laser speckle contrast waveform conduction time difference
CN112617789A (en) * 2020-07-28 2021-04-09 上海大学 Laser speckle blood flow imaging method and system
CN113197564A (en) * 2021-04-27 2021-08-03 燕山大学 Portable neurovascular coupling detection device for conscious animals
CN116718514A (en) * 2023-04-26 2023-09-08 广州医科大学附属第一医院（广州呼吸中心） Simulation detection device for depth of water laser boosting atomization steam inlet nose lung

Also Published As

Publication number Publication date
CN100502771C (en) 2009-06-24

Similar Documents

Publication Publication Date Title
CN100502771C (en) Method and equipment for transcranial cerebral blood flow high-resolution imaging
CN101485565B (en) Laser speckle blood current imaging and analyzing method
Bonin et al. Local diversity and fine-scale organization of receptive fields in mouse visual cortex
CN1208025C (en) Method and device for evaluating an kinetic pattern
CN1122682A (en) Blood analysing device which do not hurt body
WO2008149055A1 (en) Automatic characterization of cellular motion
Kato et al. A quantification of goldfish behavior by an image processing system
CN109752377B (en) Spectroscopic bimodal projection tomography tissue blood vessel imaging device and method
Chemla et al. Improving voltage-sensitive dye imaging: with a little help from computational approaches
Zhu et al. A deep learning framework for inference of single-trial neural population dynamics from calcium imaging with subframe temporal resolution
US20170350874A1 (en) Optical interrogation and control of dynamic biological functions
CN1975378A (en) Gray mold early diagnosing system based on visable/near infrared multispectral imaging
CN105342597B (en) A kind of quantitative laser blood flow detection method
CN107485383A (en) A kind of speckle blood flow imaging method and apparatus based on constituent analysis
Senseman et al. High-speed VSD imaging of visually evoked cortical waves: decomposition into intra-and intercortical wave motions
Nöbauer et al. Mesoscale volumetric light-field (MesoLF) imaging of neuroactivity across cortical areas at 18 Hz
Hu et al. Deep learning based wavefront sensor for complex wavefront detection in adaptive optical microscopes
CN113298700B (en) High-resolution image reconstruction method in scattering scene
CN117974525A (en) Laser speckle contrast blood flow imaging method based on second-order autocorrelation function calculation
Cunha et al. Movement quantification during epileptic seizures: a new technical contribution to the evaluation of seizure semiology
Sokol et al. The Use of Digital Interferometry Devices to Analyze the State of Red Blood Cell Membranes
Menzel et al. Using light and X-ray scattering to untangle complex neuronal orientations and validate diffusion MRI
Kalitzin et al. Optical flow group-parameter reconstruction from multi-channel image sequences
CN209624404U (en) A kind of beam splitting type bimodal projection chromatography tissue blood vessel imaging device
Flotho et al. Semi–synthetic dataset for the evaluation of motion compensation approaches for voltage sensitive dye imaging

Legal Events

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

Granted publication date: 20090624

Termination date: 20201230