CN108805949A - A kind of modification method of blood flow velocity, device, terminal and computer storage media - Google Patents

A kind of modification method of blood flow velocity, device, terminal and computer storage media Download PDF

Info

Publication number
CN108805949A
CN108805949A CN201810517956.7A CN201810517956A CN108805949A CN 108805949 A CN108805949 A CN 108805949A CN 201810517956 A CN201810517956 A CN 201810517956A CN 108805949 A CN108805949 A CN 108805949A
Authority
CN
China
Prior art keywords
error
discrete
equation
speed
velocity component
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201810517956.7A
Other languages
Chinese (zh)
Other versions
CN108805949B (en
Inventor
魏润杰
高琪
刘星利
吴鹏
李飞
王巍
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Hangzhou Sheng Shi Technology Co Ltd
Original Assignee
Hangzhou Sheng Shi Technology Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Hangzhou Sheng Shi Technology Co Ltd filed Critical Hangzhou Sheng Shi Technology Co Ltd
Priority to CN201810517956.7A priority Critical patent/CN108805949B/en
Publication of CN108805949A publication Critical patent/CN108805949A/en
Application granted granted Critical
Publication of CN108805949B publication Critical patent/CN108805949B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00—Image analysis
    • G06T7/20—Analysis of motion
    • G06T7/262—Analysis of motion using transform domain methods, e.g. Fourier domain methods
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/14—Fourier, Walsh or analogous domain transformations, e.g. Laplace, Hilbert, Karhunen-Loeve, transforms
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T11/00—2D [Two Dimensional] image generation
    • G06T11/003—Reconstruction from projections, e.g. tomography
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/10—Image acquisition modality
    • G06T2207/10072—Tomographic images
    • G06T2207/10088—Magnetic resonance imaging [MRI]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/30—Subject of image; Context of image processing
    • G06T2207/30004—Biomedical image processing
    • G06T2207/30101—Blood vessel; Artery; Vein; Vascular
    • G06T2207/30104—Vascular flow; Blood flow; Perfusion

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Theoretical Computer Science (AREA)
  • Pure & Applied Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • General Engineering & Computer Science (AREA)
  • Algebra (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Multimedia (AREA)
  • Operations Research (AREA)
  • Computing Systems (AREA)
  • Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)
  • Ultra Sonic Daignosis Equipment (AREA)

Abstract

The embodiment of the invention discloses a kind of modification method of blood flow velocity, this method includes：Obtain the speed field areas of blood flow, determine the convection equation of each error velocity component of blood flow on grid node, the convection equation of each error velocity component is carried out respectively discrete, obtain the discrete form of convection equation, determine the continuity equation of the error speed of blood flow on grid node, continuity equation is carried out discrete, obtain the discrete form of continuity equation, using the discrete form of the discrete form of convection equation and continuity equation as constraint equation, determine the minimum value of goal-selling function, by the corresponding each error velocity component of minimum value, it is determined as each error velocity component of blood flow on the grid node of speed field areas, according to each error velocity component of blood flow, each velocity component of blood flow to measuring is modified.The embodiment of the present invention further simultaneously discloses a kind of correcting device of blood flow velocity, terminal and computer storage media.

Description

A kind of modification method of blood flow velocity, device, terminal and computer storage media
Technical field
The present invention relates to the correction technique of blood flow velocity in blood flow measurement more particularly to a kind of amendment sides of blood flow velocity Method, device, terminal and computer storage media.
Background technology
In clinical treatment diagnosis, the imaging of color ultrasound angiocarpy, Magnetic resonance imaging (MRI, Magnetic Resonance Imaging) and Technology of Digital Subtraction Angiography (DSA, Digital Subtraction Angiography) etc. is now Clinically main medical imaging system can not only provide human organ image, can also pass through the Hemodynamic environment of time series The variation that haemodynamics is embodied in field is spent, although the development of these technologies greatly improves the efficiency and accuracy of diagnosis, Be image noise presence to velocity field extraction cause prodigious interference；Such as the speed to contrast agent flow in two-dimentional DSA Reduction, flowing (MRI 4D Flow, the Magnetic Resonance that Magnetic resonance imaging changes over time on three dimensions Imaging 4D Flow) in 3D velocity field extraction etc., wherein noise so that image is fuzzy, or even floods image, shows For with randomness, between image generally with correlation and with additivity, it is seen then that the presence of noise seriously affects The quality of image, and the correlation with the time is shown in time series data, therefore, for blood flow velocity field, in velocity field The existing error speed with time correlation, which is modified, becomes particularly important.
Image particle velocity measuring technique (PIV, Particle Image Velocimetry) is experiment Hemodynamics field One of most important quantitative measurment means, the technology obtain measurement zone by analyzing imaging of the flow field trace particle on camera Velocity field in domain, since the velocity vector of PIV is obtained by the cross-correlation analysis of two different moments, particle The movement of image has temporal correlation in different moments.When particle imaging because light intensity is excessive overexposure occurs, because of optics Aperture is too small to there is diffraction or because when the factors such as laser polarization form particle picture intensity profile specific form, this A little image error information for being unfavorable for accurately obtaining velocity vector can have a correlation in the imaging of different moments, and this and when Between relevant error speed can be embodied in velocity field, the spatially-correlated errors of it and velocity field and random error three deposit jointly It is in speed field error.
The processing for the imaging of color ultrasound angiocarpy, MRI, DSA and experiment measuring speed field is mainly filtering and noise reduction at present Method, such method is simple and effective, but the problem of bringing is quite a few, and mainly filtering and noise reduction method can draw to some extent Enter non-physical error, secondary pollution is caused to blood flow velocity field, and the filtering and noise reduction method based on statistical property is to random The inhibiting effect of error is unfavorable for handling the error of space correlation and time correlation, by three kinds of error solutions than more significant It is very difficult that coupling, which carries out processing,.
However, at present since the difficulty for obtaining the velocity field based on time resolution is larger, so not over temporal information The technology of error inhibition is directly carried out, and Eulerian configuration is all based on generally directed to the analysis and post-processing of velocity field, It is not related to temporal information, it can thus be seen that existing without for the modification method with the error speed of time correlation.
Invention content
In view of this, an embodiment of the present invention is intended to provide a kind of modification method of blood flow velocity, device, terminal organ and calculating Machine storage medium, it is intended to realize for the amendment with the error speed of time correlation.
In order to achieve the above objectives, the technical proposal of the invention is realized in this way：
In a first aspect, the embodiment of the present invention provides a kind of modification method of blood flow velocity, including：Obtain the velocity field of blood flow Region；Wherein, grid node is distributed in the speed field areas, determines blood flow on the grid node of the speed field areas The convection equation of each error velocity component, it is discrete to the progress of the convection equation of each error velocity component respectively, obtain each error speed The discrete form for spending the convection equation of component determines the continuous of the error speed of blood flow on the grid node of the speed field areas Property equation, the continuity equation of the error speed is carried out discrete, obtains the discrete form of the continuity equation of error speed, With the discrete form of the discrete form of the convection equation of each error velocity component and the continuity equation of the error speed For constraint equation, the minimum value of goal-selling function is determined, by the corresponding each error velocity component of the minimum value, be determined as institute Each error velocity component for stating blood flow on the grid node of speed field areas, it is right according to each error velocity component of the blood flow Each velocity component of blood flow is modified on the grid node of the speed field areas measured.
Further, described discrete to the progress of the convection equation of each error velocity component respectively, obtain each error speed point The discrete form of the convection equation of amount, including：It is several to the local derviation of time in the convection equation of each error velocity component carry out from It dissipates, and discrete to the several progress of the local derviation of each coordinate direction of grid node, obtains the convection equation of each error velocity component Discrete form.
Further, the several progress of the local derviation of time are discrete and right in the convection equation to each error velocity component The several progress of local derviation of each coordinate direction of grid node are discrete, obtain the discrete shape of the convection equation of each error velocity component Formula, including：According to preset time point, progress several to the local derviation of time in the convection equation of each error velocity component are discrete, Use fourth central difference scheme for the node inside grid, for the node in net boundary using forward difference format or Person's backward difference format, progress several to the local derviation of each coordinate direction in the convection equation of each error velocity component are discrete, obtain The discrete form of the convection equation of each error velocity component.
Further, the continuity equation to the error speed carries out discrete, obtains the continuity of error speed The discrete form of equation, including：Progress several to the local derviation of each error velocity component in the continuity equation of error speed are discrete, And it is discrete to the several progress of the local derviation of each measuring speed component, obtain the discrete form of the continuity equation of the error speed.
Further, in the continuity equation to error speed the local derviation of each error velocity component it is several carry out from It dissipates, and discrete to the several progress of the local derviation of each measuring speed, obtains the discrete form of the continuity equation of the error speed, wrap It includes：Fourth central difference scheme is used for the node inside grid, forward difference lattice are used for the node in net boundary Formula or backward difference format, progress several to the local derviation of each error velocity component in the continuity equation of error speed are discrete, And it is discrete to the several progress of the local derviation of each measuring speed component, obtain the discrete form of the continuity equation of the error speed.
Further, using each error velocity component as independent variable, accordingly, the goal-selling function is：It is each to become certainly The evolution of the quadratic sum of two norms of amount.
Second aspect, the embodiment of the present invention provide a kind of correcting device of blood flow velocity, including：Acquisition module, for obtaining Take the speed field areas of blood flow；Wherein, grid node, the first determining module, for determining is distributed in the speed field areas The convection equation of each error velocity component of blood flow on the grid node of the speed field areas, the first discrete block, for dividing The other convection equation progress to each error velocity component is discrete, obtains the discrete form of the convection equation of each error velocity component, Second determining module, the continuity equation of the error speed of blood flow on the grid node for determining the speed field areas, Two discrete blocks, it is discrete for the continuity equation progress to the error speed, obtain the continuity equation of error speed Discrete form, third determining module, for the discrete form of the convection equation of each error velocity component and the error The discrete form of the continuity equation of speed is constraint equation, the minimum value of goal-selling function is determined, by the minimum value pair Each error velocity component answered is determined as each error velocity component of blood flow on the grid node of the speed field areas, corrects Module, for each error velocity component according to the blood flow, to blood on the grid node of the speed field areas measured Each velocity component of stream is modified.
Further, first discrete block is specifically used for：To the time in the convection equation of each error velocity component The several progress of local derviation are discrete, and discrete to the several progress of the local derviation of each coordinate direction of grid node, obtain each error speed The discrete form of the convection equation of component.
Further, first discrete block it is several to the local derviation of time in the convection equation of each error velocity component into Row is discrete, and discrete to the several progress of the local derviation of each coordinate direction of grid node, obtains the convection current of each error velocity component The discrete form of equation, including：According to preset time point, to the partial derivative of time in the convection equation of each error velocity component Item progress is discrete, fourth central difference scheme is used for the node inside grid, before being used for the node in net boundary To difference scheme or backward difference format, it is several to the local derviation of each coordinate direction in the convection equation of each error velocity component into Row is discrete, obtains the discrete form of the convection equation of each error velocity component.
Further, second discrete block is specifically used for：To each error speed in the continuity equation of error speed The several progress of local derviation of component are discrete, and discrete to the several progress of the local derviation of each measuring speed component, obtain the error speed Continuity equation discrete form.
Further, local derviation of second discrete block to each error velocity component in the continuity equation of error speed Several progress are discrete, and discrete to the several progress of the local derviation of each measuring speed component, obtain the continuity side of the error speed The discrete form of journey, including：Fourth central difference scheme is used for the node inside grid, for the node in net boundary Using forward difference format or backward difference format, to the local derviation of each error velocity component in the continuity equation of error speed Several progress are discrete, and discrete to the several progress of the local derviation of each measuring speed component, obtain the continuity side of the error speed The discrete form of journey.
Further, using each error velocity component as independent variable, accordingly, the goal-selling function is：It is each to become certainly The evolution of the quadratic sum of two norms of amount.
The third aspect, the embodiment of the present invention provide a kind of terminal, and the terminal includes that processor, memory and communication are total Line；The communication bus is for realizing the connection communication between processor and memory；The processor is for executing memory The revision program of the blood flow velocity of middle storage, to realize following steps：
Obtain the speed field areas of blood flow；Wherein, grid node is distributed in the speed field areas, determines the speed The convection equation of each error velocity component of blood flow on the grid node of field areas, respectively to the convection current side of each error velocity component Cheng Jinhang is discrete, obtains the discrete form of the convection equation of each error velocity component, determines the grid section of the speed field areas The continuity equation of the error speed of blood flow on point, it is discrete to the continuity equation progress of the error speed, obtain error speed The discrete form of the continuity equation of degree, with the discrete form of the convection equation of each error velocity component and error speed The discrete form of the continuity equation of degree is constraint equation, determines the minimum value of goal-selling function, the minimum value is corresponded to Each error velocity component, be determined as each error velocity component of blood flow on the grid node of the speed field areas, according to institute Each error velocity component for stating blood flow, to each velocity component of blood flow on the grid node of the speed field areas measured into Row is corrected.
Further, described discrete to the progress of the convection equation of each error velocity component respectively, obtain each error speed point In the discrete form of the convection equation of amount, the processor is used to execute the revision program of the error speed stored in memory, To realize following steps：Progress several to the local derviation of time in the convection equation of each error velocity component are discrete, and to grid section The several progress of local derviation of each coordinate direction of point are discrete, obtain the discrete form of the convection equation of each error velocity component.
Further, the several progress of the local derviation of time are discrete and right in the convection equation to each error velocity component The several progress of local derviation of each coordinate direction of grid node are discrete, obtain the discrete shape of the convection equation of each error velocity component In formula, the processor is used to execute the revision program of the error speed stored in memory, to realize following steps：According to pre- If time point, it is several to the local derviation of time in the convection equation of each error velocity component carry out it is discrete, inside grid Node uses fourth central difference scheme, and forward difference format or backward difference lattice are used for the node in net boundary Formula, progress several to the local derviation of each coordinate direction in the convection equation of each error velocity component are discrete, obtain each error speed Spend the discrete form of the convection equation of component.
Further, the continuity equation to the error speed carries out discrete, obtains the continuity of error speed In the discrete form of equation, the processor is used to execute the revision program of the error speed stored in memory, with realize with Lower step：Progress several to the local derviation of each error velocity component in the continuity equation of error speed are discrete, and to each measurement speed The several progress of local derviation for spending component are discrete, obtain the discrete form of the continuity equation of the error speed.
Further, in the continuity equation to error speed the local derviation of each error velocity component it is several carry out from It dissipates, and discrete to the several progress of the local derviation of each measuring speed, obtains in the discrete form of the continuity equation of the error speed, The processor is used to execute the revision program of the error speed stored in memory, to realize following steps：For in grid The node in portion uses fourth central difference scheme, and forward difference format or backward difference are used for the node in net boundary Format, progress several to the local derviation of each error velocity component in the continuity equation of error speed are discrete, and to each measuring speed The several progress of local derviation of component are discrete, obtain the discrete form of the continuity equation of the error speed.
Further, using each error velocity component as independent variable, accordingly, the goal-selling function is：It is each to become certainly The evolution of the quadratic sum of two norms of amount.
Fourth aspect, the embodiment of the present invention provide a kind of computer storage media, are stored in the computer storage media There are computer executable instructions, which is configured to execute the blood provided in said one or multiple embodiments The modification method of flow velocity degree.
Modification method, device, terminal and the computer storage media for a kind of blood flow velocity that the embodiment of the present invention is provided, This method includes：First, the speed field areas of blood flow is obtained, grid node is distributed in speed field areas, determines speed place The convection equation of each error velocity component of blood flow on the grid node in domain, respectively to the convection equation of each error velocity component into Row is discrete, obtains the discrete form of the convection equation of each error velocity component, similarly, determines the grid section of speed field areas The continuity equation of the error speed of blood flow on point, it is discrete to the continuity equation progress of error speed, obtain error speed The discrete form of continuity equation, in this way, with the company of the discrete form and error speed of the convection equation of each error velocity component The discrete form of continuous property equation is constraint equation, determines the minimum value of goal-selling function, by the corresponding each error speed of minimum value Component is spent, is determined as each error velocity component of blood flow on the grid node of speed field areas, finally, according to each error of blood flow Each velocity component of blood flow on grid node of the velocity component to correct the speed field areas measured；That is, in this hair In bright embodiment, by respectively to the continuity equation of the convection equation of each error velocity component and error speed determined into Row is discrete, obtains the discrete shape of the discrete form of the convection equation of each error velocity component and the continuity equation of error speed Formula, in order to determine each error velocity component of blood flow, with the discrete form and error of the convection equation of each error velocity component The discrete form of the continuity equation of speed is constraint equation, first determines the minimum value of goal-selling function, then by minimum value pair Each error velocity component answered is determined as each error velocity component of blood flow, to be repaiied with each error velocity component of blood flow Each velocity component of positive blood flow improves in this way, being corrected for the error velocity interpolation of time correlation in speed field areas The accuracy of blood flow velocity.
Description of the drawings
Fig. 1 is the flow diagram of the modification method of the blood flow velocity in the embodiment of the present invention；
Fig. 2 is the distribution schematic diagram of blood flow in the raw velocity field in the embodiment of the present invention；
Fig. 3 is the distribution schematic diagram of blood flow in the revised velocity field in the embodiment of the present invention；
Fig. 4 is the structural schematic diagram of the correcting device of the blood flow velocity in the embodiment of the present invention；
Fig. 5 is the structural schematic diagram of the terminal in the embodiment of the present invention；
Fig. 6 is the structural schematic diagram of the computer storage media in the embodiment of the present invention.
Specific implementation mode
Following will be combined with the drawings in the embodiments of the present invention, and technical solution in the embodiment of the present invention carries out clear, complete Site preparation describes.
The embodiment of the present invention provides a kind of modification method of blood flow velocity, and this method can be applied in terminal, and Fig. 1 is this The flow diagram of the modification method of blood flow velocity in inventive embodiments, as shown in Figure 1, the modification method of the blood flow velocity, May include：
S101：Obtain the speed field areas of blood flow；
Wherein, grid node is distributed in speed field areas；
In above-mentioned S101, the speed field areas of blood flow can be obtained from Medical Devices, since in MRI, electronics calculates Machine tomoscan (CT, Computed Tomography) etc. and PIV measurement results are distributed in two dimensional surface or three dimensions On interior structured grid node, so, grid node is distributed in speed field areas, wherein the grid node can be according to two Dimensional plane is distributed, can also be according to three-dimensional spatial distribution, and here, the embodiment of the present invention is not especially limited.
For by taking two dimensional surface is distributed as an example, Fig. 2 is the distribution signal of blood flow in original velocity field in the embodiment of the present invention Figure is illustrated in figure 2 the distribution schematic diagram of the blood flow before not being modified in speed field areas；Speed field areas in image Include abscissa x and ordinate y, velocity field area distribution has grid node, and the unit of the spacing between grid is millimeter, by The region that legend on the right side of in Fig. 2 can be seen that different gray scales indicates different speed, wherein the unit of speed is that millimeter is every Second；For the PIV of two-dimensional incompressible contracting flowing, the speed field areas of blood flow is got, is distributed in the velocity field of blood flow The grid node of mx × my structures, wherein mx, my are respectively number of nodes of the grid node on x, y both directions, Δ x, Δ y It is x, the grid spacing of y both directions, for the speed field data of time resolution, time interval is Δ t.Wherein, uo, voPoint The raw velocity field measured Biao Shi not tested along x, the measuring speed component of y both directions, ue, veIndicate the error of both direction Velocity component, ur, vrIndicate the true velocity component of both direction.
S102：Determine the convection equation of each error velocity component of blood flow on the grid node of speed field areas；
After the grid node for obtaining speed field areas, according to fluid dynamic theory, the grid section of speed field areas Each error component of blood flow is satisfied by convection equation on point, is based on this, by taking the velocity component in the directions x as an example, error velocity component ueConvection equation it is as follows：
Wherein, t is the time,For differential sign, according to uo=ur+ue, vo=vr+ve, above-mentioned convection equation can convert For：
This method is related to error velocity component u it can be seen from above-mentioned formula (2)eTo the time and along two coordinate directions Partial derivative.
So far, error velocity component u has just been obtainedeConvection equation, can also similarly obtain error velocity component vePair Flow equation.
S103：It is discrete to the progress of the convection equation of each error velocity component respectively, obtain the convection current of each error velocity component The discrete form of equation；
For two-dimensional space, it is above-mentioned the convection equation of each error velocity component is carried out respectively discrete may include： To error velocity component ueConvection equation carry out it is discrete, to error velocity component veConvection equation carry out it is discrete, to obtain Error velocity component ueConvection equation discrete form, error velocity component veConvection equation discrete form.
In specific implementation process, S103 may include：
Progress several to the local derviation of time in the convection equation of each error velocity component are discrete, and to each coordinate of grid node The several progress of local derviation in direction are discrete, obtain the discrete form of the convection equation of each error velocity component.
Wherein, in the above-mentioned convection equation to each error velocity component the local derviation of time it is several carry out discrete use from Format is dissipated, it can be poor for order central with the several discrete discrete schemes used of progress of local derviation to each coordinate direction of grid node Cellular and upstreame scheme etc., here, the embodiment of the present invention is not especially limited.
In a kind of optional embodiment, in the above-mentioned convection equation to each error velocity component the local derviation of time it is several into Row is discrete, and discrete to the several progress of the local derviation of each coordinate direction of grid node, obtains the convection equation of each error velocity component Discrete form, may include：
According to preset time point, progress several to the local derviation of time in the convection equation of each error velocity component are discrete；
Fourth central difference scheme is used for the node inside grid, for the node in net boundary using preceding to difference Cellular or backward difference format, it is several to the local derviation of each coordinate direction in the convection equation of each error velocity component carry out from It dissipates；
Obtain the discrete form of the convection equation of each error velocity component.
Wherein, for error velocity component ueConvection equation, in practical applications, it is assumed that t-2, t-1, t, t+1, t+2 Moment, x, the error velocity component u on the directions ye, veValue on grid node (i, j) is respectively ue ij,t-2,ue ij,t-1, ue ij,t,ue ij,t+1,ue ij,t+2, ve ij,t-2,ve ij,t-1,ve ij,t,ve ij,t+1,ve ij,t+2, then x, the error speed in the directions y in t-2, T-1, t, t+1, the t+2 moment can be expressed as the form of column vector：
[ue ij,t-2；ue ij,t-1；ue ij,t；ue ij,t+1；ue ij,t+2；ve ij,t-2；ve ij,t-1；ve ij,t；ve ij,t+1；ve ij,t+2]10M×1 (3)
Wherein, M=mx × my indicates x, and grid number of nodes summation in y both directions, node first increases along the directions y (j), Secondly increase along the directions x (i), traverse the node of entire computational domain.
To error velocity component ueConvection equation carry out it is discrete in, it is several for the local derviation to the time Phase The value answered isFourth central difference scheme, corresponding partial derivative is used to calculate public the time in all grid nodes Formula is：
The coefficient of dispersion matrix Α several to the local derviation of time can be calculated by above-mentioned formula (4)uetIt can be expressed as：
Therefore, error velocity component ueThe discrete form several to the local derviation of time t：
It is several to the local derviation in the directions grid node xCorrespondingly value isGrid internal node is adopted With fourth central difference scheme, forward difference format or backward difference format are used to the node in net boundary, correspondingly, Can be to the several calculation formula of the local derviation in the directions x：
The coefficient of dispersion matrix several to the local derviation in the directions x can be calculated by above-mentioned formula (7)It can be expressed as：
Wherein, Ο is null matrix.Matrix C1,C2,C3,C4,C5It is expressed as：
Then error velocity component ueThe several discrete form of local derviation in the directions x is：
Similarly, the discrete form several to the local derviation in the directions y is：
Wherein, by calculating,It can be expressed as：
Wherein, the C in above-mentioned formula (16)6It is expressed as：
By above-mentioned formula (2) it is found that error velocity component ueConvection equation in there is also nonlinear terms Value is accordinglyFourth central difference scheme is used to grid internal node, the node in net boundary is used Forward difference format or backward difference format, accordinglyCan be：
It can be obtained by above-mentioned formula (15)Coefficient of dispersion matrixIt can be expressed as：
Wherein, Matrix C7,C8For：
ThenDiscrete form be：
Wherein, " " is interior product code.
Similarly, for nonlinear termsDiscrete form be：
Wherein,It can be expressed as
Wherein, Matrix C9It is expressed as：
Wherein, Matrix C91,C92It is expressed as：
Then, velocity error component ueThe discrete form of convection equation is：
Wherein, uo ij,t,vo ij,tRespectively represent uo, voValue on grid node (i, j).
Similarly, according to error velocity component ueThe discrete form of convection equation can obtain error velocity component veConvection current The discrete form of equation is：
Wherein,With the local derviation of time field and velocity field Number calculating is related, referred to herein as derivation operator matrix.Using different discrete schemes, the concrete form of derivation operator matrix is also each It is different, cell centered scheme or upstreame scheme for example, by using higher order precision etc..
So far, the discrete form of the convection equation of each error velocity component has just been obtained.
S104：Determine the continuity equation of the error speed of blood flow on the grid node of speed field areas；
According to fluid dynamic theory, ur, vrMeet following continuity equation：
Since experimental data has residual epsilon on each grid node1, then uo, voMeet following continuity equation：
Then ue, veMeet following continuity equation：
The continuity equation that error speed can be obtained according to formula (31) and formula (32) is as follows：
In this way, just defining the continuity equation of error speed.
S105：It is discrete to the continuity equation progress of error speed, obtain the discrete shape of the continuity equation of error speed Formula；
For formula (33), in discrete to the continuity equation progress of error speed, including：Error velocity component ueTo x The local derviation in direction is several, error velocity component ve, measuring speed component u several to the local derviation in the directions yoIt is several to the local derviation in the directions x With measuring speed component voIt is several to the local derviation in the directions y.
In a kind of optional embodiment, S105 may include：
Progress several to the local derviation of each error velocity component in the continuity equation of error speed are discrete, and to each measurement speed The several progress of local derviation for spending component are discrete, obtain the discrete form of the continuity equation of error speed.
Here, it should be noted that the partial derivative of each error velocity component in the above-mentioned continuity equation to error speed Item progress is discrete, similar with the method that the several progress of the local derviation to each measuring speed component are discrete.
In a kind of optional embodiment, the local derviation of each error velocity component in the above-mentioned continuity equation to error speed Several progress are discrete, and discrete to the several progress of the local derviation of each measuring speed component, obtain the continuity equation of error speed Discrete form may include：
Fourth central difference scheme is used for the node inside grid, for the node in net boundary using preceding to difference Cellular or backward difference format, it is several to the local derviation of each error velocity component in the continuity equation of error speed carry out from It dissipates, and discrete to the several progress of the local derviation of each measuring speed component, obtains the discrete form of the continuity equation of error speed.
In specific implementation process, according toDiscrete method similar method may be used obtain uo, voIt is inclined DerivativeCorresponding coefficient matrix is respectivelyThen error speed continuity equation from Scattered form is：
Here, it should be noted that different discrete schemes equally may be used in the continuity equation of above-mentioned error speed, The embodiment of the present invention is not especially limited this.
So far, the discrete form of the convection equation of each error velocity component has just been obtained by S105.
S106：With the discrete of the discrete form of the convection equation of each error velocity component and the continuity equation of error speed Form is constraint equation, determines the minimum value of goal-selling function, by the corresponding each error velocity component of minimum value, is determined as speed Spend each error velocity component of blood flow on the grid node of field areas；
Wherein, using each error velocity component as independent variable, accordingly, above-mentioned goal-selling function is：Each independent variable The evolution of the quadratic sum of two norms.
For by taking two-dimensional space as an example, the expression formula of goal-selling function is：
Wherein, f is goal-selling function, | | ue ij| |, | | ve ij| | indicate independent variable ue ij,ve ijTwo norms.
Determining for constraint equation solves when in order to solve goal-selling function f minimums, and therefore, which, which is converted into solution, has about The linear multivariate function minimum problems of beam, you can the nonlinear multielement minimum value function fmincon provided using matlab It is solved.
For goal-selling function (35), variable u is takene ij,ve ijInitial value be 0, be not present inequality linear restriction, take For the minimum value and maximum value of true velocity component respectively as inequality linear restriction lower bound and the upper bound, S105 obtains error speed Continuity equation discrete form, formula (34) is used as equation linear restriction, and S103 obtains the convection current side of each error velocity component The discrete form of journey, formula (28) and (29) are used as equation nonlinear restriction, inequality nonlinear restriction, last basis are also not present The parameter that options in fmincon functions is specified optimizes solution, when obtaining goal-selling function minimum, constraint equation Determine solution ue ij,ve ij。
S107：According to each error velocity component of blood flow, to blood flow on the grid node of the speed field areas measured Each velocity component is modified.
After each error velocity component for determining blood flow in S106, blood that can be by following formula to measuring Each velocity component of stream is modified：
Wherein, uc ij,vc ijIt is each velocity component of blood flow in revised velocity field.
Fig. 3 is the distribution schematic diagram of blood flow in the revised velocity field in the embodiment of the present invention, and Fig. 3 is compared in Fig. 2 Speed field areas the variation that occurs of blood flow velocity, it is seen then that the influence with the error speed of time correlation to the speed of blood flow It is very big, have great importance it follows that pair being modified with the error speed of time correlation.
The step with S106 calculate as a result, the measuring speed component to raw velocity field is modified, the present invention Each velocity component u of blood flow in the revised velocity field that embodiment obtainsc ij,vc ijIt is inhibited to measurement error, it carries A velocity field that can more reflect true flow field structure is supplied, to provide better material to the research of flow mechanism below.
It is the example for carrying out mesh generation to velocity field for two-dimensional space above, has three below with velocity field area distribution The modification method of above-mentioned blood flow velocity is illustrated for the grid node of dimension space.
The PIV measurement result velocity field modification methods of three-dimensional potential flows are described, specific implementation step is such as Under：
First, the discrete form of the convection equation of blood flow velocity is determined；
In practical applications, the PIV measurement results of three-dimensional potential flows are distributed across nx × ny × nz regular textures The instantaneous velocity field of three velocity components in three dimensions on grid node, wherein nx, ny, nz are respectively grid node in x, y, Number of nodes on tri- directions z, Δ x, Δ y, Δ z are x, the grid spacing in tri- directions y, z.For the velocity field of time resolution Data, time interval are Δ t.If uo, vo, woIt is survey of the raw velocity field of experiment measurement on x, tri- directions y, z respectively Measure velocity component, ue, ve, weIt is the error velocity component on three directions, ur, vr, weIt is the true velocity point on three directions Amount.By taking the velocity component in the directions x as an example, error velocity component ueConvection equation form be：
Wherein, t is the time,For differential sign.According to uo=ur+ue, vo=vr+ve, wo=wr+we, by above-mentioned convection current side Journey form can be converted to：
The equation is related to error speed u it can be seen from above-mentioned formula (38)eTo the time and to the inclined of three coordinate directions Derivative term.
Assuming that t-1, t, the t+1 moment, x, y, the error velocity component u on tri- directions ze, ve, weGrid node (i, j, K) value on is respectively ue ijk,t-1,ue ijk,t,ue ijk,t+1, ve ijk,t-1,ve ijk,t,ve ijk,t+1, we ijk,t-1,we ijk,t,we ijk ,t+1.The error velocity component in tri- directions then x, y, z can be expressed as the form of column vector at t-1, t, t+1 moment：
[ue ijk,t-1；ue ijk,t；ue ijk,t+1；ve ijk,t-1；ve ijk,t；ve ijk,t+1；we ijk,t-1；we ijk,t；we ijk,t+1]9N×1 (39)
Wherein, N=nx × ny × nz indicates that x, tri- direction number of nodes summations of y, z, node first increase along the directions z (k), It is secondary to traverse the node of entire computational domain finally along the directions x (i) along y (j).
Subsequent operation for convenience, provides error velocity component u firsteConvection equation discrete form,Accordingly Value beUse Second-Order Central Difference format, the corresponding several calculating of local derviation public the time in all grid nodes Formula is as follows：
Coefficient of dispersion matrixIt can be expressed as：
Therefore, error velocity component ueThe discrete form several to the local derviation of time t：
Value is accordinglySecond-Order Central Difference format is used in grid internal node, is used on boundary Forward difference format or backward difference format, the corresponding several calculation formula of local derviation are：
Coefficient of dispersion matrixIt can be expressed as：
Wherein, Matrix C10,C11,C12It is expressed as：
Then error velocity component ueThe several discrete form of local derviation in the directions x is：
Other twoLocal derviation is several similarly to be write out：
It indicates respectively as follows：
Wherein, Matrix C13,C14,C15,C16It is expressed as：
The discrete form of the convection equation of error velocity component there are nonlinear terms, forValue is accordinglyNode inside grid equally uses Second-Order Central Difference format, and forward difference is used in borderline node Format or backward difference format, corresponding calculation formula obtain：
Coefficient of dispersion matrixIt can be expressed as：
Wherein, Matrix C17,C18It is expressed as：
ThenDiscrete form be：
Equally, for nonlinear termsMatrix form after discrete is：
Coefficient of dispersion matrixRespectively：
Wherein, Matrix C19,C20It is expressed as：
Matrix C191,C192It is expressed as：
Then, error velocity component ueThe discrete form of convection equation is：
Wherein, uo ijk,t,vo ijk,t,wo ijk,tRespectively represent uo, vo, woValue on grid node (i, j, k).
Similarly, according to error velocity component ueConvection equation discrete form, error velocity component v can be obtainede, weThe discrete form of convection equation be respectively：
It is related with time field and the partial derivative calculating of velocity field, referred to herein as derivation operator matrix.Using it is different from Format is dissipated, the concrete form of derivation operator matrix is also different, for example, by using the cell centered scheme or upstreame scheme of higher order precision Deng.
Then, it is determined that the discrete form of the continuity equation of error speed；
True flow field meets continuity equation：
Experimental data has residual epsilon on each grid node1, i.e. experimental data meets equation：
Then the residual error formula of error speed is：
According to formula (74) and formula (75), the continuity equation that can obtain error speed is：
It is provided using above-mentionedIt similar can obtain the raw velocity u that experiment measureso, vo, woLocal derviation NumberCorresponding coefficient matrix is respectivelyThe then discrete shape of the continuity equation of error speed Formula is：
Here, it should be noted that different discrete schemes may be used in the continuity equation of error speed, and the present invention is real Example is applied to be not especially limited.
Finally, it is assumed that goal-selling function is：
Wherein, f is goal-selling function, | | ue ijk| |, | | ve ijk| |, | | we ijk| | indicate variable ue ijk,ve ijk,we ijk's Two norms.
The solution of optimization problem is to solve for determining for constraint equation when goal-selling function minimum and solves.Therefore problem conversion To solve constrained linear multivariate function minimum problems, you can the nonlinear multielement minimum value letter provided using matlab Number fmincon is solved.
For goal-selling function (78), variable u is takene ijk,ve ijk,we ijkInitial value be 0, there is no inequality it is linear Constraint, take true velocity component minimum value and maximum value respectively as inequality linear restriction lower bound and the upper bound, error speed The discrete form of continuity equation, i.e. formula (77) be used as equation linear restriction, the convection equation of each error velocity component from The form of dissipating, i.e. formula (70), (71) and (72) are used as equation nonlinear restriction, inequality nonlinear restriction are also not present, finally Solution is optimized according to the options parameters specified, when obtaining object function minimum, constraint equation determines solution ue ijk, ve ijk,we ijk。
Error correction is carried out to raw velocity field according to the following formula：
Wherein, uc ijk,vc ijk,wc ijkIt is each velocity component of blood flow in revised velocity field.
Each velocity component u of blood flow in the revised velocity field that the embodiment of the present invention obtainsc ijk,vc ijk,wc ijk, eliminate The residual error of raw velocity field divergence, inhibited to measurement error, true flow field structure can more be reflected by providing one Velocity field, to provide better material to the research of flow mechanism below.
This have the advantage that：The error of temporal correlation in measurement error can be effectively eliminated, is corrected Velocity field afterwards more approaches the true velocity field of flowing, is corrected to the minimum of raw velocity field, and erection rate field height is loyal to Raw experimental data solves discrete error speed convection equation by optimization method, can obtain two model of error speed The optimal solution of number, is suitable for the velocity field amendment of two and three dimensions, can be realized by matrix operation, convenient for faster parallel It calculates, can realize the error correction of the velocity field measured experiment, the faithful to raw experimental data in erection rate field, and inhibit The velocity error of time correlation in flow field, so as to preferably react the feature in true flow field.
A kind of modification method for blood flow velocity that the embodiment of the present invention is provided, this method include：First, blood flow is obtained Speed field areas is distributed with grid node in speed field areas, determines each error of blood flow on the grid node of speed field areas The convection equation of velocity component, it is discrete to the progress of the convection equation of each error velocity component respectively, obtain each error velocity component The discrete form of convection equation similarly determine the continuous of the error speed of blood flow on the grid node of speed field areas Property equation, the continuity equation of error speed is carried out discrete, obtains the discrete form of the continuity equation of error speed, this Sample is constraint with the discrete form of the discrete form of the convection equation of each error velocity component and the continuity equation of error speed Equation determines the minimum value of goal-selling function, by the corresponding each error velocity component of minimum value, is determined as speed field areas Each error velocity component of blood flow finally corrects the speed measured according to each error velocity component of blood flow on grid node Spend each velocity component of blood flow on the grid node of field areas；That is, in embodiments of the present invention, by respectively to determination The convection equation of each error velocity component gone out and the continuity equation progress of error speed are discrete, obtain each error velocity component Convection equation discrete form and error speed continuity equation discrete form, in order to determine blood flow each error speed Component is spent, the discrete form with the discrete form of the convection equation of each error velocity component and the continuity equation of error speed is Constraint equation, first determines the minimum value of goal-selling function, then by the corresponding each error velocity component of minimum value, is determined as blood flow Each error velocity component, to correct each velocity component of blood flow with each error velocity component of blood flow, in this way, for The error velocity interpolation of time correlation is corrected, and the accuracy of blood flow velocity in speed field areas is improved.
Based on same inventive concept, the present embodiment provides a kind of correcting device of blood flow velocity, Fig. 4 is the embodiment of the present invention In blood flow velocity correcting device structural schematic diagram, as shown in figure 4, the device includes：Acquisition module 41, first determines mould Block 42, the first discrete block 43, the second determining module 44, the second discrete block 45, third determining module 46 and correcting module 47；
Wherein, grid node is distributed in speed field areas in acquisition module 41, the speed field areas for obtaining blood flow； First determining module 42, the convection equation of each error velocity component of blood flow on the grid node for determining speed field areas； First discrete block 43, for respectively to the convection equation of each error velocity component carry out it is discrete, obtain each error velocity component Convection equation discrete form；Second determining module 44, the error of blood flow on the grid node for determining speed field areas The continuity equation of speed；Second discrete block 45, it is discrete for the continuity equation progress to error speed, obtain error speed The discrete form of the continuity equation of degree；Third determining module 46, for the discrete of the convection equation of each error velocity component The discrete form of the continuity equation of form and error speed is constraint equation, determines the minimum value of goal-selling function, will most It is small to be worth corresponding each error velocity component, it is determined as each error velocity component of blood flow on the grid node of speed field areas；It repaiies Positive module 47, for each error velocity component according to blood flow, to blood flow on the grid node of the speed field areas measured Each velocity component is modified.
In a kind of optional embodiment, the first discrete block 43 is specifically used for：To the convection current side of each error velocity component The several progress of the local derviation of time are discrete in journey, and discrete to the several progress of the local derviation of each coordinate direction of grid node, obtain each mistake The discrete form of the convection equation of poor velocity component.
In a kind of optional embodiment, the first discrete block 43 is to the time in the convection equation of each error velocity component The several progress of local derviation are discrete, and discrete to the several progress of the local derviation of each coordinate direction of grid node, obtain each error velocity component Convection equation discrete form, may include：According to preset time point, when in the convection equation of each error velocity component Between local derviation it is several carry out it is discrete；Fourth central difference scheme is used for the node inside grid, in net boundary Node uses forward difference format or backward difference format, to each coordinate direction in the convection equation of each error velocity component The several progress of local derviation are discrete；Obtain the discrete form of the convection equation of each error velocity component.
In a kind of optional embodiment, the second discrete block 45 is specifically used for：To in the continuity equation of error speed The several progress of local derviation of each error velocity component are discrete, and discrete to the several progress of the local derviation of each measuring speed component, missed The discrete form of the continuity equation of differential degree.
In a kind of optional embodiment, the second discrete block 45 is to each error speed in the continuity equation of error speed The several progress of local derviation of component are discrete, and discrete to the several progress of the local derviation of each measuring speed component, obtain the company of error speed The discrete form of continuous property equation may include：Fourth central difference scheme is used for the node inside grid, for Grid Edge Node in boundary uses forward difference format or backward difference format, to each error speed in the continuity equation of error speed The several progress of local derviation of component are discrete, and discrete to the several progress of the local derviation of each measuring speed component, obtain the company of error speed The discrete form of continuous property equation.
Wherein, using each error velocity component as independent variable, accordingly, above-mentioned goal-selling function may include：Each certainly The evolution of the quadratic sum of two norms of variable.
In practical applications, above-mentioned acquisition module 41, the first determining module 42, the first discrete block 43, second determine mould Block 44, the second discrete block 45, third determining module 46 and correcting module 47 can realize by the processor on device, specifically For central processing unit (CPU, Central Processing Unit), microprocessor (MPU, Microprocessor Unit), Digital signal processor (DSP, Digital Signal Processing) or field programmable gate array (FPGA, Field Programmable Gate Array) etc. realizations.
Based on same inventive concept, the present embodiment provides a kind of terminal, Fig. 5 is the structure of the terminal in the embodiment of the present invention Schematic diagram, as shown in figure 5, the terminal includes：Processor 51, memory 52 and communication bus 53；Wherein, communication bus 53 is used for Realize the connection communication between processor 51 and memory 52；The processor 51 is for executing the blood stored in memory The revision program of flow velocity degree, to realize following steps：
The speed field areas of blood flow is obtained, grid node is distributed in speed field areas；Determine the grid of speed field areas The convection equation of each error velocity component of blood flow on node；It is discrete to the progress of the convection equation of each error velocity component respectively, Obtain the discrete form of the convection equation of each error velocity component；Determine the error speed of blood flow on the grid node of speed field areas The continuity equation of degree；It is discrete to the continuity equation progress of error speed, obtain the discrete of the continuity equation of error speed Form；Discrete form with the discrete form of the convection equation of each error velocity component and the continuity equation of error speed is about Shu Fangcheng determines the minimum value of goal-selling function, by the corresponding each error velocity component of minimum value, is determined as speed field areas Grid node on blood flow each error velocity component；According to each error velocity component of blood flow, to the speed place measured Each velocity component of blood flow is modified on the grid node in domain.
It is described discrete to the progress of the convection equation of each error velocity component respectively in a kind of optional embodiment, it obtains In the discrete form of the convection equation of each error velocity component, the processor 51 is additionally operable to execute the amendment of the blood flow velocity Program, to realize following steps：
Progress several to the local derviation of time in the convection equation of each error velocity component are discrete, and to each coordinate of grid node The several progress of local derviation in direction are discrete, obtain the discrete form of the convection equation of each error velocity component.
In a kind of optional embodiment, in the convection equation to each error velocity component the local derviation of time it is several into Row is discrete, and discrete to the several progress of the local derviation of each coordinate direction of grid node, obtains the convection equation of each error velocity component Discrete form in, the processor 51 is additionally operable to execute the revision program of the blood flow velocity, to realize following steps：
According to preset time point, progress several to the local derviation of time in the convection equation of each error velocity component are discrete； Use fourth central difference scheme for the node inside grid, for the node in net boundary using forward difference format or Person's backward difference format, progress several to the local derviation of each coordinate direction in the convection equation of each error velocity component are discrete；It obtains The discrete form of the convection equation of each error velocity component.
In a kind of optional embodiment, the continuity equation progress to error speed is discrete, obtains error speed Continuity equation discrete form in, the processor 51 is additionally operable to execute the revision program of the blood flow velocity, with realize Following steps：
Progress several to the local derviation of each error velocity component in the continuity equation of error speed are discrete, and to each measurement speed The several progress of local derviation for spending component are discrete, obtain the discrete form of the continuity equation of error speed.
In a kind of optional embodiment, the local derviation of each error velocity component in the continuity equation to error speed Several progress are discrete, and discrete to the several progress of the local derviation of each measuring speed, obtain the discrete of the continuity equation of error speed In form, the processor 41 is additionally operable to execute the revision program of the blood flow velocity, to realize following steps：
Fourth central difference scheme is used for the node inside grid, for the node in net boundary using preceding to difference Cellular or backward difference format, it is several to the local derviation of each error velocity component in the continuity equation of error speed carry out from It dissipates, and discrete to the several progress of the local derviation of each measuring speed component, obtains the discrete form of the continuity equation of error speed.
Wherein, using each error velocity component as independent variable, accordingly, above-mentioned goal-selling function may include：Each certainly The evolution of the quadratic sum of two norms of variable.
The embodiment of the present invention provides a kind of computer storage media, and Fig. 6 is that the computer storage in the embodiment of the present invention is situated between The structural schematic diagram of matter, as shown in fig. 6, computer executable instructions are stored in the computer storage media 600, the calculating Machine executable instruction is configured to execute the modification method for the blood flow velocity that other embodiments of the invention provide.
It should be understood by those skilled in the art that, the embodiment of the present invention can be provided as method, system or computer program Product.Therefore, the shape of hardware embodiment, software implementation or embodiment combining software and hardware aspects can be used in the present invention Formula.Moreover, the present invention can be used can use storage in the computer that one or more wherein includes computer usable program code The form for the computer program product implemented on medium (including but not limited to magnetic disk storage and optical memory etc.).
The present invention be with reference to according to the method for the embodiment of the present invention, the flow of equipment (system) and computer program product Figure and/or block diagram describe.It should be understood that can be realized by computer program instructions every first-class in flowchart and/or the block diagram The combination of flow and/or box in journey and/or box and flowchart and/or the block diagram.These computer programs can be provided Instruct the processor of all-purpose computer, special purpose computer, Embedded Processor or other programmable data processing devices to produce A raw machine so that the instruction executed by computer or the processor of other programmable data processing devices is generated for real The device for the function of being specified in present one flow of flow chart or one box of multiple flows and/or block diagram or multiple boxes.
These computer program instructions, which may also be stored in, can guide computer or other programmable data processing devices with spy Determine in the computer-readable memory that mode works so that instruction generation stored in the computer readable memory includes referring to Enable the manufacture of device, the command device realize in one flow of flow chart or multiple flows and/or one box of block diagram or The function of being specified in multiple boxes.
These computer program instructions also can be loaded onto a computer or other programmable data processing device so that count Series of operation steps are executed on calculation machine or other programmable devices to generate computer implemented processing, in computer or The instruction executed on other programmable devices is provided for realizing in one flow of flow chart or multiple flows and/or block diagram one The step of function of being specified in a box or multiple boxes.
The foregoing is only a preferred embodiment of the present invention, is not intended to limit the scope of the present invention.

Claims (12)

1. a kind of modification method of blood flow velocity, which is characterized in that including：
Obtain the speed field areas of blood flow；Wherein, grid node is distributed in the speed field areas；
Determine the convection equation of each error velocity component of blood flow on the grid node of the speed field areas；
It is discrete to the progress of the convection equation of each error velocity component respectively, obtain the discrete of the convection equation of each error velocity component Form；
Determine the continuity equation of the error speed of blood flow on the grid node of the speed field areas；
It is discrete to the continuity equation progress of the error speed, obtain the discrete form of the continuity equation of error speed；
With the discrete of the discrete form of the convection equation of each error velocity component and the continuity equation of the error speed Form is constraint equation, determines the minimum value of goal-selling function, and the corresponding each error velocity component of the minimum value determines For each error velocity component of blood flow on the grid node of the speed field areas；
According to each error velocity component of the blood flow, on the grid node of the speed field areas measured blood flow it is each Velocity component is modified.
2. according to the method described in claim 1, it is characterized in that, it is described respectively to the convection equation of each error velocity component into Row is discrete, obtains the discrete form of the convection equation of each error velocity component, including：
Progress several to the local derviation of time in the convection equation of each error velocity component are discrete, and to each coordinate direction of grid node Local derviation it is several carry out it is discrete, obtain the discrete form of the convection equation of each error velocity component.
3. according to the method described in claim 2, it is characterized in that, time in the convection equation to each error velocity component Local derviation it is several carry out discrete, and carry out discrete to the local derviation of each coordinate direction of grid node is several, obtain each error speed point The discrete form of the convection equation of amount, including：
According to preset time point, progress several to the local derviation of time in the convection equation of each error velocity component are discrete；
Fourth central difference scheme is used for the node inside grid, forward difference lattice are used for the node in net boundary Formula or backward difference format, progress several to the local derviation of each coordinate direction in the convection equation of each error velocity component are discrete；
Obtain the discrete form of the convection equation of each error velocity component.
4. according to the method described in claim 1, it is characterized in that, the continuity equation to the error speed carry out from It dissipates, obtains the discrete form of the continuity equation of error speed, including：
Progress several to the local derviation of each error velocity component in the continuity equation of error speed are discrete, and to each measuring speed point The several progress of local derviation of amount are discrete, obtain the discrete form of the continuity equation of the error speed.
5. according to the method described in claim 4, it is characterized in that, each error is fast in the continuity equation to error speed The several progress of local derviation for spending component are discrete, and discrete to the several progress of the local derviation of each measuring speed, obtain the error speed The discrete form of continuity equation, including：
Fourth central difference scheme is used for the node inside grid, forward difference lattice are used for the node in net boundary Formula or backward difference format, progress several to the local derviation of each error velocity component in the continuity equation of error speed are discrete, And it is discrete to the several progress of the local derviation of each measuring speed component, obtain the discrete form of the continuity equation of the error speed.
6. accordingly, described according to the method described in claim 1, it is characterized in that, using each error velocity component as independent variable Goal-selling function is：
The evolution of the quadratic sum of two norms of each independent variable.
7. a kind of correcting device of blood flow velocity, which is characterized in that including：
Acquisition module, the speed field areas for obtaining blood flow；Wherein, grid node is distributed in the speed field areas；
First determining module, the convection current of each error velocity component of blood flow on the grid node for determining the speed field areas Equation；
First discrete block, carries out the convection equation of each error velocity component for respectively discrete, obtains each error speed point The discrete form of the convection equation of amount；
Second determining module, the continuity side of the error speed of blood flow on the grid node for determining the speed field areas Journey；
Second discrete block, it is discrete for the continuity equation progress to the error speed, obtain the continuity of error speed The discrete form of equation；
Third determining module, for the discrete form of the convection equation of each error velocity component and the error speed The discrete form of continuity equation is constraint equation, determines the minimum value of goal-selling function, and the minimum value is corresponding each Error velocity component is determined as each error velocity component of blood flow on the grid node of the speed field areas；
Correcting module, for each error velocity component according to the blood flow, to the grid of the speed field areas measured Each velocity component of blood flow is modified on node.
8. device according to claim 7, which is characterized in that first discrete block is specifically used for：
Progress several to the local derviation of time in the convection equation of each error velocity component are discrete, and to each coordinate direction of grid node Local derviation it is several carry out it is discrete, obtain the discrete form of the convection equation of each error velocity component.
9. device according to claim 7, which is characterized in that second discrete block is specifically used for：
Progress several to the local derviation of each error velocity component in the continuity equation of error speed are discrete, and to each measuring speed point The several progress of local derviation of amount are discrete, obtain the discrete form of the continuity equation of the error speed.
10. device according to claim 7, which is characterized in that using each error velocity component as independent variable, accordingly, institute Stating goal-selling function is：
The evolution of the quadratic sum of two norms of each independent variable.
11. a kind of terminal, which is characterized in that the terminal includes processor, memory and communication bus；
The communication bus is for realizing the connection communication between processor and memory；
The processor is used to execute the revision program that blood flow velocity is stored in memory, to realize following steps：
Obtain the speed field areas of blood flow；Wherein, grid node is distributed in the speed field areas；
Determine the convection equation of each error velocity component of blood flow on the grid node of the speed field areas；
It is discrete to the progress of the convection equation of each error velocity component respectively, obtain the discrete of the convection equation of each error velocity component Form；
Determine the continuity equation of the error speed of blood flow on the grid node of the speed field areas；
It is discrete to the continuity equation progress of the error speed, obtain the discrete form of the continuity equation of error speed；
With the discrete of the discrete form of the convection equation of each error velocity component and the continuity equation of the error speed Form is constraint equation, determines the minimum value of goal-selling function, and the corresponding each error velocity component of the minimum value determines For each error velocity component of blood flow on the grid node of the speed field areas；
According to each error velocity component of the blood flow, on the grid node of the speed field areas measured blood flow it is each Velocity component is modified.
12. a kind of computer storage media, which is characterized in that be stored with the executable finger of computer in the computer storage media It enables, which is configured to execute the amendment side for the blood flow velocity that any one of the claims 1 to 6 provide Method.
CN201810517956.7A 2018-05-25 2018-05-25 Blood flow velocity correction method, device, terminal and computer storage medium Active CN108805949B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810517956.7A CN108805949B (en) 2018-05-25 2018-05-25 Blood flow velocity correction method, device, terminal and computer storage medium

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810517956.7A CN108805949B (en) 2018-05-25 2018-05-25 Blood flow velocity correction method, device, terminal and computer storage medium

Publications (2)

Publication Number Publication Date
CN108805949A true CN108805949A (en) 2018-11-13
CN108805949B CN108805949B (en) 2022-04-12

Family

ID=64089167

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810517956.7A Active CN108805949B (en) 2018-05-25 2018-05-25 Blood flow velocity correction method, device, terminal and computer storage medium

Country Status (1)

Country Link
CN (1) CN108805949B (en)

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20150066447A1 (en) * 2013-08-27 2015-03-05 Halliburton Energy Services, Inc. Generating a flow model grid based on truncation error threshold values
CN104636566A (en) * 2015-03-13 2015-05-20 中国矿业大学 Dynamic mesh numerical solution method based on modified fluid motion equation
CN104769593A (en) * 2012-10-31 2015-07-08 旭硝子株式会社 Simulation device, simulation method, and program
CN105095555A (en) * 2014-07-15 2015-11-25 北京航空航天大学 Non-divergence smoothing processing method and apparatus for velocity field
CN105929193A (en) * 2016-04-15 2016-09-07 北京航空航天大学 Velocity field quick correction method and device based on hydrodynamics continuity equation

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104769593A (en) * 2012-10-31 2015-07-08 旭硝子株式会社 Simulation device, simulation method, and program
US20150066447A1 (en) * 2013-08-27 2015-03-05 Halliburton Energy Services, Inc. Generating a flow model grid based on truncation error threshold values
CN105095555A (en) * 2014-07-15 2015-11-25 北京航空航天大学 Non-divergence smoothing processing method and apparatus for velocity field
CN104636566A (en) * 2015-03-13 2015-05-20 中国矿业大学 Dynamic mesh numerical solution method based on modified fluid motion equation
CN105929193A (en) * 2016-04-15 2016-09-07 北京航空航天大学 Velocity field quick correction method and device based on hydrodynamics continuity equation

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
NIU LILI ET AL.: "Algorithms for Correcting Velocity Vectors in Ultrasonic Particle Image Velocimetry", 《2009 3RD INTERNATIONAL CONFERENCE ON BIOINFORMATICS AND BIOMEDICAL ENGINEERING》 *
WANG HONGPING ET AL.: "Spatial pyramidal cross correlation for particle image velocimetry", 《SCIENCE CHINA(TECHNOLOGICAL SCIENCES)》 *
高琪等: "基于连续性条件的体PIV后处理技术", 《北京航空航天大学学报》 *

Also Published As

Publication number Publication date
CN108805949B (en) 2022-04-12

Similar Documents

Publication Publication Date Title
Guo A simple algorithm for fitting a Gaussian function [DSP tips and tricks]
Gates et al. Towards high performance digital volume correlation
US20110268334A1 (en) Apparatus for Improving Image Resolution and Apparatus for Super-Resolution Photography Using Wobble Motion and Point Spread Function (PSF), in Positron Emission Tomography
CN106659449B (en) Quantitative dark field imaging in tomography
US9263164B2 (en) Method and device for correcting artefacts during X-ray imagery, especially computer tomography, with a moving modulator field
Huhn et al. FFT integration of instantaneous 3D pressure gradient fields measured by Lagrangian particle tracking in turbulent flows
Callaghan et al. Spatial resolution and velocity field improvement of 4D‐flow MRI
Sitek Representation of photon limited data in emission tomography using origin ensembles
Schultze et al. Performance of hull-detection algorithms for proton computed tomography reconstruction
Faragó et al. Syris: a flexible and efficient framework for X-ray imaging experiments simulation
Agarwal et al. Reconstructing velocity and pressure from noisy sparse particle tracks using constrained cost minimization
US9196063B2 (en) Super-resolution apparatus and method
US9501819B2 (en) Super-resolution apparatus and method
CN101980304A (en) Three-dimensional digital volume image distortion measuring method
CN104851080A (en) TV-based 3D positron emission tomography (PET) image reconstruction method
KR20130039566A (en) Cpmpton camera and method for reconstructing image for resolution recovery thereof
Lesaint et al. GCC and FBCC for linear tomosynthesis
Wen et al. Flexible retrospective phase stepping in x-ray scatter correction and phase contrast imaging using structured illumination
Zacharopoulos et al. 3D shape reconstruction in optical tomography using spherical harmonics and BEM
Li Consistency equations in native detector coordinates and timing calibration for time-of-flight PET
CN108805949A (en) A kind of modification method of blood flow velocity, device, terminal and computer storage media
US9245359B2 (en) Apparatus and method for generating medical image using linear gamma ray source
Karimi et al. A hybrid stochastic-deterministic gradient descent algorithm for image reconstruction in cone-beam computed tomography
Fournier et al. Pinsker estimators for local helioseismology: inversion of travel times for mass-conserving flows
Kupinski et al. Scanning linear estimation: improvements over region of interest (ROI) methods

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
CB02 Change of applicant information

Address after: Room 707, building 2, Xizi international jinzuo, Nanyuan street, Yuhang District, Hangzhou City, Zhejiang Province

Applicant after: Hangzhou Sheng Shi Technology Co., Ltd.

Address before: 311100 Hangzhou, Yuhang, Zhejiang Linping new town Nanyuan street, new far CBC1 block 1502 room.

Applicant before: Hangzhou Sheng Shi Technology Co., Ltd.

CB02 Change of applicant information
GR01 Patent grant
GR01 Patent grant