CN113040812A - Medical image analysis method, medical image analysis apparatus, computer device, and storage medium - Google Patents
Medical image analysis method, medical image analysis apparatus, computer device, and storage medium Download PDFInfo
- Publication number
- CN113040812A CN113040812A CN202110148160.0A CN202110148160A CN113040812A CN 113040812 A CN113040812 A CN 113040812A CN 202110148160 A CN202110148160 A CN 202110148160A CN 113040812 A CN113040812 A CN 113040812A
- Authority
- CN
- China
- Prior art keywords
- interest
- region
- image
- blood vessel
- parameter 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.)
- Granted
Links
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/48—Diagnostic techniques
- A61B8/481—Diagnostic techniques involving the use of contrast agent, e.g. microbubbles introduced into the bloodstream
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/0059—Measuring for diagnostic purposes; Identification of persons using light, e.g. diagnosis by transillumination, diascopy, fluorescence
- A61B5/0062—Arrangements for scanning
- A61B5/0066—Optical coherence imaging
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/08—Detecting organic movements or changes, e.g. tumours, cysts, swellings
- A61B8/0891—Detecting organic movements or changes, e.g. tumours, cysts, swellings for diagnosis of blood vessels
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/52—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/5215—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/52—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/5215—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data
- A61B8/5238—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data for combining image data of patient, e.g. merging several images from different acquisition modes into one image
- A61B8/5261—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data for combining image data of patient, e.g. merging several images from different acquisition modes into one image combining images from different diagnostic modalities, e.g. ultrasound and X-ray
-
- 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
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Heart & Thoracic Surgery (AREA)
- Molecular Biology (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Pathology (AREA)
- Radiology & Medical Imaging (AREA)
- Physics & Mathematics (AREA)
- Biomedical Technology (AREA)
- Veterinary Medicine (AREA)
- Medical Informatics (AREA)
- Biophysics (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Hematology (AREA)
- Vascular Medicine (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
The application relates to a medical image analysis method, a medical image analysis device, a computer device and a storage medium. The method comprises the following steps: acquiring an angiographic image to be analyzed; obtaining a first region of interest in the angiographic image according to the angiographic image; analyzing the first region of interest to obtain first parameter information corresponding to the first region of interest; obtaining an analysis result according to the first parameter information; the analysis result is used for indicating whether the first region of interest meets a preset first condition. By adopting the method, the analysis result for indicating whether the first region of interest in the angiographic image to be analyzed meets the preset first condition can be accurately obtained, and the accuracy of the obtained analysis result of the angiographic image to be analyzed is improved.
Description
Technical Field
The present application relates to the field of medical image technology, and in particular, to a medical image analysis method, apparatus, computer device, and storage medium.
Background
The bioabsorbable stent has the advantages of avoiding coronary artery metallization caused by long-term existence of a metal stent and the like, and is widely applied to interventional therapy of coronary heart diseases. However, before the bioabsorbable stent is implanted into the human body, it is first determined whether the bioabsorbable stent can be implanted by combining the lesion site of the patient and the instruction for using the bioabsorbable stent.
In the traditional technology, a doctor manually marks blood vessels in an angiography image of a patient, a coronary angiography quantitative analysis method is used for obtaining blood vessel parameters of the patient, and whether the patient is suitable for implanting the bioabsorbable stent is judged according to the obtained blood vessel parameters.
However, the conventional determination method has a problem of low determination efficiency.
Disclosure of Invention
In view of the above, there is a need to provide a medical image analysis method, apparatus, computer device and storage medium capable of improving the efficiency of determining whether a patient is suitable for implanting a bioabsorbable stent.
A method of medical image analysis, the method comprising:
acquiring an angiographic image to be analyzed;
obtaining a first region of interest in the angiographic image according to the angiographic image;
analyzing the first region of interest to obtain first parameter information corresponding to the first region of interest;
obtaining an analysis result according to the first parameter information; the analysis result is used for indicating whether the first region of interest meets a preset first condition.
In one embodiment, the obtaining a first region of interest in the angiographic image from the angiographic image comprises:
and inputting the angiogram image into a preset positioning network, and obtaining a first region of interest in the angiogram image through the positioning network.
In one embodiment, the first parameter information includes: the position information of the first region of interest, the blood vessel distribution information in a preset range in the first region of interest, the blood vessel shape information in the preset range, and the plaque quantitative analysis information of the first region of interest.
In one embodiment, the method further comprises:
acquiring an intra-cavity image corresponding to the angiogram image;
determining a second region of interest corresponding to the first region of interest in the intracavity image;
analyzing the second region of interest to obtain second parameter information corresponding to the second region of interest;
and obtaining the analysis result according to the first parameter information and the second parameter information.
In one embodiment, the determining, in the intracavity image, a second region of interest corresponding to the first region of interest includes:
calculating a first blood vessel diameter distribution map corresponding to the angiogram image and a second blood vessel diameter distribution map corresponding to the intracavity image;
calculating a first blood vessel topological structure key point corresponding to the angiogram image and a second blood vessel topological structure key point corresponding to the intracavity image;
and determining a second region of interest corresponding to the first region of interest according to the first blood vessel diameter distribution map, the second blood vessel diameter distribution map, the first blood vessel topological structure key points and the second blood vessel topological structure key points.
In one embodiment, the second parameter information includes: the position information of the second region of interest, the blood vessel distribution information in a preset range in the second region of interest, the blood vessel shape information in the preset range, and the plaque quantitative analysis information of the second region of interest.
In one embodiment, the intraluminal images include intravascular ultrasound images and optical coherence tomography.
A medical image analysis apparatus, the apparatus comprising:
the first acquisition module is used for acquiring an angiographic image to be analyzed;
the second acquisition module is used for obtaining a first region of interest in the angiography image according to the angiography image;
the first analysis module is used for analyzing the first region of interest to obtain first parameter information corresponding to the first region of interest;
the third acquisition module is used for acquiring an analysis result according to the first parameter information; the analysis result is used for indicating whether the first region of interest meets a preset first condition.
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 an angiographic image to be analyzed;
obtaining a first region of interest in the angiographic image according to the angiographic image;
analyzing the first region of interest to obtain first parameter information corresponding to the first region of interest;
obtaining an analysis result according to the first parameter information; the analysis result is used for indicating whether the first region of interest meets a preset first condition.
A computer-readable storage medium, on which a computer program is stored which, when executed by a processor, carries out the steps of:
acquiring an angiographic image to be analyzed;
obtaining a first region of interest in the angiographic image according to the angiographic image;
analyzing the first region of interest to obtain first parameter information corresponding to the first region of interest;
obtaining an analysis result according to the first parameter information; the analysis result is used for indicating whether the first region of interest meets a preset first condition.
According to the medical image analysis method, the medical image analysis device, the computer equipment and the storage medium, the first parameter information corresponding to the first region of interest can be obtained by analyzing the first region of interest in the angiographic image to be analyzed, and then the angiographic image to be analyzed can be analyzed according to the first parameter information corresponding to the first region of interest to obtain the analysis result for indicating whether the first region of interest in the angiographic image to be analyzed meets the preset first condition, because the method is used for analyzing the first region of interest in the angiographic image to be analyzed, the angiographic image to be analyzed can be accurately analyzed according to the first parameter information corresponding to the first region of interest, so that the analysis result for indicating whether the first region of interest in the angiographic image to be analyzed meets the preset first condition can be accurately obtained, the accuracy of the analysis result of the obtained angiographic image to be analyzed is improved.
Drawings
FIG. 1 is a diagram of an exemplary medical image analysis method;
FIG. 2 is a schematic flow chart diagram of a method of medical image analysis in one embodiment;
FIG. 3 is a schematic diagram of a U-Net network in one embodiment;
FIG. 4 is a diagram illustrating a method for calculating vascularity information within a predetermined range and vascularity information within a predetermined range in a first region of interest according to an embodiment;
FIG. 5 is a schematic diagram of a method for analyzing blobs from a first region of interest in one embodiment;
FIG. 6 is a flow chart illustrating a method of medical image analysis in accordance with another embodiment;
FIG. 7 is a flow diagram illustrating a method of medical image analysis in one embodiment;
FIG. 7a is a schematic flow chart diagram illustrating a method for medical image analysis in accordance with another embodiment;
FIG. 8 is a diagram of a report format in one embodiment;
fig. 9 is a block diagram showing the configuration of a medical image analysis apparatus 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.
The medical image analysis method provided by the embodiment of the application can be applied to computer equipment shown in fig. 1. The computer device comprises a processor and a memory connected by a system bus, wherein a computer program is stored in the memory, and the steps of the method embodiments described below can be executed when the processor executes the computer program. Optionally, the computer device may further comprise a network interface, a display screen and an input device. Wherein the processor of the computer device is configured to provide computing and control capabilities. The memory of the computer device includes a nonvolatile storage medium storing an operating system and a computer program, and an internal memory. The internal memory provides an environment for the operation of an operating system and computer programs in the non-volatile storage medium. The network interface of the computer device is used for communicating with an external terminal through a network connection. Optionally, the computer device may be a server, a personal computer, a personal digital assistant, other terminal devices such as a tablet computer, a mobile phone, and the like, or a cloud or a remote server, and the specific form of the computer device is not limited in the embodiment of the present application.
In one embodiment, as shown in fig. 2, a medical image analysis method is provided, which is illustrated by applying the method to the computer device in fig. 1, and comprises the following steps:
s201, an angiographic image to be analyzed is acquired.
The angiographic image to be analyzed may be a coronary angiographic image or an angiographic image of another part. In particular, a computer device acquires an angiographic image to be analyzed. Alternatively, the computer device may acquire an angiographic image to be analyzed from a PACS (Picture Archiving and Communication Systems) server, or may acquire the angiographic image to be analyzed in real time from a corresponding medical imaging device. Optionally, after the computer device acquires the angiographic image to be analyzed, the angiographic image to be analyzed may be preprocessed to remove interference in the angiographic image to be analyzed, where the preprocessing includes resampling processing, size adjustment, gray scale normalization processing, and the like.
S202, obtaining a first region of interest in the angiography image according to the angiography image.
Specifically, the computer device obtains a first region of interest in the angiographic image according to the obtained angiographic image. Optionally, the computer device may obtain the first region of interest in the angiographic image to be analyzed by using a feature detection method, or may input the angiographic image to be analyzed into a preset neural network to obtain the first region of interest in the angiographic image to be analyzed.
S203, analyzing the first region of interest to obtain first parameter information corresponding to the first region of interest.
Specifically, the computer device analyzes a first region of interest in the angiographic image to obtain first parameter information corresponding to the first region of interest. Optionally, the computer device may analyze the position of the first region of interest, analyze the blood vessel distribution information in a preset range in the first region of interest, analyze the blood vessel morphology information in the preset range, and analyze the plaque amount in the first region of interest.
S204, obtaining an analysis result according to the first parameter information; the analysis result is used for indicating whether the first region of interest meets a preset first condition.
Specifically, the computer device obtains an analysis result for indicating whether the first region of interest meets a preset first condition according to the obtained first parameter information. Optionally, the preset first condition may be whether the first region of interest may be implanted with a bioabsorbable stent (BRS), or whether the first region of interest may be implanted with other medical devices.
In the medical image analysis method, the first parameter information corresponding to the first region of interest can be obtained by analyzing the first region of interest in the angiographic image to be analyzed, and then the angiographic image to be analyzed can be analyzed according to the first parameter information corresponding to the first region of interest to obtain an analysis result indicating whether the first region of interest in the angiographic image to be analyzed meets the preset first condition, because the method is used for analyzing the first region of interest in the angiographic image to be analyzed (coronary angiographic image), the angiographic image to be analyzed can be accurately analyzed according to the first parameter information corresponding to the first region of interest, so that the analysis result indicating whether the first region of interest in the angiographic image to be analyzed meets the preset first condition can be accurately obtained, the accuracy of the analysis result of the obtained angiographic image to be analyzed is improved.
In the obtaining the first region of interest in the angiographic image according to the angiographic image, in an embodiment, the S202 includes: inputting an angiography image (coronary angiography image) into a preset positioning network, and obtaining a first region of interest in the angiography image through the positioning network.
Specifically, the computer device inputs the angiographic image to be analyzed into a preset positioning network, and obtains a first region of interest in the angiographic image to be analyzed through the positioning network. Alternatively, the predetermined positioning network may be a U-Net network as shown in fig. 3. It can be understood that the preset positioning network is a pre-trained positioning network, and optionally, the computer device may train the preset initial positioning network through a large number of sample angiographic images (coronary angiographic images) and gold standard images of the region of interest corresponding to the sample angiographic images (coronary angiographic images) to obtain the positioning network.
In this embodiment, the computer device inputs the angiographic image to be analyzed into a preset positioning network, and since the preset positioning network is a pre-trained network, the first region of interest in the angiographic image to be analyzed can be accurately obtained through the positioning network; in addition, the first region of interest in the angiographic image to be analyzed can be quickly obtained through the preset positioning network, so that the efficiency of obtaining the first region of interest in the angiographic image to be analyzed is improved.
In a scenario where a first region of interest of the angiographic image to be analyzed is analyzed to obtain first parameter information corresponding to the first region of interest, and an analysis result indicating whether the first region of interest satisfies a preset first condition is obtained according to the first parameter information, the first parameter information includes: the position information of the first region of interest, the blood vessel distribution information in a preset range and the blood vessel shape information in the preset range in the first region of interest, and the plaque quantitative analysis information of the first region of interest.
Specifically, the first parameter information corresponding to the first region of interest includes: the position information of the first region of interest, the blood vessel distribution information in a preset range and the blood vessel shape information in the preset range in the first region of interest, and the plaque quantitative analysis information of the first region of interest. Accordingly, the computer device needs to analyze the position of the first region of interest, the blood vessel distribution within a preset range in the first region of interest and the blood vessel morphology within the preset range, and the plaque of the first region of interest, and the detailed description of the above analysis process and the specific process of obtaining the analysis result according to the first parameter information will be described below.
Specifically, analyzing the position of the first region of interest to obtain the position information of the first region of interest may include: taking the angiogram image to be analyzed and the heat map of the first region of interest as input, and transmitting the input to a trained classification network to obtain the position information of the first region of interest, wherein optionally, the trained classification network may be a classical classification network ResNet18 network. Optionally, if the position information of the first region of interest obtained by the computer device is that the first region of interest is located in the left trunk or the coronary opening, it is determined that the first region of interest does not satisfy the preset first condition, or if the position information of the first region of interest obtained by the computer device is that the restenosis/chronic total occlusion lesion in the bridged vascular lesion/Drug Eluting Stent (DES) in the first region of interest, it is determined that the first region of interest does not satisfy the preset first condition.
Specifically, the distribution of blood vessels in a preset range in the first region of interest and the morphology of the blood vessels in the preset range are analyzed, and the key is to determine the geometric morphology and the geometric parameters in the preset range in the first region of interest, including the reference diameters of the two ends of the blood vessels in the preset range in the first region of interest, whether the blood vessels are branched, the diameter of the branched blood vessels, whether the lesion is involved in branching, the degree of the involved branches, the curvature of the lesion blood vessels, the length of the lesion area, and the like. Optionally, the computer device may obtain the blood vessel distribution information within the preset range and the blood vessel morphology information within the preset range in the first region of interest by using a quantitative calculation method as shown in fig. 4, where the blood vessel extraction process may include: and constructing a graph with the image gradient of the angiogram as weight according to the preset range in the given first region of interest and the angiogram, and extracting the blood vessels in the angiogram by using the graph cut for blood vessel segmentation. The implementation process of quantitative analysis may include: calculating a blood vessel central line tree according to the blood vessel segmentation result, and then calculating a distribution diagram of the blood vessel diameter along the blood vessel central line tree; calculating reference diameters at two ends of the blood vessel, the length of the pathological changes of the blood vessel, the main branch stenosis degree, the blood vessel reconstruction line, the branch stenosis degree and the branch diameter according to the blood vessel diameter distribution map; and calculating the curvature of the main branch vessel according to the vessel centerline tree. If the blood vessel distribution information in the preset range and the blood vessel shape information in the preset range in the first region of interest obtained by the computer equipment meet the following conditions: the reference diameter of the blood vessel is beyond 2.75 mm-3.75 mm; bifurcated diseased vessels <2.0mm in diameter; bifurcated diseased vessels >2.0mm in diameter and > 50% stenosis; the lesion length is more than or equal to 20 mm; the difference between the reference diameters of the near end and the far end is greater than 0.25 mm; if any one of the blood vessel height meanders (the curvature of the blood vessel is larger than the artificially set threshold), the computer device determines that the first region of interest does not satisfy the preset first condition.
Specifically, the key to analyzing the plaque of the first region of interest is to determine the distribution of calcifications and the components of the plaque around the first region of interest, optionally, the computer device may analyze the plaque of the first region of interest through an analysis process as shown in fig. 5, optionally, the segmentation network in fig. 5 may be a U-Net network, and the gold standard image when the segmentation network is trained may be a mark of the sample image where various plaques are located. Optionally, the specific details of the quantitative analysis in fig. 5 are as follows: for the calcified area, calculating the distribution range of the calcified area around the lumen and the calcified thickness; for lipid and fibrous cap regions, the size is calculated to determine the plaque type. If the plaque quantitative analysis information of the first region of interest obtained by the computer equipment meets the following requirements: severe calcification (coverage >270 ° thickness >500 μm), lipid plaque or calcified plaque, the computer device determines that the first region of interest does not satisfy the predetermined first condition.
In this embodiment, the first parameter information corresponding to the first region of interest includes the position information of the first region of interest, the blood vessel distribution information within a preset range in the first region of interest, the blood vessel morphology information within the preset range, and the plaque quantitative analysis information of the first region of interest, so that the computer device can accurately analyze the first region of interest according to the first parameter information corresponding to the first region of interest, and the accuracy of an analysis result obtained by the computer device and used for indicating whether the first region of interest meets a preset first condition is improved.
In some scenarios, other types of intra-cavity imaging techniques besides angiograms may be used to assist in determining the nature of the current lesion, such as intravascular ultrasound (IVUS) and Optical Coherence Tomography (OCT), and in one embodiment, as shown in fig. 6, the method further comprises:
s601, acquiring an intra-cavity image corresponding to the angiogram image.
Specifically, the computer device obtains an intra-cavity image corresponding to the angiogram image to be analyzed. Optionally, the intraluminal image corresponding to the angiographic image to be analyzed may include an intravascular ultrasound image and optical coherence tomography. Optionally, the computer device may obtain an intra-cavity image corresponding to the angiogram image to be analyzed from a PACS (Picture Archiving and Communication Systems) server, or may obtain an intra-cavity image corresponding to the angiogram image to be analyzed in real time from a corresponding medical imaging device.
S602, in the intracavity image, a second region of interest corresponding to the first region of interest is determined.
It will be appreciated that when using an intra-luminal image corresponding to the angiogram image to be analyzed, it is necessary to map the region of interest in the angiogram with the region of interest in the intra-luminal image to obtain an intra-luminal cross-sectional view of the lesion site. Specifically, the computer device determines a second region of interest corresponding to the first region of interest in the intra-cavity image corresponding to the angiogram image to be analyzed. Optionally, the computer device may determine a second region of interest corresponding to the first region of interest in the intra-cavity image corresponding to the angiogram image by using an image registration method.
S603, analyzing the second region of interest to obtain second parameter information corresponding to the second region of interest.
Specifically, the computer device analyzes a second region of interest in the intra-cavity image corresponding to the angiogram image to obtain second parameter information corresponding to the second region of interest. Optionally, the computer device may analyze the position of the second region of interest, analyze the blood vessel distribution information in a preset range in the second region of interest, analyze the blood vessel morphology information in the preset range, and analyze the plaque amount in the second region of interest. Optionally, the second parameter information corresponding to the second region of interest may include: the position information of the second region of interest, the blood vessel distribution information in a preset range and the blood vessel shape information in the preset range in the second region of interest, and the plaque quantitative analysis information of the second region of interest. It should be noted that, for a specific process of obtaining, by the computer device, the second parameter information corresponding to the second region of interest, reference may be made to the specific description of obtaining, in the foregoing embodiment, the first parameter information corresponding to the first region of interest, which is not described herein again.
S604, obtaining an analysis result according to the first parameter information and the second parameter information.
Specifically, the computer device obtains an analysis result for indicating whether the first region of interest meets a preset first condition according to the first parameter information corresponding to the first region of interest and the second parameter information corresponding to the second region of interest. Optionally, if both the first parameter information corresponding to the first region of interest and the second parameter information corresponding to the second region of interest obtained by the computer device satisfy the conditions described in the above embodiment, the analysis result determined by the computer device indicates that the first region of interest satisfies the preset first condition; if any one of the first parameter information corresponding to the first interest and the second parameter information corresponding to the second interest, which are obtained by the computer device, does not satisfy the condition described in the above embodiment, the analysis result determined by the computer device indicates that the first interest does not satisfy the preset first condition.
In this embodiment, the computer device obtains an intra-cavity image corresponding to the angiographic image to be analyzed, determines a second region of interest corresponding to a first region of interest in the angiographic image to be analyzed in the intra-cavity image, so that the second region of interest can be analyzed to obtain second parameter information corresponding to the second region of interest, and further obtains an analysis result indicating whether the first region of interest in the angiographic image to be analyzed meets a preset first condition according to the first parameter information corresponding to the first region of interest and the second parameter information corresponding to the second region of interest, and therefore, according to the first parameter information corresponding to the first region of interest in the angiographic image to be analyzed and the second parameter information corresponding to the second region of interest in the intra-cavity image corresponding to the angiographic image to be analyzed Moreover, the analysis result for indicating whether the first region of interest in the angiographic image to be analyzed meets the preset first condition can be accurately obtained, and the accuracy of the obtained analysis result is improved.
In the above-mentioned scene of determining the second region of interest corresponding to the first region of interest in the angiographic image in the intra-cavity image corresponding to the angiographic image, the computer device may determine the second region of interest according to the blood vessel diameter distribution map, the blood vessel topological structure key points in the angiographic image, and the blood vessel diameter distribution map and the blood vessel topological structure key points in the intra-cavity image, as shown in fig. 7 and 7a, in an embodiment, the step S602 includes:
s701, calculating a first blood vessel diameter distribution map corresponding to the angiogram image and a second blood vessel diameter distribution map corresponding to the intracavity image.
Specifically, the computer device calculates a first blood vessel diameter distribution map corresponding to the angiogram image and a second blood vessel diameter distribution map corresponding to the intracavity image. Illustratively, the acquisition of the first vessel diameter distribution map may be obtained by a simple method exemplified by: the computer equipment calculates the center line of the target blood vessel, then makes the perpendicular line of the center line to obtain the intersection point position of the perpendicular line and the edge of the blood vessel, thereby obtaining the diameter value of the current center point, and draws the diameter value of the corresponding center point along the center line according to the physical length of the blood vessel from the near end to the far end, thus obtaining the diameter distribution map. Illustratively, the acquisition of the second vessel diameter map may be obtained by a simple method exemplified by: the computer device firstly obtains all diameters in each blood vessel section, then uses the minimum diameter as the diameter value of the current section position, and draws the diameter value of the corresponding central point along the central line according to the physical length of the blood vessel from the near end to the far end, thus obtaining the diameter distribution map.
S702, calculating a first blood vessel topological structure key point corresponding to the angiogram image and a second blood vessel topological structure key point corresponding to the intracavity image.
Specifically, the computer device calculates a first blood vessel topological structure key point corresponding to the angiogram image and a second blood vessel topological structure key point corresponding to the intra-cavity image. Optionally, there are many available methods for extracting key points of the topological structure of the blood vessel, and here, a simple and feasible scheme is exemplified: the extraction of the first vessel topological keypoints can be obtained using intersection point detection because it is a two-dimensional image. The extraction of the second blood vessel topological structure key point is the blood vessel section along the central line, so that whether the current section is a bifurcation or other structure can be judged through a simple classification network, and the topological information is obtained.
And S703, determining a second region of interest corresponding to the first region of interest according to the first blood vessel diameter distribution map, the second blood vessel diameter distribution map, the first blood vessel topological structure key points and the second blood vessel topological structure key points.
Specifically, the computer device determines a second region of interest corresponding to the first region of interest according to the obtained first blood vessel diameter distribution map, the obtained second blood vessel diameter distribution map, the obtained first blood vessel topological structure key points and the obtained second blood vessel topological structure key points. Optionally, the computer device may register the first blood vessel diameter distribution map and the second blood vessel diameter distribution map, adjust the first blood vessel topological structure key point and the second blood vessel topological structure key point until the first blood vessel topological structure key point and the second blood vessel topological structure key point coincide with each other when the first blood vessel diameter distribution map and the second blood vessel diameter distribution map coincide with each other, and map the first region of interest to the intra-cavity image when the first blood vessel topological structure key point and the second blood vessel topological structure key point coincide with each other, so as to obtain a second region of interest corresponding to the first region of interest.
In this embodiment, the computer device calculates the first blood vessel diameter distribution map corresponding to the angiogram image and the second blood vessel diameter distribution map corresponding to the intraluminal image, and the first blood vessel topological structure key point corresponding to the angiogram image and the second blood vessel topological structure key point corresponding to the intraluminal image, so that the second region of interest corresponding to the first region of interest can be accurately determined according to the first blood vessel diameter distribution map, the second blood vessel diameter distribution map, the first blood vessel topological structure key point and the second blood vessel topological structure key point, and the accuracy of the determined second region of interest corresponding to the first region of interest is improved.
In some scenarios, after obtaining an analysis result indicating whether the first region of interest in the angiographic image satisfies the preset first condition, the computer device may display the obtained analysis result in a report format as shown in fig. 8 in order to visually display the analysis result and facilitate the user to view the analysis result. It should be noted that the report format in fig. 8 is only an example, and the format of the report is not limited in this embodiment.
It should be understood that although the various steps in the flow charts of fig. 2-8 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. 2-8 may include multiple steps or multiple stages, which are not necessarily performed at the same time, but may be performed at different times, which are not necessarily performed in sequence, but may be performed in turn or alternately with other steps or at least some of the other steps.
The application of the degradable stent in the coronary artery is still a relatively new application at present, and at the present stage, clinicians mainly rely on experience for whether the degradable stent can be used for patients. Meanwhile, a large amount of manual measurement operation is required to measure coronary artery disease. This is a high requirement for the clinical experience, qualification and measurement technique of the doctor. And experience judges that reproducibility and consistency are generally poor. At present, a complete and mature tool is not available for assisting and supporting doctors to judge the implantation performance of the degradable stent. In order to solve the problems, the scheme in the embodiment of the application realizes the analysis of an angiography image (coronary angiography image), and facilitates the subsequent quick and reproducible evaluation of the implantation performance of the coronary degradable stent. The designed work flow reduces the requirements on the experience of doctors, and relatively comprehensive evaluation results can be obtained only by executing limited operations, so that the resources of the doctors are saved, and meanwhile, relatively stable services are provided for patients. The automatic/semi-automatic method allows the measurement results to be reproducible, and also provides accurate quantitative indicators for the doctor for further evaluation.
In one embodiment, as shown in fig. 9, there is provided a medical image analysis apparatus including: first acquisition module, second acquisition module, first analysis module and third acquisition module, wherein:
the first acquisition module is used for acquiring an angiographic image to be analyzed.
And the second acquisition module is used for obtaining a first region of interest in the angiography image according to the angiography image.
And the first analysis module is used for analyzing the first region of interest to obtain first parameter information corresponding to the first region of interest.
The third acquisition module is used for acquiring an analysis result according to the first parameter information; the analysis result is used for indicating whether the first region of interest meets a preset first condition.
Optionally, the first parameter information includes: the position information of the first region of interest, the blood vessel distribution information in a preset range and the blood vessel shape information in the preset range in the first region of interest, and the plaque quantitative analysis information of the first region of interest.
The medical image analysis apparatus provided in this embodiment may implement the method embodiments described above, and the implementation principle and the technical effect are similar, which are not described herein again.
On the basis of the foregoing embodiment, optionally, the second obtaining module includes: an acquisition unit, wherein:
the acquisition unit is used for inputting the angiography image into a preset positioning network and obtaining a first region of interest in the angiography image through the positioning network.
The medical image analysis apparatus provided in this embodiment may implement the method embodiments described above, and the implementation principle and the technical effect are similar, which are not described herein again.
On the basis of the foregoing embodiment, optionally, the apparatus further includes: the fourth acquisition module, the determination module, the second analysis module and the fifth acquisition module, wherein:
and the fourth acquisition module is used for acquiring the intracavity image corresponding to the angiogram image.
And the determining module is used for determining a second interested area corresponding to the first interested area in the intracavity image.
And the second analysis module is used for analyzing the second region of interest to obtain second parameter information corresponding to the second region of interest.
And the fifth acquisition module is used for acquiring an analysis result according to the first parameter information and the second parameter information.
Optionally, the second parameter information includes: the position information of the second region of interest, the blood vessel distribution information in a preset range and the blood vessel shape information in the preset range in the second region of interest, and the plaque quantitative analysis information of the second region of interest.
Optionally, the intraluminal images include intravascular ultrasound images and optical coherence tomography.
The medical image analysis apparatus provided in this embodiment may implement the method embodiments described above, and the implementation principle and the technical effect are similar, which are not described herein again.
On the basis of the foregoing embodiment, optionally, the determining module includes: a first calculation unit, a second calculation unit, and a determination unit, wherein:
the first calculating unit is used for calculating a first blood vessel diameter distribution map corresponding to the angiogram image and a second blood vessel diameter distribution map corresponding to the intracavity image.
And the second calculating unit is used for calculating a first blood vessel topological structure key point corresponding to the angiogram image and a second blood vessel topological structure key point corresponding to the intracavity image.
And the determining unit is used for determining a second region of interest corresponding to the first region of interest according to the first blood vessel diameter distribution map, the second blood vessel diameter distribution map, the first blood vessel topological structure key points and the second blood vessel topological structure key points.
The medical image analysis apparatus provided in this embodiment may implement the method embodiments described above, and the implementation principle and the technical effect are similar, which are not described herein again.
For specific limitations of the medical image analysis apparatus, reference may be made to the above limitations of the medical image analysis method, which are not described herein again. The modules in the medical image analysis apparatus may be implemented in whole or in part 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, 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 an angiographic image to be analyzed;
obtaining a first region of interest in the angiographic image according to the angiographic image;
analyzing the first region of interest to obtain first parameter information corresponding to the first region of interest;
obtaining an analysis result according to the first parameter information; the analysis result is used for indicating whether the first region of interest meets a preset first condition.
The implementation principle and technical effect of the computer device provided by the above embodiment are similar to those of the above method embodiment, and are not described herein again.
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 an angiographic image to be analyzed;
obtaining a first region of interest in the angiographic image according to the angiographic image;
analyzing the first region of interest to obtain first parameter information corresponding to the first region of interest;
obtaining an analysis result according to the first parameter information; the analysis result is used for indicating whether the first region of interest meets a preset first condition.
The implementation principle and technical effect of the computer-readable storage medium provided by the above embodiments are similar to those of the above method embodiments, and are not described herein again.
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 can include at least one of non-volatile and volatile memory. Non-volatile Memory may include Read-Only Memory (ROM), magnetic tape, floppy disk, flash Memory, optical storage, or the like. Volatile Memory can include Random Access Memory (RAM) or external cache Memory. By way of illustration and not limitation, RAM can take many forms, such as Static Random Access Memory (SRAM) or Dynamic Random Access Memory (DRAM), among others.
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 method of medical image analysis, the method comprising:
acquiring an angiographic image to be analyzed;
obtaining a first region of interest in the angiographic image according to the angiographic image;
analyzing the first region of interest to obtain first parameter information corresponding to the first region of interest;
obtaining an analysis result according to the first parameter information; the analysis result is used for indicating whether the first region of interest meets a preset first condition.
2. The method of claim 1, wherein obtaining the first region of interest in the angiographic image from the angiographic image comprises:
and inputting the angiogram image into a preset positioning network, and obtaining a first region of interest in the angiogram image through the positioning network.
3. The method according to claim 1 or 2, wherein the first parameter information comprises: the position information of the first region of interest, the blood vessel distribution information in a preset range in the first region of interest, the blood vessel shape information in the preset range, and the plaque quantitative analysis information of the first region of interest.
4. The method of claim 1, further comprising:
acquiring an intra-cavity image corresponding to the angiogram image;
determining a second region of interest corresponding to the first region of interest in the intracavity image;
analyzing the second region of interest to obtain second parameter information corresponding to the second region of interest;
and obtaining the analysis result according to the first parameter information and the second parameter information.
5. The method of claim 4, wherein said determining a second region of interest corresponding to said first region of interest in said intra-cavity image comprises:
calculating a first blood vessel diameter distribution map corresponding to the angiogram image and a second blood vessel diameter distribution map corresponding to the intracavity image;
calculating a first blood vessel topological structure key point corresponding to the angiogram image and a second blood vessel topological structure key point corresponding to the intracavity image;
and determining a second region of interest corresponding to the first region of interest according to the first blood vessel diameter distribution map, the second blood vessel diameter distribution map, the first blood vessel topological structure key points and the second blood vessel topological structure key points.
6. The method of claim 5, wherein the second parameter information comprises: the position information of the second region of interest, the blood vessel distribution information in a preset range in the second region of interest, the blood vessel shape information in the preset range, and the plaque quantitative analysis information of the second region of interest.
7. The method of claim 4, wherein the intraluminal images comprise intravascular ultrasound images and optical coherence tomography.
8. A medical image analysis apparatus, characterized in that the apparatus comprises:
the first acquisition module is used for acquiring an angiographic image to be analyzed;
the second acquisition module is used for obtaining a first region of interest in the angiography image according to the angiography image;
the first analysis module is used for analyzing the first region of interest to obtain first parameter information corresponding to the first region of interest;
the third acquisition module is used for acquiring an analysis result according to the first parameter information; the analysis result is used for indicating whether the first region of interest meets a preset first condition.
9. A computer device comprising a memory and a processor, the memory storing a computer program, characterized in that the processor, when executing the computer program, implements the steps of the method of any of claims 1 to 7.
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 |
|---|---|---|---|
| CN202110148160.0A CN113040812B (en) | 2021-02-03 | 2021-02-03 | Medical image analysis method, medical image analysis device, computer equipment and storage medium |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110148160.0A CN113040812B (en) | 2021-02-03 | 2021-02-03 | Medical image analysis method, medical image analysis device, computer equipment and storage medium |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113040812A true CN113040812A (en) | 2021-06-29 |
| CN113040812B CN113040812B (en) | 2023-08-22 |
Family
ID=76508596
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110148160.0A Active CN113040812B (en) | 2021-02-03 | 2021-02-03 | Medical image analysis method, medical image analysis device, computer equipment and storage medium |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113040812B (en) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113723406A (en) * | 2021-09-03 | 2021-11-30 | 乐普（北京）医疗器械股份有限公司 | Processing method and device for positioning bracket of coronary angiography image |
| CN113989231A (en) * | 2021-10-28 | 2022-01-28 | 上海联影医疗科技股份有限公司 | Method and device for determining kinetic parameters, computer equipment and storage medium |
| CN114638878A (en) * | 2022-03-18 | 2022-06-17 | 北京安德医智科技有限公司 | Two-dimensional echocardiogram pipe diameter detection method and device based on deep learning |
Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2004023102A2 (en) * | 2002-09-09 | 2004-03-18 | The Board Of Regents Of The University Of Texas System | Methods of detection of progression of stenosis |
| US20110071404A1 (en) * | 2009-09-23 | 2011-03-24 | Lightlab Imaging, Inc. | Lumen Morphology and Vascular Resistance Measurements Data Collection Systems, Apparatus and Methods |
| CN103462696A (en) * | 2013-09-17 | 2013-12-25 | 浙江大学 | Integrated intravascular OCT (optical coherence tomography) image and DSA (digital subtraction angiography) integrating online real-time treatment device |
| CN109716446A (en) * | 2016-09-28 | 2019-05-03 | 光学实验室成像公司 | Utilize the stent planning system and method for blood vessel presentation |
| CN110613441A (en) * | 2019-10-09 | 2019-12-27 | 南京沃福曼医疗科技有限公司 | Pressure or blood flow obtaining method based on hemodynamics principle |
| TWI698225B (en) * | 2019-06-11 | 2020-07-11 | 宏碁股份有限公司 | Blood vessel status evaluation method and blood vessel status evaluation device |
| CN111445449A (en) * | 2020-03-19 | 2020-07-24 | 上海联影智能医疗科技有限公司 | Region-of-interest classification method and device, computer equipment and storage medium |
| CN111523538A (en) * | 2020-04-14 | 2020-08-11 | 博动医学影像科技（上海）有限公司 | Blood vessel image processing method and system, computing device and storage medium |
| CN112107362A (en) * | 2020-08-24 | 2020-12-22 | 江苏大学 | Computer-assisted surgery design system for coronary heart disease |
-
2021
- 2021-02-03 CN CN202110148160.0A patent/CN113040812B/en active Active
Patent Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2004023102A2 (en) * | 2002-09-09 | 2004-03-18 | The Board Of Regents Of The University Of Texas System | Methods of detection of progression of stenosis |
| US20110071404A1 (en) * | 2009-09-23 | 2011-03-24 | Lightlab Imaging, Inc. | Lumen Morphology and Vascular Resistance Measurements Data Collection Systems, Apparatus and Methods |
| CN103462696A (en) * | 2013-09-17 | 2013-12-25 | 浙江大学 | Integrated intravascular OCT (optical coherence tomography) image and DSA (digital subtraction angiography) integrating online real-time treatment device |
| CN109716446A (en) * | 2016-09-28 | 2019-05-03 | 光学实验室成像公司 | Utilize the stent planning system and method for blood vessel presentation |
| TWI698225B (en) * | 2019-06-11 | 2020-07-11 | 宏碁股份有限公司 | Blood vessel status evaluation method and blood vessel status evaluation device |
| CN110613441A (en) * | 2019-10-09 | 2019-12-27 | 南京沃福曼医疗科技有限公司 | Pressure or blood flow obtaining method based on hemodynamics principle |
| CN111445449A (en) * | 2020-03-19 | 2020-07-24 | 上海联影智能医疗科技有限公司 | Region-of-interest classification method and device, computer equipment and storage medium |
| CN111523538A (en) * | 2020-04-14 | 2020-08-11 | 博动医学影像科技（上海）有限公司 | Blood vessel image processing method and system, computing device and storage medium |
| CN112107362A (en) * | 2020-08-24 | 2020-12-22 | 江苏大学 | Computer-assisted surgery design system for coronary heart disease |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113723406A (en) * | 2021-09-03 | 2021-11-30 | 乐普（北京）医疗器械股份有限公司 | Processing method and device for positioning bracket of coronary angiography image |
| CN113723406B (en) * | 2021-09-03 | 2023-07-18 | 乐普（北京）医疗器械股份有限公司 | Method and device for processing support positioning of coronary angiography image |
| CN113989231A (en) * | 2021-10-28 | 2022-01-28 | 上海联影医疗科技股份有限公司 | Method and device for determining kinetic parameters, computer equipment and storage medium |
| CN114638878A (en) * | 2022-03-18 | 2022-06-17 | 北京安德医智科技有限公司 | Two-dimensional echocardiogram pipe diameter detection method and device based on deep learning |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113040812B (en) | 2023-08-22 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN113040812B (en) | Medical image analysis method, medical image analysis device, computer equipment and storage medium | |
| US10803994B2 (en) | Vascular flow assessment | |
| JP6667999B2 (en) | Image processing apparatus, image processing method, and program | |
| US10052032B2 (en) | Stenosis therapy planning | |
| CN112967220B (en) | Computer-implemented method of evaluating CT data sets relating to perivascular tissue | |
| CN110889896A (en) | Method, device and system for obtaining angiostenosis lesion interval and three-dimensional synthesis | |
| JP2013154192A (en) | Lumen morphology and vascular resistance measurements data collection system, apparatus and method | |
| EP3125764A1 (en) | Processing apparatus and method for processing cardiac data of a living being | |
| CN110458837B (en) | Image post-processing method and device, electronic equipment and storage medium | |
| CN114119602B (en) | Method, apparatus and storage medium for object analysis of medical images | |
| US20220319004A1 (en) | Automatic vessel analysis from 2d images | |
| CN115565667A (en) | Prediction method, device, equipment and medium for success rate of CTO interventional therapy | |
| CN112472112B (en) | Method, system and storage medium for regulating vascular stenosis | |
| CN111627023B (en) | Method and device for generating coronary artery projection image and computer readable medium | |
| WO2021048838A2 (en) | Automated analysis of image data to determine fractional flow reserve | |
| CN116091444A (en) | Side branch evaluation method and device, storage medium and terminal | |
| US10182782B2 (en) | Evaluation apparatus, evaluation method, and evaluation program | |
| JP2024535030A (en) | Method, device, equipment and medium for determining FFR based on multimodal medical images | |
| CN114549487A (en) | Image analysis method, image analysis system, image analysis apparatus, and image analysis medium | |
| Adams et al. | Algorithm for quantifying advanced carotid artery atherosclerosis in humans using MRI and active contours | |
| CN115880310B (en) | Retina OCT fault segmentation method, device and equipment | |
| JP7560103B2 (en) | Method and apparatus for inspecting small deformation of lumen | |
| CN117476238B (en) | Computational fluid dynamics simulation analysis method and system for evaluating myocardial ischemia | |
| CN114266759A (en) | Image analysis method and device and computer equipment | |
| JP2022128173A (en) | Image processing device, image processing program and image 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 | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |