[go: up one dir, main page]

CN117409151B - Three-dimensional surface reconstruction method and device for two-dimensional intracardiac ultrasonic catheter image - Google Patents

Three-dimensional surface reconstruction method and device for two-dimensional intracardiac ultrasonic catheter image Download PDF

Info

Publication number
CN117409151B
CN117409151B CN202311713620.5A CN202311713620A CN117409151B CN 117409151 B CN117409151 B CN 117409151B CN 202311713620 A CN202311713620 A CN 202311713620A CN 117409151 B CN117409151 B CN 117409151B
Authority
CN
China
Prior art keywords
image
dimensional
ultrasonic
images
interpolation
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202311713620.5A
Other languages
Chinese (zh)
Other versions
CN117409151A (en
Inventor
王泽辉
徐乃旻
范鲁燕
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Jiangsu Tingsn Technology Co ltd
Original Assignee
Jiangsu Tingsn Technology Co ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Jiangsu Tingsn Technology Co ltd filed Critical Jiangsu Tingsn Technology Co ltd
Priority to CN202311713620.5A priority Critical patent/CN117409151B/en
Publication of CN117409151A publication Critical patent/CN117409151A/en
Application granted granted Critical
Publication of CN117409151B publication Critical patent/CN117409151B/en
Priority to PCT/CN2024/114827 priority patent/WO2025123781A1/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T17/00Three dimensional [3D] modelling, e.g. data description of 3D objects
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/02Neural networks
    • G06N3/04Architecture, e.g. interconnection topology
    • G06N3/0464Convolutional networks [CNN, ConvNet]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/02Neural networks
    • G06N3/08Learning methods
    • G06N3/084Backpropagation, e.g. using gradient descent
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/10Segmentation; Edge detection
    • G06T7/11Region-based segmentation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/10Segmentation; Edge detection
    • G06T7/12Edge-based segmentation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10028Range image; Depth image; 3D point clouds
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10132Ultrasound image
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/20Special algorithmic details
    • G06T2207/20081Training; Learning
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/20Special algorithmic details
    • G06T2207/20084Artificial neural networks [ANN]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/30Subject of image; Context of image processing
    • G06T2207/30004Biomedical image processing
    • G06T2207/30048Heart; Cardiac

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Software Systems (AREA)
  • Artificial Intelligence (AREA)
  • Computing Systems (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Mathematical Physics (AREA)
  • Biomedical Technology (AREA)
  • Biophysics (AREA)
  • Computational Linguistics (AREA)
  • Data Mining & Analysis (AREA)
  • Evolutionary Computation (AREA)
  • General Health & Medical Sciences (AREA)
  • Molecular Biology (AREA)
  • General Engineering & Computer Science (AREA)
  • Geometry (AREA)
  • Computer Graphics (AREA)
  • Ultra Sonic Daignosis Equipment (AREA)

Abstract

本发明属于心超图像处理技术技术领域,提供了一种二维心内超声导管图像的三维表面重建方法及装置。本重建方法包括以下步骤:S1:在心脏内部采集多个连续的超声图像;S2:对各超声图像进行插值处理,获得多个插值图像;S3:将各超声图像输入预训练的Unet模型中,获得各超声图像的心脏腔室分隔图像;S4:提取各心脏腔室分隔图像的各腔室表面的边缘信息,获得边缘图像;S5:根据各边缘图像和插值图像,生成对应的点云数据;S6:通过Ashape算法对生成的点云数据进行表面重建,获得三维模型。本三维表面重建方法是基于二维图像和单一超声探头的超声图像建模方法,支持实时更新的心脏腔室重建,实现心腔的位置自动分割及建模,速度更快,准确性更高。

The present invention belongs to the technical field of cardiac ultrasound image processing technology, and provides a three-dimensional surface reconstruction method and device for a two-dimensional intracardiac ultrasound catheter image. The reconstruction method comprises the following steps: S1: collecting multiple continuous ultrasound images inside the heart; S2: interpolating each ultrasound image to obtain multiple interpolated images; S3: inputting each ultrasound image into a pre-trained Unet model to obtain a cardiac chamber partition image of each ultrasound image; S4: extracting edge information of each chamber surface of each cardiac chamber partition image to obtain an edge image; S5: generating corresponding point cloud data according to each edge image and the interpolated image; S6: performing surface reconstruction on the generated point cloud data through the Ashape algorithm to obtain a three-dimensional model. The three-dimensional surface reconstruction method is an ultrasound image modeling method based on two-dimensional images and a single ultrasound probe, supports real-time updated cardiac chamber reconstruction, realizes automatic segmentation and modeling of cardiac cavity positions, and has faster speed and higher accuracy.

Description

Three-dimensional surface reconstruction method and device for two-dimensional intracardiac ultrasonic catheter image
Technical Field
The invention relates to the technical field of cardiac ultrasonic image processing, in particular to a three-dimensional surface reconstruction method and device for a two-dimensional intracardiac ultrasonic catheter image.
Background
Intracardiac ultrasound catheter images are widely used in cardiac disease diagnosis and treatment procedures. However, the reconstruction of these images remains a challenging task due to image noise and resolution limitations.
In complex scenarios, model reconstruction often depends on labeling by professionals, and the accuracy depends on experience and ability of labeling staff, which is time-consuming and laborious, resulting in missing window periods for optimal treatment of patients, and thus is extremely important for automatic and accurate modeling of models.
In recent years, there have been proposed methods for reconstructing an ultrasound image, but most have limitations such as: the method for constructing the three-dimensional ultrasonic imaging catheter and the system in the heart has the problems that the constructed model is rough, the chamber in the heart cannot be reconstructed and the like. It is difficult to provide more accurate and comprehensive information in the diagnosis and treatment of cardiovascular disease.
Disclosure of Invention
In view of the above, the embodiments of the present invention provide a three-dimensional surface reconstruction method, apparatus and electronic device for two-dimensional intracardiac ultrasound catheter images, which are used for solving or partially solving the above-mentioned problems.
In a first aspect, an embodiment of the present invention provides a three-dimensional surface reconstruction method for a two-dimensional intracardiac ultrasound catheter image, including the steps of:
s1, acquiring a plurality of continuous ultrasonic images inside a heart through an intracardiac ultrasonic catheter;
s2, carrying out interpolation processing on blank areas among the ultrasonic images to obtain a plurality of interpolation images;
s3, inputting each ultrasonic image into a pre-trained Unet model to obtain a heart chamber separation image of each output ultrasonic image;
s4, extracting edge information of each chamber surface of each heart chamber separation image to obtain an edge image of each chamber;
s5, generating corresponding point cloud data according to each edge image and each interpolation image;
And S6, carrying out surface reconstruction on the generated point cloud data through Ashape algorithm to obtain a three-dimensional model.
According to a specific implementation manner of the embodiment of the present invention, in the step S1, the ultrasound image includes rotation angle information corresponding to an intracardiac ultrasound catheter.
According to a specific implementation manner of the embodiment of the present invention, in the step S2, the interpolation process is linear interpolation, and the interpolated image includes pixel information and angle information.
According to a specific implementation manner of the embodiment of the present invention, in the step S3, the Unet model training manner is as follows:
s3.1: preparing an intra-cardiac ultrasound image dataset with corresponding cardiac chamber annotations;
s3.2: preprocessing an intra-cardiac ultrasound image within the intra-cardiac ultrasound image dataset to remove noise, enhance contrast, and unify the dimensions of the image;
S3.3: dividing the preprocessed intracardiac ultrasonic image data set into a training set, a verification set and a test set;
s3.4: training the Unet model using the intra-cardiac ultrasound image data and corresponding cardiac chamber annotations in the training set;
S3.5: verifying the Unet model in the training process by using the verification set to evaluate the performance of the model and perform tuning;
s3.6: the test set is used to make a final assessment of the Unet model that was trained and validated.
According to a specific implementation manner of the embodiment of the present invention, in the step S3.4, the difference between the prediction result and the true label is measured using the cross entropy loss function, and the training process uses the gradient descent algorithm to minimize the loss function by back-propagating the parameters of the optimization model.
According to a specific implementation manner of the embodiment of the present invention, in the step S3, the method further includes: and the Unet model judges each output heart chamber separation image, and if the heart chamber separation image is judged to be undetermined, the corresponding heart chamber separation image is marked secondarily by manual intervention.
According to a specific implementation manner of the embodiment of the present invention, in the step S4, the edge information is extracted based on the Sobel operator in the following manner:
for each pixel location in the heart chamber separation image Amplitude of gradientAnd direction ofThe calculation can be made by the following formula:
Wherein, AndRepresenting the gradients of the pixels in the image in the x-direction and the y-direction, respectively.
According to a specific implementation manner of the embodiment of the present invention, in the step S5, the point cloud data generation manner is as follows:
Wherein, For the point coordinates in the point cloud, (x, y) for each pixel position in the edge image or the interpolated image,And c is an offset constant for angle information corresponding to the edge image or the interpolation image.
The embodiment of the method has at least the following technical effects:
The first and the current conventional heart chamber modeling depends on the second catheter, and the three-dimensional surface reconstruction method is an ultrasonic image modeling method based on a two-dimensional image and a single ultrasonic probe and based on a two-dimensional image turning point cloud, so that the operation flow and difficulty are reduced, and simultaneously, the heart chamber reconstruction updated in real time is supported, so that the automatic segmentation of the heart chamber position can be realized, the speed is higher, and the accuracy is higher.
And secondly, in the three-dimensional surface reconstruction method, the automatic segmentation result is adjusted in a man-machine interaction mode, and intervention is performed at any time to ensure accuracy and reliability.
Thirdly, in the three-dimensional surface reconstruction method, boundary contours are extracted through edges after the heart chamber separation image is segmented, cloud data are turned on, so that the calculated amount is greatly reduced, and the efficiency is improved.
In a second aspect, an embodiment of the present invention provides a three-dimensional surface reconstruction apparatus for two-dimensional intracardiac ultrasound catheter images, including:
The data acquisition module is used for acquiring a plurality of continuous ultrasonic images inside the heart through an intracardiac ultrasonic catheter;
the image interpolation module is used for carrying out interpolation processing on blank areas among the ultrasonic images to obtain a plurality of interpolation images;
The region separation module is used for inputting each ultrasonic image into a pre-trained Unet model to obtain a heart chamber separation image of each output ultrasonic image;
The edge extraction module is used for extracting edge information of each chamber surface of each heart chamber separation image to obtain an edge image of each chamber;
the point cloud generation module is used for generating corresponding point cloud data according to the edge images and the interpolation images;
and the surface reconstruction module is used for carrying out surface reconstruction on the generated point cloud data through Ashape algorithm to obtain a three-dimensional model.
In a third aspect, an embodiment of the present invention further provides an electronic device, where the electronic device includes a memory, a processor, and a computer program stored on the memory and executable on the processor, where the processor implements the steps of the three-dimensional surface reconstruction method in the foregoing first aspect or any implementation manner of the first aspect when the processor executes the program.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings used in the description of the embodiments or the prior art will be briefly described below. Like elements or portions are generally identified by like reference numerals throughout the several figures. In the drawings, elements or portions thereof are not necessarily drawn to scale.
FIG. 1 is a flow chart showing the steps of a method for reconstructing a three-dimensional surface of a two-dimensional intracardiac ultrasound catheter image according to an embodiment of the present invention;
FIG. 2 shows a schematic view of ultrasound image acquisition in step S1 according to an embodiment of the present invention, wherein A is a schematic view of a sequence of acquiring continuous intra-cardiac ultrasound images, and B is a schematic view of a single intra-cardiac ultrasound image;
FIG. 3 is a schematic view of image interpolation in step S2 according to an embodiment of the present invention, wherein A, B is two adjacent two-dimensional ultrasound images, and C is an interpolated image after A, B is linearly interpolated;
Fig. 4 shows a schematic diagram of a segmentation result of Unet in step S3 according to an embodiment of the present invention, wherein A, B, C is a schematic diagram of a left heart chamber at different angles;
fig. 5 is a schematic diagram showing the edge extraction result in step S4 according to the embodiment of the present invention, wherein A, B, C is the boundary contour information of the left chamber of the heart at different angles corresponding to fig. 4;
Fig. 6 shows a schematic diagram of point cloud generation in step S5 according to an embodiment of the present invention, where a is a schematic diagram of original point cloud data, and B is a schematic diagram after FPS sampling;
FIG. 7 is a schematic diagram showing the effect of the Ashape algorithm surface reconstruction in step S6 according to an embodiment of the present invention;
FIG. 8 shows an ultrasound device connection schematic;
fig. 9 is a schematic structural diagram of an electronic device according to an embodiment of the present invention;
FIG. 10 is a block diagram of a three-dimensional surface reconstruction apparatus for two-dimensional intra-cardiac ultrasound catheter images according to an embodiment of the present invention;
In fig. 8, 1, ultrasonic transducer, 2, ultrasonic catheter, 3, ultrasonic handle, 4, rotary motor, 5, ultrasonic connector, 6, ultrasonic mainframe.
Detailed Description
Embodiments of the technical scheme of the present invention will be described in detail below with reference to the accompanying drawings. The following examples are only for more clearly illustrating the technical aspects of the present invention, and thus are merely examples, which should not be construed as limiting the scope of the present invention.
It is noted that unless otherwise indicated, technical or scientific terms used herein should be given the ordinary meaning as understood by one of ordinary skill in the art to which this application belongs.
Fig. 1 is a flowchart of steps of a three-dimensional surface reconstruction method for a two-dimensional intracardiac ultrasound catheter image according to an embodiment of the present invention, referring to fig. 1, the method includes the following steps:
s1, data acquisition: a plurality of sequential ultrasound images are acquired inside the heart through an intracardiac ultrasound catheter.
Continuous ultrasound images are acquired using an intracardiac ultrasound catheter technique, specifically by placing the catheter inside the heart and rotating it through a range of angles along the ultrasound catheter head. As shown in A in fig. 2, a continuous intracardiac ultrasonic image sequence is acquired, in the process, the motor rotates step by step, and waits for sampling instructions sent by an electrocardiogram to sample, and meanwhile, the rotation angles of all moments are recorded, and the real-time internal 360-degree high-speed rotation can be supported during rotation. In order to ensure that the acquired image sequence covers all view angles and positions of the heart chamber and has good image quality, the sampling frequency in the acquisition process needs to be sampled according to the diastolic points of the electrocardiogram of the patient so as to ensure that the three-dimensional model drawn after sampling is accurate, and meanwhile, the change of the whole process can be dynamically refreshed through real-time model reconstruction. The individual acquisition results are shown as B in fig. 2.
S2, image interpolation: and carrying out interpolation processing on blank areas among the ultrasonic images to obtain a plurality of interpolation images.
For the acquired two-dimensional ultrasonic image data, since the motor rotation angle cannot meet the minimum modeling requirement, interpolation is required for a blank area between images. And interpolating the edge image to obtain an interpolation image, wherein the interpolation image is used for filling a blank area which possibly appears, and comprises pixel information and angles.
Interpolation is to ensure not only continuity of the image but also reliability of the whole image, so that the pixel value thereof is calculated using linear interpolation. The linear interpolation calculates the pixel value of the interpolation point according to the distance ratio by using the pixel values of two adjacent points on the boundary of the known region. Let the point to be interpolated be P and its neighboring known points be a and B. The distance between the point A and the point P is d1, and the distance between the point B and the point P is d2. The interpolated pixel value for the P point can be calculated by the following formula:
P = A + (B - A) * (d1 / (d1 + d2))(1)
this formula is actually based on the linear relationship of the line segment AB, and the value of P is calculated according to the distance scale.
Fig. 3 is an effect diagram of interpolation by the above method, A, B in fig. 3 is two adjacent diagrams obtained by ultrasound, and C in fig. 3 is a result after linear interpolation from A, B in fig. 3.
S3, target area separation: each ultrasound image is input into a pre-trained Unet model, and a heart chamber separation image of each ultrasound image is obtained.
The Unet model is a common convolutional neural network architecture, and is particularly suitable for pixel-level segmentation tasks. The design inspiration comes from the self-encoder, and the local characteristics and the context information can be learned at the same time by connecting the encoder and the decoder, so that an accurate segmentation result is realized.
The Unet model training mode is as follows:
s3.1 data preparation:
An intracardiac ultrasound image dataset with corresponding cardiac chamber annotations is first prepared. Typically, these images will contain individual slices and views of the heart, as well as corresponding heart chamber boundary annotations.
S3.2, preprocessing data:
an intra-cardiac ultrasound image within the image dataset is preprocessed to remove noise, enhance contrast, and unify the dimensions of the image. Common preprocessing operations include filtering, histogram equalization, normalization, and the like.
S3.3, data division:
the preprocessed image dataset is divided into a training set, a validation set and a test set. The training set is used for training of the model, the validation set is used for adjusting the hyper-parameters of the model and monitoring the training progress, and the test set is used for evaluating the performance of the model.
S3.4 training Unet model:
Unet, consisting of an encoder and decoder, can learn to perform pixel-level cardiac chamber segmentation. The encoder portion is used to extract image features and the decoder portion is used to remap features back to the size of the input image and generate segmentation results. The Unet model is trained using intra-cardiac ultrasound image data and corresponding cardiac chamber labeling in a training set. The difference between the predicted outcome and the true annotation is typically measured using a cross entropy loss function and the parameters of the model are optimized by back propagation. The training process may use a gradient descent algorithm or variants thereof to minimize the loss function.
S3.5, verification and optimization:
And verifying the Unet model in the training process by using a verification set to evaluate the performance of the model and perform tuning. The best model may be selected based on segmentation results and metrics (e.g., accuracy, recall, F1 score, etc.) on the validation set.
S3.6 test and evaluation:
the trained and validated Unet model was finally evaluated using the test set. And inputting the intracardiac ultrasonic image on the test set into the trained model to obtain a heart chamber segmentation result generated by the model. Various evaluation indices (e.g., dice coefficients, jaccard indices, etc.) may be used to evaluate the performance and accuracy of the model.
The heart chamber separation image after Unet outputs is shown in fig. 4, wherein the three images A, B, C in fig. 4 are respectively left heart chambers at different angles.
The steps also include:
man-machine interaction: the part sets a manual labeling button in the interface aiming at the heart chamber image with deviation of labeling so as to prevent errors of the machine.
The method comprises the specific operations that a result of automatic segmentation is presented in a screen, unet can score each picture, the scoring is based on the internal consistency of two indexes 1), the consistency and smoothness of the interior of a segmented image are evaluated, the quality of the segmented result can be evaluated by calculating the similarity or difference between adjacent pixels, for example, the difference of color, texture or gradient between pixels can be calculated, and the segmented image is scored according to the size of the difference; 2) Regional connectivity: whether the areas in the segmented image have connectivity or not is evaluated, the smoothness of the area, perimeter or boundary of each area can be calculated, and different scores can be given to the segmented results according to the characteristics. After the evaluation indexes are weighted and averaged, the images are arranged in order from small to large, and a corresponding user is prompted to check the images under the condition that the images are lower than a set threshold value, and if the images are inaccurate, the corresponding images can be marked secondarily in a manual intervention mode.
S4, edge extraction: and extracting edge information of each chamber surface of each heart chamber separation image, and obtaining an edge image of each chamber.
The edge extraction method is used for obtaining edge information of the surface of the cavity. The boundary contours of the chamber and other regions are obtained by edge extraction and used for further processing and analysis. Edge extraction can help determine the pixel location of the object boundary, enabling accurate segmentation and reconstruction of the chamber region.
In this embodiment, a method based on the Sobel operator is used, which is based on the principle that the gradient of pixel intensity in an image is calculated to detect the edge, so as to calculate the image positionAndGradient in direction. Assuming that in the Sobel operator, each pixel position in the input imageThen the amplitude of the image gradientAnd direction ofThe calculation can be made by the following formula:
(2)
(3)
Wherein, AndRepresenting the gradient of the image in the x-direction and the y-direction, respectively, the larger the G value is, the more significant the pixel value near the edge will be, i.e. there is edge information.
Unet, and calculating boundary contour information of the heart chamber separation diagram through a Sobel operator, as shown in fig. 5.
S5, generating point cloud: and generating corresponding point cloud data according to each edge image and the interpolation image.
And according to each edge image and the angle obtained during scanning (corresponding to the ultrasonic image), generating corresponding point cloud data by combining the pixel information and the angle of the interpolation image in the step S2.
Assuming that the original coordinates are (x, y) and the pixel position on the nth picture is (x, y) and the angle isAssume that the transformed point cloud coordinates areThe transformation relationship is shown in formula (4).
(4)
Where c is the offset constant.
The point cloud generated by the method has the condition of excessive point number, which causes the calculation amount of the subsequent steps to be increased, and in order to solve the problem, the furthest point sampling (Farthest Point Sampling, FPS) is introduced for selecting a group of representative sampling points from large-scale point cloud data. The furthest point sampling method is based on the distance measurement between points, and covers the whole point cloud data as much as possible while keeping the uniform distribution among the sampling points. In this embodiment, the furthest point sampling is adopted to keep uniformity between data on one hand, and on the other hand, to avoid the problem that the triangular surface difference is too large during the later surface reconstruction, as shown in fig. 6, a in fig. 6 is original point cloud data, the number of point clouds is large, so that the calculation complexity of the surface reconstruction is increased, and B in fig. 6 is obtained after sampling the FPS for 10000 points, so that the improvement of the calculation performance is obtained while the integral structure is not lost.
S6, surface reconstruction: and carrying out surface reconstruction on the generated point cloud data through Ashape algorithm to obtain a three-dimensional model.
The Ashape algorithm is a commonly used surface reconstruction algorithm that can recover a smooth chamber surface based on the geometry and topology of the point cloud data. Three-dimensional shape information of the chamber can be extracted from the point cloud data through Ashape algorithm.
The generated point cloud data is surface reconstructed using Ashape algorithm. And visually displaying the chamber data subjected to Ashape algorithm surface reconstruction. Three-dimensional visualization software or libraries may be used to present the chamber reconstruction results so that doctors and researchers can more intuitively observe and analyze the heart chamber structure, as shown in fig. 7.
It should be noted that, the modules are arranged according to a streaming layout, which is only one embodiment of the present invention, and may be arranged in other manners, which is not limited by the present invention.
Different from other intracardiac catheter image modeling, the three-dimensional surface reconstruction method is a reconstruction method based on an ultrasonic image completely, firstly, an image of the heart inner cavity is acquired through the ultrasonic catheter, then a cavity boundary is determined through interpolation and automatic segmentation methods, then three-dimensional coordinate mapping is carried out through the image after edge extraction, and finally, surface reconstruction is carried out through AShape algorithm. The embodiment of the invention has at least the following technical effects:
The first and the current conventional heart chamber modeling depends on the second catheter, and the three-dimensional surface reconstruction method is an ultrasonic image modeling method based on a two-dimensional image and a single ultrasonic probe and based on a two-dimensional image turning point cloud, so that the operation flow and difficulty are reduced, and simultaneously, the heart chamber reconstruction updated in real time is supported, so that the automatic segmentation of the heart chamber position can be realized, the speed is higher, and the accuracy is higher.
And secondly, in the three-dimensional surface reconstruction method, the automatic segmentation result is adjusted in a man-machine interaction mode, and intervention is performed at any time to ensure accuracy and reliability.
Thirdly, in the three-dimensional surface reconstruction method, boundary contours are extracted through edges after the heart chamber separation image is segmented, cloud data are turned on, so that the calculated amount is greatly reduced, and the efficiency is improved.
The embodiment is based on an intracardiac ultrasonic catheter technology, the clinical medical diagnosis technology of ultrasonic waves is an ultrasonic detection technology based on echo scanning, ultrasonic diagnostic imaging adopts an array transducer with multiple array elements to emit ultrasonic waves into a human body, changes the relative delay and amplitude of excitation of each array element, and can form a focused acoustic beam emitted in a certain direction. When the acoustic beam encounters the interface of different organs and tissues in the body, a reflected echo is generated and received by the array transducer. The above process requires multiple component connections, the device connections being as shown in fig. 8. The device specifically comprises the following components:
The ultrasonic transducer 1 is composed of 64 array elements, completes the functions of ultrasonic wave emission and receiving, and can be sent to a designated position of the heart for data acquisition.
An ultrasound catheter 2, in which various connection lines for the transducer and a flexible material for bending are contained.
An ultrasonic handle 3, which is operated by a doctor to control the depth and direction of the catheter, and a sensor such as a rotary motor is arranged in the ultrasonic handle.
And a rotary motor 4 connected to the catheter for adjusting the direction of the ultrasonic transducer.
An ultrasonic connector 5, which converts between analog and digital signals.
The ultrasonic host computer 6 sends or receives the instruction, and the two-dimensional image sequence is obtained through the ultrasonic image acquisition equipment, so that faster visual feedback is provided for doctors.
Fig. 9 shows a schematic structural diagram of an electronic device 90 according to an embodiment of the present invention, where the electronic device 90 includes at least one processor 901 (e.g. a CPU), at least one input/output interface 904, a memory 902, and at least one communication bus 903 for enabling connection communication between these components. At least one processor 901 is configured to execute computer instructions stored in memory 902 to enable the at least one processor 901 to perform any one of the embodiments of the three-dimensional surface reconstruction method described above. The memory 902 is a non-transitory memory (non-transitory memory) that may include volatile memory, such as high-speed random access memory (RAM: random Access Memory), or may include non-volatile memory, such as at least one disk memory. Communication connection(s) with at least one other device or unit is effected through at least one input output interface 904 (which may be a wired or wireless communication interface).
In some embodiments, the memory 902 stores a program 9021, and the processor 901 executes the program 9021 to perform any of the foregoing sub-table method embodiments.
The electronic device may exist in a variety of forms including, but not limited to:
(1) A mobile communication device: such devices are characterized by mobile communication capabilities and are primarily aimed at providing voice, data communications. Such terminals include: smart phones, multimedia phones, functional phones, low-end phones, etc.
(2) Ultra mobile personal computer device: such devices are in the category of personal computers, having computing and processing functions, and generally also having mobile internet access characteristics. Such terminals include: PDA, MID, and UMPC devices, etc., such as iPad.
(3) Portable entertainment device: such devices may display and play multimedia content. The device comprises: audio, video players (e.g., iPod), palm game consoles, electronic books, and smart toys and portable car navigation devices.
(4) Specific server: the configuration of the server includes a processor, a hard disk, a memory, a system bus, and the like, and the server is similar to a general computer architecture, but is required to provide highly reliable services, and thus has high requirements in terms of processing capacity, stability, reliability, security, scalability, manageability, and the like.
(5) Other electronic devices with data interaction functions.
Fig. 10 is a block diagram of a three-dimensional surface reconstruction apparatus for two-dimensional intracardiac ultrasound catheter images according to an embodiment of the present invention, the apparatus comprising:
The data acquisition module is used for acquiring a plurality of continuous ultrasonic images in the heart through the intracardiac ultrasonic catheter;
the image interpolation module is used for carrying out interpolation processing on blank areas among the ultrasonic images to obtain a plurality of interpolation images;
The region separation module is used for inputting each ultrasonic image into the pre-trained Unet model to obtain a heart chamber separation image of each output ultrasonic image;
The edge extraction module is used for extracting edge information of the surfaces of the chambers of the heart chamber separation images to obtain edge images of the chambers;
the point cloud generation module is used for generating corresponding point cloud data according to each edge image and the interpolation image;
And the surface reconstruction module is used for carrying out surface reconstruction on the generated point cloud data through Ashape algorithm to obtain a three-dimensional model.
The functions of each module in the embodiment of fig. 10 correspond to the content in the corresponding method embodiment, and are not described herein.
The foregoing is merely illustrative of the present invention, and the present invention is not limited thereto, and any changes or substitutions easily contemplated by those skilled in the art within the scope of the present invention should be included in the present invention. Therefore, the protection scope of the invention is subject to the protection scope of the claims.

Claims (7)

1. A method for reconstructing a three-dimensional surface of a two-dimensional intracardiac ultrasound catheter image, comprising the steps of:
s1, acquiring a plurality of continuous ultrasonic images through the intracardiac ultrasonic catheter in a rotating way, wherein the ultrasonic images comprise rotation angle information corresponding to the intracardiac ultrasonic catheter;
S2, carrying out interpolation processing on blank areas among the ultrasonic images to obtain a plurality of interpolation images, wherein the interpolation processing is linear interpolation, and the interpolation images comprise pixel information and angle information;
S3, inputting each ultrasonic image into a pre-trained Unet model to obtain a heart chamber separation image of each ultrasonic image output by the Unet model;
s4, extracting edge information of each chamber surface of each heart chamber separation image to obtain an edge image of each chamber;
s5, generating corresponding point cloud data according to the edge images and the interpolation images, wherein the point cloud data is generated in the following mode:
Wherein, For the point coordinates in the point cloud, (x, y) for each pixel position in the edge image or the interpolated image,C is an offset constant for angle information corresponding to the edge image or the interpolation image;
And S6, carrying out surface reconstruction on the generated point cloud data through Ashape algorithm to obtain a three-dimensional model.
2. The three-dimensional surface reconstruction method according to claim 1, wherein in the step S3, the Unet model training method is as follows:
s3.1: preparing an intra-cardiac ultrasound image dataset with corresponding cardiac chamber annotations;
s3.2: preprocessing an intra-cardiac ultrasound image within the intra-cardiac ultrasound image dataset to remove noise, enhance contrast, and unify the dimensions of the image;
S3.3: dividing the preprocessed intracardiac ultrasonic image data set into a training set, a verification set and a test set;
s3.4: training the Unet model using the intra-cardiac ultrasound image data and corresponding cardiac chamber annotations in the training set;
S3.5: verifying the Unet model in the training process by using the verification set to evaluate the performance of the model and perform tuning;
s3.6: the test set is used to make a final assessment of the Unet model that was trained and validated.
3. The three-dimensional surface reconstruction method according to claim 2, wherein in step S3.4, the difference between the prediction result and the true annotation is measured using a cross entropy loss function, and the training process uses a gradient descent algorithm to minimize the loss function by back-propagating the parameters of the optimization model.
4. The three-dimensional surface reconstruction method according to claim 2, wherein in the step S3, further comprising: and the Unet model judges each output heart chamber separation image, and if the heart chamber separation image is judged to be undetermined, the corresponding heart chamber separation image is marked secondarily by manual intervention.
5. The three-dimensional surface reconstruction method according to claim 1, wherein in the step S4, the edge information is extracted based on a Sobel operator in the following manner:
for each pixel location in the heart chamber separation image Amplitude of gradientAnd direction ofThe calculation can be made by the following formula:
Wherein, AndRepresenting the gradients of the pixels in the image in the x-direction and the y-direction, respectively.
6. A three-dimensional surface reconstruction device for two-dimensional intracardiac ultrasound catheter images, comprising:
The data acquisition module is used for acquiring a plurality of continuous ultrasonic images in the heart through the intracardiac ultrasonic catheter, and the ultrasonic images comprise rotation angle information corresponding to the intracardiac ultrasonic catheter;
the image interpolation module is used for carrying out interpolation processing on blank areas among the ultrasonic images to obtain a plurality of interpolation images, wherein the interpolation processing is linear interpolation, and the interpolation images comprise pixel information and angle information;
The region separation module is used for inputting each ultrasonic image into a pre-trained Unet model to obtain a heart chamber separation image of each output ultrasonic image;
The edge extraction module is used for extracting edge information of each chamber surface of each heart chamber separation image to obtain an edge image of each chamber;
the point cloud generation module is used for generating corresponding point cloud data according to each edge image and the interpolation image, and the point cloud data generation mode is as follows:
Wherein, For the point coordinates in the point cloud, (x, y) for each pixel position in the edge image or the interpolated image,C is an offset constant for angle information corresponding to the edge image or the interpolation image;
and the surface reconstruction module is used for carrying out surface reconstruction on the generated point cloud data through Ashape algorithm to obtain a three-dimensional model.
7. An electronic device comprising a memory, a processor and a computer program stored on the memory and executable on the processor, characterized in that the processor implements the steps of the three-dimensional surface reconstruction method according to any one of claims 1 to 5 when the program is executed by the processor.
CN202311713620.5A 2023-12-14 2023-12-14 Three-dimensional surface reconstruction method and device for two-dimensional intracardiac ultrasonic catheter image Active CN117409151B (en)

Priority Applications (2)

Application Number Priority Date Filing Date Title
CN202311713620.5A CN117409151B (en) 2023-12-14 2023-12-14 Three-dimensional surface reconstruction method and device for two-dimensional intracardiac ultrasonic catheter image
PCT/CN2024/114827 WO2025123781A1 (en) 2023-12-14 2024-08-27 Three-dimensional surface reconstruction method and apparatus for two-dimensional intracardiac ultrasound catheter image

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202311713620.5A CN117409151B (en) 2023-12-14 2023-12-14 Three-dimensional surface reconstruction method and device for two-dimensional intracardiac ultrasonic catheter image

Publications (2)

Publication Number Publication Date
CN117409151A CN117409151A (en) 2024-01-16
CN117409151B true CN117409151B (en) 2024-06-28

Family

ID=89491153

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202311713620.5A Active CN117409151B (en) 2023-12-14 2023-12-14 Three-dimensional surface reconstruction method and device for two-dimensional intracardiac ultrasonic catheter image

Country Status (2)

Country Link
CN (1) CN117409151B (en)
WO (1) WO2025123781A1 (en)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117409151B (en) * 2023-12-14 2024-06-28 江苏霆升科技有限公司 Three-dimensional surface reconstruction method and device for two-dimensional intracardiac ultrasonic catheter image
CN119131298B (en) * 2024-09-09 2025-06-24 中国科学技术大学 Three-dimensional reconstruction method of cardiac chambers based on echocardiography and related devices

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111166388A (en) * 2020-02-28 2020-05-19 合肥凯碧尔高新技术有限公司 Method and device for building three-dimensional model based on two-dimensional ultrasonic image cognition
CN115619943A (en) * 2022-11-02 2023-01-17 蓝影(重庆)数字科技有限公司 Three-dimensional reconstruction method and system for surgical images

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110136264B (en) * 2019-05-30 2021-02-19 北京中盛博方环保工程技术有限公司 Three-dimensional laser scanning-based stock ground material modeling method and system
CN117078695B (en) * 2023-08-18 2024-09-03 内蒙古工业大学 A method for carotid artery plaque ultrasound image recognition and segmentation based on deep learning
CN117409151B (en) * 2023-12-14 2024-06-28 江苏霆升科技有限公司 Three-dimensional surface reconstruction method and device for two-dimensional intracardiac ultrasonic catheter image

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111166388A (en) * 2020-02-28 2020-05-19 合肥凯碧尔高新技术有限公司 Method and device for building three-dimensional model based on two-dimensional ultrasonic image cognition
CN115619943A (en) * 2022-11-02 2023-01-17 蓝影(重庆)数字科技有限公司 Three-dimensional reconstruction method and system for surgical images

Also Published As

Publication number Publication date
WO2025123781A1 (en) 2025-06-19
CN117409151A (en) 2024-01-16

Similar Documents

Publication Publication Date Title
JP7407790B2 (en) Ultrasound system with artificial neural network for guided liver imaging
CN117409151B (en) Three-dimensional surface reconstruction method and device for two-dimensional intracardiac ultrasonic catheter image
CN110807829B (en) Method for constructing three-dimensional heart model based on ultrasonic imaging
US8483488B2 (en) Method and system for stabilizing a series of intravascular ultrasound images and extracting vessel lumen from the images
EP1690230B1 (en) Automatic multi-dimensional intravascular ultrasound image segmentation method
US20120065499A1 (en) Medical image diagnosis device and region-of-interest setting method therefore
CN101849813A (en) 3D Cardiac Ultrasound Virtual Endoscopy System
JP2011131062A (en) Fast anatomical mapping using ultrasound image
WO2022213654A1 (en) Ultrasonic image segmentation method and apparatus, terminal device, and storage medium
CN112641464B (en) Method and system for enabling context-aware ultrasound scanning
JP2002224116A (en) Ultrasound diagnostic device and image processing device
KR20110128197A (en) Automatic Analysis of Cardiac M-Mode Views
CN110313941B (en) Data processing method, device, equipment and storage medium
CN115861132B (en) Blood vessel image correction method, device, medium and equipment
WO2020217860A1 (en) Diagnostic assistance device and diagnostic assistance method
CN116486019B (en) A method and system for three-dimensional cardiac modeling based on three-dimensional cardiac mapping
US20220273267A1 (en) Ultrasonic imaging method and ultrasonic imaging system
CN117115355A (en) Three-dimensional ultrasonic modeling method, system, electronic device and readable storage medium
US20230133103A1 (en) Learning model generation method, image processing apparatus, program, and training data generation method
CN114098795B (en) System and method for generating ultrasound probe guidance instructions
JP2022179433A (en) Image processing device and image processing method
CN111476764B (en) A method for three-dimensional reconstruction of motion-blurred CT images
WO2021099171A1 (en) Systems and methods for imaging screening
US12283048B2 (en) Diagnosis support device, diagnosis support system, and diagnosis support method
CN116363038A (en) Ultrasonic image fusion method, device, computer equipment and storage medium

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
CB02 Change of applicant information
CB02 Change of applicant information

Country or region after: China

Address after: 21st Floor, Block B, Phase I, Zhongdan Ecological Life Science Industrial Park, No. 3-1, Xinjinhu Road, Jiangbei New District, Nanjing City, Jiangsu Province 211899

Applicant after: JIANGSU TINGSN TECHNOLOGY Co.,Ltd.

Address before: Room 2105-2111, block B, phase I, Zhongdan ecological life science and Technology Industrial Park, No. 3-1, xinjinhu Road, Jiangbei new area, Nanjing, Jiangsu 211899

Applicant before: JIANGSU TINGSN TECHNOLOGY Co.,Ltd.

Country or region before: China

GR01 Patent grant
GR01 Patent grant