CN100594843C - Image displaying apparatus, image displaying method, storage medium and program - Google Patents
Image displaying apparatus, image displaying method, storage medium and program Download PDFInfo
- Publication number
- CN100594843C CN100594843C CN200510091028A CN200510091028A CN100594843C CN 100594843 C CN100594843 C CN 100594843C CN 200510091028 A CN200510091028 A CN 200510091028A CN 200510091028 A CN200510091028 A CN 200510091028A CN 100594843 C CN100594843 C CN 100594843C
- Authority
- CN
- China
- Prior art keywords
- mentioned
- breast
- image
- zone
- mammary gland
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Landscapes
- Apparatus For Radiation Diagnosis (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
To provide an image displaying apparatus capable of displaying an image having no overlapping mammary gland portion to make it easy to understand the mammary gland structure, an image displaying method, storage medium, and program.This image displaying apparatus is equipped with a storing unit 10 which stores data associated with an imaging area including breasts of a subject, a processing unit 13which generates a plurality of images respectively corresponding to a plurality of divided areas in the imaging area divided on the basis of positional information of pressure plates from the data, and a displaying unit 16 which displays the generated images.
Description
Technical field
The present invention relates to show graphics device, method for displaying image, storage medium and the program of the image relevant with the breast that between plate for forcing, is pressed into flat.
Background technology
As the conventional art of taking breast, there is the compressing breast to enlarge mammary gland, carry out the special-purpose X ray shooting of breast of 2 dimension image collection from 1 direction.At the three-dimensional structure that the mammary gland in the breast is made up of the cell as Fructus Vitis viniferae, each chamber was overlapping a lot of when a former state was taken, and can not carry out detailed evaluation.By clipping breast its attenuation is extended for this reason,, take with X ray still less in the eclipsed while that reduces mammary gland.
But, be difficult to differentiate and diagnose the mammary gland of pipe ball front when shooting and the eclipsed tumor of mammary gland of detector side.
[non-patent literature 1] Tao Wu and other 10 people work " Medical Physics; Vol.30No.3 " Tomographic mammography using a limited number oflow-dose conebeam projection images " " American Association ofPhysicists in Medicine distribution, in March, 2003, p.365-380
Summary of the invention
The object of the present invention is to provide a kind of demonstration to remove the overlapping of mammary gland and grasp image display device, method for displaying image, storage medium and the program of the image of mammary gland structure easily.
The 1st aspect of the present invention, image display device comprises: the parts of storing the data relevant with the camera watch region of the breast that comprises subject; Generate the parts of a plurality of images corresponding with a plurality of cut zone difference in the above-mentioned camera watch region of cutting apart according to the positional information of plate for forcing by above-mentioned data; The parts that show the image of above-mentioned generation.
The 2nd aspect of the present invention in the image display device of the image that generates breast, comprising: the parts of the regional area image of the above-mentioned image of a plurality of generations; When being presented on 1 picture, show the display unit of the mark of the subregion that is associated with the structure of breast to the subregion image of above-mentioned generation whole.
The 3rd aspect of the present invention in the image display device of the image that generates breast, comprising: the parts of calculating the Region Segmentation position; Generate the parts of the image of the cut zone of cutting apart with above-mentioned Region Segmentation position; The parts that show the image of above-mentioned generation, the parts of calculating the above-mentioned zone split position detect the positional information of the nipple of above-mentioned breast based on above-mentioned image, calculate the above-mentioned zone split position according to the positional information of the above-mentioned nipple that is detected.
The 4th aspect of the present invention in the image display device of the image that generates breast, comprising: the parts of calculating the Region Segmentation position; Generate the parts of the image of the cut zone of cutting apart with above-mentioned Region Segmentation position of calculating, the parts that show the image of above-mentioned generation, above-mentioned calculating apparatus extracts the general structure of breast inside based on above-mentioned image, calculates split position according to the position of this general structure.
The 5th aspect of the present invention, method for displaying image comprises: from the data relevant with the camera watch region of the breast that comprises subject, generate the step of a plurality of images corresponding with a plurality of cut zone difference in the above-mentioned camera watch region of cutting apart according to the positional information of plate for forcing; The step that shows the image of above-mentioned generation.
The 6th aspect of the present invention, be to record to be used to allow the machine of calculating realize the storage medium that the computer of the program of following component function can read, promptly, from with the data of the camera watch region of the breast that is included in subject, generate with the above-mentioned camera watch region of cutting apart according to the positional information of plate for forcing in a plurality of cut zone distinguish the parts of corresponding a plurality of images; The parts that show the image of above-mentioned generation.
The 7th aspect of the present invention, be to allow the program of the following component function of computer realization, promptly, from the data relevant, generate the parts of a plurality of images corresponding with a plurality of cut zone difference in the above-mentioned camera watch region of cutting apart according to the positional information of plate for forcing with the camera watch region of the breast that is included in subject; The parts that show the image of above-mentioned generation.
If employing the present invention then can provide to show image display device, method for displaying image, storage medium and the program of removing the overlapping of mammary gland and grasping the image of mammary gland structure easily.
Description of drawings
Fig. 1 is a structure chart of showing the X ray breast-image capturing device of embodiments of the present invention.
Fig. 2 is the figure that shows the display frame example on the display unit that is displayed on Fig. 1.
Fig. 3 is the oblique view of major part structure of X ray camera-shooting table of the X ray breast-image capturing device of exploded view 1.
Fig. 4 is the figure of the action of the X-ray tube ball when showing X ray shooting in the present embodiment.
Fig. 5 is the figure of the imageable areas (camera watch region) of the image forming process unit of exploded view 1.
Fig. 6 is used to supply the figure that the unitary part projected image of projection process of Fig. 1 generates.
Fig. 7 is the figure of the Region Segmentation face (median plane) of exploded view 6.
Fig. 8 is the figure that the Region Segmentation face that is used to supply Fig. 1 is set definite processing of the Region Segmentation face (median plane) that the Flame Image Process of processing unit produces.
Fig. 9 is the figure of the Region Segmentation face (predetermined surface) of exploded view 6.
Figure 10 is the figure of the Region Segmentation face (nipple face) of exploded view 6.
Figure 11 is used to supply the figure that is set definite processing of the Region Segmentation face (nipple face) that the Flame Image Process of processing unit carries out by the Region Segmentation face of Fig. 1.
Figure 12 is the figure that example is formed on the lip-deep ditch of plate for forcing in the present embodiment.
Figure 13 sets in the processing unit 12 at the Region Segmentation face of Fig. 1, and nipple extracts the figure that supplies that handles.
Figure 14 sets in the processing unit 12 at the Region Segmentation face of Fig. 1, and another nipple extracts the figure that supplies that handles.
Figure 15 sets in the processing unit 12 at the Region Segmentation face of Fig. 1, and divisional plane is set the figure that supplies that handles.
Figure 16 sets in the processing unit 12 at the Region Segmentation face of Fig. 1, and threshold value is determined the figure that supplies of processing.
Figure 17 sets in the processing unit 12 at the Region Segmentation face of Fig. 1, is used to supply the figure that the extraction in mammary gland zone is handled.
Figure 18 shows the relevant rectangular histogram of volume data of breast that is replaced into the state of fat deposit with the mammary gland layer.
Figure 19 is the figure that supplies that interior histogrammic curve approximation processing is handled in the extraction in the mammary gland zone of Figure 18.
Figure 20 is the relevant rectangular histogram of volume data of showing with the breast of inhomogeneous high concentration.
Figure 21 is the figure that supplies that interior histogrammic curve approximation processing is handled in the extraction in the mammary gland zone of Figure 20.
Figure 22 shows and has received the edge and emphasize the relevant histogrammic figure of volume data that handles.
Figure 23 is the relevant histogrammic figure of volume data that shows with the high state of breast density.
Figure 24 shows in the present embodiment, the figure of the summary of the device of usefulness ultrasonic scanning breast area.
Figure 25 sets in the processing unit 12 at the Region Segmentation face of Fig. 1, nipple and marker (marker) localization process supply figure.
Figure 26 sets in the processing unit 12 at the Region Segmentation face of Fig. 1, and the tracking of latex dust is handled supplies figure.
Figure 27 is in the present embodiment, shows the figure of corresponding relation of the dividing method of relative camera method etc.
The specific embodiment
Following with reference to the description of drawings embodiments of the present invention.If the employing present embodiment, then the present invention X ray stereo photographic device as the radiodiagnosis device can be provided, be loaded into image display device in the X ray stereo photographic device, method for displaying image, with the program of the combination of computer realization image display part, write down the storage medium that the computer of this program can read.
Adopt X ray stereo camera shooting method in the present embodiment.X ray stereo camera shooting method according to see through by a plurality of X ray of multi-direction shooting image reconstruct cross-section image (for the difference that reconstructs of X ray CT, hereinafter referred to as " imaging ").As the processing of reality, see through image according to a plurality of X ray of taking by a plurality of directions, the each point on specified section (pixel) is gone up many lonizing radiation of selecting by this point, selects the pixel value corresponding with these lonizing radiation to add.By specifying the position distinctness on the section to appear before one's eyes on this handling principle, leave the position of the position of section and blur.In such X ray stereo camera shooting method, in the X ray awl overlapping areas of group shot, can carry out the imaging of cross-section image (with reference to Fig. 4, Fig. 5) as direction.
Fig. 1 shows the structure of the X ray breast-image capturing device of embodiments of the present invention.X ray breast-image capturing device has X ray camera-shooting table 1.X ray camera-shooting table 1 as shown in Figure 3, equipment X-ray tube ball 2 on a side's of C type arm 17 end.The high-tension generation X ray that applies that X-ray tube ball 2 receives from high-voltage generator 3.On the opposing party's of C type arm 17 end, mounting plane detector (also being called flat-panel detector) 5 on the direction relative with X-ray tube ball 2.It is the Direct Transform type of electric signal that area detector 5 has the incident X-rays Direct Transform, perhaps incident X-rays is transformed to light with fluorophor, and this light is transformed to a plurality of quasiconductor detecting elements that the indirect variable of electric signal is remodeled.A plurality of quasiconductor detecting elements are aligned to 2 dimension clathrates.In a plurality of quasiconductor detecting elements, follow signal charge that X ray incident takes place to read as digital signal via data collection module 6.Area detector 5 can the alternative image booster and the combination of optical photographing unit.
The X ray breast-image capturing device of present embodiment also possesses when possessing X ray camera-shooting table 1: system control unit 7, shooting control unit 8, operating unit 9, image storage unit 10, image forming process unit 11, Region Segmentation face are set processing unit 12, projection process unit 13, section conversion process unit 14, picture formation unit 15, display unit 16.Shooting control unit 8 is carried out the multi-direction shooting (stereo camera shooting) (with reference to Fig. 4 (a), Fig. 4 (b)) of the breast of the subject that is pressed into flat on plate for forcing 18,19 as object by according to imaging conditions control high voltage generating unit 3, arm rotating mechanism 4, the data collection module 6 set in advance.The breast on right side and the breast in left side are carried out stereo camera shooting separately, the relevant a plurality of X ray of the right breast different with the shooting direction are seen through the data of image and see through the data storage of image on image storage unit 10 with the relevant a plurality of X ray of the different left breast of direction of making a video recording.Image forming process unit 11 as shown in Figure 6, for the equally spaced a plurality of sections in the imageable areas (abbreviating camera watch region as) of the roughly square shape of X beam and a pair of plate for forcing 18,19 regulations, a plurality of cross-section images that imaging is corresponding respectively.The set of being made up of a plurality of cross-section images is called volumetric image or only calls volume.The data of a plurality of cross-section images, promptly the data of volumetric image are stored in the image storage unit 10.
Region Segmentation face setting processing unit 12 is set at camera watch region a plurality of, is typically for thickness direction (Y direction), sets the Region Segmentation face that is used to be divided into 2 subregions.The Region Segmentation face can manually set arbitrarily in manipulation operations unit 9, and this initial position is automatically set by Region Segmentation face setting processing unit 12.Its initial Region Segmentation face is selected from 3 kinds.Selection operation is carrying out in advance or in the processing procedure.
The 1st Region Segmentation looks like shown in Figure 7, is the median plane of 2 plate for forcing 18,19 up and down that clips breast.Specifically, distance (thickness of oppressed breast) between the top of downside plate for forcing 18 and upside plate for forcing 19 following is set to T, plate for forcing 18,19 is parallel, sets the 1st Region Segmentation face at the face (median plane) to the distance that is positioned at T/2 above downside plate for forcing 18.Because downside plate for forcing 18 is fixed on the C type arm 17, so, then can set median plane if the lower position of upside plate for forcing 19 is determined.The lower position of upside plate for forcing 19 can detect by the position sensor of rotary encoder etc.And can determine the lower position of upside plate for forcing 19 by Flame Image Process.The lower position of upside plate for forcing 19 is because be the border of oppressed breast and upside plate for forcing 19, so can determine from the section that soft tissue looks like to change to the even picture of plate for forcing 19 as shown in Figure 8, by near a plurality of sections it, extracting cross-section image.
At this moment, as shown in figure 12,, then from the statistical picture data that generated, determine the surface location of plate for forcing 18,19 easily if on each surface of plate for forcing 18,19, the rill 21 of diving is very much formed clathrate.Because know the figure of the ditch 21 of setting, so can in 2 dimension images, calculate the index what degree this figure comprises by images match.For each layered image of image, investigate its index, the layering that index is the highest is the surface location of plate for forcing 19.If the ditch of different images is set, then can independently determine to go up the lower compression plate surface on last lower compression plate 18,19.Also the not marker of transmitted X-rays can be set without ditch 21 on plate for forcing 18,19.
The 2nd Region Segmentation looks like shown in Figure 9, is set at above the downside plate for forcing 18 of 2 plate for forcing 18,19 up and down of clipping breast beginning and only leaves on the face of position of distance D of regulation.As mentioned above, because select plate for forcing 18 to be fixed on the C type arm 17, so no matter the following position of upside plate for forcing 19 how, all is set on the certain position all the time.
The position of downside or upside plate for forcing 18,19 is determined as described below.Each layering hypothesis of volumetric image is and plate for forcing 18,19 almost parallels.The direction of plate for forcing is because determined by the condition of imaging processing, so generate layering as described above easily.
The layering number of plies of volumetric image is set to N, and initial stratified position is set to z1, and last layering position is set to z2, and the surface of the plate for forcing 18,19 of downside or upside is set to z3.Z1, z2 are the values of the condition indication during as imaging.Existing z3 is the situation of the constant in the device design, and as mentioned above, also having is the situation of using the measured value of sensor measurement.In addition, also having is the situation of the value determined according to image as mentioned above.
The Region Segmentation face for example begins to leave under the locational situation of distance L in the positive direction of z axle being designated as from the surface of plate for forcing 18,19, and the position of Region Segmentation face becomes the layering position of a=1+N (z3+L-z1)/(z2-z1) layer.General this value is not an integer, is illustrated between the integer of the integer of the maximum that is no more than a and the minimum bigger than a to have the Region Segmentation face.
As Figure 10, shown in Figure 11, that threshold process and figure matching treatment can be extracted accurately from utilizing according to cross-section image, by nipple and parallel with plate for forcing 18,19 or tilt a plurality of, select the highest section of integrated value of pixel value in each face, should be set at the 3rd Region Segmentation face by selecteed face.And then, replacing face by nipple, also can set the face of the center of gravity of passing through the pixel value in the camera watch region as the 3rd Region Segmentation face.
And then relevant nipple extracts the details of handling and narrates later on.
Camera watch region is divided into roughly first subregion (A) of upside plate for forcing 19 sides and roughly the latter half zone (B) of downside plate for forcing 18 sides by the Region Segmentation face.A plurality of cross-section images of projection process unit 13 in first subregion roughly,, generate and the corresponding projected image in first subregion roughly by maximum projection process (MIP).In addition, by the maximum projection process, generate and the regional corresponding projected image of the latter half roughly a plurality of cross-section images of projection process unit 13 in the latter half zone roughly.MIP handles as known in the art, sets the face parallel with plate for forcing 18,19 as the perspective plane at this, with this perspective plane vertical (lonizing radiation) on the maximum of pixel value as the value of the intersection point pixel on these lonizing radiation and perspective plane.That is, projection process is independently implemented in each subregion for the thickness direction of oppressed breast.As projection process, be not limited to MIP, also can be the processing of the pixel value on the integration lonizing radiation.Projection process unit 13 to each self-generating of left and right sides breast haply half subregion the part projected image and roughly in the subregional part projected image of Lower Half, for the projected image of each whole view field of self-generating of left and right sides breast.The data storage of the projected image in the part projected image of left and right sides breast roughly first subregion separately and subregional part projected image of left and right sides breast roughly Lower Half separately and left and right sides breast group shot picture zone separately in memory element 10.
Picture constitutes the display frame that unit 15 constitutes display unit 17.As shown in Figure 2, clip the Central Line of the picture that dots,, be configured to symmetry with the part projected image of roughly the latter half of relevant left breast the part projected image of roughly the first half of relevant right breast.On the downside of the part projected image of the first half roughly, the part projected image about roughly the latter half of the part projected image of roughly the latter half of right breast and relevant left breast is clipped the picture Central Line be configured to symmetry.On the outside of the part projected image of the first half roughly clipping the balanced configuration of picture Central Line about the region-wide projected image of right breast and the region-wide projected image of relevant left breast.And then, the region-wide projected image of relevant left and right sides breast, corresponding with the suitable operation of operating unit 9, can switch to the cross-section image relevant that generates by section conversion process unit 14 with the Region Segmentation face.On the outside of the part projected image of the latter half roughly, looking like to clip the roughly balanced configuration of picture Central Line about the skiagraph picture of right breast and the skiagraph of relevant left breast.In each skiagraph picture, corresponding with the Region Segmentation face respectively linear light mark overlaps each other.The suitable operation of linear light target display position and operating unit 9 is corresponding can to move up and down.When online cursor moves, follow this mobile update part projected image.
And then, in Figure 27, show mammillary structure, corresponding image capture method and display packing with it.In Figure 27, " outside ", " inboard " are the marks corresponding with the structure of breast.According to divisional plane volumetric region is divided into 2, generates the image in zone separately, and this moment, the method in the face side zone that generates divisional plane as method A, the method in the rear side zone of generation divisional plane as method B.Each mark corresponding with method A, method B can be determined with the such table of Figure 27.
As mentioned above, in the present embodiment, make the image of the volumetric image (multilamellar cross-section image) of having cut apart the compressing breast, the image at the display surface and the back side.As definite method of the Region Segmentation face of cutting apart, use the central section of plate for forcing as the method for Region Segmentation face and by Flame Image Process extraction fat deposit as the method for Region Segmentation face etc.Because a plurality of mammary gland are overlapping, so might be diagnosed as tumor mistakenly, when not having the eclipsed image of mammary gland by making the breast split image, can providing in the present invention, diagnostic accuracy improves in common X ray mammography method.
And then the 3rd above-mentioned Region Segmentation face is the face by nipple.The nipple that following declare area divisional plane is set processing unit 12 extracts processing.
As shown in figure 13,3 axles of volume data are set to X, Y, Z.Volume data hypothesis is at X, and Y has n, m, l[pixel respectively on the Z direction] size.At first, from volume data, remove the image disruption (speckle noise etc.) of distribution by Filtering Processing.Below, the volume data passing threshold after filtered is handled by 2 values.Threshold value is in order to distinguish air section and the zone based on breast and nipple in addition is set at suitable value.For the voxel (voxel) that has above the value of threshold value, the interior voxel in zone that promptly mainly comprises breast and nipple for example provides " 1 ", and for the voxel with the value below the threshold value, promptly the voxel in the air section for example provides " 0 ".
At this, suppose that vector f is the vector of the high direction of the probability that exists of expression nipple.Calculate this vector f and the voxel inner product separately that has given " 1 ".From a plurality of voxels that given " 1 ", select the voxel of inner product maximum.The centre coordinate value of the voxel of selecting is defined as teat placement.
Nipple extracting method as other has semi-automatic method.At first, the operator manually specifies teat placement on 2 dimension images of the Projection Display that generates according to volume data, MIP demonstration (maximum projection) etc.In this specified phases, just determined the position on 2 dimension pictures, and then needs are determined the position of depth direction (projecting direction).
The operator determines teat placement along projection straight line (being typically and the orthogonal direction of picture) the retrieval pixel value of having passed through specified point on 2 dimension images with the method for one of following (a) and (b).
(a) select voxel a plurality of voxels on projection straight line with maximum pixel value, the position of this voxel of selecting as teat placement.
(b) a plurality of voxels on projection straight line, extract the voxel that has above the pixel value of defined threshold.The middle position in the successive zone of voxel (voxel sequence) of extracting is arranged to teat placement.When the voxel sequence exists under a plurality of situations, select the longest voxel sequence.And then threshold value for example can be determined according to the display condition (grey level and window width) of MIP.
And then in other nipple extracting method, handling (section conversion process: multiplanar reconstruction processing) on the arbitrary cross-section of Sheng Chenging from volume data by MPR, the operator specifies teat placement, appointed positions as teat placement.
And then in other nipple extracting method, on the image that is shown by volume composition (VolumeRendering), the operator specifies teat placement.From this position, determine 3 dimension coordinates of nipple as described below.Along the projection straight line of the point that in volume composition image, passes through to click, the opacity that retrieval is used in the volume composition is made.If this opacity more than or equal to the 1st threshold value, then adds up opacity, this aggregate-value is reached the point of the 2nd threshold value as teat placement.
And then other nipple extracting method is to implement Flame Image Process for each section of MPR pictorial display, automatically calculates the image processing method of teat placement.In this image processing method, mainly comprise 2 kinds of image processing methods.The 1st kind of method is to spread the 2 dimension masks of differentiating the nipple shape, carries out the method for figure coupling (Pattern Matching Algorithm).The 2nd image processing method is as follows.At first, from volume data, remove the picture noise (speckle noise etc.) of distribution by Filtering Processing.Below, according to the volume data of pixel value 2 values through Filtering Processing.The threshold value that the is used for 2 values threshold value with the zone usefulness that is used for separation of air zone and in addition breast etc. as mentioned above is identical.For the Filtering Processing of being implemented a general subdifferential, 2 subdifferentials etc. by the volume data of 2 values, extract profile.The contour line that arrangement is extracted.
Multiple spot on the contour line is asked separately the inclination of tangent line as Figure 14 example, calculating inclination is zero position, promptly is typically use and calculates teat placement with the connecting point position of the tangent line B of the health axle almost parallel of subject.Perhaps with the position of differential value maximum as teat placement.
And then in other nipple extracting method, at first, the suitable high frequency of the volume data of nipple is removed filter process, remove the spatial radio-frequency component that is included in the volume data.It is air section and the zone that comprises breast and nipple in addition that passing threshold is handled volume data 2 values.If desired, then as processings of filling a vacancy, to comprise air section in addition breast and the zone (" 1 " zone) of nipple, after the expansion process (dilation) of small distance, carry out the shrink process (erotion) of same distance.
Only shrink process is carried out in " 1 " zone that has received the volume data A that handles of filling a vacancy, only carry out expansion process with distance n2 for this result with distance n1.At this moment, the n2 hypothesis is more than or equal to n1.The volume data that has received this processing is called B.This computing has the effect on the surface in smoothing " 1 " zone, its result, and volume data B is corresponding with the volume data A that has removed jut.By from volume data A, deducting the volume data B that has removed jut, generate the subregional volume data C of residual a plurality of jut with jut.One of a plurality of juts zone in the volume data C is a nipple area.
Labelling is carried out in connection zone in the volume data C (to be called: connectedcomponent labeling).Ask each regional volume and barycentric coodinates value.In advance serve as 1 zone of selection, basis to the definite alternative condition of every kind of image capture method (several method of standard of compressing breast).For example, select the zone of the inner product maximum of the vector f of barycentric coodinates value in the zone of volume below the above V2 of V1.Vector f is the vector of the high direction of the probability of expression nipple existence.Ask the center of gravity in the zone of selecting, it as teat placement.
Below, set in the processing unit 12 at divisional plane, also have the function of determining divisional plane according to the changes in distribution in mammary gland zone.With reference to Fig. 5, the median plane that separates 2 mammary gland zones as divisional plane.At first, the method for extracting mammary gland from volume data is described.The zone of mammary gland is the pixel value bigger than other many structures of breast on the x-ray tomography image.Thereby, if determine a certain threshold value, then can handle the zone of roughly trying to achieve mammary gland by 2 values.Threshold value determination method is narrated later on.In camera watch region, set a plurality of minutes aspects with plate for forcing 18,19 approximate horizontal.In other words, on oppressed breast thickness direction, branch aspect (XZ face) is set with certain interval.Ask the basal area in mammary gland zone of every extraction to each minute aspect.If basal area is asked in each layering position, then shown in Figure 15 (a), can obtain representing the curve chart of the variation of the mammary gland zone basal area relevant with thickness direction.On this curve chart, ask scope with zone also bigger (scope of ordering to B from the T point) than certain certain area (dotted line), the mid point M1 of this scope as the layering position of determining.
As additive method, shown in Figure 15 (b), ask the layering position P1 that occurs peak value from mid point M1 to the below basal area.Equally, from the peak value that M1 begins to ask the top, determine this layering position P2.Be set to P2.The mid point of P1 and P2 as M2.
As divisional plane, in the present example, also can determine T, P1, M1, P2, these a plurality of layerings of B, it be made the projection image that is segmented on a plurality of zones as divisional plane show to the layering position of trying to achieve like this (M1 or M2).
And then, try to achieve the display condition of the pictorial display of MIP that threshold value can be set from the operator and MPR etc.In this case, wish to constitute in the image of MIP and MPR etc., (with for example redness etc.) the overlapping threshold value of extracting in the fixed threshold value that is presented at, under the situation of display condition change,, determine threshold value once more according to this operator's operation, again extract the zone, upgrade superimposed images.In addition, the known several automatic method that has optimal threshold, big Tianjin method and percentage ratio to piece definite threshold value such as method together, these methods can be used.
Another that the following describes threshold value in the processing unit 12 determined method, and processing unit 12 is at first made the rectangular histogram (with reference to Figure 16) of pixel value from volume data.Judge zone beyond the breast according to rectangular histogram, this zone as α.Be set to β, γ according to the intramammary zone of rectangular histogram.Value in the black of this Figure 16 part is set to φ, and the value of white portion is set to ψ (φ＜ψ).
Region alpha: the fat deposit that is equivalent to breast.
Region beta: the main structure with breast is not corresponding.
Zone γ: be equivalent to the mammary gland zone of breast etc.
The pixel value in background (air) zone is also littler than φ, does not show in the rectangular histogram of Figure 16.
Set different values on region alpha and the region beta with on the regional γ, generating 2 value images, extracting the mammary gland structure.
Below, determine threshold value.Breast generally has 3 kinds of mammary gland structural regimes.State depends on age, suckling experience etc.
State a: the mammary gland layer all is replaced into the breast (with reference to Figure 17 (a)) of fat deposit
State b: the state of inhomogeneous high concentration breast (with reference to Figure 17 (b))
State c: the breast that breast density is high (with reference to Figure 17 (c))
At this, the expression breast has 3 states.Rectangular histogram when Figure 18 example mammary gland layer all is replaced into the breast of fat deposit.All be replaced in the breast of fat deposit at the mammary gland layer, mammary gland structure realm γ is divided into mammary gland zone and these 2 zones, non-mammary gland zone easily from rectangular histogram.
If pixel value is set to X, the pixel sum is set to Y, then can show with Figure 19 with curve of approximation f (X) performance rectangular histogram the time.
The boundary position φ of operator appointed area α and β, and the boundary position ψ of region beta and γ.The differential value of supposing Y=f (X) becomes positive X and is set to X1 (n is an integer) from negative.As the center, the ε % of voxel width η is set to region beta X1.And then, there are situation that ε sets by the operator and the rectangular histogram that uses a model to set the situation of default value.
Rectangular histogram when the breast of inhomogeneous high concentration is illustrated in Figure 20.If pixel value is set to X, the pixel sum is set to Y, then can represent with Figure 21 with curve of approximation f (X) performance rectangular histogram the time.But, do not see significant variation.
Thereby, if be image transformation the information of frequency field, be reduced on the frequency space under the low frequency via Fourier transform, then the edge is emphasized.In inverse-Fourier transform, return real space.In the rectangular histogram after frequency processing, as shown in figure 22, the edge is emphasized.
If pixel value is set to X, the pixel sum is set to Y, with the approximate rectangular histogram of curve of approximation f (X), then can obtain the rectangular histogram that the edge is emphasized.On the rectangular histogram of being emphasized at the edge, boundary position φ and the ψ of manual appointed area α of operator and β.
In the high breast of breast density, as Figure 23 example, nearly all be the mammary gland zone.The boundary position φ of operator appointed area α and γ.
For each breast, the threshold value of determining pixel value in advance with default value is also arranged, select the kind of breast, set the method for histogrammic threshold value.Which of the breast that the image that the operator selects to be prompted is the mammary gland layer breast that all is replaced into fat deposit, the breast of inhomogeneous high concentration, breast density is high.Pre-determine and each breast image corresponding threshold, use this threshold value to carry out threshold process.
In addition, have automatically determine the mammary gland layer all be replaced into the high breast of breast, the breast density of the breast of fat deposit, inhomogeneous high concentration which method and select after the establishing method of threshold value.The image that at first is prompted the mammary gland layer all is replaced into the high breast of breast, the breast density of the breast of fat deposit, inhomogeneous high concentration which get relevantly with modeled rectangular histogram in advance, selecting from correlation is the breast of which class.
Determine the threshold value of each breast image in advance.The threshold value of the kind of the selected relatively breast that goes out is carried out threshold process as default value (reference value).
Below, the method for setting divisional plane from volume data is described.This method from the volume data of the breast of being tried to achieve by ultrasonic scanning, is semi-automatically set divisional plane as shown in figure 24.At first, shooting person specifies by nipple or wants to cut apart the section (can be a plurality of) of demonstration.The collected volume data are come in the zone that comprises breast with ultrasound wave 3D scanning.Pictorial display is carried out as divisional plane in the plane of setting.
In addition, as shown below, also can set divisional plane according to volume data.Be that a kind of the use as the position relation of the position of knowing nipple and image can localized jig, collect the volume data in the zone that comprises breast, mechanically set the method for section.As shown in figure 25, a side of marker and nipple is alignd with the opposing party.
Collect the divisional plane of image setting by nipple.For example, the distance of service range plate for forcing 18,19 is set divisional plane.Nipple (marker) as center of rotation, is made the plane rotation, and the operator sets divisional plane.
The following describes the other method of from volume data, extracting the mammary gland structure.This method is that the operator selects to follow the tracks of starting point, uses general tracking to extract the method for mammary gland structure from this starting point.2 dimension pictorial display (projection, MIP, MPR) volume datas.As shown in figure 26, the operator clicks the position of latex dust.Use known blood vessel tracing method to follow the tracks of latex dust.
And then the present invention is not limited to above-mentioned embodiment itself, the implementation phase can the deformation construction key element in the scope that does not break away from its purport and specialize.In addition, by the suitable combination of disclosed a plurality of elements in the above-described embodiment, can form various inventions.For example, can from whole elements that present embodiment is showed, delete several elements.And then, also can suit to make up the element of different embodiments.
Claims (6)
1. image display device comprises:
Store the parts of the data relevant with the camera watch region of the breast that comprises subject;
Generate the parts of a plurality of images corresponding with a plurality of cut zone difference in the above-mentioned camera watch region of cutting apart according to the positional information of plate for forcing according to above-mentioned data;
The parts that show the image of above-mentioned generation;
The borderline divisional plane setup unit that is used for the divisional plane that above-mentioned camera watch region is divided into a plurality of cut zone is set in a plurality of mammary gland zone.
2. the described image display device of claim 1 is characterized in that: above-mentioned divisional plane setup unit is determined the border in above-mentioned mammary gland zone according to the variation of the cross-sectional area in the above-mentioned mammary gland zone relevant with the thickness direction of oppressed above-mentioned breast.
3. the described image display device of claim 2, above-mentioned divisional plane setup unit are handled according to the data relevant with above-mentioned camera watch region, passing threshold and are extracted above-mentioned mammary gland zone.
4. the described image display device of claim 3 is characterized in that: above-mentioned divisional plane setup unit is determined threshold value according to the histogrammic shape of the relation of expression pixel value relevant with the data of above-mentioned camera watch region and number of picture elements.
5. the described image display device of claim 3 is characterized in that: above-mentioned divisional plane setup unit is instructing specified value to be set at threshold value by the operator in the rectangular histogram of the relation of expression pixel value relevant with the data of above-mentioned camera watch region and number of picture elements.
6. the described image display device of claim 3 is characterized in that: above-mentioned divisional plane setup unit is implemented the edge to the data of above-mentioned camera watch region and is emphasized to handle.
Applications Claiming Priority (4)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2004226942 | 2004-08-03 | ||
| JP2004226942 | 2004-08-03 | ||
| JP2005216587 | 2005-07-27 | ||
| JP2005216587A JP4709604B2 (en) | 2004-08-03 | 2005-07-27 | Image display device, image display method, storage medium, and program |
Related Child Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN 200910160229 Division CN101612047B (en) | 2004-08-03 | 2005-08-03 | Image displaying apparatus |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN1739452A CN1739452A (en) | 2006-03-01 |
| CN100594843C true CN100594843C (en) | 2010-03-24 |
Family
ID=36092084
Family Applications (2)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN200510091028A Active CN100594843C (en) | 2004-08-03 | 2005-08-03 | Image displaying apparatus, image displaying method, storage medium and program |
| CN 200910160229 Active CN101612047B (en) | 2004-08-03 | 2005-08-03 | Image displaying apparatus |
Family Applications After (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN 200910160229 Active CN101612047B (en) | 2004-08-03 | 2005-08-03 | Image displaying apparatus |
Country Status (1)
| Country | Link |
|---|---|
| CN (2) | CN100594843C (en) |
Families Citing this family (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP5528518B2 (en) * | 2012-09-28 | 2014-06-25 | 富士フイルム株式会社 | Radiation image generating apparatus and method |
| WO2014157794A1 (en) * | 2013-03-29 | 2014-10-02 | 주식회사 레이언스 | Mammography device and subject shooting method using same |
| CN104127202B (en) * | 2013-05-03 | 2016-12-28 | 上海西门子医疗器械有限公司 | The localization method of compressor and positioner and gastrointestinal X-ray machine |
| CN104207797B (en) * | 2013-06-03 | 2018-06-19 | 上海联影医疗科技有限公司 | Determine that patient sizes, position and X ray pass through the method for patient's attenuation based on CT locating plates |
| TWI511085B (en) * | 2013-08-08 | 2015-12-01 | Univ Chien Hsin Sci & Tech | Nipple image segmentation method |
| CN104182965B (en) * | 2014-01-22 | 2016-03-30 | 上海联影医疗科技有限公司 | The method of chest muscle is split in a kind of galactophore image |
| WO2017054775A1 (en) * | 2015-09-30 | 2017-04-06 | Shanghai United Imaging Healthcare Co., Ltd. | System and method for determining a breast region in a medical image |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5050197A (en) * | 1988-05-26 | 1991-09-17 | Automed Oy | Mammographic methods and apparatuses |
| US5579360A (en) * | 1994-12-30 | 1996-11-26 | Philips Electronics North America Corporation | Mass detection by computer using digital mammograms of the same breast taken from different viewing directions |
| CN1145674A (en) * | 1994-12-30 | 1997-03-19 | 菲利浦电子有限公司 | Automatic segmentation, skinline and nipple detection in digital mammograms |
| US5640956A (en) * | 1995-06-07 | 1997-06-24 | Neovision Corporation | Methods and apparatus for correlating ultrasonic image data and radiographic image data |
-
2005
- 2005-08-03 CN CN200510091028A patent/CN100594843C/en active Active
- 2005-08-03 CN CN 200910160229 patent/CN101612047B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5050197A (en) * | 1988-05-26 | 1991-09-17 | Automed Oy | Mammographic methods and apparatuses |
| US5579360A (en) * | 1994-12-30 | 1996-11-26 | Philips Electronics North America Corporation | Mass detection by computer using digital mammograms of the same breast taken from different viewing directions |
| CN1145674A (en) * | 1994-12-30 | 1997-03-19 | 菲利浦电子有限公司 | Automatic segmentation, skinline and nipple detection in digital mammograms |
| US5640956A (en) * | 1995-06-07 | 1997-06-24 | Neovision Corporation | Methods and apparatus for correlating ultrasonic image data and radiographic image data |
Non-Patent Citations (2)
| Title |
|---|
| Tomographic mammography using a limited number oflow-dose cone-beam projection images. Tao Wu et al.Med. Phys.,Vol.30 No.3. 2003 |
| Tomographic mammography using a limited number oflow-dose cone-beam projection images. Tao Wu et al.Med.Phys,Vol.30 No.3. 2003 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN101612047A (en) | 2009-12-30 |
| CN1739452A (en) | 2006-03-01 |
| CN101612047B (en) | 2013-03-06 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP4709604B2 (en) | Image display device, image display method, storage medium, and program | |
| CN105705096B (en) | Method and apparatus for displaying medical images | |
| KR102109588B1 (en) | Methods for processing, displaying and navigating breast tissue images | |
| JP6392309B2 (en) | A system for navigating the tomosynthesis stack, including automatic focusing | |
| US8184890B2 (en) | Computer-aided diagnosis and visualization of tomosynthesis mammography data | |
| US8165365B2 (en) | Method and system for displaying tomosynthesis images | |
| US9401019B2 (en) | Imaging tomosynthesis system, in particular mammography system | |
| CN107405126B (en) | Retrieving corresponding structures of pairs of medical images | |
| US6487432B2 (en) | Method and system for selecting and displaying medical image data | |
| EP2212859B1 (en) | Method and apparatus for volume rendering of data sets | |
| CN102763137B (en) | For being generated the system and method for 2D image by chromatographic data collection | |
| CN102656608B (en) | Rebuild object of interest | |
| CN100594843C (en) | Image displaying apparatus, image displaying method, storage medium and program | |
| US20120157819A1 (en) | Imaging method and imaging device for displaying decompressed views of a tissue region | |
| CN110621233A (en) | System and method for synthesizing low-dimensional image data from high-dimensional image data using object mesh enhancement | |
| US8611492B2 (en) | Imaging method for rotating a tissue region | |
| CN101005803B (en) | Method for flexible 3dra-ct fusion | |
| CN101443815A (en) | Method and apparatus for reconstructing an image | |
| EP3560424A1 (en) | Multimodality 2d to 3d imaging navigation | |
| US10973485B1 (en) | Enhanced volume viewing | |
| KR20160117807A (en) | Method and apparatus for breast lesion diagnosis | |
| CN102548480B (en) | Device and method for displaying medical image and program | |
| Kita et al. | Correspondence between different view breast X-rays using a simulation of breast deformation | |
| US11158045B2 (en) | Method and apparatus for performing 3D imaging examinations of a structure under differing configurations and analyzing morphologic changes | |
| CN100583161C (en) | Method for depicting an object displayed in a volume data set |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C14 | Grant of patent or utility model | ||
| GR01 | Patent grant | ||
| C41 | Transfer of patent application or patent right or utility model | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20160720 Address after: Japan Tochigi Patentee after: Toshiba Medical System Co., Ltd. Address before: Tokyo, Japan, Japan Patentee before: Toshiba Corp Patentee before: Toshiba Medical System Co., Ltd. |