CN113925471A - Skin microcirculation function evaluation and visualization method - Google Patents
Skin microcirculation function evaluation and visualization method Download PDFInfo
- Publication number
- CN113925471A CN113925471A CN202111410230.1A CN202111410230A CN113925471A CN 113925471 A CN113925471 A CN 113925471A CN 202111410230 A CN202111410230 A CN 202111410230A CN 113925471 A CN113925471 A CN 113925471A
- Authority
- CN
- China
- Prior art keywords
- microcirculation
- skin
- function
- skin microcirculation
- microcirculation 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 194
- 238000011156 evaluation Methods 0.000 title claims abstract description 19
- 238000007794 visualization technique Methods 0.000 title abstract description 7
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 claims abstract description 26
- 230000017531 blood circulation Effects 0.000 claims abstract description 26
- 229910052760 oxygen Inorganic materials 0.000 claims abstract description 26
- 239000001301 oxygen Substances 0.000 claims abstract description 26
- 230000008081 blood perfusion Effects 0.000 claims abstract description 22
- 238000012800 visualization Methods 0.000 claims abstract description 18
- 230000000004 hemodynamic effect Effects 0.000 claims abstract description 16
- 238000004140 cleaning Methods 0.000 claims abstract description 13
- 238000010183 spectrum analysis Methods 0.000 claims abstract description 4
- 210000001519 tissue Anatomy 0.000 claims description 18
- 238000000034 method Methods 0.000 claims description 15
- 230000010412 perfusion Effects 0.000 claims description 12
- MWUXSHHQAYIFBG-UHFFFAOYSA-N Nitric oxide Chemical compound O=[N] MWUXSHHQAYIFBG-UHFFFAOYSA-N 0.000 claims description 10
- 238000012545 processing Methods 0.000 claims description 7
- 210000003743 erythrocyte Anatomy 0.000 claims description 6
- 238000005259 measurement Methods 0.000 claims description 6
- 210000002889 endothelial cell Anatomy 0.000 claims description 5
- 230000003287 optical effect Effects 0.000 claims description 4
- 230000000747 cardiac effect Effects 0.000 claims description 3
- 238000010276 construction Methods 0.000 claims description 3
- 238000010219 correlation analysis Methods 0.000 claims description 3
- 230000001419 dependent effect Effects 0.000 claims description 3
- 239000006185 dispersion Substances 0.000 claims description 3
- AQCHWTWZEMGIFD-UHFFFAOYSA-N metolazone Chemical compound CC1NC2=CC(Cl)=C(S(N)(=O)=O)C=C2C(=O)N1C1=CC=CC=C1C AQCHWTWZEMGIFD-UHFFFAOYSA-N 0.000 claims description 3
- 210000003205 muscle Anatomy 0.000 claims description 3
- 210000005036 nerve Anatomy 0.000 claims description 3
- 230000000241 respiratory effect Effects 0.000 claims description 3
- 239000000523 sample Substances 0.000 claims description 3
- 238000011425 standardization method Methods 0.000 claims description 3
- 238000010586 diagram Methods 0.000 abstract description 13
- 238000004458 analytical method Methods 0.000 abstract description 4
- 230000000007 visual effect Effects 0.000 abstract description 2
- 210000000056 organ Anatomy 0.000 description 6
- 230000008859 change Effects 0.000 description 4
- 201000010099 disease Diseases 0.000 description 4
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 4
- 230000002159 abnormal effect Effects 0.000 description 3
- 210000000115 thoracic cavity Anatomy 0.000 description 3
- 210000000683 abdominal cavity Anatomy 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 2
- 238000003759 clinical diagnosis Methods 0.000 description 2
- 230000003902 lesion Effects 0.000 description 2
- 230000000877 morphologic effect Effects 0.000 description 2
- 238000001228 spectrum Methods 0.000 description 2
- 239000002344 surface layer Substances 0.000 description 2
- 230000009885 systemic effect Effects 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 230000003187 abdominal effect Effects 0.000 description 1
- 210000003484 anatomy Anatomy 0.000 description 1
- 238000009412 basement excavation Methods 0.000 description 1
- 238000007405 data analysis Methods 0.000 description 1
- 238000013481 data capture Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000004064 dysfunction Effects 0.000 description 1
- 230000003511 endothelial effect Effects 0.000 description 1
- 210000003734 kidney Anatomy 0.000 description 1
- 210000000265 leukocyte Anatomy 0.000 description 1
- 210000004072 lung Anatomy 0.000 description 1
- 238000013508 migration Methods 0.000 description 1
- 230000005012 migration Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000009257 reactivity Effects 0.000 description 1
- 230000008326 skin blood flow Effects 0.000 description 1
- 210000000813 small intestine Anatomy 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 230000001457 vasomotor Effects 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)
- Measuring And Recording Apparatus For Diagnosis (AREA)
- Investigating Or Analysing Biological Materials (AREA)
Abstract
The invention provides a skin microcirculation function evaluation and visualization method, which comprises the steps of measuring skin microcirculation oxygen parameters and microcirculation hemodynamic parameters, cleaning skin microcirculation function data, performing visualization engineering on skin microcirculation function, constructing a three-dimensional visualization model, a skin microcirculation function uniaxial bubble diagram and performing wavelet transform spectrum analysis. The invention realizes accurate and comprehensive description of the skin microcirculation function characteristics through the capture and analysis of skin microcirculation hemodynamics (microcirculation blood perfusion and blood flow velocity) and microcirculation oxygen (microcirculation oxygen partial pressure). Meanwhile, the invention realizes the visual engineering of the skin microcirculation function through the universal microcirculation framework and improves the scientificity of the feature display of the skin microcirculation function.
Description
Technical Field
The invention belongs to the technical field of microcirculation data analysis, and particularly relates to a skin microcirculation function evaluation and visualization method.
Background
The skin covers the whole body, is superficial in position, and is the largest organ in surface area in the human body. The reactive change of skin microvascular is often appeared in the early stage of certain diseases, can reflect the function of the whole body microvascular to a certain extent, and is an ideal site for evaluating the dysfunction of the microvascular. The observation of clinical living skin microcirculation can provide the change index of microvascular for systemic diseases and local lesions, and is helpful for clinical diagnosis and treatment guidance.
As shown in fig. 1, the commonly used skin microvascular evaluation parameters are mainly divided into two major categories, namely morphological indexes and functional indexes, and the morphological evaluation parameters include microvascular density, vessel diameter, vasomotor motion, blood flow velocity, erythrocyte and leukocyte migration, and the like; the functional evaluation parameters comprise microvascular pressure, elasticity, release of endothelial derived vaso-and vaso-factors, skin temperature, local skin blood flow perfusion and the like.
The existing methods and techniques for skin microcirculation function suffer from the following objective disadvantages:
(1) index stability. The capillary microscope commonly used in clinic can visually observe the state of the capillary vessels and the blood flow condition on the surface layer of the skin, but the observation result is easily influenced by local and external factors and is limited by the magnification of the microscope device;
(2) index singleness and scientificity. The prior art can only measure a single index of skin microcirculation function, and cannot realize the comprehensive attention to skin microcirculation blood perfusion and skin microcirculation oxygen indexes. Furthermore, there are mostly abnormal values in cutaneous microvascular perfusion and blood flow velocity measured by the prior art. The prior art does not have strategies such as abnormal value processing, boundary value adjustment and the like, and the objective defects further influence the scientificity and accuracy of the evaluation of the skin microcirculation function.
Therefore, how to realize scientific evaluation of skin microcirculation function is a problem to be solved.
Disclosure of Invention
Aiming at the problems in the prior art, the invention aims to provide a skin microcirculation function evaluation and visualization method, which realizes accurate and comprehensive description of skin microcirculation function characteristics through capture and analysis of skin microcirculation hemodynamics (microcirculation blood perfusion and blood flow velocity) and microcirculation oxygen (microcirculation oxygen partial pressure). Meanwhile, the invention realizes the visual engineering of the skin microcirculation function through the universal microcirculation framework and improves the scientificity of the feature display of the skin microcirculation function.
The purpose of the invention is realized by the following technical scheme:
a skin microcirculation function evaluation and visualization method comprises the following steps:
step 1, measurement of skin microcirculation oxygen parameter and microcirculation hemodynamics parameter:
measuring skin microcirculation function parameter data including skin microcirculation oxygen parameters and hemodynamic parameters by using a skin microcirculation function evaluation device, wherein the skin microcirculation function evaluation device is composed of a microcirculation oxygen partial pressure monitor Microx TX3 and a dual-channel laser Doppler monitor VMS, guiding an optical probe to a position 1mm above the skin during measurement, and using two light sources of white light and laser;
importing the skin microcirculation function parameter data measured in the step 1 into a data cleaning module, and processing the outlier skin microcirculation function parameter data through a boxplot algorithm;
step 3, visualization engineering of skin microcirculation function:
step 3.1, three-dimensional visualization model construction
Importing the skin microcirculation function parameter data after the cleaning treatment in the step 2 for dimensionless treatment: processing the skin microcirculation function parameter data by a dispersion standardization method, eliminating the dimension of the multi-parameter microcirculation function data, uniformly projecting the dimension in a [0, 1] interval, and realizing the visualization of the skin microcirculation function under the same coordinate system frame;
step 3.2, skin microcirculation function uniaxial bubble chart
Importing the skin microcirculation function parameter data after the cleaning treatment in the step 2, and generating a skin integral microcirculation function uniaxial bubble chart by using EChats to realize skin microcirculation function parameter distribution and weight correlation analysis;
step 3.3, wavelet transform spectrum analysis:
further analyzing and visualizing the skin microcirculation blood flow signal by using wavelet transform, importing the hemodynamic parameter data cleaned and processed in the step 2, dividing the hemodynamic 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 skin microcirculation blood flow perfusion, blood flow velocity and erythrocyte density signals in a characteristic frequency band through the n frequency band source signal amplitudes; drawing a three-dimensional time-frequency diagram of the skin microcirculation blood flow perfusion signal through three dimensional indexes of time (sec), frequency (Hz) and microcirculation characteristic source Amplitude (AU); the distribution pattern of the amplitude related to the function of the skin microcirculation is shown by radar map.
Further, the skin microcirculation oxygen parameter includes skin microcirculation tissue oxygen partial Pressure (PO)2) The hemodynamic parameters include a level of cutaneous microcirculation Blood Perfusion (BP) and a speed of cutaneous microcirculation blood perfusion (V).
Further, when treating the outlier skin microcirculation function parameter in step 2, defining Q125% maximum, Q3Is 75% of maximum, Q3And Q1The difference between the values is the four-bit distance IQR, and (Q) is set1-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, step 3.1 uses Python and ECharts to generate a skin microcirculation function three-dimensional visualization module, wherein time, skin microcirculation function parameter variables and microcirculation function parameter variable values are respectively defined as the module X axis, Y axis and Z axis.
Further, the amplitude values related to the skin microcirculation function in step 3.3 include cardiac source amplitude values, respiratory source amplitude values, muscle source amplitude values, nerve source amplitude values, nitric oxide dependent endothelial cell source amplitude values and nitric oxide independent endothelial cell source amplitude values.
Compared with the prior art, the invention has the beneficial effects that:
1. compared with the prior art that the capillary microscope is used for directly observing the state and the blood flow condition of the skin surface layer capillary, the skin microcirculation function in the scheme of the invention comprises multiple dimensions of skin microcirculation blood perfusion, blood flow speed and skin microcirculation oxygen partial pressure, and simultaneously breaks through the limitation of the magnification of a microscope; in addition, compared with the microcirculation function index in the prior art, the method can evaluate the skin microcirculation function more stably and comprehensively;
2. in the scheme of the invention, the skin microcirculation function eliminates abnormal data from instrument firmware and optical components through an algorithm, so that the scientificity of the skin microcirculation function data is improved; in addition, the scheme of the invention can visualize the skin microcirculation function data by combining a universal microcirculation framework and a single-axis bubble diagram, a two-bit frequency spectrum diagram and a three-dimensional time-frequency diagram, and is beneficial to deep analysis and excavation of the skin microcirculation function state; the invention provides a brand new method for comprehensively and scientifically evaluating the microcirculation function of the skin.
Drawings
The invention is further illustrated by the following examples in conjunction with the accompanying drawings:
FIG. 1 is a schematic diagram of the evaluation parameters of skin microcirculation in the prior art;
FIG. 2 is a schematic flow chart of a method for evaluating and visualizing the microcirculation function of skin according to the present invention;
FIG. 3 is a schematic diagram of the skin microcirculation function evaluation and visualization method according to the present invention;
FIG. 4 is a functional three-dimensional visualization view of skin microcirculation in example 1; in the figure, (a) a three-dimensional visualization histogram of skin microcirculation, (b) a three-dimensional visualization scatter diagram of skin microcirculation;PO2: microcirculation tissue oxygen partial pressure, BP: microcirculation blood perfusion level, V: perfusion rate of microcirculation blood flow;
FIG. 5 is a graph of uniaxial bubbles produced using test data;
FIG. 6 is a diagram of the skin microcirculation function uniaxial air bubbles of 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 diagram of the skin microvascular perfusion, blood flow velocity and red blood cell density signals at a characteristic frequency band; the black curve is the blood perfusion of skin microvascular blood flow, the green curve is the microcirculation blood flow velocity, the red curve is erythrocyte density signals, and the spectrum distribution diagram of the microcirculation blood flow signals represents the frequency interval of the amplitude and the characteristic peak of each frequency band;
FIG. 8 is a three-dimensional time-frequency diagram of a skin microcirculation blood perfusion signal; in the figure, the scale of blue (low) to orange (high) represents the continuous wavelet coefficients;
FIG. 9 is a radar plot of amplitude related to skin microcirculation function.
Detailed Description
Example 1
As shown in fig. 2 and fig. 3, the embodiment provides a method for evaluating and visualizing a skin microcirculation function, the method relies on a computer algorithm to evaluate the skin microcirculation function, and relies on a general microcirculation framework to establish a method for visualizing the skin microcirculation function, compared with other tissues:
1) the skin tissue is positioned on the superficial surface of the organism, compared with the microcirculation observation of thoracic cavity organs and abdominal cavity organs and tissues, the microcirculation function of the living skin is easier to observe, so the skin tissue is an ideal site for clinically evaluating the functional state of the microvasculature of the organism;
2) the reactivity and the functional change of skin microvascular are often shown in the early stage of the disease, so the microcirculation-level index can be provided for the systemic disease and the local lesion, and the clinical diagnosis and the guidance treatment are facilitated;
3) because the skin tissue is shallow in position, the microcirculation function data capture access and the light source path are different from thoracic organs and tissues (such as lungs) and abdominal organs and tissues (such as small intestines, kidneys and the like);
4) because of the uniqueness of the anatomical structure, the interference of skin tissue microcirculation function data and the source of the interference are different from those of organs and tissues of the abdominal cavity and the thoracic cavity.
Specifically, the method for visualizing the microcirculation function of the skin comprises the following steps:
step 1, measurement of skin microcirculation oxygen parameter and microcirculation hemodynamics parameter:
measuring skin microcirculation function parameter data including skin microcirculation tissue oxygen partial Pressure (PO) by using skin microcirculation function evaluation device2) The skin microcirculation function evaluation device comprises a microcirculation oxygen partial pressure monitor Microx TX3 (Presens company in Germany) and a two-channel laser Doppler monitor VMS (Moor company in UK), an optical probe is guided to a position 1mm above the skin during measurement, and two light sources of white light and laser are used.
and (3) importing the skin microcirculation function parameter data captured in the step (1) into a data cleaning module, and processing the outlier microcirculation function parameter data through a boxplot 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 defined as outlier and adjusted to the normal range boundary value.
Step 3, visualization engineering of skin microcirculation function:
step 3.1, three-dimensional visualization model construction
Importing the skin microcirculation function parameter data after the cleaning treatment in the step 2 for dimensionless treatment: processing the skin microcirculation function parameter data by a dispersion standardization method, eliminating the multi-parameter microcirculation function data dimension and uniformly projecting the multi-parameter microcirculation function data dimension on [0, 1]]Within the interval, realize the same coordinate system frameVisualization of subclavian skin microcirculation function. Specifically, the skin microcirculation functional parameter data after dimensionless treatment, namely the oxygen partial Pressure (PO) of skin microcirculation tissue is imported2) Python and ECharts are used for generating a skin microcirculation function three-dimensional visualization module, wherein time, microcirculation function parameter variables and microcirculation function parameter variable values are respectively defined as an X axis, a Y axis and a Z axis of the module, and a skin microcirculation function three-dimensional visualization view shown in figure 4 is obtained.
Step 3.2, skin microcirculation function uniaxial bubble chart
And (3) importing the skin microcirculation function parameter data after the cleaning treatment in the step (2), and generating a skin integral microcirculation function uniaxial bubble chart by using EChats to realize skin microcirculation function parameter distribution and weight correlation analysis.
In the embodiment, a uniaxial bubble chart as shown in fig. 5 is generated by using part of test data, wherein the horizontal axis in the bubble chart represents the value range of the data distribution of the skin microcirculation function parameter, the value range is divided into 7 intervals (interval 0-2, interval 2-4, interval 4-6, interval 6-8, interval 8-10, interval 10-12 and interval 12-14), and the distribution weight correlation of the data of the skin microcirculation function parameter is reflected; 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, the more the microcirculation function data fall in the corresponding section of the circle. Therefore, the distribution and weight association condition of the skin microcirculation function data in each value interval can be visually shown in the uniaxial bubble chart.
The skin microcirculation tissue oxygen partial Pressure (PO) as shown in FIG. 62) In the distribution weight correlation bubble chart, the circular area corresponding to the interval scale 50 hPa is larger, and the circular area corresponding to the interval 0-40 hPa is smaller, which indicates that more skin microcirculation tissue oxygen partial pressure data are intensively distributed on the interval scale 50 hPa, and less data are distributed on the interval scale 0-40 hPa, namely the skin microcirculation tissue oxygen partial pressure is mainly at the level of 50 hPa, and the microcirculation blood flow perfusion speed at the level of 50 hPa is the dominant weight of the microcirculation function.
Step 3.3, wavelet transform spectrum analysis:
the present embodiment uses wavelet transformation for further analysis and visualization of skin microcirculation blood flow signals. The hemodynamic parameter data (skin microcirculation blood perfusion level (BP) and skin microcirculation blood perfusion velocity (V)) processed in step 2 are imported, and as shown in table 1, the hemodynamic parameter data can be divided into six frequency band source signal amplitudes according to the microcirculation blood perfusion signal frequency.
Through the six frequency band source signal amplitudes, a two-dimensional distribution diagram of the skin microvascular blood perfusion, blood flow velocity and red blood cell density signals in the characteristic frequency band is generated as shown in fig. 7.
A three-dimensional time-frequency graph of the skin microcirculation blood perfusion signal shown in figure 8 is drawn through three-dimensional indexes of time (sec), frequency (Hz) and microcirculation characteristic source Amplitude (AU), and represents distribution characteristics and change rules of 6 source amplitudes related to skin microcirculation functions along with time progress.
By means of a radar map as shown in fig. 9, six patterns of distribution of amplitude related to the function of the skin microcirculation are shown, including cardiac source amplitude, respiratory source amplitude, muscle source amplitude, nerve source amplitude, nitric oxide dependent and independent endothelial cell source amplitude.
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 (5)
1. A method for assessing and visualizing the microcirculation function of the skin, said method comprising the following steps:
step 1, measurement of skin microcirculation oxygen parameter and microcirculation hemodynamics parameter:
measuring skin microcirculation function parameter data including skin microcirculation oxygen parameters and hemodynamic parameters by using a skin microcirculation function evaluation device, wherein the skin microcirculation function evaluation device is composed of a microcirculation oxygen partial pressure monitor Microx TX3 and a dual-channel laser Doppler monitor VMS, guiding an optical probe to a position 1mm above the skin during measurement, and using two light sources of white light and laser;
step 2, cleaning skin microcirculation function data:
importing the skin microcirculation function parameter data measured in the step 1 into a data cleaning module, and processing the outlier skin microcirculation function parameter data through a boxplot algorithm;
step 3, visualization engineering of skin microcirculation function:
step 3.1, three-dimensional visualization model construction
Importing the skin microcirculation function parameter data after the cleaning treatment in the step 2 for dimensionless treatment: processing the skin microcirculation function parameter data by a dispersion standardization method, eliminating the dimension of the multi-parameter microcirculation function data, uniformly projecting the dimension in a [0, 1] interval, and realizing the visualization of the skin microcirculation function under the same coordinate system frame;
step 3.2, skin microcirculation function uniaxial bubble chart
Importing the skin microcirculation function parameter data after the cleaning treatment in the step 2, and generating a skin integral microcirculation function uniaxial bubble chart by using EChats to realize skin microcirculation function parameter distribution and weight correlation analysis;
step 3.3, wavelet transform spectrum analysis:
further analyzing and visualizing the skin microcirculation blood flow signal by using wavelet transform, importing the hemodynamic parameter data cleaned and processed in the step 2, dividing the hemodynamic 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 skin microcirculation blood flow perfusion, blood flow velocity and 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 skin microcirculation blood flow perfusion signal through three dimensional indexes of time, frequency and microcirculation characteristic source amplitude; the distribution pattern of the amplitude related to the function of the skin microcirculation is shown by radar map.
2. The method of assessing and visualizing skin microcirculation function of claim 1, wherein said skin microcirculation oxygen parameter of step 1 includes skin microcirculation tissue oxygen partial Pressure (PO)2) The hemodynamic parameters include a level of cutaneous microcirculation Blood Perfusion (BP) and a speed of cutaneous microcirculation blood perfusion (V).
3. The method for assessing and visualizing microcirculation function of skin according to claim 1, wherein Q is defined when treating the parameter of skin microcirculation function of outliers in step 2125% maximum, Q3Is 75% of maximum, Q3And Q1The difference between the values is the four-bit distance IQR, and (Q) is set1-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 of claim 1, wherein step 3.1 is performed by using Python and ECharts to generate a three-dimensional visualization module of skin microcirculation function, wherein time, variables of skin microcirculation function parameters and variables of microcirculation function parameters are respectively defined as X-axis, Y-axis and Z-axis of the module.
5. The method of claim 1, wherein the skin microcirculation function related amplitudes in step 3.3 include cardiac source amplitude, respiratory source amplitude, muscle source amplitude, nerve source amplitude, nitric oxide dependent endothelial cell source amplitude and nitric oxide independent endothelial cell source amplitude.
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 |
|---|---|
| CN113925471A true CN113925471A (en) | 2022-01-14 |
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 |
| 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 |
| CN202111410210.4A Pending CN113907724A (en) | 2021-11-09 | 2021-11-25 | Pancreas overall microcirculation function evaluation and visualization method |
Family Applications After (3)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| 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 |
| CN202111410210.4A Pending CN113907724A (en) | 2021-11-09 | 2021-11-25 | Pancreas overall microcirculation function evaluation and visualization method |
Country Status (1)
| Country | Link |
|---|---|
| CN (4) | CN113925471A (en) |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20120269420A1 (en) * | 2009-04-07 | 2012-10-25 | Kayvan Najarian | Image Processing and Machine Learning for Diagnostic Analysis of Microcirculation |
| AU2017221811A1 (en) * | 2010-08-12 | 2017-09-21 | Heartflow, Inc. | Method and system for patient-specific modeling of blood flow |
| CN109662735A (en) * | 2019-02-18 | 2019-04-23 | 亿慈(上海)智能科技有限公司 | The measurement method of SkBF groundwater increment |
| 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 |
|---|---|---|---|---|
| US8021303B2 (en) * | 2003-06-12 | 2011-09-20 | Bracco Research Sa | System for extracting morphological information through a perfusion assessment process |
| RU2357656C1 (en) * | 2007-12-24 | 2009-06-10 | Государственное образовательное учреждение высшего профессионального образования "Астраханская государственная медицинская академия Федерального агентства по здравоохранению и социальному развитию" (ГОУ ВПО АГМА Росздрава) | Method of diagnostics of suffered chronic antenatal hypoxia in newborn infants |
| 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 |
| CN105636512B (en) * | 2013-08-14 | 2020-11-24 | 佩德拉科技私人有限公司 | System and method for assessing vascular remodeling |
| 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 |
| WO2016187136A1 (en) * | 2015-05-15 | 2016-11-24 | Veriskin, Inc. | Cutaneous blood flow monitoring device |
| 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 | 中山大学附属第医院 | Intelligent aid decision-making system based on multi-mode ultrasound omics |
| 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 |
| 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 CN202111410368.1A patent/CN114159026A/en active Pending
- 2021-11-25 CN CN202111410229.9A patent/CN113925470A/en active Pending
- 2021-11-25 CN CN202111410210.4A patent/CN113907724A/en active Pending
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20120269420A1 (en) * | 2009-04-07 | 2012-10-25 | Kayvan Najarian | Image Processing and Machine Learning for Diagnostic Analysis of Microcirculation |
| AU2017221811A1 (en) * | 2010-08-12 | 2017-09-21 | Heartflow, Inc. | Method and system for patient-specific modeling of blood flow |
| CN109662735A (en) * | 2019-02-18 | 2019-04-23 | 亿慈(上海)智能科技有限公司 | The measurement method of SkBF groundwater increment |
| 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 |
| CN113925470A (en) | 2022-01-14 |
| CN114159026A (en) | 2022-03-11 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US11963819B2 (en) | Method and system for ultrasonic fluid spectral doppler imaging | |
| Zhong et al. | A mathematical analysis on the biological zero problem in laser Doppler flowmetry | |
| WO2022213938A1 (en) | Ultrasound elastography method and apparatus, and electronic device and storage medium | |
| CN111000544B (en) | Method and system for constructing hybrid continuous blood pressure measurement model based on PPG waveform | |
| CN113081038B (en) | Elastography method, elastography device, electronic equipment and storage medium | |
| Line et al. | Repeated measurement variation and precision of laser Doppler flowmetry measurements | |
| CN111759323B (en) | Multi-module microcirculation function evaluation device and biological tissue microcirculation visualization method based on same | |
| CN113133753A (en) | Biological tissue blood flow real-time monitoring system and simulation monitoring system based on magnetic induction phase shift | |
| US20050277817A1 (en) | Noninvasive measurement system for monitoring activity condition of living body | |
| Tazawa et al. | Noncontact measurements of avian embryo heart rate by means of the laser speckle: comparison with contact measurements | |
| CN113925471A (en) | Skin microcirculation function evaluation and visualization method | |
| CN107478414B (en) | OCT imaging back-pumping performance testing device and method | |
| RU2692220C1 (en) | Method of colour doppler mapping in endoscopic optical coherence tomography | |
| 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 | |
| CN117953145A (en) | Method for visualizing dynamic quantification of three-dimensional biological sample cells | |
| Wojtkiewicz et al. | Evaluation of algorithms for microperfusion assessment by fast simulations of laser Doppler power spectral density | |
| CN113951850A (en) | Lung tissue overall microcirculation function evaluation and visual presentation method | |
| CN109589140A (en) | A kind of ultrasonic measurement entry processing method and compuscan | |
| Rogatkin et al. | Optical Non-invasive Flowmetry without Lasers and Coherent Light. | |
| GB2595162A (en) | Biological function measurement device, and biological function measurement method, and program | |
| RU2687775C1 (en) | Device for diagnosing organ pathology | |
| Blair et al. | Transcutaneous blood flow measurements in arteries of the human hand | |
| CN202446057U (en) | Doppler laser OCT cystoscope system | |
| Belsheva et al. | Sensitivity of Frequency Domain Near Infrared Spectroscopy for Neurovascular Structure Detection in Biotissue Volume: Numerical Modeling Results | |
| McCulloch et al. | Machine Controls and Optimizing the Images |
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 |