CN110464329A - A kind of method and apparatus of velocity of blood flow distribution measuring - Google Patents
A kind of method and apparatus of velocity of blood flow distribution measuring Download PDFInfo
- Publication number
- CN110464329A CN110464329A CN201910657334.9A CN201910657334A CN110464329A CN 110464329 A CN110464329 A CN 110464329A CN 201910657334 A CN201910657334 A CN 201910657334A CN 110464329 A CN110464329 A CN 110464329A
- Authority
- CN
- China
- Prior art keywords
- reconstruction region
- unit
- electrode
- region
- several
- 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
Classifications
-
- 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/0265—Measuring blood flow using electromagnetic means, e.g. electromagnetic flowmeter
-
- 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/0285—Measuring or recording phase velocity of blood waves
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Hematology (AREA)
- Heart & Thoracic Surgery (AREA)
- Surgery (AREA)
- Biophysics (AREA)
- Pathology (AREA)
- Engineering & Computer Science (AREA)
- Biomedical Technology (AREA)
- Cardiology (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Physiology (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Electromagnetism (AREA)
- Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
The invention discloses a kind of method and apparatus of velocity of blood flow distribution measuring, wherein a kind of method of velocity of blood flow distribution measuring, comprising: in the magnetic field of the cross-section setting constant homogeneous to reconstruction region;By n electrode arrangements the cross section outside;Using one of electrode as reference electrode, remaining n-1 electrode is measuring electrode, obtains n-1 induced electricity potential difference;By the neighbouring electrode of the reference electrode clockwise or counter-clockwise, is chosen as new reference electrode, remaining n-1 electrode is measuring electrode, obtains n-1 induced electricity potential difference;Until traversing all electrodes is new reference electrode, n* (n-1) a induced electricity potential difference is finally obtained;A induced electricity potential difference of n* (n-1) calculates the blood flowrate of each reconstruction region, and then obtain the velocity of blood flow to reconstruction region and be distributed by reconstructing equation as input data.To solve the problem of that traditional measurement obtains that data volume is few to influence measurement accuracy.
Description
Technical field
The present invention relates to biomedical engineering field, specifically a kind of method and apparatus of velocity of blood flow distribution measuring.
Background technique
According to Faraday's electromagnetic induction law, when conductor passes through stationary magnetic field with a certain speed, both ends can be generated vertically
In the induced electricity potential difference of magnetic direction, induced electricity potential difference size has linear relationship with the product of magnetic field strength and conductor speed,
Here it is the basic principles based on electromagnetic induction measurement fluid velocity.
Blood flow measurement technology based on electromagnetic induction can be applied to the narrow early diagnosis and prevention of human body artery, and existing
It is few by induced potential difference data based on electromagnetic induction blood flow measurement technology and increase measuring electrode and the cost that causes increases
Contradictory limitation, the blood flow or flow velocity of regional area can only be measured at present.Although can judge the stenosis of artery,
It can not determine specific narrow location.
Application No. is 201810927376.5 patents to disclose a kind of blood of human body flow-speed measurement method and its device, bright
By the setting of multiple electrodes and multiple measured zones, the limitation of traditional multi-electrode electromagnetic flow meter convection is broken, will pass
The average speed for the entire measurement cross section that system flowmeter is had in mind is converted into the average axial velocity in each infinitesimal, solves doctor
Velocity of blood flow test problems in treatment.But the dependence there are electromagnetic induction measurement accuracy and reconstruction region to electrode number,
The artificial subjective problem of reconstruction region division mode.Such as: the above method and its device can only obtain 15 effectively
Potential difference data, at most can only 15 reconstruction regions of accurate reconstruction, and artificially selected 15 reconstruction regions, 15 reconstruct
Region includes there are the region of blood vessel and there is no the regions of blood vessel, artificially select it is subjective, and during actual measurement
The specific location of blood vessel is not known, it is difficult to be applied in practice.
Summary of the invention
In view of this, the present invention provides a kind of method and apparatus of velocity of blood flow distribution measuring, to solve traditional measurement
The problem of acquisition data volume is few, influences measurement accuracy.
In a first aspect, the present invention provides a kind of method of velocity of blood flow distribution measuring, comprising:
In the magnetic field of the cross-section setting constant homogeneous to reconstruction region;
By n electrode arrangements the cross section outside;
Using one of electrode as reference electrode, remaining n-1 electrode is measuring electrode, obtains n-1 induced electricity potential difference;
By clockwise or counter-clockwise, choosing the neighbouring electrode of the reference electrode as new reference electrode, remaining n-1
A electrode is measuring electrode, obtains n-1 induced electricity potential difference;
Until traversing all electrodes is new reference electrode, n* (n-1) a induced electricity potential difference is finally obtained;
A induced electricity potential difference of n* (n-1) calculates the blood of each reconstruction region by reconstructing equation as input data
Flow velocity value, and then obtain the velocity of blood flow to reconstruction region and be distributed;
Wherein, described to reconstruction region is region to be measured, n > 2.
Preferably, it before the blood flowrate for calculating each reconstruction region, divides, obtains to reconstruction region to described
To several regions, several regions are reconstruction region, to calculate the velocity of blood flow in each reconstruction region to the n* (n-1)
The contribution weight of a induced electricity potential difference.
Preferably, the division methods of the reconstruction region are as follows:
Selection is described to reconstruction region, and the detection shape to reconstruction region;
Find it is described to the geometric center of reconstruction region at;
To establish two-dimentional polar coordinate system as pole at the geometric center;
In the two-dimentional polar coordinate system, as pole, to be carried out to the shape by several scales at the geometric center
It reduces, obtains several scale shapes;
Using pole as initial point, using the polar axis as side, rotates, obtain in the two-dimentional polar coordinate system according to set angle
Several polar diameters;
Several scale shapes are divided into several regions respectively by several polar diameters, and several regions are reconstruct area
Domain;
Wherein, it is described to reconstruction region be need to measure velocity of blood flow distribution region.
Preferably, the cross section to reconstruction region is equivalent to justify, the shape to reconstruction region is circle;Or
The cross section to reconstruction region is equivalent to justify, the shape to reconstruction region is circle;Described two
It ties up in polar coordinate system, as pole, to be reduced to the shape by several scales at the geometric center, obtain several scales
Shape method particularly includes: the number h for determining several scale shapes, using the pole as the center of circle, using k*R/h as radius according to
It is secondary to make concentric circles outward, obtain several scale shapes;
Wherein, several scale shapes are circular shape;Radius of the R for described in reconstruction region, coefficient k=1,
2...h。
Preferably, the reconstruct equation are as follows:
Wherein, VMIndicate the average velocity of axes of m-th reconstruction region, WNMIndicate the flowing in m-th region to n-th sense
Answer the contribution weight of potential difference, AMIndicate the area of m-th reconstruction region, Δ UNIndicate that n-th measures induced electricity potential difference.
Second aspect, the present invention provide a kind of device of velocity of blood flow distribution measuring, comprising:
Magnetic field units, in the magnetic field of the cross-section setting constant homogeneous to reconstruction region;
Electrode unit, by n electrode arrangements the cross section outside;
Induced electricity potential difference acquiring unit, for using one of electrode as reference electrode, remaining n-1 electrode to be measurement electricity
Pole obtains n-1 induced electricity potential difference;
Updating unit, it is electric as new reference by the neighbouring electrode for clockwise or counter-clockwise, choosing the reference electrode
Pole, remaining n-1 electrode are measuring electrode, obtain n-1 induced electricity potential difference;
Traversal Unit finally obtains n* (n-1) a induced electricity potential difference until traversing all electrodes is new reference electrode;
Computing unit, a induced electricity potential difference of n* (n-1) are calculated each heavy as input data by reconstructing equation
The blood flowrate in structure region, and then obtain the velocity of blood flow to reconstruction region and be distributed;
Wherein, described to reconstruction region is region to be measured, n > 2.
Preferably, the device, further includes:
Reconstruction region division unit or device;
The reconstruction region division unit or device, before the blood flowrate for calculating each reconstruction region, to institute
It states and is divided to reconstruction region, obtain several regions, several regions are reconstruction region, to calculate in each reconstruction region
Velocity of blood flow to the contribution weight of a induced electricity potential difference of the n* (n-1).
Preferably, the reconstruction region division unit or device, comprising:
Detection unit, computing unit, polar coordinate system establish unit, scale shape unit polar diameter unit and division unit;Institute
It states detection unit to connect with the computing unit and the scale shape unit respectively, the computing unit and the polar coordinate system
Unit connection is established, the polar coordinate system is established unit and also connect with the scale shape unit and the polar diameter unit, described
Scale shape unit and the polar diameter unit are also connect with the division unit;
The detection unit, for selecting to reconstruction region, and the detection shape to reconstruction region;
The computing unit, for find it is described to the geometric center of reconstruction region at;
The polar coordinate system establishes unit, for establish two-dimentional polar coordinate system as pole at the geometric center;
The scale shape unit, in the two-dimentional polar coordinate system, using at the geometric center as pole, to institute
It states shape to be reduced by several scales, obtains several scale shapes;
The polar diameter unit is used for using pole as initial point, using the polar axis as side, according to set angle in the two-dimentional pole
It is rotated in coordinate system, obtains several polar diameters;
Several scale shapes are divided into several regions, institute respectively for several polar diameters by the division unit
Stating several regions is reconstruction region;
Wherein, it is described to reconstruction region be need to measure velocity of blood flow distribution region.
Preferably, the detection unit, comprising: equivalent unit;
The equivalent unit, for the cross section to reconstruction region to be equivalent to justify, the shape to reconstruction region
Shape is circle;Or
The equivalent unit, for the cross section to reconstruction region to be equivalent to justify, the scale shape unit is complete
At following operation: determine the number h of several scale shapes, it is successively outside by radius of k*R/h using the pole as the center of circle
Make concentric circles, obtains several scale shapes;
Wherein, several scale shapes are circular shape;Radius of the R for described in reconstruction region, coefficient k=1,
2...h。
Preferably, the reconstruct equation are as follows:
Wherein, VMIndicate the average velocity of axes of m-th reconstruction region, WNMIndicate the flowing in m-th region to n-th sense
Answer the contribution weight of potential difference, AMIndicate the area of m-th reconstruction region, Δ UNIndicate that n-th measures induced electricity potential difference.
The present invention at least has the following beneficial effects:
The present invention provides a kind of method and apparatus of velocity of blood flow distribution measuring, obtains data volume to solve traditional measurement
Less, the problem of affecting measurement accuracy.
Detailed description of the invention
By the description below with reference to attached drawing to the embodiment of the present invention, the above and other purposes of the present invention, feature and
Advantage is apparent, in the accompanying drawings:
Fig. 1 is the method flow diagram that a kind of reconstruction region of the embodiment of the present invention divides;
Fig. 2 is the schematic device that a kind of reconstruction region of the embodiment of the present invention divides;
Fig. 3 is that the reconstruction region that the present invention is equivalent to bowlder to the cross section of reconstruction region shape divides schematic diagram;
Fig. 4 is a kind of method flow diagram of velocity of blood flow distribution measuring of the embodiment of the present invention;
Fig. 5 is upper limb velocity of blood flow distribution reconstructing method schematic illustration of the embodiment of the present invention based on electromagnetic induction;
Specific embodiment
Below based on embodiment, present invention is described, but it is worth noting that, the present invention is not limited to these realities
Apply example.It is detailed to describe some specific detail sections below in datail description of the invention.However, for not detailed
The present invention can also be understood completely in the part described to the greatest extent, those skilled in the art.
In addition, it should be understood by one skilled in the art that provided attached drawing simply to illustrate that the purpose of the present invention,
Feature and advantage, attached drawing are not to be actually drawn to scale.
Meanwhile unless the context clearly requires otherwise, "include", "comprise" etc. otherwise throughout the specification and claims
Similar word should be construed as the meaning for including rather than exclusive or exhaustive meaning;That is, being " including but not limited to "
Meaning.
Fig. 1 is the method flow diagram that a kind of reconstruction region of the embodiment of the present invention divides.As shown in Figure 1, a kind of reconstruction region
The method of division, comprising: step 101 selection is to reconstruction region, and the detection shape to reconstruction region;Step 102 is sought
Look for it is described to the geometric center of reconstruction region at;Step 103 is to establish two-dimentional polar coordinate system as pole at the geometric center;Step
Rapid 104 in the two-dimentional polar coordinate system, as pole, to contract to the shape by several scales at the geometric center
It is small, obtain several scale shapes;Step 105 is existed as side according to set angle using pole as initial point using ray where the polar axis
It is rotated in the two dimension polar coordinate system, obtains several polar diameters;Several polar diameters described in step 106 are respectively by several scale shapes
Several regions are divided into, several regions are reconstruction region;Wherein, described to divide to reconstruction region to need to measure velocity of blood flow
The region of cloth.
Step 101 selection is to reconstruction region, and detects the shape to reconstruction region method particularly includes: as selected
Human upper limb (arm) is to reconstruction region, that is, needing to measure the region that velocity of blood flow is distributed is human upper limb, is examined using edge
Survey method detects the shape of human upper limb, detects and draws the edge lines that human upper limb waits for reconstructed sectional, obtain area to be reconstructed
Domain.
Step 102 find it is described to the geometric center of reconstruction region at method particularly includes: will be equivalent to advise to reconstruction region
Then figure such as selects human upper limb (arm) for reconstruction region, then equivalent can be round or oval its.If will be to
Reconstruction region is equivalent to circle, then to the geometric center of reconstruction region in the center point;If ellipse will be equivalent to reconstruction region
Shape, then to the geometric center of reconstruction region at the center of oval two focus lines.
Step 103 is to establish two-dimentional polar coordinate system as pole at the geometric center method particularly includes: such as with area to be reconstructed
It is pole O at the geometric center in domain, is polar axis Ox from pole to a ray is drawn to reconstruction region outside, then selected one long
The positive direction (counterclockwise or clockwise) for spending unit and angle, for any point P in plane, indicates line segment with r
The length of OP, θ indicate the angle from Ox to OP, and r is called the polar diameter of point P, and θ is called the polar angle of point P, has ordinal number to be just (r, θ)
The polar coordinates of point M, the coordinate system being built such that are called polar coordinate system, since the method for building up of two-dimentional polar coordinate system is fairly simple,
Legend is not provided herein, specifically can refer to Fig. 3 understanding.
Step 104 is in the two-dimentional polar coordinate system, as pole, to press several rulers to the shape at the geometric center
Degree is reduced, and several scale shapes are obtained method particularly includes: centered on the pole at the geometric center, in detection institute
It is done if stating and being chosen in the shape edges to reconstruction region, if several distances for doing the center are calculated, according to institute
The certain proportion coefficient (that is, several scales) for stating several distances is reduced to obtain several rulers to the shape to reconstruction region
Spend shape, wherein proportionality coefficient is less than 1.
Step 105 is using pole as initial point, using ray where the polar axis as side, sits according to set angle in the two-dimentional pole
It is rotated in mark system, obtains several polar diameters method particularly includes: using pole as initial point, using ray where the polar axis as side, according to
Set angle rotates in the two-dimentional polar coordinate system, can rotate, can also rotate clockwise counterclockwise.To rotate clockwise
For: when using pole as initial point, ray where the polar axis is side, is rotated in the two-dimentional polar coordinate system according to set angle
Afterwards, with the first polar diameter is obtained after reconstruction region edge intersects, then using pole as initial point, with ray where first polar diameter
The second polar diameter is obtained after rotating in the two-dimentional polar coordinate system according to set angle for side, and so on, obtain the pole N
Diameter.That is to say, being, the number of polar diameter is 360 °/set angle, if set angle is 20 °, the number of polar diameter is 18.
Several scale shapes are divided into several regions, several regions respectively by several polar diameters described in step 106
For reconstruction region;Wherein, it is described to reconstruction region be need to measure velocity of blood flow distribution region method are as follows: as polar diameter
Number is 3, and several scale shapes are 1, then obtains 3* (1+1)=6 reconstruction region.
It is done if being chosen in the detection shape edges to reconstruction region, if doing the center described in calculating
Several distances, according to several distances certain proportion coefficient (that is, several scales) to the shape to reconstruction region into
Row, which reduces, obtains several scale shapes, wherein proportionality coefficient is less than 1.
In a step 101, the cross section to reconstruction region is equivalent to justify, the shape to reconstruction region is circle
Shape;I.e. the cross section to reconstruction region is equivalent to justify by selection after reconstruction region, and detection is described to reconstruction region at this time
Shape is circle;Or be equivalent to justify by the cross section to reconstruction region, the shape to reconstruction region is circle;Institute
It states in two-dimentional polar coordinate system, as pole, to reduce, obtain several by several scales to the shape at the geometric center
Scale shape method particularly includes: the number h for determining several scale shapes is half with k*R/h using the pole as the center of circle
Diameter successively makees concentric circles outward, obtains several scale shapes;Wherein, several scale shapes are circular shape;R is institute
State the radius to reconstruction region, coefficient k=1,2...h.
Meanwhile before selection is to reconstruction region, the number M that setting reconstruction region divides is drawn according to the reconstruction region
The number M divided, determines the number h and 2 π of the set angle/m of several scale shapes;Wherein, m=M/ (h+1).
Fig. 2 is the schematic diagram for the device that a kind of reconstruction region of the embodiment of the present invention divides.As shown in Fig. 2, a kind of reconstruct area
The device that domain divides, comprising: detection unit 201, computing unit 202, polar coordinate system establish unit 203, scale shape unit 204
Polar diameter unit 205 and division unit 206;The detection unit 201 respectively with the computing unit 202 and the scale shape list
Member 204 connects, and the computing unit 202 is established unit 203 with the polar coordinate system and connect, and the polar coordinate system establishes unit
203 also connect with the scale shape unit 204 and the polar diameter unit 205, the scale shape unit 204 and the polar diameter
Unit 205 is also connect with the division unit 206;The detection unit 201, for selecting to reconstruction region, and detection institute
State the shape to reconstruction region;The computing unit 202, for find it is described to the geometric center of reconstruction region at;The pole
Establishment of coordinate system unit 203, for establish two-dimentional polar coordinate system as pole at the geometric center;The scale shape unit
204, it is used in the two-dimentional polar coordinate system, as pole, to be carried out to the shape by several scales at the geometric center
It reduces, obtains several scale shapes;The polar diameter unit 205, for being with ray where the polar axis using pole as initial point
Side rotates in the two-dimentional polar coordinate system according to set angle, obtains several polar diameters;The division unit 206, for described
Several scale shapes are divided into several regions respectively by several polar diameters, and several regions are reconstruction region;Wherein, described
It is the region for needing to measure velocity of blood flow distribution to reconstruction region.Concrete implementation mode can refer to the detailed description in Fig. 1.
In Fig. 2, the detection unit 201, comprising: equivalent unit;The equivalent unit, for by the area to be reconstructed
The cross section in domain is equivalent to justify, and the shape to reconstruction region is circle;Or the equivalent unit, for by described wait reconstruct
The cross section in region is equivalent to justify, and the scale shape unit 204 completes following operation: determining of several scale shapes
Number h successively makees concentric circles by radius of k*R/h using the pole as the center of circle outward, obtains several scale shapes;Wherein,
Several scale shapes are circular shape;R is the radius to reconstruction region, coefficient k=1,2...h.Concrete implementation
Mode can refer to the detailed description in Fig. 1.
In Fig. 2, a kind of device that reconstruction region divides, further includes: determination unit;The determination unit and the calculating
Unit 202 connects, and before selection is to reconstruction region, the determination unit is used to set or determine that reconstruction region divides
Number M determines the number h and the set angle 2 of several scale shapes according to the number M that the reconstruction region divides
π/m;Wherein, m=M/ (h+1).Concrete implementation mode can refer to the detailed description in Fig. 1.
Fig. 3 is that the reconstruction region that the present invention is equivalent to bowlder to the cross section of reconstruction region shape divides schematic diagram.Such as Fig. 3
It is shown, human upper limb cross section is such as divided into M reconstruction region, by taking Fig. 3 division mode as an example, will be divided into reconstruction region
64 reconstruction regions, wherein the number h=3 and the set angle of several scale shapes are 2 π/16.
In Fig. 3, before the blood flowrate for calculating each reconstruction region, divided to described to reconstruction region,
Several regions are obtained, several regions are reconstruction region, to calculate the velocity of blood flow in each reconstruction region to the n* (n-
1) the contribution weight of a induced electricity potential difference.If the quantity M of reconstruction region is 64, the label mode of each reconstruction region are as follows: polar axis,
First polar diameter and smallest dimension shape intersecting area are the first reconstruction region 1, polar axis, the first polar diameter and neighbouring next scale shape
Intersecting area is the second reconstruction region 2, third reconstructed region 3, fourth reconstructed region 4 is followed successively by from inside to outside, later according to inverse
Clocking sequence is the 5th reconstruction region 5 by the first polar diameter, the second polar diameter and smallest dimension shape intersecting area, according to above-mentioned phase
The 6th reconstruction region 6, the 7th reconstruction region 7 and the 8th reconstruction region 8 are followed successively by from inside to outside with method, according still further to the above method
Sequence counter-clockwise is arranged successively label 9,10....64 (not providing on figure) and is divided into the 9th reconstruction region and reconstructs area for the tenth
The 64th reconstruction region of domain ....Blood in 64 reconstruction regions, such as the quantity n=16 of electrode, in each reconstruction region
Contribution weight of the flow velocity to a induced electricity potential difference of the 16* (16-1).
The present invention also proposes the device that another reconstruction region divides, comprising: memory and processor and is stored in storage
On device and the computer program that can run on a processor, the computer program is such as above-mentioned method, and the processor is held
Selection is performed the steps of when row described program to reconstruction region, and the detection shape to reconstruction region;Described in searching
At the geometric center of reconstruction region;To establish two-dimentional polar coordinate system as pole at the geometric center;It is sat in the two-dimentional pole
In mark system, as pole, to be reduced to the shape by several scales at the geometric center, obtain several scale shapes;
Using pole as initial point, using ray where the polar axis as side, rotates, obtain in the two-dimentional polar coordinate system according to set angle
Several polar diameters;Several scale shapes are divided into several regions respectively by several polar diameters, and several regions are reconstruct
Region;Wherein, it is described to reconstruction region be need to measure velocity of blood flow distribution region.Concrete implementation mode can refer to Fig. 1
In detailed description.
Fig. 4 is a kind of method flow diagram of velocity of blood flow distribution measuring of the embodiment of the present invention.As shown in figure 4, a kind of blood
The method of velocity flow profile measurement, comprising: step 401 is in the magnetic field of the cross-section setting constant homogeneous to reconstruction region;Step
402 by n electrode arrangements the cross section outside;Step 403 is using one of electrode as reference electrode, and remaining n-1
Electrode is measuring electrode, obtains n-1 induced electricity potential difference;Step 404 chooses the reference electrode by clockwise or counter-clockwise
Neighbouring electrode as new reference electrode, remaining n-1 electrode is measuring electrode, obtains n-1 induced electricity potential difference;Step
405, up to traversing all electrodes as new reference electrode, finally obtain n* (n-1) a induced electricity potential difference;N* described in step 406
(n-1) a induced electricity potential difference calculates the blood flowrate of each reconstruction region, in turn by reconstructing equation as input data
The velocity of blood flow to reconstruction region is obtained to be distributed;Wherein, described to reconstruction region is region to be measured, n > 2.
In Fig. 4, be illustrated using human upper limb (arm) as to reconstruction region: step 401 is in area to be reconstructed
The magnetic field of the cross-section setting constant homogeneous in domain, i.e., be arranged a constant homogeneous in human upper limb cross-section using excitation unit
Magnetic field.For step 402 by n electrode arrangements in the outside of the cross section, the outside of the cross section is to reconstruction region
On the outside of skin, such as: when reconstruction region is human upper limb, n electrode arrangements are in the outside of human upper limb skin, i.e. skin of upper extremity
Surface layer;For step 403 using one of electrode as reference electrode, remaining n-1 electrode is measuring electrode, obtains n-1 induced electricity
Potential difference, that is, using one of electrode as reference electrode, remaining electrode is measuring electrode, obtains n-1 actual induction potential difference;
Step 404 chooses the neighbouring electrode of the reference electrode as new reference electrode, remaining n-1 by clockwise or counter-clockwise
A electrode is measuring electrode, obtains n-1 induced electricity potential difference, and such as counterclockwise, choosing next electrode is reference electrode,
Remaining all electrode is measuring electrode, then obtains n-1 actual induction potential difference;Step 405 is until it is new for traversing all electrodes
Reference electrode, finally obtain n* (n-1) a induced electricity potential difference, multiple step 404, until traversing all electrodes is reference electrode,
It is final to obtain n* (n-1) a actual induction potential difference, all actual induction potential difference N=n* (n-1);N* described in step 406
(n-1) a induced electricity potential difference calculates the blood flowrate of each reconstruction region, in turn by reconstructing equation as input data
The velocity of blood flow to reconstruction region is obtained to be distributed, that is, by reconstruct equation, the blood flowrate of each reconstruction region is calculated,
And then obtain the velocity of blood flow distribution in entire upper limb section.Wherein, described to reconstruction region is region to be measured, n > 2.
In Fig. 4, before the magnetic field to the cross-section of reconstruction region setting constant homogeneous, to described to reconstruction region
It is divided, obtains several regions, several regions are reconstruction region.Upper limb cross section described in step 101 is divided into M
A reconstruction region, wherein M≤N;Wherein, all actual induction potential difference N=n* (n-1).
In Fig. 4, the division methods of the reconstruction region are as follows: selection is described described to weight to reconstruction region, and detection
The shape in structure region;Find it is described to the geometric center of reconstruction region at;To establish two-dimentional pole as pole at the geometric center
Coordinate system;In the two-dimentional polar coordinate system, as pole, to contract to the shape by several scales at the geometric center
It is small, obtain several scale shapes;Using pole as initial point, using ray is side where the polar axis, according to set angle described two
It is rotated in dimension polar coordinate system, obtains several polar diameters;Several scale shapes are divided into several areas respectively by several polar diameters
Domain, several regions are reconstruction region;Wherein, it is described to reconstruction region be need to measure velocity of blood flow distribution region.Tool
The implementation of body can refer to the detailed description in Fig. 1 or Fig. 2.
In Fig. 4, the cross section to reconstruction region is equivalent to justify, the shape to reconstruction region is circle;
Or be equivalent to justify by the cross section to reconstruction region, the shape to reconstruction region is circle;It is sat in the two-dimentional pole
In mark system, as pole, to be reduced to the shape by several scales at the geometric center, obtain several scale shapes
Method particularly includes: determine the number h of several scale shapes, it is successively outside by radius of k*R/h using the pole as the center of circle
Make concentric circles, obtains several scale shapes;Wherein, several scale shapes are circular shape;R is the area to be reconstructed
The radius in domain, coefficient k=1,2...h.Concrete implementation mode can refer to the detailed description in Fig. 1 or Fig. 2.
In a kind of method of velocity of blood flow distribution measuring, the reconstruct equation are as follows:
Wherein, VMIndicate the average velocity of axes of m-th reconstruction region, WNMIndicate the flowing in m-th region to n-th sense
Answer the contribution weight of potential difference, AMIndicate the area of m-th reconstruction region, Δ UNIndicate that n-th measures induced electricity potential difference.
Present invention simultaneously provides a kind of devices of velocity of blood flow distribution measuring, comprising: magnetic field units, in area to be reconstructed
The magnetic field of the cross-section setting constant homogeneous in domain;Electrode unit, by n electrode arrangements the cross section outside;Induction
Potential difference acquiring unit is used for using one of electrode as reference electrode, remaining n-1 electrode is measuring electrode, obtains n-1
Induced electricity potential difference;Updating unit chooses the neighbouring electrode of the reference electrode as new reference by clockwise or counter-clockwise
Electrode, remaining n-1 electrode are measuring electrode, obtain n-1 induced electricity potential difference;Traversal Unit, until all electrodes of traversal are
New reference electrode finally obtains n* (n-1) a induced electricity potential difference;Computing unit, a induced electricity potential difference conduct of n* (n-1)
Input data calculates the blood flowrate of each reconstruction region, and then obtain the blood to reconstruction region by reconstructing equation
Flow velocity distribution;Wherein, described to reconstruction region is region to be measured, n > 2.Concrete implementation mode can refer to Fig. 4 and figure
Detailed description in 5.
A kind of device of velocity of blood flow distribution measuring, further includes: reconstruction region division unit or device;The reconstruct area
Domain division unit or device are drawn to described to reconstruction region before the blood flowrate for calculating each reconstruction region
Point, several regions are obtained, several regions are reconstruction region, to calculate the velocity of blood flow in each reconstruction region to the n*
(n-1) the contribution weight of a induced electricity potential difference.
The device of one of Fig. 2 reconstruction region division, packet can be used in the reconstruction region division unit or device
Include: detection unit 201, computing unit 202, polar coordinate system are established unit 203,204 polar diameter unit 205 of scale shape unit and are drawn
Sub-unit 206;The detection unit 201 is connect with the computing unit 202 and the scale shape unit 204 respectively, described
Computing unit 202 is established unit 203 with the polar coordinate system and is connect, the polar coordinate system establish unit 203 also with the scale
Shape unit 204 and the polar diameter unit 205 connection, the scale shape unit 204 and the polar diameter unit 205 also with it is described
Division unit 206 connects;The detection unit 201, for selecting to reconstruction region, and the detection shape to reconstruction region
Shape;The computing unit 202, for find it is described to the geometric center of reconstruction region at;The polar coordinate system establishes unit
203, for establish two-dimentional polar coordinate system as pole at the geometric center;The scale shape unit 204, for described
In two-dimentional polar coordinate system, as pole, to be reduced to the shape by several scales at the geometric center, obtain several rulers
Spend shape;The polar diameter unit 205, for being existed as side according to set angle using ray where the polar axis using pole as initial point
It is rotated in the two dimension polar coordinate system, obtains several polar diameters;The division unit 206, for several polar diameters respectively by institute
It states several scale shapes and is divided into several regions, several regions are reconstruction region;Wherein, it is described to reconstruction region be need
Measure the region of velocity of blood flow distribution.Concrete implementation mode can refer to the detailed description in Fig. 1 or Fig. 2.
Wherein, the detection unit 201, comprising: equivalent unit;The equivalent unit, being used for will be described to reconstruction region
Cross section be equivalent to justify, the shape to reconstruction region be circle;Or the equivalent unit, for by the area to be reconstructed
The cross section in domain is equivalent to justify, and the scale shape unit 204 completes following operation: determining the number of several scale shapes
H successively makees concentric circles by radius of k*R/h using the pole as the center of circle outward, obtains several scale shapes;Wherein, institute
Stating several scale shapes is circular shape;R is the radius to reconstruction region, coefficient k=1,2...h.Concrete implementation side
Formula can refer to the detailed description in Fig. 1 or Fig. 2.
A kind of reconstruct equation of the device of velocity of blood flow distribution measuring are as follows:
Wherein, V1…VMRespectively indicate the 1st to m-th reconstruction region average velocity of axes, WNMIndicate m-th region
Flow the contribution weight to n-th induced electricity potential difference, A1...AMIndicate the 1st to m-th reconstruction region area, Δ U1...Δ
UNIndicate the 1st to n-th measure induced electricity potential difference.
Wherein, different zones flowing is calculated according to the position of measuring electrode and reconstruction region and induced electricity potential difference is caused to change
Weighted value, calculation formula is as follows:
In formula, ψinIndicate the angle between pole and the line and polar axis of measuring electrode position, ψoutIndicate pole and ginseng
The angle between the line of electrode position and polar axis is examined, (r, θ) is indicated to establish as pole to reconstruction region geometric center
The coordinate of the geometric center of reconstruction region under polar coordinates, t indicate iteration coefficient, and value range is 1 to infinity, and B indicates permanent
Fixed-field size, R indicate the section radius of human upper limb.
Fig. 5 is upper limb velocity of blood flow distribution reconstructing method schematic illustration of the embodiment of the present invention based on electromagnetic induction.Such as
Shown in Figure 4 and 5, it is illustrated using human upper limb (arm) as to reconstruction region, in Fig. 5, the first hot-wire coil a,
Skin f, Fig. 5 the illustrating for Fig. 4 of second hot-wire coil b, measuring electrode c, arteriosomes d, human upper limb.
In figures 4 and 5, step 401 is in the magnetic field of the cross-section setting constant homogeneous to reconstruction region.Specifically
It says, is passed through electric current to the first hot-wire coil a and the second hot-wire coil b, generating one by the second hot-wire coil b to reconstruction region
To the stationary magnetic field B in the first direction hot-wire coil a.
In figures 4 and 5, step 402 is by n electrode arrangements in the outside of the cross section, the outside of the cross section
On the outside of the skin to reconstruction region, such as: when reconstruction region is human upper limb, n electrode arrangements skin in human upper limb
Outside.Specifically, n electrode is 16 electrodes, 16 electrode c are evenly spaced on the outside of the skin f of human upper limb, 16
The label of a electrode is respectively the tenth No. five No. zero electrode e0, number one electrode e1... electrode e15, every 2 electrode gaps
22.5 degree.When arteriosomes d is there are when blood flow, arterial blood makees cutting magnetic induction line fortune along perpendicular to flowing in terms of paper
Dynamic, according to electromagnetic induction principle, this can generate a stable induced electricity potential field in upper limb cross-section, and 16 electrodes can detect
To induced potential.Using an electrode in two electrodes as measuring electrode, another is linked into letter as reference electrode
In number acquisition device, required induced potential difference signal can be obtained.
In figures 4 and 5, for step 403 using one of electrode as reference electrode, remaining n-1 electrode is measuring electrode,
Obtain n-1 induced electricity potential difference.Specifically, n electrode is 16 electrodes, first using No. zero electrode e0 as reference electrode,
Remaining 15 electrode is measuring electrode (the tenth No. five electrode e15 of number one electrode e1...), will obtain 15 potential difference letters
Number.
In figures 4 and 5, step 404 chooses the neighbouring electrode conduct of the reference electrode by clockwise or counter-clockwise
New reference electrode, remaining n-1 electrode are measuring electrode, obtain n-1 induced electricity potential difference.Specifically, again with number one
Electrode e1 is reference electrode, remaining 15 electrode (the tenth No. five electrode e15 of No. zero electrode e0, No. second electrode e2...) is
Measuring electrode obtains 15 potential difference signals again.
In figures 4 and 5, step 405 is until traversing all electrodes is new reference electrode, and finally acquisition n* (n-1) is a
Induced electricity potential difference.Specifically, being next successively with the tenth No. five No. second electrode e2, third electrode e3... electrode e15
Reference electrode finally obtains 240 induced potential difference signals.Traditional measurement method is only with reference to electricity with one of electrode
Pole, can only obtain 15 induced potential difference signals, and this method expands than the induced potential difference signal quantity that traditional measurement mode obtains
16 times are increased.In electromagnetic induction measuring technique, effective induced potential difference data number is more, and measurement accuracy is higher.
Human upper limb cross section is such as divided into M reconstruction region, by taking Fig. 3 division mode as an example, will be drawn to reconstruction region
It is divided into 64 reconstruction regions, wherein the number h=3 and the set angle of several scale shapes are 2 π/16.As shown in figure 3,
Before the blood flowrate for calculating each reconstruction region, is divided to described to reconstruction region, obtain several regions, it is described
Several regions are reconstruction region, to calculate the velocity of blood flow in each reconstruction region to a induced electricity potential difference of the n* (n-1)
Contribute weight.The label mode of each reconstruction region are as follows: polar axis, the first polar diameter and smallest dimension shape intersecting area are the first weight
Structure region 1, polar axis, the first polar diameter and neighbouring next scale shape intersecting area are the second reconstruction region 2, are followed successively by from inside to outside
Third reconstructed region 3, fourth reconstructed region 4, later according to sequence counter-clockwise, by the first polar diameter, the second polar diameter and smallest dimension
Shape intersecting area is the 5th reconstruction region 5, is followed successively by the 6th reconstruction region the 6, the 7th from inside to outside according to above-mentioned same procedure
Reconstruction region 7 and the 8th reconstruction region 8 are arranged successively label 9 according still further to above method sequence counter-clockwise, and 10....64 is (on figure
Do not provide) it is divided into the 9th reconstruction region and for the 64th reconstruction region of the tenth reconstruction region ....In 64 reconstruction regions, often
Velocity of blood flow in a reconstruction region is to the contribution weight of a induced electricity potential difference of the 16* (16-1), with the quantity n of above-mentioned electrode
For=16.
In figures 4 and 5, a induced electricity potential difference of n* described in step 406 (n-1) is as input data, by reconstructing equation,
The blood flowrate of each reconstruction region is calculated, and then obtains the velocity of blood flow to reconstruction region and is distributed;Wherein, it is described to
Reconstruction region is region to be measured, n > 2.Specifically, 240 induced electricity potential differences that measurement is obtained are as input data
Following reconstruct equations are imported, the blood flowrate of each reconstruction region is calculated, and then obtain the blood at entire human upper limb section
Velocity flow profile.
Wherein, V1...VMRespectively indicate the 1st to m-th reconstruction region average velocity of axes, WNMIndicate m-th region
Contribution weight of the flowing to n-th induced electricity potential difference, A1...AMIndicate the 1st to m-th reconstruction region area, Δ U1...
ΔUNIndicate the 1st to n-th measure induced electricity potential difference.
Wherein, different zones flowing is calculated according to the position of measuring electrode and reconstruction region and induced electricity potential difference is caused to change
Weighted value, calculation formula is as follows:
In formula, ψinIndicate the angle between pole and the line and polar axis of measuring electrode position, ψoutIndicate pole and ginseng
The angle between the line of electrode position and polar axis is examined, (r, θ) is indicated to establish as pole to reconstruction region geometric center
The coordinate of the geometric center of reconstruction region under polar coordinates, t indicate iteration coefficient, and value range is 1 to infinity, and B indicates permanent
Fixed-field size;The cross section to reconstruction region is equivalent to justify, the shape to reconstruction region be circle, R indicate to
The section radius of the cross section of reconstruction region, that is to say, that be the section radius of human upper limb.
Meanwhile the cross section to reconstruction region is equivalent to justify, the shape to reconstruction region is circle;M weight
Structure region area AMCalculation method is as follows:
Induced electricity potential difference Δ UNIt is obtained by actual measurement, the W being finally calculatedNMAnd AMIt is updated in reconstruct equation,
The velocity amplitude V of each reconstruction region can be obtainedM, finally obtain the velocity of blood flow distribution to reconstruction region.In Fig. 5, first is logical
Electric coil a, the second hot-wire coil b, measuring electrode c, arteriosomes d, human upper limb skin f.
Embodiment described above is only to express embodiments of the present invention, and the description thereof is more specific and detailed, but can not
Therefore limitations on the scope of the patent of the present invention are interpreted as.It should be pointed out that those skilled in the art, In
Under the premise of not departing from present inventive concept, several deformations, same replacement can also be made, improved etc., these belong to the present invention
Protection scope.Therefore, the scope of protection of the patent of the invention shall be subject to the appended claims.
Claims (10)
1. a kind of method of velocity of blood flow distribution measuring characterized by comprising
In the magnetic field of the cross-section setting constant homogeneous to reconstruction region;
By n electrode arrangements the cross section outside;
Using one of electrode as reference electrode, remaining n-1 electrode is measuring electrode, obtains n-1 induced electricity potential difference;
By clockwise or counter-clockwise, choosing the neighbouring electrode of the reference electrode as new reference electrode, remaining n-1 electricity
Extremely measuring electrode obtains n-1 induced electricity potential difference;
Until traversing all electrodes is new reference electrode, n* (n-1) a induced electricity potential difference is finally obtained;
A induced electricity potential difference of n* (n-1) calculates the blood stream of each reconstruction region by reconstructing equation as input data
Speed value, and then obtain the velocity of blood flow to reconstruction region and be distributed;
Wherein, described to reconstruction region is region to be measured, n > 2.
2. according to the method described in claim 1, it is characterized by:
Before the blood flowrate for calculating each reconstruction region, is divided to described to reconstruction region, obtains several regions,
Several regions are reconstruction region, to calculate the velocity of blood flow in each reconstruction region to the n* (n-1) a induced potential
The contribution weight of difference.
3. according to the method described in claim 2, it is characterized by:
The division methods of the reconstruction region are as follows:
Selection is described to reconstruction region, and the detection shape to reconstruction region;
Find it is described to the geometric center of reconstruction region at;
To establish two-dimentional polar coordinate system as pole at the geometric center;
In the two-dimentional polar coordinate system, as pole, to be reduced to the shape by several scales at the geometric center,
Obtain several scale shapes;
Using pole as initial point, using the polar axis as side, rotates, obtain several in the two-dimentional polar coordinate system according to set angle
Polar diameter;
Several scale shapes are divided into several regions respectively by several polar diameters, and several regions are reconstruction region;
Wherein, it is described to reconstruction region be need to measure velocity of blood flow distribution region.
4. according to the method described in claim 3, it is characterized by:
The cross section to reconstruction region is equivalent to justify, the shape to reconstruction region is circle;Or
The cross section to reconstruction region is equivalent to justify, the shape to reconstruction region is circle;In the two-dimentional pole
In coordinate system, as pole, to be reduced to the shape by several scales at the geometric center, obtain several scale shapes
Method particularly includes: the number h for determining several scale shapes, using the pole as the center of circle, using k*R/h as radius successively to
It is outer to make concentric circles, obtain several scale shapes;
Wherein, several scale shapes are circular shape;R is the radius to reconstruction region, coefficient k=1,2...h.
5. according to the described in any item methods of claim 2-4, which is characterized in that the reconstruct equation are as follows:
Wherein, VMIndicate the average velocity of axes of m-th reconstruction region, WNMIndicate the flowing in m-th region to n-th induced electricity
The contribution weight of potential difference, AMIndicate the area of m-th reconstruction region, Δ UNIndicate that n-th measures induced electricity potential difference.
6. a kind of device of velocity of blood flow distribution measuring characterized by comprising
Magnetic field units, in the magnetic field of the cross-section setting constant homogeneous to reconstruction region;
Electrode unit, by n electrode arrangements the cross section outside;
Induced electricity potential difference acquiring unit is used for using one of electrode as reference electrode, remaining n-1 electrode is measuring electrode,
Obtain n-1 induced electricity potential difference;
Updating unit chooses the neighbouring electrode of the reference electrode as new reference electrode by clockwise or counter-clockwise,
Remaining n-1 electrode is measuring electrode, obtains n-1 induced electricity potential difference;
Traversal Unit finally obtains n* (n-1) a induced electricity potential difference until traversing all electrodes is new reference electrode;
Computing unit, a induced electricity potential difference of n* (n-1) calculate each reconstruct area by reconstructing equation as input data
The blood flowrate in domain, and then obtain the velocity of blood flow to reconstruction region and be distributed;
Wherein, described to reconstruction region is region to be measured, n > 2.
7. device according to claim 6, which is characterized in that further include:
Reconstruction region division unit or device;
The reconstruction region division unit or device, before the blood flowrate for calculating each reconstruction region, to it is described to
Reconstruction region is divided, and obtains several regions, several regions are reconstruction region, to calculate the blood in each reconstruction region
Contribution weight of the flow velocity to a induced electricity potential difference of the n* (n-1).
8. device according to claim 7, it is characterised in that:
The reconstruction region division unit or device, comprising:
Detection unit (201), computing unit (202), polar coordinate system establish unit (203), scale shape unit (204) polar diameter list
First (205) and division unit (206);The detection unit (201) respectively with the computing unit (202) and the scale shape
Unit (204) connection, the computing unit (202) are established unit (203) with the polar coordinate system and are connect, and the polar coordinate system is built
Vertical unit (203) are also connect with the scale shape unit (204) and the polar diameter unit (205), the scale shape unit
(204) it is also connect with the division unit (206) with the polar diameter unit (205);
The detection unit (201), for selecting to reconstruction region, and the detection shape to reconstruction region;
The computing unit (202), for find it is described to the geometric center of reconstruction region at;
The polar coordinate system establishes unit (203), for establish two-dimentional polar coordinate system as pole at the geometric center;
The scale shape unit (204), in the two-dimentional polar coordinate system, using at the geometric center as pole, it is right
The shape is reduced by several scales, obtains several scale shapes;
The polar diameter unit (205) is used for using pole as initial point, using the polar axis as side, according to set angle in the two dimension
It is rotated in polar coordinate system, obtains several polar diameters;
Several scale shapes are divided into several regions, institute respectively for several polar diameters by the division unit (206)
Stating several regions is reconstruction region;
Wherein, it is described to reconstruction region be need to measure velocity of blood flow distribution region.
9. device according to claim 8, it is characterised in that:
The detection unit (201), comprising: equivalent unit;
The equivalent unit, for the cross section to reconstruction region to be equivalent to justify, the shape to reconstruction region is
It is round;Or
The equivalent unit, for the cross section to reconstruction region to be equivalent to justify, the scale shape unit (204) is complete
At following operation: determine the number h of several scale shapes, it is successively outside by radius of k*R/h using the pole as the center of circle
Make concentric circles, obtains several scale shapes;
Wherein, several scale shapes are circular shape;R is the radius to reconstruction region, coefficient k=1,2...h.
10. according to the described in any item devices of claim 7-9, which is characterized in that the reconstruct equation are as follows:
Wherein, VMIndicate the average velocity of axes of m-th reconstruction region, WNMIndicate the flowing in m-th region to n-th induced electricity
The contribution weight of potential difference, AMIndicate the area of m-th reconstruction region, Δ UNIndicate that n-th measures induced electricity potential difference.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910657334.9A CN110464329B (en) | 2019-07-19 | 2019-07-19 | Method and device for measuring blood flow velocity distribution |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910657334.9A CN110464329B (en) | 2019-07-19 | 2019-07-19 | Method and device for measuring blood flow velocity distribution |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110464329A true CN110464329A (en) | 2019-11-19 |
CN110464329B CN110464329B (en) | 2020-06-02 |
Family
ID=68509568
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910657334.9A Active CN110464329B (en) | 2019-07-19 | 2019-07-19 | Method and device for measuring blood flow velocity distribution |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110464329B (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112120686A (en) * | 2020-09-21 | 2020-12-25 | 东北大学 | Non-contact artery imaging method and device, electronic device and storage medium |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN201408037Y (en) * | 2009-02-26 | 2010-02-17 | 江苏横河自控设备有限公司 | Independent excitation multi-electrode electromagnetic flowmeter |
CN105169526A (en) * | 2015-09-29 | 2015-12-23 | 青岛歌尔声学科技有限公司 | Infusion control device and method |
CN105344012A (en) * | 2015-11-30 | 2016-02-24 | 燕山大学 | Transcranial magnetic stimulation and laser doppler velocity measurement integrated apparatus |
CN105701824A (en) * | 2016-01-14 | 2016-06-22 | 天津大学 | Electrical tomographic image reconstruction method based on circular-arc grid subdivision |
CN106039441A (en) * | 2016-05-12 | 2016-10-26 | 北京精密机电控制设备研究所 | Sensor-free flow rate measuring method and measuring device for ventricular assisting device |
CN107028614A (en) * | 2017-05-24 | 2017-08-11 | 吴佩萱 | A kind of wearable anaemia device for fast detecting based on electromagnetic induction effect |
CN109106360A (en) * | 2018-08-15 | 2019-01-01 | 河北科技大学 | Blood of human body flow-speed measurement method and its device |
CN109416268A (en) * | 2016-06-28 | 2019-03-01 | 西门子股份公司 | Magnetic induction flowmeter |
CN109700452A (en) * | 2018-12-06 | 2019-05-03 | 河北科技大学 | A kind of velocity of blood flow distribution measurement method of limbs |
-
2019
- 2019-07-19 CN CN201910657334.9A patent/CN110464329B/en active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN201408037Y (en) * | 2009-02-26 | 2010-02-17 | 江苏横河自控设备有限公司 | Independent excitation multi-electrode electromagnetic flowmeter |
CN105169526A (en) * | 2015-09-29 | 2015-12-23 | 青岛歌尔声学科技有限公司 | Infusion control device and method |
CN105344012A (en) * | 2015-11-30 | 2016-02-24 | 燕山大学 | Transcranial magnetic stimulation and laser doppler velocity measurement integrated apparatus |
CN105701824A (en) * | 2016-01-14 | 2016-06-22 | 天津大学 | Electrical tomographic image reconstruction method based on circular-arc grid subdivision |
CN106039441A (en) * | 2016-05-12 | 2016-10-26 | 北京精密机电控制设备研究所 | Sensor-free flow rate measuring method and measuring device for ventricular assisting device |
CN109416268A (en) * | 2016-06-28 | 2019-03-01 | 西门子股份公司 | Magnetic induction flowmeter |
CN107028614A (en) * | 2017-05-24 | 2017-08-11 | 吴佩萱 | A kind of wearable anaemia device for fast detecting based on electromagnetic induction effect |
CN109106360A (en) * | 2018-08-15 | 2019-01-01 | 河北科技大学 | Blood of human body flow-speed measurement method and its device |
CN109700452A (en) * | 2018-12-06 | 2019-05-03 | 河北科技大学 | A kind of velocity of blood flow distribution measurement method of limbs |
Non-Patent Citations (1)
Title |
---|
赵宇洋: "多电极电磁流量计肢体血液流速分布测量研究", 《河北科技大学学报》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112120686A (en) * | 2020-09-21 | 2020-12-25 | 东北大学 | Non-contact artery imaging method and device, electronic device and storage medium |
CN112120686B (en) * | 2020-09-21 | 2021-07-13 | 东北大学 | Non-contact artery imaging method and device, electronic device and storage medium |
Also Published As
Publication number | Publication date |
---|---|
CN110464329B (en) | 2020-06-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Harrach et al. | Resolution guarantees in electrical impedance tomography | |
CN102435637B (en) | Perform the system and method for electrical impedance tomography | |
US8963562B2 (en) | Transducer configurations and methods for transducer positioning in electrical impedance tomography | |
Kim et al. | Experimental measurements of turbulence spectra distal to stenoses | |
Nicholson et al. | Modelling and experimental measurements of idealised and light-moderate RCF cracks in rails using an ACFM sensor | |
Muñoz et al. | Railroad inspection based on ACFM employing a non-uniform B-spline approach | |
Guardo et al. | An experimental study in electrical impedance tomography using backprojection reconstruction | |
CN101794453B (en) | Reconstruction method of node mapping image based on regression analysis | |
Wei et al. | Theoretical and experimental evaluation of rotational magnetic induction tomography | |
CN110464329A (en) | A kind of method and apparatus of velocity of blood flow distribution measuring | |
CN110490925A (en) | A kind of method and apparatus that reconstruction region divides | |
US20080275657A1 (en) | Method for noninvasive determination of a distribution of electrical current and apparatus for the same | |
CN105433912A (en) | Magnetic-nano real-time non-invasive temperature measurement method | |
CN105158326B (en) | For measuring the detection method of martensite content in austenitic stainless steel | |
CN109946388B (en) | Electrical/ultrasonic bimodal inclusion boundary reconstruction method based on statistical inversion | |
Etuke et al. | Measurement of angular velocities using electrical impedance tomography | |
Schneider | Improved methods for measuring laminar-turbulent intermittency in boundary layers | |
Cao et al. | 2D image reconstruction of a human chest by using Calderon's method and the adjacent current pattern | |
Shi et al. | Analytical investigation of an inductive flow sensor with arc-shaped electrodes for water velocity measurement in two-phase flows | |
Bakhtiari et al. | Modeling of eddy current probe response for steam generator tubes | |
CN108254418A (en) | The three-dimensional variational method of capacitance and ultrasound tomography multi-modal fusion | |
Yang et al. | Real-time Monitoring Multi-phase Flow using Electrical Impedance Tomography | |
CN116744851A (en) | Quality metric of mapping function | |
Zhu et al. | Simultaneous reconstruction of conductivity and permittivity in electrical impedance tomography | |
JPH0221289A (en) | Detection of buried pipe |
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 |