CN116458865A - Electrical impedance imaging fault boundary extraction method based on multi-origin transformation - Google Patents
Electrical impedance imaging fault boundary extraction method based on multi-origin transformation Download PDFInfo
- Publication number
- CN116458865A CN116458865A CN202310402082.1A CN202310402082A CN116458865A CN 116458865 A CN116458865 A CN 116458865A CN 202310402082 A CN202310402082 A CN 202310402082A CN 116458865 A CN116458865 A CN 116458865A
- Authority
- CN
- China
- Prior art keywords
- boundary
- field
- fault
- transformation
- coordinate system
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 230000009466 transformation Effects 0.000 title claims abstract description 62
- 238000003384 imaging method Methods 0.000 title claims abstract description 21
- 238000000605 extraction Methods 0.000 title claims description 16
- 239000013598 vector Substances 0.000 claims abstract description 19
- 238000000034 method Methods 0.000 claims abstract description 14
- 238000003325 tomography Methods 0.000 claims abstract description 3
- 238000002593 electrical impedance tomography Methods 0.000 claims description 12
- 238000005259 measurement Methods 0.000 claims description 8
- 230000008859 change Effects 0.000 claims description 2
- 238000004422 calculation algorithm Methods 0.000 abstract description 4
- 238000010586 diagram Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000005284 excitation Effects 0.000 description 2
- 230000001788 irregular Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 230000001131 transforming effect Effects 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 238000002059 diagnostic imaging Methods 0.000 description 1
- 201000010099 disease Diseases 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000009545 invasion Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 230000001575 pathological effect Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/05—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves
- A61B5/053—Measuring electrical impedance or conductance of a portion of the body
- A61B5/0536—Impedance imaging, e.g. by tomography
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/0033—Features or image-related aspects of imaging apparatus classified in A61B5/00, e.g. for MRI, optical tomography or impedance tomography apparatus; arrangements of imaging apparatus in a room
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
- Y04S—SYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
- Y04S10/00—Systems supporting electrical power generation, transmission or distribution
- Y04S10/50—Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications
- Y04S10/52—Outage or fault management, e.g. fault detection or location
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Biomedical Technology (AREA)
- Medical Informatics (AREA)
- Biophysics (AREA)
- Pathology (AREA)
- Engineering & Computer Science (AREA)
- Radiology & Medical Imaging (AREA)
- Heart & Thoracic Surgery (AREA)
- Physics & Mathematics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
The method establishes a reference coordinate system on a two-dimensional plane where a target field fault is located, then determines a plurality of transformation origins located outside the target field fault boundary on the two-dimensional plane and coordinates thereof in the reference coordinate system to obtain transformation vectors from any transformation origin to corresponding arbitrary field boundary points, and then combines the coordinates of the transformation origin in the reference coordinate system to obtain the coordinates of the field boundary points in the reference coordinate system, thereby being capable of extracting the fault boundary of the target field fault. The method has simple mathematical algorithm and strong operability, and can accurately establish the two-dimensional field boundary of the object tomography in the three-dimensional object by the electrical impedance imaging, thereby improving the quality of the reconstructed image.
Description
Technical Field
The application relates to the technical field of electrical impedance imaging, in particular to an electrical impedance imaging fault boundary extraction method based on multi-origin transformation.
Background
The electrical impedance imaging technology (Electrical Impedance Tomography, EIT) is an emerging medical imaging technology, has the advantages of safety, no invasion, simple system, low cost and the like, and has wide application prospects in aspects of disease diagnosis, pathological state monitoring and the like. In recent years, EIT application research has also been expanded to some industrial detection fields. And attaching an EIT electrode to the boundary of the target field, injecting a certain current into the boundary of the target field through the EIT electrode to serve as excitation, measuring the boundary voltage response, and processing excitation and measurement data by using an image reconstruction algorithm to reconstruct the resistivity distribution imaging of the target field. The structure and the function of the target can be evaluated and judged through the comparison analysis of the resistivity distribution image and the historical data.
In electrical impedance imaging, the accuracy of a field physical model is an important factor affecting imaging quality, and in order to improve imaging accuracy, the boundary of an imaging field should be accurately extracted. The field of two-dimensional electrical impedance imaging is a fault plane (cross section) of a three-dimensional object of a target, in simulation research, a regular circle or ellipse is generally used for simulating an imaging field, but the actual imaging target cross section is often irregular, such as a trunk cross section of a human body, so that the extraction of boundaries of the irregular imaging fields is difficult to realize at present, and because the fault field cannot be acquired by transversely cutting the three-dimensional object in practice, an accurate field physical model is difficult to establish.
Disclosure of Invention
Aiming at the problems and the technical requirements, the applicant provides an electrical impedance imaging fault boundary extraction method based on multi-origin transformation, and the technical scheme of the application is as follows:
the electrical impedance imaging fault boundary extraction method based on multi-origin transformation comprises the following steps:
establishing a reference coordinate system on a two-dimensional plane where the target field fault is located;
determining coordinates of a plurality of transformation origins which are positioned on a two-dimensional plane and outside a fault boundary of the target field fault and are in a reference coordinate system, wherein the plurality of transformation origins comprise origins of the reference coordinate system;
determining an arbitrary nth transformation origin O Tn To the corresponding arbitrary mth field boundary point B n m Is of the transform vector of (a)Any field boundary point is positioned on the fault boundary of the target field fault;
according to arbitrary transformation vectorsCombining with transformation origin O Tn Coordinates in the reference coordinate system, resulting in a field boundary point +.>And extracting coordinates in a reference coordinate system to obtain a fault boundary of the target field fault.
The further technical proposal is that a connecting line positioned outside the fault boundary is arranged between each transformation origin and the corresponding field boundary point, and the formed connecting line does not pass through other field boundary points.
The further technical proposal is that the fault boundary of the target field fault comprises a plurality of field boundary points at discrete intervals, and each field boundary point has a corresponding transformation origin.
Further technical proposal is that the transformation vector is determinedThe method of (1) comprises:
measurement transformation origin O Tn To the boundary point of the fieldStraight line distance +.>
Measurement transformation origin O Tn To the boundary point of the fieldAn included angle theta between the connecting line of the (E) and the positive direction of the x-axis of the reference coordinate system;
determining a transformation vectorExpressed as +.>
The further technical proposal is that the boundary point of the field is obtainedThe coordinates in the reference coordinate system areWherein, (x) Tn ,y Tn ) Is to change the origin O Tn Coordinates in a reference coordinate system.
The further technical proposal is that the method for extracting the fault boundary of the target field fault comprises the following steps:
and sequentially connecting all the field boundary points according to the obtained coordinates of all the field boundary points in a reference coordinate system to obtain the fault boundary of the target field fault on the two-dimensional plane.
The beneficial technical effects of this application are:
the method is simple in implementation method, can accurately establish a two-dimensional field boundary of target tomography in a three-dimensional object by electrical impedance imaging, can reduce imaging algorithm errors caused by field subdivision coordinate deviation based on the accurate fault boundary, reduces artifact generation and improves quality of reconstructed images.
When the method is realized, the transformation origin is flexible to select, only the principle that the boundary points of the field can be reached once and each boundary point of the field has a corresponding transformation origin is met, the transformation origin can be selected at any position meeting the condition without influencing the result, in the field boundary extraction, the module length and the included angle with the OX direction of the transformation vector of the origin and the transformation vector of the boundary point are directly measured, each vector can be obtained, and the extraction of the coordinate data of the boundary point is completed through a mathematical algorithm, so that the operation is easy to implement.
Drawings
FIG. 1 is a method flow diagram of a method of electrical impedance tomography boundary extraction in accordance with one embodiment of the present application.
Fig. 2 is a perspective view and a plan cross-sectional view of a target field slice in one example of the present application.
Fig. 3 is a vector diagram of the calculation process of obtaining the coordinates of a corresponding field boundary point based on a transformation origin in fig. 2.
Detailed Description
The following describes the embodiments of the present application further with reference to the accompanying drawings.
The application discloses an electrical impedance tomography fault boundary extraction method based on multi-origin transformation, please refer to a method flow chart shown in fig. 1, the method comprises the following steps:
step 1, a reference coordinate system is established on a two-dimensional plane alpha where a target field fault is located, wherein the target field fault is a field fault needing to be subjected to electrical impedance imaging, and is generally a fault plane of a three-dimensional target body, as shown in fig. 2, one fault plane of the three-dimensional target body 1 forms a target field fault 2, and a fault boundary of the target field fault 2 comprises a plurality of field boundary points 3 at discrete intervals.
And taking any point outside the fault boundary of the target field fault as an original point O (0, 0) on the two-dimensional plane, and establishing a reference coordinate system XOY.
And 2, determining coordinates of a plurality of transformation origins which are positioned on the two-dimensional plane alpha and outside the fault boundary of the target field fault and are in a reference coordinate system, wherein the plurality of transformation origins comprise origins O (0, 0) of the reference coordinate system. Any nth transformation origin O obtained by determination of the application Tn The coordinates in the reference coordinate system are (x Tn ,y Tn ) N is a parameter with a starting value of 1, FIG. 3 uses the origin O (0, 0) of the reference coordinate system as the first transformed origin O T1 。
The number of transformation origins and the positions of the transformation origins are adjusted according to the form and the complex condition of the fault boundary of the target field fault, and the principle of determining the transformation origins comprises the following steps:
the connecting line positioned outside the fault boundary is arranged between each transformation origin and the corresponding field boundary point, and the formed connecting line does not pass through other field boundary points, namely, the connecting line can be reached once. In addition, the fault boundary of the target field fault comprises a plurality of field boundary points at discrete intervals, and each field boundary point has a corresponding transformation origin. Any one field boundary point has only one corresponding transformation origin or a plurality of corresponding transformation origins.
Based on the principle of determining the transformation origins, in the example of fig. 2, 5 transformation origins are determined and obtained, respectively, as transformation origins O T1 、O T2 、O T3 、O T4 And O T5 Transform origin O T1 Corresponding to 6 field boundary points, transform origin O T2 Corresponding to 5 field boundary points, transforming origin O T3 Corresponding to 3 field boundary points, transform origin O T4 Corresponding to 4 field boundary points, transforming origin O T5 Corresponding to 3 field boundary points. Of the 16 field boundary points on the target field fault 2, 11 field boundary points each have one corresponding transformation origin, and the other 5 field boundary points each have two corresponding transformation origins.
Step 3, determining any oneMeaning nth transform origin O Tn To the corresponding arbitrary mth field boundary pointIs of the transform vector of (a)m is a parameter with a start value of 1. Comprising the following steps:
(1) Measurement transformation origin O Tn To the boundary point of the fieldStraight line distance +.>
(2) Measurement transformation origin O Tn To the boundary point of the fieldAn included angle theta between the line of the reference frame and the positive x-axis direction.
(3) Determining a transformation vectorExpressed as +.>
Step 4, according to any transformation vectorCombining with transformation origin O Tn Coordinates in the reference coordinate system, resulting in a field boundary point +.>Coordinates in a reference coordinate system.
At the determination of the transformation origin O Tn Coordinates in a reference coordinate system (x Tn ,y Tn ) Then, the origin O (0, 0) of the obtained reference coordinate system is equal toTransform origin O Tn Vector of (3)By the principle of vector addition, the origin O (0, 0) of the reference coordinate system can be determined to the field boundary point +.>Vector of->The boundary point of the field is obtainedThe coordinates in the reference coordinate system are +.>
Such as in fig. 3, to transform the origin O T2 Corresponding field boundary pointsFor example, the measurement transform origin O T2 To the boundary point of the field->And an angle theta, thereby determining a transformation vector +.>Combining with transformation origin O T2 Coordinates in the reference coordinate system can be used to obtain the field boundary point +.>Coordinates in a reference coordinate system. The coordinates of other field boundary points in the reference coordinate system are equally available.
And step 5, after the coordinates of all the field boundary points in the reference coordinate system are obtained, the fault boundary of the target field fault can be extracted. And smoothly connecting all the field boundary points in sequence according to the obtained coordinates of all the field boundary points in a reference coordinate system to obtain the fault boundary of the target field fault on the two-dimensional plane.
What has been described above is only a preferred embodiment of the present application, which is not limited to the above examples. It is to be understood that other modifications and variations which may be directly derived or contemplated by those skilled in the art without departing from the spirit and concepts of the present application are to be considered as being included within the scope of the present application.
Claims (6)
1. The electrical impedance imaging fault boundary extraction method based on multi-origin transformation is characterized by comprising the following steps of:
establishing a reference coordinate system on a two-dimensional plane where the target field fault is located;
determining coordinates of a plurality of transformation origins which are positioned outside a fault boundary of the target field fault and on the two-dimensional plane in the reference coordinate system, wherein the plurality of transformation origins comprise origins of the reference coordinate system;
determining an arbitrary nth transformation origin O Tn To the corresponding arbitrary mth field boundary pointIs>Any field boundary point is positioned on the fault boundary of the target field fault;
according to arbitrary transformation vectorsCombining with transformation origin O Tn Coordinates in the reference coordinate system, a field boundary point +.>Coordinates in the reference coordinate system and extracting the target fieldFault boundaries of faults.
2. The electrical impedance tomography boundary extraction method of claim 1 wherein each transformation origin and its corresponding field boundary point have a connection line located outside the tomography boundary, and the connection line is formed without passing other field boundary points, referred to as single pass.
3. The electrical impedance tomography boundary extraction method of claim 1 wherein the target field fault comprises a plurality of discrete spaced field boundary points on the fault boundary, each field boundary point having a corresponding transformation origin.
4. The electrical impedance tomography boundary extraction method of claim 1 wherein a transformation vector is determinedThe method of (1) comprises:
measurement transformation origin O Tn To the boundary point of the fieldStraight line distance +.>
Measurement transformation origin O Tn To the boundary point of the fieldAn included angle theta between the connecting line of the reference coordinate system and the positive direction of the x-axis;
determining a transformation vectorExpressed as +.>
5. The electrical impedance tomography boundary extraction method of claim 4, wherein a field boundary point is obtainedThe coordinates in the reference coordinate system are +.>Wherein, (x) Tn ,y Tn ) Is to change the origin O Tn Coordinates in the reference coordinate system.
6. The electrical impedance tomography fault boundary extraction method of claim 1, wherein the method of extracting the fault boundary of the target field fault comprises:
and sequentially connecting all the field boundary points according to the obtained coordinates of all the field boundary points in the reference coordinate system to obtain the fault boundary of the target field fault on the two-dimensional plane.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310402082.1A CN116458865B (en) | 2023-04-14 | 2023-04-14 | Electrical impedance imaging fault boundary extraction method based on multi-origin transformation |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310402082.1A CN116458865B (en) | 2023-04-14 | 2023-04-14 | Electrical impedance imaging fault boundary extraction method based on multi-origin transformation |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN116458865A true CN116458865A (en) | 2023-07-21 |
| CN116458865B CN116458865B (en) | 2024-02-02 |
Family
ID=87180149
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202310402082.1A Active CN116458865B (en) | 2023-04-14 | 2023-04-14 | Electrical impedance imaging fault boundary extraction method based on multi-origin transformation |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN116458865B (en) |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20170172451A1 (en) * | 2015-12-16 | 2017-06-22 | General Electric Company | System and method for enhanced electrical impedance tomography |
| CN107224287A (en) * | 2017-03-31 | 2017-10-03 | 重庆邮电大学 | Flexible sensor electrode position optimization method based on Review of Electrical Impedance Tomography |
| CN112450910A (en) * | 2020-12-04 | 2021-03-09 | 桂林电子科技大学 | Boundary measuring device and method for electrical impedance tomography |
| CN114847913A (en) * | 2022-04-14 | 2022-08-05 | 四川大学华西医院 | Bioelectrical impedance tomography device and method |
-
2023
- 2023-04-14 CN CN202310402082.1A patent/CN116458865B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20170172451A1 (en) * | 2015-12-16 | 2017-06-22 | General Electric Company | System and method for enhanced electrical impedance tomography |
| CN107224287A (en) * | 2017-03-31 | 2017-10-03 | 重庆邮电大学 | Flexible sensor electrode position optimization method based on Review of Electrical Impedance Tomography |
| CN112450910A (en) * | 2020-12-04 | 2021-03-09 | 桂林电子科技大学 | Boundary measuring device and method for electrical impedance tomography |
| CN114847913A (en) * | 2022-04-14 | 2022-08-05 | 四川大学华西医院 | Bioelectrical impedance tomography device and method |
Non-Patent Citations (1)
| Title |
|---|
| 王琦 等: "基于人体结构先验信息的胸部电阻抗成像方法", 《中国生物医学工程学报》, vol. 38, no. 1, pages 35 - 43 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN116458865B (en) | 2024-02-02 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| WO2020134826A1 (en) | Parallel magnetic resonance imaging method and related equipment | |
| CN110251130B (en) | Electrical impedance imaging method, device and system based on movable deformation component | |
| CN112381862B (en) | Full-automatic registration method and device for CAD (computer-aided design) model and triangular mesh | |
| CN111753947B (en) | Resting brain network construction method, device, equipment and computer storage medium | |
| CN112006672A (en) | Determining occurrence of focal and/or trochanteric arrhythmogenic activity in a cardiac tissue region | |
| CN111640143A (en) | Nerve navigation rapid surface registration method and system based on PointNet | |
| CN110414481A (en) | A kind of identification of 3D medical image and dividing method based on Unet and LSTM | |
| CN1679030A (en) | Method and system for displaying confidence intervals for source reconstruction | |
| CN116458865B (en) | Electrical impedance imaging fault boundary extraction method based on multi-origin transformation | |
| JP2007167374A (en) | Medical image processor | |
| CN112669386A (en) | Magnetoencephalogram automatic positioning and registering method and device based on three-dimensional optical scanning | |
| CN113948189B (en) | MEG source positioning method based on GRU neural network | |
| CN107102056B (en) | Unit expansion construction method of defect magnetic leakage signal | |
| CN114748099B (en) | Offset detection method, apparatus, device, and computer-readable storage medium | |
| US20220328195A1 (en) | System and method for heart strain determination | |
| CN115830156A (en) | Accurate electrical impedance tomography method, apparatus, system, medium and device | |
| CN111681292B (en) | Task fMRI brain decoding and visualization method based on convolutional neural network | |
| CN114755010A (en) | Rotary machine vibration fault diagnosis method and system | |
| CN112766313A (en) | Crystal segmentation and positioning method, device, equipment and medium based on U-net structure | |
| Tang et al. | Measurement of displacement in isolated heart muscle cells using markerless subpixel image registration | |
| Krinidis et al. | Reconstruction of serially acquired slices using physics-based modeling | |
| CN114757951B (en) | Sign data fusion method, data fusion equipment and readable storage medium | |
| CN101794373B (en) | Application method of rotating and sub-pixel matching algorithm to machine vision system | |
| CN117274413B (en) | EIT-based conductivity image reconstruction method, system and equipment | |
| CN116309591B (en) | Medical image 3D key point detection method, model training method and device |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |