CN114159026A - Evaluation method for microcirculation function of kidney - Google Patents
Evaluation method for microcirculation function of kidney Download PDFInfo
- Publication number
- CN114159026A CN114159026A CN202111410368.1A CN202111410368A CN114159026A CN 114159026 A CN114159026 A CN 114159026A CN 202111410368 A CN202111410368 A CN 202111410368A CN 114159026 A CN114159026 A CN 114159026A
- Authority
- CN
- China
- Prior art keywords
- microcirculation
- renal
- microcirculation function
- kidney
- function
- 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.)
- Pending
Links
- 230000004089 microcirculation Effects 0.000 title claims abstract description 202
- 210000003734 kidney Anatomy 0.000 title claims abstract description 61
- 238000011156 evaluation Methods 0.000 title claims description 9
- 238000000034 method Methods 0.000 claims abstract description 24
- 230000017531 blood circulation Effects 0.000 claims abstract description 22
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 claims abstract description 18
- 229910052760 oxygen Inorganic materials 0.000 claims abstract description 18
- 239000001301 oxygen Substances 0.000 claims abstract description 18
- 238000012545 processing Methods 0.000 claims abstract description 17
- 230000010412 perfusion Effects 0.000 claims abstract description 15
- 238000012800 visualization Methods 0.000 claims abstract description 15
- 238000007781 pre-processing Methods 0.000 claims abstract description 4
- 238000010183 spectrum analysis Methods 0.000 claims abstract description 4
- 230000008081 blood perfusion Effects 0.000 claims description 18
- 210000001519 tissue Anatomy 0.000 claims description 13
- MWUXSHHQAYIFBG-UHFFFAOYSA-N Nitric oxide Chemical compound O=[N] MWUXSHHQAYIFBG-UHFFFAOYSA-N 0.000 claims description 12
- 210000002889 endothelial cell Anatomy 0.000 claims description 7
- 239000000523 sample Substances 0.000 claims description 6
- 210000003743 erythrocyte Anatomy 0.000 claims description 5
- 238000004140 cleaning Methods 0.000 claims description 3
- 238000010276 construction Methods 0.000 claims description 3
- 230000001419 dependent effect Effects 0.000 claims description 3
- 239000006185 dispersion Substances 0.000 claims description 3
- 238000005259 measurement Methods 0.000 claims description 3
- 210000003205 muscle Anatomy 0.000 claims description 3
- 210000005036 nerve Anatomy 0.000 claims description 3
- 230000003287 optical effect Effects 0.000 claims description 3
- 238000013139 quantization Methods 0.000 claims description 3
- 238000011425 standardization method Methods 0.000 claims description 3
- 238000010219 correlation analysis Methods 0.000 claims description 2
- 230000000241 respiratory effect Effects 0.000 claims description 2
- 230000000747 cardiac effect Effects 0.000 claims 1
- 238000010586 diagram Methods 0.000 abstract description 6
- 210000004369 blood Anatomy 0.000 abstract description 2
- 239000008280 blood Substances 0.000 abstract description 2
- 238000002562 urinalysis Methods 0.000 abstract 1
- 230000006870 function Effects 0.000 description 96
- 210000002700 urine Anatomy 0.000 description 7
- 238000004458 analytical method Methods 0.000 description 3
- 208000002352 blister Diseases 0.000 description 2
- 230000004064 dysfunction Effects 0.000 description 2
- 210000004072 lung Anatomy 0.000 description 2
- 210000000056 organ Anatomy 0.000 description 2
- 230000035790 physiological processes and functions Effects 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 206010023421 Kidney fibrosis Diseases 0.000 description 1
- 230000002159 abnormal effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- -1 blood flow speed Substances 0.000 description 1
- 230000036770 blood supply Effects 0.000 description 1
- 210000004556 brain Anatomy 0.000 description 1
- 230000001413 cellular effect Effects 0.000 description 1
- 238000009535 clinical urine test Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 206010061989 glomerulosclerosis Diseases 0.000 description 1
- 230000000004 hemodynamic effect Effects 0.000 description 1
- 230000003907 kidney function Effects 0.000 description 1
- 230000002503 metabolic effect Effects 0.000 description 1
- 210000004925 microvascular endothelial cell Anatomy 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000001575 pathological effect Effects 0.000 description 1
- 210000002381 plasma Anatomy 0.000 description 1
- 230000008085 renal dysfunction Effects 0.000 description 1
- 230000029058 respiratory gaseous exchange Effects 0.000 description 1
- 238000010008 shearing Methods 0.000 description 1
- 210000003491 skin Anatomy 0.000 description 1
- 210000000813 small intestine Anatomy 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 210000000115 thoracic cavity Anatomy 0.000 description 1
- 210000005239 tubule Anatomy 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
- 238000007794 visualization technique Methods 0.000 description 1
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/0059—Measuring for diagnostic purposes; Identification of persons using light, e.g. diagnosis by transillumination, diascopy, fluorescence
- A61B5/0082—Measuring for diagnostic purposes; Identification of persons using light, e.g. diagnosis by transillumination, diascopy, fluorescence adapted for particular medical purposes
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/0059—Measuring for diagnostic purposes; Identification of persons using light, e.g. diagnosis by transillumination, diascopy, fluorescence
- A61B5/0062—Arrangements for scanning
- A61B5/0064—Body surface scanning
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/02—Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
- A61B5/026—Measuring blood flow
- A61B5/0261—Measuring blood flow using optical means, e.g. infrared light
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/20—Measuring for diagnostic purposes; Identification of persons for measuring urological functions restricted to the evaluation of the urinary system
- A61B5/201—Assessing renal or kidney functions
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7235—Details of waveform analysis
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7235—Details of waveform analysis
- A61B5/7253—Details of waveform analysis characterised by using transforms
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Molecular Biology (AREA)
- Animal Behavior & Ethology (AREA)
- Veterinary Medicine (AREA)
- Public Health (AREA)
- Physics & Mathematics (AREA)
- General Health & Medical Sciences (AREA)
- Biophysics (AREA)
- Pathology (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Surgery (AREA)
- Physiology (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Signal Processing (AREA)
- Psychiatry (AREA)
- Hematology (AREA)
- Urology & Nephrology (AREA)
- Cardiology (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Radiology & Medical Imaging (AREA)
- Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)
- Investigating Or Analysing Biological Materials (AREA)
- Measuring And Recording Apparatus For Diagnosis (AREA)
Abstract
The invention belongs to the technical field of medical data processing, and provides a method for evaluating the microcirculation function of kidneys, which comprises the following steps: the method comprises the steps of measuring the renal microcirculation function parameters, preprocessing renal microcirculation function data and visualizing the renal microcirculation function, and analyzing the renal microcirculation function from the depth of the whole layer by constructing a three-dimensional visualization model, a renal microcirculation function bubble diagram and wavelet transform spectrum analysis. Compared with the prior art, the kidney microcirculation function in the scheme of the invention comprises the indexes of the perfusion, the blood flow speed and the microcirculation oxygen partial pressure of the kidney microcirculation blood. Compared with the prior art of indirectly evaluating the renal microcirculation function through histology and urinalysis indexes, the method can evaluate the renal microcirculation function more directly and comprehensively.
Description
Technical Field
The invention belongs to the technical field of medical data processing, and particularly relates to a method for evaluating the microcirculation function of kidneys.
Background
Different from the parallel distribution of the heart microcirculation series structure and the brain microcirculation, the kidney microcirculation structure is unique and is in a saccule-shaped reticular structure. The renal microcirculation is involved in maintaining the normal physiological function of the renal tubules. Renal microcirculation dysfunction can directly affect the blood supply and distribution of glomeruli and renal interstitium, leading to glomerulosclerosis and renal interstitial fibrosis. Therefore, accurate assessment of renal microcirculatory dysfunction is helpful for the diagnosis and treatment of renal dysfunction in early stages of clinical practice.
As shown in fig. 1, in the prior art, a renal cortex or renal medulla specimen is mainly prepared, and the renal microcirculation is observed by using a microscope, a scanning electron microscope or a transmission electron microscope; or the renal microcirculation function is indirectly reflected by urine examination indexes such as urine microalbumin, urine microglobulin and the like.
The existing technology for evaluating the function of the microcirculation of the kidney has the following objective defects:
(1) the functional index is absent. The prior art can only measure the kidney histological structure, and can not evaluate microcirculation function indexes such as kidney microcirculation micro-hemodynamics, kidney microcirculation oxygen and the like.
(2) And (4) indirection. In the prior art, urine detection indexes such as urine microalbumin and urine microglobulin only can indirectly reflect the microcirculation function of the kidney, but cannot directly reflect functional indexes such as microcirculation blood perfusion, microcirculation oxygen distribution and the like of the kidney.
Therefore, the prior art method is difficult to comprehensively understand the microvascular and microvascular dynamics and functions of the living kidney. How to comprehensively and scientifically show the microcirculation function characteristics of the kidney is a problem to be solved.
Disclosure of Invention
Aiming at the problems in the prior art, the invention aims to provide a method for evaluating the function of the microcirculation of the kidney, which realizes the description and evaluation of the function characteristics of the microcirculation of the kidney by capturing and analyzing functional parameters such as perfusion of the microcirculation blood, blood flow speed, oxygen partial pressure of the microcirculation and the like.
The purpose of the invention is realized by the following technical scheme:
a method for evaluating the microcirculation function of kidney includes such steps as:
utilize the kidney microcirculation function evaluation device who comprises microcirculation oxygen partial pressure monitor and binary channels laser Doppler monitor, gather multiple kidney microcirculation function parameter simultaneously, kidney microcirculation function parameter includes kidney microcirculation tissue oxygen partial Pressure (PO)2) Renal microcirculation blood perfusion level (BP) and blood perfusion rate (V);
importing the data of the renal microcirculation function parameters measured in the step 1 into a data cleaning module, and processing the outlier renal microcirculation function parameters through a boxplot algorithm;
step 3.1, three-dimensional visualization model construction
Introducing the renal microcirculation function parameter data preprocessed in the step 2, performing dimensionless processing, and then generating a renal microcirculation function three-dimensional visualization module by applying Python and ECharts to the renal microcirculation function parameter data subjected to dimensionless quantization processing, wherein time, a microcirculation function parameter variable and a microcirculation function parameter variable value are respectively defined as an X axis, a Y axis and a Z axis of the module;
step 3.2, renal microcirculation function bubble chart
Introducing the renal microcirculation function parameter data pretreated in the step 2, and generating a renal microcirculation function uniaxial bubble chart by using ECharts to realize renal microcirculation function parameter distribution and weight correlation analysis;
step 3.3, wavelet transform Spectrum analysis
Importing the renal microcirculation function parameter data preprocessed in the step 2, dividing the renal microcirculation function parameter data into n frequency band source signal amplitudes according to the microcirculation blood flow perfusion signal frequency, and generating a two-dimensional distribution map of the renal microvascular blood flow perfusion, the blood flow velocity and the erythrocyte density signals in a characteristic frequency band through the n frequency band source signal amplitudes; drawing a three-dimensional time-frequency graph of the renal microcirculation blood perfusion signal according to three dimensional indexes of time, frequency and microcirculation characteristic source amplitude; the distribution pattern of the relevant magnitude of the renal microcirculation function is shown by radar maps.
Further, step 1 utilizes when kidney microcirculation function evaluation device carries out the collection of kidney microcirculation function parameter, integrate the light probe of two equipment to micron stereotaxic instrument, utilize micron stereotaxic instrument to guide the light probe of two equipment to kidney top 1mm department, use two kinds of light sources of white light, laser, gather multiple kidney microcirculation function parameter simultaneously.
Further, the specific method for processing the outlier renal microcirculation function parameter in the step 2 is as follows: definition of Q125% maximum, Q3Is 75% of maximum, Q3And Q1The difference between is the interquartile range (IQR) and is set to (Q)1-1.5 XIQR) and (Q)3+1.5 × IQR) is the microcirculation function parameter boundary value, and the microcirculation function data beyond the boundary value is regarded as outlier and adjusted to the normal range boundary value.
Further, the concrete method of the dimensionless processing in step 3.1 is as follows: and processing the renal microcirculation function parameter data by a dispersion standardization method, eliminating dimensional dimensions of multi-parameter data, and projecting the renal microcirculation function parameter data in a [0,1] interval to realize the visualization of the renal microcirculation function of the same universal coordinate system frame.
Further, the horizontal axis of the uniaxial bubble chart in step 3.2 represents the value range of the distribution of the renal microcirculation function parameter data, the value range is divided into m sections, continuous sections are marked by using numbers, and the distribution weight correlation of the renal microcirculation function parameter data is reflected.
Further, in step 3.3, the microcirculation blood flow perfusion signal frequency can be divided into 6 frequency range source signal amplitudes, including a heart source amplitude of 2-5Hz, a respiration source amplitude of 0.4-2Hz, a muscle source amplitude of 0.15-0.4Hz, a nerve source amplitude of 0.04-0.15Hz, a nitric oxide dependent endothelial cell source amplitude of 0.0095-0.04Hz and a nitric oxide independent endothelial cell source amplitude of 0.005-0.0095 Hz.
Compared with the prior art, the invention has the beneficial effects that:
1. in the method for evaluating the renal microcirculation function, the renal microcirculation function comprises the indexes of renal microcirculation blood perfusion, blood flow speed and microcirculation oxygen partial pressure, and compared with the prior art of indirectly evaluating the renal microcirculation function through histology and urine test indexes, the method can more directly and comprehensively evaluate the renal microcirculation function;
2. according to the method for evaluating the renal microcirculation function, null data, abnormal data and outliers from instrument firmware and optical components are removed from the microcirculation function data through a computer algorithm, so that the scientificity and the accuracy of the renal microcirculation function data are improved;
3. the study on the microcirculation function of the kidney in the prior art does not contain a visualization scheme; the method for evaluating the renal microcirculation function disclosed by the invention is used for analyzing the renal microcirculation function in a deep manner from the whole aspect by using a general microcirculation framework and wavelet analysis, realizes the visual presentation of the renal microcirculation function and provides a new method for completely, comprehensively, scientifically and accurately evaluating the renal microcirculation function.
Drawings
The invention is further illustrated by the following examples in conjunction with the accompanying drawings:
FIG. 1 is a prior art technical scheme for a study of renal microcirculation;
FIG. 2 is a flow chart of a method for assessing renal microcirculation function according to the present invention;
FIG. 3 is a schematic diagram of a method for assessing renal microcirculation function according to the present invention;
FIG. 4 is a three-dimensional visualization view of the function of the microcirculation of the kidney according to example 1; in the figure, (a) a three-dimensional visualization histogram of the renal microcirculation, (b) a three-dimensional visualization scatter diagram of the renal microcirculation; PO (PO)2: microcirculation tissue oxygen partial pressure, BP: microcirculation bloodPerfusion level, V: perfusion rate of microcirculation blood flow;
FIG. 5 is a graph of uniaxial bubbles generated from test data;
FIG. 6 is a graph of uniaxial blebs of renal microcirculation function in example 1; in the figure, PO2: microcirculation tissue oxygen partial pressure, BP: microcirculation blood perfusion level, V: perfusion rate of microcirculation blood flow;
FIG. 7 is a two-dimensional distribution graph of the renal microvascular perfusion, blood flow velocity and red blood cell density signals in the characteristic frequency band of example 1; in the figure, the black curve is the renal microvascular blood perfusion, the green curve is the microcirculation blood flow velocity, and the red curve is the red blood cell density signal;
FIG. 8 is a three-dimensional time-frequency diagram of a perfusion signal of the microcirculation of the kidney according to example 1; in the figure, the scale of blue (low) to orange (high) represents the continuous wavelet coefficients;
FIG. 9 is a radar chart of the amplitude associated with the microcirculation function of the kidney according to example 1.
Detailed Description
Example 1
As shown in fig. 2 and fig. 3, the present embodiment provides a method for assessing the microcirculation function of kidney, which has the following characteristics compared with the microcirculation function assessment of kidney of other tissues:
1) the kidney is an important organ of the human body, and the main physiological functions are to filter blood plasma, reabsorb useful substances and discharge metabolic wastes in the body in the form of urine. When the kidney is damaged acutely and chronically due to a plurality of factors, the structure of the kidney changes, and the renal function is further influenced;
2) compared with organs and tissues such as lung, small intestine, skin and the like, the kidney microcirculation structure is unique and is of a saccule-shaped reticular structure, and endothelial cells of the kidney are flatly coated on the luminal side of the capillary wall and are contacted with blood flow, so that the kidney microcirculation structure is easily influenced by pathological factors such as shearing force and the like;
3) the micro-circulation function data acquisition access, the light source path and the interference and interference source of the micro-circulation function data are different from superficial tissues (such as skin) and thoracic tissues (such as lung).
The embodiment realizes evaluation of renal microcirculation hemodynamics and microcirculation oxygen by relying on a computer algorithm, establishes a visualization method of renal microcirculation function by relying on a general microcirculation framework and wavelet analysis, and specifically comprises the following steps:
utilize the kidney microcirculation function evaluation device who comprises microcirculation oxygen partial pressure monitor micron TX3 (German Presens company) and binary channels laser Doppler monitor VMS (British Moor company), integrate above-mentioned two equipment probe to micron stereotaxic apparatus, utilize micron stereotaxic apparatus to guide the light probe of two equipment to kidney top 1mm department, with two kinds of light sources of white light and laser, gather multiple kidney microcirculation function parameter (S1) simultaneously, kidney microcirculation function parameter includes kidney microcirculation tissue oxygen partial pressure (PO 1)2) Renal microcirculation blood perfusion level (BP) and blood perfusion rate (V).
and (3) importing the data of the renal microcirculation function parameters measured in the step (1) into a data cleaning module (module), and processing the outlier renal microcirculation function parameters through a boxcar graph algorithm. Definition of Q125% maximum, Q3Is 75% of maximum, Q3And Q1The difference between is the interquartile range (IQR) and is set to (Q)1-1.5 XIQR) and (Q)3+1.5 × IQR) is the microcirculation function parameter boundary value, and the microcirculation function data exceeding the boundary value is regarded as an outlier and adjusted to the normal range boundary value (S2).
step 3.1, three-dimensional visualization model construction
Importing the renal microcirculation function parameter data preprocessed in the step 2, and performing dimensionless processing (S3): and processing the renal microcirculation function parameter data by a dispersion standardization method, eliminating dimensional dimensions of multi-parameter data, projecting the microcirculation function parameter data in a [0,1] interval, and realizing renal microcirculation function visualization of the same universal coordinate system frame (S4).
Specifically, Python and ECharts are applied to the kidney microcirculation function parameter data subjected to dimensionless quantization processing to generate a kidney microcirculation function three-dimensional visualization module, wherein time, a microcirculation function parameter variable and a microcirculation function parameter variable value are respectively defined as an X axis, a Y axis and a Z axis of the module, and a kidney microcirculation function three-dimensional visualization view as shown in fig. 4 is formed.
3.2 uniaxial bleb plot of renal microcirculation function
Introducing the renal microcirculation function parameter data (the renal microcirculation blood perfusion level BP and the blood perfusion speed V) pretreated in the step 2, generating a renal microcirculation function uniaxial bubble map by using ECharts, realizing the analysis of the distribution and weight association of the renal microcirculation function parameters, arranging the data in a form of 'bubbles' according to a horizontal axis, and displaying the distribution weight and the association of the renal microcirculation function data by the area size of the 'bubbles'.
As shown in fig. 5, in the present embodiment, a uniaxial bubble chart is generated by using part of the test data, wherein the horizontal axis in the bubble chart represents the value range of the distribution of the renal microcirculation function parameter data, the value range is divided into 5 sections, and the continuous sections are marked by using numbers, and are respectively the section 0-1, the section 1-2, the section 2-3, the section 3-4 and the section 4-5, which reflect the distribution weight correlation of the renal microcirculation function parameter data; the data distribution quantity of the microcirculation function parameter is in direct proportion to the area of the 'air bubble': the larger the circular area is, the more the microcirculation function data fall in the corresponding interval of the circle; conversely, a smaller circular area indicates less microcirculatory function data in this interval (fig. 5). Therefore, the distribution and weight correlation of the renal microcirculation function data in each value interval can be visually shown in the uniaxial bubble chart (S5).
The renal microcirculation tissue oxygen partial Pressure (PO) is shown in FIG. 62) In the distribution weight correlation bubble chart, the circular areas corresponding to the interval scales of 40-50hPa are larger, the circular areas corresponding to the interval scales of 0-40hPa and 50-60hPa are smaller, which indicates that the oxygen partial pressure data of the biological renal microcirculation tissue are more and more intensively distributed on the interval scales of 40-50hPa, and less data are distributed on the interval scales of 0-40hPa and 50-60hPa, namely the oxygen partial pressure of the renal microcirculation tissue is mainly at the level of 40-50hPa, and the oxygen partial pressure level of the microcirculation tissue at the level of 40-50hPa is the microcirculation functionDominant weighting of energy.
Step 3.3, wavelet transform spectrum analysis:
the present embodiment analyzes the renal microcirculation blood flow signal using wavelet transform. Introducing the renal microcirculation blood flow perfusion data (the renal microcirculation blood flow perfusion level BP and the blood flow perfusion speed V) preprocessed in the step 2, and dividing the data into six frequency band source signal amplitudes according to the microcirculation blood flow perfusion signal frequency, wherein the six frequency band source signal amplitudes comprise a heart source amplitude of 2-5Hz, a respiratory source amplitude of 0.4-2Hz, a muscle source amplitude of 0.15-0.4Hz, a nerve source amplitude of 0.04-0.15Hz, a nitric oxide dependent endothelial cell source amplitude of 0.0095-0.04Hz and a nitric oxide independent endothelial cell source amplitude of 0.005-0.0095Hz (S6).
The kidney microvascular endothelial cells are the most important cellular components of the kidney microcirculation and are the basic functional units of the kidney microcirculation. Through the six frequency band source signal amplitudes, a two-dimensional distribution diagram of the renal microvascular blood perfusion, blood flow velocity and red blood cell density signals in the characteristic frequency band is generated as shown in fig. 7. The two-dimensional distribution diagram of the microcirculation blood flow signal frequency spectrum represents the frequency interval of the amplitude and the characteristic peak of each frequency band.
A three-dimensional time-frequency graph (figure 8) of the renal microcirculation blood perfusion signal is drawn through three-dimensional indexes of time (sec), frequency (Hz) and microcirculation characteristic source Amplitude (AU), and distribution characteristics and change rules of 6 source amplitudes related to the renal microcirculation function along with the time progress are represented.
Six patterns of distribution of amplitude values associated with the function of the renal microcirculation are shown by radar maps as shown in fig. 9.
Finally, it should be noted that the above only illustrates the technical solution of the present invention, but not limited thereto, and although the present invention has been described in detail with reference to the preferred arrangement, it will be understood by those skilled in the art that modifications or equivalent substitutions may be made thereto without departing from the spirit and scope of the technical solution of the present invention.
Claims (6)
1. A method for assessing the function of the microcirculation of the kidney, which comprises the following steps:
step 1, measurement of renal microcirculation function parameters:
utilize the kidney microcirculation function evaluation device who comprises microcirculation oxygen partial pressure monitor and binary channels laser Doppler monitor, gather multiple kidney microcirculation function parameter simultaneously, kidney microcirculation function parameter includes kidney microcirculation tissue oxygen partial Pressure (PO)2) Renal microcirculation blood perfusion level (BP) and blood perfusion rate (V);
step 2, preprocessing the kidney microcirculation function data:
importing the data of the renal microcirculation function parameters measured in the step 1 into a data cleaning module, and processing the outlier renal microcirculation function parameters through a boxplot algorithm;
step 3, visualizing the microcirculation function of the kidney:
step 3.1, three-dimensional visualization model construction
Introducing the renal microcirculation function parameter data preprocessed in the step 2, performing dimensionless processing, and then generating a renal microcirculation function three-dimensional visualization module by applying Python and ECharts to the renal microcirculation function parameter data subjected to dimensionless quantization processing, wherein time, a microcirculation function parameter variable and a microcirculation function parameter variable value are respectively defined as an X axis, a Y axis and a Z axis of the module;
step 3.2, renal microcirculation function uniaxial bubble chart
Introducing the renal microcirculation function parameter data pretreated in the step 2, and generating a renal microcirculation function uniaxial bubble chart by using ECharts to realize renal microcirculation function parameter distribution and weight correlation analysis;
step 3.3, wavelet transform Spectrum analysis
Importing the renal microcirculation function parameter data preprocessed in the step 2, dividing the renal microcirculation function parameter data into n frequency band source signal amplitudes according to the microcirculation blood flow perfusion signal frequency, and generating a two-dimensional distribution map of the renal microvascular blood flow perfusion, the blood flow velocity and the erythrocyte density signals in a characteristic frequency band through the n frequency band source signal amplitudes; drawing a three-dimensional time-frequency graph of the renal microcirculation blood perfusion signal according to three dimensional indexes of time, frequency and microcirculation characteristic source amplitude; the distribution pattern of the relevant magnitude of the renal microcirculation function is shown by radar maps.
2. The method according to claim 1, wherein in the step 1, when the renal microcirculation function parameter is collected by using the renal microcirculation function evaluation device, the optical probes of the two devices are integrated into a micro-positioning machine, the optical probes of the two devices are guided to a position 1mm above the kidney by using the micro-positioning machine, and a plurality of renal microcirculation function parameters are collected simultaneously by using two light sources, namely white light and laser.
3. The method of claim 1, wherein the specific method of processing the outlier renal microcirculation function parameter in step 2 is: definition of Q125% maximum, Q3Is 75% of maximum, Q3And Q1The difference between is the interquartile range (IQR) and is set to (Q)1-1.5 XIQR) and (Q)3+1.5 × IQR) is the microcirculation function parameter boundary value, and the microcirculation function data beyond the boundary value is regarded as outlier and adjusted to the normal range boundary value.
4. The method for assessing renal microcirculation function according to claim 1, wherein the non-dimensionalization in step 3.1 is performed by: and processing the renal microcirculation function parameter data by a dispersion standardization method, eliminating dimensional dimensions of multi-parameter data, and projecting the renal microcirculation function parameter data in a [0,1] interval to realize the visualization of the renal microcirculation function of the same universal coordinate system frame.
5. The method of claim 1, wherein the horizontal axis of the uniaxial bubbly chart in step 3.2 represents the distribution range of the renal microcirculation function parameter data, the distribution range is divided into m sections, and the continuous sections are marked by numbers to reflect the distribution weight correlation of the renal microcirculation function parameter data.
6. The method of claim 1, wherein in step 3.3, the microcirculation blood perfusion signal frequency is divided into 6 frequency range source signal amplitudes, including a cardiac source amplitude of 2-5Hz, a respiratory source amplitude of 0.4-2Hz, a muscle source amplitude of 0.15-0.4Hz, a nerve source amplitude of 0.04-0.15Hz, a nitric oxide dependent endothelial cell source amplitude of 0.0095-0.04Hz and a nitric oxide independent endothelial cell source amplitude of 0.005-0.0095 Hz.
Applications Claiming Priority (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111321092 | 2021-11-09 | ||
| CN202111321092X | 2021-11-09 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN114159026A true CN114159026A (en) | 2022-03-11 |
Family
ID=79248111
Family Applications (4)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202111410230.1A Pending CN113925471A (en) | 2021-11-09 | 2021-11-25 | Skin microcirculation function evaluation and visualization method |
| CN202111410210.4A Pending CN113907724A (en) | 2021-11-09 | 2021-11-25 | Pancreas overall microcirculation function evaluation and visualization method |
| CN202111410368.1A Pending CN114159026A (en) | 2021-11-09 | 2021-11-25 | Evaluation method for microcirculation function of kidney |
| CN202111410229.9A Pending CN113925470A (en) | 2021-11-09 | 2021-11-25 | Method for evaluating microcirculation function of small intestine |
Family Applications Before (2)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202111410230.1A Pending CN113925471A (en) | 2021-11-09 | 2021-11-25 | Skin microcirculation function evaluation and visualization method |
| CN202111410210.4A Pending CN113907724A (en) | 2021-11-09 | 2021-11-25 | Pancreas overall microcirculation function evaluation and visualization method |
Family Applications After (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202111410229.9A Pending CN113925470A (en) | 2021-11-09 | 2021-11-25 | Method for evaluating microcirculation function of small intestine |
Country Status (1)
| Country | Link |
|---|---|
| CN (4) | CN113925471A (en) |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20080294027A1 (en) * | 2003-06-12 | 2008-11-27 | Bracco Research Sa | System for extracting morphological information through a perfusion assessment process |
| CN105636512A (en) * | 2013-08-14 | 2016-06-01 | 南洋理工大学 | Systems and methods for revascularization assessment |
| US20180146866A1 (en) * | 2015-05-15 | 2018-05-31 | Veriskin, Inc. | Cutaneous blood flow monitoring device |
| CN111759323A (en) * | 2020-07-08 | 2020-10-13 | 刘明明 | Multi-module microcirculation function evaluation device and biological tissue microcirculation visualization method based on same |
Family Cites Families (16)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| RU2357656C1 (en) * | 2007-12-24 | 2009-06-10 | Государственное образовательное учреждение высшего профессионального образования "Астраханская государственная медицинская академия Федерального агентства по здравоохранению и социальному развитию" (ГОУ ВПО АГМА Росздрава) | Method of diagnostics of suffered chronic antenatal hypoxia in newborn infants |
| WO2010117576A2 (en) * | 2009-04-07 | 2010-10-14 | Virginia Commonwealth University | Image processing and machine learning for diagnostic analysis of microcirculation |
| US8315812B2 (en) * | 2010-08-12 | 2012-11-20 | Heartflow, Inc. | Method and system for patient-specific modeling of blood flow |
| WO2013177420A2 (en) * | 2012-05-23 | 2013-11-28 | St. Jude Children's Research Hospital | Methods and compositions for the treatment of bcr-abl positive lymphoblastic leukemias |
| CN103527184B (en) * | 2013-10-28 | 2016-04-20 | 北京大学 | A kind of Forecasting Methodology of dolostone reservoirs and system |
| US10602933B2 (en) * | 2014-10-13 | 2020-03-31 | Sergei Yurievich PODTAEV | Diagnosing disorders of microvascular tone regulation mechanisms |
| US10482241B2 (en) * | 2016-08-24 | 2019-11-19 | Sap Se | Visualization of data distributed in multiple dimensions |
| CN109215095B (en) * | 2017-07-07 | 2023-07-04 | 北京国双科技有限公司 | Data display method, device, storage medium and processor |
| CN107582097A (en) * | 2017-07-18 | 2018-01-16 | 中山大学附属第医院 | A kind of Aided intelligent decision-making learned based on multi-modal ultrasound group |
| CN108753682B (en) * | 2018-05-03 | 2019-08-02 | 中国医学科学院微循环研究所 | Promote endothelial cell into the excretion body active ingredient and its preparation method and application of blood vessel |
| CN109662735B (en) * | 2019-02-18 | 2021-07-27 | 亿慈(上海)智能科技有限公司 | Method for measuring skin blood perfusion |
| CN111796095A (en) * | 2019-04-09 | 2020-10-20 | 苏州扇贝生物科技有限公司 | Proteome mass spectrum data processing method and device |
| CN109939282A (en) * | 2019-04-23 | 2019-06-28 | 四川大学 | A kind of percutaneous Left heart assistanee system |
| CN111949842A (en) * | 2019-05-15 | 2020-11-17 | 株式会社日立制作所 | Bubble chart generation method and device |
| CN111642564A (en) * | 2020-05-11 | 2020-09-11 | 天津科技大学 | Preparation method of yoghourt |
| CN112269871A (en) * | 2020-10-12 | 2021-01-26 | 国网新疆电力有限公司信息通信公司 | Data visualization analysis method and device based on LDA topic generation model |
-
2021
- 2021-11-25 CN CN202111410230.1A patent/CN113925471A/en active Pending
- 2021-11-25 CN CN202111410210.4A patent/CN113907724A/en active Pending
- 2021-11-25 CN CN202111410368.1A patent/CN114159026A/en active Pending
- 2021-11-25 CN CN202111410229.9A patent/CN113925470A/en active Pending
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20080294027A1 (en) * | 2003-06-12 | 2008-11-27 | Bracco Research Sa | System for extracting morphological information through a perfusion assessment process |
| CN105636512A (en) * | 2013-08-14 | 2016-06-01 | 南洋理工大学 | Systems and methods for revascularization assessment |
| CN112716471A (en) * | 2013-08-14 | 2021-04-30 | 佩德拉科技私人有限公司 | System and method for assessing vascular remodeling |
| US20180146866A1 (en) * | 2015-05-15 | 2018-05-31 | Veriskin, Inc. | Cutaneous blood flow monitoring device |
| CN111759323A (en) * | 2020-07-08 | 2020-10-13 | 刘明明 | Multi-module microcirculation function evaluation device and biological tissue microcirculation visualization method based on same |
Non-Patent Citations (1)
| Title |
|---|
| 宋晓红 等: "自发性高血压大鼠肾脏和胰腺微循环血流灌注和自律运动变化", 《微循环学杂志》 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113907724A (en) | 2022-01-11 |
| CN113925471A (en) | 2022-01-14 |
| CN113925470A (en) | 2022-01-14 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| KR100794721B1 (en) | Real-time diagnostic system employing a non-invasive method to analyze the electro-magnetic field radiated from a subject and the variation thereof | |
| US9027412B2 (en) | Method and system for non-invasively monitoring fluid flow in a subject | |
| CN109259822A (en) | The guidance of focused ultrasound therapy three-dimensional and dynamic real-time monitor system and method based on the detection imaging of passive cavitation | |
| JP2010071692A (en) | Measurement apparatus and measurement method | |
| CN104840219B (en) | A kind of pen type ultrasonic color imager | |
| US10105057B2 (en) | Apparatus for optical analysis of an associated tissue | |
| Zhong et al. | A mathematical analysis on the biological zero problem in laser Doppler flowmetry | |
| CN110338843A (en) | Tissue-estimating method, apparatus, equipment and computer readable storage medium | |
| CN111759323B (en) | Multi-module microcirculation function evaluation device and biological tissue microcirculation visualization method based on same | |
| CN110361357A (en) | A kind of single array element photoacoustic spectrum signal acquisition system and method for skin detection | |
| US20050277817A1 (en) | Noninvasive measurement system for monitoring activity condition of living body | |
| Zhang et al. | Photoacoustic power azimuth spectrum for microvascular evaluation | |
| Savchenko et al. | The use of speckle correlation analysis to determine blood flow velocity | |
| US20160051149A1 (en) | Photoacoustic Probe for Burn Injury Diagnosis | |
| CN114159026A (en) | Evaluation method for microcirculation function of kidney | |
| Lindén et al. | Evaluation of Enhanced High-Resolution Laser Doppler Imaging in anin VitroTube Model with the Aim of Assessing Blood Flow in Separate Microvessels | |
| US20040234187A1 (en) | Method of and apparatus for measuring nonstationary oscillatory motion | |
| DE3838396C2 (en) | ||
| CN113786211B (en) | Ultrasonic erythrocyte aggregation degree detection and regulation system | |
| TWI588492B (en) | Near-field array detection method for detecting optically high scatter material | |
| CN101889862A (en) | Three-dimensional scanning body electrographic apparatus used for detecting human body functions and detection method thereof | |
| CN111007062B (en) | OCT real-time nondestructive monitoring method in tissue engineering skin construction process | |
| Orlova et al. | Algorithm for recognition of vascular structures in the biotissue volume | |
| US20010023391A1 (en) | Mold-in method and apparatus | |
| GB2595162A (en) | Biological function measurement device, and biological function measurement method, and program |
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 |