WO2010050031A1 - 放射線撮影装置 - Google Patents

放射線撮影装置 Download PDF

Info

Publication number
WO2010050031A1
WO2010050031A1 PCT/JP2008/069781 JP2008069781W WO2010050031A1 WO 2010050031 A1 WO2010050031 A1 WO 2010050031A1 JP 2008069781 W JP2008069781 W JP 2008069781W WO 2010050031 A1 WO2010050031 A1 WO 2010050031A1
Authority
WO
WIPO (PCT)
Prior art keywords
image
intermediate image
pixel value
coefficient
histogram
Prior art date
Application number
PCT/JP2008/069781
Other languages
English (en)
French (fr)
Inventor
滝人 酒井
孝之 佐野
Original Assignee
株式会社島津製作所
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by 株式会社島津製作所 filed Critical 株式会社島津製作所
Priority to PCT/JP2008/069781 priority Critical patent/WO2010050031A1/ja
Publication of WO2010050031A1 publication Critical patent/WO2010050031A1/ja

Links

Images

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment

Definitions

  • the present invention irradiates a subject with two types of radiation beams, generates two types of fluoroscopic images, and subtracts the other fluoroscopic image from one radioscopic image while performing predetermined weighting on the images.
  • the present invention also relates to a radiation imaging apparatus capable of performing subtraction imaging for generating a composite fluoroscopic image suited to the purpose of diagnosis.
  • the radiation imaging apparatus irradiates a subject with a radiation beam, detects the radiation transmitted therethrough, and generates a radiographic image.
  • the soft tissue of the subject may be reflected by simply taking a radioscopic image, and detailed diagnosis may not be possible.
  • Radiation includes those with wavelengths that are likely to pass through bones and soft tissues, and those that do not. That is, the ease of transmission of radiation differs between the bone portion and the soft tissue depending on the wavelength. Further, when the output of the radiation beam is changed, the wavelength spectrum included in the radiation beam is changed. By utilizing this, it is possible to synthesize a synthetic fluoroscopic image (subtraction image) in which the bone portion of the subject is emphasized (see, for example, Patent Document 1).
  • the subject is irradiated with a high-power radiation beam.
  • generated at this time be a high output image.
  • the subject is irradiated with a low-power radiation beam.
  • generated at this time be a low output image.
  • the pixel values of the high output image and the low output image are multiplied by a predetermined coefficient, and the low output image is subtracted from the high output image, for example, radiation having a wavelength suitable for contrasting the bone of the subject is enhanced.
  • a synthetic fluoroscopic image in which the bone part of the specimen is clearly reflected is obtained.
  • the pixel values of the high-output image and the low-output image are multiplied by other coefficients and the above subtraction is performed, for example, radiation having a wavelength suitable for contrasting the soft tissue of the subject is emphasized, and the subject A synthetic fluoroscopic image in which the soft tissue is clearly reflected can be obtained.
  • the composite perspective image is displayed on a monitor attached to the radiation imaging apparatus.
  • the conventional example having such a configuration has the following problems. That is, according to the conventional configuration, after the first synthetic fluoroscopic image is displayed, when the coefficient is reselected and the second synthetic fluoroscopic image is generated, the diagnosis using the second synthetic fluoroscopic image starts immediately. Can not. It is assumed that the first composite perspective image is reflected on the monitor. The adjustment of the brightness of the radiographic image displayed on the monitor is suitable for visualizing the first synthetic fluoroscopic image, for example, when the operator operates a knob attached to the monitor. When the second synthetic fluoroscopic image is displayed on the monitor in this state, the second synthetic fluoroscopic image becomes extremely dark or extremely bright, and the operator must adjust the brightness of the radioscopic fluoroscopic image again. I must.
  • the operator since it is not guaranteed that the luminance distribution of the first synthetic fluoroscopic image and the luminance distribution of the second synthetic fluoroscopic image are the same, the operator is required to display the synthetic fluoroscopic image every time the synthetic fluoroscopic image is displayed on the monitor. It is necessary to adjust the brightness.
  • the coefficients multiplied by the high output image and the low output image are changed, the luminance distribution of the combined fluoroscopic image is naturally changed. This is because when the coefficient is changed, the pixel value subtracted from the high-output image is changed accordingly.
  • the present invention has been made in view of such circumstances, and an object of the present invention is to provide a radiation imaging apparatus that does not require adjustment of the brightness of the monitor each time a synthesized fluoroscopic image (subtraction image) is updated and is easy to operate. Is to provide.
  • the present invention has the following configuration in order to achieve such an object. That is, the present invention provides a radiation detection apparatus comprising: a radiation source that irradiates a radiation beam with a predetermined output; a radiation detection means that detects the radiation beam; The radiographic image generating means for converting the detection data output by the radiographic image, the coefficient storage means for storing the first coefficient and the second coefficient, and multiplying one of the radiographic images by the first coefficient
  • a first intermediate image generating unit configured to generate a first intermediate image and multiply the other radioscopic image by a second coefficient to generate the other first intermediate image; and from the first intermediate image to the other first
  • a second intermediate image generating unit that generates a second intermediate image by subtracting the intermediate image
  • a luminance adjusting unit that adjusts the luminance of the second intermediate image to generate a composite perspective image. It is an butterfly.
  • the configuration of the present invention includes a radiographic image generation unit that converts detection data output from the radiation detection unit into a radiographic image, and a coefficient storage unit that stores a first coefficient and a second coefficient.
  • a first intermediate image that multiplies the first fluoroscopic image by a first coefficient and generates one first intermediate image, and multiplies the second radiographic image by a second coefficient to generate the other first intermediate image.
  • a generating unit; and a second intermediate image generating unit configured to generate a second intermediate image by subtracting the other first intermediate image from the first intermediate image.
  • the second intermediate image is subjected to brightness correction before being used for diagnosis. Since this brightness correction is performed each time the second intermediate image is generated, even if the second intermediate image is regenerated by changing the first coefficient and the second coefficient, the composite fluoroscopy used for diagnosis is used. The brightness of the image has already been adjusted. Therefore, the operator does not need to adjust the brightness of the monitor each time a composite perspective image is taken. Therefore, a radiographic apparatus that can be easily operated by the operator can be provided.
  • the luminance adjusting means includes a histogram generating means for generating a histogram of the second intermediate image, a dividing means for dividing the entire pixel value range of the histogram into a main pixel value range and an unnecessary pixel value range, and a second intermediate Conversion means for generating a composite perspective image by successively converting each of the pixel values of the pixels constituting the image into new pixel values, and by the conversion means, among the pixel values in the second intermediate image,
  • the pixel values belonging to the main pixel value range can be assigned to the pixel values belonging to the target pixel value range in the composite perspective image.
  • the above configuration shows one of the specific modes of brightness adjustment in the present invention. That is, the above configuration generates a histogram of the second intermediate image, and divides the entire pixel value range of the histogram into a main pixel value range and an unnecessary pixel value range. Then, each of the pixel values of the pixels constituting the second intermediate image is successively converted into new pixel values, thereby generating a composite perspective image. According to the above configuration, since a histogram can be generated for each second intermediate image, luminance correction using the histogram is ensured every time no matter how many times the generation of the composite perspective image is repeated. Become.
  • pixel values belonging to the main pixel value range are allocated to pixel values belonging to the target pixel value range in the composite perspective image.
  • the above-described histogram generation means can be configured to generate a histogram based on pixels constituting a partial region of the second intermediate image.
  • the coefficient storage means further stores a pixel value correspondence table in which the pixel values in the second intermediate image and the pixel values in the composite perspective image correspond to each other in association with the first coefficient and the second coefficient, and brightness adjustment means Obtains the luminance correspondence table associated with the first coefficient and the second coefficient, and refers to the luminance correspondence table to make each of the pixel values of the pixels constituting the second intermediate image one after another as a new pixel value. It can also be set as the structure which produces
  • the above configuration shows one of the specific modes of brightness adjustment in the present invention. That is, a synthesized perspective image is generated by successively setting each pixel value of the pixels constituting the second intermediate image to a new pixel value. At this time, a luminance correspondence table is used. By doing so, a configuration for generating a histogram can be omitted, and a radiation imaging apparatus with improved image processing speed can be provided.
  • the present invention may further include a combining unit that combines a plurality of synthetic fluoroscopic images to generate a moving image.
  • a moving image suitable for diagnosis can be generated.
  • a histogram may be generated every time the second intermediate image is acquired to generate a combined perspective image.
  • the luminance of the new synthetic fluoroscopic image can be adjusted again.
  • the composition of the present invention can acquire a synthetic perspective image (subtraction image).
  • the second intermediate image is subjected to brightness correction before being used for diagnosis. Since this brightness correction is performed every time the second intermediate image is generated, even if the setting is changed and the second intermediate image is regenerated, the brightness adjustment of the composite fluoroscopic image used for diagnosis is Already done. Therefore, the operator does not need to adjust the brightness of the monitor each time a composite perspective image is taken. Therefore, a radiographic apparatus that can be easily operated by the operator can be provided.
  • FIG. 1 is a functional block diagram illustrating a configuration of an X-ray imaging apparatus according to Embodiment 1.
  • FIG. 3 is a flowchart for explaining the operation of the X-ray imaging apparatus according to Embodiment 1; 6 is a schematic diagram for explaining a histogram according to Embodiment 1.
  • FIG. FIG. 6 is a schematic diagram for explaining the operation of the conversion unit according to the first embodiment. It is a functional block diagram explaining the structure of the X-ray imaging apparatus which concerns on 1 modification of this invention. It is a schematic diagram explaining the structure of the X-ray imaging apparatus which concerns on 1 modification of this invention.
  • FIG. 1 is a functional block diagram illustrating the configuration of the X-ray imaging apparatus 1 according to the first embodiment.
  • the X-ray imaging apparatus 1 according to the first embodiment includes a top plate 2 on which a subject M is placed and an X-ray beam that irradiates a pulsed X-ray beam adjacent to the top plate 2.
  • An X-ray tube 3, a flat panel detector (FPD) 4 for detecting X-rays transmitted through the subject, and an X-ray grid 5 for removing scattered X-rays incident on the FPD 4 are provided.
  • FPD flat panel detector
  • the configuration of the first embodiment includes an X-ray tube control unit 6 that controls the tube voltage, tube current, and X-ray beam pulse width of the X-ray tube 3.
  • the X-ray imaging apparatus 1 includes an X-ray fluoroscopic image generation unit 11 that converts detection data output from the FPD 4 into an X-ray fluoroscopic image, a coefficient storage unit 12 that stores various coefficients described later, and an X-ray fluoroscopy.
  • a log conversion processing unit 13 that performs log conversion processing on an image, a weighting unit 14 that multiplies the X-ray fluoroscopic image by a coefficient, a subtraction processing unit 15 that subtracts the other X-ray fluoroscopic image from one X-ray fluoroscopic image, Histogram generation unit 16 that generates a histogram S of an X-ray fluoroscopic image, a division unit 17 that divides a pixel value range of the histogram S, a conversion unit 18 that converts a pixel value to generate a composite fluoroscopic image, and a synthetic fluoroscopic image Are combined with each other to generate a moving image.
  • the X-ray tube corresponds to the radiation source of the present invention
  • the FPD corresponds to the radiation detection means of the present invention
  • the X-ray tube control unit corresponds to the radiation source control unit of the present invention
  • the X-ray fluoroscopic image generation unit corresponds to the radiographic image generation unit of the present invention.
  • each of the coefficient storage unit, the weighting unit, and the subtraction processing unit corresponds to a coefficient storage unit, a first intermediate image generation unit, and a second intermediate image generation unit of the present invention.
  • each of the histogram generation unit, division unit, conversion unit, and combination unit corresponds to the histogram generation unit, division unit, conversion unit, and combination unit of the present invention.
  • the histogram generating means, the dividing means, and the converting means correspond to specific modes of the luminance correcting means of the present invention.
  • the X-ray imaging apparatus 1 includes an operation console 23 that receives an instruction from an operator, and a display unit 22 that displays an X-ray fluoroscopic image or a moving image.
  • the X-ray imaging apparatus 1 includes an X-ray tube control unit 6, an X-ray fluoroscopic image generation unit 11, a coefficient storage unit 12, a log conversion processing unit 13, a weighting unit 14, a subtraction processing unit 15, a histogram generation unit 16,
  • a main control unit 25 that controls the division unit 17, the conversion unit 18, and the coupling unit 19 is provided.
  • the main control unit 25 is constituted by a CPU, and realizes each unit by executing various programs. Further, each of the above-described units may be divided and executed by an arithmetic device that takes charge of them.
  • FIG. 2 is a flowchart for explaining the operation of the X-ray imaging apparatus 1 according to the first embodiment.
  • a placement step S1 for placing the subject on the top board, and high-output radiation on the subject.
  • a high output image generation step S2 for generating a high output image H by irradiating a beam; a low output image generation step S3 for generating a low output image L by irradiating the subject with a low output radiation beam;
  • a log conversion step S4 for log-converting the output image H and the low output image L, a high coefficient image H, and a low output image L are multiplied by a predetermined coefficient to generate a first intermediate image P1.
  • these steps will be described in order.
  • ⁇ Installation step S1, high output image generation step S2, low output image generation step S3> First, the subject M is placed on the top 2 and preparation for X-ray fluoroscopic imaging of the subject M is performed.
  • the X-ray tube control unit 6 determines the tube voltage, tube current, and X-ray beam irradiation time of the X-ray tube 3 so as to irradiate a high-power X-ray beam. Irradiation with a high-power X-ray beam.
  • the X-ray tube control unit 6 determines the tube voltage, tube current, and X-ray beam irradiation time of the X-ray tube 3 so as to irradiate the low-power X-ray beam, and the X-ray tube 3 Irradiate a low-power X-ray beam.
  • the X-ray beam at this time has a pulse shape.
  • the FPD 4 detects the X-ray beam that has passed through the subject M, and outputs a detection signal to the X-ray fluoroscopic image generation unit 11.
  • the X-ray fluoroscopic image generator 11 generates an X-ray fluoroscopic image based on the detection signal.
  • the X-ray fluoroscopic image generated based on the detection signal obtained by irradiating the subject M with the high output X-ray beam is the high output image H, and the subject M is irradiated with the low output X-ray beam.
  • the X-ray fluoroscopic image generated based on the detection signal obtained in this manner is a low output image L.
  • the high output image H and the low output image L are sent to the log conversion processing unit 13 at the subsequent stage.
  • One of the high output image and the low output image corresponds to one X-ray fluoroscopic image of the present invention, and the other of either the high output image or the low output image is the book. This corresponds to the other X-ray fluoroscopic image of the invention.
  • the high output X-ray beam includes the X-rays of various wavelengths.
  • a high-power X-ray beam includes, for example, more X-rays having a specific wavelength than a low-power X-ray beam.
  • the transmission characteristics when the high-power X-ray beam and the low-power X-ray beam pass through the subject M are different from each other. Therefore, in order to obtain a fluoroscopic image based on the transmittance of X-rays having a specific wavelength contained in a high-power X-ray beam, the low-power image L may be subtracted from the high-power image H.
  • the wavelength of X-rays to be emphasized varies depending on the purpose of inspection. Therefore, by multiplying the high output image H and the low output image L by a predetermined coefficient and then subtracting the low output image L from the high output image H, the wavelength of the X-ray to be enhanced can be selected. . This time, such a subtraction process is performed.
  • ⁇ Log conversion step S4> Prior to the subtraction process described above, the high output image H and the low output image L are sent to the log conversion processing unit 13.
  • the high output image H is converted into log conversion data logH
  • the low output image L is converted into log conversion data logL.
  • the reason for this processing is that an exponential relationship is superimposed on the X-ray wavelength and the transmittance of the X-ray subject. That is, if log conversion processing is performed once before subtracting both images H and L, the exponential relationship that has a strong tendency is eliminated in advance in both log conversion data logH and logL. Therefore, if the log conversion data is subtracted, the difference between the high output image H and the low output image L becomes more obvious.
  • the log conversion data logH is multiplied by a high output coefficient, and the log conversion data logL is multiplied by a low output coefficient.
  • the log conversion data logH is converted into the first intermediate image H1
  • the log conversion data logL is converted into the first intermediate image L1.
  • the first intermediate images H1 and L1 are still in a state where log conversion has been performed.
  • the first intermediate image H1 and the first intermediate image L1 are sent to the subtraction processing unit 15.
  • One of the low output coefficient and the high output coefficient corresponds to the first coefficient of the present invention, and the other corresponds to the second coefficient of the present invention.
  • the subtraction processing unit 15 performs a subtraction process for subtracting the first intermediate image L1 from the first intermediate image H1, thereby generating a second intermediate image P2.
  • the second intermediate image P2 is obtained by performing a log reverse conversion process after the subtraction process in the subtraction processing unit 15.
  • ⁇ Division step S6> characteristic brightness correction processing in the present invention will be described. That is, a new image is generated by converting the pixel values of the pixels constituting the second intermediate image P2.
  • the second intermediate image P2 is sent to the histogram generation unit 16 to generate the histogram S of the second intermediate image P2.
  • the histogram S can be represented by a graph in which the pixel value a and the frequency f are related as shown in FIG.
  • the pixel value a takes a value from 0 to 4095, and the frequency f represents how often each pixel value a appears in the second intermediate image P2.
  • the frequency f is biased only to the region where the pixel value is 0, and the frequency f of the pixel values from 1 to 4095 is 0.
  • the frequency f is within the entire range (the total pixel value range r0) of the pixel value a as shown in FIG. It is biased to a certain area (main pixel value range r1).
  • the generated histogram S is sent to the dividing unit 17.
  • the entire pixel value range r0 is divided into a main pixel value range r1 and unnecessary pixel value ranges r2 and r3. This actual processing will be described.
  • the dividing unit 17 reads the frequency f when the pixel value is 0 from the histogram S. Then, the frequency f at that time is read while increasing the pixel value by one. The frequency f read at this time is added up and expressed as a cumulative frequency F. Then, the ratio of the cumulative frequency F to the total frequency f in the entire pixel value range r0 is calculated.
  • the dividing unit 17 finishes reading the frequency f when the ratio of the cumulative frequency F reaches a predetermined set value, and reads the pixel value a at this time.
  • the pixel value at this time is A1.
  • an area between pixel values 0 and A1 is determined as an unnecessary pixel value range r2. Since the set value is set low, the frequency with which the pixel values belonging to the unnecessary pixel value range r2 appear in the second intermediate image is low.
  • the cumulative frequency F is also obtained for the high pixel value side. That is, the frequency f when the pixel value a is 4095 is read from the histogram S, and thereafter, the cumulative frequency F is obtained while decreasing the pixel value by one.
  • the dividing unit 17 finishes reading the frequency f, and sets the pixel value at this time as A2.
  • an area where the pixel value is between A2 and 4095 is determined as the unnecessary pixel value range r3.
  • the pixel values A1 and A2 are sent to the conversion unit 18.
  • the second intermediate image P ⁇ b> 2 is sent from the subtraction processing unit 15 to the conversion unit 18.
  • ⁇ Pixel value changing step S7> The conversion unit 18 recognizes the main pixel value range and the unnecessary pixel value range in the histogram S from the transmitted pixel values A1 and A2. Then, the pixel values of the pixels constituting the second intermediate image P2 are converted based on the main pixel value range r1 and the unnecessary pixel value ranges r2 and r3. Specifically, among the pixel values in the second intermediate image P2, pixel values belonging to the main pixel value range r1 are allocated by the conversion unit 18 to pixel values belonging to the target pixel value range t in the composite perspective image. .
  • FIG. 4 is a schematic diagram illustrating the operation of the conversion unit according to the first embodiment.
  • the upper array in FIG. 4 represents the entire pixel value range r0 in the second intermediate image P2, and pixel values from 0 to 4,095 are prepared.
  • the pixel value A1 is assumed to be 1,000, and the pixel value A2 is assumed to be 2,000. Therefore, the width of the main pixel value range r1 is 1,000 pixel values.
  • the lower array in FIG. 4 represents the entire pixel value range in the combined perspective image to be generated by the conversion unit 18, and the pixel values are from 0 as in the second intermediate image P2 for the sake of brevity. Assume that 4,095 have been prepared. However, actually, the configuration is not limited to this.
  • a predetermined pixel value range is set in the entire pixel value range in the combined perspective image. That is, the target pixel value range t sandwiched between the lower limit value B1 and the upper limit value B2.
  • the lower limit value B1 is 50
  • the upper limit value B2 is 4,050. Therefore, the width of the target pixel value range t is 4,000 pixel values.
  • the conversion unit 18 reads the image value corresponding to each pixel constituting the second intermediate image P2, and determines whether it is in the main pixel value range r1 or in the unnecessary pixel value ranges r2 and r3. Since the pixel values belonging to the main pixel value range r1 are allocated to the pixel values belonging to the target pixel value range t in the synthesized perspective image, for example, 1,000 pixel values in the second intermediate image P2 are synthesized. It is 50 in a fluoroscopic image. That is, all of the 2,000 pixel values in the second intermediate image P2 are converted to 50 by the conversion unit 18. Similarly, the pixel value of 1,000 in the second intermediate image P2 is the pixel value of 4,050 in the synthesized perspective image.
  • the pixel values from 1,001 to 1,999 in the second intermediate image P2 are linearly scaled and allocated within the target pixel value range t. That is, since the width of the target pixel value range t is four times the width of the main pixel value range r1, when the main pixel value range r1 is expanded four times, it matches the target pixel value range.
  • the pixel value of 1,000 in the second intermediate image P2 is set to 50 pixel value, but the pixel value of 1,001 in the second intermediate image P2 increased by 1 from the pixel value of 1,000 is 50. That is, the pixel value is increased by 4 from the pixel value of 54. Then, all the pixel values of 1,001 in the second intermediate image P2 are converted to 54 by the conversion unit 18.
  • the upper limit value of the main pixel value range r1 is assigned to the upper limit value of the target pixel value range t
  • the lower limit value of the main pixel value range r1 is assigned to the lower limit value of the target pixel value range t.
  • the pixel values in the unnecessary pixel value ranges r2 and r3 are also converted by the conversion unit 18, but the specific mode is not particularly limited.
  • all the pixel values in the unnecessary pixel value range r2 may be set to a predetermined pixel value (for example, 0) smaller than the lower limit value B1.
  • all the pixel values in the unnecessary pixel value range r3 may be set to predetermined pixel values (for example, 4,095) that are larger than the upper limit value B2.
  • the pixel values in the second intermediate image P2 are converted one after another by the conversion unit 18 to generate a composite perspective image.
  • This synthetic fluoroscopic image is displayed on the display unit 22, and the generation of the synthetic fluoroscopic image using the X-ray imaging apparatus 1 according to the first embodiment is finished.
  • the X-ray imaging apparatus 1 can generate a moving image by combining a plurality of synthetic fluoroscopic images over time. That is, the high output pulse and the low output pulse are alternately irradiated toward the subject M. This moving image is also displayed on the display unit 22. At this time, since luminance correction is performed again on each of the combined fluoroscopic images, this will be described.
  • the combining unit 19 calculates the brightness of each composite perspective image sent from the conversion unit 18. As a specific method, the pixel values of the pixels constituting the composite perspective image are averaged, and this value is used as the brightness.
  • the synthesized fluoroscopic images are sequentially sent to the connecting unit 19 over time, and the lightness is determined each time.
  • the combining unit 19 changes the pixel value of the newer synthetic fluoroscopic image to adjust the lightness.
  • a predetermined bias value is added to or subtracted from each pixel value of the composite perspective image to bring the brightness in the new composite perspective image closer to the brightness in the old composite perspective image.
  • the above-described configuration is not necessarily employed for the generation of a moving image in the X-ray imaging apparatus 1 of the first embodiment. That is, the combined fluoroscopic images may be simply combined without generating luminance correction again, and a moving image may be generated. Alternatively, the histogram S of the second intermediate image P2 acquired first is shared, and the subsequent synthetic fluoroscopic images are used. It is good also as a structure which produces
  • the configuration of the first embodiment can acquire a composite perspective image (subtraction image).
  • the second intermediate image P2 is subjected to brightness correction before being used for diagnosis. Since this brightness correction is performed every time the second intermediate image P2 is generated, even if the second output image P2 is regenerated by changing the high output coefficient and the low output coefficient, the diagnosis is performed.
  • the brightness adjustment of the synthetic perspective image to be used has already been made. Therefore, the operator does not need to adjust the brightness of the display unit 22 every time a composite perspective image is taken. Therefore, the X-ray imaging apparatus 1 that can be easily operated by the operator can be provided.
  • the configuration of the first embodiment generates a histogram S of the second intermediate image P2, and divides the entire pixel value range r0 of the histogram S into a main pixel value range r1 and unnecessary pixel value ranges r2 and r3. Then, each of the pixel values a of the pixels constituting the second intermediate image P2 is successively converted into new pixel values, thereby generating a composite perspective image. According to the configuration of the first embodiment, since the histogram S can be generated for each second intermediate image P2, no matter how many times the generation of the synthesized perspective image is repeated, the luminance correction using the histogram S is performed every time. Will be ensured.
  • the pixel value a belonging to the main pixel value range is changed to the pixel value a belonging to the target pixel value range t in the composite perspective image. Allocated.
  • the main pixel value range r1 where frequent pixel values are concentrated among the pixel values a in the second intermediate image P2 to the target pixel value range t. If the target pixel value range t is wider than the main pixel value range r1, it is possible to widen the difference in pixel values in pixels having similar pixel values. Therefore, the contrast of the synthesized fluoroscopic image is high, and a clear fluoroscopic image suitable for diagnosis can be acquired.
  • the present invention is not limited to the configuration of the above-described embodiment, and can be modified as follows.
  • the second intermediate image P2 is sent to the histogram generation unit 16, but the present invention is not limited to this configuration.
  • a configuration without a histogram generation unit may be employed. That is, although the high-output coefficient and the low-output coefficient are associated with the coefficient storage unit 12, the luminance correspondence table T1 may be associated with the coefficient.
  • the luminance correspondence table T1 is stored in the coefficient storage unit 12.
  • the luminance correspondence table T1 is a table in which the pixel values of the second intermediate image P2 correspond to the pixel values of the composite perspective image.
  • the luminance correction unit 20 reads the pixel values of the pixels constituting the second intermediate image P2, and searches for the same pixel value in the column of the second intermediate image P2 in the luminance correspondence table T1. And the pixel value of the synthetic
  • the plurality of luminance correspondence tables T1 are associated with the high output coefficient and the low output coefficient corresponding to them, a suitable luminance correspondence table T1 according to the change of the high output coefficient and the low output coefficient. Can be used, and the luminance distribution of the combined fluoroscopic image can be made constant regardless of the change of each coefficient. Note that the luminance correction unit in the present modification corresponds to the luminance adjusting means of the present invention.
  • the histogram S is generated using the entire area of the second intermediate image P2, but the present invention is not limited to this. That is, the histogram S may be generated using a part of the second intermediate image P2. By doing so, the time for generating the histogram S can be shortened, and the X-ray imaging apparatus 1 with improved image processing speed can be provided. Furthermore, by allowing the console 23 to select a region to be used for generating the histogram S in the second intermediate image P2, the luminance correction of the second intermediate image P2 based on the region of interest of the subject M can be performed. It becomes possible.
  • the present invention is suitable for medical devices.

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Medical Informatics (AREA)
  • Engineering & Computer Science (AREA)
  • Radiology & Medical Imaging (AREA)
  • Biomedical Technology (AREA)
  • Biophysics (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Optics & Photonics (AREA)
  • Pathology (AREA)
  • Physics & Mathematics (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • General Health & Medical Sciences (AREA)
  • Public Health (AREA)
  • Veterinary Medicine (AREA)
  • Apparatus For Radiation Diagnosis (AREA)

Abstract

 本発明は、FPD4が出力した検出データをX線透視画像に変換するX線透視画像生成部11と、X線透視画像に係数を乗算して、第1中間画像を生成する重み付け部14と、第1中間画像から第1中間画像を減算することで、第2中間画像P2を生成する減算処理部15と、第2中間画像の輝度調節を行って、合成透視画像を生成するヒストグラム生成部16と、分割部17と、変換部18とを備えている。これにより、第2中間画像は、生成されるたびに輝度調節がされるので、X線透視画像に乗算される係数を変更して、新たに第2中間画像を生成したとしても、表示される合成透視画像は、診断に好適なものとなっている。

Description

放射線撮影装置
 本発明は、二種類の放射線ビームを被検体に照射して、二種類の放射線透視像を生成して、それらに所定の重み付けを行いつつ、一方の放射線透視画像から他方の放射線透視画像を減算して、診断の目的に合わせた合成透視画像を生成するサブトラクション撮影を行うことができる放射線撮影装置に関する。
 放射線撮影装置は、被検体に放射線ビームを照射し、それを透過した放射線を検出して、放射線透視画像を生成する。しかしながら、例えば、被検体の骨部を診断する場合などにおいては、単に放射線透視画像を撮影するだけでは、被検体の軟組織も写りこんでしまい、詳細な診断ができない場合がある。
 放射線には、骨部、軟組織を透過しやすい波長のものと、そうでない波長のものがある。つまり、骨部と軟組織で放射線の透過しやすさが波長で異なっている。また、放射線ビームの出力を変化させると、放射線ビームに含まれる波長スペクトルが変化する。これを利用して、被検体の骨部を強調した合成透視画像(サブトラクション画像)を合成することが可能である(例えば、特許文献1参照)。
 まず、高出力の放射線ビームを被検体に照射する。このとき生成された放射線透視画像を高出力画像とする。次に、低出力の放射線ビームを被検体に照射する。このとき生成された放射線透視画像を低出力画像とする。高出力画像、および低出力画像の画素値に所定の係数を乗算し、高出力画像から低出力画像を減算すると、たとえば、被検体の骨部の造影に適した波長の放射線が強調され、被検体の骨部が鮮明に写りこんだ合成透視画像が得られる。また、高出力画像、および低出力画像の画素値に他の係数を乗算し、上述のような減算を行うと、たとえば、被検体の軟組織の造影に適した波長の放射線が強調され、被検体の軟組織が鮮明に写りこんだ合成透視画像が得られる。合成透視画像は、放射線撮影装置に付属のモニタに表示される。
 この様に、従来の構成では、診断の目的に合わせて、係数を選択し、これを高出力画像、および低出力画像に乗算することで、診断の目的に合致した合成透視画像が生成できるようになっている。
特開2007-236766号公報
 しかしながら、このような構成を有する従来例には、次のような問題がある。
 すなわち、従来の構成によれば、第1の合成透視画像が表示された後で、係数を再選択し、第2の合成透視画像を生成すると、直ちには第2の合成透視画像による診断を開始できない。仮に、モニタに第1の合成透視画像が映っているものとする。モニタが表示する放射線透視画像の輝度の調節は、例えば、オペレータがモニタに付属のツマミを操作することによって、第1の合成透視画像の視認に好適なものとされている。この状態で第2の合成透視画像がモニタに表示されると、第2の合成透視画像は、極端に暗い、または、極端に明るいものとなり、オペレータは、再び放射線透視画像の輝度の調節しなければならない。
 つまり、第1の合成透視画像の輝度分布と第2の合成透視画像の輝度分布とが同一であることが保証されていないので、オペレータは合成透視画像がモニタに表示される度に合成透視画像の輝度調節を行う必要があるのである。しかも、高出力画像、および低出力画像に乗算される係数を変更すれば、合成透視画像の輝度分布も自ずと変更してしまう。係数が変更させると、それに応じて、高出力画像から減算される画素値が変動するからである。
 この様に、従来の構成によれば、合成透視画像を生成する度に輝度の調節を行わなければならず、オペレータの操作が煩雑となるのである。
 本発明は、この様な事情に鑑みてなされたものであって、その目的は、合成透視画像(サブトラクション画像)を更新する度にモニタの輝度調節を必要とせず、操作が簡単な放射線撮影装置を提供することにある。
 本発明は、この様な目的を達成するために、次に様な構成をとる。
 すなわち、本発明は、所定の出力で放射線ビームを照射する放射線源と、放射線ビームを検出する放射線検出手段と、放射線源を制御する放射線源制御手段とを備えた放射線撮影装置において、放射線検出手段が出力した検出データを放射線透視画像に変換する放射線透視画像生成手段と、第1係数、および第2係数を記憶する係数記憶手段と、一方の放射線透視画像に第1係数を乗算して一方の第1中間画像を生成するとともに、他方の放射線透視画像に第2係数を乗算して他方の第1中間画像を生成する第1中間画像生成手段と、一方の第1中間画像から他方の第1中間画像を減算することで、第2中間画像を生成する第2中間画像生成手段と、第2中間画像の輝度調節を行って、合成透視画像を生成する輝度調節手段とを備えることを特徴とするものである。
 [作用・効果]本発明の構成は、放射線検出手段が出力した検出データを放射線透視画像に変換する放射線透視画像生成手段と、第1係数、および第2係数を記憶する係数記憶手段と、一方の放射線透視画像に第1係数を乗算するとともに一方の第1中間画像を生成するとともに、他方の放射線透視画像に第2係数を乗算して、他方の第1中間画像を生成する第1中間画像生成手段と、一方の第1中間画像から他方の第1中間画像を減算することで、第2中間画像を生成する第2中間画像生成手段とを備えている。これにより、一方の放射線透視画像から他方の放射線透視画像を減算することになり、合成透視画像(サブトラクション画像)が取得できる。そして、第2中間画像は、診断に利用される前に、輝度補正が行われる。この輝度補正は、第2中間画像が生成される度に行われるので、例え、第1係数、第2係数を変更させて第2中間画像を再生成したとしても、診断に利用される合成透視画像の輝度調節は、既になされている。したがって、オペレータは、合成透視画像を撮影するたびにモニタの輝度を調節する必要がない。したがってオペレータの操作が容易な放射線撮影装置が提供できることとなる。
 また、上述の輝度調節手段は、第2中間画像のヒストグラムを生成するヒストグラム生成手段と、ヒストグラムの全画素値範囲を主要画素値範囲と不要画素値範囲とに分割する分割手段と、第2中間画像を構成する画素の画素値の各々を新たな画素値に次々と変換していくことで合成透視画像を生成する変換手段とを備え、変換手段により、第2中間画像における画素値のうち、主要画素値範囲内に属する画素値は、合成透視画像における目標画素値範囲内に属する画素値に割り振られる構成とすることができる。
 [作用・効果]上記構成は、本発明における輝度調節の具体的な態様の1つを示している。すなわち、上記構成は、第2中間画像のヒストグラムを生成し、ヒストグラムの全画素値範囲を主要画素値範囲と不要画素値範囲とに分割する。そして、第2中間画像を構成する画素の画素値の各々を新たな画素値に次々と変換していくことで合成透視画像を生成するのである。上記構成によれば、第2中間画像毎にヒストグラムを生成することができるので、合成透視画像の生成を幾度繰返したとしても、その度ごとにヒストグラムを用いた輝度補正が確実にされることになる。
 また、上記構成によれば、第2中間画像における画素値のうち、主要画素値範囲内に属する画素値は、合成透視画像における目標画素値範囲内に属する画素値に割り振られる。この様な構成とすることで、第2中間画像における画素値のうち、頻出の画素値が集中している主要画素値範囲を目標画素値範囲に移行させることができる。目標画素値範囲を主要画素値範囲よりも広いものとすれば、強度的に互いに隣接する画素値の強度の差を広げることができる。したがって、合成透視画像のコントラストは、高いものとなり、鮮明で診断に好適な合成透視画像が取得できるのである。
 また、上述のヒストグラム生成手段は、第2中間画像の一部の領域を構成する画素を基にヒストグラムを生成する構成とすることもできる。
 [作用・効果]上記構成によれば、ヒストグラムを生成する時間を短縮でき、画像の処理速度が向上した放射線撮影装置が提供できる。さらに、第2中間画像におけるヒストグラムの生成に使用する領域を選択できるようにすることで、被検体の関心部位を基準とした第2中間画像の輝度補正が可能となる。
 また、上述の係数記憶手段は、さらに第2中間画像における画素値と合成透視画像における画素値が対応した画素値対応表を第1係数、および第2係数と連関させて記憶し、輝度調節手段は、第1係数、および第2係数に連関した輝度対応表を取得するとともに、輝度対応表を参照して、第2中間画像を構成する画素の画素値の各々を次々と新たな画素値としていくことで合成透視画像を生成する構成とすることもできる。
 [作用・効果]上記構成は、本発明における輝度調節の具体的な態様の1つを示している。すなわち、第2中間画像を構成する画素の画素値の各々を次々と新たな画素値としていくことで合成透視画像が生成されるが、このとき輝度対応表を利用する。こうすることで、ヒストグラムを生成する構成を省略でき、画像の処理速度が向上した放射線撮影装置が提供できる。
 また、本発明は、複数の合成透視画像を結合して動画を生成する結合手段を更に備えることもできる。
 [作用・効果]上記構成によれば、診断に好適な動画を生成することができる。上記構成によれば、例えば、第2中間画像を取得する毎にヒストグラムを生成して、合成透視画像を生成する構成とすることもできる。また、経時的に隣接した合成透視画像の明度があまりにも違いすぎる場合、新しい合成透視画像の輝度を再び調節する構成とすることもできる。さらには、始めに取得された第2中間画像のヒストグラムを共用のものとして、以降の合成透視画像を生成し、これを結合する構成としてもよい。いずれの構成によっても、診断に好適な動画が取得できる。
 本発明の構成は、合成透視画像(サブトラクション画像)が取得できる。そして、第2中間画像は、診断に利用される前に、輝度補正が行われる。この輝度補正は、第2中間画像が生成される度に行われるので、例え、設定を変更させて第2中間画像を再生成したとしても、診断に利用される合成透視画像の輝度調節は、既になされている。したがって、オペレータは、合成透視画像を撮影するたびにモニタの輝度を調節する必要がない。したがってオペレータの操作が容易な放射線撮影装置が提供できることとなる。
実施例1に係るX線撮影装置の構成を説明する機能ブロック図である。 実施例1に係るX線撮影装置の動作を説明するフローチャートである。 実施例1に係るヒストグラムを説明する模式図である。 実施例1に係る変換部の動作を説明する模式図である。 本発明の1変形例に係るX線撮影装置の構成を説明する機能ブロック図である。 本発明の1変形例に係るX線撮影装置の構成を説明する模式図である。
符号の説明
3  X線管(放射線源)
4  FPD(放射線検出手段)
6  X線管制御部(放射線源制御手段)
11 X線透視画像生成部(放射線透視画像生成手段)
12 係数記憶部(係数記憶手段)
14 重み付け部(第1中間画像生成手段)
15 減算処理部(第2中間画像生成手段)
16 ヒストグラム生成部(ヒストグラム生成手段)
17 分割部(分割手段)
18 変換部(変換手段)
19 結合部(結合手段)
 以下、本発明に係る放射線撮影装置の実施例について図面を参照しながら説明する。なお、実施例においては、X線を用いたX線撮影装置1について説明する。
 まず、実施例1に係るX線撮影装置1の構成について説明する。図1は、実施例1に係るX線撮影装置1の構成を説明する機能ブロック図である。図1に示すように、実施例1に係るX線撮影装置1には、被検体Mを載置する天板2と、その天板2に隣接してパルス状のX線ビームを照射するX線管3と、被検体を透過したX線を検出するフラット・パネル・ディテクタ(FPD)4と、FPD4に入射する散乱X線を除去するX線グリッド5とが設けられている。また、実施例1の構成は、X線管3の管電圧、管電流やX線ビームのパルス幅を制御するX線管制御部6とを備えている。そして、X線撮影装置1は、FPD4から出力された検出データをX線透視画像に変換するX線透視画像生成部11と、後述の種々の係数を記憶する係数記憶部12と、X線透視画像にログ変換処理を行うログ変換処理部13と、X線透視画像に係数を乗算する重み付け部14と、一方のX線透視画像から他方のX線透視画像を減算する減算処理部15と、X線透視画像のヒストグラムSを生成するヒストグラム生成部16と、ヒストグラムSの画素値範囲を分割する分割部17と、画素値を変換して合成透視画像を生成する変換部18と、合成透視画像を結合して動画を生成する結合部19とを備えている。なお、X線管は、本発明の放射線源に相当し、FPDは、本発明の放射線検出手段に相当する。また、X線管制御部は、本発明の放射線源制御手段に相当し、X線透視画像生成部は、本発明の放射線透視画像生成手段に相当する。さらに、係数記憶部、重み付け部、および減算処理部の各々は、本発明の係数記憶手段、第1中間画像生成手段、第2中間画像生成手段のそれぞれに相当する。さらにまた、ヒストグラム生成部、分割部、変換部、および結合部の各々は、本発明のヒストグラム生成手段、分割手段、変換手段、結合手段のそれぞれに相当する。また、ヒストグラム生成手段、分割手段、およ変換手段は、本発明の輝度補正手段の具体的な態様に相当する。
 また、X線撮影装置1は、オペレータの指示を受け付ける操作卓23と、X線透視画像、または動画が表示される表示部22とを備えている。
 さらにまた、X線撮影装置1は、X線管制御部6,X線透視画像生成部11,係数記憶部12,ログ変換処理部13,重み付け部14,減算処理部15,ヒストグラム生成部16,分割部17,変換部18,および結合部19を統括的に制御する主制御部25を備えている。この主制御部25は、CPUによって構成され、種々のプログラムを実行することにより、各部を実現している。また、上述の各部は、それらを担当する演算装置に分割されて実行されてもよい。
 次に、実施例1に係るX線撮影装置1の動作について説明する。図2は、実施例1に係るX線撮影装置1の動作を説明するフローチャートである。図2に示すように、実施例1に係るX線撮影装置1にて合成透視画像を生成するには、被検体を天板に載置する載置ステップS1と、被検体に高出力の放射線ビームを照射して、高出力画像Hを生成する高出力画像生成ステップS2と、被検体に低出力の放射線ビームを照射して、低出力画像Lを生成する低出力画像生成ステップS3と、高出力画像H,および低出力画像Lをログ変換するログ変換ステップS4と、高出力画像H,および低出力画像Lの各々に所定の係数を乗算して第1中間画像P1を生成し、第1中間画像P1同士を減算して第2中間画像P2を生成する中間画像生成ステップS5と、第2中間画像P2のヒストグラムS解析を実行し、ヒストグラムSの全画素値範囲を分割する分割ステップS6と、第2中間画像P2のヒストグラムSに基づいて、第2中間画像P2を構成する画素の画素値を新たな画素値に変更して、合成透視画像を生成する画素値変更ステップS7とを備えている。以降、これらの各ステップについて、順を追って説明する。
 <載置ステップS1,高出力画像生成ステップS2,低出力画像生成ステップS3>
 まず、天板2に被検体Mが載置され、被検体MのX線透視撮影の準備を行う。そして、X線管制御部6は、高出力のX線ビームを照射するように、X線管3の管電圧、管電流、X線ビームの照射時間を決定し、それにしたがってX線管3は、高出力のX線ビームを照射する。その後、X線管制御部6は、低出力のX線ビームを照射するように、X線管3の管電圧、管電流、X線ビームの照射時間を決定し、それにしたがってX線管3は、低出力のX線ビームを照射する。なお、このときのX線ビームは、パルス状となっている。
 高出力のX線ビームが照射された時点で、FPD4は、被検体Mを透過したX線ビームを検出し、検出信号をX線透視画像生成部11に向けて出力する。X線透視画像生成部11では、検出信号を基にX線透視画像を生成する。高出力のX線ビームを被検体Mに照射して得られた検出信号を基に生成されたX線透視画像は、高出力画像Hであり、低出力のX線ビームを被検体Mに照射して得られた検出信号を基に生成されたX線透視画像は、低出力画像Lである。この高出力画像H,および低出力画像Lは、後段のログ変換処理部13に送出される。なお、高出力画像、および低出力画像のいずれかのうちの一方は、本発明の一方のX線透視画像に相当し、高出力画像、および低出力画像のいずれかのうちの他方は、本発明の他方のX線透視画像に相当する。
 高出力のX線ビームは、低い出力のX線ビームは、様々な波長のX線を含んでいる。高出力のX線ビームは、例えば、低い出力のX線ビームと比べて、より多くの特定波長のX線を含んでいることになる。このように、高出力のX線ビームと、低出力のX線ビームとが被検体Mを透過するときの透過特性は、互いに異なっている。したがって、高出力のX線ビームに含まれる特定波長のX線の透過性に基づいた透視画像を得たい場合には、高出力画像Hから低出力画像Lを減算すればよいことになる。
 検査の目的によって、強調したいX線の波長は変動する。そこで、高出力画像H,および低出力画像Lに所定の係数を乗算した後、高出力画像Hから低出力画像Lを減算すれば、強調したいX線の波長を選択することができるようになる。そこで、今度は、この様な減算処理を行う。
 <ログ変換ステップS4>
 上述の減算処理に先立って、高出力画像H,および低出力画像Lは、ログ変換処理部13に送出される。そして、高出力画像Hは、ログ変換データlogHに変換され、低出力画像Lは、ログ変換データlogLに変換される。この処理の存在理由は、X線の波長と、X線の被検体の透過度には、指数関数的な関係が重畳していることによる。つまり、両画像H,Lを減算処理する前に、いったんログ変換処理を行っておけば、両ログ変換データlogH,logLには、予め強い傾向である指数関数的な関係が消去されている。したがって、ログ変換データ同士を減算処理すれば、高出力画像Hと低出力画像Lとの差異がより判然とするのである。
 <中間画像生成ステップS5>
 ログ変換処理によって生成された両ログ変換データlogH,logLは、重み付け部14に送出される。重み付け部14は、係数記憶部12で記憶された係数を読み出す。このとき、2つの係数がペアで読み出される。そのうちの一方がログ変換データlogHに乗算される高出力用係数であり、他方がログ変換データlogLに乗算される低出力用係数である。高出力用係数、および低出力用係数は、関連付けられて、係数のセットとなっている。係数記憶部12においては、この係数のセットを複数個記憶している。どの係数のセットが読み出されるかについては、検査の目的による。具体的には、操作卓23を通じてオペレータが検査目的を入力し、それにふさわしい係数のセットが重み付け部14に送出されるのである。
 そして、ログ変換データlogHには、高出力用係数が乗算され、ログ変換データlogLには、低出力用係数が乗算される。こうして、ログ変換データlogHは、第1中間画像H1に変換され、ログ変換データlogLは、第1中間画像L1に変換される。このとき、第1中間画像H1,L1は、未だログ変換が施された状態となっている。そして、第1中間画像H1と第1中間画像L1とが、減算処理部15に送出される。なお、低出力用係数、高出力用係数のうちの一方が本発明の第1係数に相当し、他方が本発明の第2係数に相当する。
 減算処理部15は、第1中間画像H1から第1中間画像L1を減算する減算処理が行われ、第2中間画像P2が生成される。この第2中間画像P2は、減算処理部15において、減算処理の後にログ逆変換処理も行われたものとなっている。
 <分割ステップS6>
 次に、本発明における特徴的な輝度補正の処理について説明する。すなわち、第2中間画像P2を構成する画素の画素値を変換して新たな画像を生成するのである。まず、第2中間画像P2は、ヒストグラム生成部16に送出され、第2中間画像P2のヒストグラムSを生成する。このヒストグラムSは、図3に示すように画素値aと、頻度fとが関連したグラフで表現することができる。画素値aは、0から4095までの値をとり、頻度fは、画素値aの各々が第2中間画像P2において、どの位の頻度で出現するかを表している。例えば、第2中間画像P2の全面が完全に黒く塗りつぶされている場合、頻度fは、画素値が0の領域のみに偏り、1から4095までの画素値の頻度fは0となる。実際は、第2中間画像P2には、被検体Mの透視画像が写りこんでいるので、頻度fは、図3に示すように、画素値aにおける全範囲(全画素値範囲r0)の内のある領域(主要画素値範囲r1)に偏っている。
 生成されたヒストグラムSは、分割部17に送出される。分割部17において、ヒストグラムSを参照することで、全画素値範囲r0は、主要画素値範囲r1と、不要画素値範囲r2,r3とに分割される。この実際の処理について説明する。分割部17は、画素値が0であるときの頻度fをヒストグラムSから読み取る。そして、画素値を1ずつ増加させながら、そのときの頻度fを読み取っていく。このとき読み取られた頻度fは、合算され、累積頻度Fとして表す。そして、累積頻度Fが、全画素値範囲r0における頻度fの合計のうちのどの位の割合を占めるかを求める。分割部17は、累積頻度Fの割合が所定の設定値となったところで、頻度fを読み取りを終了し、このときの画素値aを読み取る。このときの画素値をA1とする。分割部17においては、画素値が0からA1の間の領域を不要画素値範囲r2と判別する。設定値は、低く設定されるので、不要画素値範囲r2に属する画素値が第2中間画像に表れる頻度は、低いものとなっている。
 また、累積頻度Fは、高画素値側についても求められる。すなわち、画素値aが4095であるときの頻度fをヒストグラムSから読み取って、以降、画素値を1ずつ減少させながら累積頻度Fを求める。分割部17は、累積頻度Fの割合が所定の設定値となったところで、頻度fを読み取りを終了し、このときの画素値をA2とする。分割部17においては、画素値がA2から4095の間の領域を不要画素値範囲r3と判別する。この画素値A1,A2は、変換部18に送出される。これとは別に、減算処理部15から第2中間画像P2が変換部18に送出される。
 <画素値変更ステップS7>
 変換部18は、送出された画素値A1,A2からヒストグラムSにおける主要画素値範囲と不要画素値範囲を認識する。そして、この主要画素値範囲r1と不要画素値範囲r2,r3を基に第2中間画像P2を構成する画素の画素値を変換するのである。具体的には、第2中間画像P2における画素値のうち、主要画素値範囲r1内に属する画素値は、変換部18により、合成透視画像における目標画素値範囲t内に属する画素値に割り振られる。
 図4は、実施例1に係る変換部の動作を説明する模式図である。図4の上側の配列は、第2中間画像P2における全画素値範囲r0を表しており、画素値は0から4,095まで用意されている。また、画素値A1は、1,000であるものとし、画素値A2は、2,000であるものする。したがって、主要画素値範囲r1の幅は、1,000個分の画素値である。
 図4の下側の配列は、変換部18により、これから生成される合成透視画像における全画素値範囲を表しおり、画素値は、簡潔な説明のため、第2中間画像P2と同様に0から4,095まで用意されているものとする。しかし、実際は、この構成に限られない。合成透視画像における全画素値範囲には、所定の画素値範囲が設定されている。すなわち、下限値B1と上限値B2とに挟まれた目標画素値範囲tである。また、下限値B1は、50であるものとし、上限値B2は、4,050であるものとする。したがって、目標画素値範囲tの幅は、4,000個分の画素値である。
 変換部18は、第2中間画像P2を構成する各々の画素に対応する画像値を読み出し、これが主要画素値範囲r1にあるか、それとも不要画素値範囲r2,r3にあるかを判断する。主要画素値範囲r1内に属する画素値は、合成透視画像における目標画素値範囲t内に属する画素値に割り振られるのであるから、例えば、第2中間画像P2における1,000の画素値は、合成透視画像における50とされる。すなわち、第2中間画像P2における2,000の画素値は、変換部18によって全て50に変換される。同様に、第2中間画像P2における1,000の画素値は、合成透視画像における4,050の画素値とされる。
 第2中間画像P2における1,001から1,999までの画素値は、線形的なスケーリングが施されて、目標画素値範囲t内に割り振られる。すなわち、目標画素値範囲tの幅は、主要画素値範囲r1の幅の4倍となっているわけであるから、主要画素値範囲r1を4倍に拡大すると目標画素値範囲に一致する。第2中間画像P2における1,000の画素値は、50の画素値とされたが、1,000の画素値よりも1だけ増加した第2中間画像P2における1,001の画素値は、50の画素値よりも4だけ増加した54の画素値とされるわけである。そして、第2中間画像P2における1,001の画素値は、変換部18によって全て54に変換される。つまり、主要画素値範囲r1の上限値は、目標画素値範囲tの上限値に割り振られ、主要画素値範囲r1の下限値は、目標画素値範囲tの下限値に割り振られる。
 不要画素値範囲r2,r3における画素値も変換部18によって変換されるが、その具体的な態様は特に限定されない。たとえば、不要画素値範囲r2の画素値を全て下限値B1よりも小さい所定の画素値(例えば0)としてもよい。同様に、不要画素値範囲r3の画素値を全て上限値B2よりも大きい所定の画素値(例えば4,095)としてもよい。
 この様に、第2中間画像P2における画素値は次々と変換部18によって変換され、合成透視画像が生成される。この合成透視画像は、表示部22に表示され、実施例1に係るX線撮影装置1を用いた合成透視画像の生成は終了となる。
 なお、実施例1に係るX線撮影装置1は複数の合成透視画像を経時的に結合して動画を生成することができる。つまり、高出力のパルスと低出力のパルスとが交互に被検体Mに向けて照射されることになる。この動画も表示部22にて表示される。このとき、合成透視画像の各々について、再び輝度補正が施されるので、これについて説明する。結合部19は、変換部18から送出された合成透視画像の各々の明度を割り出す。この具体的な方法としては、合成透視画像を構成する画素の画素値を平均してこの値を明度とする。合成透視画像は、経時的に次々と結合部19に送出され、その度ごとに明度が割り出されるのである。経時的に隣接した合成透視画像の明度があまりにも違いすぎると、動画の視認性の悪化に繋がる。本発明によれば、経時的に隣接した合成透視画像の明度が5%以上相違している場合、結合部19は、新しい方の合成透視画像の画素値を変更して明度を調節する。明度の調節方法としては、合成透視画像の画素値の各々に対して所定のバイアス値を加算・減算することで、新しい方の合成透視画像における明度を古い方の合成透視画像における明度に近づける。
 なお、実施例1のX線撮影装置1における動画の生成については、必ずしも上述の構成を採用しなくともよい。すなわち、輝度補正を再び行わずに合成透視画像を単に結合して動画を生成してもよいし、始めに取得された第2中間画像P2のヒストグラムSを共用のものとして、以降の合成透視画像を生成し、これを結合する構成としてもよい。つまり、X線撮影装置1に搭載されるCPUの処理能力と、診断の種別に応じて、上述のような種々の構成が選択可能である。
 以上のように、実施例1の構成は、合成透視画像(サブトラクション画像)が取得できる。そして、第2中間画像P2は、診断に利用される前に、輝度補正が行われる。この輝度補正は、第2中間画像P2が生成される度に行われるので、例え、高出力用係数、および低出力用係数を変更させて第2中間画像P2を再生成したとしても、診断に利用される合成透視画像の輝度調節は、既になされている。したがって、オペレータは、合成透視画像を撮影するたびに表示部22の輝度を調節する必要がない。したがってオペレータの操作が容易なX線撮影装置1が提供できることとなる。
 実施例1の構成は、第2中間画像P2のヒストグラムSを生成し、ヒストグラムSの全画素値範囲r0を主要画素値範囲r1と不要画素値範囲r2,r3とに分割する。そして、第2中間画像P2を構成する画素の画素値aの各々を新たな画素値に次々と変換していくことで合成透視画像を生成するのである。実施例1の構成によれば、第2中間画像P2毎にヒストグラムSを生成することができるので、合成透視画像の生成を幾度繰返したとしても、その度ごとにヒストグラムSを用いた輝度補正が確実にされることになる。
 また、実施例1の構成によれば、第2中間画像P2における画素値のうち、主要画素値範囲内に属する画素値aは、合成透視画像における目標画素値範囲t内に属する画素値aに割り振られる。この様な構成とすることで、第2中間画像P2における画素値aのうち、頻出の画素値が集中している主要画素値範囲r1を目標画素値範囲tに移行させることができる。目標画素値範囲tを主要画素値範囲r1よりも広いものとすれば、画素値が互いに似通った画素において、画素値の差を広げることができる。したがって、合成透視画像のコントラストは、高いものとなり、鮮明で診断に好適な合成透視画像が取得できるのである。
 本発明は、上述の実施例の構成に限られず、下記のような変形実施が可能である。
 (1)上述した実施例においては、第2中間画像P2は、ヒストグラム生成部16に送出されていたが、本発明は、この構成に限らない。図5に示すように、ヒストグラム生成部を有しない構成とすることもできる。すなわち、係数記憶部12には、高出力用係数と低出力用係数とが関連付けられていたが、これらに加えて輝度対応表T1をも連関させる構成としてもよい。輝度対応表T1は、係数記憶部12において記憶されている。この輝度対応表T1は、図6に示すように、第2中間画像P2の画素値と合成透視画像の画素値が対応した表となっている。輝度補正部20は、第2中間画像P2を構成する画素の画素値を読み出して、輝度対応表T1における第2中間画像P2の欄において、これと同一の画素値を検索する。そして、これに対応する合成透視画像の画素値を取得するわけである。最後に、第2中間画像P2における画素値を取得された合成透視画像の画素値に変換する。こうすることで、ヒストグラムSを生成する構成を省略でき、画像の処理速度が向上したX線撮影装置1が提供できる。しかも、複数の輝度対応表T1がそれらに対応する高出力用係数、および低出力用係数に関連付けられているので、高出力用係数と低出力用係数の変更に応じて好適な輝度対応表T1を使用することができ、各係数の変更に係らず合成透視画像の輝度分布を一定とすることができる。なお、本変形例における輝度補正部は、本発明の輝度調節手段に相当する。
 (2)上述した実施例において、第2中間画像P2の全域を用いてヒストグラムSを生成していたが、本発明はこれに限らない。すなわち、第2中間画像P2の一部を用いてヒストグラムSを生成する構成としてもよい。こうすることで、ヒストグラムSを生成する時間を短縮でき、画像の処理速度が向上したX線撮影装置1が提供できる。さらに、操作卓23にて、第2中間画像P2におけるヒストグラムSの生成に使用する領域を選択できるようにすることで、被検体Mの関心部位を基準とした第2中間画像P2の輝度補正が可能となる。
 (3)上述した実施例は、医用の装置であったが、この発明は、工業用や、原子力用の装置にも適応することができる。
 以上のように、本発明は、医療機器に適している。

Claims (5)

  1.  所定の出力で放射線ビームを照射する放射線源と、前記放射線ビームを検出する放射線検出手段と、前記放射線源を制御する放射線源制御手段とを備えた放射線撮影装置において、
     前記放射線検出手段が出力した検出データを放射線透視画像に変換する放射線透視画像生成手段と、
     第1係数、および第2係数を記憶する係数記憶手段と、
     一方の放射線透視画像に前記第1係数を乗算して一方の第1中間画像を生成するとともに、他方の放射線透視画像に前記第2係数を乗算して他方の第1中間画像を生成する第1中間画像生成手段と、
     前記一方の第1中間画像から前記他方の第1中間画像を減算することで、第2中間画像を生成する第2中間画像生成手段と、
     前記第2中間画像の輝度調節を行って、合成透視画像を生成する輝度調節手段とを備えることを特徴とする放射線撮影装置。
  2.  請求項1に記載の放射線撮影装置において、
     前記輝度調節手段は、
     前記第2中間画像のヒストグラムを生成するヒストグラム生成手段と、
     前記ヒストグラムの全画素値範囲を主要画素値範囲と不要画素値範囲とに分割する分割手段と、
     前記第2中間画像を構成する画素の画素値の各々を新たな画素値に次々と変換していくことで合成透視画像を生成する変換手段とを備え、
     前記変換手段により、前記第2中間画像における画素値のうち、前記主要画素値範囲内に属する画素値は、前記合成透視画像における目標画素値範囲内に属する画素値に割り振られることを特徴とする放射線撮影装置。
  3.  請求項2に記載の放射線撮影装置において、
     前記ヒストグラム生成手段は、前記第2中間画像の一部の領域を構成する画素を基に前記ヒストグラムを生成することを特徴とする放射線撮影装置。
  4.  請求項1に記載の放射線撮影装置において、
     前記係数記憶手段は、さらに前記第2中間画像における画素値と前記合成透視画像における画素値が対応した画素値対応表を前記第1係数、および前記第2係数と連関させて記憶し、
     前記輝度調節手段は、前記第1係数、および前記第2係数に連関した輝度対応表を取得するとともに、前記輝度対応表を参照して、前記第2中間画像を構成する画素の画素値の各々を次々と新たな画素値としていくことで前記合成透視画像を生成することを特徴とする放射線撮影装置。
  5.  請求項1ないし請求項4のいずれかに記載の放射線撮影装置において、
     複数の合成透視画像を結合して動画を生成する結合手段を更に備えていることを特徴とする放射線撮影装置。
PCT/JP2008/069781 2008-10-30 2008-10-30 放射線撮影装置 WO2010050031A1 (ja)

Priority Applications (1)

Application Number Priority Date Filing Date Title
PCT/JP2008/069781 WO2010050031A1 (ja) 2008-10-30 2008-10-30 放射線撮影装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
PCT/JP2008/069781 WO2010050031A1 (ja) 2008-10-30 2008-10-30 放射線撮影装置

Publications (1)

Publication Number Publication Date
WO2010050031A1 true WO2010050031A1 (ja) 2010-05-06

Family

ID=42128409

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2008/069781 WO2010050031A1 (ja) 2008-10-30 2008-10-30 放射線撮影装置

Country Status (1)

Country Link
WO (1) WO2010050031A1 (ja)

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS60250792A (ja) * 1984-05-28 1985-12-11 Fuji Photo Film Co Ltd サブトラクシヨン画像の階調処理最適化方法および装置
JPS6331641A (ja) * 1986-07-24 1988-02-10 富士写真フイルム株式会社 放射線画像情報処理方法
JP2004023408A (ja) * 2002-06-17 2004-01-22 Hitachi Medical Corp X線画像診断装置
JP2004321310A (ja) * 2003-04-22 2004-11-18 Canon Inc 放射線撮像装置
JP2005006918A (ja) * 2003-06-19 2005-01-13 Hitachi Medical Corp 放射線画像撮像装置

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS60250792A (ja) * 1984-05-28 1985-12-11 Fuji Photo Film Co Ltd サブトラクシヨン画像の階調処理最適化方法および装置
JPS6331641A (ja) * 1986-07-24 1988-02-10 富士写真フイルム株式会社 放射線画像情報処理方法
JP2004023408A (ja) * 2002-06-17 2004-01-22 Hitachi Medical Corp X線画像診断装置
JP2004321310A (ja) * 2003-04-22 2004-11-18 Canon Inc 放射線撮像装置
JP2005006918A (ja) * 2003-06-19 2005-01-13 Hitachi Medical Corp 放射線画像撮像装置

Similar Documents

Publication Publication Date Title
JP5669828B2 (ja) 画像処理装置および蛍光観察装置
JP2008194374A (ja) 放射線画像投影装置および方法
JP2006312040A (ja) 医用イメージング装置のための露出を制御する方法及びシステム
JP2011245117A (ja) X線画像診断装置
JP3667317B2 (ja) 放射線断層撮影装置
JP5422171B2 (ja) X線画像診断装置
US7760856B2 (en) Radiographic apparatus and radiation detection signal processing method
JP2005296332A (ja) X線診断装置、画像生成装置及び画像生成方法
JP5244164B2 (ja) 内視鏡装置
JPWO2013118193A1 (ja) 放射線撮影装置
JP2015024068A (ja) 医用画像処理装置
US7469034B2 (en) Method for analyzing and representing x-ray projection images and x-ray examination unit
JP2011019712A (ja) X線透視撮影装置
US20080075379A1 (en) Image processing device and image processing method
US8121372B2 (en) Method for reducing image noise in the context of capturing an image using two different radiation spectra
JP2014064608A (ja) 放射線撮影装置、その制御方法及びプログラム
WO2010050031A1 (ja) 放射線撮影装置
JP2011239804A (ja) X線検査装置
JP2010075555A (ja) X線画像診断装置
US9959614B2 (en) Image processing device and method
JP2005034539A (ja) 骨密度分布測定機能付きx線画像診断装置
US7949174B2 (en) System and method for calibrating an X-ray detector
US10159457B2 (en) X-ray diagnostic apparatus
JP2007244484A (ja) X線診断装置
JP5044457B2 (ja) 放射線画像撮影装置及び撮影条件設定方法

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 08877742

Country of ref document: EP

Kind code of ref document: A1

NENP Non-entry into the national phase

Ref country code: DE

122 Ep: pct application non-entry in european phase

Ref document number: 08877742

Country of ref document: EP

Kind code of ref document: A1

NENP Non-entry into the national phase

Ref country code: JP