CN100589760C - X-ray CT apparatus and image processing apparatus - Google Patents
X-ray CT apparatus and image processing apparatus Download PDFInfo
- Publication number
- CN100589760C CN100589760C CN200580040425A CN200580040425A CN100589760C CN 100589760 C CN100589760 C CN 100589760C CN 200580040425 A CN200580040425 A CN 200580040425A CN 200580040425 A CN200580040425 A CN 200580040425A CN 100589760 C CN100589760 C CN 100589760C
- Authority
- CN
- China
- Prior art keywords
- mentioned
- unit
- image
- ray
- value
- 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.)
- Expired - Fee Related
Links
Images
Landscapes
- Apparatus For Radiation Diagnosis (AREA)
Abstract
An X-ray CT apparatus comprises a unit (20) for extracting a high-contrast area of a relatively high X-ray attenuation coefficient from an image, units (21, 22) for generating a blurred image on the high-contrast area on the basis of both the position of the high-contrast area extracted and a point-image intensity distribution function intrinsic to the apparatus, and a unit (23) for subtracting the blurred one from the image so as to generate a low-contrast image on a low-contrast area of a relatively low X-ray attenuation coefficient.
Description
Technical field
The present invention relates to generate the X ray CT device and the image processing apparatus of the image of subject inside according to subject is shone the data for projection that X ray obtains.
Background technology
The progress of X ray CT device is very remarkable, with hope higher meticulous (high-resolution) and in wide scope, photograph from the on-the-spot strong desire of medical treatment correspondingly, in recent years, developed many sections X ray CT device, it is considerably popularized.Should possess by many section X ray CT devices: irradiation tangential section direction (length direction of Wo Tai) expansion and have the x-ray source of the fan beam X ray of width; Having the 2 dimension detectors of structure of having arranged the detecting element row of multiple row (4 row, 8 row, 16 row etc.) at slice direction, is the scanning device that makes its action by many scannings or helical scanning.Thus, compare with single section X ray CT device, can high accuracy and in the short time, obtain volume data (volume data) in the wide scope of subject.
Have near at hand under the situation of high-contrast material (material with high X ray attenuation quotient) by the blood vessel in the CT image of such X ray CT device generation, be subjected to its fuzzy influence, the appearance of intravascular space is not distinct.Now, improve the method for the image quality of region-of-interest as the influence that is not subjected to high-contrast area, known following method: be set to carry out Flame Image Process after the suitable value consistent in the CT of high-contrast area value with on every side CT value, then, the value of recovery high-contrast area.But, in the method, can't remove high-contrast area fuzzy (opening flat 10-40372 communique) on every side well with reference to the spy.
Summary of the invention
The objective of the invention is to: absorb at the X ray that shows the subject inside obtain by X ray CT device under the situation of distributed image, be reduced in produce around the object of high X ray attenuation quotient fuzzy.
A first aspect of the present invention provides a kind of X ray CT device, be to collect the data for projection relevant with subject,, possess in the X ray CT device of the image of the inside of the above-mentioned subject of reconstruct according to the above-mentioned data for projection of collecting: constitute and from above-mentioned image, extract the unit that X ray declines into the higher high-contrast area of coefficient ratio out; Constitute according to the position and the device of the high-contrast area of above-mentioned extraction inherently, produce the unit of the broad image relevant with above-mentioned high-contrast area as intensity distribution function; Constitute in order to produce the low contrast regions relevant soft image lower, from the unit of the above-mentioned broad image of above-mentioned figure image subtraction with the X ray attenuation quotient.
A second aspect of the present invention provides a kind of X ray CT device, be to collect the data for projection relevant with subject, according to the above-mentioned data for projection of collecting and in the X ray CT device of the image of the inside of the above-mentioned subject of reconstruct, possess: constitute in order to produce the relevant soft image in zone lower, from the unit of the above-mentioned figure image subtraction broad image more relevant than higher zone with the X ray attenuation quotient being included in above-mentioned image with the X ray attenuation quotient; Accordingly above-mentioned soft image is categorized as the unit in a plurality of zones with the CT value; The pixel value of above-mentioned soft image is replaced into the unit of inherent value separately, above-mentioned a plurality of zones of having classified.
A third aspect of the present invention provides a kind of X ray CT device, be to collect the data for projection relevant with subject, according to the above-mentioned data for projection of collecting and in the X ray CT device of the image of the inside of the above-mentioned subject of reconstruct, possess: from above-mentioned image, according to a plurality of threshold values, sequentially extract a plurality of regional candidates' that vary in size that comprised unit out; Calculate the unit of a plurality of regional candidate center of gravity separately of above-mentioned extraction; According to the distance between the above-mentioned center of gravity, from a plurality of regional candidate of above-mentioned extraction, select the unit in specific zone.
Description of drawings
Fig. 1 is the structure chart of X ray CT device of the present invention.
Fig. 2 is the structure chart of the Flame Image Process parts of embodiments of the invention 1.
Fig. 3 is the key diagram of transfer function.
Fig. 4 is the structure chart of the Flame Image Process parts of embodiments of the invention 2.
Fig. 5 is the structure chart of the classification processing unit of embodiments of the invention 2.
Fig. 6 is 3 zones cutting apart are handled in expression by classification figure.
Fig. 7 is the flow chart that the zone of the zone decision parts of presentation graphs 5 determines treatment step.
Fig. 8 be Fig. 7 S13 supply figure.
Fig. 9 is the figure that is illustrated in region R 1 of extracting out among the S14 of Fig. 7 and the position of centre of gravity B1 that calculates in S15.
Figure 10 is the figure that is illustrated in region R 2 of extracting out among the S14 of Fig. 7 and the position of centre of gravity B2 that calculates in S15.
Figure 11 is the figure that is illustrated in region R 3 of extracting out among the S14 of Fig. 7 and the position of centre of gravity B3 that calculates in S15.
The specific embodiment
(embodiment 1)
Below, embodiments of the invention 1 are described.
Fig. 1 is the structure chart of the X ray CT device of embodiment 1.X ray CT device 1 possesses: constitute the pallet 2 of collecting the data for projection relevant with subject; The platform for sleeping in 3 that is used for mounting subject P and moves it; Be used for the input of operation of X ray CT device and the operating console of pictorial display.
Between the negative electrode anode of X-ray tube 5, apply tube voltage from high voltage generator 10.Supply with heater current from high voltage generator 10 to the filament of X-ray tube 5.The supply with heater current of applying by tube voltage produces X ray.For X-ray detector 6, adopt 1 dimension array-type detector or 2 dimension array-type detectors (being also referred to as many sections type detector).X-ray detection device for example has the foursquare sensitive surface of 0.5mm * 0.5mm.For example on the channel direction, arrange 916 x-ray detection devices.For example being set up in parallel 40 these row on slice direction, then is 2 dimension array-type detectors.By single row constitute be 1 the dimension array-type detector.
Transacter 11 is commonly referred to as DAS (data acquisition system).Transacter 11 will be transformed to voltage signal and amplify from the signal that detector 6 is exported each channel, and then will be transformed to digital signal.These data (initial data) are supplied with the operating console 4 of pallet outside.
12 pairs of data (initial data) from transacter 11 outputs of the pre-treatment parts of operating console 4 are implemented correcting process such as sensitivity correction, the output data for projection.This data for projection is sent to reconstruction processing parts 13.Reconstruction processing parts 13 for example according to by helical scanning, used the swept-volume of conical beam X ray or used the data for projection that they collect simultaneously, reconstructed image data, and storing in the image storage part spare 14.Flame Image Process parts 15 generate display image according to the view data that is stored in the image storage part spare 14.According to the input of operator, be used for being provided with etc. of condition setting, region-of-interest of display image to input equipment 16.To describe the image processing unit in the back in detail.Display device 17 is presented at the image that generates in the Flame Image Process parts 15.In addition, the scan control parts 18 of operating console 4 are according to operator's input, and control high pressure generating apparatus 10, pallet driving device 7, transacter 11, top board driving device 9 make scannings such as carrying out helical scanning.In addition, operating console 4 can constitute with specialized hardware, also can use a computer to realize same function with software.
Fig. 2 is the figure of structure of the Flame Image Process parts 15 of key diagram 1.In addition, Fig. 2 (A) is the ideograph as the blood vessel cross-section structure of an example of Flame Image Process object.Comprise blood vessel wall P1, support P2, contrast agent (blood flow part) P3, fatty P4.Below, the situation of having carried out handling with the image to such entity is that example describes.
Flame Image Process parts 15 possess region-of-interest image extraction parts 19, high-contrast extraction parts 20, convolution parts 21, CT value proportional zoom (scaling) parts 22, subtractor 23, adder 24 are set.
Region-of-interest is provided with image and extracts parts 19 out according to the input to input equipment 16, and region-of-interest is set, and extracts the view data and the output of this region-of-interest out from image storage part spare 14.(B) expression of Fig. 2 is provided with the image that image is extracted the region-of-interest of parts 19 extractions out by region-of-interest.
High-contrast is extracted parts 20 out from the image by the region-of-interest of extracting parts 19 extractions out, extracts the very high high-contrast area of X ray attenuation quotient out, the pixel group that promptly has very high CT value.As the high-contrast material, have imbed metallicity utensils such as tested intravital support, calcification for treatment calcium etc.In addition, the high-contrast material as in the soft image of having removed support etc. described later is not support etc., but is object with contrast agent mainly.Optionally use contrast agent with the CT value between support and the blood vessel wall.
Handle as concrete extraction, difference between the meansigma methods of the value by obtaining concerned pixel and the value of a plurality of pixels around this concerned pixel, obtain the differential map picture of the image of input, by using this differential map picture is implemented threshold process, obtain binary image with the edge corresponding threshold of support.This binary image is the different image of pixel value in the part that may have the high-contrast material and in addition part.That is, this binary picture similarly is the image (hereinafter referred to as high-contrast location drawing picture) that comprises the information of the position that may have the high-contrast material.(C) of Fig. 2 represents this high-contrast location drawing picture.In its differential, binary conversion treatment, can easily obtain high-contrast location drawing picture with fewer computing.
CT value proportional zoom parts 22 be standard CT value as contrast agent at this promptly according to as CT value value especially, and the image of exporting from convolution parts 21 (broad image) is carried out standardization (proportional zoom).This standardization is the processing that is used for the pixel value level of subtraction process described later.(D) expression of Fig. 2 is from the image (hereinafter referred to as high-contrast image) of CT value proportional zoom parts 22 outputs.
Adding unit 24 adds the positional information of high-contrast material by the high-contrast location drawing is looked like to add this soft image.Fig. 2 (F) represents this output image.This output image is that the positional information with the high-contrast material is added to from the image of original region-of-interest and has reduced image on the image of fuzzy (pseudomorphism) that the high-contrast material causes, is presented in the display device 17 via the processing of back level.In this output image, can observe well in original image because fuzzy influence and fat around the impalpable high-contrast material.
According to such embodiment 1, can be reduced in the fuzzy composition that produces of high-contrast material on every side, therefore can observe high-contrast material material on every side well.Have under support or the calcific situation in the blood capillary of the diameter of 5mm at the such 3mm that has of the coronary artery of heart especially, can observe the blood vessel wall around it, the savings state of fat, the state of contrast agent etc. well.
Also can carry out various distortion and implement this embodiment 1.For example, in the above-described embodiments, added high-contrast location drawing picture, but also can not addition, and only carry out the demonstration of soft image.In addition, the different sign of color of the position of all right overlapping expression high-contrast material.
And then, in the above-described embodiments, X ray CT device has been described, but also can the medical image-processing apparatus of the demonstration of carrying out image according to the data for projection of exporting from X ray CT device, CT image etc., have implemented.And then, in the above-described embodiments, 2 dimension treatment of picture have been described, but also can have carried out the processing of 3 d image by a plurality of 2 dimension images are carried out above-mentioned processing.
And then in the above-described embodiments, the pixel value that image has been described is the situation of CT value, but if expression X ray attenuation quotient then also can be the value beyond the CT value.In addition, in the above-described embodiments, obtained high-contrast location drawing picture by differential, binary conversion treatment, but also can use following method:, obtained the method for high-contrast location drawing picture the fuzzy transfer function that the contrary Using Convolution high-contrast material of image causes; The information of the position of storage support, shape is obtained the method for positional information etc. according to this information in advance.
(embodiment 2)
Below, embodiments of the invention 2 are described.In addition, omit explanation to the part identical with embodiment 1.
Embodiment 2 handles by classification (clustering), automatically classifies to make similar substance in the easy recognition image.Known under ideal condition, the CT value of platelet or fat expression-100 to 50, blood vessel wall is represented 50 to 129 CT value, is included in the CT value that contrast agent in the blood flow is represented 130-350.But, in the inspection of reality, being subjected to the influence of physique, beam hardening, the reconstruction of function of subject, the size of material, the state of the material beyond the region-of-interest, the concentration of contrast agent etc., CT has change.Therefore, in the method for carrying out zone extraction (Region Segmentation), the change of CT value that can't be corresponding such can't be classified well.In the present embodiment, from soft image, extract the main regional A corresponding, main area B, main with the contrast agent corresponding zone C corresponding accurately out with blood vessel wall with platelet.
Fig. 4 is the structure chart of the X ray CT device 1 of embodiment 2.In Fig. 4, between subtraction parts 23 and adding unit 24, has classification processing unit 25.This classification processing unit 25 is obtained a plurality of zones of regulation according to the CT value, and carries out the conversion process of CT value, makes to show each zone with same brightness or color.In addition, in the present embodiment, the situation that the CT value is divided into 3 zones has been described, but number that should the zone also can be other numbers.In addition, also can change number of regions accordingly with operator's input.
Fig. 5 is the structure chart of classification processing unit 25.Classification processing unit 25 possesses histogram generating means 26, zone decision parts 27, CT value transform parts 28.Histogram generating means 26 is obtained the rectangular histogram of the soft image of being imported.Rectangular histogram represents that the frequency that the pixel of each CT value occurs distributes.At this, the CT value is distributed in the scope between-100 to 537.
Fig. 7 represents to be extracted out from soft image by zone decision parts 27 step of the processing of 3 regional A, B, C.At first, from soft image, extract vascular wall area B (S11) out.From approximate centre a plurality of tracker wires are set radially by the specified blood vessel of operator.Follow the tracks of the CT value along each tracker wire.The position that the change of CT value is bigger, i.e. expression is defined as the interior point and the exterior point of blood vessel wall above 2 positions of the differential value of setting.Then, from soft image, extract platelet zone A (S12) out.As platelet zone A, extract the pixel group with CT value lower of the inboard that is positioned at blood vessel wall out than blood vessel wall.
Then, extract the zone C (S13～S19) of contrast agent out.At first, example is such as shown in Figure 8, and by the k-means method, the scope (in Fig. 6, the scope for-100 to 537) that the CT value on the rectangular histogram is distributed is divided into N section SEG
1～SEG
N(S13).If the section that the CT value is the highest is SEG
1N is set to surpass the integer that final classified regions is counted n (at this n=3).It is desirable to N is set to final classified regions and counts 2 times or 3 times of n.M is a treatment variable.M is initialized as 1.
At first, as shown in Figure 9, from soft image, extract out and section SEG
1Corresponding regional candidate R
1(S14).With section SEG
1Minimum as threshold value, from soft image, extract region R out
1The regional candidate R that calculating has been extracted out
1Position of centre of gravity B
1(S15).In S17, variable m is added 1 (S17).In addition, in S18, as the stop condition that is used to prevent spread, judge the minima of " CTmin seg m＜CTmax/2 " CTmin seg m section of being m, CTmax is the maximum CT value of soft image.That is, in the minima of section m 1/2 o'clock, stop area processing and amplifying less than the maximum CT value of soft image.During No in S18, carry out S19.During Yes in S18, turn back to S14.
As shown in Figure 10, from soft image, extract out and section SEG
1With section SEG
2Corresponding regional candidate R
2(S14).In fact, with section SEG
2Minimum as threshold value, from soft image, extract regional candidate R out
2The regional candidate R that calculating has been extracted out
2Position of centre of gravity B
2(S15).Calculate this zone candidate R
2Position of centre of gravity B
2With a regional candidate R before
1Position of centre of gravity B
1Between distance D
2, and with the regulation threshold value Th compare (S16).
Addendum modification D up to position of centre of gravity
mSurpass till the threshold value Th, the processing of S14～S18 is carried out in circulation.That is, along with increasing progressively of variable m, regional candidate R enlarges bit by bit.Along with the expansion of regional candidate R, center of gravity B displacement.When the extraction object was identical, the addendum modification of center of gravity B was fewer.Enlarging when make extracting object out and comprising other materials, the addendum modification of center of gravity B is bigger.That is, extract out to as if during contrast agent, even the enlarged area candidate R, position of centre of gravity does not roughly change yet.But when having added blood vessel wall etc. in extracting object out, regional candidate R becomes greatly by leaps and bounds, and also great changes have taken place for position of centre of gravity.The search position of centre of gravity section that great changes have taken place.
When m=3, as shown in figure 11, at regional candidate R
3Position of centre of gravity B
3With a regional candidate R before
2Position of centre of gravity B
2Between distance D
3When having surpassed threshold value Th (S16), with a regional candidate R before
2Be defined as the zone C (S19) of contrast agent.That is, the addendum modification of center of gravity is converged to regional candidate R less than the maximum of threshold value Th
2Be chosen as specific zone C.
Like this, by being conceived to the displacement that the zone enlarges the position of centre of gravity that causes, can be subjected to the influence of CT value variable hardly, the zone C of extracting contrast agent accurately out.In addition, for the extraction of blood vessel wall too, change, can be subjected to the influence of CT value variable hardly, extract the area B of blood vessel wall accurately out by being conceived to the CT value.And then, for hematoblastic extraction too,, can be subjected to the influence of CT value variable hardly by the inboard such condition of blood vessel wall is set, extract hematoblastic regional A accurately out.
In addition, zone decision parts 27 also can use the k-mean method, and this rectangular histogram is categorized as 3 zones.By the k-mean method, obtain and be used to obtain each regional threshold value Th, make that the dispersion value of the CT value in each zone is equal.By determining each zone to make that the dispersion value of the CT value in each zone is equal, thereby the high material of similarity is included in the same area.
Fig. 6 represents to carry out sectional example by present embodiment, and the CT value with-100 to 537 is categorized as 3 zones.The pixel of CT value of scope with regional A is main corresponding with fatty part.The pixel of CT value of scope with area B is main corresponding with the part of blood vessel wall.The pixel of CT value of scope with zone C is main corresponding with hematoblastic part.
CT value transform parts 28 will determine the CT value of the pixel that 3 zones of parts 27 extractions or classification comprise respectively to be replaced into each the regional eigenvalue that sets in advance by the zone.Thus, on display image, show the same area part, show can grasp each regional form with same brightness or color.
According to such present embodiment and since with the CT value show subregion accordingly image, the image of different materials such as the tissue that can grasp subject inside well, fat, contrast agent can be provided.In addition, automatically obtain the threshold value that is used to determine the zone, even therefore have under the situation of change, also can provide the image of having classified well in the CT value according to rectangular histogram.
In addition, in this embodiment, reducing behind the fuzzy composition (pseudomorphism composition) that produces around the high-contrast material, as object, therefore the processing of classifying can observe high-contrast material material on every side well with soft image.Have under support or the calcific situation in the blood capillary of the diameter of 5mm at the such 3mm that has of the coronary artery of heart especially, can observe the blood vessel wall around it, the savings state of fat, the state of contrast agent etc. well.
In addition, owing to the processing of classifying at the image that does not comprise the high-contrast material, add the information of high-contrast material then, so can carry out subregion to part well with low CT value.
In addition, the present invention has more than and is limited to the foregoing description, the implementation phase, in the scope that does not break away from its aim, can be out of shape and specialize element.In addition, a plurality of elements that disclosed by suitably making up the foregoing description can form various inventions.For example, can from the whole elements shown in the embodiment, delete several elements.And then, also can suitably make up the element among the different embodiment.
For example, in the above-described embodiments, use the k-mean method to obtain threshold value, but also can use group average method, other sorting techniques such as Ward method obtain threshold value.In addition, in the above-described embodiments, at having removed processings of classifying of high-contrast material and fuzzy image on every side thereof, but also can not such processing, and at the original image of region-of-interest, the processing of directly classifying, generation display image.
And then, in the above-described embodiments, the high-contrast location drawing is looked like to carry out addition, but also can not carry out addition, and only show soft image.In addition, the different labelling of color of also can overlapping demonstration representing the position of high-contrast material.
And then, in the above-described embodiments, X ray CT device has been described, but also can according to the medical image-processing apparatus that carries out pictorial display from the data for projection CT image of X ray CT device output etc., have implemented.And then, in the above-described embodiments, 2 dimension treatment of picture have been described, but also can have carried out the processing of 3 d image by a plurality of 2 dimension images are carried out above-mentioned processing.And then in the above-described embodiments, the pixel value that image has been described is the situation of CT value, but if the value of expression X ray attenuation quotient then also can be the value beyond the CT value.
According to the present invention, absorb at the X ray that shows the subject inside obtain by X ray CT device under the situation of distributed image, can be reduced in produce around the object of high X ray attenuation quotient fuzzy.
Claims (17)
1. X ray CT device is collected the data for projection relevant with subject in above-mentioned X ray CT device, according to the above-mentioned data for projection of collecting and the image of the inside of the above-mentioned subject of reconstruct, this X ray CT device is characterised in that and comprises:
Constitute and from above-mentioned image, extract the unit of X ray attenuation quotient out than higher high-contrast area;
Constitute according to the position and the device of the high-contrast area of above-mentioned extraction inherently, produce the unit of the broad image relevant with above-mentioned high-contrast area as intensity distribution function;
Constitute in order to produce the low contrast regions relevant soft image lower with the X ray attenuation quotient, from the unit of the above-mentioned broad image of above-mentioned figure image subtraction,
Accordingly above-mentioned soft image is categorized as the unit in a plurality of zones with the CT value; And
The pixel value of above-mentioned soft image is replaced into the unit of inherent value separately, above-mentioned a plurality of zones of having classified,
Above-mentioned unit of classifying possesses:
According to a plurality of threshold values, from above-mentioned soft image, extract the unit in a plurality of zones that vary in size out;
Calculate the unit of a plurality of zones center of gravity separately of above-mentioned extraction;
According to the distance between the above-mentioned center of gravity, from a plurality of zones of above-mentioned extraction, select the unit in the specific zone corresponding with extracting object out,
The above-mentioned unit of selecting selects the addendum modification of above-mentioned center of gravity to converge to zone less than the maximum of ormal weight as above-mentioned specific zone.
2. X ray CT device according to claim 1 is characterized in that:
The unit of above-mentioned extraction high-contrast area produces the differential map picture according to above-mentioned image, by above-mentioned differential map is looked like to carry out threshold process, and extracts above-mentioned high-contrast area out.
3. X ray CT device according to claim 1 is characterized in that:
The unit of above-mentioned generation broad image possesses:
To the location drawing picture of the high-contrast area of above-mentioned extraction, the above-mentioned unit of convolution as intensity distribution function;
According to the standard CT value of above-mentioned high-contrast area, the location drawing of the high-contrast area of having carried out convolution is looked like to carry out standardized unit.
4. X ray CT device according to claim 1 is characterized in that:
Above-mentioned high-contrast area be embedded to tested intravital metallic utensil or calcification for treatment calcium.
5. X ray CT device according to claim 1 is characterized in that also comprising:
With the location drawing picture of the high-contrast area of above-mentioned extraction and the unit of above-mentioned soft image addition.
6. X ray CT device according to claim 1 is characterized in that also comprising:
With the location drawing picture of the high-contrast area of above-mentioned extraction and the unit of the soft image addition after the above-mentioned displacement.
7. X ray CT device according to claim 1 is characterized in that:
Above-mentioned a plurality of zones of having classified comprise mainly be the zone of fat, mainly be the zone of blood vessel wall, mainly be the zone of contrast agent.
8. X ray CT device according to claim 1 is characterized in that:
The frequency that occurs according to the pixel of each CT value relevant with above-mentioned soft image distributes, and determines above-mentioned a plurality of threshold value.
9. X ray CT device is collected the data for projection relevant with subject in above-mentioned X ray CT device, according to the above-mentioned data for projection of collecting and the image of the inside of the above-mentioned subject of reconstruct, this X ray CT device is characterised in that and comprises:
Constitute in order to produce the zone relevant soft image lower, from the unit of the above-mentioned figure image subtraction broad image more relevant than higher zone with the X ray attenuation quotient being included in above-mentioned image with the X ray attenuation quotient;
Accordingly above-mentioned soft image is categorized as the unit in a plurality of zones with the CT value;
The pixel value of above-mentioned soft image is replaced into the unit of inherent value separately, above-mentioned a plurality of zones of having classified,
Above-mentioned unit of classifying is categorized as contrast medium areas, vascular wall area and platelet zone with above-mentioned soft image,
Above-mentioned unit of classifying possesses in order to determine above-mentioned contrast medium areas:
According to a plurality of threshold values, order is extracted the inner a plurality of regional candidates' that vary in size that comprise unit out from above-mentioned soft image;
Calculate the unit of a plurality of regional candidate center of gravity separately of above-mentioned extraction;
According to the distance between the above-mentioned center of gravity, from a plurality of regional candidate of above-mentioned extraction, select the unit of above-mentioned contrast medium areas,
The above-mentioned unit of selecting selects the addendum modification of above-mentioned center of gravity to converge to zone less than the maximum of ormal weight as above-mentioned contrast medium areas.
10. X ray CT device according to claim 9 is characterized in that:
The frequency that occurs according to the pixel of each CT value relevant with above-mentioned soft image distributes, and determines above-mentioned a plurality of threshold value.
11. X ray CT device according to claim 9 is characterized in that:
Above-mentioned unit of classifying is followed the tracks of the CT value in order to determine above-mentioned vascular wall area radially from the blood vessel approximate centre of above-mentioned soft image, the variation of above-mentioned CT value is defined as the inner and outer wall of blood vessel above 2 positions of setting.
12. X ray CT device according to claim 11 is characterized in that:
Above-mentioned unit of classifying is extracted pixel group above-mentioned blood vessel wall inboard, that have the CT value lower than above-mentioned blood vessel wall out in order to determine above-mentioned platelet zone.
13. an X ray CT device is collected the data for projection relevant with subject in above-mentioned X ray CT device, according to the above-mentioned data for projection of collecting and the image of the inside of the above-mentioned subject of reconstruct, this X ray CT device is characterised in that and comprises:
From above-mentioned image,, sequentially extract the inner a plurality of regional candidates' that vary in size that comprise unit out according to a plurality of threshold values;
Calculate the unit of a plurality of regional candidate center of gravity separately of above-mentioned extraction;
According to the distance between the above-mentioned center of gravity, from a plurality of regional candidate of above-mentioned extraction, select the unit in specific zone,
The above-mentioned unit of selecting selects the addendum modification of above-mentioned center of gravity to converge to zone less than the maximum of ormal weight as above-mentioned specific zone.
14. X ray CT device according to claim 13 is characterized in that:
The frequency that occurs according to the pixel with above-mentioned image-related each CT value distributes, and determines above-mentioned a plurality of threshold value.
15. an image processing apparatus is characterized in that comprising:
Constitute and from radioscopic image, extract the unit of X ray attenuation quotient out than higher high-contrast area;
Constitute according to the position and the device of the high-contrast area of above-mentioned extraction inherently, produce the unit of the broad image relevant with above-mentioned high-contrast area as intensity distribution function;
Constitute in order to produce the low contrast regions relevant soft image lower, deduct the unit of above-mentioned broad image from above-mentioned radioscopic image with the X ray attenuation quotient,
Accordingly above-mentioned soft image is categorized as the unit in a plurality of zones with the CT value; And
The pixel value of above-mentioned soft image is replaced into the unit of inherent value separately, above-mentioned a plurality of zones of having classified,
Above-mentioned unit of classifying possesses:
According to a plurality of threshold values, from above-mentioned soft image, extract the unit in a plurality of zones that vary in size out;
Calculate the unit of a plurality of zones center of gravity separately of above-mentioned extraction;
According to the distance between the above-mentioned center of gravity, from a plurality of zones of above-mentioned extraction, select the unit in the specific zone corresponding with extracting object out,
The above-mentioned unit of selecting selects the addendum modification of above-mentioned center of gravity to converge to zone less than the maximum of ormal weight as above-mentioned specific zone.
16. an image processing apparatus is characterized in that comprising:
Constitute from radioscopic image deduct be included in above-mentioned radioscopic image the X ray attenuation quotient broad image more relevant than higher zone, produce the unit of the soft image relevant with the lower zone of X ray attenuation quotient;
Accordingly above-mentioned soft image is categorized as the unit in a plurality of zones with pixel value;
The pixel value of above-mentioned soft image is replaced into the unit of inherent value separately, above-mentioned a plurality of zones of having classified,
Above-mentioned unit of classifying possesses:
According to a plurality of threshold values, from above-mentioned soft image, extract the unit in a plurality of zones that vary in size out;
Calculate the unit of a plurality of zones center of gravity separately of above-mentioned extraction;
According to the distance between the above-mentioned center of gravity, from a plurality of zones of above-mentioned extraction, select the unit in the specific zone corresponding with extracting object out,
The above-mentioned unit of selecting selects the addendum modification of above-mentioned center of gravity to converge to zone less than the maximum of ormal weight as above-mentioned specific zone.
17. an image processing apparatus is characterized in that comprising:
According to a plurality of threshold values, order is extracted the inner a plurality of regional candidates' that vary in size that comprise unit out from radioscopic image;
Calculate the unit of a plurality of regional candidate center of gravity separately of above-mentioned extraction;
According to the distance between the above-mentioned center of gravity, from a plurality of regional candidate of above-mentioned extraction, select the unit in specific zone,
The above-mentioned unit of selecting selects the addendum modification of above-mentioned center of gravity to converge to zone less than the maximum of ormal weight as above-mentioned specific zone.
Applications Claiming Priority (3)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2004342272 | 2004-11-26 | ||
| JP342272/2004 | 2004-11-26 | ||
| JP342273/2004 | 2004-11-26 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN101065065A CN101065065A (en) | 2007-10-31 |
| CN100589760C true CN100589760C (en) | 2010-02-17 |
Family
ID=38965623
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN200580040425A Expired - Fee Related CN100589760C (en) | 2004-11-26 | 2005-11-24 | X-ray CT apparatus and image processing apparatus |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN100589760C (en) |
Families Citing this family (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP5305821B2 (en) * | 2008-10-10 | 2013-10-02 | 株式会社東芝 | Medical image processing apparatus and medical image diagnostic apparatus |
| JP2012200285A (en) * | 2011-03-23 | 2012-10-22 | Toshiba Corp | Image processing apparatus, x-ray ct apparatus and image processing method |
| US9486178B2 (en) * | 2012-08-31 | 2016-11-08 | Shimadzu Corporation | Radiation tomographic image generating apparatus, and radiation tomographic image generating method |
| JP5971411B2 (en) * | 2013-05-10 | 2016-08-17 | 株式会社島津製作所 | Image processing device |
| CN104517263B (en) * | 2013-09-30 | 2019-06-14 | Ge医疗系统环球技术有限公司 | The method and apparatus for reducing pseudomorphism in computed tomography images reconstruct |
| CN107623816A (en) * | 2013-10-09 | 2018-01-23 | 佳能株式会社 | Image processing equipment, image pick up equipment and image processing method |
| CN107106109B (en) * | 2014-11-06 | 2020-11-06 | 皇家飞利浦有限公司 | Computed tomography system |
| EP3633603A4 (en) | 2017-06-30 | 2020-06-17 | Shanghai United Imaging Healthcare Co., Ltd. | Tissue density analysis method and system |
| CN107392893A (en) * | 2017-06-30 | 2017-11-24 | 上海联影医疗科技有限公司 | Tissue density's analysis method and system |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6266388B1 (en) * | 1999-09-07 | 2001-07-24 | General Electric Company | Methods and apparatus for two-pass cone beam image reconstruction |
| CN1547161A (en) * | 2003-12-08 | 2004-11-17 | 西安理工大学 | Automatic generating method for colour multi-window CT image |
-
2005
- 2005-11-24 CN CN200580040425A patent/CN100589760C/en not_active Expired - Fee Related
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6266388B1 (en) * | 1999-09-07 | 2001-07-24 | General Electric Company | Methods and apparatus for two-pass cone beam image reconstruction |
| CN1547161A (en) * | 2003-12-08 | 2004-11-17 | 西安理工大学 | Automatic generating method for colour multi-window CT image |
Also Published As
| Publication number | Publication date |
|---|---|
| CN101065065A (en) | 2007-10-31 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN100589760C (en) | X-ray CT apparatus and image processing apparatus | |
| JP4874810B2 (en) | X-ray CT system | |
| US7835562B2 (en) | Methods and apparatus for noise reduction filtering of images | |
| JP4751495B2 (en) | Method and apparatus for identifying structural image components | |
| JP4149598B2 (en) | Method for automatically setting collimator of X-ray imaging system during image acquisition, and X-ray imaging system | |
| US6754376B1 (en) | Method for automatic segmentation of medical images | |
| EP2891455B1 (en) | Radiation tomographic image-generating apparatus and radiation tomographic image-generating method | |
| US9659390B2 (en) | Tomosynthesis reconstruction with rib suppression | |
| CN100350434C (en) | Image data area extracting system and image data area extracting method | |
| JP5559375B2 (en) | X-ray computed tomography apparatus, reconstruction processing apparatus, and image processing apparatus | |
| US8483470B2 (en) | Radiological image area extracting apparatus, radiological image area extraction program, and radiographic apparatus | |
| KR101493375B1 (en) | Image processing apparatus, image processing method, and computer-readable storage medium | |
| JP2000011170A (en) | Method for emphasizing discrete pixel image and device therefor | |
| CN105321155A (en) | Ring artifact elimination method for CBCT image | |
| WO2007047059A2 (en) | Image processing method for boundary extraction between at least two tissues, the boundary being found as the cost minimum path connecting a start point and an end point using the fast marching algorithm with the cost function decreases the more likely a pixel is not on the boudary | |
| JP4236405B2 (en) | Tomographic segmentation | |
| JP2002099896A (en) | Abnormal shadow candidate detecting method and apparatus | |
| JP5329103B2 (en) | Image processing apparatus and X-ray CT apparatus | |
| US7751626B2 (en) | Method and apparatus for detection using gradient-weighted and/or distance-weighted graph cuts | |
| WO1985004090A1 (en) | X-ray ct image processor | |
| JP3229200B2 (en) | Bone measurement method | |
| JP2013119021A (en) | X-ray ct device and image processing method | |
| WO2009004571A1 (en) | Method and apparatus for image reconstruction | |
| JPH08287230A (en) | Computer-aided image diagnostic device | |
| JP2017131454A (en) | Medical image display system, medical image display program and medical image display method |
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: 20160706 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. |
|
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20100217 Termination date: 20181124 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |