CN109242840A

CN109242840A - Detect breast image in beam-defining clipper region method, in breast image determine boundary method and Medical Devices

Info

Publication number: CN109242840A
Application number: CN201810996464.0A
Authority: CN
Inventors: 闫晶; 王汉禹; 牛杰
Original assignee: Shanghai United Imaging Healthcare Co Ltd
Current assignee: Shanghai United Imaging Healthcare Co Ltd
Priority date: 2018-08-29
Filing date: 2018-08-29
Publication date: 2019-01-18
Anticipated expiration: 2038-08-29
Also published as: CN109242840B

Abstract

The present invention relates to technical field of medical image processing, the method and Medical Devices on boundary are determined more particularly to a kind of method in beam-defining clipper region in detection breast image, in breast image.A kind of method in beam-defining clipper region in detection breast image, by carrying out the operation such as gradient transformation and binaryzation to breast image, and then the prime area in the breast image is determined based on the direct exposure area template of the binaryzation breast gradient image and binaryzation that obtain after binaryzation, and the linear feature in combination binaryzation breast gradient image, the precise positioning to beam-defining clipper zone boundary in breast image can be realized, can ensure that in turn will not remain beam-defining clipper region based on the breast image after the border cuts positioned, the complexity and difficulty of breast image post-processing is effectively reduced, promote the post-processing effect of breast image.

Description

It detects the method in beam-defining clipper region in breast image, determine boundary in breast image Method and Medical Devices

Technical field

The present invention relates to technical field of medical image processing, more particularly to beam-defining clipper region in a kind of detection breast image Method, in breast image determine boundary method and Medical Devices.

Background technique

It is clapped using digital galactophore X-ray shooting system (Digital Mammography System, abbreviation DMG) scanning When taking the photograph breast image, according to clinical demand, by selecting different plate for forcing to oppress corresponding shooting tissue, and pass through The size of mammary gland beam-defining clipper opening is adjusted, so that the size of mammary gland beam-defining clipper opening and used plate for forcing size dimension one It causes, it is subsequent to be cut according to the breast image that the coordinate pair of machine feedback is shot, to remove the beam-defining clipper in breast image Region.

But the defects of because of such as mechanical position tolerances, it will lead in the breast image after cutting, still retain limited Beam device region, and since beam-defining clipper is to belong to highly attenuating substance, so that remaining limit on the breast image after cutting can be made Beam device region is highly attenuating region, which will increase the complexity and difficulty of breast image post-processing, and then reduce The post-processing effect of breast image.

Summary of the invention

Based on this, it is necessary to provide it is a kind of detection breast image in beam-defining clipper region method, in breast image determine The method and Medical Devices on boundary can be applied to obtain captured by such as digital galactophore X-ray shooting system DMG scanning imagery On the breast image taken, accurately detected by the boundary to the beam-defining clipper region in breast image, in order to subsequent defeated Accurately beam-defining clipper region template out, so based on the beam-defining clipper region template can to the beam-defining clipper region in breast image into The operations such as row cutting can remain beam-defining clipper region in the breast image after effectively avoiding cutting and be formed by highly attenuating region While, the post-processing complexity and difficulty of breast image are reduced, the post-processing effect of breast image is promoted.

In view of the above technical problems, beam-defining clipper region in a kind of detection breast image is provided in embodiments herein Method, the method can include:

The breast image direct exposure area template of binaryzation corresponding with the breast image is obtained, the binaryzation directly exposes Light region template includes direct exposure area；

Gradient transformation and binaryzation are carried out to obtain binaryzation breast gradient image to the breast image；

Prime area is determined based on the binaryzation breast gradient image and the direct exposure area template of the binaryzation；

Obtain the linear feature on the binaryzation breast gradient image；And

Within the scope of the prime area on the breast image, linear feature is met into predetermined length threshold value, and position Region corresponding to pixel column and pixel column except the direct exposure area is as boundary corresponding to the beam-defining clipper.

In above-mentioned detection breast image in the method in beam-defining clipper region, by carrying out gradient transformation and two-value to breast image The operation such as change, and then based on the binaryzation breast gradient image obtained after binaryzation and the direct exposure area template of binaryzation come really Prime area in the fixed breast image, and the linear feature in binaryzation breast gradient image is combined, it can be realized to breast The precise positioning of beam-defining clipper zone boundary in image, and then can ensure that will not based on the breast image after the border cuts positioned Remain beam-defining clipper region, the complexity and difficulty of breast image post-processing is effectively reduced, promotes the post-processing of breast image Effect.

In an alternative embodiment, the step for obtaining the linear feature on the binaryzation breast gradient image Suddenly, it may include:

The linear feature on the binaryzation breast gradient image is obtained based on Hough transformation；Or

Based on the binaryzation breast gradient image, by being obtained respectively in pixel row direction and in pixel column side Upward projection value, to obtain the linear feature on the binaryzation breast gradient image.

In an alternative embodiment, the breast image includes breastwork side, opposite side, upside far from the breastwork side The downside and；

Wherein, when based on the binaryzation breast gradient image, by obtained in pixel row direction respectively and Projection value in pixel column direction, come when obtaining the linear feature on the binaryzation breast gradient image, in the pixel column Projection value on direction includes the first projection value corresponding with the upside and corresponding second projection value with the downside, in pixel Projection value on column direction includes third projection value corresponding with the opposite side, and the predetermined length threshold value includes the first pre- fixed length Spend threshold value, the second predetermined length threshold value and third predetermined length threshold value；And

Within the scope of the prime area on the breast image, linear feature is met into predetermined length threshold value, And it is located at region corresponding to pixel column and pixel column except the direct exposure area as corresponding to the beam-defining clipper The step of boundary, it may include:

In within the scope of the prime area on the breast image, first projection value is greater than described first and is made a reservation for Length threshold and be located at the direct exposure area except pixel column, second projection value be greater than second predetermined length Threshold value and be located at the direct exposure area except pixel column and the third projection value be greater than third predetermined length threshold value And it is located at the pixel column except the direct exposure area as boundary corresponding to the beam-defining clipper region.

In an alternative embodiment, the binaryzation breast gradient image may include the first subgraph, the second subgraph Picture, third subgraph and the 4th subgraph, the pixel column direction are consistent with the extending direction of the breastwork side；

Wherein, described that gradient transformation and binaryzation are carried out to obtain binaryzation breast gradient image to the breast image Step, comprising:

Gradient transformation and binaryzation are carried out to the breast image, and first subgraph is obtained based on first gradient threshold value Picture, and

It is obtained respectively based on the second Grads threshold and second subgraph corresponding to the upside and the downside institute The corresponding third subgraph and with the 4th subgraph corresponding to the opposite side；

Wherein, the first gradient threshold value is greater than second Grads threshold, and first subgraph is described for characterizing Breast image contour feature.

In an alternative embodiment, described directly to be exposed based on the binaryzation breast gradient image and the binaryzation The step of light region template determines prime area, it may include:

The prime area is determined based on first subgraph and the direct exposure area template of the binaryzation.

In an alternative embodiment, described to be based on the binaryzation breast gradient image, by obtaining respectively Projection value in pixel row direction and in pixel column direction, to obtain the spy of the straight line on the binaryzation breast gradient image The step of sign, it may include:

In acquisition on second subgraph in first projection value of pixel row direction；

In acquisition on the third subgraph in second projection value of pixel row direction；And

In acquisition on the 4th subgraph in the third projection value of pixel column direction；

Wherein, first projection value is on second subgraph in projection maximum in the pixel row direction Value, second projection value are on the third subgraph in projection value maximum in the pixel row direction, the third Projection value is on the 4th subgraph in projection value maximum in the pixel column direction.

In an alternative embodiment, described directly to be exposed based on the binaryzation breast gradient image and the binaryzation The step of light region template determines prime area, may also include that

The prime area is determined based on first subgraph and the direct exposure area template of the binaryzation and is somebody's turn to do The central point of prime area；

Wherein, central point and second Grads threshold based on the prime area, obtain second subgraph respectively Picture and the third subgraph；Second subgraph is used to characterize straight above the central point in the breast image Line feature, the third subgraph are used to characterize the linear feature being located at below the central point in the breast image.

In an alternative embodiment, first projection value is on second subgraph in the pixel column side Upward maximum projection value, second projection value are on the third subgraph in throwing maximum in the pixel row direction Shadow value, the third projection value are on the 4th subgraph in projection value maximum in the pixel column direction.

In an alternative embodiment, the first predetermined length threshold value is the prime area of the first predetermined ratio The value of upside or downside, the third predetermined length threshold value are the value of the prime area breastwork opposite side of the second predetermined ratio；

Wherein, the second predetermined length threshold value is equal with the first predetermined length threshold value.

In an alternative embodiment, above-mentioned method may also include that

Obtain the physical location of the beam-defining clipper；

When can not determine the boundary of any side directions in beam-defining clipper region based on the binaryzation breast gradient image, Then using the beam-defining clipper the side directions projected position of the physical location in breast image as boundary.

In an alternative embodiment, above-mentioned method may also include that

All boundaries based on acquisition cut the breast image, described in removing in the breast image Beam-defining clipper region.

In an alternative embodiment, above-mentioned method may also include that

According to projected position of the physical location of the beam-defining clipper on initial breast image, to the initial breast image It is cut to generate the breast image after cutting；

Wherein, the breast image that obtains is the breast image obtained after the cutting.

In an alternative embodiment, above-mentioned method may also include that

To the breast image carry out gradient transformation before, to the breast image carry out gray scale linearization process and/ Or median filter process.

In an alternative embodiment, above-mentioned method may also include that

Before carrying out gradient transformation to the breast image, bilinear interpolation processing is carried out to contract to the breast image Small image.

In an alternative embodiment, above-mentioned method may also include that

Arest neighbors interpolation processing is carried out with downscaled images to the direct exposure area template of the binaryzation；

Wherein, the direct exposure area template reduced scale of the binaryzation and the breast image reduced scale phase Together.

In an alternative embodiment, it is directly exposed based on the binaryzation breast gradient image and the binaryzation After region template determines the step of prime area, the method may also include that

In the prime area, benefit line is carried out to the linear feature partly covered by highly attenuating organization chart picture.

In an alternative embodiment, it is described to the linear feature partly covered by highly attenuating organization chart picture into Row mends the step of line, it may include:

Based on the breast image and scheduled first gray threshold, the first low ash degree template is obtained, wherein the height declines Subtract organization chart picture and is presented as first area in the first low ash degree template；

The first area of the first low ash degree template is expanded with predetermined expansion core to become second area and shape At the second low gray level areas template, wherein the area of first area described in the area ratio of the second area is big；

The part other than the prime area is removed in the second low ash degree template to form third low ash degree template, Correspondingly, the second area becomes third region, wherein the area in the third region is less than the face of the second area Product；And

When the extended line wait mend the linear feature of line is fallen in the third region, extend the linear feature.

In an alternative embodiment, it is directly exposed based on the binaryzation breast gradient image and the binaryzation Before region template determines the step of prime area, the method may also include that

Direct exposure area in the direct exposure area template of the binaryzation is corroded；

Wherein, corroded based on the first core or the second verification direct exposure area, and second core is greater than First core.

In an alternative embodiment, the direct exposure area in the direct exposure area template of the binaryzation The step of being corroded, it may include:

Based on the second predetermined gray threshold, the high gray scale template of the breast image is extracted；

It is straight in the direct exposure area template of the binaryzation to determine that the first gray areas in the high gray scale template accounts for The ratio of exposure area is connect, the gray scale of first gray areas is not less than the second predetermined gray threshold；

When the ratio is greater than predetermined ratio, to the direct exposure area in the direct exposure area template of the binaryzation Carry out the corrosion based on second core；Otherwise, the corrosion based on first core is carried out, or is not corroded.

The union of the binaryzation breast gradient image Yu the direct exposure area template of the binaryzation is sought, and simultaneously with this The boundary characteristic of collection determines the prime area；Or

The first boundary characteristic in the binaryzation breast gradient image is obtained, the direct exposure area of the binaryzation is obtained The second boundary feature of template, and the original area is determined with the union of first boundary characteristic and the second boundary feature Domain.

In another optional embodiment of the application, a kind of side that boundary is determined in breast image is additionally provided Method, the breast image are shot via the mammary machine for including beam-defining clipper, the method can include:

Obtain the physical location of beam-defining clipper when shooting；

Based on the breast image, binaryzation breast image and the direct exposure area template of binaryzation are obtained；The two-value Changing direct exposure area template includes direct exposure area；

Prime area is determined based on the binaryzation breast image and the direct exposure area template of the binaryzation；

Meet predetermined length along the extending direction acquisition on the boundary of the beam-defining clipper in the binaryzation breast image The outside at threshold value and the center far from the breast image does not include the linear feature of the direct exposure area；And

In the range of the prime area limits, according to the physical location of the linear feature and the beam-defining clipper in institute The relationship for stating the correspondence projected position of breast image, determines the boundary of the breast image.

In the above-mentioned method for determining boundary in breast image, by carrying out gradient transformation and binaryzation to breast image Deng operation, and then determined based on the direct exposure area template of the binaryzation breast gradient image and binaryzation that are obtained after binaryzation Prime area in the breast image, and the linear feature in binaryzation breast gradient image is combined, it can be realized to breast figure As the precise positioning on boundary, and then it can ensure that beam-defining clipper area will not be remained based on the breast image after the border cuts positioned Domain promotes the post-processing effect of breast image so that the complexity and difficulty of breast image post-processing is effectively reduced.

In an alternative embodiment, the binaryzation breast image is to pass through: first carrying out ladder to the breast image Degree transformation carries out binaryzation to the transformed image of gradient again and obtains.

In an alternative embodiment, the gradient transformation is the feature extracting method based on Sobel operator.

In an alternative embodiment, before gradient transformation further include:

The breast image is normalized, gradient transformation is carried out on image after normalization.

In an alternative embodiment, before the normalization further include:

Bilinear interpolation is carried out with downscaled images to the breast image, the normalization is that the image after diminution is enterprising Capable.

In an alternative embodiment, it is directly exposed described based on the binaryzation breast image and the binaryzation Before region template determines the step of prime area, the method also includes: in the direct exposure area template of the binaryzation Direct exposure area corroded；

Wherein, the corrosion is carried out based on the first core or the second core, and the second core is greater than the first core.

In an alternative embodiment, it is directly exposed described based on the binaryzation breast image and the binaryzation After region template determines the step of prime area, the method may also include that

In the prime area, to the institute partly covered on the binaryzation breast image by highly attenuating organization chart picture It states linear feature and carries out benefit line.

In an alternative embodiment, it is described on the binaryzation breast image by highly attenuating organization chart picture partly The step of linear feature of masking carries out benefit line, it may include:

Based on the breast image and the first predetermined gray threshold, the first low ash degree template is obtained, wherein described highly attenuating Organization chart picture is presented as first area in the first low ash degree template；

In an alternative embodiment, described in the range of the prime area limits, according to the linear feature With the physical location of the beam-defining clipper in the relationship of the corresponding projected position of the breast image, the side of the breast image is determined The step of boundary, comprising:

Compared with correspondence projected position of the physical location of the beam-defining clipper in the breast image, when the linear feature When closer to the center of the breast image, the linear feature is determined as to the boundary of the breast image；It otherwise, will be described The physical location of beam-defining clipper is determined as the boundary of the breast image in the correspondence projected position of the breast image.

In an alternative embodiment, it is described in the binaryzation breast image along the boundary of the beam-defining clipper It does not include the direct exposure that extending direction acquisition, which meets predetermined length threshold value and the outside at the center far from the breast image, In the step of linear feature in region:

The linear feature is obtained based on Hough transformation；Or

Based on the binaryzation breast image, by being obtained in pixel row direction and in pixel column direction respectively Maximal projection value, to obtain the linear feature.

One of the application it is optional be to additionally provide a kind of method that boundary is determined in breast image in embodiment, The breast image is shot via the mammary machine for including beam-defining clipper, the method can include:

Obtain the direct exposure area template of binaryzation of breast image and the breast image；

The physical location of beam-defining clipper when based on the shooting breast image, to the breast image and the breast image The direct exposure area template of binaryzation cut；

Based on the breast image after cutting, binaryzation breast image is obtained；

Based on the binaryzation breast image and the direct exposure area template of the binaryzation, prime area is determined；

In the range of the prime area limits, along the first extending direction of the beam-defining clipper, in the binaryzation The linear feature for meeting the first predetermined length threshold value is obtained in breast image, and the linear feature is determined as the beam-defining clipper edge First extending direction corresponding boundary on the breast image.

In the above-mentioned method for determining boundary in breast image, the physical location of acquired beam-defining clipper is first passed through to cream Room image and the direct exposure area template of binaryzation are cut, then again to after cutting breast image carry out gradient transformation and The operation such as binaryzation, based on the binaryzation breast gradient image obtained after binaryzation and the direct exposure area template of binaryzation come really Prime area in the fixed breast image, and the linear feature in binaryzation breast gradient image is combined, it can be realized to breast The precise positioning of image boundary, and then can effectively lower the difficulty and complexity of image procossing, while can ensure that again and being based on determining Breast image after the border cuts of position will not remain beam-defining clipper region, effectively to promote the post-processing effect of breast image.

In an alternative embodiment, the beam-defining clipper is also in the second extension side parallel with first extending direction To extension, the method may also include that

In the range of the prime area limits, along second extending direction of the beam-defining clipper, described two The linear feature for meeting the second predetermined length threshold value is obtained in value breast image, and the linear feature is determined as the limit beam Device corresponding boundary on the breast image along second extending direction.

In an alternative embodiment, the beam-defining clipper is also in the third extension side vertical with first extending direction To extension, the method may also include that

In the range of the prime area limits, along the third extending direction of the beam-defining clipper, described two The linear feature for meeting third predetermined length threshold value is obtained in value breast image, and the linear feature is determined as the limit beam Device corresponding boundary on the breast image along the third extending direction；

Wherein, the first predetermined length threshold value and the second predetermined length threshold value are equal.

In an alternative embodiment, it is directly exposed described based on the binaryzation breast image and the binaryzation Region template, before the step of determining prime area, the method may also include that

Bilinear interpolation is carried out with downscaled images to the breast image after cutting；And

Accordingly with the breast image, arest neighbors is carried out to the direct exposure area template of the binaryzation after cutting to insert Value is with downscaled images.

Global offset parameter based on detector carries out gray scale linearisation to the breast image.

The direct exposure area template of binaryzation of breast image and the breast image is obtained, the breast image has initial Boundary；The direct exposure area template of binaryzation includes direct exposure area；

Based on breast image, binaryzation breast image is obtained；

In the prime area, efficiency frontier linear feature of the breast image in all directions and described is obtained Boundary corresponding to initial boundary；Wherein, effective linear feature refers to obtaining from the binaryzation breast image, long The outside that degree is greater than the center of predetermined length threshold value and separate breast image does not have the straight line of the direct exposure area special Sign；

If the efficiency frontier linear feature has not been obtained in either direction, the breast image is in this direction Boundary corresponding to initial boundary is as efficiency frontier linear feature；

Based on all efficiency frontier linear features, the boundary of the breast image is determined.

In the above-mentioned method for determining boundary in breast image, by carrying out gradient transformation and binaryzation to breast image Deng operation, which is determined based on the direct exposure area template of the binaryzation breast gradient image and binaryzation that obtain after binaryzation Prime area in room image, and the initial boundary of the linear feature and breast image in binaryzation breast gradient image is combined, Comprehensive precise positioning to breast image boundary can be realized, and then can effectively lower the difficulty and complexity of image procossing, together When can ensure that again beam-defining clipper region will not be remained based on the breast image after the border cuts positioned, effectively to promote breast The post-processing effect of image.

In the alternative embodiment of the application, present invention also provides a kind of to obtain beam-defining clipper in breast image The method on boundary, it may include:

The direct exposure area template of binaryzation of breast image and the breast image is obtained, the breast image includes limit beam Device region and the highly attenuating low gray level areas lower than predetermined gray threshold organized the formation of, the direct exposure area of binaryzation Template includes direct exposure area；

Change of gradient and binary conversion treatment are carried out to the breast image, to obtain binaryzation breast gradient image；

Based on the binaryzation breast gradient image, the linear feature on each length direction of the beam-defining clipper is extracted；

When the outside at center of the linear feature far from the breast image, there is no the direct exposure areas, described Linear feature is unsatisfactory for predetermined length threshold value again, and when one end of the linear feature terminates at the low gray level areas, to this Linear feature is extended, and using the linear feature after the extension for meeting the predetermined length threshold value as the side of the beam-defining clipper Boundary.

In an alternative embodiment, the outside when center of the linear feature far from the breast image is not There are the direct exposure area, the linear features to be unsatisfactory for predetermined length threshold value again, and one end of the linear feature is whole When terminating in the low gray level areas, the step of extension to the linear feature, comprising:

Based on the breast image and scheduled first gray threshold, the first low ash degree template is obtained；Wherein, the height declines Subtract organization chart picture and is presented as first area in the first low ash degree template；

The first area of the first low ash degree template is expanded with predetermined expansion core to become second area and shape At the second low gray level areas template；Wherein, the second area is bigger than the first area；

When the direct exposure area, and institute is not present in the outside at center of the linear feature far from the breast image When stating the extended line of linear feature and falling in the second area, extend the linear feature.

In an alternative embodiment, described when the linear feature is unsatisfactory for predetermined length threshold value, to described straight The step of line feature is extended, it may include:

Prime area is limited in the second low ash degree template, wherein the prime area is by the binaryzation breast What gradient image and the direct exposure area template of binaryzation determined；

After the second low ash degree template is defined the prime area, institute is removed in the second low ash degree template The part other than prime area is stated to form third low ash degree template, correspondingly, the second area becomes third region, In, the area in the third region is less than the area of the second area；

When the direct exposure area and institute is not present in the outside at center of the linear feature far from the breast image When stating the extended line of linear feature and falling in the third region, extend the linear feature.

In an optional embodiment of the application, a kind of computer readable storage medium is additionally provided, is stored thereon There is computer program, realizes when which is executed by processor such as the step of above-mentioned any one the method.

In an optional embodiment of the application, a kind of Medical Devices are additionally provided, it may include:

Radiographic source includes the beam-defining clipper for limiting the X-ray boundary in the radiographic source for issuing X-ray；

Detector receives the X-ray through the breast being positioned between the radiographic source and detector, to generate breast figure Picture；

Memory is directly exposed for storing computer program, the breast image, the corresponding binaryzation of the breast image Region template and when shooting the breast image beam-defining clipper physical location the breast image projected position；With And

Processor, for according to the breast image, the direct exposure area template of the binaryzation and the beam-defining clipper Projected position, realize when executing the computer program such as the step in above-mentioned any one the method.

In above-mentioned Medical Devices, using processor to image acquisition units breast image collected, it is based on breast image And the direct exposure area template of the corresponding binaryzation of breast image, binaryzation is determined by gradient transformation and binarization operation Breast image, then determine prime area, and combine the linear feature in binary image, it can be realized to limiting beam in breast image The precise positioning of device zone boundary, and then can ensure that the breast image after cutting will not remain beam-defining clipper region, effectively to drop The complexity and difficulty of low breast image post-processing, promote the post-processing effect of breast image.

Detailed description of the invention

Fig. 1 is the flow diagram that the method in beam-defining clipper region in breast image is detected in one embodiment；

Fig. 2 a is the schematic diagram of beam-defining clipper zone boundary in one embodiment；

Fig. 2 b is the grey value characteristics figure of beam-defining clipper zone boundary shown in Fig. 2 a；

Fig. 3 a is to carry out pretreated breast image；

Fig. 3 b is the overall profile figure based on the extracted beam-defining clipper region of image shown in Fig. 3 a；

Fig. 3 c is the coboundary figure based on the extracted beam-defining clipper region of image shown in Fig. 3 a；

Fig. 3 d is the lower boundary figure based on the extracted beam-defining clipper region of image shown in Fig. 3 a；

Fig. 3 e is the left margin figure based on the extracted beam-defining clipper region of image shown in Fig. 3 a；

Fig. 4 a is the image in one embodiment before median filtering；

Fig. 4 b is the image carried out after median filtering to Fig. 4 a described image；

Fig. 5 a is the breast image there are overexposure phenomenon；

Fig. 5 b is based on the direct exposure area template of the extracted binaryzation of image shown in Fig. 5 a；

Fig. 5 c is based on the extracted high gray scale template of image shown in Fig. 5 b；

Fig. 6 a is original breast image after pretreatment in one embodiment；

Fig. 6 b is based on the extracted binaryzation profile diagram of image shown in Fig. 6 a；

Fig. 6 c is the direct exposure area template of binaryzation based on breast image shown in Fig. 6 a；

Fig. 6 d is acquired prime area after the union for seeking Fig. 6 b and Fig. 6 c；

Fig. 7 a is breast image after pretreatment in one embodiment；

Fig. 7 b is based on the extracted binaryzation profile diagram of image shown in Fig. 7 a；

Fig. 7 c is the direct exposure area template of binaryzation based on breast image shown in Fig. 7 a；

Fig. 7 d is acquired prime area after the union for seeking Fig. 7 b and Fig. 7 c；

Fig. 8 a is containing highly attenuating tissue and with the breast image of low gray level areas；

Fig. 8 b is the lower boundary straight line characteristic extracted based on image shown in Fig. 8 a；

Fig. 8 c is the low ash degree Prototype drawing extracted based on image shown in Fig. 8 a；

Fig. 8 d is the low ash degree Prototype drawing after being expanded based on image shown in Fig. 8 c；

Fig. 8 e be based on image shown in Fig. 8 d remove initial range outside low gray level areas after be formed by low ash degree template Figure；

Fig. 8 f is to be formed by lower boundary after carrying out benefit line to image shown in Fig. 8 b based on low ash degree template shown in Fig. 8 e Characteristic pattern；

Fig. 9 is the flow diagram that the method in beam-defining clipper region in breast image is detected in another embodiment；

Figure 10 is the flow diagram for determining the method on boundary in one embodiment in breast image；

Figure 11 is the flow diagram for determining the method on boundary in another embodiment in breast image；

Figure 12 is a kind of flow diagram for the method that boundary is determined in breast image in one embodiment；

Figure 13 is the flow diagram for obtaining the method on beam-defining clipper boundary in one embodiment in breast image；

Figure 14 is the structural schematic diagram of Medical Devices in one embodiment.

Specific embodiment

In order to make the objectives, technical solutions, and advantages of the present invention clearer, with reference to the accompanying drawings and embodiments, right The present invention is further elaborated.It should be appreciated that the specific embodiments described herein are merely illustrative of the present invention, and It is not used in the restriction present invention.

During carrying out image capture operations to the breast tissue of compressing, in captured breast image there may be Beam-defining clipper region, and since above-mentioned beam-defining clipper is highly attenuating substance, so corresponding model in above-mentioned breast image can be made Enclose it is interior there are highly attenuating region, the highly attenuating region will increase breast image post-processing complexity and difficulty, and then reduce cream The post-processing effect of room image.Traditional is cut by the breast image shot according to the coordinate pair of machine feedback, is come Highly attenuating region corresponding to the beam-defining clipper in breast image is removed, but the defects of because of such as mechanical position tolerances, will lead to In breast image after cutting, beam-defining clipper region is still remained with.

In order to solve the deficiency that tradition is cut according to machine feedback coordinate, by the breast figure in embodiments herein Beam-defining clipper region as in is detected and is positioned, and is cut the beam-defining clipper region in breast image according to location information Removal, can completely remove the beam-defining clipper region in breast image, and then effectively avoid because remaining beam-defining clipper in breast image Region and for subsequent image post-process caused by adverse effect.

Fig. 1 is the flow diagram that the method in beam-defining clipper region in breast image is detected in one embodiment.Such as Fig. 1 institute Show, a method of beam-defining clipper region in detection breast image, it may include following steps:

Step S1 carries out gradient transformation and binaryzation to breast image to obtain binaryzation breast gradient image.

For example, can select corresponding plate for forcing first according to clinical scene to squeeze the breast tissue of subject, such as may be selected Big plate for forcing, small plate for forcing press plate for forcing or oxter plate for forcing etc. to squeeze the breast tissue of subject, in order to subsequent right The opening of the mammary gland beam-defining clipper of x-ray imaging equipment should be adjusted, in the opening and selected pressure for making the mammary gland beam-defining clipper After compeling board size matching, image capture operations are carried out to the breast tissue of compressing, to obtain the breast image of subject.

After obtaining above-mentioned breast image, gradient transformation and binary conversion treatment can be successively carried out to the breast image, To obtain binaryzation breast gradient image, and then the behaviour such as beam-defining clipper region detection and positioning is carried out to breast image convenient for subsequent Make.

Step S2 obtains the direct exposure area template of binaryzation corresponding with above-mentioned breast image.

For example, the direct exposure region of binaryzation of the breast image can be obtained based on breast image acquired in step S1 Domain template.Wherein, the direct exposure area template of above-mentioned binaryzation includes direct exposure area.

It should be noted that above-mentioned, for the transformation of breast image progress gradient, simultaneously binaryzation obtains binaryzation breast gradient map It the step of picture, and the step of exposure area template direct based on breast image acquisition binaryzation, can carry out simultaneously, it can also be according to need It asks and successively carries out in a predetermined sequence, the embodiment of the present application has no specific limitation to two step along orders.

Step S3 determines prime area based on binaryzation breast gradient image and the direct exposure area template of binaryzation.

For example, can be based on acquired in binaryzation breast gradient image obtained in above-mentioned steps S1 and step S2 The direct exposure area template of binaryzation, to can determine in above-mentioned breast image, the region where beam-defining clipper region, i.e., initially Region.That is, the beam-defining clipper region on breast image is located in the prime area, which can be used for determining beam-defining clipper area Approximate range locating for domain is generally greater than beam-defining clipper region.

Step S4 obtains the linear feature on binaryzation breast gradient image.

For example, can be obtained using Hough transformation above-mentioned based on binaryzation breast gradient image acquired in step S1 Linear feature can also be projected by respectively in the pixel column direction and pixel row direction of the binaryzation breast gradient image It calculates, to obtain the projection value of the binaryzation breast gradient image on pixel column and pixel column, and then obtains above-mentioned straight line Feature.

Linear feature is met preset condition, and be located at straight in the range of the prime area on breast image by step S5 Region corresponding to the pixel column and pixel column except the range of exposure area is connect, as the corresponding boundary in beam-defining clipper region.

For example, can lead to based on identified prime area in step S3, and within the scope of the prime area of breast image Whether the projection value crossed on binaryzation breast gradient image acquired in judgment step S4 is greater than corresponding predetermined length threshold value, and after Continuous judgement is greater than pixel column or pixel column corresponding to the projection value of corresponding predetermined length threshold value, if positioned at above-mentioned direct exposure Except light region, or the outside by judging the rows and columns whether there is the modes such as direct exposure area, to carry out beam-defining clipper The detection of zone boundary, i.e., projection value is greater than corresponding predetermined length threshold value and the pixel column that is located at except direct exposure area and Region corresponding to pixel column, can also be within the scope of the prime area of breast image as the corresponding boundary in beam-defining clipper region, base The binaryzation breast gradient map of direct exposure area is not present in central point in the prime area, the side that will be far from the central point The region as corresponding to upper pixel column and pixel column is as the corresponding boundary in beam-defining clipper region；Correspondingly, may be based on Hough The acquired linear feature of transformation carries out the judgement of above-mentioned predetermined length threshold value, and then to obtain the corresponding side in beam-defining clipper region Boundary.

In an alternative embodiment, since when carrying out galactophore scanning, the side of neighbouring patient body generally will not Beam-defining clipper is opened, therefore there is no beam-defining clipper images in the breastwork side of scanning the formed breast image of mammary gland, if i.e. breast figure As being regular quadrangle, and it is generally breastwork side adjacent to patient body side (right side i.e. in present specification attached drawing), and The pixel column direction of the breast image is consistent with the extending direction of the breastwork side, i.e., left side refers to one opposite with its breastwork side Side, and the upper side and lower side is then remaining opposite two sides is then only needed to the upside, downside and left side of the breast image three totally The detection and positioning of side side extending direction progress beam-defining clipper zone boundary.

In an alternative embodiment, in above-mentioned binaryzation breast gradient image, the projection in pixel row direction Value can include: and corresponding first projection value in upside and corresponding with the downside of above-mentioned breast image second of above-mentioned breast image Projection value；Projection value in pixel column direction is then can include: third projection value corresponding with the left side of breast image.Accordingly , above-mentioned predetermined length threshold value then may include the first predetermined length threshold value corresponding with the upside of above-mentioned breast image, with it is upper State the corresponding second predetermined length threshold value in downside and the pre- fixed length of third corresponding with the left side of above-mentioned breast image of breast image Threshold value is spent, i.e. the first predetermined length threshold value is for judging whether the first projection value meets preset condition, the second predetermined length threshold value For judging whether the second projection value meets preset condition and third predetermined length threshold value for whether judging third projection value Meet preset condition.

In an alternative embodiment, can by judge pixel column and pixel column whether be located at direct exposure area it In, it is illustrated for the boundary in beam-defining clipper region to determine in breast image, i.e., as shown in Figure 1, above-mentioned detection breast image Step S5 in the method in middle beam-defining clipper region can include:

Within the scope of the prime area on above-mentioned breast image, the first projection value is greater than the first predetermined length threshold value and remote It is greater than the second predetermined length threshold there is no the pixel column of direct exposure area, by the second projection value from prime area central point side Value and the pixel column that direct exposure area is not present far from prime area central point side, and third projection value is greater than third Predetermined length threshold value and the pixel column that direct exposure area is not present far from prime area central point side, it is right on breast image Boundary of the region answered as corresponding side direction upper limit beam device region, it may be assumed that the first projection value is greater than the first predetermined length threshold value, And the pixel column of direct exposure area is not present far from prime area central point side, the corresponding region on breast image, Boundary as beam-defining clipper region on the breast image；Second projection value is greater than the second predetermined length threshold value, and separate The pixel column of direct exposure area is not present in prime area central point side, the corresponding region on breast image, as this The downside boundary in beam-defining clipper region on breast image；Third projection value is greater than third predetermined length threshold value, and far from original area The pixel column of direct exposure area is not present in domain central point side, the corresponding region on breast image, as the breast figure As the left border in upper beam-defining clipper region.In addition, subsequent can be considered to limit by the region corresponded on breast image other than boundary Beam device region, and then beam-defining clipper region template is inputted, the subsequent beam-defining clipper region template that is based on is to cut out breast image It the processing operations such as cuts, and then eliminates the adverse effect that beam-defining clipper region post-processes breast image.

In another alternative embodiment, since breast image under normal circumstances is positive square, i.e. breast image Upper and lower side edge length is identical, and above-mentioned the first predetermined length threshold value, the second predetermined length threshold value and third predetermined length threshold value The generally length value of its corresponding side predetermined ratio on prime area or breast image, so can be by the first predetermined length Threshold value setting is identical as the second predetermined length threshold value.In addition, the first above-mentioned predetermined length threshold value, the second predetermined length threshold value and Third predetermined length threshold value etc., may be defined as the length value of side predetermined ratio corresponding to above-mentioned prime area, also with further Promote the precision of beam-defining clipper region detection and positioning；Wherein, above-mentioned predetermined ratio can for such as 90%, 80%, 70% or 60% etc., and in order to further enhance the precision for detecting and positioning, generally above-mentioned predetermined ratio can be greater than 80%.

In an alternative embodiment, pretreatment operation first can be carried out to breast image, then carries out gradient transformation again And the operation such as binary conversion treatment, with the binaryzation breast gradient image in obtaining step S1.The binaryzation breast gradient image can Including the first subgraph, the second subgraph, third subgraph and the 4th subgraph, and first subgraph may include indicating above-mentioned The overall profile linear feature in beam-defining clipper region in breast image, the second subgraph may include indicating to limit beam in above-mentioned breast image The upper side edge boundary straight line feature in device region, third subgraph may include the downside for indicating beam-defining clipper region in above-mentioned breast image Side boundary straight line feature, the 4th subgraph may include indicating that the left side boundary straight line in beam-defining clipper region in above-mentioned breast image is special Sign；Wherein, it can be inserted by carrying out such as gray scale linearization process, median filtering, bilinearity to the acquired breast image of scanning Above-mentioned initial pictures are obtained in the pretreatment operations such as value diminution and/or normalized.

For example, as shown in Figure 1, step S1 in above-mentioned detection breast image in the method in beam-defining clipper region can include:

Based on breast image acquired in scanning mammary gland, Gradient Features and binaryzation behaviour are extracted based on different Grads threshold Make, to obtain the first above-mentioned subgraph, the second subgraph, third subgraph and the 4th subgraph respectively.Wherein, due to One subgraph is the overall profile linear feature that can indicate beam-defining clipper region and/or breast image, and the second subgraph, third Subgraph and the 4th subgraph can then indicate the linear feature of beam-defining clipper region respective side, therefore same Grads threshold can be used (such as the second Grads threshold) carries out gradient transformation to obtain above-mentioned second subgraph, third subgraph and the 4th subgraph, and can Above-mentioned first is obtained to carry out gradient transformation using the Grads threshold (such as first gradient threshold value) bigger than the second Grads threshold Subgraph.In an alternative embodiment, ladder can be normalized and then carried out to the breast image of acquisition Feature extraction and binarization operation are spent, the calculation amount of subsequent image processing is effectively reduced.

In an alternative embodiment, as shown in Figure 1, in above-mentioned detection breast image in the method in beam-defining clipper region Step S3 can include: above-mentioned prime area is determined based on the first subgraph and the direct exposure area template of binaryzation.And another It, can also be by obtaining the maximal projection value in pixel row direction on the second subgraph as above-mentioned in an optional embodiment The first projection value, obtain on third subgraph maximal projection value in pixel row direction as above-mentioned the second projection value and The maximal projection value in pixel column direction is obtained on the 4th subgraph as above-mentioned third projection value.

In an alternative embodiment, for such as because dosage is excessively high, and image is caused to be more than the reasons such as saturation gray scale, The caused unconspicuous breast image of boundary straight line feature, can be before above-mentioned determining prime area range, first for satiety Direct exposure area in the direct exposure area template of binaryzation corresponding to the breast image of sum carries out big nuclear corrosion, and is directed to Direct exposure area in the direct exposure area template of binaryzation corresponding to non-oversaturated breast image carries out small nuclear corrosion, With the further accuracy for extracting the boundary characteristic in beam-defining clipper region in breast image.Meanwhile for including that height declines Subtract the breast image of tissue, can first be directed to and be obtained after carrying out gradient to breast image and converting determining prime area range Direct exposure area in the direct exposure area template of the binaryzation taken carries out expansion process, based on acquired binaryzation breast Low gray level areas in gradient image carries out mending line processing for the linear feature partly covered by highly attenuating organization chart picture, Then judge whether linear feature meets preset condition in binaryzation breast gradient image again.

Fig. 2 a is the schematic diagram of beam-defining clipper zone boundary in one embodiment, and Fig. 2 b is beam-defining clipper regional edge shown in Fig. 2 a The grey value characteristics figure on boundary；Wherein, the horizontal axis in Fig. 2 b indicates the distance between breast image left side, and the longitudinal axis is expressed as ash Angle value.Referring to fig. 2 shown in a~2b, using the left side of breast image in Fig. 2 a as starting point, edge is parallel to breast image side up and down Extending direction on any straight line in breast image, is between the gray value G of the distance between left side D and current location The characteristic relation of curve L shown in Fig. 2 b.Referring to the extension trend of figure b curve L it is found that the gray value G in breast image exists Slope trend is nearby presented in 10 transition region of beam-defining clipper zone boundary in breast image, i.e. its first derivative can regard constant as, therefore Linear feature can be obtained for example, by using Sobel operator or Hough transformation etc. by First-order Gradient operator, carry out beam-defining clipper region The detection on boundary, specific:

Above-mentioned Sobel operator may include the operator of horizontal direction and the operator of vertical direction, and can be based on shown below Operator matrix table, horizontal direction formula and vertical direction formula is respectively adopted, to the horizontal direction and Vertical Square of breast image The extraction of straight line of beam-defining clipper zone boundary is carried out upwards.

Horizontal direction formula:

Vertical direction formula:

The amplitude formula of integral gradient are as follows: M (x, y) ≈ | g_x|+|g_y|

In the linear feature using the extracted horizontal direction of horizontal direction formula of above-mentioned Sobel operator, beam-defining clipper The gradient value of the coboundary (i.e. upper side edge) in region is positive value, and the gradient value of lower boundary (i.e. the lower side of breast image) is negative Value.In the linear feature using the extracted horizontal direction of vertical direction formula of above-mentioned Sobel operator, beam-defining clipper region The linear feature value of left margin (i.e. left side) is positive value；I.e. based on the feature for calculating each boundary using above-mentioned Sobel operator, energy Enough boundary characteristics for accurately extracting each side in beam-defining clipper region, further decrease the complexity and difficulty of image procossing, promote limit Beam device region detection and the precision of positioning.

3a~3e is directed to reference to the accompanying drawing retouches the linear feature on breast image using Sobel operator in detail It states:

Fig. 3 a is to carry out pretreated breast image；It includes beam-defining clipper that Fig. 3 b, which is extracted based on image shown in Fig. 3 a, The figure of the overall profile in region；Fig. 3 c is the coboundary figure based on the extracted beam-defining clipper region of image shown in Fig. 3 a；Fig. 3 d is Figure based on the extracted lower boundary including beam-defining clipper region of image shown in Fig. 3 a；Fig. 3 e is is mentioned based on image shown in Fig. 3 a The figure of the left margin including beam-defining clipper region taken；Wherein, Fig. 3 a is through normalizing pretreated breast image, Fig. 3 b~3e For binaryzation breast gradient image.As shown in Fig. 3 a~3e, it can be scanned first according to clinical demand using galactophore scanning equipment It is imaged and is pre-processed (such as normalized) and form breast image 21 as shown in Figure 3a afterwards；Secondly, using above-mentioned Sobel operator after being based on overall profile and the binarization operation that first gradient threshold value (T1) extracts the breast image 21, is formed Overall profile (i.e. entirety Sobel linear feature) 22 as shown in Figure 3b；Then, using the level side of above-mentioned Sobel operator To formula, after horizontal linear feature and binarization operation based on breast image 21 in the second Grads threshold (T2) extraction Fig. 3 a, And then coboundary linear feature 23 as shown in Figure 3c and lower boundary linear feature 25 as shown in Figure 3d are obtained respectively；Equally Ground breast image 21 in Fig. 3 a is extracted based on the second Grads threshold (T2) using the vertical direction formula of above-mentioned Sobel operator Vertical line feature and binarization operation after, and then obtain Fig. 3 e shown in left margin linear feature 26；Wherein, above-mentioned One Grads threshold can be greater than the second Grads threshold (i.e. T1 > T2), in order to obtain more accurately contour line, and it is specific Value can be set according to actual gray value of image and demand.In addition, in the straight line characteristic shown in above-mentioned Fig. 3 b~3e, When obtaining corresponding linear feature or profile, it can be possible to extract interfering line 24 as illustrated in figure 3 c simultaneously and do accordingly Disturb point (point of white i.e. shown in figure).

In an alternative embodiment, as shown in Fig. 3 a~3e, straight line is being carried out to breast image 21 shown in Fig. 3 a Before feature extraction, first the breast image can be normalized so that the equal normalizing of the pixel value on its image to 0~1 it Between in the range of, breast image 21 after then obtaining normalized shown in Fig. 3 a obtains linear feature in order to subsequent Threshold value reduces subsequent required process resource；Wherein, above-mentioned normalized may include first obtaining maximum in breast image Pixel value P_maxAnd minimum pixel value P_min, and by each pixel value P in the breast image_nSubtract above-mentioned minimum pixel value P_minAfterwards, divided by max pixel value P_maxWith minimum pixel value P_minDifference, that is, normalize after pixel value P_n'=(P_n-P_min)/ (P_max-P_min), n is positive integer, so that range of each pixel value between 0~1 on obtained image after normalization It is interior.

In an alternative embodiment, the above-mentioned method for detecting beam-defining clipper region in breast image, can also be based on upper The overall profile figure (that is, first subgraph) and the direct exposure area template of binaryzation in beam-defining clipper region shown in Fig. 3 b stated, To determine prime area of the beam-defining clipper region on breast image.Meanwhile the can be obtained based on the second subgraph shown in Fig. 3 c One projection value is obtained the second projection value based on third subgraph shown in Fig. 3 d and is obtained based on the 4th subgraph shown in Fig. 3 e Third projection value.For example, can using on the second subgraph in the maximum projection value of pixel row direction as the first above-mentioned projection value, Using on third subgraph in the maximum projection value of pixel row direction as the second above-mentioned projection value, by the 4th subgraph in picture The plain maximum projection value of column direction is as above-mentioned third projection value.

It in an alternative embodiment, can also be based on two corresponding to the first subgraph shown in above-mentioned Fig. 3 b and Fig. 3 a The direct exposure area template (not shown) of value first determines the central point of above-mentioned prime area, and is based on the central point and figure Breast image shown in 3a obtain respectively in turn the second subgraph shown in above-mentioned Fig. 3 c, third subgraph shown in Fig. 3 d and 4th subgraph shown in Fig. 3 e, and then it is convenient for the generation of each subgraph, reduce the difficulty and complexity of image procossing.Its In, the second above-mentioned subgraph can be used for characterizing the linear feature in breast image above the central point of prime area, on The third subgraph stated can be used for characterizing the linear feature in breast image below the central point of prime area, and above-mentioned the Four subgraphs can be used for characterizing the linear feature in breast image on the left of the central point of prime area.

For example, gradient change can be carried out just for the region between the upper side edge and central point on breast image shown in Fig. 3 a It changes and binary conversion treatment, and then obtains the second subgraph shown in Fig. 3 c, so that second subgraph only shows above-mentioned breast Linear feature on image near neighbouring upper side edge；Same lower side and central point on breast image shown in Fig. 3 a it Between region, carry out gradient transformation and binary conversion treatment to obtain third subgraph shown in above-mentioned Fig. 3 d；For Fig. 3 a institute Show the region between the left side and central point on breast image, carries out gradient transformation and binary conversion treatment to obtain above-mentioned figure 4th subgraph shown in 3e；It ensures that in this way, based on the maximum throwing in pixel row direction acquired in the second subgraph Pixel column corresponding to shadow value (i.e. the first projection value) is the coboundary in beam-defining clipper region the most doubtful, is based on third subgraph institute The pixel column corresponding to maximum projection value (i.e. the second projection value) in pixel row direction obtained is beam-defining clipper the most doubtful The lower boundary in region, and based on acquired in the 4th subgraph in pixel column direction maximum projection value (i.e. third projection value) Corresponding pixel column is the left margin in beam-defining clipper region the most doubtful；Namely the first above-mentioned projection value is in the second subgraph On in the maximum projection value of central point neighbouring in pixel row direction, the second projection value is on third subgraph in pixel row direction The upper neighbouring maximum projection value of central point, third projection value be on the 4th subgraph in central point neighbouring in pixel column direction most Big projection value.

In addition, after the boundary straight line feature for obtaining above-mentioned beam-defining clipper region, it can also be special to the boundary straight line extracted Sign is verified.For example, by judging that the outside of every extracted boundary straight line feature is directly exposed with the presence or absence of binaryzation Direct exposure area included by region template judges the straight line as boundary far from above-mentioned central point side with the presence or absence of straight Connect exposure area.If illustrating that extracted boundary is not true there are direct exposure area in the outside of boundary straight line feature Real beam-defining clipper boundary, then cast out the extracted boundary straight line feature, and can obtain corresponding side based on mechanical coordinate Boundary's feature；Conversely, then illustrating that extracted boundary straight line feature is true beam-defining clipper boundary.

In an alternative embodiment, the method in beam-defining clipper region may also include that in above-mentioned detection breast image

Beam-defining clipper mechanical coordinate is obtained from galactophore scanning equipment, and can not be obtained based on above-mentioned binaryzation breast gradient image When taking the boundary of any side in beam-defining clipper region, then position that can be corresponding on breast image by the beam-defining clipper mechanical coordinate As boundary, can based on the physical location of acquired beam-defining clipper, by beam-defining clipper the side directions physical location in cream Projected position in room image obtains the boundary on the upside, downside and three, left side direction of breast image as boundary, To guarantee to carry out accurate beam-defining clipper region detection and positioning on breast image；Finally owned again based on above-mentioned acquired Boundary (i.e. beam-defining clipper region is in upside, downside and three, left side direction boundary) ask union as boundary line to breast image It is cut, to ensure completely to wipe out beam-defining clipper region on breast image.

In another alternative embodiment, in order to reduce the difficulty and complexity of image procossing, promoted detection, positioning and The speed of image cropping, can be first according to above-mentioned beam-defining clipper mechanical coordinate to the above-mentioned direct exposure region of breast image and binaryzation Domain template is cut (i.e. for the first time cut) for the first time, to remove the beam-defining clipper region of mechanical coordinate feedback, and herein basis into The hand-manipulating of needle is detected and is positioned to the beam-defining clipper region for causing the image retention after cutting because of mechanical coordinate error, i.e., to cut Breast image afterwards is as the first image, the direct exposure area template of binaryzation after cutting as the first exposure area.So Afterwards, then to first image above-mentioned change of gradient and binary conversion treatment are carried out to obtain above-mentioned binaryzation breast gradient map Picture, and then continue above-mentioned every image processing operations, on the breast image after cutting to remaining beam-defining clipper region The operation such as detected, positioned and cut.It wherein, usually will be direct for the above-mentioned direct exposure area template of binaryzation The pixel value of exposure area is set to " 1 ", and the pixel value in other regions is set to " 0 ".

Further, beam-defining clipper region is not remained on the breast image after cutting if detecting；For example, cutting On breast image afterwards, the boundary in the residual beam-defining clipper region detected be boundary in received mechanical coordinate (examine The boundary in the beam-defining clipper region measured with the mechanical coordinate based on acquisition and the boundary of determination is consistent) when, illustrate it is above-mentioned just During secondary cutting, the beam-defining clipper region in breast image is wiped out completely, then can not had at this time after carrying out Continuous positioning and cutting etc. operation, and then promoted breast image in beam-defining clipper region detection, positioning and cut whole process into Degree and precision.

Fig. 4 a is the image in one embodiment before median filtering, and Fig. 4 b is after carrying out median filtering to Fig. 4 a described image Image.It, in an alternative embodiment, can be in order to promote the precision of beam-defining clipper region detection as shown in Fig. 4 a-4b Before carrying out gradient transformation to breast image 11 shown in Fig. 4 a, such as gray scale is carried out to breast image 11 shown in Fig. 4 a and is linearized The operation such as processing and median filter process, to reject bad line 12 in breast image 11 shown in Fig. 4 a, to obtain such as Fig. 4 b institute Show pretreated breast image 13；It is subsequent to continue the behaviour such as gradient transformation according to the breast image 13 after the optimization Make, then effectively avoids in breast image, interference the defects of because of such as bad line for beam-defining clipper region detection.In an example In, it can first pass through the whole global offset amount (Global for subtracting flat panel detector of breast image 11 shown in Fig. 4 a Offset after), the part less than 0 is set to 0, and then realize the gray scale linearization process carried out for breast image 11.

In another alternative embodiment, in order to further decrease process resource spent by image procossing, can into Before the transformation of row gradient, processing is zoomed in and out to above-mentioned breast image, and gradient is carried out based on the breast image after scaling processing Transformation and binarization operation, and then above-mentioned binaryzation breast gradient image is obtained, in order to carry out above-mentioned beam-defining clipper regional edge The detection and positioning on boundary, and the amplifying operation of corresponding ratio is carried out after determining beam-defining clipper zone boundary, to obtain breast figure As the actual boundary (i.e. beam-defining clipper template) in upper beam-defining clipper region.

For example, bilinear interpolation diminution processing can be carried out to above-mentioned breast image, by the breast image according to first Scale smaller is third image, and the breast image is such as reduced into its size 1/10th are used as third image；Meanwhile it can be right The direct exposure area template of the corresponding binaryzation of above-mentioned breast image carries out arest neighbors interpolation diminution processing, by the binaryzation Direct exposure area template also according to above-mentioned first scale smaller is and the direct exposure region of third image size matched second Domain template.It is subsequent then to continue the operations such as above-mentioned gradient transformation and binaryzation, and then obtain according to above-mentioned third image Above-mentioned binaryzation breast gradient image.Due to carrying out the acquired third image of above-mentioned image down processing, the second directly exposure Breast image more corresponding than its and the direct exposure area template of binaryzation are much smaller respectively for light region template, so on carrying out Many process resources will be saved by stating many image procossings, while can also greatly improve the speed of processing.

It in another alternative embodiment, can be to the limit of subsequent progress since breast image is there are when overexposure region The inspection and positioning in beam device region adversely affect, so settable one second picture bigger than above-mentioned first pixel threshold Plain threshold value, and high gray scale region template, and the high gray scale region template can be obtained based on second pixel threshold and breast image It is subsequent directly to be exposed by judging that the high gray areas in acquired high gray scale template accounts for binaryzation including high gray areas Whether the ratio of direct exposure area is greater than default ratio in region template, to determine that breast image is existing with the presence or absence of overexposure As.Wherein, the second pixel threshold can be substantially larger than the first pixel threshold；For example, the first pixel threshold can be set as 50~100, and Two pixel thresholds then may be configured as 40000~60000.

Further, if ratio between above-mentioned high gray areas and direct exposure area is greater than default ratio, to can Big nuclear corrosion is carried out to the direct exposure area in the direct exposure area template of binaryzation, and if high gray areas and directly exposure Ratio between region is less than or equal to default ratio, then can be to the direct exposure area in the direct exposure area template of binaryzation Carry out small nuclear corrosion (can also be without processing).It then, then can be based on the direct exposure area of binaryzation after above-mentioned corrosion treatment Template continues the determination of subsequent prime area range.Wherein, the kernel matrix of above-mentioned big nuclear corrosion is greater than the interior of small nuclear corrosion Nuclear matrix.

Fig. 5 a is the breast image there are overexposure phenomenon, and Fig. 5 b is straight based on the extracted binaryzation of image shown in Fig. 5 a Exposure area template is connect, Fig. 5 c is based on the extracted high gray scale template of image shown in Fig. 5 a；Wherein, Fig. 5 a can be at through diminution Breast image after reason, and Fig. 5 b can be to be reduced treated the direct exposure area template of binaryzation, which can also be base In the direct exposure area template of binaryzation that Fig. 5 a is obtained.The case where there are overexposures for breast image, due to overexposure meeting It causes direct exposure area excessive, or even beam-defining clipper region can be covered, and then will lead to and detected based on linear feature Cast out on beam-defining clipper boundary out.As shown in Fig. 5 a~5c, in an alternative embodiment, for breast figure shown in Fig. 5 a As due to causing the area of direct exposure area 31 excessive there are overexposure, it is possible to influence whether beam-defining clipper zone boundary It extracts and verifies；It is obtained shown in Fig. 5 b so can first extract based on the direct exposure area 31 on breast image shown in Fig. 5 a High gray scale template 33 shown in the direct exposure area template 32 of binaryzation and Fig. 5 c, i.e., reach saturation by gray scale in breast image The number of pixels in the region of gray scale, and the number of pixels of direct exposure area in breast image is extracted, then pass through judgement ash Degree reaches the number of pixels in region of saturation gray scale and whether the ratio between the number of pixels of direct exposure area reaches the default of exposure Fig. 5 b (is compared) by value with Fig. 5 c；For example judge that gray scale reaches the number of pixels and directly exposure in the region of saturation gray scale Whether the ratio between the number of pixels in region reaches 90%, illustrates that there are overexposure phenomenons in the breast image if reaching, otherwise There is no overexposure phenomenons.Then, when in breast image there are when overexposure phenomenon, due to overexposure meeting so that subsequent be extracted Direct exposure area it is more larger than actual direct exposure area, that is, exist in the direct exposure area template of binaryzation of overexposure Direct exposure area can become larger compared to actual direct exposure area, thus can be to subsequent beam-defining clipper boundary characteristic Extraction adversely affects.For example, extracted correct beam-defining clipper can be made since extracted direct exposure area is excessive The outside of the boundary straight line in region can still have direct exposure area, so that the boundary straight line can be verified subsequent It is rejected in journey.So at this point, can there are the direct exposure area templates of the binaryzation of overexposure to carry out big nuclear corrosion for above-mentioned, It will not be cast out with the boundary in the correct beam-defining clipper region for ensuring to be detected.

In another alternative embodiment, in order to promote the accuracy that beam-defining clipper zone boundary linear feature verifies, also Small nuclear corrosion can be carried out for the direct exposure area template of binaryzation there is no overexposure, be further continued for prime area range really It is fixed；Wherein, " core " in above-mentioned big nuclear corrosion and small nuclear corrosion can be any shape and size, possesses one and individually determines The reference point (anchorpoint, i.e. anchor point) that justice comes out.For example, the core K2 of core K1 and small nuclear corrosion in big nuclear corrosion are One centre is with reference point and filled square table or discoid table, and the matrix of core K1 is greater than the square of core K2 Battle array；If core K1 can be the matrix of 5*5, core K2 can be the matrix of 3*3.

In an alternative embodiment, as shown in Figure 1, in above-mentioned detection breast image in the method in beam-defining clipper region Other modes can be used also to determine prime area in step S3；For example, above-mentioned binaryzation breast gradient image can be obtained first First boundary characteristic, and the second boundary feature of the direct exposure area template of binaryzation is obtained, and with first boundary characteristic Above-mentioned prime area is determined with the union of the second boundary feature；Alternatively, can also be by first seeking binaryzation breast gradient image With the union between the direct exposure area template of binaryzation, such as Fig. 7 b and Fig. 7 c are sought into union and obtain Fig. 7 d, it is subsequent can base Boundary characteristic is being sought in the range that Fig. 7 d is limited, the prime area of above-mentioned beam-defining clipper is then determined with boundary characteristic.

In an alternative embodiment, in order to enable the positioning in beam-defining clipper region is more accurate, above-mentioned Sobel can be used Operator integral gradient amplitude characteristic and the direct exposure area template of binaryzation determine the prime area in beam-defining clipper region.Specifically , it can first determine the integral gradient characteristic image of breast image and the direct respective boundary position of exposure area template of binaryzation, Such as can respectively it be carried out horizontal and perpendicular by the integral gradient characteristic image and the direct exposure area template of binaryzation to breast image Histogram to projection, and in the vertical direction obtain floor projection the top first projection value be " 1 " position as upper Boundary, and first projection value for obtaining floor projection bottom in the vertical direction is the position of " 1 " as lower boundary, simultaneously Obtaining the first projection value of the vertical projection value leftmost side in the horizontal direction is the position of " 1 " as left border, and in level The initial position on the right side (i.e. breastwork side) of figure is as right side boundary on direction, and then obtains beam-defining clipper in breast image respectively The profile in region and the profile of direct exposure area are then based on the profile in above-mentioned beam-defining clipper region and the wheel of direct exposure area Wide union is as final prime area.The process of union is sought as in first throwing of all directions of the above two images Shadow value be " 1 " position in, take far from image center first projection value be " 1 " position be the side it is corresponding initially Range boundary.

Fig. 6 a is breast image after pretreatment in one embodiment, and Fig. 6 b is extracted based on image shown in Fig. 6 a Binaryzation profile diagram, Fig. 6 c are the direct exposure area template of binaryzation based on breast image shown in Fig. 6 a, and Fig. 6 d is to seek figure Acquired prime area after the union of 6b and Fig. 6 c.It, can be based on cream after pretreatment shown in Fig. 6 a as shown in Fig. 6 a~6d Room image 41, the overall profile 42 for converting linear feature by obtaining gradient shown in Fig. 6 b after gradient transformation and binary conversion treatment (i.e. the first subgraph)；The profile 43 that the direct exposure area of binaryzation shown in Fig. 6 c is obtained based on Fig. 6 a, then seeks Fig. 6 b Prime area 44 shown in Fig. 6 d is finally determined with the union of Fig. 6 c.Wherein, by being formed after being cut using mechanical coordinate Beam-defining clipper region is not present in Fig. 6 a, so obtaining final prime area is entire image.

Fig. 7 a is breast image after pretreatment in one embodiment, and Fig. 7 b is extracted based on image shown in Fig. 7 a Binaryzation profile diagram, Fig. 7 c are the direct exposure area template of binaryzation based on breast image shown in Fig. 7 a, and Fig. 7 d is to seek figure Acquired prime area after the union of 7b and Fig. 7 c.It, can be based on cream after pretreatment shown in Fig. 7 a as shown in Fig. 7 a~7d Room image 61 obtains the monolithic wheel of binaryzation breast gradient breast image shown in Fig. 7 b by gradient transformation and binary conversion treatment Exterior feature 62；The profile 63 that direct exposure area shown in Fig. 7 c is obtained based on Fig. 7 a, the union for then seeking Fig. 7 b and Fig. 7 c are come most Prime area 64 shown in Fig. 7 d is determined eventually.Wherein, since there are beam-defining clipper regions in breast image 61, so obtain final Prime area 64 be not entire image, but the region in white box as shown in Figure 7 d.

In an alternative embodiment, in order to eliminate low gray level areas middle short line for above-mentioned beam-defining clipper region detection and The adverse effect of positioning can obtain before judging whether the linear feature on binaryzation breast gradient image meets preset condition Take the low gray level areas on above-mentioned breast image, and based on acquired prime area range, to pass through the low gray level areas and Line segment within the scope of above-mentioned prime area carry out mend line operation, with based on mend line after binaryzation breast gradient image after The continuous judgement for carrying out linear feature.Wherein, above-mentioned low gray level areas generally refer in breast image pixel value be less than it is preset The region of first pixel threshold, and first pixel threshold can be for according to the pixel value progress of each region in practical breast image It is arranged and is adaptively adjusted, as long as it can be used to distinguish the region for causing broken string on the breast image.

It is formed by breast image in tissue for example, covering based on beam-defining clipper, due to beam-defining clipper area in breast image Domain and the tissue region of high density low ash degree are there are Chong Die, so that the boundary in the beam-defining clipper region of overlapping can not extract Come, and then the linear feature on beam-defining clipper boundary corresponding to the side can be made shorter, to be unable to reach threshold range, causes The linear feature may be rejected.It can be by based on above-mentioned breast figure to above-mentioned breast image for this situation As extracting low gray level areas, then judge whether the extending direction of the boundary sections in above extracted beam-defining clipper region passes through this Low gray level areas, and direction is extended along in prime area range to the boundary sections for passing through the low gray level areas and carries out benefit line, To extend the length of above-mentioned boundary sections.

In addition, in order to further enhance the efficiency of image procossing, it is only low in prime area to being located at when carrying out benefit line Gray areas carries out benefit line, and for the region being located at other than prime area, even if it is low gray level areas, also without mending line Operation.Wherein, above-mentioned low gray level areas refers to grey scale pixel value lower than preset threshold (for example, first predetermined gray threshold) Pixel region, and the preset threshold is less than the preset threshold of high gray areas signified in the application other embodiments.It is specific:

Fig. 8 a is containing highly attenuating tissue and with the breast image of low gray level areas, and Fig. 8 b is based on image shown in Fig. 8 a And the lower boundary straight line characteristic extracted, Fig. 8 c are the low ash degree Prototype drawing extracted based on image shown in Fig. 8 a, Fig. 8 d is base The image shown in Fig. 8 c expanded after low ash degree Prototype drawing, Fig. 8 e be based on image shown in Fig. 8 d remove initial range outside Low ash degree Prototype drawing is formed by after low gray level areas, Fig. 8 f is based on low ash degree template shown in Fig. 8 e to image shown in Fig. 8 b Lower boundary characteristic pattern is formed by after carrying out benefit line.As shown in Fig. 8 a~8e, in the present embodiment, contain for shown in Fig. 8 a The breast image 71 of highly attenuating tissue, since there are highdensity low gray level areas for its lower side, so the lower boundary extracted is straight Line feature obtains in lower boundary linear feature Figure 72 shown in Fig. 8 b, and lower boundary linear feature is partly hidden by low gray level areas It covers, so need to carry out Fig. 8 b to mend line operation.Specifically, low ash degree mould first can be obtained based on breast image 71 shown in Fig. 8 a Plate 73 (i.e. image shown in Fig. 8 c), and after carrying out expansion process to above-mentioned low ash degree template 73, it is swollen described in Fig. 8 d to obtain Low ash degree template 74 after swollen, the low ash degree template 74 is more mellow and fuller relative to above-mentioned low ash degree template 73, and region area Can be larger, so that the operation of subsequent benefit line is more accurate；Then, prime area range in above-mentioned low ash degree template 74 is removed Behind region in addition, low ash degree template 75 final shown in Fig. 8 e can be formed, and then the low gray area in prime area has been determined Domain (i.e. white area 77 shown in Fig. 8 e), in order to guarantee to be mended just for the low gray level areas within the scope of prime area Low ash degree template outside the range of prime area can be set to " 0 " before carrying out benefit line by line；Finally, based on final shown in Fig. 8 e After low ash degree template 75 carries out benefit line to image 72 shown in Fig. 8 b, lower boundary linear feature Figure 76 as illustrated in fig. 8f is formed. Wherein, comparison diagram 8b and Fig. 8 f is it is found that the line segment for partly being covered in image 72 and being formed by low gray level areas carries out After mending line, the line segment (i.e. white line segment 78) being formed by image 7f, length is obviously increased, and then can effectively be kept away Beam-defining clipper zone boundary line segment correct boundary that is too short, and causing caused by exempting from because of the highly attenuating low gray level areas for organizing to generate Line segment is rejected.

In another alternative embodiment, for shown in Fig. 8 a, there are the images of highly attenuating substance, can first carry out base It is straight to go out coboundary linear feature, lower boundary linear feature and left margin using Sobel operator extraction for central point in prime area After line feature and binarization operation, line operation is mended carrying out the low ash degree as shown in Fig. 8 b~8f, and after mending line operation, into The judgement of row beam-defining clipper zone boundary (for details, reference can be made to Fig. 3 a~3e and related elaborations).

Further, after getting boundary corresponding to above-mentioned each side, the direct exposure area of binaryzation can be based on Template judges the position of acquired maximal projection value, that is, judge pixel column where the maximal projection value of the acquisition or Pixel column whether there is direct exposure region in the side far from above-mentioned prime area central point；Direct exposure area if it does not exist, Then illustrate that acquired boundary is real boundary, it is subsequent that coordinate conversion is carried out based on the real boundary, to breast figure As being cut, and then remove the beam-defining clipper region in the side directions.And if the picture where the maximal projection value in the acquisition Plain row or pixel column, in the side far from above-mentioned prime area central point, there are when direct exposure region, then illustrate above-mentioned boundary inspection There may be mistakes for survey, therefore the pixel column where the maximal projection value of the acquisition is given up, and can will correspond to the direction On mechanical coordinate corresponding to boundary or the side of breast image as real boundary, carry out subsequent trimming operation, into And area-of-interest is extracted, convenient for carrying out post-processing operation to the breast image after cutting.In addition, in an optional implementation In example, if the projection value of edge lengths 85% on the upside of all such as larger than breast images can not be found in above-mentioned all pixel columns, It can be directly using the corresponding side of breast image as the corresponding boundary in beam-defining clipper region.

Further, when the side using breast image is as real boundary, it can first judge above-mentioned mechanical coordinate In corresponding side whether be located on breast image, if side corresponding in above-mentioned mechanical coordinate is located at breast image On, then using side corresponding in mechanical coordinate as true boundary；Otherwise, then using the side of breast image as real side Boundary.Meanwhile in same side directions, when using true side acquired in above-mentioned Sobel operator (such as Fig. 3 c~3e institute shape At side border) or mechanical coordinate in side can be only formed a part on respective side boundary when, in combination with above-mentioned true Side and mechanical coordinate form the above-mentioned same side border, also in combination with true side, mechanical coordinate and breast image side While to form the above-mentioned same side border, equally also formed in combination with mechanical coordinate and breast image side above-mentioned same One side border, i.e., the same side border are that different modes is combined to be determined；Wherein, it at same position, uses The side of above-mentioned true side, the side in mechanical coordinate and breast image is successively reduced as the priority of side border.

Fig. 9 is the flow diagram that the method in beam-defining clipper region in breast image is detected in another embodiment.Such as Fig. 9 institute It states, in conjunction with real image process flow, the method in beam-defining clipper region in breast image is detected in the present embodiment can include:

Step S02 obtains breast image, and can obtain the images such as the direct exposure area template of binaryzation based on breast image Parameter, and step S04 can be continued.Wherein, the direct exposure area template of the binaryzation may include direct exposure area.

Step S04 obtains beam-defining clipper mechanical coordinate corresponding with breast image in step S02, in order to reduce subsequent image The calculation amount and complexity of processing, can also according to the beam-defining clipper mechanical coordinate in step S02 breast image and binaryzation it is straight Connect the progress first time trimming operation such as exposure area template.

Step S06 can be to step in order to remove the bad line in breast image and reduce the calculation amount of subsequent image processing Breast image after cutting in S04 carries out such as gray scale linearisation, median filtering (removing bad line) and/or bilinearity difference and reduces Pretreatment operations such as (reduce calculation amount), and step S10, step S12 and step can be carried out simultaneously or successively according to predesigned order S14.Wherein, can be after the global offset amount (nGlobal offset) by the way that breast image integrally to be subtracted to flat panel detector, it will Part of the pixel value less than 0 is set to 0, to realize the gray scale linearization process to breast image.

Step S10, acquired breast image, extracts the low ash degree of the breast image after being pre-processed based on step S06 Template, and step S24 can be continued, it carries out mending line operation in order to subsequent.

Step S12, acquired breast image after being pre-processed based on step S06 can carry out gradient to the breast image Transformation and binary conversion treatment to obtain binaryzation breast gradient image, and can proceed with step S22.Wherein, in the step 12 Binaryzation breast gradient image may include whole boundary characteristic binaryzation breast gradient image and each side border feature binaryzation Breast gradient image.For example, Sobel operator can be used for the breast image with breastwork side, gradient transformation and two-value are carried out Change operation, and then obtain Sobel entirety boundary characteristic binaryzation breast gradient image (i.e. the first subgraph), and is upper and lower, left The Sobel boundary characteristic binaryzation breast gradient image (i.e. second and third, four subgraphs) of three sides.

Step S14, acquired breast image, extracts high gray scale template, and can continue after being pre-processed based on step S06 Step S16, in order to subsequent progress overexposure light processing.

Step S16, based on height acquired in the direct exposure area template of binaryzation acquired in step S02 and step S14 Whether gray scale template, the ratio that high gray areas of the judgement at this time in breast image accounts for direct exposure area are greater than predetermined value (such as 90%)；If being greater than predetermined value, that is, there is supersaturation, therefore step S18 can be continued, otherwise, then can continue step 20.

Step S18 carries out big nuclear corrosion to the direct exposure area template of binaryzation acquired in step S02, and can continue Step S22.

Step S20 carries out small nuclear corrosion to the direct exposure area template of binaryzation acquired in step S02, and can continue Step S22.Wherein, the nuclear matrix of big nuclear corrosion is greater than the nuclear matrix of small nuclear corrosion.

Step S22, based on the direct exposure area template of binaryzation after above-mentioned steps S18 or step S20 corrosion treatment, with And acquired whole boundary characteristic binaryzation breast gradient image in step S12, to determine the breast figure after step S04 is cut The central point of the prime area range in beam-defining clipper region and the prime area as in, and step S24 can be continued.

Step S24 is right based on the prime area range in the extracted low ash degree template of step S10 and step S22 Acquired each side border feature binaryzation breast gradient image carries out mending line operation in step S12, and can continue step S26。

Step S26 is obtained based on each side border feature binaryzation breast gradient image acquired after step S24 benefit line It takes on each side border feature binaryzation breast gradient image, the maximal projection value on pixel column, column direction, and can continue to walk Rapid S28.

Step S28 can be sentenced for the boundary characteristic of any side directions based on maximal projection value acquired in step S26 Whether the maximal projection value of breaking is greater than corresponding predetermined length threshold value, and to pixel column or pixel column far from acquired in step S22 Central point side whether there is direct exposure area；If maximal projection value is greater than corresponding predetermined length threshold value, and to picture There is no direct exposure areas for the side of plain row or pixel column far from central point acquired in step S22, then can continue step S30；Otherwise, then it can proceed with step S32.Wherein, above-mentioned predetermined length threshold value can be the length of corresponding initial areas side The certain percentage length of degree, such as edge direction on the upside of prime area, predetermined length threshold value can be 85%*L, and L is original area For domain upper side edge in the projection value of line direction, L is positive integer.

Step S30 will be located within the scope of prime area acquired in step S22, and projection value is greater than predetermined length threshold value, Direct exposure area is not present in (side i.e. far from above-mentioned central point) on the outside simultaneously, with pixel column corresponding to the projection value Or pixel column, as on respective side direction, boundary corresponding to beam-defining clipper region, and step S34 can be continued.

Step S32, by mechanical coordinate acquired in step S04, the boundary on respective side direction is as respective side Boundary corresponding to beam-defining clipper region on direction, and step S34 can be continued.

Step S34 seeks the union on boundary acquired in step S30 and step S32, and using the union as beam-defining clipper template It is exported, and second can be carried out to the breast image after cutting in step S04 based on the beam-defining clipper template and cut, and then gone Except first time cuts remaining beam-defining clipper region in rear udder attachment image.

It in an alternative embodiment, then will be in step when step S06 carries out the diminution of bilinearity difference to breast image The arest neighbors diminution for carrying out corresponding ratio to the direct exposure area template of binaryzation after rapid S02 and before step S16 is handled；Together When, in step S34, after obtaining beam-defining clipper template, breast image is carried out second before cutting, it be to beam-defining clipper template The amplification of corresponding ratio is carried out, so that acquired beam-defining clipper template is suitable for the breast image after cutting for the first time.

It in another alternative embodiment, can also be by the image of breast image between above-mentioned step S06 and step S12 Tonal range normalizes between 0~1, i.e., after breast image being normalized, is carrying out gradient transformation and binaryzation Operation, to reduce the complexity of subsequent image processing.

The method that beam-defining clipper region in breast image is detected in the present embodiment, can be applied to the image procossing stream of breast image During journey (Pipeline), i.e., after the detection of the direct exposure area of breast image, first according to the mechanical coordinate of feedback The first cutting in beam-defining clipper region is carried out to breast image, and after the detection of the direct exposure area of breast image, based on this The relevant technologies content in breast image in the method in beam-defining clipper region is detected in embodiment, to the breast image after first cut The detection and accurate positioning in the remaining beam-defining clipper region of further progress, and when detecting remaining beam-defining clipper region to first Breast image after cutting is cut again, and with the mechanical coordinate that effectively prevent due to above-mentioned feedback, there are errors, is caused Remain beam-defining clipper region on breast image after first cutting, and influences the post-processing operation of breast image, and then promotion feeling The accuracy rate of interest extracted region and the quality for post-processing image.

Figure 10 is the flow diagram for determining the method on boundary in one embodiment in breast image.As shown in Figure 10, In an alternative embodiment, the method that present invention also provides a kind of determines boundary in breast image, the breast image It can be for via the medical image for the mammary machine shooting for including beam-defining clipper, this method be can comprise the following steps that

Step S111, the physical location of beam-defining clipper when obtaining shooting.

Step S112 is based on breast image, obtains binaryzation breast image and the direct exposure area template of binaryzation respectively, And the direct exposure area template of binaryzation may include direct exposure area.

Specifically, can first to breast image carry out the diminution of bilinearity difference processing after, then to the breast image after diminution into Row normalized, and gradient transformation is carried out to breast image using Sobel operator, continue to the transformed breast image of gradient Binaryzation is carried out, above-mentioned binaryzation breast image is then obtained.Meanwhile it can be to the above-mentioned direct exposure area template of binaryzation After the closest interpolation diminution processing for carrying out ratio corresponding with above-mentioned bilinearity difference diminution, and then the binaryzation reduced is straight Connect exposure area template.

Step S113 determines prime area based on binaryzation breast image and the direct exposure area template of binaryzation.

Specifically, the second predetermined gray threshold can be first based on, after the high gray scale template for extracting above-mentioned breast image, determining should The first gray areas in high gray scale template accounts for the ratio of the direct exposure area in the direct exposure area template of above-mentioned binaryzation, And when the ratio is greater than predetermined ratio, base is carried out to the direct exposure area in the above-mentioned direct exposure area template of binaryzation In the corrosion of the second core；Otherwise, the direct exposure area in the above-mentioned direct exposure area template of binaryzation is carried out based on the The corrosion of one core, or do not corrode；It, can be first to straight in the direct exposure area template of binaryzation i.e. before carrying out step S113 Connect exposure area, corroded based on the first core or the second core, then based on after corrosion the direct exposure area template of binaryzation and Binaryzation breast image determines above-mentioned prime area.Wherein, the gray scale of the first gray areas is not less than the second predetermined gray scale threshold Value, and the second core is greater than the first core.

In an alternative embodiment, based on the direct exposure area template of binaryzation and binaryzation breast after corrosion After image determines above-mentioned prime area, it can be also first based on breast image and the first predetermined gray threshold, obtain the first low ash Template is spent, and highly attenuating organization chart picture can be presented as first area in the first low ash degree template；Then, to the first low ash degree mould The first area of plate is expanded with predetermined expansion core to become second area and form the second low gray level areas template, and second The area of the area ratio first area in region is bigger；Then, the portion other than prime area can be removed in the second low ash degree template Point to form third low ash degree template, correspondingly, second area becomes third region, and the area in the third region is less than the The area in two regions；Finally, extend the linear feature when the extended line of the linear feature wait mend line is fallen in third region, To obtain mending line linear feature, and the extended length of above-mentioned linear feature is less than or equal on the linear feature extending direction The length of the internal diameter in third region；I.e. after carrying out step S113, first in prime area, to quilt on binaryzation breast image After the linear feature that highly attenuating organization chart picture is partly covered carries out benefit line, it is further continued for subsequent step S114.

Step S114 meets predetermined length along the extending direction acquisition on the boundary of beam-defining clipper in binaryzation breast image The outside at threshold value and the center far from breast image does not include the linear feature of direct exposure area.

Specifically, the linear feature can be obtained based on Hough transformation；Or based on binaryzation breast image, pass through The maximal projection value in pixel row direction and in pixel column direction is obtained, respectively to obtain above-mentioned linear feature.

Step S115, in the range of prime area limits, according to the physical location of linear feature and beam-defining clipper in breast The relationship of the correspondence projected position of image, determines the boundary of breast image.

Specifically, the corresponding projected position phase of the linear feature that can be will acquire and the physical location of beam-defining clipper in breast image Than when linear feature is closer to the center of breast image, then linear feature to be determined as to the boundary of breast image；Otherwise, then The physical location of beam-defining clipper can be determined as to the boundary of breast image in the correspondence projected position of breast image.

Figure 11 is the flow diagram for determining the method on boundary in another embodiment in breast image.Such as Figure 11 institute Show, the method that present invention also provides a kind of to determine boundary in breast image, which can be for via including beam-defining clipper Mammary machine captured by medical image, this method can comprise the following steps that

Step S121 obtains the direct exposure area template of binaryzation of breast image and the breast image.

Step S122, the physical location of used beam-defining clipper when based on the above-mentioned breast image of shooting, to the breast figure of acquisition The direct exposure area template of binaryzation of picture and the breast image is cut.

Step S123 obtains binaryzation breast image based on the breast image after cutting.

Step S124, the direct exposure area template of binaryzation based on above-mentioned binaryzation breast image and after cutting, determines Prime area.

Step S125, in the range of prime area limits, along the first extending direction of beam-defining clipper, in binaryzation breast The linear feature for meeting the first predetermined length threshold value is obtained in image, and the linear feature is determined as above-mentioned beam-defining clipper along first Extending direction corresponding boundary on breast image.

Step S126, in the range of prime area limits, along the second extending direction of beam-defining clipper, in binaryzation breast The linear feature for meeting the second predetermined length threshold value is obtained in image, and the linear feature is determined as beam-defining clipper and is extended along second Direction corresponding boundary on breast image；Wherein, the first extending direction is parallel to the second extending direction, and the first predetermined length Threshold value and the second predetermined length threshold value are equal.

Step S127, in the range of prime area limits, along the third extending direction of beam-defining clipper, in the binaryzation The linear feature for meeting third predetermined length threshold value is obtained in breast image, and the linear feature is determined as beam-defining clipper along third Extending direction corresponding boundary on breast image；Wherein, first third extending direction is extended perpendicularly to.

It should be noted that can be carried out simultaneously between above-mentioned step S125, step S126 and step S127, can also press It successively carries out according to scheduled sequence, can specifically be set according to actual demand.

It in an alternative embodiment, can be based on the global offset parameter of detector, first before carrying out step S124 Gray scale linearization process is carried out to breast image (or breast image after cutting)；Simultaneously can also to the breast image after cutting into The processing of row bilinear interpolation downscaled images, and the contracting of arest neighbors interpolation is carried out to the direct exposure area template of binaryzation after cutting Small image procossing, and the ratio of the ratio of bilinear interpolation downscaled images and arest neighbors interpolation downscaled images is identical.In addition, subsequent After obtaining beam-defining clipper region template, the amplifying operation of corresponding ratio is carried out, so that final beam-defining clipper region template and sanction The size of breast image after cutting is adapted.

Figure 12 is a kind of flow diagram for the method that boundary is determined in breast image in one embodiment.Such as Figure 12 institute Show, in an alternative embodiment, the method that present invention also provides a kind of to determine boundary in breast image, the breast figure As that can be via medical image captured by the mammary machine including beam-defining clipper, this method be can comprise the following steps that

Step S131 obtains the direct exposure area template of binaryzation of breast image and the breast image；Wherein, breast figure As having initial boundary, and the direct exposure area template of binaryzation includes direct exposure area.

Step S132 is based on breast image, obtains binaryzation breast image.

Step S133 is based on binaryzation breast image and the direct exposure area template of binaryzation, determines prime area.

Step S134 obtains efficiency frontier linear feature of the breast image in all directions and just in prime area Boundary corresponding to initial line circle；If efficiency frontier linear feature has not been obtained in either direction, breast image is in the direction On initial boundary corresponding to boundary as efficiency frontier linear feature；Wherein, effective linear feature is referred to from binaryzation The outside that breast image obtains, length is greater than the center of predetermined length threshold value and separate breast image does not expose directly The linear feature in region.

Step S135 determines the boundary of breast image based on all efficiency frontier linear features.

Figure 13 is the flow diagram for obtaining the method on beam-defining clipper boundary in one embodiment in breast image.Such as Figure 13 It is shown, a method of obtaining the boundary of beam-defining clipper in breast image, it may include following steps:

Step S231, obtains breast image and binaryzation directly includes region template.

Specifically, obtaining breast image using medical imaging device, and direct based on breast image acquisition binaryzation Exposure area template；Wherein, above-mentioned breast image may include beam-defining clipper region and it is highly attenuating organize the formation of lower than predetermined The low gray level areas of gray threshold, and the direct exposure area template of binaryzation then may include direct exposure area etc..

Step S232 can carry out change of gradient and binary conversion treatment to breast image, to obtain binaryzation breast gradient map Picture.

Step S233 can be based on above-mentioned binaryzation breast gradient image, straight on each length direction to extract beam-defining clipper Line feature.

Step S234 is not present direct exposure area to the center outside far from breast image, is unsatisfactory for predetermined length threshold The linear feature that value and one end terminate at low gray level areas is extended.

Specifically, direct exposure area is not present in the outside when center of the above-mentioned linear feature far from breast image, and it is somebody's turn to do Linear feature is unsatisfactory for predetermined length threshold value again, and when one end of the linear feature terminates at low gray level areas again, can be straight to this Line feature is extended.It wherein, can be right on the extending direction of the linear feature when extending to above-mentioned linear feature The linear feature is extended in above-mentioned low gray level areas, and extended length may be less than or equal to above-mentioned low gray level areas and exist Interior electrical path length on the extending direction of above-mentioned linear feature.

Step S235, using the linear feature after the extension for meeting predetermined length threshold value as the boundary of above-mentioned beam-defining clipper.

In an alternative embodiment, above-mentioned step S234 can include:

Firstly, breast image and scheduled first gray threshold can be based on, to obtain the first low ash degree template；Wherein, on The highly attenuating organization chart picture stated is presented as first area in the first low ash degree template.

Secondly, can first area to the first low ash degree template expanded with predetermined expansion core, to obtain second area With the second low gray level areas template, i.e. second area is the first area after expansion, and the second low gray level areas template is that progress is swollen The first low ash degree template after swollen；Wherein, the area of above-mentioned second area is bigger than the area of above-mentioned first area.

Finally, when direct exposure area is not present in the outside at center of any linear feature far from breast image, and this is straight (i.e. benefit line) operation when the extended line of line feature is fallen in above-mentioned second area, can be extended to the linear feature.

In another alternative embodiment, above-mentioned step S234 can include:

Firstly, breast image and scheduled first gray threshold can be based on, to obtain the first low ash degree template；Wherein, on It states highly attenuating organization chart picture and is presented as first area in the first low ash degree template.

Secondly, can first area to the first low ash degree template expanded with predetermined expansion core to become second area simultaneously Form the second low gray level areas template；Wherein, the area of first area described in the area ratio of above-mentioned second area is big；

Then, prime area is limited in the second low ash degree template；Wherein, above-mentioned prime area can be according to binaryzation cream Room gradient image and the direct exposure area template of binaryzation are determined.

Later, above-mentioned by being removed in the second low ash degree template after the second low ash degree template is limited at prime area Part other than prime area, to form third low ash degree template；Correspondingly, above-mentioned second area is transformed into third region, And the area in the third region is less than the area of above-mentioned second area；

Finally, when direct exposure area is not present in the outside at center of the above-mentioned linear feature far from breast image, and should (i.e. benefit line) operation when the extended line of linear feature is fallen in above-mentioned third region, can be extended to the linear feature.

It should be noted that the method for determining boundary in breast image in above-mentioned implementation, obtaining in breast image Each embodiment such as the method in beam-defining clipper region, between each other may be used in the method on beam-defining clipper boundary, and detection breast image To be combined with each other, and identical or approximate technical characteristic and technology contents, it can be replaced mutually under the premise of not occurring contradictory And it combines.

In another alternative embodiment, present invention also provides a kind of computer readable storage mediums, and the calculating It can be stored with computer program on machine readable storage medium storing program for executing, while being that can be realized when the computer program is executed by processor In the application in any one embodiment the step of institute's procedures set forth.Wherein, above-mentioned computer readable storage medium may include Semiconductor memory, core memory, magnetic drum memory and/or magnetic disk storage etc..

Figure 14 is the structural schematic diagram of Medical Devices in one embodiment.As shown in figure 14, a kind of Medical Devices 5 may include The components such as radiographic source 51, detector 52, memory 53 and processor 54；Above-mentioned radiographic source 51 can be used for issuing X-ray, and should It may include the beam-defining clipper for limiting above-mentioned X-ray boundary in radiographic source 51；Detector 52, which then can be used for receiving penetrating to be positioned at, to be penetrated The X-ray of breast between line source 51 and detector 52, to generate breast image；Memory 53 can be used for storing computer journey The direct exposure area template of sequence, above-mentioned breast image, the corresponding binaryzation of the breast image and when shooting above-mentioned breast image Projected position etc. of the physical location of beam-defining clipper in breast image；Processor 54 then can be used for according to above-mentioned breast image, two The direct exposure area template of value and the projected position of beam-defining clipper etc. are realized when executing computer program as any in the application Each step of institute's procedures set forth in one embodiment.

Each technical characteristic of embodiment described above can be combined arbitrarily, for simplicity of description, not to above-mentioned reality It applies all possible combination of each technical characteristic in example to be all described, as long as however, the combination of these technical characteristics is not deposited In contradiction, all should be considered as described in this specification.

The embodiments described above only express several embodiments of the present invention, and the description thereof is more specific and detailed, but simultaneously It cannot therefore be construed as limiting the scope of the patent.It should be pointed out that coming for those of ordinary skill in the art It says, without departing from the inventive concept of the premise, various modifications and improvements can be made, these belong to protection of the invention Range.Therefore, the scope of protection of the patent of the invention shall be subject to the appended claims.

Claims

1. a kind of method in beam-defining clipper region in detection breast image, which is characterized in that the described method includes:

Obtain the breast image direct exposure area template of binaryzation corresponding with the breast image, the direct exposure region of binaryzation Domain template includes direct exposure area；

Obtain the linear feature on the binaryzation breast gradient image；And

Within the scope of the prime area on the breast image, linear feature is met into predetermined length threshold value, and be located at institute Region corresponding to the pixel column and pixel column except direct exposure area is stated as boundary corresponding to the beam-defining clipper.

2. the method according to claim 1, wherein straight on the acquisition binaryzation breast gradient image The step of line feature, comprising:

Based on the binaryzation breast gradient image, by being obtained in pixel row direction and in pixel column direction respectively Projection value, to obtain the linear feature on the binaryzation breast gradient image.

3. according to the method described in claim 2, it is characterized in that, the breast image includes breastwork side, far from the breastwork Opposite side, the upper side and lower side of side；

Wherein, when based on the binaryzation breast gradient image, by being obtained respectively in pixel row direction and in pixel Projection value on column direction, come when obtaining the linear feature on the binaryzation breast gradient image, in the pixel row direction On projection value include the first projection value corresponding with the upside and corresponding second projection value with the downside, in pixel column side Upward projection value includes third projection value corresponding with the opposite side, and the predetermined length threshold value includes the first predetermined length threshold Value, the second predetermined length threshold value and third predetermined length threshold value；And

Within the scope of the prime area on the breast image, linear feature is met into predetermined length threshold value, and position Region corresponding to pixel column and pixel column except the direct exposure area is as boundary corresponding to the beam-defining clipper The step of, comprising:

In within the scope of the prime area on the breast image, first projection value is greater than first predetermined length Threshold value and be located at the direct exposure area except pixel column, second projection value be greater than the second predetermined length threshold value And it is located at the pixel column except the direct exposure area and the third projection value greater than third predetermined length threshold value and position Pixel column except the direct exposure area is as boundary corresponding to the beam-defining clipper region.

4. according to the method described in claim 3, it is characterized in that, the binaryzation breast gradient image includes the first subgraph Picture, the second subgraph, third subgraph and the 4th subgraph, the extending direction one of the pixel column direction and the breastwork side It causes；

Wherein, described that gradient transformation and binaryzation are carried out to obtain the step of binaryzation breast gradient image to the breast image Suddenly, comprising:

Gradient transformation and binaryzation are carried out to the breast image, and first subgraph is obtained based on first gradient threshold value, And

Obtained respectively based on the second Grads threshold with the upside corresponding to second subgraph, with the downside corresponding to The third subgraph and with the 4th subgraph corresponding to the opposite side；

Wherein, the first gradient threshold value is greater than second Grads threshold, and first subgraph is for characterizing the breast Image outline feature.

5. according to the method described in claim 4, it is characterized in that, described based on the binaryzation breast gradient image and described The step of direct exposure area template of binaryzation determines prime area, comprising:

6. according to the method described in claim 5, it is characterized in that, described based on the binaryzation breast gradient image, By obtaining the projection value in pixel row direction and in pixel column direction respectively, to obtain the binaryzation breast gradient map As the step of upper linear feature, comprising:

Wherein, first projection value is on second subgraph in projection value maximum in the pixel row direction, institute Stating the second projection value is on the third subgraph in projection value maximum in the pixel row direction, the third projection value For on the 4th subgraph in projection value maximum in the pixel column direction.

7. according to the method described in claim 6, it is characterized in that, described based on the binaryzation breast gradient image and described The step of direct exposure area template of binaryzation determines prime area, further includes:

The prime area is determined based on first subgraph and the direct exposure area template of the binaryzation and this is initial The central point in region；

Wherein, central point and second Grads threshold based on the prime area, obtain respectively second subgraph and The third subgraph；It is special that second subgraph is used to characterize the straight line being located above the central point in the breast image Sign, the third subgraph are used to characterize the linear feature being located at below the central point in the breast image.

8. the method according to the description of claim 7 is characterized in that first projection value be on second subgraph in Maximum projection value in the pixel row direction, second projection value are on the third subgraph in the pixel column side Upward maximum projection value, the third projection value are on the 4th subgraph in throwing maximum in the pixel column direction Shadow value.

9. the method according to any one of claim 3-8, which is characterized in that the first predetermined length threshold value is the On the upside of the prime area of one predetermined ratio or the value of downside, the third predetermined length threshold value for the second predetermined ratio institute State the value of prime area breastwork opposite side；

10. the method according to claim 1, wherein the method also includes:

Obtain the physical location of the beam-defining clipper；

When can not determine the boundary of any side directions in beam-defining clipper region based on the binaryzation breast gradient image, then will The beam-defining clipper the side directions projected position of the physical location in breast image as boundary.

11. according to the method described in claim 10, it is characterized in that, the method also includes:

All boundaries based on acquisition cut the breast image, to remove the limit beam in the breast image Device region.

12. the method according to claim 1, wherein the method also includes:

According to projected position of the physical location of the beam-defining clipper on initial breast image, the initial breast image is carried out It cuts to generate the breast image after cutting；

13. the method according to claim 1, wherein the method also includes:

To the breast image carry out gradient transformation before, to the breast image carry out gray scale linearization process and/or in Value filtering processing.

14. the method according to claim 1, wherein the method also includes:

Before carrying out gradient transformation to the breast image, bilinear interpolation processing is carried out to reduce figure to the breast image Picture.

15. according to the method for claim 14, which is characterized in that the method also includes:

Wherein, the direct exposure area template reduced scale of the binaryzation is identical as the breast image reduced scale.

16. the method according to claim 1, wherein based on the binaryzation breast gradient image and described After the step of direct exposure area template of binaryzation determines prime area, the method also includes:

17. according to the method for claim 16, which is characterized in that described to partly being covered by highly attenuating organization chart picture The linear feature carries out the step of benefit line, comprising:

Based on the breast image and the first predetermined gray threshold, the first low ash degree template is obtained, wherein the highly attenuating tissue Image is presented as first area in the first low ash degree template；

The first area of the first low ash degree template is expanded with predetermined expansion core to become second area and form the Two low gray level areas templates, wherein the area of first area described in the area ratio of the second area is big；

The part other than the prime area is removed to form third low ash degree template, accordingly in the second low ash degree template Ground, the second area become third region, wherein the area in the third region is less than the area of the second area；With And

18. the method according to claim 1, wherein based on the binaryzation breast gradient image and described Before the step of direct exposure area template of binaryzation determines prime area, the method also includes:

19. according to the method for claim 18, which is characterized in that described in the direct exposure area template of the binaryzation Direct exposure area the step of being corroded, comprising:

Determine the direct exposure that the first gray areas in the high gray scale template accounts in the direct exposure area template of the binaryzation The gray scale of the ratio in light region, first gray areas is not less than the second predetermined gray threshold；

When the ratio is greater than predetermined ratio, the direct exposure area in the direct exposure area template of the binaryzation is carried out Corrosion based on second core；Otherwise, the corrosion based on first core is carried out, or is not corroded.

20. the method according to claim 1, wherein described be based on the binaryzation breast gradient image and institute State the step of direct exposure area template of binaryzation determines prime area, comprising:

The union of the binaryzation breast gradient image Yu the direct exposure area template of the binaryzation is sought, and with the union Boundary characteristic determines the prime area；Or

The first boundary characteristic in the binaryzation breast gradient image is obtained, the direct exposure area template of the binaryzation is obtained The second boundary feature, and the prime area is determined with the union of first boundary characteristic and the second boundary feature.

21. a kind of method for determining boundary in breast image, the breast image is shot via the mammary machine for including beam-defining clipper, It is characterized in that, which comprises

Obtain the physical location of beam-defining clipper when shooting；

Based on the breast image, binaryzation breast image and the direct exposure area template of binaryzation are obtained respectively；The two-value Changing direct exposure area template includes direct exposure area；

In the binaryzation breast image along the boundary of the beam-defining clipper extending direction acquisition meet predetermined length threshold value, And the outside at the center far from the breast image does not include the linear feature of the direct exposure area；And

In the range of the prime area limits, according to the physical location of the linear feature and the beam-defining clipper in the cream The relationship of the correspondence projected position of room image, determines the boundary of the breast image.

22. according to the method for claim 21, which is characterized in that the binaryzation breast image is to pass through: first to described Breast image carries out gradient transformation, carries out binaryzation to the transformed image of gradient again and obtain.

23. according to the method for claim 22, which is characterized in that the gradient transformation is the feature based on Sobel operator Extracting method.

24. according to the method for claim 22, which is characterized in that before gradient transformation further include:

25. according to the method for claim 24, which is characterized in that before the normalization further include:

Bilinear interpolation is carried out with downscaled images to the breast image, the normalization is carried out on image after diminution 's.

26. according to the method for claim 21, which is characterized in that described based on the binaryzation breast image and described Before the step of direct exposure area template of binaryzation determines prime area, the method also includes: it is direct to the binaryzation Corroded direct exposure area in the template of exposure area；

27. according to the method for claim 26, which is characterized in that described in the direct exposure area template of the binaryzation Direct exposure area the step of being corroded, comprising:

28. according to the method for claim 21, which is characterized in that described based on the binaryzation breast image and described After the step of direct exposure area template of binaryzation determines prime area, the method also includes:

It is described straight to partly being covered on the binaryzation breast image by highly attenuating organization chart picture in the prime area Line feature carries out benefit line.

29. according to the method for claim 28, which is characterized in that described to highly attenuating on the binaryzation breast image The step of linear feature that organization chart picture is partly covered carries out benefit line, comprising:

Based on the breast image and scheduled first predetermined gray threshold, the first low ash degree template is obtained, wherein the height declines Subtract organization chart picture and is presented as first area in the first low ash degree template；

30. according to the method for claim 21, which is characterized in that described in the range of the prime area limits, root Institute is determined in the relationship of the corresponding projected position of the breast image with the physical location of the beam-defining clipper according to the linear feature The step of stating the boundary of breast image, comprising:

Compared with correspondence projected position of the physical location of the beam-defining clipper in the breast image, when the linear feature more leans on When the center of the nearly breast image, the linear feature is determined as to the boundary of the breast image；Otherwise, by the limit beam The physical location of device is determined as the boundary of the breast image in the correspondence projected position of the breast image.

31. according to the method for claim 21, which is characterized in that it is described in the binaryzation breast image along described The extending direction acquisition on the boundary of beam-defining clipper meets predetermined length threshold value and the outside at the center far from the breast image is not wrapped In the step of including the linear feature of the direct exposure area:

The linear feature is obtained based on Hough transformation；Or

Based on the binaryzation breast image, by obtained in pixel row direction respectively and in pixel column direction most Big projection value, to obtain the linear feature.

32. a kind of method for determining boundary in breast image, the breast image is shot via the mammary machine for including beam-defining clipper, It is characterized in that, which comprises

Based on the physical location of beam-defining clipper when shooting the breast image, to the two of the breast image and the breast image The direct exposure area template of value is cut；

In the range of the prime area limits, along the first extending direction of the beam-defining clipper, in the binaryzation breast The linear feature for meeting the first predetermined length threshold value is obtained in image, and the linear feature is determined as the beam-defining clipper described in First extending direction corresponding boundary on the breast image.

33. according to the method for claim 32, which is characterized in that the beam-defining clipper is also flat with first extending direction The second capable extending direction extends, the method also includes:

In the range of the prime area limits, along second extending direction of the beam-defining clipper, in the binaryzation The linear feature for meeting the second predetermined length threshold value is obtained in breast image, and the linear feature is determined as the beam-defining clipper edge Second extending direction corresponding boundary on the breast image.

34. the method according to claim 32 or 33, which is characterized in that the beam-defining clipper also with the first extension side Extend to vertical third extending direction, the method also includes:

In the range of the prime area limits, along the third extending direction of the beam-defining clipper, in the binaryzation The linear feature for meeting third predetermined length threshold value is obtained in breast image, and the linear feature is determined as the beam-defining clipper edge The third extending direction corresponding boundary on the breast image；

35. according to the method for claim 32, which is characterized in that described based on the binaryzation breast image and described The direct exposure area template of binaryzation, before the step of determining prime area, the method also includes:

Accordingly with the breast image, to the direct exposure area template of the binaryzation after cutting carry out arest neighbors interpolation with Downscaled images.

36. according to the method for claim 32, which is characterized in that described based on the binaryzation breast image and described The direct exposure area template of binaryzation, before the step of determining prime area, the method also includes:

37. a kind of method for determining boundary in breast image, the breast image is shot via the mammary machine for including beam-defining clipper, It is characterized in that, which comprises

The direct exposure area template of binaryzation of breast image and the breast image is obtained, the breast image has initial edge Boundary；The direct exposure area template of binaryzation includes direct exposure area；

Based on breast image, binaryzation breast image is obtained；

In the prime area, efficiency frontier linear feature of the breast image in all directions and described initial is obtained Boundary corresponding to boundary；Wherein, effective linear feature refers to that obtaining from the binaryzation breast image, length is big There is no the linear feature of the direct exposure area in the outside at the center of predetermined length threshold value and separate breast image；

If the efficiency frontier linear feature has not been obtained in either direction, the breast image in this direction initial Boundary corresponding to boundary is as efficiency frontier linear feature；

38. a kind of method on the boundary for obtaining beam-defining clipper in breast image, which is characterized in that the described method includes:

The direct exposure area template of binaryzation of breast image and the breast image is obtained, the breast image includes beam-defining clipper area Domain and the highly attenuating low gray level areas lower than predetermined gray threshold organized the formation of, the direct exposure area template of binaryzation Including direct exposure area；

When the direct exposure area, the straight line is not present in the outside at center of the linear feature far from the breast image Feature is unsatisfactory for predetermined length threshold value again, and when one end of the linear feature terminates at the low gray level areas, to the straight line Feature is extended, and using the linear feature after the extension for meeting the predetermined length threshold value as the boundary of the beam-defining clipper.

39. according to the method for claim 38, which is characterized in that described when the linear feature is far from the breast image The outside at center be unsatisfactory for predetermined length threshold value again there is no the direct exposure area, the linear feature, and it is described directly When one end of line feature terminates at the low gray level areas, the step of extension to the linear feature, comprising:

Based on the breast image and scheduled first gray threshold, the first low ash degree template is obtained；Wherein, described highly attenuating group It knits image and is presented as first area in the first low ash degree template；

The first area of the first low ash degree template is expanded with predetermined expansion core to become second area and form the Two low gray level areas templates；Wherein, the second area is bigger than the first area；

When the direct exposure area is not present in the outside at center of the linear feature far from the breast image, and it is described straight When the extended line of line feature is fallen in the second area, extend the linear feature.

40. a kind of Medical Devices characterized by comprising

Detector receives the X-ray through the breast being positioned between the radiographic source and detector, to generate breast image；

Memory, for storing computer program, the breast image, the direct exposure area of the corresponding binaryzation of the breast image Template and when shooting the breast image beam-defining clipper physical location the breast image projected position；And

Processor, for the throwing according to the breast image, the direct exposure area template of the binaryzation and the beam-defining clipper Shadow position is realized when executing the computer program such as the step in any one of the claims 1-39 the method.