WO2020077961A1 - Image-based breast lesion identification method and device - Google Patents

Image-based breast lesion identification method and device Download PDF

Info

Publication number
WO2020077961A1
WO2020077961A1 PCT/CN2019/082687 CN2019082687W WO2020077961A1 WO 2020077961 A1 WO2020077961 A1 WO 2020077961A1 CN 2019082687 W CN2019082687 W CN 2019082687W WO 2020077961 A1 WO2020077961 A1 WO 2020077961A1
Authority
WO
WIPO (PCT)
Prior art keywords
breast
image
feature
convolution
images
Prior art date
Application number
PCT/CN2019/082687
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 WO2020077961A1 publication Critical patent/WO2020077961A1/en

Links

Images

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
    • A61B6/52Devices using data or image processing specially adapted for radiation diagnosis
    • A61B6/5211Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data
    • A61B6/5217Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data extracting a diagnostic or physiological parameter from medical diagnostic data
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
    • A61B6/50Clinical applications
    • A61B6/502Clinical applications involving diagnosis of breast, i.e. mammography

Definitions

  • Embodiments of the present invention relate to the technical field of machine learning, and in particular, to a method and device for identifying breast imaging lesions.
  • breast imaging can use low-dose X-rays to examine human breasts. It can detect various breast tumors, cysts and other lesions, which helps to detect breast cancer early and reduce its mortality.
  • Breast imaging is an effective detection method that can be used to diagnose a variety of female breast-related diseases. Of course, the most important use is breast cancer, especially early breast cancer screening. Therefore, if you can effectively detect the early manifestations of various breast cancers on the breast image, it will be of great help to the doctor.
  • Embodiments of the present invention provide a method and device for identifying breast imaging lesions, which are used to solve the problem of low efficiency of the method for judging breast lesions in breast imaging based on doctor experience in the prior art.
  • Embodiments of the present invention provide a method for identifying breast imaging lesions, including:
  • the breast image is input into a feature extraction module to obtain feature images of different sizes of the breast image;
  • the feature extraction module includes N convolution modules;
  • the N convolution modules are down-sampling convolution blocks and / or Or up-sampling convolution block; the size of the feature image extracted by each down-sampling convolution block or up-sampling convolution block is different, each of the N convolution modules includes a first convolution layer, The second convolutional layer; the number of feature images output by the first convolution layer is less than the number of feature images input by the first convolution layer; the number of feature images output by the second convolution layer Is greater than the number of feature images output by the first convolutional layer; N is greater than 0;
  • the breast lesion of the breast image is determined.
  • the acquiring feature images of different sizes of the breast image includes:
  • the first feature image output from the N / 2th down-sampling convolution block is used to sequentially extract N / 2 second feature images of the mammary gland image through the N / 2 up-sampling convolution block, each up-sampling convolution block
  • the sizes of the extracted second feature images are different;
  • N feature images of different sizes of the breast images are determined.
  • a possible implementation manner before the feature processing module, further includes a feature preprocessing module; before the input of the breast image to the feature extraction module, further includes:
  • the feature preprocessing module includes a convolution layer, a BN layer, a Relu layer and a pooling layer; the convolution kernel of the feature preprocessing module The size is larger than the size of the convolution kernel in the N convolution modules;
  • the feature pre-processing module includes multiple consecutive convolutional layers, a BN layer, a Relu layer, and a pooling layer; the size of the convolution kernel of the feature pre-processing module and the N convolution modules The size of the largest convolution kernel is equal.
  • the method before inputting the breast image to the feature extraction module, the method further includes:
  • the mammary gland image according to the picture format corresponding to the at least one set of window width and window level is used as the mammary gland image input to the feature extraction module.
  • the mammary gland image includes mammary gland images with different projection positions of different breasts; the inputting the mammary gland image to a feature extraction module includes:
  • the breast lesion identification frame is determined from the feature image; including:
  • the first breast lesion identification frame is deleted.
  • An embodiment of the present invention provides a device for identifying breast imaging lesions, including:
  • Acquisition unit for acquiring mammary gland image
  • the feature extraction module includes N convolution modules;
  • the N convolution modules are for downsampling Convolution block or up-sampling convolution block; the size of the feature image extracted by each down-sampling convolution block or up-sampling convolution block is different, and each of the N convolution modules includes a first volume Multilayer, second convolutional layer; the number of feature images output by the first convolutional layer is less than the number of feature images input by the first convolutional layer; the feature images output by the second convolutional layer The number of is greater than the number of feature images output by the first convolutional layer; N is greater than 0; for any one of the feature images of different sizes of the breast image, the breast is determined from the feature image Lesion recognition frame; according to the breast lesion recognition frame determined from each characteristic image, determine the breast lesion of the breast image.
  • the processing unit is specifically used to:
  • the breast image is sequentially passed through N / 2 down-sampling convolution blocks to extract N / 2 first feature images of the breast image; the first feature images output from the N / 2 down-sampling convolution block are sequentially passed N / 2 up-sampling convolution blocks extract N / 2 second feature images of the mammography image, the size of the second feature images extracted by each up-sampling convolution block are different; the first feature images of the same size After merging with the second feature image, N feature images of different sizes of the breast images are determined.
  • the mammary gland image includes mammary gland images with different projection positions of different breasts; the processing unit is specifically used for:
  • the feature extraction module Take the breast image of the other breast of the same projection position of the breast image as the reference image of the breast image, and input it to the feature extraction module to obtain a reference feature image; determine the first breast in the feature image The lesion identification frame and the second breast lesion identification frame in the reference feature image; if it is determined that the positions and / or sizes of the first breast lesion identification frame and the second breast lesion identification frame are the same, delete the The first breast lesion identification frame.
  • an embodiment of the present invention provides a computer device including at least one processing unit and at least one storage unit, wherein the storage unit stores a computer program, and when the program is executed by the processing unit, The processing unit executes the steps of the method for identifying breast imaging lesions.
  • an embodiment of the present invention provides a computer-readable storage medium that stores a computer program executable by a computer device, and when the program runs on the computer device, causes the computer device to execute the mammary gland Steps of the method of image lesion identification.
  • an embodiment of the present invention further provides a computer program product
  • the computer program product includes a computer program stored on a computer-readable storage medium
  • the computer program includes program instructions, and when the program instructions are When the device is executed, the computer device is caused to perform the steps of the method for identifying breast imaging lesions.
  • the lesion of the mammary gland can be quickly identified, and the efficiency of identifying the breast lesion is improved.
  • the number of channels output by the first convolution layer is reduced, and the number of channels output by the second convolution layer is increased, so that the effective information in the image is effectively retained during the convolution process. While reducing the amount of parameters, the effectiveness of feature image extraction is improved, thereby improving the accuracy of breast lesion recognition in breast images.
  • FIG. 1a is a schematic diagram of a breast image provided by an embodiment of the present invention.
  • FIG. 1b is a schematic diagram of a breast image provided by an embodiment of the present invention.
  • 1c is a schematic diagram of a breast image provided by an embodiment of the present invention.
  • 1d is a schematic diagram of a breast image provided by an embodiment of the present invention.
  • FIG. 2 is a schematic flowchart of a method for identifying a breast imaging lesion according to an embodiment of the present invention
  • 3a is a schematic structural diagram of a feature extraction module provided by an embodiment of the present invention.
  • 3b is a schematic structural diagram of a feature extraction module provided by an embodiment of the present invention.
  • 3c is a schematic structural diagram of a feature extraction module provided by an embodiment of the present invention.
  • FIG. 3 is a schematic flowchart of a breast imaging lesion recognition provided by an embodiment of the present invention.
  • FIG. 4 is a schematic flowchart of a breast imaging lesion recognition provided by an embodiment of the present invention.
  • FIG. 5 is a schematic structural diagram of a device for identifying breast imaging lesions according to an embodiment of the present invention.
  • FIG. 6 is a schematic structural diagram of a computer device according to an embodiment of the present invention.
  • the breast X-ray image is taken as an example for an exemplary description, and other images will not be repeated here.
  • Breast X-ray images can be used to examine the breasts of humans (mainly women) using low-dose (about 0.7 mSv) X-rays. It can detect various breast tumors, cysts and other lesions, which helps to detect breast cancer early. And reduce its mortality. Some countries encourage older women (generally over 45 years old) to perform mammography regularly (with intervals ranging from one year to five years) to screen for early breast cancer.
  • the mammary gland image generally includes four X-ray images, which are four mammary gland images of the two projection positions of the two breasts (head and tail CC, medial and lateral oblique MLO), as shown in Figure 1a, Figure 1b, Figure 1c, 1d.
  • the prior art often only detects a single type of lesions such as calcifications or masses, and cannot simultaneously detect multiple lesions, and the application range is narrow. At the same time, for many types of lesions such as calcifications, masses, asymmetry, and structural distortion, the accuracy of detection is poor and cannot meet the application requirements.
  • an embodiment of the present invention provides a method for identifying breast imaging lesions, as shown in FIG. 2, including:
  • Step 201 Obtain a breast image
  • Step 202 Input the mammary gland image into a feature extraction module to obtain feature images of different sizes of the mammary gland image;
  • the feature extraction module includes N convolution modules; the N convolution modules are down-sampling convolution blocks and / or up-sampling convolution blocks; each down-sampling convolution block or up-sampling convolution block extraction The size of the feature images of the two are different.
  • Each of the N convolution modules includes a first convolution layer and a second convolution layer; the number of feature images output by the first convolution layer is less than The number of feature images input by the first convolution layer; the number of feature images output by the second convolution layer is greater than the number of feature images output by the first convolution layer; N is greater than 0;
  • the feature extraction module may include three down-sampling convolution blocks.
  • Each convolution module may include a first convolution layer and a second convolution layer.
  • the first convolution layer includes a convolution layer, a normalization (BN) layer connected to the convolution layer, and a connection to the BN layer
  • the activation function layer of Fig. 3a includes a first convolution layer and a second convolution layer.
  • the step of the feature image passing through the convolution module may include:
  • Step 1 input the feature image input by the convolution module to the first convolution layer to obtain the first feature image;
  • the convolution kernel of the first convolution layer may be N1 * m * m * N2;
  • N1 is Describe the number of channels of the feature image input by the convolution module,
  • N2 is the number of channels of the first feature image; N1> N2;
  • Step 2 Input the first feature image into the second convolution layer to obtain the second feature image;
  • the convolution kernel of the first convolution layer may be N2 * m * m * N3;
  • N3 is the channel of the second feature image Number; N3> N2;
  • Step 3 After combining the feature image input by the convolution module and the second feature image, it is determined as the feature image output by the convolution module.
  • the method for determining the feature image corresponding to the breast image described above is only one possible implementation manner. In other possible implementation manners, the feature image corresponding to the breast image may also be determined by other methods, which is not specifically limited.
  • the activation function in the embodiment of the present invention may be multiple types of activation functions, for example, it may be a linear rectification function (Rectified Linear Unit, ReLU), which is not specifically limited;
  • ReLU Rectified Linear Unit
  • the feature extraction module in the embodiment of the present invention may be a feature extraction module in a (2Dimensions, 2D) convolutional neural network.
  • the first convolution layer The size of the convolution kernel can be m * m, and the size of the second convolution layer can be n * n; m and n can be the same or different, which is not limited here; where, m, n is greater than or An integer equal to 1.
  • the number of feature images output by the first convolution layer is less than the number of feature images input by the first convolution layer; the number of feature images output by the second convolution layer is greater than the first convolution layer The number of output feature images.
  • a possible implementation manner as shown in FIG. 3c, further includes a third convolution layer between the first convolution layer and the second convolution layer;
  • the feature image input by the three convolution layers is the image output by the first convolution layer, and the feature image output by the third convolution layer is the image input by the second convolution layer.
  • the size of the convolution kernel of the third convolutional layer may be k * k, and k may be the same as m or n, or may be different, which is not limited herein.
  • the size of the convolution kernel of the first convolution layer is 3 * 3; the size of the convolution kernel of the second convolution layer is 3 * 3; the third convolution layer The size of the convolution kernel is 1 * 1.
  • the perception field of feature extraction can be effectively improved, which is beneficial to improve the accuracy of breast lesion recognition.
  • the feature images of different sizes may be feature images of different pixels, for example, the feature image with pixels 500 ⁇ 500 and the feature image with pixels 1000 ⁇ 1000 are feature images with different sizes.
  • a pre-trained breast lesion detection model is used to extract feature images of different sizes of breast images.
  • the model is determined after training a plurality of labeled breast images using a 2D convolutional neural network.
  • the image is scaled to a specific size so that the scale of the pixels in each direction is the same as the actual length.
  • the feature extraction module includes N / 2 down-sampling convolution blocks and N / 2 up-sampling convolution blocks; and acquiring feature images of different sizes of the breast image includes:
  • the first feature image output from the N / 2th down-sampling convolution block is used to sequentially extract N / 2 second feature images of the mammary gland image through the N / 2 up-sampling convolution block, each up-sampling convolution block
  • the sizes of the extracted second feature images are different;
  • N feature images of different sizes of the breast images are determined.
  • the feature extraction module also includes a feature preprocessing module before; the feature preprocessing module includes a convolution layer and a BN layer, One Relu layer and one pooling layer; the size of the convolution kernel of the feature preprocessing module is larger than that of any of the N convolution modules.
  • the size of the convolution kernel of the convolution layer may be 7 * 7, and the interval is 2 pixels.
  • the pooling layer is 2 * 2 maximum pooling.
  • the feature preprocessing module includes a plurality of continuous convolutional layers, a BN layer, a Relu layer, and a pooling layer; the size of the convolution kernel of the feature preprocessing module and the N The largest convolution kernel in each convolution module has the same size.
  • the step of the feature image passing through the feature preprocessing module may include: inputting the breast image to the feature preprocessing module to obtain a preprocessed feature image; and using the preprocessed feature image as an input of the feature extraction module.
  • Step 203 For any one of the feature images of different sizes of the breast image, determine a breast lesion recognition frame from the feature image.
  • a pre-trained breast lesion detection model is used to determine the breast lesion recognition frame from the feature image.
  • the breast lesion detection model is determined after training multiple breast images of the marked breast lesion using a 2D convolutional neural network. .
  • the area framed by the breast lesion identification frame determined from the feature image does not necessarily contain breast lesions, so each breast lesion identification frame needs to be screened according to the breast lesion probability of the breast lesion identification frame, and the breast lesion probability is less than the preset threshold
  • the breast lesion identification frame is deleted, where the breast lesion probability is the probability that the area framed by the breast lesion identification frame is the breast lesion.
  • Step 204 Determine the breast lesion of the breast image according to the breast lesion identification frame determined from each feature image.
  • the recognition frame is output as the breast lesion in the breast image
  • the output breast lesion parameters include the central coordinates of the breast lesion and the diameter of the breast lesion, wherein the central coordinates of the breast lesion are the breast lesion identification
  • the center coordinate of the frame, the diameter of the breast lesion is the distance from the center of the breast lesion identification frame to one of the faces.
  • both large-sized breast lesions and small-sized breast lesions can be detected, which improves the detection of breast lesions Precision.
  • the method of automatically detecting the breast lesion in the present application effectively improves the recognition efficiency of the breast lesion.
  • the breast lesion identification frame determined from each feature image may have multiple identification frames corresponding to one breast lesion, if the number of breast lesions in the breast image is directly determined according to the number of breast lesion identification frames, the number of detected breast lesions will result There is a large deviation, so it is necessary to convert each feature image into a feature image of the same size and align it, and then screen the breast lesion recognition frame determined from each feature image, and determine the screened breast lesion recognition frame as the breast Breast lesions in the image.
  • the breast image includes breast images of different breasts with different projection positions;
  • the input of the breast image to the feature extraction module includes:
  • the breast lesion identification frame is determined from the feature image; including:
  • the first breast lesion identification frame is deleted.
  • the screening process of the breast lesion identification frame includes the following steps, as shown in Figure 3:
  • Step 301 Determine the breast lesion identification frame with the highest probability of breast lesion from the breast lesion identification frame of each feature image.
  • Step 302 Calculate the intersection ratio of the breast lesion identification frame with the largest breast lesion probability and the other breast lesion identification frames.
  • Step 303 Delete the identification frame of other breast lesions whose cross-combination ratio is greater than a preset threshold.
  • Step 304 Determine the breast lesion identification frame with the highest probability of breast lesions from the remaining other breast lesion identification frames, and repeat the screening process of the breast lesion identification frame until there are no other remaining breast lesion identification frames.
  • the screening process of the breast lesion identification frame described above will be described in conjunction with specific examples, and the breast lesion identification frames identified in each feature image are set to A, B, C, D, E, and F, and the above breast lesion identification frames
  • the above breast lesion recognition frame is sorted according to the probability of breast lesion from large to small: E, C, A, B, D, F, after sorting, we can know the breast lesion with the highest probability of breast lesion in the breast lesion recognition frame of each feature image
  • the recognition frame is E, and then calculate the intersection ratio IOU between the breast lesion recognition frame E and each other breast lesion recognition frame.
  • the calculation method of the intersection ratio is shown in equation (1):
  • m is the breast lesion identification frame with the highest probability of breast lesions
  • n is the breast lesion identification frame compared with the breast lesion identification frame m
  • IOU is the intersection ratio between the breast lesion identification frame m and the breast lesion identification frame n.
  • the preset threshold to 0.5, if the intersection ratio between the breast lesion identification frame C and the breast lesion identification frame E is greater than 0.5, the intersection ratio between the breast lesion identification frame A and the breast lesion identification frame E is greater than 0.5, the breast If the intersection ratio of the lesion identification frame B, the breast lesion identification frame D, the breast lesion identification frame F and the breast lesion identification frame E is less than 0.5, delete the breast lesion identification frame C and the breast lesion identification frame A, and identify the breast lesion Box E is identified as a breast lesion in the breast image.
  • the remaining other breast lesion recognition frames B, D, F are sorted according to the breast lesion probability, and the breast lesion recognition frame with the highest breast lesion probability is determined as the breast lesion recognition frame B, and then the breast lesion recognition frame B and the breast lesion are calculated The cross ratio between the recognition frames D and the cross ratio between the breast lesion recognition frame B and the breast lesion recognition frame F. If the intersection ratio between the breast lesion identification frame B and the breast lesion identification frame D is greater than 0.5, and the intersection ratio between the breast lesion identification frame B and the breast lesion identification frame F is less than 0.5, then delete the breast lesion identification frame D, replace The breast lesion recognition frame B and the breast lesion recognition frame F are determined as breast lesions in the breast image.
  • the breast lesion identification frame identified in each feature image is screened based on the breast lesion probability of the breast lesion identification frame and the intersection and comparison between the breast lesion identification frames, avoiding repeated detection and output of the same breast lesion in the breast image, Improve the accuracy of detecting the number of breast lesions in breast images.
  • Step 401 Obtain a breast image as a training sample.
  • the acquired multiple breast images can be used directly as training samples, or the acquired multiple breast images can be enhanced to expand the data volume of the training samples. Enhancement operations include, but are not limited to: randomly setting pixels up, down, left, and right (Such as 0-20 pixels), random rotation setting angle (such as -15-15 degrees), random zoom setting multiple (such as 0.85-1.15 times).
  • Step 402 Manually mark the breast lesions in the training sample.
  • the breast lesions in the training sample can be marked by doctors and other professionals, and the content of the marking includes the central coordinates of the breast lesions and the diameter of the breast lesions. Specifically, multiple doctors can mark the breast lesions, and determine the final breast lesions and the parameters of the breast lesions through a multiple-vote synthesis method, and the results are saved in the form of a mask.
  • the manual labeling of the breast lesions in the training sample and the training sample are in no particular order of enhancement. You can manually mark the breast lesions in the training sample, and then the training sample of the labeled breast lesions can be enhanced. The training samples are enhanced, and then the training samples after the enhancement operations are manually marked.
  • step 403 the training samples are input to the convolutional neural network for training to determine the breast lesion recognition model.
  • the structure of the convolutional neural network includes an input layer, a down-sampling convolution block, an up-sampling convolution block, a target detection network, and an output layer. Pre-process the training samples and input them into the convolutional neural network, calculate the loss function of the output breast lesions and the mask image of the pre-labeled training samples, and then iterate iteratively using the back propagation algorithm and the sgd optimization algorithm to determine the breast lesions Detection model.
  • the process of extracting feature images of different sizes of breast images using the breast lesion detection model determined by the above training includes the following steps:
  • the mammary gland image is successively passed through N / 2 down-sampling convolution blocks to extract the first feature images of the N mammary gland images.
  • the size of the first feature image extracted by each down-sampling convolution block is different, and N / 2 is greater than 0.
  • the down-sampling convolution block includes a first convolution layer and a second convolution layer, a group connection layer, a front-back connection layer, and a down-sampling layer.
  • Step 2 The first feature image output from the N / 2 down-sampling convolution block is sequentially used to extract the second feature image of N / 2 mammary gland images through the N / 2 up-sampling convolution block.
  • the size of the second feature image extracted by each up-sampling convolution block is different.
  • the up-sampling convolution block includes a convolution layer, a group connection layer, a front-back connection layer, an up-sampling layer, and a synthesis connection layer.
  • Convolution layer includes convolution operation, batch normalization layer and RELU layer.
  • step three after combining the first feature image and the second feature image with the same size, feature images of different sizes of N / 2 breast images are determined.
  • the first feature image and the second feature image of the same size are combined through the up-sampling convolution block to determine feature images of different sizes.
  • the number of channels of the first feature image and the second feature image are combined, and the size of the feature image obtained after the merging is the same as the size of the first feature image and the second feature image.
  • the process of determining the breast lesion recognition frame from the feature image using the breast lesion detection model determined by the above training includes the following steps:
  • Step 1 For any pixel in the feature image, with the pixel as the center, diffuse to the surroundings to determine the first area.
  • Step 2 Set multiple preset frames in the first area according to preset rules.
  • the preset frame can be set to various shapes.
  • the preset rule may be that the center of the preset frame coincides with the center of the first area, or that the corner of the preset frame coincides with the angle of the first area, and so on.
  • the way to select the preset frame of the breast lesion is that, for each pixel of each feature map, it is considered as an anchor point. Set multiple preset frames with different aspect ratios on each anchor point.
  • the convolution of the feature map predicts a coordinate and size offset and confidence, and the preset frame is determined based on the coordinate and size offset and confidence.
  • Step 3 For any preset frame, predict the position deviation of the preset frame from the first area.
  • Step 4 Adjust the preset frame according to the position deviation to determine the breast lesion identification frame, and predict the breast lesion probability of the breast lesion identification frame.
  • the breast lesion probability is the probability that the area selected by the breast lesion identification frame is the breast lesion.
  • the specific training process may include: inputting the training data image to the above-mentioned convolutional neural network for calculation.
  • multiple images of different window widths and positions of the lesion are introduced.
  • the prediction frame set with the highest confidence and the prediction frame set with the maximum coincidence with the training sample are selected.
  • the cross-entropy of the confidence of the prediction frame and the labeling of the sample and the cross-entropy of the labeled lesion of the training sample and the offset of the prediction frame are used as the loss function.
  • the training optimization algorithm uses the sgd algorithm with momentum and step attenuation.
  • the input image is preprocessed to improve the effect of feature extraction.
  • the acquiring breast image includes:
  • Step 1 Determine the binary image of the breast image according to Gaussian filtering
  • Step 2 Obtain the connected area of the binarized image, and use the area with the largest area in the connected area corresponding to the breast image as the segmented breast image;
  • Step 3 Add the segmented breast image to a preset image template to generate a pre-processed breast image; and use the pre-processed breast image as the breast image input to the feature extraction module.
  • the input of the preprocessing module is a breast image saved in Dicom format.
  • Preprocessing can include gland segmentation and image normalization; the main purpose of gland segmentation is to extract the breast part of the input breast image to remove other unrelated interference images; image normalization is to normalize the image into Unified format images, specifically, include:
  • the specific binarized threshold can be obtained by finding the maximum class interval of the grayscale histogram of the image.
  • the binarized result can be obtained by flooding to obtain independent regional blocks, and the area of each regional block is counted; the area on the image corresponding to the largest regional block is used as Segmented breast image.
  • the preset image template may be a square image of a black bottom plate; specifically, the obtained divided breast image may be expanded into a 1: 1 square image by adding a black border.
  • the output breast image can be scaled by pixels, for example, the image difference can be scaled to 4096 pixels ⁇ 4096 pixels.
  • the window width and position of the mammary gland can be adjusted to obtain a better identification effect of breast lesion identification.
  • the method before inputting the breast image to the feature extraction module, the method further includes:
  • the mammary gland image according to the picture format corresponding to the at least one set of window width and window level is used as the mammary gland image input to the feature extraction module.
  • the dicom image can be converted into a png image through three sets of window width and window levels.
  • the first set of window width is 4000 and the window level is 2000; the second set of window width is 1000; the window level is 2000 ;
  • the third group has a window width of 1500 and a window level of 1500.
  • an embodiment of the present invention provides a device for identifying breast lesions. As shown in FIG. 5, the device can perform the flow of a method for identifying breast lesions.
  • the device includes an acquiring unit 501 and a processing unit 502.
  • the obtaining unit 501 is used to obtain breast images
  • the processing unit 502 is configured to input the breast image into a feature extraction module to obtain feature images of different sizes of the breast image;
  • the feature extraction module includes N convolution modules;
  • the N convolution modules are as follows Sampling convolutional blocks and / or upsampling convolutional blocks; the size of the feature image extracted by each downsampling convolutional block or upsampling convolutional block is different, and each of the N convolutional modules includes The first convolution layer and the second convolution layer; the number of feature images output by the first convolution layer is less than the number of feature images input by the first convolution layer; the second convolution layer outputs The number of feature images of is greater than the number of feature images output by the first convolutional layer; N is greater than 0; for any one of the feature images of different sizes of the breast image, from the feature image Identify the breast lesion recognition frame; determine the breast lesion of the breast image according to the breast lesion recognition frame determined from each feature image.
  • processing unit 502 is specifically configured to:
  • the breast image is sequentially passed through N / 2 down-sampling convolution blocks to extract N / 2 first feature images of the breast image; the first feature images output from the N / 2 down-sampling convolution block are sequentially passed N / 2 up-sampling convolution blocks extract N / 2 second feature images of the mammography image, the size of the second feature images extracted by each up-sampling convolution block are different; the first feature images of the same size After merging with the second feature image, N feature images of different sizes of the breast images are determined.
  • the processing unit 502 is specifically configured to:
  • the feature preprocessing module includes a convolution layer, a BN layer, a Relu layer and a pooling layer; the convolution kernel of the feature preprocessing module The size is larger than the size of the convolution kernel in the N convolution modules;
  • the feature pre-processing module includes multiple consecutive convolutional layers, a BN layer, a Relu layer, and a pooling layer; the size of the convolution kernel of the feature pre-processing module and the N convolution modules The size of the largest convolution kernel is equal.
  • the obtaining unit 501 is configured to:
  • the processing unit 502 is specifically used for:
  • the mammary gland image includes mammary gland images with different projection positions of different breasts; the acquiring unit 501 is configured to:
  • the feature extraction module Take the breast image of the other breast of the same projection position of the breast image as the reference image of the breast image, and input it to the feature extraction module to obtain a reference feature image; determine the first breast in the feature image The lesion identification frame and the second breast lesion identification frame in the reference feature image; if it is determined that the positions and / or sizes of the first breast lesion identification frame and the second breast lesion identification frame are the same, delete the The first breast lesion identification frame.
  • an embodiment of the present invention provides a computer device including at least one processing unit and at least one storage unit, wherein the storage unit stores a computer program when the program is executed by the processing unit So that the processing unit executes the steps of the method for identifying breast lesions.
  • FIG. 6 it is a schematic diagram of the hardware structure of the computer device described in the embodiment of the present invention.
  • the computer device may specifically be a desktop computer, a portable computer, a smart phone, a tablet computer, or the like.
  • the computer device may include a memory 801, a processor 802, and a computer program stored on the memory.
  • the memory 801 may include a read-only memory (ROM) and a random access memory (RAM), and provide the processor 802 with program instructions and data stored in the memory 801.
  • the computer equipment described in the embodiments of the present application may further include an input device 803 and an output device 804.
  • the input device 803 may include a keyboard, a mouse, a touch screen, etc .
  • the output device 804 may include a display device, such as a liquid crystal display (Liquid Crystal Display, LCD), a cathode ray tube (Cathode Ray Tube, CRT), a touch screen, and the like.
  • the memory 801, the processor 802, the input device 803, and the output device 804 may be connected through a bus or in other ways. In FIG. 6, connection through a bus is used as an example.
  • the processor 802 calls the program instructions stored in the memory 801 and executes the method for identifying breast lesions provided in the foregoing embodiments according to the obtained program instructions.
  • an embodiment of the present invention also provides a computer-readable storage medium that stores a computer program executable by a computer device, and when the program runs on the computer device, causes the computer device to execute a breast Steps of the method of lesion identification.
  • an embodiment of the present invention also provides a computer program product, the computer program product includes a computer program stored on a computer-readable storage medium, the computer program includes program instructions, when the program instructions When executed by a computer device, the computer device is caused to perform the steps of the method for identifying breast imaging lesions.
  • the embodiments of the present invention may be provided as methods or computer program products. Therefore, the present invention may take the form of an entirely hardware embodiment, an entirely software embodiment, or an embodiment combining software and hardware. Moreover, the present invention may take the form of a computer program product implemented on one or more computer usable storage media (including but not limited to disk storage, CD-ROM, optical storage, etc.) containing computer usable program code.
  • computer usable storage media including but not limited to disk storage, CD-ROM, optical storage, etc.
  • These computer program instructions may also be stored in a computer readable memory that can guide a computer or other programmable data processing device to work in a specific manner, so that the instructions stored in the computer readable memory produce an article of manufacture including an instruction device, the instructions The device implements the functions specified in one block or multiple blocks of the flowchart one flow or multiple flows and / or block diagrams.
  • These computer program instructions can also be loaded onto a computer or other programmable data processing device, so that a series of operating steps are performed on the computer or other programmable device to produce computer-implemented processing, which is executed on the computer or other programmable device
  • the instructions provide steps for implementing the functions specified in one block or multiple blocks of the flowchart one flow or multiple flows and / or block diagrams.

Abstract

An image-based breast lesion identification method and device, which relate to the field of machine learning technology. The method includes: acquiring a breast image (201), inputting the breast image into a feature extraction module to acquire feature images of different sizes of the breast image (202); with respect to any one of the feature images of different sizes of the breast image, determining a breast lesion identification frame (203) from that feature image; according to the breast lesion identification frame determined from each feature image, determining the breast lesion of the breast image ( 204).

Description

一种乳腺影像病灶识别的方法及装置Method and device for identifying breast imaging lesions
本申请要求在2018年10月16日提交中国专利局、申请号为201811201699.2、申请名称为“一种乳腺影像病灶识别的方法及装置”的中国专利申请的优先权,其全部内容通过引用结合在本申请中。This application requires the priority of the Chinese patent application filed on October 16, 2018 with the Chinese Patent Office, the application number is 201811201699.2, and the application name is "a method and device for breast imaging lesion recognition", the entire content of which is incorporated by reference in In this application.
技术领域Technical field
本发明实施例涉及机器学习技术领域,尤其涉及一种乳腺影像病灶识别的方法及装置。Embodiments of the present invention relate to the technical field of machine learning, and in particular, to a method and device for identifying breast imaging lesions.
背景技术Background technique
目前,乳腺影像可以利用低剂量的X光检查人类的乳房,它能侦测各种乳房肿瘤、囊肿等病灶,有助于早期发现乳癌,并降低其死亡率。乳腺影像是一种有效的检测方法,可以用于诊断多种女性乳腺相关的疾病。当然,其中最主要的使用还是在乳腺癌,尤其是早期乳腺癌的筛查上。因此若能有效的检测出乳腺影像上各种乳腺癌早期表现,对医生的帮助是巨大的。At present, breast imaging can use low-dose X-rays to examine human breasts. It can detect various breast tumors, cysts and other lesions, which helps to detect breast cancer early and reduce its mortality. Breast imaging is an effective detection method that can be used to diagnose a variety of female breast-related diseases. Of course, the most important use is breast cancer, especially early breast cancer screening. Therefore, if you can effectively detect the early manifestations of various breast cancers on the breast image, it will be of great help to the doctor.
当患者拍摄乳腺影像之后,医生通过个人经验判断乳腺影像中的病灶,该方法效率较低,并且存在较大的主观性。After the patient takes the breast image, the doctor judges the lesion in the breast image through personal experience. This method is inefficient and subjective.
发明内容Summary of the invention
本发明实施例提供一种乳腺影像病灶识别的方法及装置,用于解决现有技术中通过医生经验判断乳腺影像中乳腺病灶的方法效率低的问题。Embodiments of the present invention provide a method and device for identifying breast imaging lesions, which are used to solve the problem of low efficiency of the method for judging breast lesions in breast imaging based on doctor experience in the prior art.
本发明实施例提供了一种乳腺影像病灶识别的方法,包括:Embodiments of the present invention provide a method for identifying breast imaging lesions, including:
获取乳腺影像;Get breast images;
将所述乳腺影像输入至特征提取模块中,获取所述乳腺影像不同尺寸的特征图像;所述特征提取模块包括N个卷积模块;所述N个卷积模块为下采样卷积块和/或上采样卷积块;每个下采样卷积块或上采样卷积块提取的特征 图像的尺寸均不同,所述N个卷积模块的每个卷积模块中包括第一卷积层、第二卷积层;所述第一卷积层输出的特征图像的个数小于所述第一卷积层输入的特征图像的个数;所述第二卷积层输出的特征图像的个数大于所述第一卷积层输出的特征图像的个数;N大于0;The breast image is input into a feature extraction module to obtain feature images of different sizes of the breast image; the feature extraction module includes N convolution modules; the N convolution modules are down-sampling convolution blocks and / or Or up-sampling convolution block; the size of the feature image extracted by each down-sampling convolution block or up-sampling convolution block is different, each of the N convolution modules includes a first convolution layer, The second convolutional layer; the number of feature images output by the first convolution layer is less than the number of feature images input by the first convolution layer; the number of feature images output by the second convolution layer Is greater than the number of feature images output by the first convolutional layer; N is greater than 0;
针对所述乳腺影像的不同尺寸的特征图像中的任意一个特征图像,从所述特征图像中确定出乳腺病灶识别框;For any one of the feature images of different sizes of the breast image, determine the breast lesion identification frame from the feature image;
根据从各特征图像中确定出的乳腺病灶识别框,确定乳腺影像的乳腺病灶。According to the breast lesion identification frame identified from each feature image, the breast lesion of the breast image is determined.
一种可能的实现方式,所述获取所述乳腺影像的不同尺寸的特征图像,包括:A possible implementation manner, the acquiring feature images of different sizes of the breast image includes:
将所述乳腺影像依次通过N/2个下采样卷积块提取N/2个所述乳腺影像的第一特征图像;Successively extracting the first feature image of the breast image through N / 2 down-sampling convolution blocks through the breast image;
将第N/2个下采样卷积块输出的第一特征图像依次通过N/2个上采样卷积块提取N/2个所述乳腺影像的第二特征图像,每个上采样卷积块提取的第二特征图像的尺寸均不同;The first feature image output from the N / 2th down-sampling convolution block is used to sequentially extract N / 2 second feature images of the mammary gland image through the N / 2 up-sampling convolution block, each up-sampling convolution block The sizes of the extracted second feature images are different;
将尺寸相同的第一特征图像和第二特征图像合并后,确定N个所述乳腺影像的不同尺寸的特征图像。After merging the first feature image and the second feature image of the same size, N feature images of different sizes of the breast images are determined.
一种可能的实现方式,所述特征处理模块之前还包括特征预处理模块;所述将所述乳腺影像输入至特征提取模块之前,还包括:A possible implementation manner, before the feature processing module, further includes a feature preprocessing module; before the input of the breast image to the feature extraction module, further includes:
将所述乳腺影像输入至所述特征预处理模块中,所述特征预处理模块包括一个卷积层,一个BN层,一个Relu层和一个池化层;所述特征预处理模块的卷积核大小大于所述N个卷积模块中的卷积核的大小;Input the breast image into the feature preprocessing module, the feature preprocessing module includes a convolution layer, a BN layer, a Relu layer and a pooling layer; the convolution kernel of the feature preprocessing module The size is larger than the size of the convolution kernel in the N convolution modules;
或者,所述特征预处理模块包括连续的多个卷积层,一个BN层,一个Relu层和一个池化层;所述特征预处理模块的卷积核大小与所述N个卷积模块中的最大的卷积核的大小相等。Alternatively, the feature pre-processing module includes multiple consecutive convolutional layers, a BN layer, a Relu layer, and a pooling layer; the size of the convolution kernel of the feature pre-processing module and the N convolution modules The size of the largest convolution kernel is equal.
一种可能的实现方式,所述将所述乳腺影像输入至特征提取模块之前,还包括:In a possible implementation manner, before inputting the breast image to the feature extraction module, the method further includes:
获取所述乳腺影像的原始文件;Obtain the original file of the breast image;
在所述乳腺影像的原始文件中选取至少一组窗宽窗位,并获取所述至少一组窗宽窗位对应的图片格式的乳腺影像;Selecting at least one set of window width and window levels from the original file of the breast image, and acquiring a breast image in a picture format corresponding to the at least one set of window width and window levels;
根据所述至少一组窗宽窗位对应的图片格式的乳腺影像,作为输入至所述特征提取模块的乳腺影像。The mammary gland image according to the picture format corresponding to the at least one set of window width and window level is used as the mammary gland image input to the feature extraction module.
一种可能的实现方式,所述乳腺影像包括不同侧乳房的不同投照位的乳腺影像;所述将所述乳腺影像输入至特征提取模块,包括:A possible implementation manner, the mammary gland image includes mammary gland images with different projection positions of different breasts; the inputting the mammary gland image to a feature extraction module includes:
将所述乳腺影像的同一投照位的另一侧乳房的乳腺影像作为所述乳腺影像的参考影像,输入至所述特征提取模块,获得参考特征图像;Taking the breast image of the other breast in the same projection position of the breast image as the reference image of the breast image, and inputting it to the feature extraction module to obtain a reference feature image;
所述针对所述乳腺影像的不同尺寸的特征图像中的任意一个特征图像,从所述特征图像中确定出乳腺病灶识别框;包括:Any one of the feature images of different sizes of the breast image, the breast lesion identification frame is determined from the feature image; including:
确定所述特征图像中的第一乳腺病灶识别框和所述参考特征图像中的第二乳腺病灶识别框;Determining a first breast lesion recognition frame in the feature image and a second breast lesion recognition frame in the reference feature image;
若确定所述第一乳腺病灶识别框和所述第二乳腺病灶识别框的位置和/或大小都相同,则删除所述第一乳腺病灶识别框。If it is determined that the positions and / or sizes of the first breast lesion identification frame and the second breast lesion identification frame are the same, the first breast lesion identification frame is deleted.
本发明实施例提供了一种乳腺影像病灶识别的装置,包括:An embodiment of the present invention provides a device for identifying breast imaging lesions, including:
获取单元,用于获取乳腺影像;Acquisition unit for acquiring mammary gland image;
处理单元,用于将所述乳腺影像输入至特征提取模块中,获取所述乳腺影像不同尺寸的特征图像;所述特征提取模块包括N个卷积模块;所述N个卷积模块为下采样卷积块或上采样卷积块;每个下采样卷积块或上采样卷积块提取的特征图像的尺寸均不同,所述N个卷积模块的每个卷积模块中包括第一卷积层、第二卷积层;所述第一卷积层输出的特征图像的个数小于所述第一卷积层输入的特征图像的个数;所述第二卷积层输出的特征图像的个数大于所述第一卷积层输出的特征图像的个数;N大于0;针对所述乳腺影像的不同尺寸的特征图像中的任意一个特征图像,从所述特征图像中确定出乳腺病灶识别框;根据从各特征图像中确定出的乳腺病灶识别框,确定乳腺影像的乳腺病灶。A processing unit for inputting the breast image into a feature extraction module to obtain feature images of different sizes of the breast image; the feature extraction module includes N convolution modules; the N convolution modules are for downsampling Convolution block or up-sampling convolution block; the size of the feature image extracted by each down-sampling convolution block or up-sampling convolution block is different, and each of the N convolution modules includes a first volume Multilayer, second convolutional layer; the number of feature images output by the first convolutional layer is less than the number of feature images input by the first convolutional layer; the feature images output by the second convolutional layer The number of is greater than the number of feature images output by the first convolutional layer; N is greater than 0; for any one of the feature images of different sizes of the breast image, the breast is determined from the feature image Lesion recognition frame; according to the breast lesion recognition frame determined from each characteristic image, determine the breast lesion of the breast image.
一种可能的实现方式,所述处理单元,具体用于:In a possible implementation manner, the processing unit is specifically used to:
将所述乳腺影像依次通过N/2个下采样卷积块提取N/2个所述乳腺影像的第一特征图像;将第N/2个下采样卷积块输出的第一特征图像依次通过N/2个上采样卷积块提取N/2个所述乳腺影像的第二特征图像,每个上采样卷积块提取的第二特征图像的尺寸均不同;将尺寸相同的第一特征图像和第二特征图像合并后,确定N个所述乳腺影像的不同尺寸的特征图像。The breast image is sequentially passed through N / 2 down-sampling convolution blocks to extract N / 2 first feature images of the breast image; the first feature images output from the N / 2 down-sampling convolution block are sequentially passed N / 2 up-sampling convolution blocks extract N / 2 second feature images of the mammography image, the size of the second feature images extracted by each up-sampling convolution block are different; the first feature images of the same size After merging with the second feature image, N feature images of different sizes of the breast images are determined.
一种可能的实现方式,所述乳腺影像包括不同侧乳房的不同投照位的乳腺影像;所述处理单元,具体用于:A possible implementation manner, the mammary gland image includes mammary gland images with different projection positions of different breasts; the processing unit is specifically used for:
将所述乳腺影像的同一投照位的另一侧乳房的乳腺影像作为所述乳腺影像的参考影像,输入至所述特征提取模块,获得参考特征图像;确定所述特征图像中的第一乳腺病灶识别框和所述参考特征图像中的第二乳腺病灶识别框;若确定所述第一乳腺病灶识别框和所述第二乳腺病灶识别框的位置和/或大小都相同,则删除所述第一乳腺病灶识别框。Take the breast image of the other breast of the same projection position of the breast image as the reference image of the breast image, and input it to the feature extraction module to obtain a reference feature image; determine the first breast in the feature image The lesion identification frame and the second breast lesion identification frame in the reference feature image; if it is determined that the positions and / or sizes of the first breast lesion identification frame and the second breast lesion identification frame are the same, delete the The first breast lesion identification frame.
另一方面,本发明实施例提供了一种计算机设备,包括至少一个处理单元以及至少一个存储单元,其中,所述存储单元存储有计算机程序,当所述程序被所述处理单元执行时,使得所述处理单元执行上述乳腺影像病灶识别的方法的步骤。On the other hand, an embodiment of the present invention provides a computer device including at least one processing unit and at least one storage unit, wherein the storage unit stores a computer program, and when the program is executed by the processing unit, The processing unit executes the steps of the method for identifying breast imaging lesions.
又一方面,本发明实施例提供了一种计算机可读存储介质,其存储有可由计算机设备执行的计算机程序,当所述程序在所述计算机设备上运行时,使得所述计算机设备执行上述乳腺影像病灶识别的方法的步骤。In still another aspect, an embodiment of the present invention provides a computer-readable storage medium that stores a computer program executable by a computer device, and when the program runs on the computer device, causes the computer device to execute the mammary gland Steps of the method of image lesion identification.
又一方面,本发明实施例还提供了一种计算机程序产品,所述计算机程序产品包括存储在计算机可读存储介质上的计算机程序,所述计算机程序包括程序指令,当所述程序指令被计算机设备执行时,使所述计算机设备执行上述乳腺影像病灶识别的方法的步骤。In still another aspect, an embodiment of the present invention further provides a computer program product, the computer program product includes a computer program stored on a computer-readable storage medium, the computer program includes program instructions, and when the program instructions are When the device is executed, the computer device is caused to perform the steps of the method for identifying breast imaging lesions.
本发明实施例中,由于提取乳腺影像的特征图像,并识别每一个特征图像中的乳腺,可以快速识别乳腺的病灶,提高了乳腺病灶识别的效率。另外,通过在特征提取模块中,设置第一卷积层输出的通道数减少,且第二卷积层 输出的通道数增加,使得卷积过程中,有效的保留了图像中的有效信息,在减少参数量的同时,提高了特征图像的提取的有效性,进而提高了乳腺影像中乳腺病灶识别的准确性。In the embodiment of the present invention, since the feature image of the mammary gland image is extracted and the mammary gland in each feature image is identified, the lesion of the mammary gland can be quickly identified, and the efficiency of identifying the breast lesion is improved. In addition, in the feature extraction module, the number of channels output by the first convolution layer is reduced, and the number of channels output by the second convolution layer is increased, so that the effective information in the image is effectively retained during the convolution process. While reducing the amount of parameters, the effectiveness of feature image extraction is improved, thereby improving the accuracy of breast lesion recognition in breast images.
附图说明BRIEF DESCRIPTION
为了更清楚地说明本发明实施例中的技术方案,下面将对实施例描述中所需要使用的附图作简要介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域的普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。In order to more clearly explain the technical solutions in the embodiments of the present invention, the drawings required in the description of the embodiments will be briefly introduced below. Obviously, the drawings in the following description are only some embodiments of the present invention. Those of ordinary skill in the art can obtain other drawings based on these drawings without paying any creative labor.
图1a为本发明实施例提供的一种乳腺影像的示意图;FIG. 1a is a schematic diagram of a breast image provided by an embodiment of the present invention;
图1b为本发明实施例提供的一种乳腺影像的示意图;1b is a schematic diagram of a breast image provided by an embodiment of the present invention;
图1c为本发明实施例提供的一种乳腺影像的示意图;1c is a schematic diagram of a breast image provided by an embodiment of the present invention;
图1d为本发明实施例提供的一种乳腺影像的示意图;1d is a schematic diagram of a breast image provided by an embodiment of the present invention;
图2为本发明实施例提供的一种乳腺影像病灶识别的方法的流程示意图;2 is a schematic flowchart of a method for identifying a breast imaging lesion according to an embodiment of the present invention;
图3a为本发明实施例提供的一种特征提取模块的结构示意图;3a is a schematic structural diagram of a feature extraction module provided by an embodiment of the present invention;
图3b为本发明实施例提供的一种特征提取模块的结构示意图;3b is a schematic structural diagram of a feature extraction module provided by an embodiment of the present invention;
图3c为本发明实施例提供的一种特征提取模块的结构示意图;3c is a schematic structural diagram of a feature extraction module provided by an embodiment of the present invention;
图3为本发明实施例提供的一种乳腺影像病灶识别的流程示意图;FIG. 3 is a schematic flowchart of a breast imaging lesion recognition provided by an embodiment of the present invention;
图4为本发明实施例提供的一种乳腺影像病灶识别的流程示意图;4 is a schematic flowchart of a breast imaging lesion recognition provided by an embodiment of the present invention;
图5为本发明实施例提供的一种乳腺影像病灶识别的装置的结构示意图;5 is a schematic structural diagram of a device for identifying breast imaging lesions according to an embodiment of the present invention;
图6为本发明实施例提供的一种计算机设备的结构示意图。6 is a schematic structural diagram of a computer device according to an embodiment of the present invention.
具体实施方式detailed description
为了使本发明的目的、技术方案及有益效果更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。In order to make the objectives, technical solutions and beneficial effects of the present invention clearer, the present invention will be further described in detail below in conjunction with the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are only used to explain the present invention, and are not intended to limit the present invention.
本发明实施例中,以乳腺X射线影像为例,进行示例性的描述,其他影 像在此不再赘述。乳腺X射线影像可以是利用低剂量(约为0.7毫西弗)的X光检查人类(主要是女性)的乳房,它能侦测各种乳房肿瘤、囊肿等病灶,有助于早期发现乳癌,并降低其死亡率。有一些国家提倡年长(一般为45周岁以上)的女性定期(间隔从一年到五年不等)进行乳腺摄影,以筛检出早期的乳腺癌。乳腺影像一般包含四份X光摄像,分别为2侧乳房的2种投照位(头尾位CC,内外侧斜位MLO)的四份乳腺影像,如图1a、图1b、图1c、图1d所示。In the embodiment of the present invention, the breast X-ray image is taken as an example for an exemplary description, and other images will not be repeated here. Breast X-ray images can be used to examine the breasts of humans (mainly women) using low-dose (about 0.7 mSv) X-rays. It can detect various breast tumors, cysts and other lesions, which helps to detect breast cancer early. And reduce its mortality. Some countries encourage older women (generally over 45 years old) to perform mammography regularly (with intervals ranging from one year to five years) to screen for early breast cancer. The mammary gland image generally includes four X-ray images, which are four mammary gland images of the two projection positions of the two breasts (head and tail CC, medial and lateral oblique MLO), as shown in Figure 1a, Figure 1b, Figure 1c, 1d.
现有技术往往只检测钙化或者肿块这样单独类型的病灶,不能同时对多种病灶同时进行检出,应用范围狭窄。同时针对钙化、肿块、不对称性、结构扭曲等多种类型病灶,检测的准确性较差,无法满足应用要求。The prior art often only detects a single type of lesions such as calcifications or masses, and cannot simultaneously detect multiple lesions, and the application range is narrow. At the same time, for many types of lesions such as calcifications, masses, asymmetry, and structural distortion, the accuracy of detection is poor and cannot meet the application requirements.
针对上述问题,本发明实施例提供一种乳腺影像病灶识别的方法,如图2所示,包括:In response to the above problems, an embodiment of the present invention provides a method for identifying breast imaging lesions, as shown in FIG. 2, including:
步骤201:获取乳腺影像;Step 201: Obtain a breast image;
步骤202:将所述乳腺影像输入至特征提取模块中,获取所述乳腺影像不同尺寸的特征图像;Step 202: Input the mammary gland image into a feature extraction module to obtain feature images of different sizes of the mammary gland image;
其中,所述特征提取模块包括N个卷积模块;所述N个卷积模块为下采样卷积块和/或上采样卷积块;每个下采样卷积块或上采样卷积块提取的特征图像的尺寸均不同,所述N个卷积模块的每个卷积模块中包括第一卷积层、第二卷积层;所述第一卷积层输出的特征图像的个数小于所述第一卷积层输入的特征图像的个数;所述第二卷积层输出的特征图像的个数大于所述第一卷积层输出的特征图像的个数;N大于0;Wherein, the feature extraction module includes N convolution modules; the N convolution modules are down-sampling convolution blocks and / or up-sampling convolution blocks; each down-sampling convolution block or up-sampling convolution block extraction The size of the feature images of the two are different. Each of the N convolution modules includes a first convolution layer and a second convolution layer; the number of feature images output by the first convolution layer is less than The number of feature images input by the first convolution layer; the number of feature images output by the second convolution layer is greater than the number of feature images output by the first convolution layer; N is greater than 0;
举例来说,该特征提取模块可以包括三个下采样卷积块。每个卷积模块可以包括第一卷积层和第二卷积层,第一卷积层包括卷积层,与卷积层连接的归一化(Batch Normalization,BN)层、与BN层连接的激活函数层,如图3a示出的卷积模块包括第一卷积层和第二卷积层。For example, the feature extraction module may include three down-sampling convolution blocks. Each convolution module may include a first convolution layer and a second convolution layer. The first convolution layer includes a convolution layer, a normalization (BN) layer connected to the convolution layer, and a connection to the BN layer The activation function layer of Fig. 3a includes a first convolution layer and a second convolution layer.
为增加特征提取模块的深度,如图3b所示,一种可能的实现方式,特征图像经过卷积模块的步骤可以包括:In order to increase the depth of the feature extraction module, as shown in FIG. 3b, a possible implementation manner, the step of the feature image passing through the convolution module may include:
步骤一:将所述卷积模块输入的特征图像输入至所述第一卷积层获得第一特征图像;第一卷积层的卷积核可以为N1*m*m*N2;N1为所述卷积模块输入的特征图像的通道数,N2为第一特征图像的通道数;N1>N2;Step 1: input the feature image input by the convolution module to the first convolution layer to obtain the first feature image; the convolution kernel of the first convolution layer may be N1 * m * m * N2; N1 is Describe the number of channels of the feature image input by the convolution module, N2 is the number of channels of the first feature image; N1> N2;
步骤二:将第一特征图像输入至所述第二卷积层获得第二特征图像;第一卷积层的卷积核可以为N2*m*m*N3;N3为第二特征图像的通道数;N3>N2;Step 2: Input the first feature image into the second convolution layer to obtain the second feature image; the convolution kernel of the first convolution layer may be N2 * m * m * N3; N3 is the channel of the second feature image Number; N3> N2;
步骤三:将所述卷积模块输入的特征图像和所述第二特征图像合并后,确定为所述卷积模块输出的特征图像。Step 3: After combining the feature image input by the convolution module and the second feature image, it is determined as the feature image output by the convolution module.
在一种具体的实施例中,第二卷积层输出的特征图像的个数可以与第一卷积层输入的特征图像的个数相等。即,N1=N2。In a specific embodiment, the number of feature images output by the second convolution layer may be equal to the number of feature images input by the first convolution layer. That is, N1 = N2.
上文所描述的乳腺影像对应的特征图像的确定方式仅为一种可能的实现方式,在其它可能的实现方式中,也可以通过其它方式确定乳腺影像对应的特征图像,具体不做限定。The method for determining the feature image corresponding to the breast image described above is only one possible implementation manner. In other possible implementation manners, the feature image corresponding to the breast image may also be determined by other methods, which is not specifically limited.
需要说明的是:本发明实施例中的激活函数可以为多种类型的激活函数,比如,可以为线性整流函数(Rectified Linear Unit,ReLU),具体不做限定;It should be noted that the activation function in the embodiment of the present invention may be multiple types of activation functions, for example, it may be a linear rectification function (Rectified Linear Unit, ReLU), which is not specifically limited;
由于本发明实施例中输入的图像为二维图像,因此,本发明实施例中的特征提取模块可以为(2Dimensions,2D)卷积神经网络中的特征提取模块,相应地,第一卷积层的卷积核大小可以为m*m、第二卷积层的卷积核大小可以为n*n;m和n可以相同也可以不同,在此不做限定;其中,m,n为大于或等于1的整数。第一卷积层输出的特征图像的个数小于所述第一卷积层输入的特征图像的个数;所述第二卷积层输出的特征图像的个数大于所述第一卷积层输出的特征图像的个数。Since the input image in the embodiment of the present invention is a two-dimensional image, the feature extraction module in the embodiment of the present invention may be a feature extraction module in a (2Dimensions, 2D) convolutional neural network. Accordingly, the first convolution layer The size of the convolution kernel can be m * m, and the size of the second convolution layer can be n * n; m and n can be the same or different, which is not limited here; where, m, n is greater than or An integer equal to 1. The number of feature images output by the first convolution layer is less than the number of feature images input by the first convolution layer; the number of feature images output by the second convolution layer is greater than the first convolution layer The number of output feature images.
进一步的,为优化特征提取模块,一种可能的实现方式,如图3c所示,所述第一卷积层和所述第二卷积层之间还包括第三卷积层;所述第三卷积层输入的特征图像为所述第一卷积层输出的图像,所述第三卷积层输出的特征图像为所述第二卷积层输入的图像。Further, to optimize the feature extraction module, a possible implementation manner, as shown in FIG. 3c, further includes a third convolution layer between the first convolution layer and the second convolution layer; The feature image input by the three convolution layers is the image output by the first convolution layer, and the feature image output by the third convolution layer is the image input by the second convolution layer.
其中,第三卷积层的卷积核大小可以为k*k,k与m,n可以相同,也可以不同,在此不做限定。The size of the convolution kernel of the third convolutional layer may be k * k, and k may be the same as m or n, or may be different, which is not limited herein.
一个具体的实施例中,所述第一卷积层的卷积核的大小为3*3;所述第二卷积层的卷积核的大小为3*3;所述第三卷积层的卷积核的大小为1*1。In a specific embodiment, the size of the convolution kernel of the first convolution layer is 3 * 3; the size of the convolution kernel of the second convolution layer is 3 * 3; the third convolution layer The size of the convolution kernel is 1 * 1.
通过上述卷积核的设置方式,可以有效的提高特征提取的感知野,有利于提高乳腺病灶识别的准确度。Through the above-mentioned setting method of the convolution kernel, the perception field of feature extraction can be effectively improved, which is beneficial to improve the accuracy of breast lesion recognition.
不同尺寸的特征图像可以为不同像素的特征图像,比如像素为500×500的特征图像与像素为1000×1000的特征图像为不同尺寸的特征图像。The feature images of different sizes may be feature images of different pixels, for example, the feature image with pixels 500 × 500 and the feature image with pixels 1000 × 1000 are feature images with different sizes.
可选地,采用预先训练好的乳腺病灶检测模型提取乳腺影像的不同尺寸的特征图像,模型是采用2D卷积神经网络对已标记的多个乳腺影像进行训练后确定的。Optionally, a pre-trained breast lesion detection model is used to extract feature images of different sizes of breast images. The model is determined after training a plurality of labeled breast images using a 2D convolutional neural network.
可选地,在提取乳腺影像的不同尺寸的特征图像之前,将图像缩放到特定尺寸,使各方向上像素与实际长度的比例尺一定。Optionally, before extracting feature images of different sizes of the mammary gland image, the image is scaled to a specific size so that the scale of the pixels in each direction is the same as the actual length.
另一种可能的实现方式,所述特征提取模块包括N/2个下采样卷积块和N/2个上采样卷积块;所述获取所述乳腺影像的不同尺寸的特征图像,包括:In another possible implementation manner, the feature extraction module includes N / 2 down-sampling convolution blocks and N / 2 up-sampling convolution blocks; and acquiring feature images of different sizes of the breast image includes:
将所述乳腺影像依次通过N/2个下采样卷积块提取N/2个所述乳腺影像的第一特征图像;Successively extracting the first feature image of the breast image through N / 2 down-sampling convolution blocks through the breast image;
将第N/2个下采样卷积块输出的第一特征图像依次通过N/2个上采样卷积块提取N/2个所述乳腺影像的第二特征图像,每个上采样卷积块提取的第二特征图像的尺寸均不同;The first feature image output from the N / 2th down-sampling convolution block is used to sequentially extract N / 2 second feature images of the mammary gland image through the N / 2 up-sampling convolution block, each up-sampling convolution block The sizes of the extracted second feature images are different;
将尺寸相同的第一特征图像和第二特征图像合并后,确定N个所述乳腺影像的不同尺寸的特征图像。After merging the first feature image and the second feature image of the same size, N feature images of different sizes of the breast images are determined.
为提高特征提取的感知野,提高特征提取的性能,一种可能的实现方式,所述特征提取模块之前还包括特征预处理模块;所述特征预处理模块包括一个卷积层,一个BN层,一个Relu层和一个池化层;所述特征预处理模块的卷积核大小大于所述N个卷积模块中任一卷积模块的卷积核的大小。In order to improve the perception field of feature extraction and improve the performance of feature extraction, a possible implementation manner, the feature extraction module also includes a feature preprocessing module before; the feature preprocessing module includes a convolution layer and a BN layer, One Relu layer and one pooling layer; the size of the convolution kernel of the feature preprocessing module is larger than that of any of the N convolution modules.
优选的,所述卷积层的卷积核大小可以为7*7,间隔为2个像素。池化层为2*2的最大值池化。通过特征预处理模块,可以将图像面积迅速缩小,边长变为原有1/4,有效的提高特征图像的感知野,快速的提取浅层特征,有效 的减少原始信息的损失。Preferably, the size of the convolution kernel of the convolution layer may be 7 * 7, and the interval is 2 pixels. The pooling layer is 2 * 2 maximum pooling. Through the feature preprocessing module, the image area can be rapidly reduced, and the side length becomes 1/4, effectively improving the perception field of the feature image, quickly extracting shallow features, and effectively reducing the loss of original information.
一种可能的实现方式,所述特征预处理模块包括连续的多个卷积层,一个BN层,一个Relu层和一个池化层;所述特征预处理模块的卷积核大小与所述N个卷积模块中的最大的卷积核的大小相等。A possible implementation manner, the feature preprocessing module includes a plurality of continuous convolutional layers, a BN layer, a Relu layer, and a pooling layer; the size of the convolution kernel of the feature preprocessing module and the N The largest convolution kernel in each convolution module has the same size.
特征图像经过特征预处理模块的步骤可以包括:将所述乳腺影像输入至特征预处理模块,获得预处理的特征图像;将所述预处理的特征图像作为所述特征提取模块的输入。The step of the feature image passing through the feature preprocessing module may include: inputting the breast image to the feature preprocessing module to obtain a preprocessed feature image; and using the preprocessed feature image as an input of the feature extraction module.
步骤203:针对所述乳腺影像的不同尺寸的特征图像中的任意一个特征图像,从所述特征图像中确定出乳腺病灶识别框。Step 203: For any one of the feature images of different sizes of the breast image, determine a breast lesion recognition frame from the feature image.
可选地,采用预先训练好的乳腺病灶检测模型从特征图像中确定出乳腺病灶识别框,乳腺病灶检测模型是采用2D卷积神经网络对已标记乳腺病灶的多个乳腺影像进行训练后确定的。从特征图像中确定出的乳腺病灶识别框框选的区域并不一定都包含乳腺病灶,故需要根据乳腺病灶识别框的乳腺病灶概率对各乳腺病灶识别框进行筛选,将乳腺病灶概率小于预设阈值的乳腺病灶识别框删除,其中,乳腺病灶概率为乳腺病灶识别框框选的区域为乳腺病灶的概率。Optionally, a pre-trained breast lesion detection model is used to determine the breast lesion recognition frame from the feature image. The breast lesion detection model is determined after training multiple breast images of the marked breast lesion using a 2D convolutional neural network. . The area framed by the breast lesion identification frame determined from the feature image does not necessarily contain breast lesions, so each breast lesion identification frame needs to be screened according to the breast lesion probability of the breast lesion identification frame, and the breast lesion probability is less than the preset threshold The breast lesion identification frame is deleted, where the breast lesion probability is the probability that the area framed by the breast lesion identification frame is the breast lesion.
步骤204:根据从各特征图像中确定出的乳腺病灶识别框,确定乳腺影像的乳腺病灶。Step 204: Determine the breast lesion of the breast image according to the breast lesion identification frame determined from each feature image.
具体的,确定出乳腺病灶识别框之后,将识别框作为乳腺影像中的乳腺病灶输出,输出的乳腺病灶参数包括乳腺病灶的中心坐标以及乳腺病灶的直径,其中乳腺病灶的中心坐标为乳腺病灶识别框的中心坐标,乳腺病灶的直径为乳腺病灶识别框的中心至其中一个面的距离。Specifically, after identifying the breast lesion recognition frame, the recognition frame is output as the breast lesion in the breast image, and the output breast lesion parameters include the central coordinates of the breast lesion and the diameter of the breast lesion, wherein the central coordinates of the breast lesion are the breast lesion identification The center coordinate of the frame, the diameter of the breast lesion is the distance from the center of the breast lesion identification frame to one of the faces.
由于提取乳腺影像的不同尺寸的特征图像,并识别每一个特征图像中的乳腺病灶,故既能检测到大尺寸的乳腺病灶,同时也能检测到小尺寸的乳腺病灶,提高了乳腺病灶检测的精度。其次,相较于人工判断乳腺影像中是否存在乳腺病灶的方法,本申请中自动检测乳腺病灶的方法有效地提高了乳腺病灶识别效率。Because feature images of different sizes of breast images are extracted and the breast lesions in each feature image are identified, both large-sized breast lesions and small-sized breast lesions can be detected, which improves the detection of breast lesions Precision. Secondly, compared with the method of manually judging whether there is a breast lesion in the breast image, the method of automatically detecting the breast lesion in the present application effectively improves the recognition efficiency of the breast lesion.
由于从各个特征图像中确定出的乳腺病灶识别框可能存在多个识别框对应一个乳腺病灶,若直接根据乳腺病灶识别框的数量确定乳腺影像中乳腺病灶的数量,将导致检测得到的乳腺病灶数量存在很大偏差,故需要将各特征图像转化为同一尺寸的特征图像并对齐,然后将从各特征图像中确定出的乳腺病灶识别框进行筛选,并将筛选后的乳腺病灶识别框确定为乳腺影像中的乳腺病灶。Since the breast lesion identification frame determined from each feature image may have multiple identification frames corresponding to one breast lesion, if the number of breast lesions in the breast image is directly determined according to the number of breast lesion identification frames, the number of detected breast lesions will result There is a large deviation, so it is necessary to convert each feature image into a feature image of the same size and align it, and then screen the breast lesion recognition frame determined from each feature image, and determine the screened breast lesion recognition frame as the breast Breast lesions in the image.
为进一步提高乳腺病灶的识别准确率,一种可能的实现方式,所述乳腺影像包括不同侧乳房的不同投照位的乳腺影像;所述将所述乳腺影像输入至特征提取模块,包括:In order to further improve the recognition accuracy of breast lesions, a possible implementation manner, the breast image includes breast images of different breasts with different projection positions; the input of the breast image to the feature extraction module includes:
将所述乳腺影像的同一投照位的另一侧乳房的乳腺影像作为所述乳腺影像的参考影像,输入至所述特征提取模块,获得参考特征图像;Taking the breast image of the other breast in the same projection position of the breast image as the reference image of the breast image, and inputting it to the feature extraction module to obtain a reference feature image;
所述针对所述乳腺影像的不同尺寸的特征图像中的任意一个特征图像,从所述特征图像中确定出乳腺病灶识别框;包括:Any one of the feature images of different sizes of the breast image, the breast lesion identification frame is determined from the feature image; including:
确定所述特征图像中的第一乳腺病灶识别框和所述参考特征图像中的第二乳腺病灶识别框;Determining a first breast lesion recognition frame in the feature image and a second breast lesion recognition frame in the reference feature image;
若确定所述第一乳腺病灶识别框和所述第二乳腺病灶识别框的位置和/或大小都相同,则删除所述第一乳腺病灶识别框。If it is determined that the positions and / or sizes of the first breast lesion identification frame and the second breast lesion identification frame are the same, the first breast lesion identification frame is deleted.
可选地,乳腺病灶识别框的筛选过程包括以下步骤,如图3所示:Optionally, the screening process of the breast lesion identification frame includes the following steps, as shown in Figure 3:
步骤301,从各特征图像的乳腺病灶识别框中确定乳腺病灶概率最大的乳腺病灶识别框。Step 301: Determine the breast lesion identification frame with the highest probability of breast lesion from the breast lesion identification frame of each feature image.
步骤302,计算乳腺病灶概率最大的乳腺病灶识别框与其他乳腺病灶识别框的交并比。Step 302: Calculate the intersection ratio of the breast lesion identification frame with the largest breast lesion probability and the other breast lesion identification frames.
步骤303,将交并比大于预设阈值的其他乳腺病灶识别框删除。Step 303: Delete the identification frame of other breast lesions whose cross-combination ratio is greater than a preset threshold.
步骤304,从剩余的其他乳腺病灶识别框中确定乳腺病灶概率最大的乳腺病灶识别框,重复执行乳腺病灶识别框的筛选过程,直到没有剩余的其他乳腺病灶识别框。Step 304: Determine the breast lesion identification frame with the highest probability of breast lesions from the remaining other breast lesion identification frames, and repeat the screening process of the breast lesion identification frame until there are no other remaining breast lesion identification frames.
下面结合具体的例子对上述乳腺病灶识别框的筛选过程进行说明,设定 各特征图像中确定出的乳腺病灶识别框分别为A、B、C、D、E、F,上述各乳腺病灶识别框的乳腺病灶概率分别为:P(A)=0.9、P(B)=0.85、P(C)=0.95、P(D)=0.75、P(E)=0.96、P(F)=0.65。将上述乳腺病灶识别框按照乳腺病灶概率从大到小进行排序后为:E、C、A、B、D、F,排序后可知各特征图像的乳腺病灶识别框中乳腺病灶概率最大的乳腺病灶识别框为E,然后分别计算乳腺病灶识别框E与其他各乳腺病灶识别框之间的交并比IOU,其中交并比的计算方式如式(1)所示:The screening process of the breast lesion identification frame described above will be described in conjunction with specific examples, and the breast lesion identification frames identified in each feature image are set to A, B, C, D, E, and F, and the above breast lesion identification frames The probabilities of breast lesions are: P (A) = 0.9, P (B) = 0.85, P (C) = 0.95, P (D) = 0.75, P (E) = 0.96, P (F) = 0.65. The above breast lesion recognition frame is sorted according to the probability of breast lesion from large to small: E, C, A, B, D, F, after sorting, we can know the breast lesion with the highest probability of breast lesion in the breast lesion recognition frame of each feature image The recognition frame is E, and then calculate the intersection ratio IOU between the breast lesion recognition frame E and each other breast lesion recognition frame. The calculation method of the intersection ratio is shown in equation (1):
Figure PCTCN2019082687-appb-000001
Figure PCTCN2019082687-appb-000001
其中,m为乳腺病灶概率最大的乳腺病灶识别框,n为与乳腺病灶识别框m比较的乳腺病灶识别框,IOU为乳腺病灶识别框m与乳腺病灶识别框n之间的交并比。Among them, m is the breast lesion identification frame with the highest probability of breast lesions, n is the breast lesion identification frame compared with the breast lesion identification frame m, and IOU is the intersection ratio between the breast lesion identification frame m and the breast lesion identification frame n.
设定预设阈值为0.5,若乳腺病灶识别框C与乳腺病灶识别框E之间的交并比大于0.5,乳腺病灶识别框A与乳腺病灶识别框E之间的交并比大于0.5,乳腺病灶识别框B、乳腺病灶识别框D、乳腺病灶识别框F与乳腺病灶识别框E之间的交并比均小于0.5,则删除乳腺病灶识别框C和乳腺病灶识别框A,将乳腺病灶识别框E确定为乳腺影像中的乳腺病灶。Set the preset threshold to 0.5, if the intersection ratio between the breast lesion identification frame C and the breast lesion identification frame E is greater than 0.5, the intersection ratio between the breast lesion identification frame A and the breast lesion identification frame E is greater than 0.5, the breast If the intersection ratio of the lesion identification frame B, the breast lesion identification frame D, the breast lesion identification frame F and the breast lesion identification frame E is less than 0.5, delete the breast lesion identification frame C and the breast lesion identification frame A, and identify the breast lesion Box E is identified as a breast lesion in the breast image.
进一步地,将剩余的其他乳腺病灶识别框B、D、F根据乳腺病灶概率进行排序,确定乳腺病灶概率最大的乳腺病灶识别框为乳腺病灶识别框B,然后计算乳腺病灶识别框B与乳腺病灶识别框D之间的交并比以及乳腺病灶识别框B与乳腺病灶识别框F之间的交并比。若乳腺病灶识别框B与乳腺病灶识别框D之间的交并比大于0.5,乳腺病灶识别框B与乳腺病灶识别框F之间的交并比小于0.5,则删除乳腺病灶识别框D,将乳腺病灶识别框B以及乳腺病灶识别框F确定为乳腺影像中的乳腺病灶。由于根据乳腺病灶识别框的乳腺病灶概率以及乳腺病灶识别框之间的交并比对各特征图像中确定出的乳腺病灶识别框进行筛选,避免重复检测乳腺影像中的同一个乳腺病灶并输出,提高检测乳腺影像中乳腺病灶数量的准确性。Further, the remaining other breast lesion recognition frames B, D, F are sorted according to the breast lesion probability, and the breast lesion recognition frame with the highest breast lesion probability is determined as the breast lesion recognition frame B, and then the breast lesion recognition frame B and the breast lesion are calculated The cross ratio between the recognition frames D and the cross ratio between the breast lesion recognition frame B and the breast lesion recognition frame F. If the intersection ratio between the breast lesion identification frame B and the breast lesion identification frame D is greater than 0.5, and the intersection ratio between the breast lesion identification frame B and the breast lesion identification frame F is less than 0.5, then delete the breast lesion identification frame D, replace The breast lesion recognition frame B and the breast lesion recognition frame F are determined as breast lesions in the breast image. Because the breast lesion identification frame identified in each feature image is screened based on the breast lesion probability of the breast lesion identification frame and the intersection and comparison between the breast lesion identification frames, avoiding repeated detection and output of the same breast lesion in the breast image, Improve the accuracy of detecting the number of breast lesions in breast images.
下面具体介绍一下通过2D卷积神经网络对已标记乳腺病灶的多个乳腺影像进行训练确定乳腺病灶检测模型过程,如图4所示,包括以下步骤:The following is a detailed introduction to the process of training the multiple breast images of the marked breast lesions through the 2D convolutional neural network to determine the detection model of the breast lesions, as shown in Figure 4, including the following steps:
步骤401,获取乳腺影像作为训练样本。Step 401: Obtain a breast image as a training sample.
具体地,可以将获取的多幅乳腺影像直接作为训练样本,也可以对获取的多幅乳腺影像进行增强操作,扩大训练样本的数据量,增强操作包括但不限于:随机上下左右平移设定像素(比如0~20像素)、随机旋转设定角度(比如-15~15度)、随机缩放设定倍数(比如0.85~1.15倍)。Specifically, the acquired multiple breast images can be used directly as training samples, or the acquired multiple breast images can be enhanced to expand the data volume of the training samples. Enhancement operations include, but are not limited to: randomly setting pixels up, down, left, and right (Such as 0-20 pixels), random rotation setting angle (such as -15-15 degrees), random zoom setting multiple (such as 0.85-1.15 times).
步骤402,人工标记训练样本中的乳腺病灶。Step 402: Manually mark the breast lesions in the training sample.
可以通过医生等专业人员对训练样本中的乳腺病灶进行标记,标记的内容包括乳腺病灶的中心坐标以及乳腺病灶的直径。具体地,可以由多名医生对乳腺病灶进行标注,并通过多人投票合成的方式确定最终的乳腺病灶以及乳腺病灶参数,结果用掩码图的方式保存。需要说明的是,人工标记训练样本中乳腺病灶与训练样本的增强操作不分先后,可以先人工标记训练样本中的乳腺病灶,然后再将标记乳腺病灶的训练样本进行增强操作,也可以先将训练样本进行增强操作,然后人工对增强操作后的训练样本进行标记。The breast lesions in the training sample can be marked by doctors and other professionals, and the content of the marking includes the central coordinates of the breast lesions and the diameter of the breast lesions. Specifically, multiple doctors can mark the breast lesions, and determine the final breast lesions and the parameters of the breast lesions through a multiple-vote synthesis method, and the results are saved in the form of a mask. It should be noted that the manual labeling of the breast lesions in the training sample and the training sample are in no particular order of enhancement. You can manually mark the breast lesions in the training sample, and then the training sample of the labeled breast lesions can be enhanced. The training samples are enhanced, and then the training samples after the enhancement operations are manually marked.
步骤403,将训练样本输入卷积神经网络进行训练,确定乳腺病灶识别模型。In step 403, the training samples are input to the convolutional neural network for training to determine the breast lesion recognition model.
该卷积神经网络的结构包括输入层、下采样卷积块、上采样卷积块、目标检测网络以及输出层。将训练样本进行预处理后输入上述卷积神经网络,将输出的乳腺病灶与预先标记的训练样本的掩码图进行损失函数计算,然后采用反向传播算法以及sgd优化算法反复迭代,确定乳腺病灶检测模型。The structure of the convolutional neural network includes an input layer, a down-sampling convolution block, an up-sampling convolution block, a target detection network, and an output layer. Pre-process the training samples and input them into the convolutional neural network, calculate the loss function of the output breast lesions and the mask image of the pre-labeled training samples, and then iterate iteratively using the back propagation algorithm and the sgd optimization algorithm to determine the breast lesions Detection model.
进一步地,采用上述训练确定的乳腺病灶检测模型提取乳腺影像的不同尺寸的特征图像的过程,包括以下步骤:Further, the process of extracting feature images of different sizes of breast images using the breast lesion detection model determined by the above training includes the following steps:
步骤一,将乳腺影像依次通过N/2个下采样卷积块提取N个乳腺影像的第一特征图像。In the first step, the mammary gland image is successively passed through N / 2 down-sampling convolution blocks to extract the first feature images of the N mammary gland images.
每个下采样卷积块提取的第一特征图像的尺寸均不同,N/2大于0。The size of the first feature image extracted by each down-sampling convolution block is different, and N / 2 is greater than 0.
可选地,下采样卷积块包括第一卷积层和第二卷积层、组连接层、前后 连接层、下采样层。Optionally, the down-sampling convolution block includes a first convolution layer and a second convolution layer, a group connection layer, a front-back connection layer, and a down-sampling layer.
步骤二,将第N/2个下采样卷积块输出的第一特征图像依次通过N/2个上采样卷积块提取N/2个乳腺影像的第二特征图像。Step 2: The first feature image output from the N / 2 down-sampling convolution block is sequentially used to extract the second feature image of N / 2 mammary gland images through the N / 2 up-sampling convolution block.
每个上采样卷积块提取的第二特征图像的尺寸均不同。The size of the second feature image extracted by each up-sampling convolution block is different.
可选地,上采样卷积块包括卷积层、组连接层、前后连接层、上采样层以及合成连接层。卷积层包括卷积运算,batch normalization层和RELU层。Optionally, the up-sampling convolution block includes a convolution layer, a group connection layer, a front-back connection layer, an up-sampling layer, and a synthesis connection layer. Convolution layer includes convolution operation, batch normalization layer and RELU layer.
步骤三,将尺寸相同的第一特征图像和第二特征图像合并后,确定N/2个乳腺影像的不同尺寸的特征图像。In step three, after combining the first feature image and the second feature image with the same size, feature images of different sizes of N / 2 breast images are determined.
通过上采样卷积块中的合成连接层将尺寸相同的第一特征图像和第二特征图像合并确定不同尺寸的特征图像。可选地,在合并时,是将第一特征图像和第二特征图像的通道数进行合并,合并后得到的特征图像的尺寸与第一特征图像和第二特征图像的尺寸相同。The first feature image and the second feature image of the same size are combined through the up-sampling convolution block to determine feature images of different sizes. Optionally, when merging, the number of channels of the first feature image and the second feature image are combined, and the size of the feature image obtained after the merging is the same as the size of the first feature image and the second feature image.
进一步地,采用上述训练确定的乳腺病灶检测模型从特征图像中确定出乳腺病灶识别框的过程,包括以下步骤:Further, the process of determining the breast lesion recognition frame from the feature image using the breast lesion detection model determined by the above training includes the following steps:
步骤一,针对特征图像中任意一个像素,以像素为中心,向四周扩散确定第一区域。Step 1: For any pixel in the feature image, with the pixel as the center, diffuse to the surroundings to determine the first area.
步骤二,在第一区域中根据预设规则设置多个预设框。Step 2: Set multiple preset frames in the first area according to preset rules.
由于乳腺病灶的形状不一,故可以将预设框设置为多种形状。预设规则可以是将预设框中心与第一区域的中心重合,也可以是预设框的角与第一区域的角重合等等。Due to the different shapes of breast lesions, the preset frame can be set to various shapes. The preset rule may be that the center of the preset frame coincides with the center of the first area, or that the corner of the preset frame coincides with the angle of the first area, and so on.
在一个具体的实施例中,乳腺病灶预设框选取的方式为,对于每个特征图的每个像素,认为其为一个锚点。在每个锚点上设置多个长宽比不一的预设框。In a specific embodiment, the way to select the preset frame of the breast lesion is that, for each pixel of each feature map, it is considered as an anchor point. Set multiple preset frames with different aspect ratios on each anchor point.
对于每个预设框,通过对特征图进行卷积,预测一个坐标和尺寸的偏移,以及置信度,根据坐标和尺寸的偏移,以及置信度,确定预设框。For each preset frame, the convolution of the feature map predicts a coordinate and size offset and confidence, and the preset frame is determined based on the coordinate and size offset and confidence.
步骤三,针对任意一个预设框,预测预设框与第一区域的位置偏差。Step 3: For any preset frame, predict the position deviation of the preset frame from the first area.
步骤四,根据位置偏差调整预设框后确定乳腺病灶识别框,并预测乳腺 病灶识别框的乳腺病灶概率。Step 4: Adjust the preset frame according to the position deviation to determine the breast lesion identification frame, and predict the breast lesion probability of the breast lesion identification frame.
其中,乳腺病灶概率为乳腺病灶识别框框选的区域为乳腺病灶的概率。通过预测预设框与第一区域的位置偏差,然后采用位置偏差调整预设框确定识别框,以使识别框更多地框选特征图中的乳腺病灶区域,提高乳腺病灶检测的准确性。The breast lesion probability is the probability that the area selected by the breast lesion identification frame is the breast lesion. By predicting the position deviation between the preset frame and the first area, and then using the position deviation to adjust the preset frame to determine the recognition frame, so that the recognition frame can more frame the breast lesion area in the feature map, and improve the accuracy of breast lesion detection.
具体的训练过程可以包括:将训练数据影像输入上述的卷积神经网络进行计算。传入时,将病灶不同窗宽窗位的多张影像传入。训练时,在卷积神经网络输出的预测框中,选取置信度最高的预测框集和与训练样本重合最大的预测框集合。将预测框置信度和样本标注的交叉熵,与训练样本的标注病灶和预测框的偏移的交叉熵,两者的加权和作为loss函数。通过反向传播的方法训练,训练的优化算法使用带有动量和阶梯衰减的sgd算法。The specific training process may include: inputting the training data image to the above-mentioned convolutional neural network for calculation. During the introduction, multiple images of different window widths and positions of the lesion are introduced. During training, in the prediction frame output by the convolutional neural network, the prediction frame set with the highest confidence and the prediction frame set with the maximum coincidence with the training sample are selected. The cross-entropy of the confidence of the prediction frame and the labeling of the sample and the cross-entropy of the labeled lesion of the training sample and the offset of the prediction frame are used as the loss function. Trained by the method of back propagation, the training optimization algorithm uses the sgd algorithm with momentum and step attenuation.
在算法使用过程中,通过预处理模块,将输入图像预处理,以提高特征提取的效果。In the process of using the algorithm, through the preprocessing module, the input image is preprocessed to improve the effect of feature extraction.
一种可能的实现方式,所述获取乳腺影像,包括:A possible implementation manner, the acquiring breast image includes:
步骤一、将拍摄的乳腺影像图像,根据高斯滤波,确定所述乳腺影像图像的二值化图像;Step 1: Determine the binary image of the breast image according to Gaussian filtering;
步骤二、获取所述二值化图像的连通区域,将连通区域中最大的区域对应于所述乳腺影像图像的区域作为分割出的乳腺图像;Step 2: Obtain the connected area of the binarized image, and use the area with the largest area in the connected area corresponding to the breast image as the segmented breast image;
步骤三、将所述分割出的乳腺图像添加至预设的图像模板中,生成预处理后的乳腺图像;并将所述预处理后的乳腺图像作为输入至所述特征提取模块的乳腺影像。Step 3: Add the segmented breast image to a preset image template to generate a pre-processed breast image; and use the pre-processed breast image as the breast image input to the feature extraction module.
具体的,预处理模块的输入为以Dicom格式形式保存的乳腺影像。预处理可以包括腺体分割和图像归一化;腺体分割的主要目的是将输入的乳腺影像中的乳腺部分提取出,剔除其他无关的干扰的图像;图像归一化是将图像化归为统一格式图像,具体的,包括:Specifically, the input of the preprocessing module is a breast image saved in Dicom format. Preprocessing can include gland segmentation and image normalization; the main purpose of gland segmentation is to extract the breast part of the input breast image to remove other unrelated interference images; image normalization is to normalize the image into Unified format images, specifically, include:
在步骤一中,具体的二值化的阈值可以通过求图像灰度直方图的最大类间距方法获得。In step one, the specific binarized threshold can be obtained by finding the maximum class interval of the grayscale histogram of the image.
在步骤二中,可以将二值化的结果,通过漫水法(flood fill)获得独立的区域块,并统计每个区域块的面积;将面积最大的区域块对应的图像上的区域,作为分割出来的乳腺图像。In step two, the binarized result can be obtained by flooding to obtain independent regional blocks, and the area of each regional block is counted; the area on the image corresponding to the largest regional block is used as Segmented breast image.
在步骤三中,预设的图像模板可以为黑色底板的正方形图像;具体的,可以将获得的分割出来的乳腺图像,通过加黑边填充的方式扩充为1:1的正方形图像。In step three, the preset image template may be a square image of a black bottom plate; specifically, the obtained divided breast image may be expanded into a 1: 1 square image by adding a black border.
另外,输出的乳腺影像可以通过像素缩放,例如,可以将图像差值缩放到4096像素×4096像素大小。In addition, the output breast image can be scaled by pixels, for example, the image difference can be scaled to 4096 pixels × 4096 pixels.
针对乳腺,由于乳腺照射剂量以及拍摄的外界因素等原因,可以通过调整乳腺的窗宽窗位,以获得更好的乳腺病灶识别的识别效果。一种可能的实现方式,所述将所述乳腺影像输入至特征提取模块之前,还包括:For the mammary gland, due to the dose of the mammary gland and external factors of the shooting, the window width and position of the mammary gland can be adjusted to obtain a better identification effect of breast lesion identification. In a possible implementation manner, before inputting the breast image to the feature extraction module, the method further includes:
获取所述乳腺影像的原始文件;Obtain the original file of the breast image;
在所述乳腺影像的原始文件中选取至少一组窗宽窗位,并获取所述至少一组窗宽窗位对应的图片格式的乳腺影像;Selecting at least one set of window width and window levels from the original file of the breast image, and acquiring a breast image in a picture format corresponding to the at least one set of window width and window levels;
根据所述至少一组窗宽窗位对应的图片格式的乳腺影像,作为输入至所述特征提取模块的乳腺影像。The mammary gland image according to the picture format corresponding to the at least one set of window width and window level is used as the mammary gland image input to the feature extraction module.
在一个具体实施例中,可以通过三组窗宽窗位,将dicom图像转换为png图像,例如,第一组窗宽为4000,窗位2000;第二组窗宽为1000;窗位为2000;第三组窗宽为1500,窗位为1500。In a specific embodiment, the dicom image can be converted into a png image through three sets of window width and window levels. For example, the first set of window width is 4000 and the window level is 2000; the second set of window width is 1000; the window level is 2000 ; The third group has a window width of 1500 and a window level of 1500.
基于相同的技术构思,本发明实施例提供了一种乳腺病灶识别的装置,如图5所示,该装置可以执行乳腺病灶识别的方法的流程,该装置包括获取单元501、处理单元502。Based on the same technical concept, an embodiment of the present invention provides a device for identifying breast lesions. As shown in FIG. 5, the device can perform the flow of a method for identifying breast lesions. The device includes an acquiring unit 501 and a processing unit 502.
获取单元501,用于获取乳腺影像;The obtaining unit 501 is used to obtain breast images;
处理单元502,用于将所述乳腺影像输入至特征提取模块中,获取所述乳腺影像不同尺寸的特征图像;所述特征提取模块包括N个卷积模块;所述N个卷积模块为下采样卷积块和/或上采样卷积块;每个下采样卷积块或上采样卷积块提取的特征图像的尺寸均不同,所述N个卷积模块的每个卷积模块中 包括第一卷积层、第二卷积层;所述第一卷积层输出的特征图像的个数小于所述第一卷积层输入的特征图像的个数;所述第二卷积层输出的特征图像的个数大于所述第一卷积层输出的特征图像的个数;N大于0;针对所述乳腺影像的不同尺寸的特征图像中的任意一个特征图像,从所述特征图像中确定出乳腺病灶识别框;根据从各特征图像中确定出的乳腺病灶识别框,确定乳腺影像的乳腺病灶。The processing unit 502 is configured to input the breast image into a feature extraction module to obtain feature images of different sizes of the breast image; the feature extraction module includes N convolution modules; the N convolution modules are as follows Sampling convolutional blocks and / or upsampling convolutional blocks; the size of the feature image extracted by each downsampling convolutional block or upsampling convolutional block is different, and each of the N convolutional modules includes The first convolution layer and the second convolution layer; the number of feature images output by the first convolution layer is less than the number of feature images input by the first convolution layer; the second convolution layer outputs The number of feature images of is greater than the number of feature images output by the first convolutional layer; N is greater than 0; for any one of the feature images of different sizes of the breast image, from the feature image Identify the breast lesion recognition frame; determine the breast lesion of the breast image according to the breast lesion recognition frame determined from each feature image.
一种可能的实现方式,所述处理单元502,具体用于:In a possible implementation manner, the processing unit 502 is specifically configured to:
将所述乳腺影像依次通过N/2个下采样卷积块提取N/2个所述乳腺影像的第一特征图像;将第N/2个下采样卷积块输出的第一特征图像依次通过N/2个上采样卷积块提取N/2个所述乳腺影像的第二特征图像,每个上采样卷积块提取的第二特征图像的尺寸均不同;将尺寸相同的第一特征图像和第二特征图像合并后,确定N个所述乳腺影像的不同尺寸的特征图像。The breast image is sequentially passed through N / 2 down-sampling convolution blocks to extract N / 2 first feature images of the breast image; the first feature images output from the N / 2 down-sampling convolution block are sequentially passed N / 2 up-sampling convolution blocks extract N / 2 second feature images of the mammography image, the size of the second feature images extracted by each up-sampling convolution block are different; the first feature images of the same size After merging with the second feature image, N feature images of different sizes of the breast images are determined.
一种可能的实现方式,所述特征处理模块之前还包括特征预处理模块;所述处理单元502,具体用于:A possible implementation manner, before the feature processing module further includes a feature preprocessing module; the processing unit 502 is specifically configured to:
将所述乳腺影像输入至所述特征预处理模块中,所述特征预处理模块包括一个卷积层,一个BN层,一个Relu层和一个池化层;所述特征预处理模块的卷积核大小大于所述N个卷积模块中的卷积核的大小;Input the breast image into the feature preprocessing module, the feature preprocessing module includes a convolution layer, a BN layer, a Relu layer and a pooling layer; the convolution kernel of the feature preprocessing module The size is larger than the size of the convolution kernel in the N convolution modules;
或者,所述特征预处理模块包括连续的多个卷积层,一个BN层,一个Relu层和一个池化层;所述特征预处理模块的卷积核大小与所述N个卷积模块中的最大的卷积核的大小相等。Alternatively, the feature pre-processing module includes multiple consecutive convolutional layers, a BN layer, a Relu layer, and a pooling layer; the size of the convolution kernel of the feature pre-processing module and the N convolution modules The size of the largest convolution kernel is equal.
一种可能的实现方式,所述获取单元501,用于:In a possible implementation manner, the obtaining unit 501 is configured to:
获取所述乳腺影像的原始文件;Obtain the original file of the breast image;
所述处理单元502,具体用于:The processing unit 502 is specifically used for:
在所述乳腺影像的原始文件中选取至少一组窗宽窗位,并获取所述至少一组窗宽窗位对应的图片格式的乳腺影像;根据所述至少一组窗宽窗位对应的图片格式的乳腺影像,作为输入至所述特征提取模块的乳腺影像。Selecting at least one group of window width and window levels from the original file of the breast image, and obtaining a breast image in a picture format corresponding to the at least one group of window width and window levels; according to the pictures corresponding to the at least one group of window width and window levels The breast image in the format is used as the breast image input to the feature extraction module.
一种可能的实现方式,所述乳腺影像包括不同侧乳房的不同投照位的乳 腺影像;所述获取单元501,用于:In a possible implementation manner, the mammary gland image includes mammary gland images with different projection positions of different breasts; the acquiring unit 501 is configured to:
将所述乳腺影像的同一投照位的另一侧乳房的乳腺影像作为所述乳腺影像的参考影像,输入至所述特征提取模块,获得参考特征图像;确定所述特征图像中的第一乳腺病灶识别框和所述参考特征图像中的第二乳腺病灶识别框;若确定所述第一乳腺病灶识别框和所述第二乳腺病灶识别框的位置和/或大小都相同,则删除所述第一乳腺病灶识别框。Take the breast image of the other breast of the same projection position of the breast image as the reference image of the breast image, and input it to the feature extraction module to obtain a reference feature image; determine the first breast in the feature image The lesion identification frame and the second breast lesion identification frame in the reference feature image; if it is determined that the positions and / or sizes of the first breast lesion identification frame and the second breast lesion identification frame are the same, delete the The first breast lesion identification frame.
基于相同的技术构思,本发明实施例提供了一种计算机设备,包括至少一个处理单元以及至少一个存储单元,其中,所述存储单元存储有计算机程序,当所述程序被所述处理单元执行时,使得所述处理单元执行乳腺病灶识别的方法的步骤。如图6所示,为本发明实施例中所述的计算机设备的硬件结构示意图,该计算机设备具体可以为台式计算机、便携式计算机、智能手机、平板电脑等。具体地,该计算机设备可以包括存储器801、处理器802及存储在存储器上的计算机程序,所述处理器802执行所述程序时实现上述实施例中的任一乳腺病灶识别的方法的步骤。其中,存储器801可以包括只读存储器(ROM)和随机存取存储器(RAM),并向处理器802提供存储器801中存储的程序指令和数据。Based on the same technical concept, an embodiment of the present invention provides a computer device including at least one processing unit and at least one storage unit, wherein the storage unit stores a computer program when the program is executed by the processing unit So that the processing unit executes the steps of the method for identifying breast lesions. As shown in FIG. 6, it is a schematic diagram of the hardware structure of the computer device described in the embodiment of the present invention. The computer device may specifically be a desktop computer, a portable computer, a smart phone, a tablet computer, or the like. Specifically, the computer device may include a memory 801, a processor 802, and a computer program stored on the memory. When the processor 802 executes the program, any steps of the method for identifying a breast lesion in the foregoing embodiments are implemented. The memory 801 may include a read-only memory (ROM) and a random access memory (RAM), and provide the processor 802 with program instructions and data stored in the memory 801.
进一步地,本申请实施例中所述的计算机设备还可以包括输入装置803以及输出装置804等。输入装置803可以包括键盘、鼠标、触摸屏等;输出装置804可以包括显示设备,如液晶显示器(Liquid Crystal Display,LCD)、阴极射线管(Cathode Ray Tube,CRT),触摸屏等。存储器801,处理器802、输入装置803和输出装置804可以通过总线或者其他方式连接,图6中以通过总线连接为例。处理器802调用存储器801存储的程序指令并按照获得的程序指令执行上述实施例提供的乳腺病灶识别的方法。Further, the computer equipment described in the embodiments of the present application may further include an input device 803 and an output device 804. The input device 803 may include a keyboard, a mouse, a touch screen, etc .; the output device 804 may include a display device, such as a liquid crystal display (Liquid Crystal Display, LCD), a cathode ray tube (Cathode Ray Tube, CRT), a touch screen, and the like. The memory 801, the processor 802, the input device 803, and the output device 804 may be connected through a bus or in other ways. In FIG. 6, connection through a bus is used as an example. The processor 802 calls the program instructions stored in the memory 801 and executes the method for identifying breast lesions provided in the foregoing embodiments according to the obtained program instructions.
基于相同的技术构思,本发明实施例还提供了一种计算机可读存储介质,其存储有可由计算机设备执行的计算机程序,当所述程序在计算机设备上运行时,使得所述计算机设备执行乳腺病灶识别的方法的步骤。Based on the same technical concept, an embodiment of the present invention also provides a computer-readable storage medium that stores a computer program executable by a computer device, and when the program runs on the computer device, causes the computer device to execute a breast Steps of the method of lesion identification.
基于相同的技术构思,本发明实施例还提供了一种计算机程序产品,所 述计算机程序产品包括存储在计算机可读存储介质上的计算机程序,所述计算机程序包括程序指令,当所述程序指令被计算机设备执行时,使所述计算机设备执行上述乳腺影像病灶识别的方法的步骤。Based on the same technical concept, an embodiment of the present invention also provides a computer program product, the computer program product includes a computer program stored on a computer-readable storage medium, the computer program includes program instructions, when the program instructions When executed by a computer device, the computer device is caused to perform the steps of the method for identifying breast imaging lesions.
本领域内的技术人员应明白,本发明的实施例可提供为方法、或计算机程序产品。因此,本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且,本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。Those skilled in the art should understand that the embodiments of the present invention may be provided as methods or computer program products. Therefore, the present invention may take the form of an entirely hardware embodiment, an entirely software embodiment, or an embodiment combining software and hardware. Moreover, the present invention may take the form of a computer program product implemented on one or more computer usable storage media (including but not limited to disk storage, CD-ROM, optical storage, etc.) containing computer usable program code.
本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。The present invention is described with reference to flowcharts and / or block diagrams of methods, devices (systems), and computer program products according to embodiments of the present invention. It should be understood that each flow and / or block in the flowchart and / or block diagram and a combination of the flow and / or block in the flowchart and / or block diagram may be implemented by computer program instructions. These computer program instructions can be provided to the processor of a general-purpose computer, special-purpose computer, embedded processing machine, or other programmable data processing device to produce a machine that enables the generation of instructions executed by the processor of the computer or other programmable data processing device A device for realizing the functions specified in one block or multiple blocks of one flow or multiple blocks of a flowchart.
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。These computer program instructions may also be stored in a computer readable memory that can guide a computer or other programmable data processing device to work in a specific manner, so that the instructions stored in the computer readable memory produce an article of manufacture including an instruction device, the instructions The device implements the functions specified in one block or multiple blocks of the flowchart one flow or multiple flows and / or block diagrams.
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上,使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。These computer program instructions can also be loaded onto a computer or other programmable data processing device, so that a series of operating steps are performed on the computer or other programmable device to produce computer-implemented processing, which is executed on the computer or other programmable device The instructions provide steps for implementing the functions specified in one block or multiple blocks of the flowchart one flow or multiple flows and / or block diagrams.
尽管已描述了本发明的优选实施例,但本领域内的技术人员一旦得知了基本创造性概念,则可对这些实施例作出另外的变更和修改。所以,所附权利要求意欲解释为包括优选实施例以及落入本发明范围的所有变更和修改。 显然,本领域的技术人员可以对本发明进行各种改动和变型而不脱离本发明的精神和范围。这样,倘若本发明的这些修改和变型属于本发明权利要求及其等同技术的范围之内,则本发明也意图包含这些改动和变型在内。Although the preferred embodiments of the present invention have been described, those skilled in the art can make additional changes and modifications to these embodiments once they learn the basic inventive concept. Therefore, the appended claims are intended to be interpreted as including the preferred embodiments and all changes and modifications falling within the scope of the present invention. Obviously, those skilled in the art can make various modifications and variations to the present invention without departing from the spirit and scope of the present invention. In this way, if these modifications and variations of the present invention fall within the scope of the claims of the present invention and their equivalent technologies, the present invention is also intended to include these modifications and variations.

Claims (11)

  1. 一种乳腺影像病灶识别的方法,其特征在于,包括:A method for identifying breast imaging lesions, characterized in that it includes:
    获取乳腺影像;Get breast images;
    将所述乳腺影像输入至特征提取模块中,获取所述乳腺影像不同尺寸的特征图像;所述特征提取模块包括N个卷积模块;所述N个卷积模块为下采样卷积块和/或上采样卷积块;每个下采样卷积块或上采样卷积块提取的特征图像的尺寸均不同,所述N个卷积模块的每个卷积模块中包括第一卷积层、第二卷积层;所述第一卷积层输出的特征图像的个数小于所述第一卷积层输入的特征图像的个数;所述第二卷积层输出的特征图像的个数大于所述第一卷积层输出的特征图像的个数;N大于0;The breast image is input into a feature extraction module to obtain feature images of different sizes of the breast image; the feature extraction module includes N convolution modules; the N convolution modules are down-sampling convolution blocks and / or Or up-sampling convolution block; the size of the feature image extracted by each down-sampling convolution block or up-sampling convolution block is different, each of the N convolution modules includes a first convolution layer, The second convolutional layer; the number of feature images output by the first convolution layer is less than the number of feature images input by the first convolution layer; the number of feature images output by the second convolution layer Is greater than the number of feature images output by the first convolutional layer; N is greater than 0;
    针对所述乳腺影像的不同尺寸的特征图像中的任意一个特征图像,从所述特征图像中确定出乳腺病灶识别框;For any one of the feature images of different sizes of the breast image, determine the breast lesion identification frame from the feature image;
    根据从各特征图像中确定出的乳腺病灶识别框,确定乳腺影像的乳腺病灶。According to the breast lesion identification frame identified from each feature image, the breast lesion of the breast image is determined.
  2. 如权利要求1所述的方法,其特征在于,所述获取所述乳腺影像的不同尺寸的特征图像,包括:The method according to claim 1, wherein the acquiring feature images of different sizes of the breast image include:
    将所述乳腺影像依次通过N/2个下采样卷积块提取N/2个所述乳腺影像的第一特征图像;Successively extracting the first feature image of the breast image through N / 2 down-sampling convolution blocks through the breast image;
    将第N/2个下采样卷积块输出的第一特征图像依次通过N/2个上采样卷积块提取N/2个所述乳腺影像的第二特征图像,每个上采样卷积块提取的第二特征图像的尺寸均不同;The first feature image output from the N / 2th down-sampling convolution block is used to sequentially extract N / 2 second feature images of the mammary gland image through the N / 2 up-sampling convolution block, each up-sampling convolution block The sizes of the extracted second feature images are different;
    将尺寸相同的第一特征图像和第二特征图像合并后,确定N个所述乳腺影像的不同尺寸的特征图像。After merging the first feature image and the second feature image of the same size, N feature images of different sizes of the breast images are determined.
  3. 如权利要求1所述的方法,其特征在于,所述特征处理模块之前还包括特征预处理模块;所述将所述乳腺影像输入至特征提取模块之前,还包括:The method according to claim 1, wherein before the feature processing module, a feature preprocessing module is further included; before the breast image is input to the feature extraction module, the method further includes:
    将所述乳腺影像输入至所述特征预处理模块中,所述特征预处理模块包 括一个卷积层,一个BN层,一个Relu层和一个池化层;所述特征预处理模块的卷积核大小大于所述N个卷积模块中的卷积核的大小;Input the breast image into the feature preprocessing module, the feature preprocessing module includes a convolution layer, a BN layer, a Relu layer and a pooling layer; the convolution kernel of the feature preprocessing module The size is larger than the size of the convolution kernel in the N convolution modules;
    或者,所述特征预处理模块包括连续的多个卷积层,一个BN层,一个Relu层和一个池化层;所述特征预处理模块的卷积核大小与所述N个卷积模块中的最大的卷积核的大小相等。Alternatively, the feature pre-processing module includes multiple consecutive convolutional layers, a BN layer, a Relu layer, and a pooling layer; the size of the convolution kernel of the feature pre-processing module and the N convolution modules The size of the largest convolution kernel is equal.
  4. 如权利要求1所述的方法,其特征在于,所述将所述乳腺影像输入至特征提取模块之前,还包括:The method according to claim 1, wherein before the inputting the breast image to the feature extraction module, the method further comprises:
    获取所述乳腺影像的原始文件;Obtain the original file of the breast image;
    在所述乳腺影像的原始文件中选取至少一组窗宽窗位,并获取所述至少一组窗宽窗位对应的图片格式的乳腺影像;Selecting at least one set of window width and window levels from the original file of the breast image, and acquiring a breast image in a picture format corresponding to the at least one set of window width and window levels;
    根据所述至少一组窗宽窗位对应的图片格式的乳腺影像,作为输入至所述特征提取模块的乳腺影像。The mammary gland image according to the picture format corresponding to the at least one set of window width and window level is used as the mammary gland image input to the feature extraction module.
  5. 如权利要求1-4任一项所述的方法,其特征在于,所述乳腺影像包括不同侧乳房的不同投照位的乳腺影像;所述将所述乳腺影像输入至特征提取模块,包括:The method according to any one of claims 1 to 4, wherein the mammary gland image includes mammary gland images with different projection positions of different breasts; and the inputting the mammary gland image to a feature extraction module includes:
    将所述乳腺影像的同一投照位的另一侧乳房的乳腺影像作为所述乳腺影像的参考影像,输入至所述特征提取模块,获得参考特征图像;Taking the breast image of the other breast in the same projection position of the breast image as the reference image of the breast image, and inputting it to the feature extraction module to obtain a reference feature image;
    所述针对所述乳腺影像的不同尺寸的特征图像中的任意一个特征图像,从所述特征图像中确定出乳腺病灶识别框;包括:Any one of the feature images of different sizes of the breast image, the breast lesion identification frame is determined from the feature image; including:
    确定所述特征图像中的第一乳腺病灶识别框和所述参考特征图像中的第二乳腺病灶识别框;Determining a first breast lesion recognition frame in the feature image and a second breast lesion recognition frame in the reference feature image;
    若确定所述第一乳腺病灶识别框和所述第二乳腺病灶识别框的位置和/或大小都相同,则删除所述第一乳腺病灶识别框。If it is determined that the positions and / or sizes of the first breast lesion identification frame and the second breast lesion identification frame are the same, the first breast lesion identification frame is deleted.
  6. 一种乳腺影像病灶识别的装置,其特征在于,包括:A device for identifying breast imaging lesions, characterized in that it includes:
    获取单元,用于获取乳腺影像;Acquisition unit for acquiring mammary gland image;
    处理单元,用于将所述乳腺影像输入至特征提取模块中,获取所述乳腺影像不同尺寸的特征图像;所述特征提取模块包括N个卷积模块;所述N个 卷积模块为下采样卷积块或上采样卷积块;每个下采样卷积块或上采样卷积块提取的特征图像的尺寸均不同,所述N个卷积模块的每个卷积模块中包括第一卷积层、第二卷积层;所述第一卷积层输出的特征图像的个数小于所述第一卷积层输入的特征图像的个数;所述第二卷积层输出的特征图像的个数大于所述第一卷积层输出的特征图像的个数;N大于0;针对所述乳腺影像的不同尺寸的特征图像中的任意一个特征图像,从所述特征图像中确定出乳腺病灶识别框;根据从各特征图像中确定出的乳腺病灶识别框,确定乳腺影像的乳腺病灶。A processing unit for inputting the breast image into a feature extraction module to obtain feature images of different sizes of the breast image; the feature extraction module includes N convolution modules; the N convolution modules are for downsampling Convolution block or up-sampling convolution block; the size of the feature image extracted by each down-sampling convolution block or up-sampling convolution block is different, and each of the N convolution modules includes a first volume Multilayer, second convolutional layer; the number of feature images output by the first convolutional layer is less than the number of feature images input by the first convolutional layer; the feature images output by the second convolutional layer The number of is greater than the number of feature images output by the first convolutional layer; N is greater than 0; for any one of the feature images of different sizes of the breast image, the breast is determined from the feature image Lesion recognition frame; according to the breast lesion recognition frame determined from each characteristic image, determine the breast lesion of the breast image.
  7. 如权利要求6所述的装置,其特征在于,所述处理单元,具体用于:The apparatus according to claim 6, wherein the processing unit is specifically configured to:
    将所述乳腺影像依次通过N/2个下采样卷积块提取N/2个所述乳腺影像的第一特征图像;将第N/2个下采样卷积块输出的第一特征图像依次通过N/2个上采样卷积块提取N/2个所述乳腺影像的第二特征图像,每个上采样卷积块提取的第二特征图像的尺寸均不同;将尺寸相同的第一特征图像和第二特征图像合并后,确定N个所述乳腺影像的不同尺寸的特征图像。The breast image is sequentially passed through N / 2 down-sampling convolution blocks to extract N / 2 first feature images of the breast image; the first feature images output from the N / 2 down-sampling convolution block are sequentially passed N / 2 up-sampling convolution blocks extract N / 2 second feature images of the mammography image, the size of the second feature images extracted by each up-sampling convolution block are different; the first feature images of the same size After merging with the second feature image, N feature images of different sizes of the breast images are determined.
  8. 如权利要求6或7所述的装置,其特征在于,所述乳腺影像包括不同侧乳房的不同投照位的乳腺影像;所述处理单元,具体用于:The device according to claim 6 or 7, wherein the mammary gland image comprises mammary gland images with different projection positions of different breasts; the processing unit is specifically used for:
    将所述乳腺影像的同一投照位的另一侧乳房的乳腺影像作为所述乳腺影像的参考影像,输入至所述特征提取模块,获得参考特征图像;确定所述特征图像中的第一乳腺病灶识别框和所述参考特征图像中的第二乳腺病灶识别框;若确定所述第一乳腺病灶识别框和所述第二乳腺病灶识别框的位置和/或大小都相同,则删除所述第一乳腺病灶识别框。Take the breast image of the other breast of the same projection position of the breast image as the reference image of the breast image, and input it to the feature extraction module to obtain a reference feature image; determine the first breast in the feature image The lesion identification frame and the second breast lesion identification frame in the reference feature image; if it is determined that the positions and / or sizes of the first breast lesion identification frame and the second breast lesion identification frame are the same, delete the The first breast lesion identification frame.
  9. 一种计算机设备,其特征在于,包括至少一个处理单元以及至少一个存储单元,其中,所述存储单元存储有计算机程序,当所述程序被所述处理单元执行时,使得所述处理单元执行权利要求1~5任一权利要求所述方法的步骤。A computer device, characterized by comprising at least one processing unit and at least one storage unit, wherein the storage unit stores a computer program, and when the program is executed by the processing unit, the processing unit executes the right The method of any one of claims 1 to 5 is required.
  10. 一种计算机可读存储介质,其特征在于,其存储有可由计算机设备执行的计算机程序,当所述程序在所述计算机设备上运行时,使得所述计算 机设备执行权利要求1~5任一所述方法的步骤。A computer-readable storage medium, characterized in that it stores a computer program that can be executed by a computer device, and when the program is run on the computer device, the computer device is caused to execute any of claims 1 to 5. Describe the steps of the method.
  11. 一种计算机程序产品,其特征在于,包括存储在计算机可读存储介质上的计算机程序,所述计算机程序包括程序指令,当所述程序指令被计算机设备执行时,使所述计算机设备执行权利要求1~5任一所述方法的步骤。A computer program product, characterized in that it includes a computer program stored on a computer-readable storage medium, the computer program includes program instructions, and when the program instructions are executed by a computer device, causes the computer device to execute the claims Steps 1 to 5 of the method.
PCT/CN2019/082687 2018-10-16 2019-04-15 Image-based breast lesion identification method and device WO2020077961A1 (en)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
CN201811201699.2 2018-10-16
CN201811201699.2A CN109363697B (en) 2018-10-16 2018-10-16 Method and device for identifying focus of breast image

Publications (1)

Publication Number Publication Date
WO2020077961A1 true WO2020077961A1 (en) 2020-04-23

Family

ID=65400521

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/CN2019/082687 WO2020077961A1 (en) 2018-10-16 2019-04-15 Image-based breast lesion identification method and device

Country Status (2)

Country Link
CN (1) CN109363697B (en)
WO (1) WO2020077961A1 (en)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109363697B (en) * 2018-10-16 2020-10-16 杭州依图医疗技术有限公司 Method and device for identifying focus of breast image
TWI769370B (en) * 2019-03-08 2022-07-01 太豪生醫股份有限公司 Focus detection apparatus and method thereof
CN110400302B (en) * 2019-07-25 2021-11-09 杭州依图医疗技术有限公司 Method and device for determining and displaying focus information in breast image
CN110930385A (en) * 2019-11-20 2020-03-27 北京推想科技有限公司 Breast lump detection and positioning method and device
CN111325743A (en) * 2020-03-05 2020-06-23 北京深睿博联科技有限责任公司 Mammary gland X-ray image analysis method and device based on combined signs
CN116258717B (en) * 2023-05-15 2023-09-08 广州思德医疗科技有限公司 Lesion recognition method, device, apparatus and storage medium

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106326931A (en) * 2016-08-25 2017-01-11 南京信息工程大学 Mammary gland molybdenum target image automatic classification method based on deep learning
CN106682435A (en) * 2016-12-31 2017-05-17 西安百利信息科技有限公司 System and method for automatically detecting lesions in medical image through multi-model fusion
CN107045720A (en) * 2017-05-04 2017-08-15 深圳硅基智能科技有限公司 Artificial neural network and system for recognizing eye fundus image lesion
US20170249739A1 (en) * 2016-02-26 2017-08-31 Biomediq A/S Computer analysis of mammograms
CN109363698A (en) * 2018-10-16 2019-02-22 杭州依图医疗技术有限公司 A kind of method and device of breast image sign identification
CN109363697A (en) * 2018-10-16 2019-02-22 杭州依图医疗技术有限公司 A kind of method and device of breast image lesion identification
CN109363699A (en) * 2018-10-16 2019-02-22 杭州依图医疗技术有限公司 A kind of method and device of breast image lesion identification
CN109447065A (en) * 2018-10-16 2019-03-08 杭州依图医疗技术有限公司 A kind of method and device of breast image identification

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2010003041A2 (en) * 2008-07-03 2010-01-07 Nec Laboratories America, Inc. Mitotic figure detector and counter system and method for detecting and counting mitotic figures
US9430829B2 (en) * 2014-01-30 2016-08-30 Case Western Reserve University Automatic detection of mitosis using handcrafted and convolutional neural network features
US10424069B2 (en) * 2017-04-07 2019-09-24 Nvidia Corporation System and method for optical flow estimation
CN107133933B (en) * 2017-05-10 2020-04-28 广州海兆印丰信息科技有限公司 Mammary X-ray image enhancement method based on convolutional neural network

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20170249739A1 (en) * 2016-02-26 2017-08-31 Biomediq A/S Computer analysis of mammograms
CN106326931A (en) * 2016-08-25 2017-01-11 南京信息工程大学 Mammary gland molybdenum target image automatic classification method based on deep learning
CN106682435A (en) * 2016-12-31 2017-05-17 西安百利信息科技有限公司 System and method for automatically detecting lesions in medical image through multi-model fusion
CN107045720A (en) * 2017-05-04 2017-08-15 深圳硅基智能科技有限公司 Artificial neural network and system for recognizing eye fundus image lesion
CN109363698A (en) * 2018-10-16 2019-02-22 杭州依图医疗技术有限公司 A kind of method and device of breast image sign identification
CN109363697A (en) * 2018-10-16 2019-02-22 杭州依图医疗技术有限公司 A kind of method and device of breast image lesion identification
CN109363699A (en) * 2018-10-16 2019-02-22 杭州依图医疗技术有限公司 A kind of method and device of breast image lesion identification
CN109447065A (en) * 2018-10-16 2019-03-08 杭州依图医疗技术有限公司 A kind of method and device of breast image identification

Also Published As

Publication number Publication date
CN109363697B (en) 2020-10-16
CN109363697A (en) 2019-02-22

Similar Documents

Publication Publication Date Title
WO2020077962A1 (en) Method and device for breast image recognition
CN109363699B (en) Method and device for identifying focus of breast image
CN109363698B (en) Method and device for identifying mammary gland image signs
WO2020077961A1 (en) Image-based breast lesion identification method and device
US10991093B2 (en) Systems, methods and media for automatically generating a bone age assessment from a radiograph
Li et al. DeepSEED: 3D squeeze-and-excitation encoder-decoder convolutional neural networks for pulmonary nodule detection
US10957079B2 (en) Systems and methods for automated detection of an indication of malignancy in a mammographic image
CN111325739B (en) Method and device for detecting lung focus and training method of image detection model
US10111632B2 (en) System and method for breast cancer detection in X-ray images
US9480439B2 (en) Segmentation and fracture detection in CT images
CN110046627B (en) Method and device for identifying mammary gland image
Banerjee et al. Automated 3D segmentation of brain tumor using visual saliency
US20170109880A1 (en) System and method for blood vessel analysis and quantification in highly multiplexed fluorescence imaging
JP2007307358A (en) Method, apparatus and program for image treatment
CN109461144B (en) Method and device for identifying mammary gland image
CN110689525A (en) Method and device for recognizing lymph nodes based on neural network
US11684333B2 (en) Medical image analyzing system and method thereof
CN112053325A (en) Breast mass image processing and classifying system
TWI587844B (en) Medical image processing apparatus and breast image processing method thereof
Montaha et al. A shallow deep learning approach to classify skin cancer using down-scaling method to minimize time and space complexity
US20210398282A1 (en) Digital tissue segmentation using image entropy
WO2020168647A1 (en) Image recognition method and related device
US20230334660A1 (en) Digital tissue segmentation and viewing
CN116229236A (en) Bacillus tuberculosis detection method based on improved YOLO v5 model
TW202042250A (en) Medical image analyzing system and method thereof

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: 19874375

Country of ref document: EP

Kind code of ref document: A1

NENP Non-entry into the national phase

Ref country code: DE

122 Ep: pct application non-entry in european phase

Ref document number: 19874375

Country of ref document: EP

Kind code of ref document: A1

122 Ep: pct application non-entry in european phase

Ref document number: 19874375

Country of ref document: EP

Kind code of ref document: A1