WO2023142956A1 - Total hip replacement preoperative planning system based on deep learning - Google Patents

Total hip replacement preoperative planning system based on deep learning Download PDF

Info

Publication number
WO2023142956A1
WO2023142956A1 PCT/CN2023/070788 CN2023070788W WO2023142956A1 WO 2023142956 A1 WO2023142956 A1 WO 2023142956A1 CN 2023070788 W CN2023070788 W CN 2023070788W WO 2023142956 A1 WO2023142956 A1 WO 2023142956A1
Authority
WO
WIPO (PCT)
Prior art keywords
dimensional
hip joint
total hip
femoral
image
Prior art date
Application number
PCT/CN2023/070788
Other languages
French (fr)
Chinese (zh)
Inventor
张逸凌
刘星宇
Original Assignee
北京长木谷医疗科技有限公司
张逸凌
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 北京长木谷医疗科技有限公司, 张逸凌 filed Critical 北京长木谷医疗科技有限公司
Publication of WO2023142956A1 publication Critical patent/WO2023142956A1/en

Links

Images

Classifications

    • 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/045Combinations of networks
    • 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
    • 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
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T19/00Manipulating 3D models or images for computer graphics
    • G06T19/20Editing of 3D images, e.g. changing shapes or colours, aligning objects or positioning parts
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2219/00Indexing scheme for manipulating 3D models or images for computer graphics
    • G06T2219/20Indexing scheme for editing of 3D models
    • G06T2219/2016Rotation, translation, scaling

Definitions

  • This application relates to the technical field of artificial intelligence, in particular to a preoperative planning system for total hip replacement based on deep learning.
  • Joint replacement refers to the use of materials such as metal, polymer polyethylene or ceramics to make artificial joint prostheses according to the shape, structure and function of human joints, and implant them into the human body through surgical techniques.
  • the present application provides a preoperative planning system for total hip replacement based on deep learning.
  • the present application provides a preoperative planning system for total hip joint replacement based on deep learning, including: a total hip joint image acquisition module, which is used to acquire a three-dimensional block diagram of the total hip joint to be identified, and the three-dimensional block image to be identified
  • the picture is formed by stacking multiple two-dimensional cross-sectional images of the total hip joint;
  • the total hip joint identification module is used to input the three-dimensional block diagram of the total hip joint to be identified into the trained three-dimensional segmentation neural network, Obtaining the femoral region in each two-dimensional cross-sectional image of the total hip joint output by the trained three-dimensional segmentation neural network, the trained three-dimensional segmentation neural network of the total hip joint is preset by a label marked with the femoral region
  • the 3D block diagram of the total hip joint is used as a training sample, which is obtained by training the convolutional neural network;
  • the 3D image construction module of the total hip joint is used for 3D reconstruction based on the femur in each 2D cross-
  • the total hip recognition module is also used to: acquire several two-dimensional cross-sectional images of sample hip joints; Annotate the sample femoral region in the two-dimensional cross-sectional image of the joint, and mark the femoral head region label on the femoral head pixels in the sample femoral region to obtain several first preset two-dimensional cross-sectional images of the hip joint;
  • the acquisition sequence of the two-dimensional cross-section is to stack several first preset two-dimensional cross-sectional images of the hip joint to obtain the corresponding first preset three-dimensional block image;
  • the shape image is input to the initial three-dimensional segmentation neural network for training to obtain a trained three-dimensional segmentation neural network; wherein, the initial three-dimensional segmentation neural network is constructed by the U-Net convolutional network, and the U The convolution kernel of the -Net convolutional network is a three-dimensional convolution kernel.
  • the total hip joint recognition module when inputting a plurality of the first preset three-dimensional block images into the initial The three-dimensional segmentation neural network is trained, and after the trained three-dimensional segmentation neural network is obtained, it is also used to: mark the cortical bone region labels on the cortical bone pixels of the sample femur region in the first preset two-dimensional cross-sectional images of the hip joint , obtaining several second preset two-dimensional cross-sectional images of the hip joint; stacking the several second preset two-dimensional cross-sectional images of the hip joint in order to obtain a corresponding second preset three-dimensional block diagram; A plurality of the second preset three-dimensional block diagrams are used to optimize the parameters of the trained three-dimensional segmentation neural network to obtain a hip joint recognition model.
  • the total hip recognition module after acquiring the three-dimensional block diagram of the total hip to be recognized, is further used to: if The three-dimensional block diagram of the total hip joint to be identified is a three-dimensional block diagram of the hip joint, and the three-dimensional block diagram of the hip joint is input into the hip joint identification model to obtain the femur in each two-dimensional cross-sectional image of the hip joint.
  • the three-dimensional image construction module of the total hip joint in the three-dimensional reconstruction technology based on the three-dimensional reconstruction technology, according to each two-dimensional cross-sectional image of the total hip joint
  • it is also used to: obtain the pixel coordinates of the femoral head center point of the femoral region in the three-dimensional image based on the centroid formula of the planar image; convert the pixel coordinates into an image coordinates; determine the position of the center of rotation of the femoral head; obtain first size information according to the position of the center of rotation of the femoral head, to determine second size information based on the first size information, wherein the first size information is the size corresponding to the femoral head Information, the second size information is the size information corresponding to the acetabular cup prosthesis model.
  • the preoperative planning system further includes: a position correction module, which is used to place the placed femoral stem prosthesis model The position or the placement position of the acetabular cup prosthesis model is corrected, so that the placement position of the femoral stem prosthesis model and the placement position of the acetabular cup prosthesis model meet the preset position requirements.
  • the three-dimensional image construction module of the total hip joint when acquiring the first size information according to the position of the femoral head rotation center, is specifically used to: determine The intersection point between the rotation center position of the femoral head and the edge of the femoral head area; the first size information is obtained through the length value between the intersection point and the center position of the femoral head.
  • the three-dimensional image construction module of the total hip joint determines the intersection point between the center of rotation of the femoral head and the edge of the femoral head area, and passes the intersection point and the center of the femoral head
  • the length value between the positions, when obtaining the first size information, is specifically used to determine multiple intersections between the center of rotation of the femoral head and the edge of the femoral head area at different angles; through the multiple intersections and the central position of the femoral head The average or median of multiple length values to obtain the first size information.
  • a deep learning-based preoperative planning system for total hip replacement provided by the present application, two continuous three-dimensional convolution kernels in the convolutional neural network are connected through a residual structure, and the convolutional neural network
  • the loss function of is composed of DICE loss and BCE loss.
  • the total hip joint recognition module is specifically used to determine the medullary cavity area according to the femoral area and the cortical bone area. : filter out the cortical bone area in the femoral area to obtain the medullary cavity area.
  • the preoperative planning system for total hip replacement based on deep learning stacks multiple two-dimensional cross-sectional images of the total hip joint into a three-dimensional block diagram of the total hip joint, and is based on a convolutional neural network with a three-dimensional convolution kernel structure. Identify the three-dimensional block diagram of the total hip joint, effectively extract the femoral region of each two-dimensional cross-sectional image of the total hip joint, and perform three-dimensional modeling based on the extracted femoral region to obtain a more accurate three-dimensional model of the femoral region , so as to improve the recognition accuracy of the total hip joint during preoperative planning of total hip joint replacement according to the three-dimensional model of the femoral region.
  • Fig. 1 is a schematic flow chart of the preoperative planning method for total hip arthroplasty based on deep learning provided by the present application;
  • Fig. 2 is a schematic structural diagram of the three-dimensional residual U-Net convolutional neural network provided by the present application
  • Fig. 3 is a schematic diagram of the 3D convolution process provided by the present application.
  • Fig. 4 is the femur region identification rendering based on 3D Res U-Net provided by the application;
  • Fig. 5 is the schematic diagram of a kind of acetabular cup prosthesis model plan provided by the present application.
  • Fig. 6 is a schematic diagram of a femoral stem prosthesis model plan provided by the present application.
  • Fig. 7 is the effect diagram of the bone cortex region recognition based on 3D Res U-Net provided by the present application.
  • FIG. 8 is a schematic flowchart of another deep learning-based preoperative planning method for total hip arthroplasty provided by the present application.
  • FIG. 9 is a schematic structural diagram of a preoperative planning system for total hip arthroplasty based on deep learning provided by the present application.
  • FIG. 10 is a schematic structural diagram of an electronic device provided by the present application.
  • Figure 1 is a schematic flowchart of the preoperative planning method for total hip arthroplasty based on deep learning provided by this application. As shown in Figure 1, this application provides a preoperative planning method for total hip arthroplasty based on deep learning, including:
  • Step 101 acquiring a three-dimensional block diagram of the total hip joint to be identified, where the three-dimensional block diagram of the total hip joint to be identified is formed by stacking multiple two-dimensional cross-sectional images of the total hip joint.
  • the medical image dataset of the total hip joint to be identified is constructed.
  • the image is a two-dimensional cross-sectional image of the total hip joint. Different two-dimensional cross-sectional images of the total hip joint correspond to different bone regions.
  • the bone area corresponding to the image is the tibia area (including a part of the femur area).
  • Image of the right femur then, convert the 2D cross-sectional DICOM data of the hip joint into a JPG image, and convert multiple converted hip joint
  • the cross-sectional images are stacked to form a three-dimensional block map of the total hip joint to be identified, and the three-dimensional block map of the hip joint is obtained.
  • Step 102 input the three-dimensional block diagram of the total hip joint to be identified into the trained three-dimensional segmentation neural network, and obtain each two-dimensional cross-sectional image of the total hip joint output by the trained three-dimensional segmentation neural network
  • the femoral region, the trained three-dimensional segmentation neural network is obtained by training the convolutional neural network from a preset three-dimensional block diagram marked with a label of the femoral region.
  • the stacked 3D block diagram of the hip joint is recognized through the trained 3D segmentation neural network.
  • the convolution kernel of the convolutional neural network used to construct the three-dimensional segmentation neural network is a three-dimensional convolution kernel, so that the three-dimensional segmentation neural network obtained based on the convolutional neural network training can be used for image segmentation.
  • Feature extraction is carried out in three dimensions of the three-dimensional block map of the hip joint, effectively extracting the feature information of the femoral region between each cross-section, and finally segmenting the three-dimensional block map corresponding to the femoral region.
  • the sample femur region of the three-dimensional block image of the sample hip joint has been marked, so that the three-dimensional segmentation neural network can directly identify the femur region on the three-dimensional hip joint image.
  • the pixels of the femoral head in the sample femoral area in the three-dimensional block diagram of the sample hip joint are also marked, so that the model can identify the entire femoral area while also being able to The femoral head region is further identified, so that the femoral region in the hip joint can be quickly identified, saving time and cost, and improving the recognition accuracy of the hip joint.
  • the knee joint in addition to obtaining the two-dimensional cross-sectional image of the hip joint provided in the above embodiment, it also includes a two-dimensional cross-sectional image of the knee joint (which can be reconstructed through three-dimensional, Construct the 3D images of the tibia region and the patella region), and then segment the femur region through the trained 3D segmentation neural network, and then obtain the end of the femur region, and the 3D images of the tibia and patella region according to the segmentation recognition, and determine the relationship between the knee joint and the tibia Prosthesis type and installation position.
  • a two-dimensional cross-sectional image of the knee joint which can be reconstructed through three-dimensional, Construct the 3D images of the tibia region and the patella region
  • segment the femur region through the trained 3D segmentation neural network
  • Step 103 based on the three-dimensional reconstruction technology, according to the femoral region in each two-dimensional cross-sectional image of the total hip joint, obtain the three-dimensional image of the femoral region.
  • the identified femoral region is formed by stacking multiple two-dimensional cross-sectional images, and then three-dimensional information is reconstructed by using key point information in multiple two-dimensional images through three-dimensional reconstruction technology, thereby obtaining a three-dimensional image of the femur.
  • the preoperative planning method for total hip arthroplasty based on deep learning stacks multiple two-dimensional cross-sectional images of the total hip joint into a three-dimensional block diagram of the total hip joint, and based on the convolutional neural network of the three-dimensional convolution kernel structure, Identify the three-dimensional block diagram of the total hip joint, effectively extract the femoral region of each two-dimensional cross-sectional image of the total hip joint, and perform three-dimensional modeling based on the extracted femoral region to obtain a more accurate three-dimensional model of the femoral region , so as to improve the recognition result of the total hip joint during the preoperative planning of the total hip joint replacement according to the three-dimensional model of the femoral region.
  • the trained three-dimensional segmentation neural network is obtained through the following steps:
  • the preset three-dimensional block image set is used as a training sample, input to the initial three-dimensional segmentation neural network for training, and a trained three-dimensional segmentation neural network is obtained;
  • the set of preset three-dimensional block images includes several preset three-dimensional block images obtained by stacking preset two-dimensional cross-sectional images.
  • the preset three-dimensional block diagram is constructed according to the actual bone region segmentation requirements.
  • the preset 3D block diagram is formed by stacking the 2D cross-sectional images of the hip joint; if the preoperative planning of the knee joint and tibial prosthesis is required, the preset In addition to the 2D cross-sectional image of the hip joint, the 3D block diagram also includes the corresponding 2D cross-sectional images of the patella region and the tibial region.
  • the trained 3D segmentation neural network After identifying the femur region through the trained 3D segmentation neural network, it is combined with the existing 3D reconstruction Technology, based on the two-dimensional cross-sectional images of the patella region and the tibia region, constructs a three-dimensional image of the patella and tibia, thereby providing more accurate image data for the preoperative planning of the knee joint and tibial prosthesis.
  • the training process specifically includes:
  • Step 201 acquiring several two-dimensional cross-sectional images of the sample hip joint
  • Step 202 labeling the sample femoral region in each sample hip joint two-dimensional cross-sectional image, and marking the femoral head region label on the femoral head pixels of the sample femoral region, and obtaining several first preset two-dimensional hip joint images. cross-sectional images;
  • Step 203 According to the acquisition sequence of the sample hip joint two-dimensional cross-section, stack several first preset two-dimensional cross-sectional images of the hip joint to obtain a corresponding first preset three-dimensional block image;
  • Step 204 inputting a plurality of the first preset three-dimensional block images to the initial three-dimensional segmentation neural network for training to obtain a trained three-dimensional segmentation neural network;
  • the initial three-dimensional segmentation neural network is constructed by the U-Net convolution network
  • the convolution kernel of the U-Net convolution network is a three-dimensional convolution kernel
  • the hip joint medical image sample data set is obtained, the sample hip joint two-dimensional cross-sectional image in the sample set is manually labeled with the sample femur region, and only the label containing the femur part is extracted as a mask mask. Since this application can also extract the femoral head part of the femoral region while extracting the femoral region, the pixels of the femoral head in the sample femoral region are also marked during labeling, and finally the marked sample hip joint two
  • the three-dimensional cross-sectional images that is, the first preset two-dimensional cross-sectional image of the hip joint
  • this application converts the DICOM data of the first preset two-dimensional cross-section of the hip joint into a picture in JPG format, and at the same time, converts the annotation file into a picture in PNG format.
  • the ratio is divided into training set, validation set and test set. Since the input image of the convolutional neural network of the present application has one more dimension than the input image of the existing 2D network, that is, the three-dimensional block diagram of the sample hip joint (the first preset three-dimensional block image) is formed by stacking multiple two-dimensional cross-sectional images Therefore, the annotation file corresponds to the 3D block diagram of the sample hip joint, which is also a block diagram.
  • the first preset three-dimensional block image is obtained by pre-extracting and stacking two-dimensional cross-sections of the femoral part, and the femoral head region is calculated on the three-dimensional block image of the femoral part stack. Labeling enables the segmentation network to more quickly segment the femoral head region in the image from the background after the training is completed.
  • Fig. 2 is a schematic structural diagram of the three-dimensional residual U-Net convolutional neural network provided by the application, as shown in Fig. 2, in this application, the three-dimensional residual U-Net convolutional neural network (referred to as 3D Res UNet ) is built based on U-Net, and also includes an Encoder (encoder) part and a Decoder (decoder) part, where the Encoder part is used to analyze the entire picture and perform feature extraction and analysis, and with it The corresponding Decoder part is the process of restoring features.
  • the obtained 3D segmentation neural network can segment the femur region in the 3D block image to be recognized, and obtain the segmented femoral region block state diagram.
  • the Encoder part is composed of the basic residual module ResBlock and the maximum pooling layer MaxPooling. Its input size (Input shape) is 1*8*256*256*1.
  • ResBlock is composed of basic convolutional blocks, including two A continuous 3D convolution kernel, each 3D convolution kernel includes two sets of operations: 2*(conv+relu+bn), that is, 3D convolution (conv), activation function (relu) and batch normalization (Batch Normalization , referred to as BN) constitute a set of operations, among which, 3D convolution can effectively extract the information between the cross-sections and reduce the false detection rate; the activation function can increase the nonlinear ability of the model and improve the feature extraction ability of the model; batch normalization
  • the optimization operation can change the distribution of data, which is conducive to the rapid convergence of network training.
  • Figure 3 is a schematic diagram of the 3D convolution process provided by this application, as shown in Figure 3, since the convolution kernel is three-dimensional, it can slide in a certain step in three dimensions, and it will surround 3*3*3 The information of the region is combined to form a point, which can extract richer feature information.
  • two consecutive three-dimensional convolution kernels in the convolutional neural network are connected through a residual structure.
  • two consecutive 3D convolution kernels can prevent network degradation through a residual connection (Skip Connection).
  • the Maxpooling operation uses the maximum value of the adjacent fixed-size region as the characteristic representation of the region, which can effectively reduce the calculation of network parameters.
  • the Decoder part is composed of upsampling and convolution blocks (deconvolution Deconv).
  • the Encoder feature map and the upsampled feature map will be channel-stacked (Concat), and then deconvolved through the convolution block of the Decoder part, and its output shape (output shape) is 1*8* 256*256*1, where the composition of the convolution block is the same as the Encoder process, and will not be repeated here.
  • the optimizer uses an Adam optimizer, and the loss function of the convolutional neural network is composed of DICE loss and BCE loss (DICE loss and BCE loss are respectively a type of loss function). Since the loss function used is a fusion of DICE loss and BCE loss, this can avoid oscillations during network training when only DICE loss is used.
  • the DICOM data and annotation files of the two-dimensional cross-section of the hip joint of the entire sample case were converted into images in JPG and PNG formats in sequence, and packaged into image blocks (ie, block diagrams), which passed the test. Get test DICE.
  • the method further includes:
  • the first size information is obtained to determine the second size information according to the first size, wherein the first size information is the size information corresponding to the femoral head, and the second size information is the hip
  • the size information corresponding to the acetabular cup prosthesis model is obtained.
  • the size information includes at least a diameter and a radius, for example, by obtaining the diameter of the femoral head, the diameter of the acetabular cup prosthesis model is determined.
  • 3D Res U-Net is used to identify each pixel area of the three-dimensional block map of the hip joint.
  • the pixel label can be divided into two attribute values, named 0 and 1 respectively, where the value 0 Represents the background pixel, 1 represents the femoral head pixel; after the labeling is completed, the marked image data is passed into the convolutional neural network (ie 3D Res UNet), and the convolution pooling sampling is used to iteratively learn and train.
  • Fig. 4 is the effect diagram of the femoral region recognition based on 3D Res U-Net provided by the present application. As shown in Fig.
  • the trained three-dimensional segmentation neural network can automatically recognize the position of the femoral head and complete the femoral head region identification Identification (marked schematically with a white line in Figure 4). Finally, through 3D reconstruction, a 3D image of the femoral head region is obtained.
  • the centroid and centroid coincide.
  • the centroid formula of the planar image the coordinates of the center point of the femoral head in the 3D image of the femoral region can be obtained, that is, the femoral head rotation in the 3D image center.
  • B[i, j] represents the pixel value of the i-th row and j-th column pixel in the binary image B, so the following formula can be used to obtain the femoral head of the femoral region in the three-dimensional image
  • the location of the center point is
  • A represents the sum of the pixel values of all pixels in the binary image
  • n represents the maximum number of rows of pixels in the binary image
  • m represents the maximum number of columns of pixels in the binary image
  • the center coordinates of the image plane coordinates are:
  • S x , S y are the row and column spacing of the image array respectively.
  • the diameter of the femoral head is based on the center of the femoral head, intersects a ray with the edge of the femoral head region obtained based on the preoperative planning method for total hip arthroplasty provided by this application, calculates the length from the intersection point to the center of the femoral head, and then calculates the length from the intersection point to the center of the femoral head
  • One degree is the step length, and it is calculated once per one degree of rotation.
  • the method After inputting a plurality of the first preset three-dimensional block images into the initial three-dimensional segmentation neural network for training and obtaining a trained three-dimensional segmentation neural network, the method also includes:
  • the trained three-dimensional segmentation neural network is fine-tuned to obtain a hip joint recognition model.
  • the parameters of the three-dimensional segmentation neural network are optimized according to the obtained new training set, so that the model recognizes the femoral head region while , but also to identify cortical areas.
  • the method further includes:
  • FIG. 6 provides a schematic diagram of a femoral stem prosthesis model plan.
  • the three-dimensional segmentation neural network after parameter optimization that is, the hip joint recognition model
  • the pixel label is divided into three attribute values, named 0, 1 and 2 respectively, where the value 0 represents background pixels, 1 represents femoral head pixels, and 2 represents cortical bone.
  • Figure 7 is the effect diagram of the recognition of the cortical bone region based on 3D Res U-Net provided by the present application. As shown in Figure 7, after optimizing the parameters of the three-dimensional segmentation neural network, the femoral head in the block diagram can be simultaneously identified. and identification of cortical bone (schematically marked with a black line in Figure 7).
  • the present application intercepts images from the end of the lesser trochanter to the end of the femur from the recognition results output by the hip joint recognition model, and subtracts the cortical bone region from the femoral region in the image to obtain the medullary cavity region.
  • each row i.e., the level of the medullary canal is identified from the two-dimensional cross-section of each hip joint, and each row refers to the two-dimensional cross-sectional image of the hip joint segmented by this application, after three-dimensional Reconstruct the simulated X-ray projection effect map, and then start below the end position of the lesser trochanter, draw a horizontal line on the image at certain preset positions, and when the horizontal line intersects with the edges of the two femoral medullary canals, four points will be obtained) and
  • the intersection of the medullary cavity is four coordinates, which are named A1, A2, B1, and B2 from left to right; the midpoint can be obtained based on two points, namely A1 (X 1 , Y 1 ), A2 (X 2 , Y 2 )
  • the midpoint coordinates of can be obtained by the following formula:
  • B1 and B2 can be calculated in the same way.
  • the coordinates of the midpoint of the medullary cavity are calculated in turn for each row, and these points are fitted into a straight line to determine the anatomical axis of the medullary cavity.
  • the angle value of the femoral neck-shaft angle was calculated, combined with the shape of the medullary cavity and the position of the center of rotation of the femoral head, the model and placement position of the femoral stem prosthesis model could be determined together.
  • the femoral neck-shaft angle is the angle between the anatomical axis of the medullary canal and the femoral neck axis.
  • the preoperative planning plan for the total hip joint can be output.
  • the preset position requirements of the acetabular cup prosthesis can be, for example, that after the acetabular cup prosthesis is placed in the acetabular socket, the coverage rate covering the acetabular socket is greater than 75%, and the preset position requirements of the femoral stem prosthesis can be, for example, After the femoral stem prosthesis is placed in the medullary cavity, the angle between the long axis of the femoral stem prosthesis and the long axis of the femur is less than or equal to 3°.
  • the preoperative planning method for total hip replacement based on deep learning includes:
  • Step 801 Obtain the three-dimensional block diagram of the total hip joint to be identified.
  • Step 802 Input the 3D block image of the total hip joint to be identified into the trained 3D segmentation neural network to obtain the femur region in each 2D cross-sectional image of the total hip joint.
  • Step 803 Based on the three-dimensional reconstruction technology, according to the femoral region in each two-dimensional cross-sectional image of the total hip joint, a three-dimensional image of the femoral region is obtained.
  • Step 804 Determine the size information of the acetabular cup prosthesis, as well as the type and placement position of the femoral stem prosthesis.
  • Step 805 Put the determined posterior acetabular cup prosthesis and femoral stem prosthesis into the acetabular fossa and medullary cavity respectively.
  • Step 806 Correcting the position or type of the inserted acetabular cup prosthesis and the position or type of the femoral stem prosthesis.
  • Step 807 Output the preoperative planning scheme for total hip replacement.
  • the preoperative planning system for total hip replacement based on deep learning provided by this application.
  • the preoperative planning system for total hip replacement based on deep learning described below is the same as the preoperative planning system for total hip replacement based on deep learning described above Planning methods can be cross-referenced.
  • Fig. 9 is a schematic structural diagram of the preoperative planning system for total hip replacement based on deep learning provided by the present application.
  • the present application provides a preoperative planning system for total hip replacement based on deep learning, including A hip joint image collection module 901, a total hip joint recognition module 902 and a total hip joint three-dimensional image construction module 903, wherein the total hip joint image collection module 901 is used to obtain a three-dimensional block diagram of the total hip joint to be identified, and the total hip joint
  • the three-dimensional block diagram of the joint is formed by stacking multiple two-dimensional cross-sectional images of the total hip joint; the total hip joint identification module 902 is used to input the three-dimensional block diagram of the total hip joint to be identified into the trained three-dimensional segmentation
  • the femur region in each two-dimensional cross-sectional image of the total hip joint output by the trained three-dimensional segmentation neural network is obtained, and the trained three-dimensional segmentation neural network is a pre-prepared image of a label marked with the f
  • a two-dimensional cross-sectional image of the hip joint is used for illustration.
  • the two-dimensional cross-sectional image is DICOM data
  • each two-dimensional cross-sectional image of the hip joint includes an image of the pelvic region, an image of the left femur, and an image of the right femur. image.
  • the total hip joint image acquisition module 901 converts the DICOM data of the two-dimensional cross-section of the hip joint into a picture in JPG format, and according to the acquisition order of the cross-section (for example, from the femoral head to the end of the femur), multiple converted hip joints
  • the two-dimensional cross-sectional images of the joints are stacked to generate the three-dimensional block diagram of the total hip joint to be identified, that is, the three-dimensional block diagram of the hip joint is obtained.
  • the trained three-dimensional segmentation neural network is configured in the total hip joint identification module 902, which can identify the stacked three-dimensional block diagrams of the hip joint.
  • the convolution kernel of the convolutional neural network used to construct the three-dimensional segmentation neural network is a three-dimensional convolution kernel, so that the three-dimensional segmentation neural network obtained based on the convolutional neural network training can be used for image segmentation.
  • Feature extraction is carried out in three dimensions of the three-dimensional block map of the hip joint, effectively extracting the feature information of the femoral region between each cross-section, and finally segmenting the three-dimensional block map corresponding to the femoral region.
  • two continuous three-dimensional convolution kernels in the convolutional neural network are connected through a residual structure; and the loss function of the convolutional neural network is composed of DICE loss and BCE loss.
  • the femur region of the sample hip joint three-dimensional block diagram has been marked, so that the three-dimensional segmentation neural network in the total hip joint recognition module 902 can directly perform femur region analysis on the three-dimensional hip joint image. recognition.
  • the pixels of the femoral head in the sample femoral area in the three-dimensional block diagram of the sample hip joint are also marked, so that the total hip joint identification module 802 can identify the entire femoral area At the same time, it can further identify the femoral head region, so that the femoral region in the hip joint can be quickly identified, saving time and cost, and improving the recognition accuracy of the hip joint.
  • the key point information of the identified femoral region is used to perform three-dimensional reconstruction to obtain a three-dimensional image of the femur, so that in the subsequent preoperative planning of total hip joint replacement, the three-dimensional image of the femur can be Determine the size, type and placement of the prosthesis relative to the total hip.
  • the preoperative planning system for total hip replacement based on deep learning stacks multiple two-dimensional cross-sectional images of the total hip joint into a three-dimensional block diagram of the total hip joint, and is based on a convolutional neural network with a three-dimensional convolution kernel structure. Identify the three-dimensional block diagram of the total hip joint, effectively extract the femoral region of each two-dimensional cross-sectional image of the total hip joint, and perform three-dimensional modeling according to the extracted femoral region, so that according to the three-dimensional model of the femoral region, in the Improve the accuracy of total hip joint recognition when performing preoperative planning for total hip arthroplasty.
  • the system also includes a training module, which is used to input the preset three-dimensional block image set as a training sample into the initial three-dimensional segmentation neural network for training, and obtain a trained three-dimensional segmentation neural network;
  • the set of preset three-dimensional block images includes several preset three-dimensional block images obtained by stacking preset two-dimensional cross-sectional images.
  • the training module includes a sample two-dimensional cross-sectional image acquisition unit, a first labeling unit, a block diagram first construction unit and a first training unit, wherein the sample two-dimensional cross-sectional image acquisition unit It is used to obtain several two-dimensional cross-sectional images of the sample hip joint;
  • the first labeling unit is used to mark the sample femoral area in each sample hip joint two-dimensional cross-sectional image, and to mark the femoral head of the sample femoral area
  • the pixels are labeled with the femoral head area label to obtain several first preset two-dimensional cross-sectional images of the hip joint;
  • the preset two-dimensional cross-sectional images of the hip joint are stacked to obtain a corresponding first preset three-dimensional block image;
  • the first training unit is configured to input a plurality of the first preset three-dimensional block images into the initial three-dimensional Segment the neural network for training to obtain a trained three-dimensional segmentation neural network;
  • the initial three-dimensional segmentation neural network is constructed by the U-Net convolution network
  • the convolution kernel of the U-Net convolution network is a three-dimensional convolution kernel
  • the system further includes a second labeling unit, a second construction unit of the block diagram, and a second training unit, wherein the second labeling unit is used to identify several sheets of the first preset hip joint
  • the cortical bone pixels of the sample femoral region are labeled with the cortical bone region label to obtain several second preset two-dimensional cross-sectional images of the hip joint
  • the preset two-dimensional cross-sectional images of the hip joint are stacked in order to obtain the corresponding second preset three-dimensional block diagram
  • the second training unit is used to train the training
  • the parameters of a good three-dimensional segmentation neural network are optimized to obtain a hip joint recognition model.
  • the system also includes a cortical bone area identification module, a medullary cavity area determination module, a medullary cavity anatomical axis determination module, a femoral neck-shaft angle calculation module and a femoral stem prosthesis determination module, wherein the cortical bone
  • the area recognition module is used to input the three-dimensional block diagram of the hip joint into the hip joint recognition model if the three-dimensional block diagram of the total hip joint to be identified is a three-dimensional block diagram of the hip joint, and obtain two dimensions of each hip joint.
  • the medullary cavity region determination module is used to determine the medullary cavity region according to the femoral region and the cortical bone region;
  • the medullary cavity anatomical axis determination module is used to calculate the medullary cavity
  • the midpoint coordinates of each medullary cavity level in the region, and according to the midpoint coordinates, all the central points are fitted with a straight line to determine the anatomical axis of the medullary cavity;
  • the femoral neck shaft angle calculation module is used to calculate the anatomical axis of the medullary cavity and femoral neck axis to calculate the angle value of the femoral neck shaft angle;
  • the femoral stem prosthesis determination module is used to determine the type of femoral stem prosthesis model and Placement.
  • the system also includes a femoral head center point pixel coordinate calculation module, a coordinate conversion module, a femoral head rotation center determination module, and an acetabular cup prosthesis determination module, wherein the bone center point pixel coordinate calculation module It is used to obtain the pixel coordinates of the center point of the femoral head in the femoral region in the three-dimensional image based on the centroid formula of the planar image; the coordinate transformation module is used to convert the pixel coordinates into image coordinates; the femoral head rotation center determination module is used to determine The position of the center of rotation of the femoral head; the acetabular cup prosthesis determination module is used to obtain the first size information according to the position of the center of rotation of the femoral head, so as to determine the second size information according to the first size information, wherein the first size information is the size information corresponding to the femoral head, and the second size information is the size information corresponding to the acetabular cup prosthesis
  • the preoperative planning system also includes: a position correction module, which is used to adjust the placement position of the placed femoral stem prosthesis model or the placement position of the acetabular cup prosthesis model Correction, so that the placement position of the femoral stem prosthesis model and the placement position of the acetabular cup prosthesis model meet the preset position requirements.
  • a position correction module which is used to adjust the placement position of the placed femoral stem prosthesis model or the placement position of the acetabular cup prosthesis model Correction, so that the placement position of the femoral stem prosthesis model and the placement position of the acetabular cup prosthesis model meet the preset position requirements.
  • the three-dimensional image construction module of the total hip joint is specifically used for:
  • the first size information is obtained through the length value between the intersection point and the center position of the femoral head.
  • the three-dimensional image construction module of the total hip joint after determining the intersection point between the center of rotation of the femoral head and the edge of the femoral head region, obtains the first size information through the length value between the intersection point and the center position of the femoral head when, specifically for;
  • the first dimension information is obtained by using the average value or the median value of the multiple length values between the multiple intersection points and the central position of the femoral head.
  • the total hip joint recognition module is specifically used for:
  • the cortical bone area in the femoral area was filtered out to obtain the medullary cavity area.
  • Fig. 10 is a schematic structural diagram of an electronic device provided by the present application.
  • the electronic device may include: a processor (Processor) 1001, a communication interface (Communications Interface) 1002, a memory (Memory) 1003 and a communication bus 1004, Wherein, the processor 1001 , the communication interface 1002 , and the memory 1003 communicate with each other through the communication bus 1004 .
  • the processor 1001 can call the logic instructions in the memory 1003 to execute a deep learning-based image planning method for preoperative total hip arthroplasty.
  • the processor 1001 is configured to obtain a three-dimensional block diagram of the total hip joint to be identified, and the three-dimensional block diagram of the total hip joint to be identified is obtained by stacking multiple two-dimensional cross-sectional images of the total hip joint
  • the three-dimensional block diagram of the total hip joint to be identified is input into the trained three-dimensional segmentation neural network, and each two-dimensional cross-sectional image of the total hip joint output by the trained three-dimensional segmentation neural network is obtained.
  • the trained three-dimensional segmentation neural network is obtained by training the convolutional neural network using the preset three-dimensional block diagram marked with the label of the femoral region as a training sample; based on the three-dimensional reconstruction technology, according to each The femoral region in the two-dimensional cross-sectional image of the total hip joint is obtained to obtain a three-dimensional image of the femoral region.
  • the above logic instructions in the memory 1003 may be implemented in the form of software function units and when sold or used as an independent product, may be stored in a computer-readable storage medium.
  • the technical solution of the present application is essentially or the part that contributes to the prior art or the part of the technical solution can be embodied in the form of a software product, and the computer software product is stored in a storage medium, including Several instructions are used to make a computer device (which may be a personal computer, a server, or a network device, etc.) execute all or part of the steps of the methods described in the various embodiments of the present application.
  • the aforementioned storage media include: U disk, mobile hard disk, read-only memory (ROM, Read-Only Memory), random access memory (RAM, Random Access Memory), magnetic disk or optical disc, etc., which can store program codes. .
  • the present application also provides a computer program product
  • the computer program product includes a computer program stored on a non-transitory computer-readable storage medium
  • the computer program includes program instructions, and when the program instructions are executed by a computer During execution, the computer can execute the deep learning-based preoperative planning method for total hip replacement provided by the above methods.
  • An example is to obtain the three-dimensional block diagram of the total hip joint to be identified, and the three-dimensional block diagram of the total hip joint to be identified is formed by stacking multiple two-dimensional cross-sectional images of the total hip joint;
  • the three-dimensional block diagram of the total hip joint is input into the trained three-dimensional segmentation neural network to obtain the femur region in each two-dimensional cross-sectional image of the total hip joint output by the three-dimensional segmentation neural network trained well.
  • the 3D segmentation neural network is obtained by training the convolutional neural network with the preset 3D block diagram marked with the label of the femoral region as the training sample; based on the 3D reconstruction technology, according to each 2D cross-sectional image of the total hip joint In the femoral region, a three-dimensional image of the femoral region is obtained.
  • the present application also provides a non-transitory computer-readable storage medium, on which a computer program is stored, and when the computer program is executed by a processor, it is implemented to perform the deep learning-based total hip joint provided by the above-mentioned embodiments. Preoperative planning methods for replacement surgery.
  • An example is to obtain the three-dimensional block diagram of the total hip joint to be identified, and the three-dimensional block diagram of the total hip joint to be identified is formed by stacking multiple two-dimensional cross-sectional images of the total hip joint;
  • the three-dimensional block diagram of the total hip joint is input into the trained three-dimensional segmentation neural network to obtain the femur region in each two-dimensional cross-sectional image of the total hip joint output by the three-dimensional segmentation neural network trained well.
  • the 3D segmentation neural network is obtained by training the convolutional neural network with the preset 3D block diagram marked with the label of the femoral region as the training sample; based on the 3D reconstruction technology, according to each 2D cross-sectional image of the total hip joint In the femoral region, a three-dimensional image of the femoral region is obtained.
  • the device embodiments described above are only illustrative, and the units described as separate components may or may not be physically separated, and the components shown as units may or may not be physical units, that is, they may be located in One place, or it can be distributed to multiple network elements. Part or all of the modules can be selected according to actual needs to achieve the purpose of the solution of this embodiment. It can be understood and implemented by those skilled in the art without any creative effort.
  • each implementation can be implemented by means of software plus a necessary general-purpose hardware platform, and of course also by hardware.
  • the essence of the above technical solution or the part that contributes to the prior art can be embodied in the form of software products, and the computer software products can be stored in computer-readable storage media, such as ROM/RAM, magnetic discs, optical discs, etc., including several instructions to make a computer device (which may be a personal computer, server, or network device, etc.) execute the methods described in various embodiments or some parts of the embodiments.

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Software Systems (AREA)
  • General Physics & Mathematics (AREA)
  • General Engineering & Computer Science (AREA)
  • Artificial Intelligence (AREA)
  • Computational Linguistics (AREA)
  • Mathematical Physics (AREA)
  • Computing Systems (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Computer Graphics (AREA)
  • Biomedical Technology (AREA)
  • Biophysics (AREA)
  • Molecular Biology (AREA)
  • Data Mining & Analysis (AREA)
  • Evolutionary Computation (AREA)
  • General Health & Medical Sciences (AREA)
  • Architecture (AREA)
  • Geometry (AREA)
  • Computer Hardware Design (AREA)
  • Image Analysis (AREA)
  • Image Processing (AREA)

Abstract

The present application provides a total hip replacement preoperative planning system based on deep learning, comprising: a total hip joint image acquisition module which is used for obtaining a total hip joint three-dimensional block diagram to be identified; a total hip joint identification module which is used for inputting the total hip joint three-dimensional block diagram to be identified into a trained three-dimensional segmentation neural network to obtain a femur area in each total hip joint two-dimensional cross-sectional image, wherein the trained three-dimensional segmentation neural network is obtained by training a convolutional neural network by means of preset three-dimensional block diagrams labeled with labels of femur areas; and a total hip joint three-dimensional image construction module which is used for, according to the femur area in each total hip joint two-dimensional cross-sectional image and on the basis of three-dimensional reconstruction technology, obtaining a three-dimensional image of the femur area. According to the present application, total hip joint three-dimensional block diagrams are identified, three-dimensional modeling is performed on extracted femur areas on the basis of the three-dimensional block diagrams, and thus, the total hip joint identification precision is improved according to the femur area three-dimensional model.

Description

基于深度学习的全髋关节置换术前规划系统Preoperative planning system for total hip replacement based on deep learning
相关申请的交叉引用Cross References to Related Applications
本申请要求于2022年01月27日提交的申请号为202210101412.9,名称为“基于深度学习的全髋关节置换术前规划系统”的中国专利申请的优先权,其通过引用方式全部并入本文。This application claims the priority of the Chinese patent application with application number 202210101412.9 and titled "Deep Learning-Based Preoperative Planning System for Total Hip Joint Replacement" filed on January 27, 2022, which is incorporated herein by reference in its entirety.
技术领域technical field
本申请涉及人工智能技术领域,尤其涉及一种基于深度学习的全髋关节置换术前规划系统。This application relates to the technical field of artificial intelligence, in particular to a preoperative planning system for total hip replacement based on deep learning.
背景技术Background technique
关节置换术是指采用金属、高分子聚乙烯或陶瓷等材料,根据人体关节的形态、构造及功能制成人工关节假体,通过外科技术植入人体内。Joint replacement refers to the use of materials such as metal, polymer polyethylene or ceramics to make artificial joint prostheses according to the shape, structure and function of human joints, and implant them into the human body through surgical techniques.
在全髋关节置换的术前规划中,医师需要凭借自身的经验,从髋关节医学图像中判断整个股骨区域,确定需要置换的假体尺寸和型号等。但是传统的全髋关节置换技术依赖于医师经验,不同经验的医师会出现不同识别结果,难以保证结果的统一性。近年来,随着医疗水平的提高,利用二维图像分割神经网络识别可以消除这一缺点。In the preoperative planning of total hip replacement, physicians need to rely on their own experience to judge the entire femoral region from medical images of the hip joint, and determine the size and type of the prosthesis that needs to be replaced. However, traditional total hip replacement technology relies on physician experience, and physicians with different experience will have different identification results, making it difficult to guarantee the uniformity of results. In recent years, with the improvement of medical level, the use of two-dimensional image segmentation neural network recognition can eliminate this shortcoming.
由于髋关节形状是三维结构的,使用二维图像分割神经网络进行全髋关节分割,会丢失掉关节连续切片层之间的特征信息,导致全髋关节的识别精度较低。Since the shape of the hip joint is a three-dimensional structure, using a two-dimensional image segmentation neural network for total hip joint segmentation will lose the feature information between the continuous slice layers of the joint, resulting in low recognition accuracy of the total hip joint.
发明内容Contents of the invention
针对现有技术存在的问题,本申请提供一种基于深度学习的全髋关节置换术前规划系统。In view of the problems existing in the prior art, the present application provides a preoperative planning system for total hip replacement based on deep learning.
本申请提供一种基于深度学习的全髋关节置换置换术前规划系统,包括:全髋关节图像采集模块,用于获取待识别的全髋关节三维块状图,所述待识别的三维块状图是由多张全髋关节二维横断面图像堆叠而成的;全髋关节识 别模块,用于将所述待识别的全髋关节三维块状图输入到训练好的三维分割神经网络中,获取所述训练好的三维分割神经网络输出的每张全髋关节二维横断面图像中的股骨区域,所述训练好的全髋关节三维分割神经网络是由标注有股骨区域的标签的预设全髋关节三维块状图作为训练样本,对卷积神经网络进行训练得到的;全髋关节三维图像构建模块,用于基于三维重建技术,根据每张全髋关节二维横断面图像中的股骨区域,得到所述股骨区域的三维图像。The present application provides a preoperative planning system for total hip joint replacement based on deep learning, including: a total hip joint image acquisition module, which is used to acquire a three-dimensional block diagram of the total hip joint to be identified, and the three-dimensional block image to be identified The picture is formed by stacking multiple two-dimensional cross-sectional images of the total hip joint; the total hip joint identification module is used to input the three-dimensional block diagram of the total hip joint to be identified into the trained three-dimensional segmentation neural network, Obtaining the femoral region in each two-dimensional cross-sectional image of the total hip joint output by the trained three-dimensional segmentation neural network, the trained three-dimensional segmentation neural network of the total hip joint is preset by a label marked with the femoral region The 3D block diagram of the total hip joint is used as a training sample, which is obtained by training the convolutional neural network; the 3D image construction module of the total hip joint is used for 3D reconstruction based on the femur in each 2D cross-sectional image of the total hip joint. region to obtain a three-dimensional image of the femoral region.
根据本申请提供的一种基于深度学习的全髋关节置换术前规划系统,所述全髋关节识别模块,还用于:获取若干张样本髋关节二维横断面图像;对每一张样本髋关节二维横断面图像中的样本股骨区域进行标注,并对所述样本股骨区域的股骨头像素标注股骨头区域标签,得到若干张第一预设髋关节二维横断面图像;根据样本髋关节二维横断面的采集顺序,将若干张所述第一预设髋关节二维横断面图像进行堆叠,得到对应的第一预设三维块状图像;将多个所述第一预设三维块状图像输入到所述初始的三维分割神经网络进行训练,得到训练好的三维分割神经网络;其中,所述初始的三维分割神经网络是由U-Net卷积网络构建得到的,且所述U-Net卷积网络的卷积核为三维卷积核。According to a deep learning-based preoperative planning system for total hip replacement provided by the present application, the total hip recognition module is also used to: acquire several two-dimensional cross-sectional images of sample hip joints; Annotate the sample femoral region in the two-dimensional cross-sectional image of the joint, and mark the femoral head region label on the femoral head pixels in the sample femoral region to obtain several first preset two-dimensional cross-sectional images of the hip joint; The acquisition sequence of the two-dimensional cross-section is to stack several first preset two-dimensional cross-sectional images of the hip joint to obtain the corresponding first preset three-dimensional block image; The shape image is input to the initial three-dimensional segmentation neural network for training to obtain a trained three-dimensional segmentation neural network; wherein, the initial three-dimensional segmentation neural network is constructed by the U-Net convolutional network, and the U The convolution kernel of the -Net convolutional network is a three-dimensional convolution kernel.
根据本申请提供的一种基于深度学习的全髋关节置换术前规划系统,所述全髋关节识别模块,在所述将多个所述第一预设三维块状图像输入到所述初始的三维分割神经网络进行训练,得到训练好的三维分割神经网络之后,还用于:对若干张所述第一预设髋关节二维横断面图像中样本股骨区域的骨皮质像素标注骨皮质区域标签,得到若干张第二预设髋关节二维横断面图像;将若干张所述第二预设髋关节二维横断面图像按顺序进行堆叠,得到对应的第二预设三维块状图;通过多个所述第二预设三维块状图,对所述训练好的三维分割神经网络的参数进行优化,得到髋关节识别模型。According to a deep learning-based preoperative planning system for total hip joint replacement provided by the present application, the total hip joint recognition module, when inputting a plurality of the first preset three-dimensional block images into the initial The three-dimensional segmentation neural network is trained, and after the trained three-dimensional segmentation neural network is obtained, it is also used to: mark the cortical bone region labels on the cortical bone pixels of the sample femur region in the first preset two-dimensional cross-sectional images of the hip joint , obtaining several second preset two-dimensional cross-sectional images of the hip joint; stacking the several second preset two-dimensional cross-sectional images of the hip joint in order to obtain a corresponding second preset three-dimensional block diagram; A plurality of the second preset three-dimensional block diagrams are used to optimize the parameters of the trained three-dimensional segmentation neural network to obtain a hip joint recognition model.
根据本申请提供的一种基于深度学习的全髋关节置换术前规划系统,所述全髋关节识别模块,在所述获取待识别的全髋关节三维块状图之后,还用于:若所述待识别的全髋关节三维块状图为髋关节三维块状图,将所述髋关节三维块状图输入到所述髋关节识别模型,得到每张髋关节二维横断面图像中的股骨区域和骨皮质区域;根据所述股骨区域和所述骨皮质区域,确定髓 腔区域;计算所述髓腔区域中每个髓腔层面的中点坐标,并根据所述中点坐标,将所有中心点进行直线拟合,确定髓腔解剖轴线;根据所述髓腔解剖轴线和股骨颈轴线,计算得到股骨颈干角的角度值;根据所述角度值、所述髓腔区域和股骨头旋转中心位置,确定股骨柄假体模型的类型和放置位置。According to a deep learning-based preoperative planning system for total hip replacement provided by the present application, the total hip recognition module, after acquiring the three-dimensional block diagram of the total hip to be recognized, is further used to: if The three-dimensional block diagram of the total hip joint to be identified is a three-dimensional block diagram of the hip joint, and the three-dimensional block diagram of the hip joint is input into the hip joint identification model to obtain the femur in each two-dimensional cross-sectional image of the hip joint. area and cortical bone area; according to the femoral area and the cortical bone area, determine the medullary canal area; calculate the midpoint coordinates of each medullary canal level in the medullary canal area, and according to the midpoint coordinates, all Carry out straight line fitting at the center point to determine the anatomical axis of the medullary cavity; calculate the angle value of the femoral neck-shaft angle according to the anatomical axis of the medullary cavity and the axis of the femoral neck; Center position, to determine the type and placement of the femoral stem prosthesis model.
根据本申请提供的一种基于深度学习的全髋关节置换术前规划系统,所述全髋关节三维图像构建模块,在所述基于三维重建技术,根据每张全髋关节二维横断面图像中的股骨区域,得到所述股骨区域的三维图像之后,还用于:基于平面图像的质心公式,获取所述三维图像中股骨区域的股骨头中心点的像素坐标;将所述像素坐标转换为图像坐标;确定股骨头旋转中心位置;根据股骨头旋转中心位置,获取第一尺寸信息,以根据所述第一尺寸信息确定第二尺寸信息,其中,所述第一尺寸信息为股骨头对应的尺寸信息,所述第二尺寸信息为髋臼杯假体模型对应的尺寸信息。According to a preoperative planning system for total hip joint replacement based on deep learning provided by the present application, the three-dimensional image construction module of the total hip joint, in the three-dimensional reconstruction technology based on the three-dimensional reconstruction technology, according to each two-dimensional cross-sectional image of the total hip joint After obtaining the three-dimensional image of the femoral region, it is also used to: obtain the pixel coordinates of the femoral head center point of the femoral region in the three-dimensional image based on the centroid formula of the planar image; convert the pixel coordinates into an image coordinates; determine the position of the center of rotation of the femoral head; obtain first size information according to the position of the center of rotation of the femoral head, to determine second size information based on the first size information, wherein the first size information is the size corresponding to the femoral head Information, the second size information is the size information corresponding to the acetabular cup prosthesis model.
根据本申请提供的一种基于深度学习的全髋关节置换术前规划系统,所述术前规划系统还包括:位置矫正模块,所述位置矫正模块用于对放置的股骨柄假体模型的放置位置或者髋臼杯假体模型的放置位置进行矫正,以使所述股骨柄假体模型的放置位置以及所述髋臼杯假体模型的放置位置满足预设位置要求。According to a deep learning-based preoperative planning system for total hip arthroplasty provided by the present application, the preoperative planning system further includes: a position correction module, which is used to place the placed femoral stem prosthesis model The position or the placement position of the acetabular cup prosthesis model is corrected, so that the placement position of the femoral stem prosthesis model and the placement position of the acetabular cup prosthesis model meet the preset position requirements.
根据本申请提供的一种基于深度学习的全髋关节置换术前规划系统,所述全髋关节三维图像构建模块,在根据股骨头旋转中心位置,获取第一尺寸信息时,具体用于:确定股骨头旋转中心位置与股骨头区域边缘交点;通过交点与股骨头中心位置之间的长度值,获取第一尺寸信息。According to a preoperative planning system for total hip arthroplasty based on deep learning provided by the present application, the three-dimensional image construction module of the total hip joint, when acquiring the first size information according to the position of the femoral head rotation center, is specifically used to: determine The intersection point between the rotation center position of the femoral head and the edge of the femoral head area; the first size information is obtained through the length value between the intersection point and the center position of the femoral head.
根据本申请提供的一种基于深度学习的全髋关节置换术前规划系统,所述全髋关节三维图像构建模块,在确定股骨头旋转中心位置与股骨头区域边缘交点,通过交点与股骨头中心位置之间的长度值,获取第一尺寸信息时,具体用于;确定股骨头旋转中心位置与股骨头区域边缘在不同角度下的多个交点;通过多个交点与股骨头中心位置之间的多个长度值的平均值或中值,获取第一尺寸信息。According to a preoperative planning system for total hip arthroplasty based on deep learning provided by the present application, the three-dimensional image construction module of the total hip joint determines the intersection point between the center of rotation of the femoral head and the edge of the femoral head area, and passes the intersection point and the center of the femoral head The length value between the positions, when obtaining the first size information, is specifically used to determine multiple intersections between the center of rotation of the femoral head and the edge of the femoral head area at different angles; through the multiple intersections and the central position of the femoral head The average or median of multiple length values to obtain the first size information.
根据本申请提供的一种基于深度学习的全髋关节置换术前规划系统,所述卷积神经网络中两个连续的三维卷积核是通过残差结构连接的,且所述卷积神经网络的损失函数是由DICE loss和BCE loss组成的。According to a deep learning-based preoperative planning system for total hip replacement provided by the present application, two continuous three-dimensional convolution kernels in the convolutional neural network are connected through a residual structure, and the convolutional neural network The loss function of is composed of DICE loss and BCE loss.
根据本申请提供的一种基于深度学习的全髋关节置换术前规划系统,所述全髋关节识别模块,在根据所述股骨区域和所述骨皮质区域,确定髓腔区域时,具体用于:将所述股骨区域中骨皮质区域滤除,得到髓腔区域。According to a preoperative planning system for total hip arthroplasty based on deep learning provided by the present application, the total hip joint recognition module is specifically used to determine the medullary cavity area according to the femoral area and the cortical bone area. : filter out the cortical bone area in the femoral area to obtain the medullary cavity area.
本申请的上述技术方案至少具有如下有益效果:The above technical solution of the present application has at least the following beneficial effects:
本申请提供的基于深度学习的全髋关节置换术前规划系统,通过将多张全髋关节二维横断面图像堆叠成全髋关节三维块状图,基于三维卷积核结构的卷积神经网络,对全髋关节三维块状图进行识别,有效的提取到每张全髋关节二维横断面图像的股骨区域,并根据提取得到的股骨区域进行三维建模,得到更为精准的股骨区域三维模型,以根据该股骨区域三维模型,在进行全髋关节置换术前规划时,提高全髋关节识别精度。The preoperative planning system for total hip replacement based on deep learning provided by this application stacks multiple two-dimensional cross-sectional images of the total hip joint into a three-dimensional block diagram of the total hip joint, and is based on a convolutional neural network with a three-dimensional convolution kernel structure. Identify the three-dimensional block diagram of the total hip joint, effectively extract the femoral region of each two-dimensional cross-sectional image of the total hip joint, and perform three-dimensional modeling based on the extracted femoral region to obtain a more accurate three-dimensional model of the femoral region , so as to improve the recognition accuracy of the total hip joint during preoperative planning of total hip joint replacement according to the three-dimensional model of the femoral region.
附图说明Description of drawings
为了更清楚地说明本申请或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作一简单地介绍,显而易见地,下面描述中的附图是本申请的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。In order to more clearly illustrate the technical solutions in this application or the prior art, the accompanying drawings that need to be used in the description of the embodiments or the prior art will be briefly introduced below. Obviously, the accompanying drawings in the following description are the present For some embodiments of the application, those skilled in the art can also obtain other drawings based on these drawings without creative work.
图1为本申请提供的基于深度学习的全髋关节置换术前规划方法的流程示意图;Fig. 1 is a schematic flow chart of the preoperative planning method for total hip arthroplasty based on deep learning provided by the present application;
图2为本申请提供的三维残差U-Net卷积神经网络的结构示意图;Fig. 2 is a schematic structural diagram of the three-dimensional residual U-Net convolutional neural network provided by the present application;
图3为本申请提供的3D卷积过程的示意图;Fig. 3 is a schematic diagram of the 3D convolution process provided by the present application;
图4为本申请提供的基于3D Res U-Net的股骨区域识别效果图;Fig. 4 is the femur region identification rendering based on 3D Res U-Net provided by the application;
图5为本申请提供的一种髋臼杯假体模型计划示意图;Fig. 5 is the schematic diagram of a kind of acetabular cup prosthesis model plan provided by the present application;
图6为本申请提供的一种股骨柄假体模型计划示意图;Fig. 6 is a schematic diagram of a femoral stem prosthesis model plan provided by the present application;
图7为本申请提供的基于3D Res U-Net的骨皮质区域识别效果图;Fig. 7 is the effect diagram of the bone cortex region recognition based on 3D Res U-Net provided by the present application;
图8为本申请提供的另一种基于深度学习的全髋关节置换术前规划方法的流程示意图;FIG. 8 is a schematic flowchart of another deep learning-based preoperative planning method for total hip arthroplasty provided by the present application;
图9为本申请提供的基于深度学习的全髋关节置换术前规划系统的结构示意图;FIG. 9 is a schematic structural diagram of a preoperative planning system for total hip arthroplasty based on deep learning provided by the present application;
图10为本申请提供的电子设备的结构示意图。FIG. 10 is a schematic structural diagram of an electronic device provided by the present application.
具体实施方式Detailed ways
为使本申请的目的、技术方案和优点更加清楚,下面将结合本申请中的附图,对本申请中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本申请一部分实施例,而不是全部的实施例。基于本申请中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本申请保护的范围。In order to make the purpose, technical solutions and advantages of this application clearer, the technical solutions in this application will be clearly and completely described below in conjunction with the accompanying drawings in this application. Obviously, the described embodiments are part of the embodiments of this application , but not all examples. Based on the embodiments in this application, all other embodiments obtained by persons of ordinary skill in the art without creative efforts fall within the protection scope of this application.
图1为本申请提供的基于深度学习的全髋关节置换术前规划方法的流程示意图,如图1所示,本申请提供了一种基于深度学习的全髋关节置换术前规划方法,包括:Figure 1 is a schematic flowchart of the preoperative planning method for total hip arthroplasty based on deep learning provided by this application. As shown in Figure 1, this application provides a preoperative planning method for total hip arthroplasty based on deep learning, including:
步骤101,获取待识别的全髋关节三维块状图,所述待识别的全髋关节三维块状图是由多张全髋关节二维横断面图像堆叠而成的。 Step 101 , acquiring a three-dimensional block diagram of the total hip joint to be identified, where the three-dimensional block diagram of the total hip joint to be identified is formed by stacking multiple two-dimensional cross-sectional images of the total hip joint.
在本申请中,首先,基于全髋关节术前规划的医学数字成像和通信(Digital Imaging and Communications in Medicine,简称DICOM)数据,构建待识别的全髋关节医学图像数据集,该图像数据集中的图像为全髋关节二维横断面图像,不同的全髋关节二维横断面图像对应不同的骨头区域,例如,髋关节二维横断面图像对应的骨头区域为股骨区域,膝关节二维横断面图像对应的骨头区域为胫骨区域(也包括一部分股骨区域),本申请以髋关节二维横断面图像进行说明,在各个髋关节二维横断面图像中包括有骨盆区域图像、左侧股骨图像和右侧股骨图像;然后,将髋关节二维横断面DICOM数据转换成JPG格式的图片,按照横断面的采集顺序(例如,从股骨头到股骨末端),将多张转换后的髋关节二维横断面图像堆叠成待识别的全髋关节三维块状图,即得到髋关节三维块状图。In this application, firstly, based on the digital imaging and communications in medicine (DICOM) data of the preoperative planning of the total hip joint, the medical image dataset of the total hip joint to be identified is constructed. The image is a two-dimensional cross-sectional image of the total hip joint. Different two-dimensional cross-sectional images of the total hip joint correspond to different bone regions. The bone area corresponding to the image is the tibia area (including a part of the femur area). Image of the right femur; then, convert the 2D cross-sectional DICOM data of the hip joint into a JPG image, and convert multiple converted hip joint The cross-sectional images are stacked to form a three-dimensional block map of the total hip joint to be identified, and the three-dimensional block map of the hip joint is obtained.
步骤102,将所述待识别的全髋关节三维块状图输入到训练好的三维分割神经网络中,获取所述训练好的三维分割神经网络输出的每张全髋关节二维横断面图像中的股骨区域,所述训练好的三维分割神经网络是由标注有股骨区域的标签的预设三维块状图,对卷积神经网络进行训练得到的。 Step 102, input the three-dimensional block diagram of the total hip joint to be identified into the trained three-dimensional segmentation neural network, and obtain each two-dimensional cross-sectional image of the total hip joint output by the trained three-dimensional segmentation neural network The femoral region, the trained three-dimensional segmentation neural network is obtained by training the convolutional neural network from a preset three-dimensional block diagram marked with a label of the femoral region.
在本申请中,通过训练好的三维分割神经网络,对堆叠而成的髋关节三维块状图进行识别。在本申请中,用于构建三维分割神经网络的卷积神经网络的卷积核为三维卷积核,使得基于该卷积神经网络训练得到的三维 分割神经网络,在进行图像分割时,可以对髋关节三维块状图的三个维度上进行特征提取,有效的提取了每个横截面之间的股骨区域特征信息,最终将股骨区域对应的三维块状图分割处理。由于三维分割神经网络在训练过程中,样本髋关节三维块状图的样本股骨区域已进行了标注,使得三维分割神经网络可直接对三维的髋关节图像进行股骨区域的识别。需要说明的是,本申请在对卷积神经网络进行训练时,将样本髋关节三维块状图中样本股骨区域的股骨头像素也进行了标注,使得模型在识别整个股骨区域的同时,还能进一步对股骨头区域进行识别,从而可快速识别髋关节中的股骨区域,节约时间成本,提高髋关节识别精度。需要说明的是,若需要对膝关节二维横断面图像进行图像识别,除了获取上述实施例提供的髋关节二维横断面图像以外,还包括膝关节二维横断面图像(可通过三维重建,构建胫骨区域和髌骨区域的三维图像),进而通过训练好的三维分割神经网络将股骨区域分割出来,再根据分割识别得到股骨区域末端,与胫骨及髌骨区域的三维图像,确定膝关节和胫骨的假体型号和安装位置。In this application, the stacked 3D block diagram of the hip joint is recognized through the trained 3D segmentation neural network. In this application, the convolution kernel of the convolutional neural network used to construct the three-dimensional segmentation neural network is a three-dimensional convolution kernel, so that the three-dimensional segmentation neural network obtained based on the convolutional neural network training can be used for image segmentation. Feature extraction is carried out in three dimensions of the three-dimensional block map of the hip joint, effectively extracting the feature information of the femoral region between each cross-section, and finally segmenting the three-dimensional block map corresponding to the femoral region. During the training process of the three-dimensional segmentation neural network, the sample femur region of the three-dimensional block image of the sample hip joint has been marked, so that the three-dimensional segmentation neural network can directly identify the femur region on the three-dimensional hip joint image. It should be noted that, when the application is training the convolutional neural network, the pixels of the femoral head in the sample femoral area in the three-dimensional block diagram of the sample hip joint are also marked, so that the model can identify the entire femoral area while also being able to The femoral head region is further identified, so that the femoral region in the hip joint can be quickly identified, saving time and cost, and improving the recognition accuracy of the hip joint. It should be noted that if it is necessary to perform image recognition on the two-dimensional cross-sectional image of the knee joint, in addition to obtaining the two-dimensional cross-sectional image of the hip joint provided in the above embodiment, it also includes a two-dimensional cross-sectional image of the knee joint (which can be reconstructed through three-dimensional, Construct the 3D images of the tibia region and the patella region), and then segment the femur region through the trained 3D segmentation neural network, and then obtain the end of the femur region, and the 3D images of the tibia and patella region according to the segmentation recognition, and determine the relationship between the knee joint and the tibia Prosthesis type and installation position.
步骤103,基于三维重建技术,根据每张全髋关节二维横断面图像中的股骨区域,得到所述股骨区域的三维图像。 Step 103, based on the three-dimensional reconstruction technology, according to the femoral region in each two-dimensional cross-sectional image of the total hip joint, obtain the three-dimensional image of the femoral region.
在本申请中,识别得到的股骨区域是由多张二维横断面图像堆叠而成,进而通过三维重建技术,利用多个二维图像中的关键点信息重建出三维信息,从而得到股骨三维图像。In this application, the identified femoral region is formed by stacking multiple two-dimensional cross-sectional images, and then three-dimensional information is reconstructed by using key point information in multiple two-dimensional images through three-dimensional reconstruction technology, thereby obtaining a three-dimensional image of the femur.
本申请提供的基于深度学习的全髋关节置换术前规划方法,通过将多张全髋关节二维横断面图像堆叠成全髋关节三维块状图,基于三维卷积核结构的卷积神经网络,对全髋关节三维块状图进行识别,有效的提取到每张全髋关节二维横断面图像的股骨区域,并根据提取得到的股骨区域进行三维建模,得到更为精准的股骨区域三维模型,以根据该股骨区域三维模型,在进行全髋关节置换术前规划时,提高全髋关节识别结果。The preoperative planning method for total hip arthroplasty based on deep learning provided by this application stacks multiple two-dimensional cross-sectional images of the total hip joint into a three-dimensional block diagram of the total hip joint, and based on the convolutional neural network of the three-dimensional convolution kernel structure, Identify the three-dimensional block diagram of the total hip joint, effectively extract the femoral region of each two-dimensional cross-sectional image of the total hip joint, and perform three-dimensional modeling based on the extracted femoral region to obtain a more accurate three-dimensional model of the femoral region , so as to improve the recognition result of the total hip joint during the preoperative planning of the total hip joint replacement according to the three-dimensional model of the femoral region.
在上述实施例的基础上,所述训练好的三维分割神经网络通过以下步骤得到:On the basis of the foregoing embodiments, the trained three-dimensional segmentation neural network is obtained through the following steps:
将预设三维块状图像集作为训练样本,输入到初始的三维分割神经网络进行训练,得到训练好的三维分割神经网络;The preset three-dimensional block image set is used as a training sample, input to the initial three-dimensional segmentation neural network for training, and a trained three-dimensional segmentation neural network is obtained;
其中,所述预设三维块状图像集中包含有若干张由预设二维横断面图像进行堆叠得到的预设三维块状图像。Wherein, the set of preset three-dimensional block images includes several preset three-dimensional block images obtained by stacking preset two-dimensional cross-sectional images.
在本申请中,预设三维块状图中不同的髓腔对应不同的骨头区域,根据实际骨头区域分割需求,构建相应的预设三维块状图。例如,股骨头及股骨柄假体的术前规划,预设三维块状图是由髋关节二维横断面图像堆叠而成;若需要进行膝关节及胫骨假体的术前规划,则预设三维块状图除了包括髋关节二维横断面图像,还包括相应的髌骨区域和胫骨区域的二维横断面图像,通过训练好的三维分割神经网络对股骨区域识别之后,再配合现有三维重建技术,根据髌骨区域和胫骨区域的二维横断面图像,构建髌骨和胫骨的三维图像,从而为膝关节及胫骨假体的术前规划提供更为精确的图像数据。In this application, different medullary cavities in the preset three-dimensional block diagram correspond to different bone regions, and the corresponding preset three-dimensional block diagram is constructed according to the actual bone region segmentation requirements. For example, for the preoperative planning of the femoral head and femoral stem prosthesis, the preset 3D block diagram is formed by stacking the 2D cross-sectional images of the hip joint; if the preoperative planning of the knee joint and tibial prosthesis is required, the preset In addition to the 2D cross-sectional image of the hip joint, the 3D block diagram also includes the corresponding 2D cross-sectional images of the patella region and the tibial region. After identifying the femur region through the trained 3D segmentation neural network, it is combined with the existing 3D reconstruction Technology, based on the two-dimensional cross-sectional images of the patella region and the tibia region, constructs a three-dimensional image of the patella and tibia, thereby providing more accurate image data for the preoperative planning of the knee joint and tibial prosthesis.
在上述实施例的基础上,训练过程具体包括:On the basis of the foregoing embodiments, the training process specifically includes:
步骤201,获取若干张样本髋关节二维横断面图像;Step 201, acquiring several two-dimensional cross-sectional images of the sample hip joint;
步骤202,对每一张样本髋关节二维横断面图像中样本股骨区域进行标注,并对所述样本股骨区域的股骨头像素标注股骨头区域标签,得到若干张第一预设髋关节二维横断面图像;Step 202, labeling the sample femoral region in each sample hip joint two-dimensional cross-sectional image, and marking the femoral head region label on the femoral head pixels of the sample femoral region, and obtaining several first preset two-dimensional hip joint images. cross-sectional images;
步骤203,根据样本髋关节二维横断面的采集顺序,将若干张所述第一预设髋关节二维横断面图像进行堆叠,得到对应的第一预设三维块状图像;Step 203: According to the acquisition sequence of the sample hip joint two-dimensional cross-section, stack several first preset two-dimensional cross-sectional images of the hip joint to obtain a corresponding first preset three-dimensional block image;
步骤204,将多个所述第一预设三维块状图像输入到所述初始的三维分割神经网络进行训练,得到训练好的三维分割神经网络;Step 204, inputting a plurality of the first preset three-dimensional block images to the initial three-dimensional segmentation neural network for training to obtain a trained three-dimensional segmentation neural network;
其中,所述初始的三维分割神经网络是由U-Net卷积网络构建得到的,且所述U-Net卷积网络的卷积核为三维卷积核。Wherein, the initial three-dimensional segmentation neural network is constructed by the U-Net convolution network, and the convolution kernel of the U-Net convolution network is a three-dimensional convolution kernel.
在本申请中,获取髋关节医学图像样本数据集,对样本集中的样本髋关节二维横断面图像进行手动标注样本股骨区域,只提取含有股骨部分的标签作为掩膜mask。由于本申请在提取股骨区域的同时,还可以对股骨区域的股骨头部分进行提取,因此,在标注时也对样本股骨区域中的股骨头像素进行标注,最终通过将标注后的样本髋关节二维横断面图像(即第 一预设髋关节二维横断面图像)进行堆叠,从而构建得到第一预设三维块状图像。具体地,本申请将第一预设髋关节二维横断面DICOM数据转换成JPG格式的图片,同时,标注文件转换成PNG格式的图片,将其进行打乱顺序后,按照6:2:2的比例划分为训练集、验证集和测试集。由于本申请的卷积神经网络的输入图像比现有2D网络的输入图像多一个维度,即样本髋关节三维块状图(第一预设三维块状图像)是多张二维横截面图像堆叠而成的,因此,标注文件与该样本髋关节三维块状图对应,也是一个块状图。需要说明的是,在本申请中,第一预设三维块状图像是通过预先提取得到股骨部分的二维横断面堆叠而成,并在股骨部分堆叠的三维块状图上对股骨头区域进行标注,使得分割网络在训练完成之后,更加快速的将图像中的股骨头区域与背景分割。In this application, the hip joint medical image sample data set is obtained, the sample hip joint two-dimensional cross-sectional image in the sample set is manually labeled with the sample femur region, and only the label containing the femur part is extracted as a mask mask. Since this application can also extract the femoral head part of the femoral region while extracting the femoral region, the pixels of the femoral head in the sample femoral region are also marked during labeling, and finally the marked sample hip joint two The three-dimensional cross-sectional images (that is, the first preset two-dimensional cross-sectional image of the hip joint) are stacked to obtain the first preset three-dimensional block image. Specifically, this application converts the DICOM data of the first preset two-dimensional cross-section of the hip joint into a picture in JPG format, and at the same time, converts the annotation file into a picture in PNG format. The ratio is divided into training set, validation set and test set. Since the input image of the convolutional neural network of the present application has one more dimension than the input image of the existing 2D network, that is, the three-dimensional block diagram of the sample hip joint (the first preset three-dimensional block image) is formed by stacking multiple two-dimensional cross-sectional images Therefore, the annotation file corresponds to the 3D block diagram of the sample hip joint, which is also a block diagram. It should be noted that, in this application, the first preset three-dimensional block image is obtained by pre-extracting and stacking two-dimensional cross-sections of the femoral part, and the femoral head region is calculated on the three-dimensional block image of the femoral part stack. Labeling enables the segmentation network to more quickly segment the femoral head region in the image from the background after the training is completed.
图2为本申请提供的三维残差U-Net卷积神经网络的结构示意图,可参考图2所示,在本申请中,采用的三维残差U-Net卷积神经网络(简称3D Res UNet)是基于U-Net构建得到的,同样包含了一个Encoder(编码器)部分和一个Decoder(解码器)部分,其中,Encoder部分是用来分析整张图片并且进行特征提取与分析,而与之相对应的Decoder部分是还原特征的过程,通过对该卷积神经网络进行训练,得到的三维分割神经网络,可对待识别的三维块状图中的股骨区域进行分割,得到分割好的股骨区域块状图。Fig. 2 is a schematic structural diagram of the three-dimensional residual U-Net convolutional neural network provided by the application, as shown in Fig. 2, in this application, the three-dimensional residual U-Net convolutional neural network (referred to as 3D Res UNet ) is built based on U-Net, and also includes an Encoder (encoder) part and a Decoder (decoder) part, where the Encoder part is used to analyze the entire picture and perform feature extraction and analysis, and with it The corresponding Decoder part is the process of restoring features. By training the convolutional neural network, the obtained 3D segmentation neural network can segment the femur region in the 3D block image to be recognized, and obtain the segmented femoral region block state diagram.
具体地,Encoder部分是由基本的残差模块ResBlock和最大池化层MaxPooling构成,其输入尺寸(Input shape)为1*8*256*256*1,ResBlock由基本的卷积块组成,包括两个连续的3D卷积核,每个3D卷积核包括两组操作:2*(conv+relu+bn),即将3D卷积(conv)、激活函数(relu)和批量归一化(Batch Normalization,简称BN)构成一组操作,其中,3D卷积可以有效的提取到横截面之间的信息,减少误检率;激活函数可以增加模型的非线性能力,提升模型特征提取能力;批量归一化操作可以改变数据的分布,有利于网络训练时快速收敛。图3为本申请提供的3D卷积过程的示意图,可参考图3所示,由于卷积核是三维的,可以在三个维度上以一定的步长滑动,其将周围3*3*3区域的信息结合在一起形成一个点,能够提取到更丰富的特征信息。Specifically, the Encoder part is composed of the basic residual module ResBlock and the maximum pooling layer MaxPooling. Its input size (Input shape) is 1*8*256*256*1. ResBlock is composed of basic convolutional blocks, including two A continuous 3D convolution kernel, each 3D convolution kernel includes two sets of operations: 2*(conv+relu+bn), that is, 3D convolution (conv), activation function (relu) and batch normalization (Batch Normalization , referred to as BN) constitute a set of operations, among which, 3D convolution can effectively extract the information between the cross-sections and reduce the false detection rate; the activation function can increase the nonlinear ability of the model and improve the feature extraction ability of the model; batch normalization The optimization operation can change the distribution of data, which is conducive to the rapid convergence of network training. Figure 3 is a schematic diagram of the 3D convolution process provided by this application, as shown in Figure 3, since the convolution kernel is three-dimensional, it can slide in a certain step in three dimensions, and it will surround 3*3*3 The information of the region is combined to form a point, which can extract richer feature information.
优选地,所述卷积神经网络中两个连续的三维卷积核是通过残差结构连接的。可参考图2所示,在本申请中,两个连续的3D卷积核通过残差连接(Skip Connection),可以防止网络退化。通过残差结构连接两个3D卷积核,将网络中的特征图经过连续卷积块作用前和作用后相加(Add)在一起,这样可以使网络在训练时自行选择合适的反向传播路径。另外,Maxpooling操作将相邻的固定大小区域的最大值作为该区域特征表示,这样做可以有效减少网络参数运算。Preferably, two consecutive three-dimensional convolution kernels in the convolutional neural network are connected through a residual structure. As shown in Figure 2, in this application, two consecutive 3D convolution kernels can prevent network degradation through a residual connection (Skip Connection). Connect the two 3D convolution kernels through the residual structure, and add the feature maps in the network before and after the continuous convolution blocks, so that the network can choose the appropriate backpropagation during training. path. In addition, the Maxpooling operation uses the maximum value of the adjacent fixed-size region as the characteristic representation of the region, which can effectively reduce the calculation of network parameters.
进一步地,Decoder部分是由上采样和卷积块(反卷积Deconv)构成。在Decoder层,会先将Encoder的特征图与上采样过后的特征图进行通道堆叠(Concat),然后经过Decoder部分的卷积块进行反卷积,其输出尺寸(output shape)为1*8*256*256*1,其中,卷积块的构成与Encoder过程相同,此处不再赘述。Further, the Decoder part is composed of upsampling and convolution blocks (deconvolution Deconv). In the Decoder layer, the Encoder feature map and the upsampled feature map will be channel-stacked (Concat), and then deconvolved through the convolution block of the Decoder part, and its output shape (output shape) is 1*8* 256*256*1, where the composition of the convolution block is the same as the Encoder process, and will not be repeated here.
在本申请中,3D Res UNet在训练时,每一次喂入网络的样本数batch size为8,初始化学习率设置为1e-4,每隔5000个iteration(迭代),学习率衰减为原来的0.95。优选地,在本申请中,优化器使用Adam优化器,所述卷积神经网络的损失函数是由DICE loss和BCE loss组成的(DICE loss和BCE loss分别为一种损失函数的类型)。由于使用的损失函数为DICE loss和BCE loss的融合体,这样可以避免只使用DICE loss时造成网络训练过程中的振荡现象。In this application, when 3D Res UNet is training, the batch size of samples fed to the network each time is 8, the initial learning rate is set to 1e-4, and every 5000 iterations (iterations), the learning rate decays to the original 0.95 . Preferably, in the present application, the optimizer uses an Adam optimizer, and the loss function of the convolutional neural network is composed of DICE loss and BCE loss (DICE loss and BCE loss are respectively a type of loss function). Since the loss function used is a fusion of DICE loss and BCE loss, this can avoid oscillations during network training when only DICE loss is used.
进一步地,基于前期划分好的训练集、验证集和测试集,设置每迭代1000次,对训练集和验证集做一次验证,并测算模型的train loss(训练误差)和val loss(验证误差),以及train DICE和val DICE,通过早停法来判断网络训练停止时机,得到训练好的三维分割神经网络。另外,在测试阶段,将整例样本病例的髋关节二维横断面DICOM数据和标注文件,按照顺序分别转换为JPG、PNG格式的图片,并包装成图像块(即块状图),通过测试得到test DICE。Further, based on the training set, verification set and test set divided in the previous period, set each iteration 1000 times, do a verification on the training set and verification set, and measure the model's train loss (training error) and val loss (verification error) , as well as train DICE and val DICE, use the early stopping method to judge the timing of network training stop, and get the trained three-dimensional segmentation neural network. In addition, in the testing phase, the DICOM data and annotation files of the two-dimensional cross-section of the hip joint of the entire sample case were converted into images in JPG and PNG formats in sequence, and packaged into image blocks (ie, block diagrams), which passed the test. Get test DICE.
在上述实施例的基础上,在所述基于三维重建技术,根据每张全髋关节二维横断面图像中的股骨区域,得到所述股骨区域的三维图像之后,所述方法还包括:On the basis of the above-mentioned embodiments, after obtaining the three-dimensional image of the femoral region according to the femoral region in each two-dimensional cross-sectional image of the total hip joint based on the three-dimensional reconstruction technology, the method further includes:
基于平面图像的质心公式,获取所述三维图像中股骨区域的股骨头中心点的像素坐标;Based on the centroid formula of the planar image, obtain the pixel coordinates of the center point of the femoral head in the femoral region in the three-dimensional image;
将所述像素坐标转换为图像坐标;converting said pixel coordinates to image coordinates;
确定股骨头旋转中心位置;Determine the position of the center of rotation of the femoral head;
根据股骨头旋转中心位置,获取第一尺寸信息,以根据所述第一尺寸确定第二尺寸信息,其中,所述第一尺寸信息为股骨头对应的尺寸信息,所述第二尺寸信息为髋臼杯假体模型对应的尺寸信息。可选地,在本申请中,尺寸信息至少包括直径和半径,例如,通过获取股骨头的直径,确定髋臼杯假体模型的直径。According to the rotation center position of the femoral head, the first size information is obtained to determine the second size information according to the first size, wherein the first size information is the size information corresponding to the femoral head, and the second size information is the hip The size information corresponding to the acetabular cup prosthesis model. Optionally, in the present application, the size information includes at least a diameter and a radius, for example, by obtaining the diameter of the femoral head, the diameter of the acetabular cup prosthesis model is determined.
在本申请中,通过3D Res U-Net识别髋关节三维块状图的每个像素区域,在训练过程中,可将像素标注划分为两种属性值,分别命名0和1,其中,数值0代表背景像素,1代表股骨头像素;在完成标注之后,将标注好的图像数据传入到卷积神经网络中(即3D Res UNet),使用卷积池化采样一直迭代学习训练。图4为本申请提供的基于3D Res U-Net的股骨区域识别效果图,如图4所示,在训练完成之后,训练好的三维分割神经网络可自动识别股骨头位置,完成股骨头区域的识别(图4中用白色线框示意标出)。最后通过三维重建,得到股骨头区域的三维图像。In this application, 3D Res U-Net is used to identify each pixel area of the three-dimensional block map of the hip joint. During the training process, the pixel label can be divided into two attribute values, named 0 and 1 respectively, where the value 0 Represents the background pixel, 1 represents the femoral head pixel; after the labeling is completed, the marked image data is passed into the convolutional neural network (ie 3D Res UNet), and the convolution pooling sampling is used to iteratively learn and train. Fig. 4 is the effect diagram of the femoral region recognition based on 3D Res U-Net provided by the present application. As shown in Fig. 4, after the training is completed, the trained three-dimensional segmentation neural network can automatically recognize the position of the femoral head and complete the femoral head region identification Identification (marked schematically with a white line in Figure 4). Finally, through 3D reconstruction, a 3D image of the femoral head region is obtained.
在获取到股骨头区域的三维图像之后,由于三维分割神经网络输出的图像为二值图像,在二值图像中,仅有“0”和“1”两种像素值,其质量分布是均匀的,所以在识别得到股骨区域的三维图像中,质心和形心重合,根据平面图像的质心公式,可求得股骨区域的三维图像中股骨头的中心点坐标,即该三维图像中的股骨头旋转中心。具体地,设二值图像为B,则B[i,j]表示二值图像B中第i行,第j列像素的像素值,因此可使用下列公式求得三维图像中股骨区域的股骨头中心点的位置:After obtaining the 3D image of the femoral head region, since the image output by the 3D segmentation neural network is a binary image, in the binary image, there are only two pixel values of "0" and "1", and its mass distribution is uniform , so in the recognized 3D image of the femoral region, the centroid and centroid coincide. According to the centroid formula of the planar image, the coordinates of the center point of the femoral head in the 3D image of the femoral region can be obtained, that is, the femoral head rotation in the 3D image center. Specifically, if the binary image is B, then B[i, j] represents the pixel value of the i-th row and j-th column pixel in the binary image B, so the following formula can be used to obtain the femoral head of the femoral region in the three-dimensional image The location of the center point:
Figure PCTCN2023070788-appb-000001
Figure PCTCN2023070788-appb-000001
Figure PCTCN2023070788-appb-000002
Figure PCTCN2023070788-appb-000002
其中,A表示二值图像中所有像素的像素值之和,Among them, A represents the sum of the pixel values of all pixels in the binary image,
Figure PCTCN2023070788-appb-000003
n表示二值图像中像素的最大行数,m表示二值图像中像素的最大列数,从而可得到的股骨头中心点的像素坐标,进而将像素坐标转换为图像坐标。
Figure PCTCN2023070788-appb-000003
n represents the maximum number of rows of pixels in the binary image, and m represents the maximum number of columns of pixels in the binary image, so that the pixel coordinates of the center point of the femoral head can be obtained, and then the pixel coordinates are converted into image coordinates.
具体地,图像平面坐标中心坐标为:Specifically, the center coordinates of the image plane coordinates are:
Figure PCTCN2023070788-appb-000004
Figure PCTCN2023070788-appb-000004
则像素坐标
Figure PCTCN2023070788-appb-000005
到图像坐标(x′,y′)的变换公式为:
Then the pixel coordinates
Figure PCTCN2023070788-appb-000005
The transformation formula to image coordinates (x', y') is:
Figure PCTCN2023070788-appb-000006
Figure PCTCN2023070788-appb-000006
Figure PCTCN2023070788-appb-000007
Figure PCTCN2023070788-appb-000007
其中,S x,S y分别为图像阵列的行列间距。最后根据图像坐标,得到股骨头旋转中心在三维图像的位置,确定股骨头直径,由于股骨头直径等同于髋臼杯内圆直径,进而推算出髋臼杯假体模型的直径。示例地,股骨头直径是以股骨头中心为出发点,通过一条射线与基于本申请提供的全髋关节置换术前规划方法得到的股骨头区域边缘相交,计算交点到股骨头中心的长度,然后以一度为步长,每旋转一度计算一次,最后通过统计均值得到半径,从而得到髋臼杯假体模型的直径,保证了根据模型的识别结果确定髋臼杯假体模型的规格型号(示例的,图5提供了一种髋臼杯假体模型计划示意图)。 Wherein, S x , S y are the row and column spacing of the image array respectively. Finally, according to the image coordinates, the position of the center of rotation of the femoral head in the three-dimensional image is obtained, and the diameter of the femoral head is determined. Since the diameter of the femoral head is equal to the diameter of the inner circle of the acetabular cup, the diameter of the acetabular cup prosthesis model is calculated. Exemplarily, the diameter of the femoral head is based on the center of the femoral head, intersects a ray with the edge of the femoral head region obtained based on the preoperative planning method for total hip arthroplasty provided by this application, calculates the length from the intersection point to the center of the femoral head, and then calculates the length from the intersection point to the center of the femoral head One degree is the step length, and it is calculated once per one degree of rotation. Finally, the radius is obtained through the statistical mean value, so as to obtain the diameter of the acetabular cup prosthesis model, which ensures that the specification of the acetabular cup prosthesis model is determined according to the recognition results of the model (for example, Figure 5 provides a schematic diagram of the model plan of the acetabular cup prosthesis).
在上述实施例的基础上,在所述将多个所述第一预设三维块状图像输入到所述初始的三维分割神经网络进行训练,得到训练好的三维分割神经网络之后,所述方法还包括:On the basis of the above embodiments, after inputting a plurality of the first preset three-dimensional block images into the initial three-dimensional segmentation neural network for training and obtaining a trained three-dimensional segmentation neural network, the method Also includes:
对若干张所述第一预设髋关节二维横断面图像中样本股骨区域的骨皮质像素标注骨皮质区域标签,得到若干张第二预设髋关节二维横断面图像;Marking the cortical bone region labels on the cortical bone pixels of the sample femur region in the first preset two-dimensional cross-sectional images of the hip joint to obtain several second preset two-dimensional cross-sectional images of the hip joint;
将若干张所述第二预设髋关节二维横断面图像按顺序进行堆叠,构建对应的第二预设三维块状图;Stacking several second preset hip joint two-dimensional cross-sectional images in order to construct a corresponding second preset three-dimensional block diagram;
通过多个所述第二预设三维块状图,对所述训练好的三维分割神经网 络进行微调,得到髋关节识别模型。Through a plurality of the second preset three-dimensional block diagrams, the trained three-dimensional segmentation neural network is fine-tuned to obtain a hip joint recognition model.
在本申请中,通过在第一预设髋关节二维横断面图像中标记新的标签,从而根据得到的新训练集对三维分割神经网络的参数进行优化,使得模型在识别股骨头区域的同时,还能识别骨皮质区域。In this application, by marking a new label in the first preset two-dimensional cross-sectional image of the hip joint, the parameters of the three-dimensional segmentation neural network are optimized according to the obtained new training set, so that the model recognizes the femoral head region while , but also to identify cortical areas.
在上述实施例的基础上,在所述获取待识别的全髋关节三维块状图之后,所述方法还包括:On the basis of the above embodiments, after the acquisition of the three-dimensional block diagram of the total hip joint to be identified, the method further includes:
将所述髋关节三维块状图输入到所述髋关节识别模型,得到每张髋关节二维横断面图像中的股骨区域和骨皮质区域;Inputting the three-dimensional block diagram of the hip joint into the hip joint recognition model to obtain the femoral region and the cortical bone region in each two-dimensional cross-sectional image of the hip joint;
根据所述股骨区域和所述骨皮质区域,确定髓腔区域;determining the medullary cavity area according to the femoral area and the cortical bone area;
计算所述髓腔区域中每个髓腔层面的中点坐标,并根据所述中点坐标,将所有中心点进行直线拟合,确定髓腔解剖轴线;calculating the midpoint coordinates of each medullary canal level in the medullary cavity region, and performing straight line fitting on all the midpoints according to the midpoint coordinates to determine the anatomical axis of the medullary cavity;
根据所述髓腔解剖轴线和股骨颈轴线,计算得到股骨颈干角的角度值;Calculate the angle value of the femoral neck-shaft angle according to the anatomical axis of the medullary cavity and the femoral neck axis;
根据所述角度值、所述髓腔区域和股骨头旋转中心位置,确定股骨柄假体模型的类型和放置位置(示例的,图6提供了一种股骨柄假体模型计划示意图)。According to the angle value, the medullary cavity area and the position of the center of rotation of the femoral head, the type and placement position of the femoral stem prosthesis model are determined (for example, FIG. 6 provides a schematic diagram of a femoral stem prosthesis model plan).
在本申请中,参数优化后的三维分割神经网络,即髋关节识别模型,在识别图像的每个像素区域时,像素标注划分为三种属性值,分别命名0、1和2,其中,数值0代表背景像素,1代表股骨头像素,2代表骨皮质。图7为本申请提供的基于3D Res U-Net的骨皮质区域识别效果图,如图7所示,在对三维分割神经网络的参数进行优化之后,可实现同时对块状图中的股骨头和骨皮质的识别(图7中用黑色线条示意标出)。In this application, the three-dimensional segmentation neural network after parameter optimization, that is, the hip joint recognition model, when recognizing each pixel region of the image, the pixel label is divided into three attribute values, named 0, 1 and 2 respectively, where the value 0 represents background pixels, 1 represents femoral head pixels, and 2 represents cortical bone. Figure 7 is the effect diagram of the recognition of the cortical bone region based on 3D Res U-Net provided by the present application. As shown in Figure 7, after optimizing the parameters of the three-dimensional segmentation neural network, the femoral head in the block diagram can be simultaneously identified. and identification of cortical bone (schematically marked with a black line in Figure 7).
进一步地,本申请从髋关节识别模型输出的识别结果中,截取小转子结束处直到股骨末端部位的图像,将图像中股骨区域减去骨皮质区域,得到髓腔区域。然后,从小转子结束位置以下,每横行(即从每个髋关节二维横断面中识别得到髓腔层面,每一横行指的是通过本申请分割后的髋关节二维横截面图像,经过三维重建模拟X射线投影效果图,然后在小转子结束位置的下方开始,每隔一定预设位置在该图像上画一条横线,横线与两个股骨髓腔边缘相交则会得到四个点)与髓腔交点为四个坐标,从左至 右分别命名为A1,A2,B1,B2;依据两点可以求出中点,即A1(X 1,Y 1),A2(X 2,Y 2)的中点坐标,可通过以下公式得到: Further, the present application intercepts images from the end of the lesser trochanter to the end of the femur from the recognition results output by the hip joint recognition model, and subtracts the cortical bone region from the femoral region in the image to obtain the medullary cavity region. Then, below the end position of the lesser trochanter, each row (i.e., the level of the medullary canal is identified from the two-dimensional cross-section of each hip joint, and each row refers to the two-dimensional cross-sectional image of the hip joint segmented by this application, after three-dimensional Reconstruct the simulated X-ray projection effect map, and then start below the end position of the lesser trochanter, draw a horizontal line on the image at certain preset positions, and when the horizontal line intersects with the edges of the two femoral medullary canals, four points will be obtained) and The intersection of the medullary cavity is four coordinates, which are named A1, A2, B1, and B2 from left to right; the midpoint can be obtained based on two points, namely A1 (X 1 , Y 1 ), A2 (X 2 , Y 2 ) The midpoint coordinates of can be obtained by the following formula:
Figure PCTCN2023070788-appb-000008
Figure PCTCN2023070788-appb-000008
B1,B2同理可算得。每行依次算得髓腔的中点坐标,将这些点拟合成一条直线,从而确定为髓腔解剖轴线。最后,根据髓腔解剖轴线和股骨颈轴线,计算得到股骨颈干角的角度值,再结合髓腔形态和股骨头旋转中心位置,可共同确定股骨柄假体模型的型号与放置位置。具体地,股骨颈干角为髓腔解剖轴线和股骨颈轴线之间的夹角,以这些获得的参数为筛选条件,对已有的假体库进行过滤,通过模板匹配获取假体库中对应的最佳型号;进一步地,通过移动股骨柄假体,将股骨柄假体的旋转中心位置,与之前计算得到的髋臼杯旋转中心位置(即股骨头旋转中心位置)重合,得到股骨柄假体实际放置位置。本申请基于卷积神经网络识别得到股骨区域,可快速且准确的确定假体模型的规格型号和安放位置。B1 and B2 can be calculated in the same way. The coordinates of the midpoint of the medullary cavity are calculated in turn for each row, and these points are fitted into a straight line to determine the anatomical axis of the medullary cavity. Finally, according to the anatomical axis of the medullary cavity and the axis of the femoral neck, the angle value of the femoral neck-shaft angle was calculated, combined with the shape of the medullary cavity and the position of the center of rotation of the femoral head, the model and placement position of the femoral stem prosthesis model could be determined together. Specifically, the femoral neck-shaft angle is the angle between the anatomical axis of the medullary canal and the femoral neck axis. With these obtained parameters as the screening conditions, the existing prosthesis library is filtered, and the corresponding The optimal model; further, by moving the femoral stem prosthesis, the rotation center position of the femoral stem prosthesis coincides with the previously calculated acetabular cup rotation center position (that is, the femoral head rotation center position), and the femoral stem prosthesis is obtained. The actual location of the body. The application obtains the femoral region based on convolutional neural network identification, which can quickly and accurately determine the specification, model and placement position of the prosthesis model.
在确定出髋臼杯假体第二尺寸信息和股骨柄假体的型号和位置后,分别模拟将髋臼杯假体和股骨柄假体置入目标位置(如图5和图6所示)。若髋臼杯假体的的放置位置、股骨柄假体的放置位置均满足预设位置要求时,则可以输出全髋关节术前规划方案。其中,髋臼杯假体的预设位置要求例如可以是髋臼杯假体放入髋臼窝后,覆盖髋臼窝的覆盖率大于75%,股骨柄假体的预设位置要求例如可以是股骨柄假体放入髓腔后,股骨柄假体长轴与股骨长轴之间的角度小于或等于3°。After determining the second size information of the acetabular cup prosthesis and the type and position of the femoral stem prosthesis, simulate the placement of the acetabular cup prosthesis and the femoral stem prosthesis into the target position (as shown in Figure 5 and Figure 6) . If the placement position of the acetabular cup prosthesis and the placement position of the femoral stem prosthesis meet the preset position requirements, the preoperative planning plan for the total hip joint can be output. Among them, the preset position requirements of the acetabular cup prosthesis can be, for example, that after the acetabular cup prosthesis is placed in the acetabular socket, the coverage rate covering the acetabular socket is greater than 75%, and the preset position requirements of the femoral stem prosthesis can be, for example, After the femoral stem prosthesis is placed in the medullary cavity, the angle between the long axis of the femoral stem prosthesis and the long axis of the femur is less than or equal to 3°.
接下来,参加图8,对本申请的上述技术方案进行整体说明。本申请提供的基于深度学习的全髋关节置换术前规划方法,包括:Next, referring to FIG. 8 , the above technical solution of the present application will be described as a whole. The preoperative planning method for total hip replacement based on deep learning provided by this application includes:
步骤801:获取待识别的全髋关节三维块状图。Step 801: Obtain the three-dimensional block diagram of the total hip joint to be identified.
步骤802:将待识别的全髋关节三维块状图输入到训练好的三维分割神经网络中,得到每张全髋关节二维横断面图像中的股骨区域。Step 802: Input the 3D block image of the total hip joint to be identified into the trained 3D segmentation neural network to obtain the femur region in each 2D cross-sectional image of the total hip joint.
步骤803:基于三维重建技术,根据每张全髋关节二维横断面图像中的股骨区域,得到股骨区域的三维图像。Step 803: Based on the three-dimensional reconstruction technology, according to the femoral region in each two-dimensional cross-sectional image of the total hip joint, a three-dimensional image of the femoral region is obtained.
步骤804:确定髋臼杯假体的尺寸信息,以及股骨柄假体的类型和放 置位置。Step 804: Determine the size information of the acetabular cup prosthesis, as well as the type and placement position of the femoral stem prosthesis.
步骤805:分别将确定后髋臼杯假体和股骨柄假体放入髋臼窝和髓腔内。Step 805: Put the determined posterior acetabular cup prosthesis and femoral stem prosthesis into the acetabular fossa and medullary cavity respectively.
步骤806:对放入的髋臼杯假体的位置或者型号,以及对股骨柄假体的位置或者型号进行矫正。Step 806: Correcting the position or type of the inserted acetabular cup prosthesis and the position or type of the femoral stem prosthesis.
步骤807:输出全髋关节置换术前规划方案。Step 807: Output the preoperative planning scheme for total hip replacement.
下面对本申请提供的基于深度学习的全髋关节置换术前规划系统进行描述,下文描述的基于深度学习的全髋关节置换术前规划系统与上文描述的基于深度学习的全髋关节置换术前规划方法可相互对应参照。The following describes the preoperative planning system for total hip replacement based on deep learning provided by this application. The preoperative planning system for total hip replacement based on deep learning described below is the same as the preoperative planning system for total hip replacement based on deep learning described above Planning methods can be cross-referenced.
图9为本申请提供的基于深度学习的全髋关节置换术前规划系统的结构示意图,如图9所示,本申请提供了一种基于深度学习的全髋关节置换术前规划系统,包括全髋关节图像采集模块901、全髋关节识别模块902和全髋关节三维图像构建模块903,其中,全髋关节图像采集模块901用于获取待识别的全髋关节三维块状图,所述全髋关节三维块状图是由多张全髋关节二维横断面图像堆叠而成的;全髋关节识别模块902用于将所述待识别的全髋关节三维块状图输入到训练好的三维分割神经网络中,获取所述训练好的三维分割神经网络输出的每张全髋关节二维横断面图像中的股骨区域,所述训练好的三维分割神经网络是由标注有股骨区域的标签的预设三维块状图,对卷积神经网络进行训练得到的;全髋关节三维图像构建模块903用于基于三维重建技术,根据每张全髋关节二维横断面图像中的股骨区域,得到所述股骨区域的三维图像。Fig. 9 is a schematic structural diagram of the preoperative planning system for total hip replacement based on deep learning provided by the present application. As shown in Fig. 9, the present application provides a preoperative planning system for total hip replacement based on deep learning, including A hip joint image collection module 901, a total hip joint recognition module 902 and a total hip joint three-dimensional image construction module 903, wherein the total hip joint image collection module 901 is used to obtain a three-dimensional block diagram of the total hip joint to be identified, and the total hip joint The three-dimensional block diagram of the joint is formed by stacking multiple two-dimensional cross-sectional images of the total hip joint; the total hip joint identification module 902 is used to input the three-dimensional block diagram of the total hip joint to be identified into the trained three-dimensional segmentation In the neural network, the femur region in each two-dimensional cross-sectional image of the total hip joint output by the trained three-dimensional segmentation neural network is obtained, and the trained three-dimensional segmentation neural network is a pre-prepared image of a label marked with the femoral region. Assuming a three-dimensional block diagram, it is obtained by training a convolutional neural network; the three-dimensional image construction module 903 of the total hip joint is used to obtain the described 3D image of the femoral region.
在本申请中,以髋关节二维横断面图像进行说明,该二维横断面图像为DICOM数据,在各个髋关节二维横断面图像中包括有骨盆区域图像、左侧股骨图像和右侧股骨图像。进一步地,全髋关节图像采集模块901将髋关节二维横断面DICOM数据转换成JPG格式的图片,按照横断面的采集顺序(例如,从股骨头到股骨末端),将多张转换后的髋关节二维横断面图像进行堆叠,从而生成待识别的全髋关节三维块状图,即得到髋关节三维块状图。In this application, a two-dimensional cross-sectional image of the hip joint is used for illustration. The two-dimensional cross-sectional image is DICOM data, and each two-dimensional cross-sectional image of the hip joint includes an image of the pelvic region, an image of the left femur, and an image of the right femur. image. Further, the total hip joint image acquisition module 901 converts the DICOM data of the two-dimensional cross-section of the hip joint into a picture in JPG format, and according to the acquisition order of the cross-section (for example, from the femoral head to the end of the femur), multiple converted hip joints The two-dimensional cross-sectional images of the joints are stacked to generate the three-dimensional block diagram of the total hip joint to be identified, that is, the three-dimensional block diagram of the hip joint is obtained.
在本申请中,全髋关节识别模块902内配置有训练好的三维分割神经 网络,可对堆叠而成的髋关节三维块状图进行识别。在本申请中,用于构建三维分割神经网络的卷积神经网络的卷积核为三维卷积核,使得基于该卷积神经网络训练得到的三维分割神经网络,在进行图像分割时,可以对髋关节三维块状图的三个维度上进行特征提取,有效的提取了每个横截面之间的股骨区域特征信息,最终将股骨区域对应的三维块状图分割处理。优选地,所述卷积神经网络中两个连续的三维卷积核是通过残差结构连接的;且所述卷积神经网络的损失函数是由DICE loss和BCE loss组成的。In this application, the trained three-dimensional segmentation neural network is configured in the total hip joint identification module 902, which can identify the stacked three-dimensional block diagrams of the hip joint. In this application, the convolution kernel of the convolutional neural network used to construct the three-dimensional segmentation neural network is a three-dimensional convolution kernel, so that the three-dimensional segmentation neural network obtained based on the convolutional neural network training can be used for image segmentation. Feature extraction is carried out in three dimensions of the three-dimensional block map of the hip joint, effectively extracting the feature information of the femoral region between each cross-section, and finally segmenting the three-dimensional block map corresponding to the femoral region. Preferably, two continuous three-dimensional convolution kernels in the convolutional neural network are connected through a residual structure; and the loss function of the convolutional neural network is composed of DICE loss and BCE loss.
由于三维分割神经网络在训练过程中,样本髋关节三维块状图的股骨区域已进行了标注,使得全髋关节识别模块902中的三维分割神经网络,可直接对三维的髋关节图像进行股骨区域的识别。需要说明的是,本申请在对卷积神经网络进行训练时,将样本髋关节三维块状图中样本股骨区域的股骨头像素也进行了标注,使得全髋关节识别模块802在识别整个股骨区域的同时,还能进一步对股骨头区域进行识别,从而可快速识别髋关节中的股骨区域,节约时间成本,提高髋关节识别精度。During the training process of the three-dimensional segmentation neural network, the femur region of the sample hip joint three-dimensional block diagram has been marked, so that the three-dimensional segmentation neural network in the total hip joint recognition module 902 can directly perform femur region analysis on the three-dimensional hip joint image. recognition. It should be noted that, when the present application is training the convolutional neural network, the pixels of the femoral head in the sample femoral area in the three-dimensional block diagram of the sample hip joint are also marked, so that the total hip joint identification module 802 can identify the entire femoral area At the same time, it can further identify the femoral head region, so that the femoral region in the hip joint can be quickly identified, saving time and cost, and improving the recognition accuracy of the hip joint.
最后,通过全髋关节三维图像构建模块903,利用识别得到的股骨区域的关键点信息进行三维重建,得到股骨三维图像,从而在后续的全髋关节置换术前规划中,可基于该股骨三维图像确定全髋关节相关的假体尺寸类型和放置位置。Finally, through the total hip joint three-dimensional image construction module 903, the key point information of the identified femoral region is used to perform three-dimensional reconstruction to obtain a three-dimensional image of the femur, so that in the subsequent preoperative planning of total hip joint replacement, the three-dimensional image of the femur can be Determine the size, type and placement of the prosthesis relative to the total hip.
本申请提供的基于深度学习的全髋关节置换术前规划系统,通过将多张全髋关节二维横断面图像堆叠成全髋关节三维块状图,基于三维卷积核结构的卷积神经网络,对全髋关节三维块状图进行识别,有效的提取到每张全髋关节二维横断面图像的股骨区域,并根据提取得到的股骨区域进行三维建模,以根据该股骨区域三维模型,在进行全髋关节置换术前规划时,提高全髋关节识别精度。The preoperative planning system for total hip replacement based on deep learning provided by this application stacks multiple two-dimensional cross-sectional images of the total hip joint into a three-dimensional block diagram of the total hip joint, and is based on a convolutional neural network with a three-dimensional convolution kernel structure. Identify the three-dimensional block diagram of the total hip joint, effectively extract the femoral region of each two-dimensional cross-sectional image of the total hip joint, and perform three-dimensional modeling according to the extracted femoral region, so that according to the three-dimensional model of the femoral region, in the Improve the accuracy of total hip joint recognition when performing preoperative planning for total hip arthroplasty.
在上述实施例的基础上,所述系统还包括训练模块,用于将预设三维块状图像集作为训练样本,输入到初始的三维分割神经网络进行训练,得到训练好的三维分割神经网络;其中,所述预设三维块状图像集中包含有若干张由预设二维横断面图像进行堆叠得到的预设三维块状图像。On the basis of the above embodiments, the system also includes a training module, which is used to input the preset three-dimensional block image set as a training sample into the initial three-dimensional segmentation neural network for training, and obtain a trained three-dimensional segmentation neural network; Wherein, the set of preset three-dimensional block images includes several preset three-dimensional block images obtained by stacking preset two-dimensional cross-sectional images.
在上述实施例的基础上,所述训练模块包括样本二维横断面图像获取 单元、第一标注单元、块状图第一构建单元和第一训练单元,其中,样本二维横断面图像获取单元用于获取若干张样本髋关节二维横断面图像;第一标注单元用于对每一张样本髋关节二维横断面图像中的样本股骨区域进行标注,并对所述样本股骨区域的股骨头像素标注股骨头区域标签,得到若干张第一预设髋关节二维横断面图像;块状图第一构建单元用于根据样本髋关节二维横断面的采集顺序,将若干张所述第一预设髋关节二维横断面图像进行堆叠,得到对应的第一预设三维块状图像;第一训练单元用于将多个所述第一预设三维块状图像输入到所述初始的三维分割神经网络进行训练,得到训练好的三维分割神经网络;On the basis of the above embodiments, the training module includes a sample two-dimensional cross-sectional image acquisition unit, a first labeling unit, a block diagram first construction unit and a first training unit, wherein the sample two-dimensional cross-sectional image acquisition unit It is used to obtain several two-dimensional cross-sectional images of the sample hip joint; the first labeling unit is used to mark the sample femoral area in each sample hip joint two-dimensional cross-sectional image, and to mark the femoral head of the sample femoral area The pixels are labeled with the femoral head area label to obtain several first preset two-dimensional cross-sectional images of the hip joint; The preset two-dimensional cross-sectional images of the hip joint are stacked to obtain a corresponding first preset three-dimensional block image; the first training unit is configured to input a plurality of the first preset three-dimensional block images into the initial three-dimensional Segment the neural network for training to obtain a trained three-dimensional segmentation neural network;
其中,所述初始的三维分割神经网络是由U-Net卷积网络构建得到的,且所述U-Net卷积网络的卷积核为三维卷积核。Wherein, the initial three-dimensional segmentation neural network is constructed by the U-Net convolution network, and the convolution kernel of the U-Net convolution network is a three-dimensional convolution kernel.
在上述实施例的基础上,所述系统还包括第二标注单元、块状图第二构建单元和第二训练单元,其中,第二标注单元用于对若干张所述第一预设髋关节二维横断面图像中样本股骨区域的骨皮质像素标注骨皮质区域标签,得到若干张第二预设髋关节二维横断面图像;块状图第二构建单元用于将若干张所述第二预设髋关节二维横断面图像按顺序进行堆叠,得到对应的第二预设三维块状图;第二训练单元用于通过多个所述第二预设三维块状图,对所述训练好的三维分割神经网络的参数进行优化,得到髋关节识别模型。On the basis of the above-mentioned embodiments, the system further includes a second labeling unit, a second construction unit of the block diagram, and a second training unit, wherein the second labeling unit is used to identify several sheets of the first preset hip joint In the two-dimensional cross-sectional image, the cortical bone pixels of the sample femoral region are labeled with the cortical bone region label to obtain several second preset two-dimensional cross-sectional images of the hip joint; The preset two-dimensional cross-sectional images of the hip joint are stacked in order to obtain the corresponding second preset three-dimensional block diagram; the second training unit is used to train the training The parameters of a good three-dimensional segmentation neural network are optimized to obtain a hip joint recognition model.
在上述实施例的基础上,所述系统还包括骨皮质区域识别模块、髓腔区域确定模块、髓腔解剖轴线确定模块、股骨颈干角计算模块和股骨柄假体确定模块,其中,骨皮质区域识别模块用于若所述待识别的全髋关节三维块状图为髋关节三维块状图,将所述髋关节三维块状图输入到所述髋关节识别模型,得到每张髋关节二维横断面图像中的股骨区域和骨皮质区域;髓腔区域确定模块用于根据所述股骨区域和所述骨皮质区域,确定髓腔区域;髓腔解剖轴线确定模块用于计算所述髓腔区域中每个髓腔层面的中点坐标,并根据所述中点坐标,将所有中心点进行直线拟合,确定髓腔解剖轴线;股骨颈干角计算模块用于根据所述髓腔解剖轴线和股骨颈轴线,计算得到股骨颈干角的角度值;股骨柄假体确定模块用于根据所述角度值、所述髓腔区域和股骨头旋转中心位置,确定股骨柄假体模型的类型和放置 位置。On the basis of the above embodiments, the system also includes a cortical bone area identification module, a medullary cavity area determination module, a medullary cavity anatomical axis determination module, a femoral neck-shaft angle calculation module and a femoral stem prosthesis determination module, wherein the cortical bone The area recognition module is used to input the three-dimensional block diagram of the hip joint into the hip joint recognition model if the three-dimensional block diagram of the total hip joint to be identified is a three-dimensional block diagram of the hip joint, and obtain two dimensions of each hip joint. The femoral region and the cortical bone region in the three-dimensional cross-sectional image; the medullary cavity region determination module is used to determine the medullary cavity region according to the femoral region and the cortical bone region; the medullary cavity anatomical axis determination module is used to calculate the medullary cavity The midpoint coordinates of each medullary cavity level in the region, and according to the midpoint coordinates, all the central points are fitted with a straight line to determine the anatomical axis of the medullary cavity; the femoral neck shaft angle calculation module is used to calculate the anatomical axis of the medullary cavity and femoral neck axis to calculate the angle value of the femoral neck shaft angle; the femoral stem prosthesis determination module is used to determine the type of femoral stem prosthesis model and Placement.
在上述实施例的基础上,所述系统还包括股骨头中心点像素坐标计算模块、坐标转换模块、股骨头旋转中心确定模块和髋臼杯假体确定模块,其中,骨头中心点像素坐标计算模块用于基于平面图像的质心公式,获取所述三维图像中股骨区域的股骨头中心点的像素坐标;坐标转换模块用于将所述像素坐标转换为图像坐标;股骨头旋转中心确定模块用于确定股骨头旋转中心位置;髋臼杯假体确定模块用于根据股骨头旋转中心位置,获取第一尺寸信息,以根据所述第一尺寸信息确定第二尺寸信息,其中,所述第一尺寸信息为股骨头对应的尺寸信息,所述第二尺寸信息为髋臼杯假体模型对应的尺寸信息。On the basis of the foregoing embodiments, the system also includes a femoral head center point pixel coordinate calculation module, a coordinate conversion module, a femoral head rotation center determination module, and an acetabular cup prosthesis determination module, wherein the bone center point pixel coordinate calculation module It is used to obtain the pixel coordinates of the center point of the femoral head in the femoral region in the three-dimensional image based on the centroid formula of the planar image; the coordinate transformation module is used to convert the pixel coordinates into image coordinates; the femoral head rotation center determination module is used to determine The position of the center of rotation of the femoral head; the acetabular cup prosthesis determination module is used to obtain the first size information according to the position of the center of rotation of the femoral head, so as to determine the second size information according to the first size information, wherein the first size information is the size information corresponding to the femoral head, and the second size information is the size information corresponding to the acetabular cup prosthesis model.
在上述实施例的基础上,所述术前规划系统还包括:位置矫正模块,所述位置矫正模块用于对放置的股骨柄假体模型的放置位置或者髋臼杯假体模型的放置位置进行矫正,以使所述股骨柄假体模型的放置位置以及所述髋臼杯假体模型的放置位置满足预设位置要求。On the basis of the above-mentioned embodiments, the preoperative planning system also includes: a position correction module, which is used to adjust the placement position of the placed femoral stem prosthesis model or the placement position of the acetabular cup prosthesis model Correction, so that the placement position of the femoral stem prosthesis model and the placement position of the acetabular cup prosthesis model meet the preset position requirements.
在上述实施例的基础上,所述全髋关节三维图像构建模块,在根据股骨头旋转中心位置,获取第一尺寸信息时,具体用于:On the basis of the above-mentioned embodiments, the three-dimensional image construction module of the total hip joint is specifically used for:
确定股骨头旋转中心位置与股骨头区域边缘交点;Determine the intersection point between the center of rotation of the femoral head and the edge of the femoral head area;
通过交点与股骨头中心位置之间的长度值,获取第一尺寸信息。The first size information is obtained through the length value between the intersection point and the center position of the femoral head.
在上述实施例的基础上,所述全髋关节三维图像构建模块,在确定股骨头旋转中心位置与股骨头区域边缘交点,通过交点与股骨头中心位置之间的长度值,获取第一尺寸信息时,具体用于;On the basis of the above-described embodiments, the three-dimensional image construction module of the total hip joint, after determining the intersection point between the center of rotation of the femoral head and the edge of the femoral head region, obtains the first size information through the length value between the intersection point and the center position of the femoral head when, specifically for;
确定股骨头旋转中心位置与股骨头区域边缘在不同角度下的多个交点;Determining multiple intersections between the position of the center of rotation of the femoral head and the edge of the femoral head area at different angles;
通过多个交点与股骨头中心位置之间的多个长度值的平均值或中值,获取第一尺寸信息。The first dimension information is obtained by using the average value or the median value of the multiple length values between the multiple intersection points and the central position of the femoral head.
在上述实施例的基础上,所述全髋关节识别模块,在根据所述股骨区域和所述骨皮质区域,确定髓腔区域时,具体用于:On the basis of the above-mentioned embodiments, the total hip joint recognition module is specifically used for:
将所述股骨区域中骨皮质区域滤除,得到髓腔区域。The cortical bone area in the femoral area was filtered out to obtain the medullary cavity area.
本申请提供的系统是用于执行上述各方法实施例的,具体流程和详细内容请参照上述实施例,此处不再赘述。The system provided in the present application is used to execute the above method embodiments. For the specific process and details, please refer to the above embodiments, which will not be repeated here.
图10为本申请提供的电子设备的结构示意图,如图10所示,该电子设备可以包括:处理器(Processor)1001、通信接口(Communications Interface)1002、存储器(Memory)1003和通信总线1004,其中,处理器1001,通信接口1002,存储器1003通过通信总线1004完成相互间的通信。处理器1001可以调用存储器1003中的逻辑指令,以执行基于深度学习的全髋关节置换术前影像规划方法。Fig. 10 is a schematic structural diagram of an electronic device provided by the present application. As shown in Fig. 10, the electronic device may include: a processor (Processor) 1001, a communication interface (Communications Interface) 1002, a memory (Memory) 1003 and a communication bus 1004, Wherein, the processor 1001 , the communication interface 1002 , and the memory 1003 communicate with each other through the communication bus 1004 . The processor 1001 can call the logic instructions in the memory 1003 to execute a deep learning-based image planning method for preoperative total hip arthroplasty.
一种示例,所述处理器1001,用于获取待识别的全髋关节三维块状图,所述待识别的全髋关节三维块状图是由多张全髋关节二维横断面图像堆叠而成的;将所述待识别的全髋关节三维块状图输入到训练好的三维分割神经网络中,得到所述训练好的三维分割神经网络输出的每张全髋关节二维横断面图像中的股骨区域,所述训练好的三维分割神经网络是由标注有股骨区域的标签的预设三维块状图作为训练样本,对卷积神经网络进行训练得到的;基于三维重建技术,根据每张全髋关节二维横断面图像中的股骨区域,得到所述股骨区域的三维图像。An example, the processor 1001 is configured to obtain a three-dimensional block diagram of the total hip joint to be identified, and the three-dimensional block diagram of the total hip joint to be identified is obtained by stacking multiple two-dimensional cross-sectional images of the total hip joint The three-dimensional block diagram of the total hip joint to be identified is input into the trained three-dimensional segmentation neural network, and each two-dimensional cross-sectional image of the total hip joint output by the trained three-dimensional segmentation neural network is obtained. femoral region, the trained three-dimensional segmentation neural network is obtained by training the convolutional neural network using the preset three-dimensional block diagram marked with the label of the femoral region as a training sample; based on the three-dimensional reconstruction technology, according to each The femoral region in the two-dimensional cross-sectional image of the total hip joint is obtained to obtain a three-dimensional image of the femoral region.
此外,上述的存储器1003中的逻辑指令可以通过软件功能单元的形式实现并作为独立的产品销售或使用时,可以存储在一个计算机可读取存储介质中。基于这样的理解,本申请的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的部分可以以软件产品的形式体现出来,该计算机软件产品存储在一个存储介质中,包括若干指令用以使得一台计算机设备(可以是个人计算机,服务器,或者网络设备等)执行本申请各个实施例所述方法的全部或部分步骤。而前述的存储介质包括:U盘、移动硬盘、只读存储器(ROM,Read-Only Memory)、随机存取存储器(RAM,Random Access Memory)、磁碟或者光盘等各种可以存储程序代码的介质。In addition, the above logic instructions in the memory 1003 may be implemented in the form of software function units and when sold or used as an independent product, may be stored in a computer-readable storage medium. Based on this understanding, the technical solution of the present application is essentially or the part that contributes to the prior art or the part of the technical solution can be embodied in the form of a software product, and the computer software product is stored in a storage medium, including Several instructions are used to make a computer device (which may be a personal computer, a server, or a network device, etc.) execute all or part of the steps of the methods described in the various embodiments of the present application. The aforementioned storage media include: U disk, mobile hard disk, read-only memory (ROM, Read-Only Memory), random access memory (RAM, Random Access Memory), magnetic disk or optical disc, etc., which can store program codes. .
另一方面,本申请还提供一种计算机程序产品,所述计算机程序产品包括存储在非暂态计算机可读存储介质上的计算机程序,所述计算机程序包括程序指令,当所述程序指令被计算机执行时,计算机能够执行上述各方法所提供的基于深度学习的全髋关节置换术前规划方法。On the other hand, the present application also provides a computer program product, the computer program product includes a computer program stored on a non-transitory computer-readable storage medium, the computer program includes program instructions, and when the program instructions are executed by a computer During execution, the computer can execute the deep learning-based preoperative planning method for total hip replacement provided by the above methods.
一种示例,获取待识别的全髋关节三维块状图,所述待识别的全髋关节三维块状图是由多张全髋关节二维横断面图像堆叠而成的;将所述待识别的全髋关节三维块状图输入到训练好的三维分割神经网络中,得到所述训练好的三维分割神经网络输出的每张全髋关节二维横断面图像中的股骨区域,所述训练好的三维分割神经网络是由标注有股骨区域的标签的预设三维块状图作为训练样本,对卷积神经网络进行训练得到的;基于三维重建技术,根据每张全髋关节二维横断面图像中的股骨区域,得到所述股骨区域的三维图像。An example is to obtain the three-dimensional block diagram of the total hip joint to be identified, and the three-dimensional block diagram of the total hip joint to be identified is formed by stacking multiple two-dimensional cross-sectional images of the total hip joint; The three-dimensional block diagram of the total hip joint is input into the trained three-dimensional segmentation neural network to obtain the femur region in each two-dimensional cross-sectional image of the total hip joint output by the three-dimensional segmentation neural network trained well. The 3D segmentation neural network is obtained by training the convolutional neural network with the preset 3D block diagram marked with the label of the femoral region as the training sample; based on the 3D reconstruction technology, according to each 2D cross-sectional image of the total hip joint In the femoral region, a three-dimensional image of the femoral region is obtained.
又一方面,本申请还提供一种非暂态计算机可读存储介质,其上存储有计算机程序,该计算机程序被处理器执行时实现以执行上述各实施例提供的基于深度学习的全髋关节置换术前规划方法。In another aspect, the present application also provides a non-transitory computer-readable storage medium, on which a computer program is stored, and when the computer program is executed by a processor, it is implemented to perform the deep learning-based total hip joint provided by the above-mentioned embodiments. Preoperative planning methods for replacement surgery.
一种示例,获取待识别的全髋关节三维块状图,所述待识别的全髋关节三维块状图是由多张全髋关节二维横断面图像堆叠而成的;将所述待识别的全髋关节三维块状图输入到训练好的三维分割神经网络中,得到所述训练好的三维分割神经网络输出的每张全髋关节二维横断面图像中的股骨区域,所述训练好的三维分割神经网络是由标注有股骨区域的标签的预设三维块状图作为训练样本,对卷积神经网络进行训练得到的;基于三维重建技术,根据每张全髋关节二维横断面图像中的股骨区域,得到所述股骨区域的三维图像。An example is to obtain the three-dimensional block diagram of the total hip joint to be identified, and the three-dimensional block diagram of the total hip joint to be identified is formed by stacking multiple two-dimensional cross-sectional images of the total hip joint; The three-dimensional block diagram of the total hip joint is input into the trained three-dimensional segmentation neural network to obtain the femur region in each two-dimensional cross-sectional image of the total hip joint output by the three-dimensional segmentation neural network trained well. The 3D segmentation neural network is obtained by training the convolutional neural network with the preset 3D block diagram marked with the label of the femoral region as the training sample; based on the 3D reconstruction technology, according to each 2D cross-sectional image of the total hip joint In the femoral region, a three-dimensional image of the femoral region is obtained.
以上所描述的装置实施例仅仅是示意性的,其中所述作为分离部件说明的单元可以是或者也可以不是物理上分开的,作为单元显示的部件可以是或者也可以不是物理单元,即可以位于一个地方,或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部模块来实现本实施例方案的目的。本领域普通技术人员在不付出创造性的劳动的情况下,即可以理解并实施。The device embodiments described above are only illustrative, and the units described as separate components may or may not be physically separated, and the components shown as units may or may not be physical units, that is, they may be located in One place, or it can be distributed to multiple network elements. Part or all of the modules can be selected according to actual needs to achieve the purpose of the solution of this embodiment. It can be understood and implemented by those skilled in the art without any creative effort.
通过以上的实施方式的描述,本领域的技术人员可以清楚地了解到各实施方式可借助软件加必需的通用硬件平台的方式来实现,当然也可以通过硬件。基于这样的理解,上述技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来,该计算机软件产品可以存储在 计算机可读存储介质中,如ROM/RAM、磁碟、光盘等,包括若干指令用以使得一台计算机设备(可以是个人计算机,服务器,或者网络设备等)执行各个实施例或者实施例的某些部分所述的方法。Through the above description of the implementations, those skilled in the art can clearly understand that each implementation can be implemented by means of software plus a necessary general-purpose hardware platform, and of course also by hardware. Based on this understanding, the essence of the above technical solution or the part that contributes to the prior art can be embodied in the form of software products, and the computer software products can be stored in computer-readable storage media, such as ROM/RAM, magnetic discs, optical discs, etc., including several instructions to make a computer device (which may be a personal computer, server, or network device, etc.) execute the methods described in various embodiments or some parts of the embodiments.
最后应说明的是:以上实施例仅用以说明本申请的技术方案,而非对其限制;尽管参照前述实施例对本申请进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本申请各实施例技术方案的精神和范围。。Finally, it should be noted that: the above embodiments are only used to illustrate the technical solutions of the present application, rather than limiting them; although the present application has been described in detail with reference to the foregoing embodiments, those of ordinary skill in the art should understand that: it can still Modifications are made to the technical solutions described in the foregoing embodiments, or equivalent replacements are made to some of the technical features; and these modifications or replacements do not make the essence of the corresponding technical solutions deviate from the spirit and scope of the technical solutions of the various embodiments of the present application. .

Claims (10)

  1. 一种基于深度学习的全髋关节置换术前规划系统,包括:A preoperative planning system for total hip replacement based on deep learning, including:
    全髋关节图像采集模块,用于获取待识别的全髋关节三维块状图,所述待识别的全髋关节三维块状图是由多张全髋关节二维横断面图像堆叠而成的;A total hip joint image acquisition module, configured to acquire a three-dimensional block diagram of the total hip joint to be identified, wherein the three-dimensional block diagram of the total hip joint to be identified is formed by stacking multiple two-dimensional cross-sectional images of the total hip joint;
    全髋关节识别模块,用于将所述待识别的全髋关节三维块状图输入到训练好的三维分割神经网络中,获取所述训练好的三维分割神经网络输出的每张全髋关节二维横断面图像中的股骨区域,所述训练好的三维分割神经网络是由标注有股骨区域的标签的预设三维块状图像作为训练样本,对卷积神经网络进行训练得到的;The total hip joint identification module is used to input the three-dimensional block diagram of the total hip joint to be identified into the trained three-dimensional segmentation neural network, and obtain the two images of each total hip joint output by the trained three-dimensional segmentation neural network. The femur region in the three-dimensional cross-sectional image, the trained three-dimensional segmentation neural network is obtained by training the convolutional neural network by using the preset three-dimensional block image marked with the label of the femoral region as a training sample;
    全髋关节三维图像构建模块,用于基于三维重建技术,根据每张全髋关节二维横断面图像中的股骨区域,得到所述股骨区域的三维图像。The three-dimensional image construction module of the total hip joint is used to obtain a three-dimensional image of the femoral region based on the three-dimensional reconstruction technology based on the femoral region in each two-dimensional cross-sectional image of the total hip joint.
  2. 根据权利要求1所述的基于深度学习的全髋关节置换术前规划系统,所述全髋关节识别模块,还用于:According to the deep learning-based preoperative planning system for total hip replacement according to claim 1, the total hip recognition module is also used for:
    获取若干张样本髋关节二维横断面图像;Obtain several two-dimensional cross-sectional images of the sample hip joint;
    对每一张样本髋关节二维横断面图像中的样本股骨区域进行标注,并对所述样本股骨区域的股骨头像素标注股骨头区域标签,得到若干张第一预设髋关节二维横断面图像;Annotate the sample femoral region in each sample hip joint two-dimensional cross-sectional image, and mark the femoral head region label on the femoral head pixel of the sample femoral region, and obtain several first preset two-dimensional cross-sectional hip joints image;
    根据样本髋关节二维横断面的采集顺序,将若干张所述第一预设髋关节二维横断面图像进行堆叠,得到对应的第一预设三维块状图像;According to the acquisition order of the two-dimensional cross-section of the sample hip joint, stacking several first preset two-dimensional cross-sectional images of the hip joint to obtain a corresponding first preset three-dimensional block image;
    将多个所述第一预设三维块状图像输入到初始的三维分割神经网络进行训练,得到训练好的三维分割神经网络;Inputting a plurality of the first preset three-dimensional block images into the initial three-dimensional segmentation neural network for training to obtain a trained three-dimensional segmentation neural network;
    其中,所述初始的三维分割神经网络是由U-Net卷积网络构建得到的,且所述U-Net卷积网络的卷积核为三维卷积核。Wherein, the initial three-dimensional segmentation neural network is constructed by the U-Net convolution network, and the convolution kernel of the U-Net convolution network is a three-dimensional convolution kernel.
  3. 根据权利要求2所述的基于深度学习的全髋关节置换术前规划系统,所述全髋关节识别模块,在所述将多个所述第一预设三维块状图像输入到初始的三维分割神经网络进行训练,得到训练好的三维分割神经网络之后,还用于:According to the preoperative planning system for total hip arthroplasty based on deep learning according to claim 2, the total hip joint recognition module, when inputting a plurality of the first preset three-dimensional block images into the initial three-dimensional segmentation The neural network is trained, and after the trained three-dimensional segmentation neural network is obtained, it is also used for:
    对若干张所述第一预设髋关节二维横断面图像中样本股骨区域的骨皮质像素标注骨皮质区域标签,得到若干张第二预设髋关节二维横断面图像;Marking the cortical bone region labels on the cortical bone pixels of the sample femur region in the first preset two-dimensional cross-sectional images of the hip joint to obtain several second preset two-dimensional cross-sectional images of the hip joint;
    将若干张所述第二预设髋关节二维横断面图像按顺序进行堆叠,得到对 应的第二预设三维块状图;Stacking several second preset hip joint two-dimensional cross-sectional images in order to obtain a corresponding second preset three-dimensional block diagram;
    通过多个所述第二预设三维块状图,对所述训练好的三维分割神经网络的参数进行优化,得到髋关节识别模型。A hip joint recognition model is obtained by optimizing the parameters of the trained three-dimensional segmentation neural network through a plurality of second preset three-dimensional block diagrams.
  4. 根据权利要求3所述的基于深度学习的全髋关节置换术前规划系统,所述全髋关节识别模块,在所述获取待识别的关节三维块状图之后,还用于:According to the preoperative planning system for total hip arthroplasty based on deep learning according to claim 3, the total hip joint identification module, after acquiring the three-dimensional block diagram of the joint to be identified, is further used for:
    将所述髋关节三维块状图输入到所述髋关节识别模型,得到每张髋关节二维横断面图像中的股骨区域和骨皮质区域;Inputting the three-dimensional block diagram of the hip joint into the hip joint recognition model to obtain the femoral region and the cortical bone region in each two-dimensional cross-sectional image of the hip joint;
    根据所述股骨区域和所述骨皮质区域,确定髓腔区域;determining the medullary cavity area according to the femoral area and the cortical bone area;
    计算所述髓腔区域中每个髓腔层面的中点坐标,并根据所述中点坐标,将所有中心点进行直线拟合,确定髓腔解剖轴线;calculating the midpoint coordinates of each medullary canal level in the medullary cavity region, and performing straight line fitting on all the midpoints according to the midpoint coordinates to determine the anatomical axis of the medullary cavity;
    根据所述髓腔解剖轴线和股骨颈轴线,计算得到股骨颈干角的角度值;Calculate the angle value of the femoral neck-shaft angle according to the anatomical axis of the medullary cavity and the femoral neck axis;
    根据所述角度值、所述髓腔区域和股骨头旋转中心位置,确定股骨柄假体模型的类型和放置位置。According to the angle value, the medullary cavity area and the position of the center of rotation of the femoral head, the type and placement position of the femoral stem prosthesis model are determined.
  5. 根据权利要求1至4任一项所述的基于深度学习的全髋关节置换术前规划系统,所述全髋关节三维图像构建模块,在所述基于三维重建技术,根据每张全髋关节二维横断面图像中的股骨区域,得到所述股骨区域的三维图像之后,还用于:According to the preoperative planning system for total hip joint replacement based on deep learning according to any one of claims 1 to 4, the three-dimensional image construction module of the total hip joint, based on the three-dimensional reconstruction technology, according to two images of each total hip joint The femur region in the three-dimensional cross-sectional image, after obtaining the three-dimensional image of the femur region, is also used for:
    基于平面图像的质心公式,获取所述三维图像中股骨区域的股骨头中心点的像素坐标;Based on the centroid formula of the planar image, obtain the pixel coordinates of the center point of the femoral head in the femoral region in the three-dimensional image;
    将所述像素坐标转换为图像坐标;converting said pixel coordinates to image coordinates;
    确定股骨头旋转中心位置;Determine the position of the center of rotation of the femoral head;
    根据股骨头旋转中心位置,获取第一尺寸信息,以根据所述第一尺寸信息确定第二尺寸信息,其中,所述第一尺寸信息为股骨头对应的尺寸信息,所述第二尺寸信息为髋臼杯假体模型对应的尺寸信息。According to the rotation center position of the femoral head, the first size information is obtained to determine the second size information according to the first size information, wherein the first size information is the size information corresponding to the femoral head, and the second size information is Size information corresponding to the acetabular cup prosthesis model.
  6. 根据权利要求5所述的基于深度学习的全髋关节置换术前规划系统,所述术前规划系统还包括:位置矫正模块,所述位置矫正模块用于对放置的股骨柄假体模型的放置位置或者髋臼杯假体模型的放置位置进行矫正,以使所述股骨柄假体模型的放置位置以及所述髋臼杯假体模型的放置位置满足预设位置要求。According to the preoperative planning system for total hip arthroplasty based on deep learning according to claim 5, the preoperative planning system also includes: a position correction module, which is used to place the placed femoral stem prosthesis model The position or the placement position of the acetabular cup prosthesis model is corrected, so that the placement position of the femoral stem prosthesis model and the placement position of the acetabular cup prosthesis model meet the preset position requirements.
  7. 根据权利要求5所述的基于深度学习的全髋关节置换术前规划系统, 其中,所述全髋关节三维图像构建模块,在根据股骨头旋转中心位置,获取第一尺寸信息时,具体用于:The preoperative planning system for total hip arthroplasty based on deep learning according to claim 5, wherein the three-dimensional image construction module of the total hip joint is specifically used to obtain the first size information according to the position of the center of rotation of the femoral head. :
    确定股骨头旋转中心位置与股骨头区域边缘交点;Determine the intersection point between the center of rotation of the femoral head and the edge of the femoral head area;
    通过交点与股骨头中心位置之间的长度值,获取第一尺寸信息。The first size information is obtained through the length value between the intersection point and the center position of the femoral head.
  8. 根据权利要求7所述的基于深度学习的全髋关节置换术前规划系统,其中,所述全髋关节三维图像构建模块,在确定股骨头旋转中心位置与股骨头区域边缘交点,通过交点与股骨头中心位置之间的长度值,获取第一尺寸信息时,具体用于;The preoperative planning system for total hip arthroplasty based on deep learning according to claim 7, wherein the three-dimensional image construction module of the total hip joint determines the intersection point between the center of rotation of the femoral head and the edge of the femoral head area, and passes the intersection point and the femoral head area edge. The length value between the center positions of the bones is specifically used when obtaining the first size information;
    确定股骨头旋转中心位置与股骨头区域边缘在不同角度下的多个交点;Determining multiple intersections between the position of the center of rotation of the femoral head and the edge of the femoral head area at different angles;
    通过多个交点与股骨头中心位置之间的多个长度值的平均值或中值,获取第一尺寸信息。The first dimension information is obtained by using the average value or the median value of the multiple length values between the multiple intersection points and the central position of the femoral head.
  9. 根据权利要求2所述的基于深度学习的全髋关节置换术前规划系统,其中,所述卷积神经网络中两个连续的三维卷积核是通过残差结构连接的,且所述卷积神经网络的损失函数是由DICE loss和BCE loss组成的。The preoperative planning system for total hip arthroplasty based on deep learning according to claim 2, wherein two continuous three-dimensional convolution kernels in the convolutional neural network are connected through a residual structure, and the convolution The loss function of the neural network is composed of DICE loss and BCE loss.
  10. 根据权利要求4所述的基于深度学习的全髋关节置换术前规划系统,其中,所述全髋关节识别模块,在根据所述股骨区域和所述骨皮质区域,确定髓腔区域时,具体用于:The preoperative planning system for total hip arthroplasty based on deep learning according to claim 4, wherein the total hip joint identification module, when determining the medullary cavity area according to the femoral area and the cortical bone area, specifically Used for:
    将所述股骨区域中骨皮质区域滤除,得到髓腔区域。The cortical bone area in the femoral area was filtered out to obtain the medullary cavity area.
PCT/CN2023/070788 2022-01-27 2023-01-05 Total hip replacement preoperative planning system based on deep learning WO2023142956A1 (en)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
CN202210101412.9 2022-01-27
CN202210101412.9A CN114419618B (en) 2022-01-27 2022-01-27 Total hip replacement preoperative planning system based on deep learning

Publications (1)

Publication Number Publication Date
WO2023142956A1 true WO2023142956A1 (en) 2023-08-03

Family

ID=81279449

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/CN2023/070788 WO2023142956A1 (en) 2022-01-27 2023-01-05 Total hip replacement preoperative planning system based on deep learning

Country Status (2)

Country Link
CN (1) CN114419618B (en)
WO (1) WO2023142956A1 (en)

Families Citing this family (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114419618B (en) * 2022-01-27 2024-02-02 北京长木谷医疗科技股份有限公司 Total hip replacement preoperative planning system based on deep learning
CN115393272B (en) * 2022-07-15 2023-04-18 北京长木谷医疗科技有限公司 Knee joint patella replacement three-dimensional preoperative planning system and method based on deep learning
CN115252233B (en) * 2022-07-28 2024-04-26 大连理工大学 Automatic planning method for acetabular cup in total hip arthroplasty based on deep learning
CN116071372B (en) * 2022-12-30 2024-03-19 北京长木谷医疗科技股份有限公司 Knee joint segmentation method, knee joint segmentation device, electronic equipment and storage medium
CN116128912A (en) * 2023-02-06 2023-05-16 北京纳通医用机器人科技有限公司 Method, device, equipment and medium for segmenting hip joint image
CN116597002B (en) * 2023-05-12 2024-01-30 北京长木谷医疗科技股份有限公司 Automatic femoral stem placement method, device and equipment based on deep reinforcement learning
CN116650110B (en) * 2023-06-12 2024-05-07 北京长木谷医疗科技股份有限公司 Automatic knee joint prosthesis placement method and device based on deep reinforcement learning
CN116883326A (en) * 2023-06-21 2023-10-13 北京长木谷医疗科技股份有限公司 Knee joint anatomical site recognition method, device, equipment and readable storage medium

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110648337A (en) * 2019-09-23 2020-01-03 武汉联影医疗科技有限公司 Hip joint segmentation method, hip joint segmentation device, electronic apparatus, and storage medium
CN112150473A (en) * 2020-09-24 2020-12-29 北京羽医甘蓝信息技术有限公司 Three-dimensional jaw bone image segmentation modeling method and device based on CT and terminal equipment
CN113689402A (en) * 2021-08-24 2021-11-23 北京长木谷医疗科技有限公司 Deep learning-based femoral medullary cavity form identification method, device and storage medium
CN114419618A (en) * 2022-01-27 2022-04-29 北京长木谷医疗科技有限公司 Deep learning-based preoperative planning system for total hip replacement

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111179350B (en) * 2020-02-13 2022-04-08 张逸凌 Hip joint image processing system
CN111888059B (en) * 2020-07-06 2021-07-27 北京长木谷医疗科技有限公司 Full hip joint image processing method and device based on deep learning and X-ray
CN111938813A (en) * 2020-07-21 2020-11-17 南京市第一医院 Preoperative planning method for hip revision surgery
CN113017829B (en) * 2020-08-22 2023-08-29 张逸凌 Preoperative planning method, system, medium and device for total knee arthroplasty based on deep learning
CN112971981B (en) * 2021-03-02 2022-02-08 北京长木谷医疗科技有限公司 Deep learning-based total hip joint image processing method and equipment
CN113658142B (en) * 2021-08-19 2024-03-12 江苏金马扬名信息技术股份有限公司 Hip joint femur near-end segmentation method based on improved U-Net neural network

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110648337A (en) * 2019-09-23 2020-01-03 武汉联影医疗科技有限公司 Hip joint segmentation method, hip joint segmentation device, electronic apparatus, and storage medium
CN112150473A (en) * 2020-09-24 2020-12-29 北京羽医甘蓝信息技术有限公司 Three-dimensional jaw bone image segmentation modeling method and device based on CT and terminal equipment
CN113689402A (en) * 2021-08-24 2021-11-23 北京长木谷医疗科技有限公司 Deep learning-based femoral medullary cavity form identification method, device and storage medium
CN114419618A (en) * 2022-01-27 2022-04-29 北京长木谷医疗科技有限公司 Deep learning-based preoperative planning system for total hip replacement

Also Published As

Publication number Publication date
CN114419618A (en) 2022-04-29
CN114419618B (en) 2024-02-02

Similar Documents

Publication Publication Date Title
WO2023142956A1 (en) Total hip replacement preoperative planning system based on deep learning
WO2022170768A1 (en) Unicondylar joint image processing method and apparatus, device, and storage medium
EP2996599B1 (en) Planning methods for surgical correction of abnormal bones
AU2020101836A4 (en) A method for generating femoral x-ray films based on deep learning and digital reconstruction of radiological image
Tian et al. Efficient computer-aided design of dental inlay restoration: a deep adversarial framework
CN112037200A (en) Method for automatically identifying anatomical features and reconstructing model in medical image
CN111652888B (en) Method and device for determining medullary cavity anatomical axis based on deep learning
WO2023160272A1 (en) Deep learning-based hip replacement postoperative image evaluation method and system
CN108309334B (en) Data processing method of spine X-ray image
WO2024011943A1 (en) Deep learning-based knee joint patella resurfacing three-dimensional preoperative planning method and system
WO2024001140A1 (en) Vertebral body sub-region segmentation method and apparatus, and storage medium
Kim et al. 3D reconstruction of leg bones from X-ray images using CNN-based feature analysis
CN112785609B (en) CBCT tooth segmentation method based on deep learning
Kodym et al. Skull shape reconstruction using cascaded convolutional networks
WO2023165265A1 (en) Dynamic assessment method and apparatus for knee joint tension balance
Wang et al. A DCNN system based on an iterative method for automatic landmark detection in cephalometric X-ray images
CN114792326A (en) Surgical navigation point cloud segmentation and registration method based on structured light
CN106991720B (en) Personalized acetabulum reconstruction steel plate pre-bending method based on finite element analysis and operation
CN114119872A (en) Method for analyzing 3D printing intraspinal plants based on artificial intelligence big data
CN113077499A (en) Pelvis registration method, pelvis registration device and pelvis registration system
Liu et al. Tracking-based deep learning method for temporomandibular joint segmentation
CN115252233B (en) Automatic planning method for acetabular cup in total hip arthroplasty based on deep learning
Li et al. Malocclusion treatment planning via pointnet based spatial transformation network
Xie et al. Automatic Individual Tooth Segmentation in Cone-Beam Computed Tomography Based on Multi-Task CNN and Watershed Transform
Liu A system for computational analysis and reconstruction of 3D comminuted bone fractures

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 23745841

Country of ref document: EP

Kind code of ref document: A1