CN102448374B - Image processing apparatus, X-ray image capturing apparatus and image processing method - Google Patents
Image processing apparatus, X-ray image capturing apparatus and image processing method Download PDFInfo
- Publication number
- CN102448374B CN102448374B CN201080022803.7A CN201080022803A CN102448374B CN 102448374 B CN102448374 B CN 102448374B CN 201080022803 A CN201080022803 A CN 201080022803A CN 102448374 B CN102448374 B CN 102448374B
- Authority
- CN
- China
- Prior art keywords
- image
- mentioned
- ray
- unit
- contrast
- 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
- 238000012545 processing Methods 0.000 title claims abstract description 65
- 238000003672 processing method Methods 0.000 title claims description 4
- 238000005304 joining Methods 0.000 claims abstract description 4
- 230000009466 transformation Effects 0.000 claims description 64
- 238000010586 diagram Methods 0.000 claims description 37
- 238000013459 approach Methods 0.000 claims description 14
- 238000000034 method Methods 0.000 claims description 14
- 238000000605 extraction Methods 0.000 claims description 10
- 238000003384 imaging method Methods 0.000 claims description 7
- 238000003860 storage Methods 0.000 claims description 5
- 239000000284 extract Substances 0.000 claims description 3
- 230000008569 process Effects 0.000 claims description 3
- 238000006243 chemical reaction Methods 0.000 abstract description 5
- 238000012937 correction Methods 0.000 description 13
- 230000002441 reversible effect Effects 0.000 description 9
- 238000003745 diagnosis Methods 0.000 description 6
- 210000003050 axon Anatomy 0.000 description 4
- 238000009795 derivation Methods 0.000 description 4
- 238000001514 detection method Methods 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 4
- 238000002372 labelling Methods 0.000 description 4
- 230000006870 function Effects 0.000 description 3
- 210000000988 bone and bone Anatomy 0.000 description 2
- 210000004027 cell Anatomy 0.000 description 2
- 230000001186 cumulative effect Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 230000005855 radiation Effects 0.000 description 2
- 230000003187 abdominal effect Effects 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 230000001154 acute effect Effects 0.000 description 1
- 230000004075 alteration Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 230000002349 favourable effect Effects 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 238000009499 grossing Methods 0.000 description 1
- 210000003141 lower extremity Anatomy 0.000 description 1
- 230000014759 maintenance of location Effects 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 230000005622 photoelectricity Effects 0.000 description 1
- 238000005096 rolling process Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 210000001519 tissue Anatomy 0.000 description 1
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/52—Devices using data or image processing specially adapted for radiation diagnosis
- A61B6/5211—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data
- A61B6/5229—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data combining image data of a patient, e.g. combining a functional image with an anatomical image
- A61B6/5235—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data combining image data of a patient, e.g. combining a functional image with an anatomical image combining images from the same or different ionising radiation imaging techniques, e.g. PET and CT
- A61B6/5241—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data combining image data of a patient, e.g. combining a functional image with an anatomical image combining images from the same or different ionising radiation imaging techniques, e.g. PET and CT combining overlapping images of the same imaging modality, e.g. by stitching
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/50—Image enhancement or restoration using two or more images, e.g. averaging or subtraction
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/90—Dynamic range modification of images or parts thereof
- G06T5/92—Dynamic range modification of images or parts thereof based on global image properties
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10116—X-ray image
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Medical Informatics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Pathology (AREA)
- Molecular Biology (AREA)
- Biophysics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Optics & Photonics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Radiology & Medical Imaging (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- High Energy & Nuclear Physics (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Image Processing (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
This invention provides an image processing apparatus capable of suppressing the variation in density and contrast of a long-length image that is a joint image obtained by joining a plurality of X-ray images. An image processing apparatus (6) comprises: an image memory (62) for storing the data of a plurality of X-ray images (20, 21, 22) partially overlapping with one another; a density/contrast correcting unit (64) for evening up the densities and contrasts of the plurality of X-ray images; and a long-length image obtaining unit (65) for joining all of the X-ray images (20, 21', 22") as evened up in density and contrast by the density/contrast correcting unit (64), thereby obtaining a long-length X-ray image corresponding to a long-length image-captured area of an examined body (M). The density/contrast correcting unit (64) comprises: a first table creating unit (642) for obtaining, based on image information of the overlapping portions (20a, 21a, 21b, 22a) of two adjacent ones (20, 21; or 21, 22) of the plurality of X-ray images, a gray scale conversion table (AB, BC) with which the densities and contrasts of the two X-ray images (20, 21; or 21, 22) can be evened up; a second table creating unit (643) for obtaining, based on the gray scale conversion table (AB, BC), a second gray scale conversion table (BA, CB) with which the densities and contrasts of a reference X-ray image (20) and other remaining X-ray images (21, 22) selected from the plurality of X-ray images (20-22) can be evened up; and an evening unit (644) for evening up the densities and contrasts of the remaining X-ray images (21, 22) on basis of the second gray scale conversion table (BA, CB).
Description
Technical field
The present invention relates to a kind of image processing apparatus, x-ray imaging apparatus and image processing method.
Background technology
When use has partly overlapping multiple x-ray image and obtains the long span image as index map picture, after the display density correspondence of the each image before making to engage in conjunction with and obtain long span image, this technology is in the public eye.For example, in the disclosed technology of patent documentation 1, first the each image before engaging is carried out to logarithmic transformation.Then obtain the block diagram of each image according to the each image after logarithmic transformation.Then detect the skew of the position of centre of gravity of each block diagram of obtaining, revise this skew.And in conjunction with the each image after offset correction, obtain the long span image as index map picture.
Patent documentation 1: JP 2007-275228 communique
Summary of the invention
But existing method is such technology: consistent for the display density of the joint integral image that makes to obtain, the concentration of the each image before engaging is adjusted into setting, and how does not pay close attention to for the contrast that engages integral image.
The line matter (energy size) of the x-ray using when the contrast of the x-ray image generally, obtaining by X-ray photographic depends on photography.Therefore, in the time that shooting has partly overlapping multiple image, in each photography, use while having changed the x-ray of line matter, though the concentration of the each image before can engaging by existing method adjustment, between the each image after concentration is adjusted, contrast is different sometimes.When contrast is different between each image, for certain bone portion, can be created in certain image and can clearly confirms, but the problem that cannot confirm in other images.
The problem that the present invention will solve is, provide a kind of at least can suppress as engaged multiple x-ray images index map picture long span image contrast uneven image processing apparatus and comprise the x-ray imaging apparatus of this device.
The present invention solves above-mentioned problem by following settling mode.And in following settling mode, the additional corresponding labelling of accompanying drawing that represents working of an invention mode is described, but this labelling is only for easily understanding invention, is not used as and limits invention.
The image processing apparatus (6) that invention relates to has: memory element (62), the data of multiple x-ray images (20,21,22) that storage area repeats; Long span image obtaining section (65), engages multiple x-ray images, obtains and the long span x-ray image corresponding to long span photography region of measured body (M).
In addition also there is correcting unit (64), it makes the processing corresponding to contrast (preferred concentration and contrast) of multiple x-ray images (20,21,22), long span image is obtained unit (65) and is engaged the x-ray image (20,21 ', 22 ") that makes contrast (preferred concentration and contrast) correspondence by correcting unit (64), obtains long span x-ray image.Correcting unit (64) has: form is obtained unit (642), according to the image information of the repeating part (20a, 21a, 21b, 22a) of two x-ray images (20,21 or 21,22) adjacent in multiple x-ray images, obtain the greyscale transformation form (AB, BC) corresponding to contrast (preferred concentration and contrast) that can make two x-ray images (20,21 or 21,22); Equalization unit (644), according to greyscale transformation form (AB, BC), makes the processing corresponding to contrast (preferred concentration and contrast) of multiple x-ray images (20～22).
The x-ray imaging apparatus (100) that relates to of invention has: image photography unit (Isosorbide-5-Nitrae), obtain the data of multiple x-ray images (20～22) that part repeats; Graphics processing unit (6), carries out joining process to multiple x-ray images; Display unit (7), the result of demonstration graphics processing unit (6).Graphics processing unit (6) has the image processing apparatus (6) that foregoing invention relates to.Display unit (7) shows obtains by long span image the long span x-ray image of obtaining unit (65).
According to foregoing invention, can suppress the inequality as the contrast of the long span image of index map picture.Result can obtain the long span x-ray image that is easy to diagnosis, contributes to improve diagnosis efficiency.
And, in prior art, when the position that measured body thickness is thicker is contained in long span photography region, because of the restriction of photography conditions, photography to the thicker position of this thickness need to be sacrificed image quality, but according to foregoing invention, the line matter that can change x-ray is taken the x-ray image that engages object, therefore can avoid sacrificing image quality, and suffer the angle of radiation also comparatively favourable from measured body (patient).
Accompanying drawing explanation
Fig. 1 is the block diagram that the entirety of the x-ray imaging apparatus that represents that an embodiment of the invention relate to forms.
Fig. 2 (a)～Fig. 2 (c) is the figure representing as an example of the x-ray image of the basic joint object of long span x-ray image.
Fig. 3 is the flow chart of the sequence of movement of the image processing part of presentation graphs 1.
Fig. 4 (a) is the block diagram in the repeat region 20a of the image 20 shown in Fig. 2 (a), and Fig. 4 (b) is the block diagram in the repeat region 21a of the image 21 shown in Fig. 2 (b).
Fig. 5 (a) will accumulate the stacked bar graph ratio of the block diagram of Fig. 4 (a), and Fig. 5 (b) will accumulate the stacked bar graph ratio of the block diagram of Fig. 4 (b).
Fig. 6 (a) and Fig. 6 (b) are the figure of one example of greyscale transformation form between 2 images that represent to generate by the step 323 of Fig. 3.
Fig. 7 is the figure that represents an example of the standardization greyscale transformation form generating by the step 332 of Fig. 3.
Fig. 8 is the figure that represents an example of the standardization greyscale transformation form generating by the step 332 of Fig. 3.
Fig. 9 is the figure of an example of the concentration contrast to handle averagely of the step 34 of presentation graphs 3.
Reference numeral
100X line photographic attachment, 1X spool (image photography unit), 4X thread detector (image photography unit), M measured body, 6 image processing parts (image processing apparatus), 61 image obtaining sections, 62 image storages (memory element), 63 bonding station configuration parts (bonding station setup unit), 64 concentration contrast correction portions (correcting unit), 641 logarithmic transformation handling parts (logarithmic transformation processing unit), 642 the 1st forms are made portion's (form is obtained unit), 642a block diagram obtaining section (block diagram is obtained unit), 642b characteristic quantity obtaining section (characteristic quantity is obtained unit), 642c form generating unit, 643 the 2nd forms are made portion's (the 2nd form is obtained unit), 643a benchmark image selection portion (benchmark image selected cell), 643b form generating unit, 644 to handle averagely portions (equalization unit), 65 long span image production parts (long span image is obtained unit), 7 display parts (display unit), 20～22X line image, 20a, 21a, 21b, 22a repeat region, 21 ', image after 22 ' greyscale transformation, 22 " image of gray scale after converting again, AB, greyscale transformation form (greyscale transformation form) between 2 images of BC, BA, CB reverse search form (the 2nd greyscale transformation form)
The specific embodiment
Referring to an embodiment of accompanying drawing explanation foregoing invention.
Diagnostic device in the such as hospital of x-ray imaging apparatus 100 of the present embodiment shown in Fig. 1 etc., has X-ray tube 1 (image photography unit) that roentgen radiation x uses and as plate etc. the x-ray detector 4 that sees through two-dimentional x-ray detector that x-ray image detection uses.This X-ray tube 1 and detector 4 are synchronized with each other, can form movably along the axon direction Z of measured body M.Make under the state that measured body M is static, X-ray tube 1 and detector 4 are repeatedly photographed in the time that the axon direction Z of measured body M moves, thereby can obtain along the continuous of the axon of measured body M and multiple x-ray images that part repeats.
Detector 4 (image photography unit) in the time of each photography, detect from X-ray tube 1 irradiate, and by measured body M see through x-ray image, its testing result is transformed to electric signal (x-ray detection signal), output to image processing part 6.
Image processing part 6 (image processing apparatus) is made up of CPU, ROM, RAM etc., multiple x-ray images of obtaining one by one according to the x-ray detection signal by exporting from detector 4, carry out the image processing (greyscale transformation, joint) for showing, the corresponding long span x-ray image in long span photography region (for example, region till from abdominal part to lower limb) of generation and measured body M.
Display part 7 (display unit) is accepted the output of image processing part 6, shows long span x-ray image.
In addition also can form as follows: the long span x-ray image generating by image processing part 6, under the state of Digital Image Data, is saved in the recording equipments such as magneto-optical disc apparatus, is sent to external equipment by network.
Then according to Fig. 1 of formation and the flow chart of Fig. 3 that image processing part 6 is shown with block diagram, the action (S1～S4) of image processing part 6 is described.
(S1)
First,, in the step shown in Fig. 3 (being designated hereinafter simply as " S ") 1, the data that engage multiple x-ray images of object are provided to image processing part 6 (engaging the supply of the view data of object) from x-ray detector 4.The view data that is provided to image processing part 6 is that the x-ray detection signal of exporting from x-ray detector 4 is transformed to Digital Image Data by A/D changer (omitting diagram).
In the present embodiment, for example provide shown in Fig. 2 (a)～Fig. 2 (c), along the continuous of the axon of measured body M and 3 x-ray images 20～22 that part repeats.And a part for image 20 and image 21 repeats, a part for image 21 and image 22 repeats.The picture number providing is not limited to 3, also more than 2 or 4.
In addition, to establish at this number of x-ray image that engages object be N, to them according to photography order-assigned 1 till the cognizance code of N.In the present embodiment, the image number that engages object is 3, therefore according to the cognizance code of photography order-assigned 1 to 3.Cognizance code is that the image 20 shown in Fig. 2 (a) is that the image 21 shown in 1, Fig. 2 (b) is that the image 22 shown in 2, Fig. 2 (c) is 3.
Turn back to Fig. 1 and Fig. 3.The image 20～22 that is provided to image processing part 6 is obtained by image obtaining section 61, temporarily stores in image storage 62 (memory element).
(S2)
Then in S2, store the data reading of 3 x-ray images of memorizer 62 into bonding station configuration part 63 (bonding station setup unit), set adjacent image bonding station (setting of bonding station) separately at this.
In the present embodiment, set as the image 20 of epigraph and as in the bonding station of image 21 of image, and set image 21 and the bonding station as the image 22 of hypograph.The setting of bonding station can be determined according to user's input, or such as, be decided by known processing (the judgement processing at the disclosed junction surfaces such as JP 2008-67916 communique).
In addition, labelling 20a, 21a in Fig. 2 (a) and Fig. 2 (b) refer to, have engaged the repeat region of the image 20,21 image 20,21 on the bonding station of being set by set positions portion 63.Labelling 21b, 22a in Fig. 2 (b) and Fig. 2 (c) too, refer to the repeat region that has engaged the image image of 21,22 o'clock 21,22 on the bonding station of being set by set positions portion 63.
(S3)
The long span x-ray image obtaining by the processing of image processing part 6 resembles an image in order to look on the whole, needs connection joining part smoothly.In order to connect smoothly, need to and make the operation that display density is corresponding according to each Image Adjusting gray scale (concentration).This is because if jointing edge boundary line interior of the long span x-ray image obtaining and have concentration difference outward is difficult to carry out image interpretation, diagnosis capability may decline.
But, only make display density correspondence also not enough, adjust engage before each image contrast, also make together display comparison degree corresponding with display density, this point is very important.When contrast is different even if concentration is identical, for certain bone portion, can be created in certain image and can clearly confirms, but the probability that cannot confirm in other images becomes the long span x-ray image that is difficult to diagnosis.
Therefore, in the present embodiment, as an embodiment of the invention, in order to suppress the display density of multiple x-ray images and the inequality of contrast before joint, concentration contrast correction portion 64 is set on image processing part 6, and the processing (processing of concentration contrast correction) that the inequality of carrying out having set by S2 the concentration contrast between the image of bonding station at this suppresses (S3).In addition, in the present embodiment, the below corrected concentrations of explanation and the mode of contrast, orientates the mode of the present invention of contrast correction as.
The concentration contrast correction of present embodiment is processed and is substantially carried out following processing: take the benchmark image selected as starting point, make the display density of adjacent with it image and contrast corresponding with display density and the contrast of benchmark image.Afterwards, on the bonding station that concentration and each image corresponding to contrast are set in S2, engage, generate long span x-ray image.
One example of concentration contrast correction processing (S31～S34) is below described.
(S31)
First, in the time having set each image 20～22 of bonding station and be provided to concentration contrast correction portion 64 (concentration contrast correction unit), in S31, by logarithmic transformation handling part 641 (logarithmic transformation processing unit), each image 20～22 is carried out to logarithmic transformation processing.
Generally, x-ray amount in the photon energy band of the x-ray using in body photography, that see through (inciding the line amount of x-ray detector), has the thickness of the material relatively seeing through and the characteristic that attenuation quotient logarithm reduces.Therefore, by carrying out logarithmic transformation processing, can offset the logarithm minimizing tendency from the information of x-ray detector, substitute simply is the thickness of material and the information that adds of attenuation quotient seeing through.This is the reason of carrying out logarithmic transformation.
(S32)
Then in S32, by the 1st form make portion 642 (form is obtained unit) make of making in adjacent x-ray image 20,21 with another image color and consistent 2 images of contrast between greyscale transformation form AB (greyscale transformation form).Simultaneously in the present embodiment, also make greyscale transformation form BC (greyscale transformation form) between 2 images of adjacent x-ray image 21,22.
Below explanation is used for an example (S321～S323) of the method for making form AB (with reference to Fig. 6 (a)).In addition, form BC (with reference to Fig. 6 (b)) also can make by same step, and therefore the description thereof will be omitted.
(S321)
First in S321, obtain block diagram by block diagram obtaining section 642a (block diagram is obtained unit), this histogram graph representation has been carried out the intensity (concentration, transverse axis) of picture element signal and the corresponding relation of the frequency of occurrences (frequency, the longitudinal axis) in repeat region 20a, the 21a of image 20,21 of logarithmic transformation.Fig. 4 (a) and Fig. 4 (b) are illustrated in an example of block diagram A, the B of this derivation.
Can understand following content according to block diagram A, the B of Fig. 4 (a) and Fig. 4 (b).First, the image color of the region 21a of image 21 is greater than the region 20a of image 20.This is because the concentration range (waveform starting position (minima) is to the scope of end position (maximum)) of block diagram B is compared and kept right with block diagram A, is positioned at high concentration one side.Secondly, the contrast of the region 20a of image 20 is greater than the region 21a of image 21.This is because the amplitude W2 of the concentration range of the amplitude W1 of the concentration range of block diagram A and block diagram B compares greatly.Hence one can see that, and image 20 is compared with image 21, and image color is lower, but contrast is larger.The correction of concentration and contrast is the consistent processing of block diagram that makes image 20,21.
(S322)
Turn back to Fig. 1 and Fig. 3.Then in S322, by the 642b of Characteristic Extraction portion (Characteristic Extraction unit), extract characteristic quantity A ', B ' from the each block diagram A, the B that are derived by obtaining section 642a.
Comprise any one or their combination in the minima, area average, maximum, intermediate value, mode of block diagram as an example of characteristic quantity A ', B ' of extracting.And, obtain the stacked bar graph of having accumulated each block diagram, make this stacked bar graph ratio (for example, with reference to Fig. 5 (a) and Fig. 5 (b)), can be using the concentration corresponding with multiple ratios as characteristic quantity A ', B '.Wherein, stacked bar graph ratio is referred to, (establish the concentration of image 20 as A1, A2 take concentration value ... the concentration of An, image 21 is B1, B2 ... Bn) be transverse axis, to have accumulated the ratio of value of the frequency of occurrences (frequency) of the concentration below this concentration value, ( image 20,21 is R1, R2 ... Rn) figure representing for the longitudinal axis.
In the present embodiment, while extracting characteristic quantity A ', B ' from each block diagram A, B, the multiple above-mentioned characteristic quantities of preferred compositions (such as maximum and maximum etc.), extract this assemblage characteristic amount.So, be conducive to make the greyscale transformation form of degree of precision.
(S323)
Then in S323, by form generating unit 642c, according to the each block diagram A being extracted by extraction unit 642b, characteristic quantity A ', the B ' of B, generate make by handling part 641 carried out in the image 20,21 of logarithmic transformation one with another image color and consistent 2 images of contrast between greyscale transformation form AB.
According to the characteristic quantity A ', the B ' that obtain from 2 images 20,21, the characteristic quantity of known any one image is corresponding with the characteristic quantity of any another image, but does not know that how corresponding the pixel value beyond characteristic quantity is.Therefore, need to carry out the pixel value between interpolation characteristic quantity according to multiple characteristic quantities, make the greyscale transformation form of variable pixel value of the image pixel value that is changed to another image.
In the present embodiment, while making form AB by above-mentioned characteristic quantity, derive the straight line that approaches of multiple characteristic quantities, carry out making of greyscale transformation form according to the straight line that approaches of this derivation.Make thus greyscale transformation form AB between 2 images shown in Fig. 6 (a), based on characteristic quantity.The derivation that approaches orthoscopic can be used simultaneous equations as being only at 2, makes the characteristic quantity number of middle use when uncertain when form, can use minimum to take advantage of the derivation such as method.The reason of approaching straight line in this use is, in the scope of the line matter (tube voltage 30kV～200kV) of the x-ray using in common medical diagnosis, when to as if when tissue, photoelectricity effect and Compton effect and attenuation coefficient can approach by exponential function, therefore in the time of the image using after logarithmic transformation, exponential function partial offset, pixel value and attenuation quotient and thickness with proportional, therefore can approach the straight line (being called exponential function rule in X-ray photographic method) of expression of first degree.
Wherein, in form is made, to approach all pixel values of straight-line interpolation with one, but the scope of the pixel value of interpolation repeat region 20a, 21a beyond time, preferably make form approach interpolation with pixel region and other straight lines of high pixel region of low pixel region, repeat region.Particularly, in low pixel region, for carrying out the interpolation of low pixel region of form, only use the characteristic quantity (approaching the characteristic quantity of several percentage ratios etc. of low pixel one side of minimum, meansigma methods, block diagram cumulative percentage) that approaches low pixel region, obtain and approach straight line.In the pixel region of repeat region, for carrying out interpolation, use whole characteristic quantities or the characteristic quantity except the characteristic quantity near low pixel region and high pixel region, obtain and approach straight line.
In high pixel region, for carrying out interpolation, use the characteristic quantity (characteristic quantities of several percentage ratios of high pixel value one side of peak, meansigma methods, block diagram cumulative percentage etc.) that approaches high pixel region, obtain and approach straight line.And, approach using each region as being obtained of using with each region that straight line carrys out interpolation and the form that connects is made.Now, for making the not acute variation of inclination of form of coupling part in each region, preferably use smoothing techniques (rolling average etc.).
(S33)
Then in S33, make portion 643 (the 2nd form is obtained unit) by the 2nd form, to the specific x-ray image of the benchmark as the processing of concentration contrast correction, make the standardization greyscale transformation form ABC (2nd greyscale transformation form) consistent with the image color of other images and contrast.
Below an example (S331, S332) of the method for standardization greyscale transformation form ABC is made in explanation.
(S331)
First,, in S331, by benchmark image selection portion 643a (benchmark image selected cell), select the benchmark image as the starting point of concentration contrast correction processing.The benchmark image of selecting is as long as any one in the image 20,21,22 having existed, or also can select image 20,21,22 to carry out other images of concentration contrast correction.The selection reference of benchmark image is not particularly limited, for example can be according to greyscale transformation form AB between make by the portion of making 642 2 images (with reference to Fig. 6 (a), BC (with reference to Fig. 6 (b)), to there is minima or peaked image is chosen as benchmark image, or will be chosen as benchmark image with the image of the average minimum (contrast is high) of the inclination of benchmark image to seasonable conversion form.The situation that following example is selected image 20 as benchmark image.
(S332)
Then in S332, by form generating unit 643b, generate make to select by selection portion 643a as the image 20 of benchmark image, the standardization greyscale transformation form consistent with the image color of the image 21,22 as other images and contrast.
In the present embodiment, as the 1st stage, first make the standardization greyscale transformation form of the image 21 of relative image 20.Because image 20 is benchmark images, thus the reverse search form BA (with reference to Fig. 7) of greyscale transformation form AB between 2 images making by S32 (with reference to Fig. 6 (a)) can be made, and use it.
As the 2nd stage, make the standardization gray scale form of the image 22 of relative image 20.Make the reverse search form CB (with reference to Fig. 8) of greyscale transformation form BC between 2 images making by S32 (with reference to Fig. 6 (b)), make the synthetic reverse search form CBA (omitting diagram) of it and above-mentioned reverse search form BA, and use it.
(S34)
Then in S34, by handle averagely portion 644 (equalization unit), according to the standardization greyscale transformation form of making, carry out the to handle averagely as the image 21,22 beyond the image 20 of benchmark image.
For example as shown in Figure 9, first whole pixel values of image 21 are transformed to the gray scale of the reverse search form BA based on making by S332, as image 21 '.Make thus concentration and the contrast consistent (image 21 → image 21 ') of image 20 and the image 21 ' after the greyscale transformation adjacent with its downstream.Then, whole pixel values of image 22 are transformed to the gray scale of the synthetic reverse search form CBA based on making by S332, as image 22 ".So, make the image 22 after image 21 ' after greyscale transformation and the gray scale adjacent with its downstream convert again " concentration and contrast consistent (image 22 → image 22 ").
In addition image 22 → image 22, " greyscale transformation also can be undertaken by 2 stages.
For example, first whole pixel values of image 22 are transformed to the gray scale of the reverse search form CB based on making by S332, as image 22 '.So, make concentration and the contrast consistent (image 22 → image 22 ') of the image 22 ' after image 21 before greyscale transformation and the greyscale transformation adjacent with its downstream.Then, whole pixel values of image 22 ' are transformed to the gray scale of the reverse search form BA based on making by S332, as image 22 " (image 22 ' → image 22 ").
Through above S31～S34, can obtain image corresponding to concentration and contrast 20,21 ', 22 " (with reference to Fig. 9).
(S4)
Then in S4, by long span image production part 65 (long span image is obtained unit), on the bonding station of setting by S2, engage image corresponding to concentration and contrast 20,21 ', 22 ", make the long span image as an image.
The joint of each image as shown in known method, can the arbitrary region weighting to repeat region 20a, 21a, 21b, 22a and proofread and correct after engage, or can use adjacent image 20,21 ' or 21 ', 22 " an image.
In addition, generally, the long span image of making outputs to image output device, but this output device can be monitor, film printer, also can be used as data and is saved in server, recording medium.
In the present embodiment, the long span image of making outputed to display part 7 and show process ends.
According to present embodiment, engaging before multiple x-ray images 20～22, make display density and the contrast correspondence of each image, the long span x-ray image after therefore engaging can not produce the inequality of concentration and contrast.Consequently, can obtain the long span x-ray image of easy diagnosis.
Embodiment described above, for easily understanding foregoing invention, is not used in restriction foregoing invention.Therefore, in above-mentioned embodiment, disclosed each key element also comprises all design alterations, the equipollent in the technical scope that belongs to foregoing invention.
Claims (8)
1. an image processing apparatus, has: memory element, the data of multiple x-ray images that storage area repeats; Long span image is obtained unit, engages above-mentioned multiple x-ray image, obtains and the long span x-ray image corresponding to long span photography region of measured body, it is characterized in that,
Have correcting unit, it makes the processing corresponding to contrast of above-mentioned multiple x-ray images, and above-mentioned correcting unit has:
Form is obtained unit, and it has: logarithmic transformation processing unit, carries out logarithmic transformation processing to the repeating part of two x-ray images adjacent in above-mentioned multiple x-ray images; Block diagram is obtained unit, from the image of the repeating part of two x-ray images through logarithmic transformation, obtains the block diagram of the corresponding relation that represents pixel concentration and occurrence rate; Characteristic Extraction unit, from least one characteristic quantity of above-mentioned histogram extraction, above table is obtained unit according to above-mentioned characteristic quantity, generates the contrast and the greyscale transformation form of concentration to seasonable use that make above-mentioned two x-ray images;
Equalization unit, according to above-mentioned greyscale transformation form, makes contrast and the processing corresponding to concentration of above-mentioned multiple x-ray images,
Above-mentioned long span image is obtained unit joint and is made contrast and x-ray image corresponding to concentration by above-mentioned correcting unit, obtains above-mentioned long span x-ray image.
2. image processing apparatus according to claim 1, is characterized in that,
More than above-mentioned Characteristic Extraction unit extracts any one in the minima, area average, intermediate value, maximum, mode of above-mentioned block diagram, as above-mentioned characteristic quantity.
3. image processing apparatus according to claim 1, is characterized in that,
Above-mentioned characteristic quantity is extracted by the ratio of the accumulated value of above-mentioned block diagram.
4. image processing apparatus according to claim 1, is characterized in that,
Above table is obtained unit in the time obtaining above-mentioned greyscale transformation form, to low pixel region and high pixel region in the image beyond above-mentioned repeating part, uses and utilizes the straight line that approaches that method of least square calculates to carry out linear interpolation.
5. image processing apparatus according to claim 1, is characterized in that,
Above-mentioned correcting unit has benchmark x-ray image selected cell, and from above-mentioned multiple x-ray image selection reference x-ray images, above-mentioned equalization unit makes the contrast processing corresponding with the contrast of said reference x-ray image of the image beyond said reference x-ray image.
6. image processing apparatus according to claim 5, is characterized in that,
Above-mentioned correcting unit has the 2nd form and obtains unit, it is according to above-mentioned greyscale transformation form, obtain the 2nd greyscale transformation form corresponding to contrast that can make said reference x-ray image and said reference x-ray image image in addition, above-mentioned equalization unit, according to above-mentioned the 2nd greyscale transformation form, makes the contrast of image and the processing corresponding to contrast of said reference x-ray image in addition of said reference x-ray image.
7. an x-ray imaging apparatus, has: image photography unit, obtain the data of multiple x-ray images that part repeats; Graphics processing unit, carries out joining process to above-mentioned multiple x-ray images; Display unit, shows and it is characterized in that the result of above-mentioned graphics processing unit,
Above-mentioned graphics processing unit has:
Memory element, the data of storing above-mentioned multiple x-ray images;
Correcting unit, makes the processing corresponding to contrast of above-mentioned multiple x-ray images;
Long span image is obtained unit, engages and makes x-ray image corresponding to contrast by above-mentioned correcting unit, obtain and the long span x-ray image corresponding to long span photography region of measured body,
Above-mentioned correcting unit has:
Form is obtained unit, and it has: logarithmic transformation processing unit, carries out logarithmic transformation processing to the repeating part of two x-ray images adjacent in above-mentioned multiple x-ray images; Block diagram is obtained unit, from the image of the repeating part of two x-ray images through logarithmic transformation, obtains the block diagram of the corresponding relation that represents pixel concentration and occurrence rate; Characteristic Extraction unit, from least one characteristic quantity of above-mentioned histogram extraction, above table is obtained unit according to above-mentioned characteristic quantity, generates the contrast and the greyscale transformation form of concentration to seasonable use that make above-mentioned two x-ray images;
Equalization unit, according to above-mentioned greyscale transformation form, makes contrast and the processing corresponding to concentration of above-mentioned multiple x-ray images,
Above-mentioned display unit shows obtains by above-mentioned long span image the above-mentioned long span x-ray image that unit is obtained.
8. an image processing method, has: storing step, the data of multiple x-ray images that storage area repeats; Long span image is obtained step, engages above-mentioned multiple x-ray image, obtains and the long span x-ray image corresponding to long span photography region of measured body, it is characterized in that,
Have aligning step, it makes the processing corresponding to contrast of above-mentioned multiple x-ray images, and above-mentioned aligning step has:
Form is obtained step, and it has: logarithmic transformation treatment step, carries out logarithmic transformation processing to the repeating part of two x-ray images adjacent in above-mentioned multiple x-ray images; Block diagram is obtained step, from the image of the repeating part of two x-ray images through logarithmic transformation, obtains the block diagram of the corresponding relation that represents pixel concentration and occurrence rate; Characteristic Extraction step, from least one characteristic quantity of above-mentioned histogram extraction, above table is obtained step according to above-mentioned characteristic quantity, generates the contrast and the greyscale transformation form of concentration to seasonable use that make above-mentioned two x-ray images;
Equalization step, according to above-mentioned greyscale transformation form, makes contrast and the processing corresponding to concentration of above-mentioned multiple x-ray images,
Above-mentioned long span image is obtained step joint and is made contrast and x-ray image corresponding to concentration by above-mentioned aligning step, obtains above-mentioned long span x-ray image.
Applications Claiming Priority (3)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2009127273 | 2009-05-27 | ||
| JP2009-127273 | 2009-05-27 | ||
| PCT/JP2010/058264 WO2010137482A1 (en) | 2009-05-27 | 2010-05-17 | Image processing apparatus, x-ray image capturing apparatus and image processing method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN102448374A CN102448374A (en) | 2012-05-09 |
| CN102448374B true CN102448374B (en) | 2014-05-14 |
Family
ID=43222593
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201080022803.7A Expired - Fee Related CN102448374B (en) | 2009-05-27 | 2010-05-17 | Image processing apparatus, X-ray image capturing apparatus and image processing method |
Country Status (3)
| Country | Link |
|---|---|
| JP (1) | JP5670327B2 (en) |
| CN (1) | CN102448374B (en) |
| WO (1) | WO2010137482A1 (en) |
Families Citing this family (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP5759405B2 (en) * | 2012-03-22 | 2015-08-05 | 富士フイルム株式会社 | Luminance unevenness detection apparatus and method |
| JP2015047194A (en) * | 2013-08-30 | 2015-03-16 | 株式会社日立メディコ | X-ray diagnostic apparatus and data processing method by the same |
| JP6141995B2 (en) * | 2013-09-27 | 2017-06-07 | 富士フイルム株式会社 | Mammography apparatus, radiographic imaging method and program |
| CA3065957A1 (en) * | 2017-08-14 | 2019-02-21 | Raytheon Company | Subtraction algorithm for detection of tumors |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1278650C (en) * | 2002-07-29 | 2006-10-11 | 株式会社东芝 | X-ray image processing device and method |
Family Cites Families (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| DE10135427A1 (en) * | 2001-07-20 | 2003-02-13 | Siemens Ag | Areal image detector for electromagnetic rays, especially X-rays |
| JP4533587B2 (en) * | 2003-02-07 | 2010-09-01 | 株式会社東芝 | Medical image laminating device |
| JP2005046444A (en) * | 2003-07-30 | 2005-02-24 | Hitachi Medical Corp | Medical x-ray apparatus |
-
2010
- 2010-05-17 CN CN201080022803.7A patent/CN102448374B/en not_active Expired - Fee Related
- 2010-05-17 WO PCT/JP2010/058264 patent/WO2010137482A1/en active Application Filing
- 2010-05-17 JP JP2011515979A patent/JP5670327B2/en not_active Expired - Fee Related
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1278650C (en) * | 2002-07-29 | 2006-10-11 | 株式会社东芝 | X-ray image processing device and method |
Non-Patent Citations (1)
| Title |
|---|
| JP特开2005-46444A 2005.02.24 |
Also Published As
| Publication number | Publication date |
|---|---|
| JP5670327B2 (en) | 2015-02-18 |
| WO2010137482A1 (en) | 2010-12-02 |
| CN102448374A (en) | 2012-05-09 |
| JPWO2010137482A1 (en) | 2012-11-12 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP6071144B2 (en) | Radiation image analysis apparatus and method, and program | |
| CN102448374B (en) | Image processing apparatus, X-ray image capturing apparatus and image processing method | |
| CN102254310B (en) | For the method reducing image artifacts in CT view data, especially metal artifacts | |
| RU2589383C2 (en) | Two-mode imaging including quality metrics | |
| EP2638858B1 (en) | Diagnostic imaging device and method | |
| CN101930595A (en) | Image processing method and image processing equipment | |
| CN106683146A (en) | Image reconstruction method and parameter determining method of image reconstruction algorithm | |
| US7313225B2 (en) | Method and apparatus for generating of a digital x-ray image | |
| US10748319B1 (en) | Composite radiographic image that corrects effects of parallax distortion | |
| US9486177B2 (en) | X-ray computed tomography apparatus, medical image processing apparatus, and medical image processing method | |
| CN103156629A (en) | Image processing apparatus, image processing method, and program | |
| DE102017200282B3 (en) | Method for reducing image artifacts | |
| CN106296615A (en) | CT image corrects the method and system of metal artifacts | |
| US20160134852A1 (en) | System and method for multi-material correction of image data | |
| CN106551702A (en) | The beam hardening artifact correcting unit of CT images and method | |
| KR102126355B1 (en) | Radiography imaging apparatus and method for generating an radiographic image | |
| Shi et al. | A Virtual Monochromatic Imaging Method for Spectral CT Based on Wasserstein Generative Adversarial Network With a Hybrid Loss. | |
| CN110270015B (en) | sCT generation method based on multi-sequence MRI | |
| KR101146006B1 (en) | Method and apparatus to remove ring artifacts in x-ray CT | |
| US10650561B2 (en) | Composite radiographic image that corrects effects of parallax distortion | |
| WO2013191220A1 (en) | Diagnostic imaging device, and control method thereof | |
| US9959614B2 (en) | Image processing device and method | |
| JP5305687B2 (en) | X-ray video imaging system | |
| Luckner et al. | Towards full-body x-ray images | |
| KR101676304B1 (en) | Image refinement method and computer readable storage medium using spr |
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 | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20170306 Address after: Tokyo, Japan, Japan Patentee after: Hitachi Ltd. Address before: Tokyo, Japan Patentee before: Hitachi Medical Corporation |
|
| TR01 | Transfer of patent right | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20140514 Termination date: 20170517 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |