CN109662716B - Cartilage thickness measuring method, cartilage thickness measuring device, computer equipment and storage medium - Google Patents
Cartilage thickness measuring method, cartilage thickness measuring device, computer equipment and storage medium Download PDFInfo
- Publication number
- CN109662716B CN109662716B CN201811566467.7A CN201811566467A CN109662716B CN 109662716 B CN109662716 B CN 109662716B CN 201811566467 A CN201811566467 A CN 201811566467A CN 109662716 B CN109662716 B CN 109662716B
- Authority
- CN
- China
- Prior art keywords
- cartilage
- point
- interest point
- edge
- position information
- 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
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/103—Detecting, measuring or recording devices for testing the shape, pattern, colour, size or movement of the body or parts thereof, for diagnostic purposes
- A61B5/107—Measuring physical dimensions, e.g. size of the entire body or parts thereof
- A61B5/1075—Measuring physical dimensions, e.g. size of the entire body or parts thereof for measuring dimensions by non-invasive methods, e.g. for determining thickness of tissue layer
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/103—Detecting, measuring or recording devices for testing the shape, pattern, colour, size or movement of the body or parts thereof, for diagnostic purposes
- A61B5/107—Measuring physical dimensions, e.g. size of the entire body or parts thereof
- A61B5/1079—Measuring physical dimensions, e.g. size of the entire body or parts thereof using optical or photographic means
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/13—Edge detection
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/70—Determining position or orientation of objects or cameras
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
- G06T2207/30008—Bone
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Medical Informatics (AREA)
- Animal Behavior & Ethology (AREA)
- Pathology (AREA)
- Oral & Maxillofacial Surgery (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Dentistry (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Biophysics (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Computer Vision & Pattern Recognition (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Image Analysis (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
The application relates to a cartilage thickness measuring method, a cartilage thickness measuring device, computer equipment and a storage medium. The method comprises the following steps: acquiring cartilage surface contour information according to the 3D cartilage image data; identifying a first surface and a second surface of the cartilage surface from the cartilage surface contour information; calculating cartilage thickness of a first interest point in the first surface according to the position information of the first interest point; and/or calculating cartilage thickness of a second point of interest position in the second surface according to position information of the second point of interest. The method can be used for measuring the thickness of any point on the surface of the cartilage.
Description
Technical Field
The present application relates to the field of medical imaging technologies, and in particular, to a cartilage thickness measurement method and apparatus, a computer device, and a storage medium.
Background
Osteoarthritis is the most common rheumatic disease, and is increasingly valued by the fact that it imposes a heavy mental and economic burden on patients with a very high disability rate. By measuring the thickness of the articular cartilage and observing the change in the thickness, diagnosis of osteoarthritis, monitoring of the condition of an illness, evaluation of the effect of treatment, and the like can be assisted.
At present, the cartilage measurement method generally aims at several angles or layers of cartilage, but cannot measure the cartilage thickness at any point on the cartilage surface.
Disclosure of Invention
In view of the above, there is a need to provide a cartilage thickness method, apparatus, computer device and storage medium capable of measuring cartilage thickness at all points on the cartilage surface.
A cartilage thickness measurement method, the method comprising:
acquiring cartilage surface contour information according to the 3D cartilage image data;
identifying a first surface and a second surface of the cartilage surface from the cartilage surface contour information;
calculating cartilage thickness of a first interest point in the first surface according to the position information of the first interest point; and/or
And calculating the cartilage thickness of the second interest point position according to the position information of the second interest point in the second surface.
In one embodiment, the acquiring of the cartilage surface contour information according to the 3D cartilage image data includes: obtaining 3D cartilage model information according to the 3D cartilage image data, wherein the 3D cartilage model information comprises position information of a target point in the 3D cartilage model; judging whether a background point exists in the neighborhood of a target point according to the position information of the target point in the 3D cartilage model; and if the background points exist in the neighborhood of the target point, the target point is a point forming the contour of the cartilage surface, and the position information of the target point forms the position information of the midpoint of the contour of the cartilage surface.
In one embodiment, the calculating the cartilage thickness at the first point of interest position according to the position information of the first point of interest in the first surface includes:
and calculating the shortest distance from the first interest point to the second surface according to the position information of the first interest point in the first surface to obtain the cartilage thickness of the first interest point.
The calculating the cartilage thickness of the first interest point position according to the position information of the first interest point in the first surface comprises the following steps:
acquiring a first tangent plane of a first interest point on the cartilage surface according to the position information of the first interest point in the first surface;
obtaining a first intersection point of the first normal line and the second surface according to the first normal line of the first tangent plane on the first surface;
and calculating the distance from the first interest point to the first intersection point to obtain the cartilage thickness of the first interest point.
A cartilage thickness measurement method, the method comprising:
acquiring 2D cartilage section image data of a cartilage image;
acquiring cartilage edge information according to the 2D cartilage section image data, wherein the cartilage edge information carries position information of any point forming the edge;
identifying a first edge and a second edge of the cartilage edges according to the cartilage edge information;
calculating cartilage thickness of a first interest point in the first edge according to the position information of the first interest point; and/or
And calculating the cartilage thickness of the second interest point position according to the position information of the second interest point in the second edge.
In one embodiment, the calculating the cartilage thickness of the first interest point position according to the position information of the first interest point in the first edge includes: and calculating the shortest distance from the first interest point to the second edge according to the position information of the first interest point in the first edge to obtain the cartilage thickness of the first interest point.
In one embodiment, the calculating the cartilage thickness of the first interest point position according to the position information of the first interest point in the first edge includes:
acquiring a first tangent of a first interest point in the cartilage edge according to the position information of the first interest point in the first edge;
obtaining a first intersection point of the first normal line and the second edge according to a first normal line of the first tangent line at the first edge;
and calculating the distance from the first interest point to the first intersection point to obtain the cartilage thickness of the first interest point.
A cartilage thickness measurement device, the device comprising:
the surface contour acquisition module is used for acquiring cartilage surface contour information according to the 3D cartilage image data;
a surface classification module for identifying a first surface and a second surface in the cartilage surface according to the cartilage surface contour information;
the cartilage thickness calculation module is used for calculating the cartilage thickness of a first interest point in the first surface according to the position information of the first interest point; and/or calculating cartilage thickness of a second point of interest position in the second surface according to position information of the second point of interest.
A computer device comprising a memory and a processor, the memory storing a computer program, the processor implementing the following steps when executing the computer program:
acquiring cartilage surface contour information according to the 3D cartilage image data;
identifying a first surface and a second surface of the cartilage surface from the cartilage surface contour information;
calculating cartilage thickness of a first interest point in the first surface according to the position information of the first interest point; and/or
And calculating the cartilage thickness of the second interest point position according to the position information of the second interest point in the second surface.
A computer-readable storage medium, on which a computer program is stored which, when executed by a processor, carries out the steps of:
acquiring cartilage surface contour information according to the 3D cartilage image data;
identifying a first surface and a second surface of the cartilage surface from the cartilage surface contour information;
calculating cartilage thickness of a first interest point in the first surface according to the position information of the first interest point; and/or
And calculating the cartilage thickness of the second interest point position according to the position information of the second interest point in the second surface.
According to the cartilage thickness measuring method, the cartilage thickness measuring device, the computer equipment and the storage medium, the cartilage surface is separated from the cartilage image, the cartilage thickness of the point position is calculated according to the position information of the point on one surface of the cartilage, the cartilage thickness of any point position on the cartilage surface can be calculated, and a basis is provided for subsequent diagnosis and treatment of the cartilage.
Drawings
FIG. 1 is a schematic flow chart of a cartilage thickness measurement method according to an embodiment;
FIG. 2 is a schematic flow chart illustrating a method for obtaining a cartilage surface contour according to one embodiment;
FIG. 3 is a schematic representation of the cartilage surface contour in one embodiment;
FIG. 4 is a schematic flow chart of a cartilage thickness measurement method in another embodiment;
FIG. 5 is a schematic illustration of cartilage thickness calculation from normal in one embodiment;
FIG. 6 is a block diagram showing the structure of a cartilage thickness measuring device according to an embodiment;
FIG. 7 is a diagram illustrating an internal structure of a computer device according to an embodiment.
Detailed Description
In order to make the objects, technical solutions and advantages of the present application more apparent, the present application is described in further detail below with reference to the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are merely illustrative of the present application and are not intended to limit the present application.
In one embodiment, as shown in fig. 1, there is provided a cartilage thickness measurement method comprising the steps of:
and step S110, acquiring cartilage surface contour information according to the 3D cartilage image data.
The 3D cartilage image data may be MR medical images or CT medical images, including but not limited to knee joint cartilage, hip joint cartilage, and finger joint cartilage. In the process of capturing the 3D cartilage image, images of other bones or tissues are also collected, and therefore, the cartilage region needs to be separated from the 3D cartilage image according to the density or shape characteristics of cartilage.
And step S120, identifying a first surface and a second surface in the cartilage surface according to the cartilage surface contour information.
Wherein the cartilage is generally flat, and the first surface and the second surface can be the upper and lower surfaces of the flat cartilage. Specifically, the first surface and the second surface are surfaces connected with other bones, as shown in fig. 3, the surface of the crescent-shaped area is a cartilage surface contour, and the first surface and the second surface are respectively located on two sides of the crescent-shaped area.
The identifying the first surface and the second surface of the cartilage surface comprises adopting an image segmentation method, a differentiation method, an extreme value method, machine learning and surface fitting.
Step S130, calculating a cartilage thickness of a first point of interest in the first surface according to the position information of the first point of interest.
And calculating the cartilage thickness of the position of the first interest point through a geometric calculation formula according to the position information of the first interest point.
Step S140, calculating cartilage thickness of a second point of interest in the second surface according to the position information of the second point of interest.
Similarly, the cartilage thickness at the position of the second interest point is calculated through a geometric calculation formula according to the position information of the second interest point and the cartilage surface contour information.
In this embodiment, the execution sequence of step S130 and step S140 is not sequential, or only step S130 or step S140 may be executed, or both steps may be executed.
In one embodiment, in step S130 and step S140, the calculating the cartilage thickness of the first interest point position according to the position information of the first interest point in the first surface includes: and calculating the shortest distance from the first interest point to the second surface according to the position information of the first interest point in the first surface to obtain the cartilage thickness of the first interest point. Similarly, calculating the cartilage thickness of the second interest point position according to the position information of the second interest point in the second surface includes: and calculating the shortest distance from a second interest point to the first surface according to the position information of the second interest point in the second surface to obtain the cartilage thickness of the second interest point.
The positions of the first interest point and the second interest point are known, meanwhile, the position information of any point in the first surface and the second surface is known, and the shortest distance between the point and the surface is directly calculated. In this embodiment, the shortest distance is the thickness of the cartilage.
And calculating the shortest distance from the first interest point to the second surface, wherein the calculation method comprises a traversal method, a Lagrange multiplier method, a tracking method and a tangent plane method, and the point in the second surface corresponding to the shortest distance is called a cartilage thickness endpoint corresponding to the first interest point. The shortest distance of the second point of interest to the first surface, and its corresponding cartilage thickness end point, may also be obtained. In one embodiment, in step S130, the calculating a cartilage thickness of a first point of interest position in the first surface according to the position information of the first point of interest includes: acquiring a first tangent plane of a first interest point on the cartilage surface according to the position information of the first interest point in the first surface; obtaining a first intersection point of the first normal line and the second surface according to the first normal line of the first tangent plane on the first surface; and calculating the distance from the first interest point to the first intersection point to obtain the cartilage thickness of the first interest point. And the distance from the first interest point to the first intersection point is the cartilage thickness of the first interest point.
In one embodiment, in step S140, calculating the cartilage thickness of the second point of interest position in the second surface according to the position information of the second point of interest, includes: acquiring a second tangent plane of a second interest point on the cartilage surface according to the position information of the second interest point on the second surface; obtaining a second intersection point of the second normal line and the first surface according to a second normal line of the second tangent plane on the second surface; and calculating the distance from the second interest point to the second intersection point to obtain the cartilage thickness of the second interest point. And the distance from the second interest point to the second intersection point is the cartilage thickness of the second interest point.
In the cartilage thickness measuring method, the cartilage surface is separated from the cartilage image, and the cartilage thickness of the point position is calculated according to the position information of the point on one surface of the cartilage, so that the cartilage thickness of any point position on the cartilage surface can be calculated, and a basis is provided for subsequent diagnosis and treatment of the cartilage.
In one embodiment, in step S110, as shown in fig. 2, the acquiring cartilage surface contour information according to the 3D cartilage image data includes the steps of:
step S111, obtaining 3D cartilage model information according to the 3D cartilage image data, wherein the 3D cartilage model information comprises position information of a target point in the 3D cartilage model.
Wherein, the 3D cartilage model information is image data only including cartilage parts. The 3D cartilage model can reduce the shape and position of the cartilage.
Step S112, determining whether a background point exists in the neighborhood of the target point according to the position information of the target point in the 3D cartilage model.
The 3D cartilage model is composed of a plurality of cartilage image points, the points located in the 3D cartilage model only contain cartilage image points in the neighborhood, and the points located on the surface of the 3D cartilage model have image points of other bones or tissues in the neighborhood. The background points are image points of other bones or tissues except cartilage.
Step S113, if there is a background point in the neighborhood of the target point, the target point is a point forming a contour of the cartilage surface, and the position information of the target point forms the position information of a midpoint of the contour of the cartilage surface.
In one embodiment, as shown in fig. 4, there is provided a cartilage thickness measuring method including the steps of:
step S210, 2D cartilage section image data of the cartilage image is acquired.
The 2D cartilage section image data only comprises two-dimensional information of a cartilage image, and only a cartilage plane image can be displayed. For example, a 2D cartilage section image is taken perpendicular to the thickness direction of cartilage to fully reveal the thickness of cartilage.
Step S220, cartilage edge information is obtained according to the 2D cartilage section image data, and the cartilage edge information carries position information of any point forming the edge.
Wherein, the cartilage includes but is not limited to knee joint cartilage, hip joint cartilage and finger joint cartilage. Since the density or shape characteristics of cartilage are different from the surrounding bone and tissue, cartilage edge information can be extracted.
Step S230, identifying a first edge and a second edge of the cartilage edges according to the cartilage edge information.
Specifically, the cartilage edge encloses a crescent shape, and the first edge and the second edge are located on two sides of the crescent shape.
Step S240, calculating a cartilage thickness of a first interest point in the first edge according to the position information of the first interest point.
The first interest point is located on the first edge, the cartilage edge information is known, and the cartilage thickness of the first interest point can be calculated through a geometric calculation formula according to the position information of the first interest point.
Step S250, calculating cartilage thickness of a second interest point in the second edge according to the position information of the second interest point.
Similarly, the cartilage thickness at the position of the second interest point is calculated through a geometric calculation formula according to the position information and the cartilage edge information of the second interest point.
In this embodiment, the execution sequence of step S240 and step S250 is not sequential, or only step S240 or step S250 may be executed, or both steps may be executed.
In one embodiment, in step S240 and step S250, the calculating the cartilage thickness of the first interest point position according to the position information of the first interest point in the first edge includes: and calculating the shortest distance from the first interest point to the second edge according to the position information of the first interest point in the first edge to obtain the cartilage thickness of the first interest point. Similarly, the calculating the cartilage thickness of the second interest point position according to the position information of the second interest point in the second edge includes: and calculating the shortest distance from the second interest point to the first edge according to the position information of the second interest point in the second edge to obtain the cartilage thickness of the second interest point.
The positions of the first interest point and the second interest point are known, meanwhile, the position information of any one point of the first edge and the second edge is known, and the shortest distance between the point and the surface is directly calculated. In this embodiment, the shortest distance is the thickness of the cartilage.
And calculating the shortest distance from the first interest point to the second edge, wherein the calculation method comprises a traversal method, a Lagrangian multiplier method and a tangent method, and the point in the second edge corresponding to the shortest distance is called a cartilage thickness endpoint corresponding to the first interest point. The shortest distance from the second point of interest to the first edge, and its corresponding cartilage thickness end point, may also be obtained.
In one embodiment, in step S240, calculating the cartilage thickness of the first interest point position according to the position information of the first interest point in the first edge includes: acquiring a first tangent of a first interest point in the cartilage edge according to the position information of the first interest point in the first edge; obtaining a first intersection point of the first normal line and the second edge according to a first normal line of the first tangent line at the first edge; and calculating the distance from the first interest point to the first intersection point to obtain the cartilage thickness of the first interest point. And the distance from the first interest point to the first intersection point is the cartilage thickness of the first interest point.
In one embodiment, in step S250, calculating the cartilage thickness of the second interest point position according to the position information of the second interest point in the second edge includes: acquiring a second tangent of a second interest point at the cartilage edge according to the position information of the second interest point in the second edge; obtaining a second intersection point of the second normal line and the first edge according to a second normal line of the second tangent line at the second edge; and calculating the distance from the second interest point to the second intersection point to obtain the cartilage thickness of the second interest point. And the distance from the second interest point to the second intersection point is the cartilage thickness of the second interest point.
Specifically, as shown in fig. 5, the crescent shape is a cartilage region, the first edge and the second edge are located on two sides of the crescent shape, a normal MN perpendicular to the second edge is made at one point B in the second edge, the normal MN intersects with the first edge at a point a, the distance from the point B to the point a is calculated, and the distance between the points AB is the cartilage thickness at the point B.
According to the cartilage thickness measuring method, the cartilage edge is divided into the first edge and the second edge, and the cartilage thickness of the point position is calculated according to the position information of the point in the first edge and the second edge, so that the cartilage thickness of any point position on the cartilage surface can be calculated, and a basis is provided for subsequent diagnosis and treatment of cartilage.
It should be understood that although the steps in the flowcharts of fig. 1, 2, and 4 are shown in order as indicated by the arrows, the steps are not necessarily performed in order as indicated by the arrows. The steps are not performed in the exact order shown and described, and may be performed in other orders, unless explicitly stated otherwise. Moreover, at least some of the steps in fig. 1, 2, and 4 may include multiple sub-steps or multiple stages, which are not necessarily performed at the same time, but may be performed at different times, and the order of performing the sub-steps or stages is not necessarily sequential, but may be performed alternately or alternatingly with other steps or at least some of the sub-steps or stages of other steps.
In one embodiment, as shown in fig. 6, there is provided a cartilage thickness measuring device including: a surface contour acquisition module 310, a surface classification module 320, and a cartilage thickness calculation module 330, wherein:
the surface contour obtaining module 310 is configured to obtain cartilage surface contour information according to the 3D cartilage image data.
A surface classification module 320, configured to identify a first surface and a second surface of the cartilage surface according to the cartilage surface contour information.
A cartilage thickness calculating module 330, configured to calculate a cartilage thickness at a first point of interest in the first surface according to the position information of the first point of interest; and/or calculating cartilage thickness of a second point of interest position in the second surface according to position information of the second point of interest.
In one embodiment, the cartilage thickness calculating module 330 is specifically configured to calculate a shortest distance from a first interest point in the first surface to the second surface according to the position information of the first interest point, so as to obtain the cartilage thickness at the first interest point.
In one embodiment, the cartilage thickness calculating module 330 is specifically configured to calculate a shortest distance from a second interest point in the second surface to the first surface according to the position information of the second interest point in the second surface, so as to obtain the cartilage thickness at the second interest point.
In one embodiment, the cartilage thickness calculation module 330 includes: the first tangent plane acquisition unit is used for acquiring a first tangent plane of a first interest point on the cartilage surface according to the position information of the first interest point in the first surface; a first intersection point obtaining unit, configured to obtain a first intersection point of the first normal line and the second surface according to a first normal line of the first tangent plane on the first surface; and the thickness calculating unit is used for calculating the distance from the first interest point to the first intersection point to obtain the cartilage thickness of the first interest point. And the distance from the first interest point to the first intersection point is the cartilage thickness of the first interest point.
In one embodiment, the cartilage thickness calculation module 330 includes: a second tangent plane obtaining unit, configured to obtain a second tangent plane of a second point of interest on the cartilage surface according to position information of the second point of interest on the second surface; a second intersection point obtaining unit, configured to obtain a second intersection point of the second normal line and the first surface according to a second normal line of the second tangent plane on the second surface; and the thickness calculating unit is used for calculating the distance from the second interest point to the second intersection point to obtain the cartilage thickness of the second interest point.
In one embodiment, the surface profile obtaining module 310 includes: the cartilage model obtaining unit is used for obtaining 3D cartilage model information according to the 3D cartilage image data, wherein the 3D cartilage model information comprises position information of a target point in the 3D cartilage model; and the background point judging unit is used for judging whether a background point exists in the neighborhood of the target point according to the position information of the target point in the 3D cartilage model, if the background point exists in the neighborhood of the target point, the target point is a point forming a cartilage surface contour, and the position information of the target point forms the position information of a midpoint of the cartilage surface contour.
In one embodiment, there is provided a cartilage thickness measurement device comprising:
and the section image acquisition module is used for acquiring 2D cartilage section image data of the cartilage image.
And the cartilage edge acquisition module is used for acquiring cartilage edge information according to the 2D cartilage section image data, wherein the cartilage edge information carries position information of any point forming the edge.
And the edge classification module is used for identifying a first edge and a second edge in the cartilage edges according to the cartilage edge information.
The cartilage thickness calculation module is used for calculating the cartilage thickness of a first interest point in the first edge according to the position information of the first interest point; and/or calculating the cartilage thickness of the second interest point position according to the position information of the second interest point in the second edge.
In one embodiment, the cartilage thickness calculation module is specifically configured to calculate a shortest distance from a first interest point in the first edge to the second edge according to position information of the first interest point, so as to obtain a cartilage thickness at the first interest point; and/or calculating the shortest distance from a second interest point to the first edge according to the position information of the second interest point in the second edge to obtain the cartilage thickness of the second interest point.
In one embodiment, the cartilage thickness calculation module comprises: the first tangent line obtaining unit is used for obtaining a first tangent line of a first interest point in the cartilage edge according to the position information of the first interest point in the first edge; a first intersection point obtaining unit, configured to obtain a first intersection point of the first normal line and the second edge according to a first normal line of the first tangent line at the first edge; and the distance calculation unit is used for calculating the distance from the first interest point to the first intersection point to obtain the cartilage thickness of the first interest point. And the distance from the first interest point to the first intersection point is the cartilage thickness of the first interest point.
In one embodiment, the cartilage thickness calculation module comprises: a second tangent line obtaining unit, configured to obtain a second tangent line of a second interest point in the cartilage edge according to position information of the second interest point in the second edge; a second intersection bottom obtaining unit, configured to obtain a second intersection point of the second normal line and the first edge according to a second normal line of the second tangent line at the second edge; and the distance calculation unit is used for calculating the distance from the second interest point to the second intersection point to obtain the cartilage thickness of the second interest point. And the distance from the second interest point to the second intersection point is the cartilage thickness of the second interest point.
For specific definition of the cartilage thickness measuring device, reference may be made to the above definition of the cartilage thickness measuring method, which is not described herein again. The modules in the cartilage thickness measuring device can be wholly or partially realized by software, hardware and a combination thereof. The modules can be embedded in a hardware form or independent from a processor in the computer device, and can also be stored in a memory in the computer device in a software form, so that the processor can call and execute operations corresponding to the modules.
In one embodiment, a computer device is provided, which may be a server, the internal structure of which may be as shown in fig. 7. The computer device includes a processor, a memory, a network interface, and a database connected by a system bus. Wherein the processor of the computer device is configured to provide computing and control capabilities. The memory of the computer device comprises a nonvolatile storage medium and an internal memory. The non-volatile storage medium stores an operating system, a computer program, and a database. The internal memory provides an environment for the operation of an operating system and computer programs in the non-volatile storage medium. The database of the computer device is used for storing cartilage image data. The network interface of the computer device is used for communicating with an external terminal through a network connection. The computer program is executed by a processor to implement a cartilage thickness measurement method.
Those skilled in the art will appreciate that the architecture shown in fig. 7 is merely a block diagram of some of the structures associated with the disclosed aspects and is not intended to limit the computing devices to which the disclosed aspects apply, as particular computing devices may include more or less components than those shown, or may combine certain components, or have a different arrangement of components.
In one embodiment, a computer device is provided, comprising a memory and a processor, the memory having a computer program stored therein, the processor implementing the following steps when executing the computer program:
acquiring cartilage surface contour information according to the 3D cartilage image data;
identifying a first surface and a second surface of the cartilage surface from the cartilage surface contour information;
calculating cartilage thickness of a first interest point in the first surface according to the position information of the first interest point; and/or
And calculating the cartilage thickness of the second interest point position according to the position information of the second interest point in the second surface.
In one embodiment, a computer-readable storage medium is provided, having a computer program stored thereon, which when executed by a processor, performs the steps of:
acquiring cartilage surface contour information according to the 3D cartilage image data;
identifying a first surface and a second surface of the cartilage surface from the cartilage surface contour information;
calculating cartilage thickness of a first interest point in the first surface according to the position information of the first interest point; and/or
And calculating the cartilage thickness of the second interest point position according to the position information of the second interest point in the second surface.
It will be understood by those skilled in the art that all or part of the processes of the methods of the embodiments described above can be implemented by hardware instructions of a computer program, which can be stored in a non-volatile computer-readable storage medium, and when executed, can include the processes of the embodiments of the methods described above. Any reference to memory, storage, database, or other medium used in the embodiments provided herein may include non-volatile and/or volatile memory, among others. Non-volatile memory can include read-only memory (ROM), Programmable ROM (PROM), Electrically Programmable ROM (EPROM), Electrically Erasable Programmable ROM (EEPROM), or flash memory. Volatile memory can include Random Access Memory (RAM) or external cache memory. By way of illustration and not limitation, RAM is available in a variety of forms such as Static RAM (SRAM), Dynamic RAM (DRAM), Synchronous DRAM (SDRAM), Double Data Rate SDRAM (DDRSDRAM), Enhanced SDRAM (ESDRAM), Synchronous Link DRAM (SLDRAM), Rambus Direct RAM (RDRAM), direct bus dynamic RAM (DRDRAM), and memory bus dynamic RAM (RDRAM).
The technical features of the above embodiments can be arbitrarily combined, and for the sake of brevity, all possible combinations of the technical features in the above embodiments are not described, but should be considered as the scope of the present specification as long as there is no contradiction between the combinations of the technical features.
The above-mentioned embodiments only express several embodiments of the present application, and the description thereof is more specific and detailed, but not construed as limiting the scope of the invention. It should be noted that, for a person skilled in the art, several variations and modifications can be made without departing from the concept of the present application, which falls within the scope of protection of the present application. Therefore, the protection scope of the present patent shall be subject to the appended claims.
Claims (10)
1. A cartilage thickness measurement method, the method comprising:
acquiring cartilage surface contour information according to the 3D cartilage image data;
identifying a first surface and a second surface in the cartilage surface according to the cartilage surface contour information, wherein the first surface and the second surface are the upper surface and the lower surface of the cartilage respectively, and particularly are surfaces connected with other bones;
calculating cartilage thickness of a first interest point in the first surface according to position information of the first interest point, wherein the first interest point is any point in the first surface; and/or
Calculating cartilage thickness of a second interest point in the second surface according to position information of the second interest point in the second surface, wherein the second interest point is any point in the second surface;
the calculating the cartilage thickness of the first interest point position according to the position information of the first interest point in the first surface comprises the following steps:
and calculating the shortest distance from the first interest point to the second surface according to the position information of the first interest point in the first surface to obtain the cartilage thickness of the first interest point.
2. A cartilage thickness measurement method, the method comprising:
acquiring cartilage surface contour information according to the 3D cartilage image data;
identifying a first surface and a second surface in the cartilage surface according to the cartilage surface contour information, wherein the first surface and the second surface are the upper surface and the lower surface of the cartilage respectively, and particularly are surfaces connected with other bones;
calculating cartilage thickness of a first interest point in the first surface according to position information of the first interest point, wherein the first interest point is any point in the first surface; and/or
Calculating cartilage thickness of a second interest point in the second surface according to position information of the second interest point in the second surface, wherein the second interest point is any point in the second surface;
the calculating the cartilage thickness of the first interest point position according to the position information of the first interest point in the first surface comprises the following steps:
acquiring a first tangent plane of a first interest point on the cartilage surface according to the position information of the first interest point in the first surface;
obtaining a first intersection point of the first normal line and the second surface according to the first normal line of the first tangent plane on the first surface;
and calculating the distance from the first interest point to the first intersection point to obtain the cartilage thickness of the first interest point.
3. The method according to claim 1 or 2, wherein the obtaining of the cartilage surface contour information from the 3D cartilage image data comprises:
obtaining 3D cartilage model information according to the 3D cartilage image data, wherein the 3D cartilage model information comprises position information of a target point in the 3D cartilage model;
judging whether a background point exists in the neighborhood of a target point according to the position information of the target point in the 3D cartilage model;
and if the background points exist in the neighborhood of the target point, the target point is a point forming the contour of the cartilage surface, and the position information of the target point forms the position information of the midpoint of the contour of the cartilage surface.
4. The method of claim 1, wherein calculating the shortest distance of the first point of interest to the second surface comprises: at least one of a traversal method, a Lagrange multiplier method, a tracking method, and a tangent plane method.
5. A cartilage thickness measurement method, the method comprising:
acquiring 2D cartilage section image data of a cartilage image;
acquiring cartilage edge information according to the 2D cartilage section image data, wherein the cartilage edge information carries position information of any point forming the edge;
identifying a first edge and a second edge of the cartilage edges according to the cartilage edge information;
calculating cartilage thickness of a first interest point in the first edge according to the position information of the first interest point; and/or
Calculating cartilage thickness of a second interest point in the second edge according to the position information of the second interest point;
the calculating the cartilage thickness of the first interest point position according to the position information of the first interest point in the first edge comprises the following steps:
and calculating the shortest distance from the first interest point to the second edge according to the position information of the first interest point in the first edge to obtain the cartilage thickness of the first interest point.
6. A cartilage thickness measurement method, the method comprising:
acquiring 2D cartilage section image data of a cartilage image;
acquiring cartilage edge information according to the 2D cartilage section image data, wherein the cartilage edge information carries position information of any point forming the edge;
identifying a first edge and a second edge of the cartilage edges according to the cartilage edge information;
calculating cartilage thickness of a first interest point in the first edge according to the position information of the first interest point; and/or
Calculating cartilage thickness of a second interest point in the second edge according to the position information of the second interest point;
the calculating the cartilage thickness of the first interest point position according to the position information of the first interest point in the first edge comprises the following steps:
acquiring a first tangent of a first interest point in the cartilage edge according to the position information of the first interest point in the first edge;
obtaining a first intersection point of the first normal line and the second edge according to a first normal line of the first tangent line at the first edge;
and calculating the distance from the first interest point to the first intersection point to obtain the cartilage thickness of the first interest point.
7. The method of claim 5, wherein calculating the shortest distance from the first point of interest to the second edge comprises: at least one of a traversal method, a Lagrangian multiplier method, and a tangent method.
8. A cartilage thickness measurement device, the device comprising:
the surface contour acquisition module is used for acquiring cartilage surface contour information according to the 3D cartilage image data;
the surface classification module is used for identifying a first surface and a second surface in the cartilage surface according to the cartilage surface contour information, wherein the first surface and the second surface are the upper surface and the lower surface of the cartilage respectively, and particularly are surfaces connected with other bones;
the cartilage thickness calculating module is used for calculating the cartilage thickness of a first interest point in the first surface according to the position information of the first interest point, wherein the first interest point is any point in the first surface; and/or calculating the cartilage thickness of a second interest point in the second surface according to the position information of the second interest point in the second surface, wherein the second interest point is any point in the second surface; the calculating the cartilage thickness of the first interest point position according to the position information of the first interest point in the first surface comprises the following steps: and calculating the shortest distance from the first interest point to the second surface according to the position information of the first interest point in the first surface to obtain the cartilage thickness of the first interest point.
9. A computer device comprising a memory and a processor, the memory storing a computer program, wherein the processor implements the steps of the method of any one of claims 1 to 7 when executing the computer program.
10. A computer-readable storage medium, on which a computer program is stored, which, when being executed by a processor, carries out the steps of the method of any one of claims 1 to 7.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811566467.7A CN109662716B (en) | 2018-12-19 | 2018-12-19 | Cartilage thickness measuring method, cartilage thickness measuring device, computer equipment and storage medium |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811566467.7A CN109662716B (en) | 2018-12-19 | 2018-12-19 | Cartilage thickness measuring method, cartilage thickness measuring device, computer equipment and storage medium |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109662716A CN109662716A (en) | 2019-04-23 |
| CN109662716B true CN109662716B (en) | 2021-10-22 |
Family
ID=66144125
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811566467.7A Active CN109662716B (en) | 2018-12-19 | 2018-12-19 | Cartilage thickness measuring method, cartilage thickness measuring device, computer equipment and storage medium |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109662716B (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110827285B (en) * | 2019-11-15 | 2022-07-26 | 上海联影智能医疗科技有限公司 | Cartilage thickness detection method and device, computer equipment and readable storage medium |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1735373A (en) * | 2003-01-07 | 2006-02-15 | 成像治疗仪股份有限公司 | Methods of predicting musculoskeletal disease |
| CN1780594A (en) * | 2002-11-07 | 2006-05-31 | 康复米斯公司 | Methods for determining meniscal size and shape and for devising treatment |
| CN102438559A (en) * | 2009-02-25 | 2012-05-02 | 穆罕默德·拉什万·马赫福兹 | Customized orthopaedic implants and related methods |
| CN105361883A (en) * | 2014-08-22 | 2016-03-02 | 方学伟 | Method for determining lower limb biological force line in three-dimensional space for total knee arthroplasty |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| AU9088701A (en) * | 2000-09-14 | 2002-03-26 | Univ Leland Stanford Junior | Assessing condition of a joint and cartilage loss |
-
2018
- 2018-12-19 CN CN201811566467.7A patent/CN109662716B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1780594A (en) * | 2002-11-07 | 2006-05-31 | 康复米斯公司 | Methods for determining meniscal size and shape and for devising treatment |
| CN1735373A (en) * | 2003-01-07 | 2006-02-15 | 成像治疗仪股份有限公司 | Methods of predicting musculoskeletal disease |
| CN102438559A (en) * | 2009-02-25 | 2012-05-02 | 穆罕默德·拉什万·马赫福兹 | Customized orthopaedic implants and related methods |
| CN105361883A (en) * | 2014-08-22 | 2016-03-02 | 方学伟 | Method for determining lower limb biological force line in three-dimensional space for total knee arthroplasty |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109662716A (en) | 2019-04-23 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110766730B (en) | Image registration and follow-up evaluation method, storage medium and computer equipment | |
| CN111080573B (en) | Rib image detection method, computer device and storage medium | |
| CN110415792B (en) | Image detection method, image detection device, computer equipment and storage medium | |
| EP3724849B1 (en) | Image analysis for scoring motion of a heart wall | |
| CN110717905B (en) | Brain image detection method, computer device, and storage medium | |
| CN110599465B (en) | Image positioning method and device, computer equipment and storage medium | |
| CN110838104B (en) | Multi-time point region of interest matching method, device and storage medium | |
| JP6442530B2 (en) | Shape prediction | |
| CN110717961B (en) | Multi-modal image reconstruction method and device, computer equipment and storage medium | |
| CN110210519B (en) | Classification method, computer device, and storage medium | |
| CN110706207A (en) | Image quantization method, image quantization device, computer equipment and storage medium | |
| US20220036575A1 (en) | Method for measuring volume of organ by using artificial neural network, and apparatus therefor | |
| CN111340756A (en) | Medical image lesion detection and combination method, system, terminal and storage medium | |
| CN108765421B (en) | Breast medical image processing method and device and AEC exposure parameter acquisition method | |
| CN109662716B (en) | Cartilage thickness measuring method, cartilage thickness measuring device, computer equipment and storage medium | |
| CN110310257B (en) | Medical image processing method, apparatus, computer device and storage medium | |
| CN115063397A (en) | Computer-aided image analysis method, computer device and storage medium | |
| Queirós et al. | Assessment of aortic valve tract dynamics using automatic tracking of 3D transesophageal echocardiographic images | |
| CN116681716B (en) | Method, device, equipment and storage medium for dividing intracranial vascular region of interest | |
| CN111373438A (en) | Method and apparatus for imaging an organ | |
| CN111681205B (en) | Image analysis method, computer device, and storage medium | |
| CN111161369B (en) | Image reconstruction storage method, device, computer equipment and storage medium | |
| CN111160442B (en) | Image classification method, computer device, and storage medium | |
| KR101659056B1 (en) | Automated diagnosis system for craniosynostosis using a 2d shape descriptor and automated diagnosis method for craniosynostosis using the same | |
| CN113348485A (en) | Abnormality detection method, abnormality detection program, abnormality detection device, server device, and information processing method |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| CB02 | Change of applicant information |
Address after: 201807 Shanghai City, north of the city of Jiading District Road No. 2258 Applicant after: Shanghai Lianying Medical Technology Co., Ltd Address before: 201807 Shanghai City, north of the city of Jiading District Road No. 2258 Applicant before: SHANGHAI UNITED IMAGING HEALTHCARE Co.,Ltd. |
|
| CB02 | Change of applicant information | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |