CN114271924B - Catheter calibration method and device based on grid partition - Google Patents
Catheter calibration method and device based on grid partition Download PDFInfo
- Publication number
- CN114271924B CN114271924B CN202011035782.4A CN202011035782A CN114271924B CN 114271924 B CN114271924 B CN 114271924B CN 202011035782 A CN202011035782 A CN 202011035782A CN 114271924 B CN114271924 B CN 114271924B
- Authority
- CN
- China
- Prior art keywords
- catheter
- pressure
- matrix
- parameter
- basic
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 44
- 238000005192 partition Methods 0.000 title claims abstract description 19
- 239000011159 matrix material Substances 0.000 claims abstract description 102
- 239000013598 vector Substances 0.000 claims abstract description 44
- 238000003491 array Methods 0.000 claims description 11
- 238000004364 calculation method Methods 0.000 claims description 8
- 239000002131 composite material Substances 0.000 claims description 3
- 239000000203 mixture Substances 0.000 claims description 3
- 238000003825 pressing Methods 0.000 claims description 2
- 241000208340 Araliaceae Species 0.000 claims 2
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 claims 2
- 235000003140 Panax quinquefolius Nutrition 0.000 claims 2
- 235000008434 ginseng Nutrition 0.000 claims 2
- 238000002679 ablation Methods 0.000 abstract description 4
- 238000013507 mapping Methods 0.000 abstract description 2
- 238000010586 diagram Methods 0.000 description 6
- 230000000694 effects Effects 0.000 description 3
- 229920001971 elastomer Polymers 0.000 description 3
- 239000000806 elastomer Substances 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 238000004519 manufacturing process Methods 0.000 description 3
- 210000001519 tissue Anatomy 0.000 description 3
- 230000006793 arrhythmia Effects 0.000 description 2
- 206010003119 arrhythmia Diseases 0.000 description 2
- 201000010099 disease Diseases 0.000 description 2
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 230000002107 myocardial effect Effects 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000013500 data storage Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 210000005003 heart tissue Anatomy 0.000 description 1
- 230000006698 induction Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 210000000056 organ Anatomy 0.000 description 1
- 238000007674 radiofrequency ablation Methods 0.000 description 1
- 238000003860 storage Methods 0.000 description 1
- 230000003245 working effect Effects 0.000 description 1
Landscapes
- Media Introduction/Drainage Providing Device (AREA)
- Measuring And Recording Apparatus For Diagnosis (AREA)
Abstract
The invention relates to the field of electrophysiological ablation and mapping, in particular to a catheter calibration method and device based on grid partition. The method comprises the following steps: s1, establishing a parameter matrix at a catheter stress end according to a preset azimuth angle, an elevation angle and an acting force; vectorizing the parameter matrix, and establishing a basic parameter vector matrix; s2, applying acting force to the stress end of the catheter according to the parameter matrix, and simultaneously collecting the pressure value of the pressure sensor to construct a corresponding basic pressure matrix; and S3, establishing a pressure value calibration relation model according to the basic pressure matrix and the basic parameter vector matrix, wherein the pressure value calibration relation model is used for solving the stress magnitude and the stress direction of the catheter according to the pressure sensor acquired in real time. The pressure calibration method of the electrophysiological catheter with pressure sensing reduces the requirement on the linearity of the pressure sensor, and can quickly calibrate and output the pressure of the catheter with high precision.
Description
Technical Field
The invention relates to the field of electrophysiological ablation and mapping, in particular to a catheter calibration method and device based on grid partition.
Background
Arrhythmia is one of the common diseases of arrhythmia in the world, and clinically, the radio frequency ablation technology by using a catheter has been widely applied to treat the diseases. Radio frequency energy is delivered through the catheter to the electrodes and the site where the electrodes contact and surrounding myocardial tissue for ablation. Clinical results prove that the better treatment effect can be achieved only under the condition that the contact pressure of the electrode at the distal end of the catheter and myocardial tissue is proper. In ablation catheter treatment, the catheter is inserted into the heart with the distal end of the catheter in contact with the inner wall of the heart, where it is often important to have the distal end of the catheter in good contact with the inner wall of the heart and to determine the correct direction and location of abutment, otherwise excessive pressure or incorrect location of abutment may cause undesirable damage to the heart tissue, even perforation of the heart wall, while precise positioning of the catheter is also critical.
The prior art adopts electromagnetic or optical technology to measure the abutting pressure of the distal end and the tissue, has higher and complex equipment requirements and relatively higher manufacturing cost. As mentioned in chinese patent CN103908337a, a sensor using magnetic induction is added to the catheter to sense the contact force between the distal end of the catheter and the organ, and such a sensor is easily distorted by the interference of external magnetic field in application, and the technology requires a plurality of magnetic sensors to be installed in a very small space of the distal end of the catheter, which makes the process difficult, thereby increasing the manufacturing cost. In modern scientific development, pressure is measured by simply using a strain gauge, but the pressure is greatly influenced by temperature, so that no pressure conduit adopting a strain gauge mode exists at present.
Currently, if strain gauges are used to detect catheter pressure, a matched calibration device and method are needed, and the calibration of strain gauge sensors is extremely complex and extremely high in equipment, so that an efficient and reliable calibration method and device system are needed to meet the requirement of mass production.
Disclosure of Invention
The invention solves the problems of extremely complex calibration and extremely high equipment of strain gauge sensors in the prior art, and provides a catheter calibration method and device based on grid partition.
In order to achieve the above object, the present invention provides the following technical solutions:
a grid partition-based catheter calibration method for a catheter with a plurality of pressure sensors at a stress end, comprising the following steps:
s1, establishing a parameter matrix at a catheter stress end according to a preset azimuth angle, an elevation angle and an acting force; vectorizing the parameter matrix, and establishing a basic parameter vector matrix;
s2, applying acting force to the stress end of the catheter according to the parameter matrix, collecting the pressure value of the pressure sensor, and constructing a basic pressure matrix, wherein an array in the basic pressure matrix and a vector in the basic parameter vector matrix are in one-to-one correspondence;
and S3, establishing a pressure value calibration relation model according to the basic pressure matrix and the basic parameter vector matrix, wherein the pressure value calibration relation model is used for solving the stress magnitude and the stress direction of the catheter according to the value of the pressure sensor acquired in real time.
As a preferred scheme of the invention, the parameter matrix establishment process of the step S1 is as follows: the center of the cross section of the catheter is taken as the origin of spherical coordinates, the cross section of the catheter is equally divided into N azimuth angles, the angle from the axis of the catheter to the cross section of the catheter is equally divided into K elevation angles, and the acting force in the direction of the axis of the catheter is a scalar distance.
As a preferable scheme of the invention, the value of N is 12, and the value of azimuth angle is: 0 °,30 °,60 °,90 °,120 °,150 °,180 °,210 °,240 °,270 °,300 °,330 °,360 °.
As a preferable scheme of the invention, the value of K is 4, and the value of the elevation angle is: 0 °,30 °,60 °,90 °.
As a preferable mode of the invention, the acting force in the axial direction of the catheter takes the values of 10g,20g,30g,40g,50g,60g and 70g.
In a preferred embodiment of the present invention, in step S1, the vector in the base parameter vector matrix is represented by a first parameter δ i Second parameter epsilon i And a third parameter mu i The composition of the composite material comprises the components,
first parameter delta i The calculation formula of (2) is as follows:
second parameter epsilon i The calculation formula of (2) is as follows:
third parameter mu i The calculation formula of (2) is as follows:
wherein alpha is i Is the elevation angle, beta i The azimuth angle is F, the acting force is 1.ltoreq.i.ltoreq.n, and n is the total number of arrays formed according to the azimuth angle, the elevation angle and the acting force in the parameter matrix.
As a preferred embodiment of the present invention, the pressure value calibration relation model is:
M i ＝pinv(Q i )*R i
wherein M is i Cell relation matrix of pressure value calibration relation model, R i Is a basic parameter vector matrix, Q i The matrix is a basic pressure matrix, an array in the basic pressure matrix is a pressure value array corresponding to a vector in a basic parameter vector matrix, pinv () is an inversion formula, i is more than or equal to 1 and less than or equal to n, and n is the total number of vectors in the basic parameter vector matrix.
Based on the same conception, the invention also provides a method for detecting the stress of the catheter, which is used for the catheter with a plurality of pressure sensors distributed at the stress end, and comprises the following steps:
a1, simultaneously acquiring pressure values of a plurality of pressure sensors to form a pressure acquisition array;
a2, a pressure value calibration relation model obtained by adopting any one of the calibration methods is adopted, and a reference cell relation matrix with the shortest relative distance with the pressure acquisition array is found from the pressure value calibration relation model;
a3, according to the reference cell relation matrix and the pressure acquisition array, the stress magnitude and the stress direction of the catheter are calculated.
As a preferred embodiment of the present invention, the step A2 specifically includes the steps of:
a21, solving the relative distance between the current pressure acquisition array and the array in the basic pressure matrix;
a22, finding the index number of the corresponding array when the relative distance is minimum;
a23, finding out a corresponding reference cell relation matrix in the pressure value calibration relation model according to the index number.
Based on the same conception, the invention also provides a catheter calibration device based on grid partition, which comprises a force application device, a force application size measuring device and an azimuth angle and elevation angle control device for controlling the force application direction, wherein the force application device is used for applying pressure to the stress end of the catheter,
the force application measuring device is used for collecting the pressure value applied to the stress end of the guide tube in real time,
the azimuth angle and elevation angle control device for controlling the force application direction is used for controlling the force application device to apply pressure to the force receiving end of the catheter according to the preset azimuth angle and elevation angle,
the force application device, the force application magnitude measuring device and the azimuth angle and elevation angle control device for controlling the force application direction apply force to the force receiving end of the catheter according to any one of the calibration methods.
Compared with the prior art, the invention has the beneficial effects that:
the invention discloses a pressure calibration method and a pressure calibration device for a pressure sensing electrophysiological catheter.
Description of the drawings:
FIG. 1 is a schematic view of a pressure sensing catheter used in the method of the present invention;
FIG. 2 is a schematic view of the distal electrode layout of the catheter on the pressure sensing catheter used in the method of the present invention;
FIG. 3 is a flow chart of a grid partition catheter calibration method according to embodiment 1 of the present invention;
FIG. 4 is a schematic diagram of the present invention in example 1 for establishing a spherical coordinate system at the distal end of a catheter according to preset azimuth, elevation and force;
FIG. 5 is a schematic view of the azimuth partition in the top view of embodiment 1 of the present invention;
FIG. 6 is a schematic diagram of a side elevation view partition in embodiment 1 of the present invention;
FIG. 7 is a diagram showing a calibrating apparatus according to embodiment 2 of the present invention;
FIG. 8 is a diagram showing a calibrating apparatus according to embodiment 3 of the present invention;
FIG. 9 is a flow chart of a method for detecting stress of a catheter according to embodiment 4 of the present invention;
fig. 10 is a graph showing the calibration effect in embodiment 4 of the present invention.
Detailed Description
The present invention will be described in further detail with reference to test examples and specific embodiments. It should not be construed that the scope of the above subject matter of the present invention is limited to the following embodiments, and all techniques realized based on the present invention are within the scope of the present invention.
The pressure sensing catheter, as shown in fig. 1, can be divided into three parts, distal, proximal, and handle, wherein the distal end contains a pressure sensor. The arrangement of the pressure sensors at the distal end of the catheter has various forms, an example of arranging three and four pressure sensors at the distal end of the pressure sensing catheter is shown in fig. 1, three pressure sensors (or four pressure sensors) are arranged along the circumferential direction of the catheter when the cross section of the catheter is seen, by the calibration method of the invention, real-time calibration output of the pressure value and the pressure direction of the catheter can be realized, and fig. 2 shows a schematic diagram of the arrangement of the electrode at the distal end of the catheter, which is used for energy output, and the stress direction and the size of the distal end of the catheter directly influence the working effect of the electrode.
Specifically, the pressure sensing catheter distal electrode is fixedly connected with an elastomer, the elastomer enables the distal electrode to return to a natural state after pressure is removed, and a plurality of pressure sensors are uniformly distributed in a stress concentration area on the elastomer. Taking 3 pressure sensors as an example, the 3 pressure sensors are symmetrically arranged on the elastic body. The value reflected by each pressure sensor when sensing stretching and compressing can be regarded as a vector with directivity (namely positive and negative), and when the pressure sensing catheter is subjected to pressure with direction, the three pressure sensors can detect the corresponding vector value with uniqueness. And measuring different numerical groups of the pressure sensors (each group is provided with positive and negative numerical values of 3 pressure sensors) when various pressures and various directions are known, and sorting the obtained corresponding relations into a database storage data storage module. The method can correct the stress magnitude and direction of the catheter according to pre-stored data in the actual measurement process, so that the calculated value is more accurate.
Example 1
A flowchart of a grid partition based catheter calibration method, as shown in fig. 3, includes the steps of:
s1, a spherical coordinate system is established at a catheter stress end according to a preset azimuth angle, an elevation angle and acting force along the axial direction of the catheter, and a parameter matrix is established.
S2, vectorizing the parameter matrix, and establishing a basic parameter vector matrix.
And S3, applying acting force to the stress end of the catheter according to the parameter matrix, and simultaneously collecting the pressure value of the pressure sensor to construct a basic pressure matrix, wherein an array in the basic pressure matrix and a vector in the basic parameter vector matrix are in one-to-one correspondence.
And S4, establishing a pressure value calibration relation model according to the basic pressure matrix and the basic parameter vector matrix, wherein the pressure value calibration relation model is used for solving the stress magnitude and the stress direction of the catheter according to the pressure sensor acquired in real time.
In step S1, a spherical coordinate system is established at the distal end of the catheter according to the preset azimuth, elevation and acting force, and a schematic diagram of the establishment of the spherical coordinate system at the distal end of the catheter according to the preset azimuth, elevation and acting force is shown in fig. 4. The direction of the axis of the catheter is taken as a Z axis, the center of the cross section of the catheter is taken as an origin of spherical coordinates, the cross section of the catheter is taken as an XOY plane, a spherical coordinate system is correspondingly established, the azimuth angle is beta, the elevation angle is alpha, the cross section of the catheter is equally divided into N azimuth angles, and the angle from the axis of the catheter to the cross section of the catheter is equally divided into K elevation angles. The force F applied to the catheter acts as a scalar distance in the spherical coordinate system.
As a preferred scheme, the cross section of the catheter is equally divided into 12 azimuth angles, n=12, the azimuth angle takes the value of 0 °,30 °,60 °,90 °,120 °,150 °,180 °,210 °,240 °,270 °,300 °,330 °,360 °, and the azimuth angle subareas are seen in the top view direction as shown in fig. 5. The azimuth angle may be equally divided by other values, such as n=6, n=24, n=36, etc., as desired.
As a preferred solution, the angle from the catheter axis to the catheter cross section is equally divided into 4 elevation angles, k=4, the azimuth angle has a value of 0 °,30 °,60 °,90 °, and the elevation angle is divided as shown in fig. 6 when looking in the lateral direction. The azimuth angle may be equally divided by other values, such as n=3, n= 5,N =6, etc., as desired.
As a preferred embodiment, the applied force F to the guide tube has values of 10g,20g,30g,40g,50g,60g,70g, where g is the unit gram, which is the unit of reading of the load cell, 1g corresponds to the weight of 1g of the object, and if g takes 9.8N/Kg, 1g is the unit of reading, which means 9.8X10 -3 N。
For a clearer explanation scheme, the method of the invention will be described below with values of azimuth angle of 0 °,30 °,60 °,90 °,120 °,150 °,180 °,210 °,240 °,270 °,300 °,330 °,360 °, elevation angle of 0 °,30 °,60 °,90 °, and applied forces F by the catheter of 10g,20g,30g,40g,50g,60g,70g, respectively, but the invention is not limited to a scheme that can only be implemented in sequence, but is still within the scope of the invention based on the inventive concept with other values for catheter calibration and catheter pressure testing.
From the azimuth angle beta, the elevation angle alpha and the force F exerted by the catheter, a parameter matrix P can be established, wherein,n is the total number of arrays of parameters in the parameter matrix P based on azimuth, elevation and applied force. For example, the azimuth angle β is 12, the elevation angle α is 4, and the acting force F is 7, and after being arranged and combined, the force application schemes of 12×4×7=364, that is, the force application schemes of 364 arrays are formed, namely, the force application schemes of 364 arrays are formed>
In step S2, the parameter matrix is vectorized, and the calculation method is as follows:
setting the basic parameter vector matrix as R, wherein,the vectors in the basic parameter vector matrix are represented by a first parameter delta i Second parameter epsilon i And a third parameter mu i The composition of the composite material comprises the components,
establishing an array in the parameter matrix and a vector in the basic parameter vector matrix through the subscript iEach parameter matrix array corresponds to the vector of the basic parameter vector matrix, so the same subscript is adopted. When (when)Correspondingly, a->
In step S3, when a plurality of pressure sensors are disposed at the stress end of the catheter, when an acting force is applied to the stress end of the catheter according to the parameter matrix, pressure values of the plurality of pressure sensors can be collected to form a plurality of sets of pressure value arrays, the plurality of pressure value arrays can form a base pressure matrix, and taking a pressure sensing catheter with three pressure sensors as an example, the obtained base pressure matrix is expressed as:when->When corresponding to(α 1 ，β 1 Under the force application scheme of F), the acquired basic pressure array is (x) 1 ，y 1 ，z 1 )，(α 2 ，β 2 ，F 2 ) Under the force application scheme of (2), the acquired basic pressure array is (x) 2 ，y 2 ，z 2 ) And by analogy, establishing a one-to-one correspondence between the array in the parameter matrix and the array in the basic pressure matrix through the subscript i. Therefore, by the subscript i, a one-to-one correspondence is established between the vectors in the base parameter vector matrix and the arrays in the base pressure matrix.
In step S4, a pressure value calibration relation model is established according to the basic pressure matrix and the basic parameter vector matrix, and the method specifically includes the following steps:
according to the basic parameter vector matrix R i And a base pressure matrix Q i Can be calculatedSolving a relation matrix of the two, wherein the relation matrix has a calculation formula as follows:
M i ＝pinv(Q i )*R i
wherein pinv () is an inversion formula, i is 1-n, n is the total number of arrays formed by the parameter matrix P according to azimuth angle, elevation angle and acting force, and is the number of vectors in the basic parameter vector matrix, and is the number of arrays in the basic pressure matrix. M is M i And (5) calibrating the relation model by the established pressure value. Taking a catheter with three pressure sensors at its end as an example, a matrix of base pressuresCorrespondingly, M i Is a3 x 3 matrix.
Matrix of relations M i With base pressure sub-matrix Q i Write memory module, denoted as cell= { M i ，Q i And the output procedure call is supplied and used for solving the stress size and the stress direction of the catheter according to the pressure sensor acquired in real time, wherein i is more than or equal to 1 and less than or equal to n, and n is the total number of the array formed according to azimuth angle, elevation angle and acting force in the parameter matrix P.
Example 2
Example 1 shows a grid partition based catheter calibration method and builds a pressure value calibration relationship model. The method requires a corresponding device to carry out. One of the calibration devices is shown in fig. 7, and includes a bracket, a force applying device movable up and down on the bracket, a force applying magnitude measuring device mounted on the force applying device, and a rocker arm angle and rotation angle control device for controlling the force applying direction.
The force application device module controls the force application device to move up and down or left and right along the support so as to realize the contact and pressure application of the force application platform and the catheter head end, and the force application size measurement device feeds back the size of the applied pressure in real time. The force application device control module realizes accurate force application to the catheter head end according to feedback.
The rocker arm angle and rotation angle control device for controlling the force application direction comprises a force application angle control module, a rocker arm angle control motor and a rotation angle control motor. The force application angle control module realizes accurate control of force application direction through the rocker arm angle control motor and the rotating angle control motor.
The catheter is fixed on the fixture of the force application angle control module, the tail end electrode and the pressure sensor part are exposed, the angle relation between the catheter head end and the force application table, namely, the rocker angle (corresponding to the azimuth angle in the spherical coordinate system) is controlled under the drive of the rocker control motor, and the rotation angle of the catheter head end and the force application table (corresponding to the elevation angle in the spherical coordinate system) is controlled by the rotation angle control motor. The force application detection module is a standard pressure electronic scale and can slide up and down and left and right on the force application device so as to meet the control of the set rocker angle and the rotation angle.
By the calibration device, the front end of the catheter can be applied with force according to the set azimuth angle, elevation angle and acting force, data are collected, and a pressure value calibration relation model is established.
Example 3
Example 3 shows another calibration device, as shown in fig. 8. The device comprises an angle adjusting device, a standard electronic scale and control equipment.
The catheter is fixed on the angle adjusting device, the tail end electrode and the pressure sensor part are exposed, the contact part of the head end of the standard electronic scale device and the catheter is a flat plane, the contact part is recorded when a force is applied to the plane, the included angle (elevation angle) between the tail end electrode of the catheter and the plane of the electronic scale and the rocker angle (azimuth angle) are controlled by the angle adjusting device on the control equipment, the rotation of the catheter, namely the rotation angle between the tail end electrode of the catheter and the plane of the electronic scale, can be controlled by the rotation control device on the control equipment, and the integral movement of the tail end electrode of the catheter can be controlled by the control equipment so as to realize the control of the set angle and pressure.
Example 4
The pressure and direction of the conduit can be obtained by calibrating the relation model according to the obtained pressure value, so a method for detecting the stress of the conduit is provided, and the flow chart is shown in figure 9. The method comprises the following steps:
a1, simultaneously acquiring pressure values of a plurality of pressure sensors to form a pressure acquisition array. Wherein the pressure isThe sensors are arranged at the far end of the catheter as shown in fig. 2, and the sensors are arranged at the far end of the catheter, the embodiment is illustrated by 3 pressure sensors, and a real-time collected pressure collection array is marked as V in ＝(x now ，y now ，z now )。
A2, a pressure value calibration relation model obtained by the calibration method of the embodiment 1 is adopted, and a reference cell relation matrix with the shortest relative distance with the pressure acquisition array is found from the pressure value calibration relation model.
A3, according to the reference cell relation matrix and the pressure acquisition array, the stress magnitude and the stress direction of the catheter are calculated.
Further, step A2 includes the steps of:
a21, obtaining the relative distance between the current pressure acquisition array and the array in the basic pressure matrix.
A22, finding the index number of the corresponding array when the relative distance is minimum;
a23, finding out a corresponding cell relation matrix in the pressure value calibration relation model according to the index number, wherein the cell relation matrix corresponding to the index number is used as the reference cell relation matrix.
In step A21, L is preferably used p Distance (L) p Distance) or Minkowski Distance (Minkowski Distance) to calculate the current pressure acquisition array V in ＝(x now ，y now ，z now ) And a base pressure matrix Q i The relative distance of the array in the middle is L p Distance (Min Shi distance) example, current pressure acquisition array V in ＝(x now ，y now ，z now ) The relative distance to the ith base pressure array is calculated as follows:
wherein,is the current pressure acquisition array V in ＝(x now ，y now ，z now ) The relative distance from the i-th base pressure array is 1.ltoreq.i.ltoreq.n, where n is the base pressure matrix Q i The total number of arrays, p, is a constant, and a slice may take 1 or 2.
In step A22, the relative distances are comparedFinding out the minimum value of the relative distance and the index number i corresponding to the minimum value min 。
In step A23, index number i corresponding to the minimum relative distance min Finding a relationship matrix in a pressure value calibration relationship modelThe corrected conduit pressure gamma can be obtained by using the following formula out And direction (by rocker angle alpha) out Rotation angle beta out Representation):
wherein, gamma out Is the output conduit pressure magnitude, alpha out Is the output rocker angle beta out Is the rotation angle of the output and,minimum relative distance index number i min Corresponding relation matrix, x out 、y out And z out Is via a relation matrix->Vector values of the three corrected pressure sensors.
With the above calibration device and calibration scheme, 100 groups are randomly verified, and the effect is as shown in fig. 10, the pressure deviation is the true value minus the system output pressure value, the rocker arm angle deviation is the true value minus the system output rocker arm angle value, and the rotation angle deviation is the true value minus the system output rotation angle value.
The foregoing description of the preferred embodiments of the invention is not intended to limit the invention to the particular embodiments disclosed, but on the contrary, the intention is to cover all modifications, equivalents, and alternatives falling within the spirit and scope of the invention.
Claims (9)
1. A grid partition-based catheter calibration method for a catheter with a plurality of pressure sensors at a stress end, comprising the following steps:
s1, establishing a parameter matrix at a catheter stress end according to a preset azimuth angle, an elevation angle and an acting force; vectorizing the parameter matrix, and establishing a basic parameter vector matrix;
s2, applying acting force to the stress end of the catheter according to the parameter matrix, collecting the pressure value of the pressure sensor at the same time, and constructing a basic pressure matrix, wherein an array in the basic pressure matrix and a vector in the basic parameter vector matrix are in one-to-one correspondence;
s3, establishing a pressure value calibration relation model according to the basic pressure matrix and the basic parameter vector matrix, wherein the pressure value calibration relation model is used for solving the stress magnitude and the stress direction of the catheter according to the value of the pressure sensor acquired in real time;
the pressure value calibration relation model is as follows:
M i ＝pinv(Q i )*R i
wherein M is i Cell relation matrix of pressure value calibration relation model, R i Is a basic parameter vector matrix, Q i The system is a basic pressure matrix, wherein an array in the basic pressure matrix is a pressure value array corresponding to a vector in the basic parameter vector matrix, pinv () is an inversion function, i is more than or equal to 1 and less than or equal to n, and n is the total number of vectors in the basic parameter vector matrix.
2. The grid partition based catheter calibration method of claim 1, wherein the parameter matrix establishment process of step S1 is: the center of the cross section of the catheter is taken as the origin of spherical coordinates, the cross section of the catheter is equally divided into N azimuth angles, the angle from the axis of the catheter to the cross section of the catheter is equally divided into K elevation angles, and the acting force in the direction of the axis of the catheter is a scalar distance.
3. A grid partition based catheter calibration method as recited in claim 2, wherein N is 12 and the azimuth is: 0 °,30 °,60 °,90 °,120 °,150 °,180 °,210 °,240 °,270 °,300 °,330 °,360 °.
4. A grid partition based catheter calibration method as recited in claim 3, wherein K is 4 and said elevation angle is: 0 °,30 °,60 °,90 °.
5. The grid partition based catheter calibration method of claim 4, wherein the force applied in the axial direction of the catheter is 10g,20g,30g,40g,50g,60g,70g.
6. The grid partition based catheter calibration method of claim 1, wherein in step S1, the vectors in the base parameter vector matrix are defined by a first parameter δ i Second parameter epsilon i And a third parameter mu i The composition of the composite material comprises the components,
first ginseng (radix Ginseng)Number delta i The calculation formula of (2) is as follows:
second parameter epsilon i The calculation formula of (2) is as follows:
third parameter mu i The calculation formula of (2) is as follows:
wherein alpha is i Is the elevation angle, beta i The azimuth angle is F, the acting force is 1.ltoreq.i.ltoreq.n, and n is the total number of arrays formed according to the azimuth angle, the elevation angle and the acting force in the parameter matrix.
7. A method for detecting the stress of a catheter, which is used for the catheter with a plurality of pressure sensors distributed at the stress end, and is characterized by comprising the following steps:
a1, simultaneously acquiring pressure values of a plurality of pressure sensors to form a pressure acquisition array;
a2, a pressure value calibration relation model obtained by the calibration method according to any one of claims 1-6 is adopted, and a reference cell relation matrix with the shortest relative distance with the pressure acquisition array is found from the pressure value calibration relation model;
a3, according to the reference cell relation matrix and the pressure acquisition array, the stress magnitude and the stress direction of the catheter are calculated.
8. The method of claim 7, wherein step A2 comprises the steps of:
a21, solving the relative distance between the current pressure acquisition array and the array in the basic pressure matrix;
a22, finding the index number of the corresponding array when the relative distance is minimum;
a23, finding out a corresponding reference cell relation matrix in the pressure value calibration relation model according to the index number.
9. A catheter calibration device based on grid partition is characterized by comprising a force application device, a force application size measuring device and an azimuth angle and elevation angle control device for controlling the force application direction, wherein the force application device is used for applying pressure to a stress end of a catheter,
the force application measuring device is used for collecting the pressure value applied to the stress end of the guide tube in real time,
the azimuth angle and elevation angle control device for controlling the force application direction is used for controlling the force application device to apply pressure to the force receiving end of the catheter according to the preset azimuth angle and elevation angle,
the force application device, the force application magnitude measuring device and the azimuth and elevation control device for controlling the force application direction apply force to the force receiving end of the catheter according to the method according to any one of claims 1-6 according to the parameter matrix.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011035782.4A CN114271924B (en) | 2020-09-27 | 2020-09-27 | Catheter calibration method and device based on grid partition |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011035782.4A CN114271924B (en) | 2020-09-27 | 2020-09-27 | Catheter calibration method and device based on grid partition |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN114271924A CN114271924A (en) | 2022-04-05 |
| CN114271924B true CN114271924B (en) | 2024-03-01 |
Family
ID=80867748
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011035782.4A Active CN114271924B (en) | 2020-09-27 | 2020-09-27 | Catheter calibration method and device based on grid partition |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114271924B (en) |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6066125A (en) * | 1997-09-05 | 2000-05-23 | Cordis Webster, Inc. | Omni-directional steerable catheter |
| CN102510735A (en) * | 2009-07-17 | 2012-06-20 | 计算机心脏股份有限公司 | Heart treatment kit, system, and method for radiosurgically alleviating arrhythmia |
| CN103908337A (en) * | 2012-12-31 | 2014-07-09 | 韦伯斯特生物官能(以色列)有限公司 | Catheter with serially connected sensing structures and methods of calibration and detection |
| CN106606373A (en) * | 2016-08-30 | 2017-05-03 | 四川锦江电子科技有限公司 | Ablation catheter capable of accurate pressure measurement |
| CN111084659A (en) * | 2018-10-24 | 2020-05-01 | 韦伯斯特生物官能(以色列)有限公司 | Instant calibration for catheter position and orientation |
| CN111345901A (en) * | 2018-12-20 | 2020-06-30 | 四川锦江电子科技有限公司 | Method for measuring stress of radio frequency ablation pressure catheter |
-
2020
- 2020-09-27 CN CN202011035782.4A patent/CN114271924B/en active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6066125A (en) * | 1997-09-05 | 2000-05-23 | Cordis Webster, Inc. | Omni-directional steerable catheter |
| CN102510735A (en) * | 2009-07-17 | 2012-06-20 | 计算机心脏股份有限公司 | Heart treatment kit, system, and method for radiosurgically alleviating arrhythmia |
| CN103908337A (en) * | 2012-12-31 | 2014-07-09 | 韦伯斯特生物官能(以色列)有限公司 | Catheter with serially connected sensing structures and methods of calibration and detection |
| CN106606373A (en) * | 2016-08-30 | 2017-05-03 | 四川锦江电子科技有限公司 | Ablation catheter capable of accurate pressure measurement |
| CN111084659A (en) * | 2018-10-24 | 2020-05-01 | 韦伯斯特生物官能(以色列)有限公司 | Instant calibration for catheter position and orientation |
| CN111345901A (en) * | 2018-12-20 | 2020-06-30 | 四川锦江电子科技有限公司 | Method for measuring stress of radio frequency ablation pressure catheter |
Also Published As
| Publication number | Publication date |
|---|---|
| CN114271924A (en) | 2022-04-05 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| EP2438881B1 (en) | Calibration system for a force-sensing catheter | |
| EP2394578B1 (en) | Weight-based calibration system for a pressure sensitive catheter | |
| CN102119871B (en) | Actuator-based calibration system for a pressure-sensitive catheter | |
| CN109141697B (en) | Strain gauge, pressure sensor and interventional medical catheter | |
| EP2338411B1 (en) | Calibration system for a pressure-sensitive catheter | |
| JP2006064465A (en) | 3-d drag sensor | |
| KR20150060949A (en) | System and method for three dimensional calibration of force plates | |
| CN106606373A (en) | Ablation catheter capable of accurate pressure measurement | |
| CN114271924B (en) | Catheter calibration method and device based on grid partition | |
| CN1325134C (en) | System and method for deriving angular isokinetic measurements using a linear dynamometer | |
| CN105232141A (en) | Medical catheter head end capable of measuring contact force | |
| CN116673796A (en) | Calibration tool and calibration method for robot hole making system | |
| CN113483816B (en) | Shape-position-force composite sensing unit and measuring method thereof | |
| JP4845734B2 (en) | Coordinate measuring device, method, computer program and data carrier | |
| CN116115338A (en) | Handheld probe calibration device and calibration method | |
| CN111345901B (en) | Method for measuring stress of radio frequency ablation pressure catheter | |
| CN106859762B (en) | Electrophysiology catheter | |
| JP2021177177A (en) | Catheter shape and position detection using flexible magnetic sensor | |
| Jiang et al. | Analysis of a 3-Axis Force Sensor for Minimally Invasive Surgery | |
| CN111968187B (en) | Annular structure optical parameter calibration method | |
| CN112236100B (en) | Method for calibrating a torque sensor of an instrument drive unit of a surgical robot | |
| CN109984838B (en) | Take electrophysiology pipe of pressure perception | |
| CN114593848A (en) | Wide-range high-sensitivity resistance strain measurement system | |
| WO2021222471A2 (en) | System for substrate detection and mapping using force sensing and magnetic coupling | |
| CN118776718A (en) | Puncture needle shape and interaction force sensing sensor and calibration method thereof |
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 | ||
| CB02 | Change of applicant information |
Address after: No. 5, Wuke East 3rd Road, Wuhou District, Chengdu, Sichuan 610000 Applicant after: Sichuan Jinjiang Electronic Medical Device Technology Co.,Ltd. Address before: No.5, Wuke East 3rd road, Wuhou Science Park, Chengdu hi tech Industrial Development Zone, Sichuan 610045 Applicant before: SICHUAN JINJIANG ELECTRONIC SCIENCE AND TECHNOLOGY Co.,Ltd. |
|
| GR01 | Patent grant | ||
| GR01 | Patent grant |