CN113081258B - Optimal point calibration method for puncturing effusion drainage in joint cavity treatment - Google Patents

Optimal point calibration method for puncturing effusion drainage in joint cavity treatment Download PDF

Info

Publication number
CN113081258B
CN113081258B CN202110257430.1A CN202110257430A CN113081258B CN 113081258 B CN113081258 B CN 113081258B CN 202110257430 A CN202110257430 A CN 202110257430A CN 113081258 B CN113081258 B CN 113081258B
Authority
CN
China
Prior art keywords
regular
optimal
puncture
expressed
joint cavity
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
Application number
CN202110257430.1A
Other languages
Chinese (zh)
Other versions
CN113081258A (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.)
Yangzhou University
Original Assignee
Yangzhou University
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 Yangzhou University filed Critical Yangzhou University
Priority to CN202110257430.1A priority Critical patent/CN113081258B/en
Publication of CN113081258A publication Critical patent/CN113081258A/en
Application granted granted Critical
Publication of CN113081258B publication Critical patent/CN113081258B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B34/00—Computer-aided surgery; Manipulators or robots specially adapted for use in surgery
    • A61B34/10—Computer-aided planning, simulation or modelling of surgical operations
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B17/00—Surgical instruments, devices or methods, e.g. tourniquets
    • A61B17/34—Trocars; Puncturing needles

Abstract

The invention discloses a method for calibrating an optimal point of a puncture effusion drainage for joint cavity treatment, which comprises the following steps: s1, acquiring image information of an original lesion part; s2, decomposing mechanistic data based on the lesion part image; s3, positioning partition body selection parameters based on pathological change position data; s4, determining an optimal puncture liquid extraction point based on the selection parameters of the partition body; the method has high detection precision, improves the efficiency of extracting the effusion, and has important practical significance for the treatment of the effusion extraction by joint cavity puncture.

Description

Method for calibrating optimal point of hydrops aspiration by puncture in joint cavity treatment
Technical Field
The invention relates to a data processing method, in particular to a method for calibrating an optimal point of a puncture effusion extraction for joint cavity treatment.
Background
Nowadays, the society is more and more serious in aging, and the incidence of joint cavity effusion is increased year by year. However, the joint cavity treatment position space is small and the structure is complex, so that the needle can be quickly and accurately inserted into a narrow joint gap with great difficulty, and a doctor is exposed to CT irradiation for a long time, which inevitably causes certain influence on the health of the doctor. However, the traditional method of manual puncture by doctors under CT guidance is adopted for treating joint cavity effusion at home and abroad at present, and an optimal puncture site calibration method does not exist. Therefore, it is necessary to determine the optimal puncture site for the treatment of effusion in the joint cavity.
Disclosure of Invention
The invention aims to: the invention aims to provide a method for calibrating an optimal point for puncturing and effusion pumping in joint cavity treatment, which has high detection precision.
The technical scheme is as follows: the invention provides a method for calibrating an optimal point of a puncture effusion drainage for joint cavity treatment, which comprises the following steps:
s1, acquiring image information of an original lesion part:
selecting a puncture part commonly used for joint puncture points, carrying out axial plane and longitudinal section CT scanning on the part, reconstructing an obtained image into a three-dimensional space three-dimensional structure and determining the position of each pixel point in a coordinate system, wherein x is the horizontal direction of the axial plane and the longitudinal section, y is the vertical direction of the axial plane, z is the vertical direction of the longitudinal section, and o is an origin.
S2, decomposing mechanistic data based on the lesion site image:
and (3) dividing the space polyhedron obtained in the step (S1) along the crossed connecting line of the boundary characteristic points to obtain alpha regular tetrahedrons, beta regular hexahedrons, gamma regular octahedrons, delta regular dodecahedrons and e regular icosahedrons.
S3, positioning segmentation body selection parameters based on pathological change position data:
let four vertexes of one regular tetrahedron of the alpha regular tetrahedrons be (x) 1 ，y 1 ，z 1 )，(x 2 ，y 2 ，z 2 )，(x 3 ，y 3 ，z 3 )，(x 4 ，y 4 ，z 4 ) For this regular tetrahedron, the preferred parameter is the locus designation U, whose coordinates can be expressed as:
Figure BDA0002966859690000021
for a regular tetrahedrons, the preferred parameter is the locus designation V, whose coordinates can be expressed as:
Figure BDA0002966859690000022
wherein k is x ，k y ，k z And the scale factor is represented, and is the ratio of the x, y and z dimensional calibration factors corresponding to the optimal puncture sites of the alpha regular tetrahedrons to the dimensional calibration factor of the tetrahedron, and preferably, the value is 1.10 usually.
Similarly, for a regular hexahedron with a length a, the preferred parameter is the position index N, and the coordinates of N can be expressed as:
Figure BDA0002966859690000023
for the β regular hexahedrons, the preferred parameter is the site designation R, whose coordinates can be expressed as:
Figure BDA0002966859690000024
wherein l x ，l y ，l z The representative proportional coefficient is the ratio of the x, y and z dimensional calibration coefficients corresponding to the beta regular hexahedron optimal puncture sites to the calibration coefficient of the regular hexahedron dimension with the edge length of a, and preferably, the normal value is 1.07.
Similarly, when there are γ octahedra, for each octahedra, the preferred parameter is the locus number Q, whose coordinates can be expressed as:
Figure BDA0002966859690000031
wherein j is x ，j y ，j z Represents a proportionality coefficient, preferably, the usual value is 1.11.
Similarly, for a regular dodecahedron with a ridge length of b, the preferred parameter is the locus designation D, whose coordinates can be expressed as:
Figure BDA0002966859690000032
for η regular dodecahedrons, the preferred parameter is the site designation T, whose coordinates can be expressed as:
Figure BDA0002966859690000033
wherein s is x ，s y ，s z And the scaling factor is represented, and is the ratio of the calibration factors of x, y and z corresponding to the optimal puncture sites of the eta regular dodecahedron to the calibration factor of the dimension of the regular dodecahedron with the edge length of b, and preferably, the value is 1.19 usually.
Similarly, for a regular icosahedron with a ridge length c, the preferred parameter is the site designation E, whose coordinates can be expressed as:
Figure BDA0002966859690000034
for μ regular icosahedrons, the preferred parameter is the site designation Y, whose coordinates can be expressed as:
Figure BDA0002966859690000041
wherein, t x ，t y ，t z And the scale factor is represented, and is the ratio of the calibration coefficients of the x, y and z dimensions corresponding to the optimal puncture site of the mu regular icosahedron to the calibration coefficient of the dimension of the regular icosahedron with the edge length of c, and preferably, the value is 1.13.
S4, determining an optimal puncture liquid extraction point based on the selection parameters of the partition body:
and (4) calculating the optimal puncture site for treating the joint cavity effusion according to the optimal parameters of the five regular polyhedrons obtained in the step (S3). Through calculation, the optimal puncture site coordinate is a matrix M
Figure BDA0002966859690000042
Has the beneficial effects that: according to the invention, the acquired image information of the original lesion site can be analyzed, and then the division body is decomposed by the mechanical data, so that the optimal puncture site of the division body is obtained, and the optimal point for puncturing and effusion extraction is determined according to the optimal puncture site, so that the puncture position can be conveniently determined, and the effusion extraction efficiency is improved.
Drawings
FIG. 1 is a block diagram of the process flow of the present invention.
Detailed Description
As shown in fig. 1, the method for calibrating the optimal point of the effusion puncturing during the joint cavity treatment in the embodiment comprises the following steps:
s1, acquiring image information of an original lesion part:
selecting a puncture part commonly used for joint puncture points, carrying out axial plane and longitudinal section CT scanning on the puncture part, reconstructing an obtained image into a three-dimensional space structure and determining the position of each pixel point in a coordinate system, wherein x is the horizontal direction of an axial plane and the longitudinal section, y is the vertical direction of the axial plane, z is the vertical direction of the longitudinal section, and o is an origin.
S2, decomposing mechanistic data based on a lesion part image:
and (3) dividing the space polyhedron obtained in the step (S1) along the cross connecting line of the boundary characteristic points to obtain alpha regular tetrahedrons, beta regular hexahedrons, gamma regular octahedrons, delta regular dodecahedrons and epsilon regular icosahedrons.
S3, positioning segmentation body selection parameters based on pathological change position data:
let four vertexes of one regular tetrahedron of the alpha regular tetrahedrons be (x) 1 ，y 1 ，z 1 )，(x 2 ，y 2 ，z 2 )，(x 3 ，y 3 ，z 3 )，(x 4 ，y 4 ，z 4 ) For this regular tetrahedron, the preferred parameter is the locus designation U, whose coordinates can be expressed as:
Figure BDA0002966859690000051
for a regular tetrahedrons, the preferred parameter is the locus designation V, whose coordinates can be expressed as:
Figure BDA0002966859690000061
wherein k is x ，k y ，k z And the scale coefficient is represented, is the ratio of the calibration coefficients of the x, y and z dimensions corresponding to the optimal puncture sites of the alpha regular tetrahedrons to the calibration coefficient of the dimensionality of the tetrahedron, and has the value of 1.10.
Similarly, for a regular hexahedron with a length a, the preferred parameter is the locus number N, and the coordinates of N can be expressed as:
Figure BDA0002966859690000062
for the β regular hexahedrons, the preferred parameter is the site designation R, whose coordinates can be expressed as:
Figure BDA0002966859690000063
wherein l x ，l y ，l z The representative scaling factor is the ratio of the x, y and z dimension calibration factors corresponding to the beta regular hexahedron optimal puncture sites to the regular hexahedron dimension calibration factor with the edge length of a, and the value is 1.07.
Similarly, when there are γ octahedra, for each octahedra, the preferred parameter is the locus number Q, whose coordinates can be expressed as:
Figure BDA0002966859690000071
wherein j is x ，j y ，j z Representing a proportionality coefficient, and the value is 1.11.
Similarly, for a regular dodecahedron with a ridge length of b, the preferred parameter is the locus designation D, whose coordinates can be expressed as:
Figure BDA0002966859690000072
for η regular dodecahedrons, the preferred parameter is the site designation T, whose coordinates can be expressed as:
Figure BDA0002966859690000073
wherein s is x ，s y ，s z Representing a proportional coefficient, and calibrating coefficients for x, y and z dimensions corresponding to the eta n regular dodecahedron optimal puncture sitesThe ratio of the calibration coefficients of the dimensions of the regular dodecahedron with the edge length b is 1.19.
Similarly, for a regular icosahedron with a ridge length c, the preferred parameter is the site designation E, whose coordinates can be expressed as:
Figure BDA0002966859690000074
for the case of μ regular icosahedrons, the preferred parameter is the site designation Y, whose coordinates can be expressed as:
Figure BDA0002966859690000075
wherein, t x ，t y ，t z And the scaling coefficient is represented, and is the ratio of the calibration coefficients of the x, y and z dimensions corresponding to the mu optimal puncture sites of the regular icosahedron to the calibration coefficient of the dimension of the regular icosahedron with the edge length of c, and the value is 1.13.
S4, determining an optimal puncture liquid extraction point based on the partition body selection parameters:
and (4) calculating the optimal puncture site for treating the joint cavity effusion according to the optimal parameters of the five regular polyhedrons obtained in the step (S3).
Through calculation, the optimal puncture site coordinate is a matrix M,
Figure BDA0002966859690000081

Claims (6)

1. a method for determining an optimal point of effusion extraction by puncture for joint cavity treatment is characterized by comprising the following steps: the method comprises the following steps:
s1, acquiring image information of an original lesion part:
selecting a puncture part of a joint puncture point, carrying out axial plane and vertical section CT scanning on the part, reconstructing an obtained image into a three-dimensional space three-dimensional structure and determining the position of each pixel point in a coordinate system, wherein x is the horizontal direction of the axial plane and the vertical section, y is the vertical direction of the axial plane, z is the vertical direction of the vertical section, o is an origin point,
s2, decomposing mechanistic data based on the lesion site image:
the space polyhedron obtained in the S1 is divided along the crossed connecting line of the boundary characteristic points to obtain alpha regular tetrahedrons, beta regular hexahedrons, gamma regular octahedrons, delta regular dodecahedrons and epsilon regular icosahedrons,
s3, positioning segmentation body selection parameters based on pathological part data:
let four vertexes of one regular tetrahedron of the alpha regular tetrahedrons be (x) 1 ，y 1 ，z 1 )，(x 2 ，y 2 ，z 2 )，(x 3 ，y 3 ，z 3 )，(x 4 ，y 4 ，z 4 ) The vertex coordinates are obtained according to the image information of the original lesion site, for this regular tetrahedron, the optimal parameter is a site label U, and the coordinates of U can be expressed as:
Figure FDA0003808227240000011
for α regular tetrahedrons, the optimal parameter is the locus index V, whose coordinates can be expressed as:
Figure FDA0003808227240000012
Figure FDA0003808227240000021
wherein k is x ，k y ，k z Representing a proportionality coefficient which is the ratio of the x, y and z dimension calibration coefficients corresponding to the optimal puncture sites of alpha regular tetrahedrons to the calibration coefficient of the tetrahedron dimension,
similarly, for a regular hexahedron with a length of a, the optimal parameter is the locus index N, and the coordinates of N can be expressed as:
Figure FDA0003808227240000022
for the β regular hexahedrons, the optimal parameter is the site designation R, whose coordinates can be expressed as:
Figure FDA0003808227240000023
wherein l x ，l y ，l z Representing a proportionality coefficient which is the ratio of x, y and z dimension calibration coefficients corresponding to beta regular hexahedron optimal puncture sites to a regular hexahedron dimension calibration coefficient with the edge length a obtained from image information of an original lesion part,
similarly, when there are γ octahedrons, the optimal parameter for each octahedron is the locus number Q, whose coordinates can be expressed as:
Figure FDA0003808227240000024
wherein j is x ，j y ，j z Represents the coefficient of proportionality in the form of,
similarly, for a regular dodecahedron with a ridge length b, the optimal parameter is the locus index D, and the coordinates of D can be expressed as:
Figure FDA0003808227240000031
for δ regular dodecahedrons, the optimal parameter is the locus designation T, whose coordinates can be expressed as:
Figure FDA0003808227240000032
wherein s is x ，s y ，s z Representative ratio systemThe optimal puncture sites of the delta regular dodecahedron are corresponding to the ratio of the calibration coefficients of the three dimensions x, y and z to the calibration coefficient of the dimension of the regular dodecahedron with the edge length b, the edge length b is obtained by the image information of the original lesion part,
similarly, for a regular icosahedron with a ridge length of c, the optimal parameter is the locus index E, and the coordinates of E can be expressed as:
Figure FDA0003808227240000033
for epsilon regular icosahedrons, the optimal parameter is the locus designation Y, whose coordinates can be expressed as:
Figure FDA0003808227240000034
wherein, t x ，t y ，t z Representing a proportion coefficient which is the ratio of x, y and z dimension calibration coefficients corresponding to epsilon optimal puncture sites of the regular icosahedron to a dimension calibration coefficient of the regular icosahedron with the edge length c obtained by the image information of the original lesion part,
s4, determining an optimal puncture liquid-extracting point based on the segmentation body selection parameters:
calculating the optimal puncture site for treating the joint cavity effusion according to the optimal parameters of the five regular polyhedrons obtained in the step S3, wherein the coordinate of the optimal puncture site is a matrix M through calculation,
Figure FDA0003808227240000041
2. the method for determining the optimal point for puncturing and fluid collection for the treatment of the joint cavity according to claim 1, wherein the method comprises the following steps: k is a radical of x ，k y ，k z The values are 1.10 respectively.
3. The method of claim 1The method for determining the optimal point of the puncture effusion extraction for the joint cavity treatment is characterized by comprising the following steps: l x ，l y ，l z The values are 1.07 respectively.
4. The method for determining the optimal point for puncturing and effusion collection for joint cavity treatment according to claim 1, wherein: s x ，s y ，s z The values are 1.19 respectively.
5. The method for determining the optimal point for puncturing and effusion collection for joint cavity treatment according to claim 1, wherein: t is t x ，t y ，t z The values are 1.13 respectively.
6. The method for determining the optimal point for puncturing and effusion collection for joint cavity treatment according to claim 1, wherein: j is a unit of a group x ，j y ，j z The values are 1.11 respectively.
CN202110257430.1A 2021-03-09 2021-03-09 Optimal point calibration method for puncturing effusion drainage in joint cavity treatment Active CN113081258B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110257430.1A CN113081258B (en) 2021-03-09 2021-03-09 Optimal point calibration method for puncturing effusion drainage in joint cavity treatment

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110257430.1A CN113081258B (en) 2021-03-09 2021-03-09 Optimal point calibration method for puncturing effusion drainage in joint cavity treatment

Publications (2)

Publication Number Publication Date
CN113081258A CN113081258A (en) 2021-07-09
CN113081258B true CN113081258B (en) 2022-12-30

Family

ID=76666676

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110257430.1A Active CN113081258B (en) 2021-03-09 2021-03-09 Optimal point calibration method for puncturing effusion drainage in joint cavity treatment

Country Status (1)

Country Link
CN (1) CN113081258B (en)

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2012147652A1 (en) * 2011-04-27 2012-11-01 株式会社デージーエス・コンピュータ Puncture treatment support method, puncture treatment support device, and program for puncture treatment support device
CN110353774A (en) * 2018-12-15 2019-10-22 深圳铭杰医疗科技有限公司 Assist Needle-driven Robot and its control method, computer equipment, storage medium
CN110464459A (en) * 2019-07-10 2019-11-19 丽水市中心医院 Intervention plan navigation system and its air navigation aid based on CT-MRI fusion
CN111067597A (en) * 2019-12-10 2020-04-28 山东大学 System and method for determining puncture path according to human body posture in tumor puncture

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107296645B (en) * 2017-08-03 2020-04-14 东北大学 Optimal path planning method for lung puncture operation and lung puncture operation navigation system
WO2019129606A1 (en) * 2017-12-28 2019-07-04 Koninklijke Philips N.V. Apparatus and method for assisting puncture planning

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2012147652A1 (en) * 2011-04-27 2012-11-01 株式会社デージーエス・コンピュータ Puncture treatment support method, puncture treatment support device, and program for puncture treatment support device
CN110353774A (en) * 2018-12-15 2019-10-22 深圳铭杰医疗科技有限公司 Assist Needle-driven Robot and its control method, computer equipment, storage medium
CN110464459A (en) * 2019-07-10 2019-11-19 丽水市中心医院 Intervention plan navigation system and its air navigation aid based on CT-MRI fusion
CN111067597A (en) * 2019-12-10 2020-04-28 山东大学 System and method for determining puncture path according to human body posture in tumor puncture

Also Published As

Publication number Publication date
CN113081258A (en) 2021-07-09

Similar Documents

Publication Publication Date Title
US8199985B2 (en) Automatic interpretation of 3-D medicine images of the brain and methods for producing intermediate results
CN108573502B (en) Method for automatically measuring Cobb angle
WO2019000455A1 (en) Method and system for segmenting image
CN107665497B (en) Method for calculating cardiothoracic ratio in medical image
US8970581B2 (en) System and method for interactive contouring for 3D medical images
Pulagam et al. Automated lung segmentation from HRCT scans with diffuse parenchymal lung diseases
JP2005517249A (en) Correspondence identification of objects that do not use the entire volume
CN108830852B (en) Three-dimensional ultrasonic tumor auxiliary measurement system and method
CN1776357A (en) In-situ measuring method of plant root system geometry
CN111340825A (en) Method and system for generating mediastinal lymph node segmentation model
EP2525323A1 (en) Visualizing deformation fields for image transformations
CN113223004A (en) Liver image segmentation method based on deep learning
CN108510580A (en) A kind of vertebra CT image three-dimensional visualization methods
CN113081258B (en) Optimal point calibration method for puncturing effusion drainage in joint cavity treatment
CN113298826B (en) Image segmentation method based on LA-Net network
KR20060110456A (en) Method for analyzing hippocampus for aiding diagnosis of brain diseases and the recording media therein readable by computer
CN117237435B (en) Tumor prognosis effect evaluation method, device, electronic equipment and storage medium
CN112348826B (en) Interactive liver segmentation method based on geodesic distance and V-net
CN112258536B (en) Integrated positioning and segmentation method for calluses and cerebellum earthworm parts
CN110533667B (en) Lung tumor CT image 3D segmentation method based on image pyramid fusion
CN110033454B (en) Segmentation method of lung tumor with large-area adhesion lung boundary tissue in CT image
CN105869169A (en) Automatic dividing method of tumor issue micro array image
CN113012127A (en) Cardiothoracic ratio measuring method based on chest medical image
Ge et al. Accurate localization of cortical convolutions in MR brain images
Erdt et al. Computer aided segmentation of kidneys using locally shape constrained deformable models on CT images

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant