KR101092207B1 - Method and Apparatus for Correcting Balance of Digital Radiation Detection Image having Nonlinear Response Characteristics - Google Patents

Method and Apparatus for Correcting Balance of Digital Radiation Detection Image having Nonlinear Response Characteristics Download PDF

Info

Publication number
KR101092207B1
KR101092207B1 KR1020090116235A KR20090116235A KR101092207B1 KR 101092207 B1 KR101092207 B1 KR 101092207B1 KR 1020090116235 A KR1020090116235 A KR 1020090116235A KR 20090116235 A KR20090116235 A KR 20090116235A KR 101092207 B1 KR101092207 B1 KR 101092207B1
Authority
KR
South Korea
Prior art keywords
correction
image
section
luminance
value
Prior art date
Application number
KR1020090116235A
Other languages
Korean (ko)
Other versions
KR20110059486A (en
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 KR1020090116235A priority Critical patent/KR101092207B1/en
Publication of KR20110059486A publication Critical patent/KR20110059486A/en
Application granted granted Critical
Publication of KR101092207B1 publication Critical patent/KR101092207B1/en

Links

Images

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
    • A61B6/52Devices using data or image processing specially adapted for radiation diagnosis
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00Image enhancement or restoration
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10116X-ray image
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/20Special algorithmic details
    • G06T2207/20048Transform domain processing

Abstract

본 발명은 비선형 응답특성을 가지는 디지털 방사선 검출 이미지의 각 화소 마다 휘도 특성이 평탄화되도록 이미지 밸런스를 보정하기 위하여, 이미지 센서의 응답 특성에 따라 구간별 응답 특성에 맞게 비선형 보정 방식과 선형보정 방식을 적절히 혼합 사용함으로써, 세분화 구간의 수를 줄이면서 보정 복잡도를 줄이며 정밀한 밸런스 보정이 이루어질 수 있는 이미지의 밸런스 보정 방법 및 장치에 관한 것이다. 본 발명의 일면에 따른, 이미지의 밸런스를 위한 보정 방법은, 입력 휘도값별 출력 휘도값을 나타내는 이미지 센서의 응답 특성 데이터에 대한 2차 기울기 정보를 임계치 상하로 구분되는 복수의 구간으로 나누는 단계, 및 상기 이미지 센서로부터 수신한 화소 데이터의 상기 복수의 구간 중 해당 구간에 대응된 보정 방식으로 상기 화소 데이터를 보정하는 단계를 포함한다.According to the present invention, in order to correct an image balance such that luminance characteristics are flattened for each pixel of a digital radiation detection image having a nonlinear response characteristic, a nonlinear correction method and a linear correction method are appropriately adapted to the response characteristics of each section according to the response characteristics of the image sensor. By using a mixture, the present invention relates to a method and apparatus for correcting a balance of an image in which the number of subdivision sections is reduced while reducing the complexity of correction and precise balance correction is possible. According to an aspect of the present invention, there is provided a correction method for balancing an image, comprising: dividing second slope information of response characteristic data of an image sensor indicating an output luminance value for each input luminance value into a plurality of sections divided by upper and lower threshold values, and Correcting the pixel data in a correction scheme corresponding to a corresponding section among the plurality of sections of pixel data received from the image sensor.

방사선 이미지, 비선형 응답 특성, 평탄화, 1차 기울기, 2차 기울기 Radiographic image, nonlinear response characteristics, planarization, first-order slope, second-order slope

Description

비선형 응답특성을 가지는 디지털 방사선 검출 이미지의 밸런스 보정 방법 및 장치{Method and Apparatus for Correcting Balance of Digital Radiation Detection Image having Nonlinear Response Characteristics}Method and Apparatus for Correcting Balance of Digital Radiation Detection Image with Non-Linear Response Characteristics {Method and Apparatus for Correcting Balance of Digital Radiation Detection Image having Nonlinear Response Characteristics}

본 발명은 이미지의 밸런스 보정 방법 및 장치에 관한 것으로서, 비선형 응답특성을 가지는 디지털 방사선 검출 이미지의 각 화소 마다 휘도 특성이 평탄화되도록 이미지 밸런스를 보정하기 위하여, 이미지 센서의 응답 특성에 따라 구간별 응답 특성에 맞게 비선형 보정 방식과 선형보정 방식을 적절히 혼합 사용함으로써, 세분화 구간의 수를 줄이면서 보정 복잡도를 줄이며 정밀한 밸런스 보정이 이루어질 수 있는 이미지의 밸런스 보정 방법 및 장치에 관한 것이다.The present invention relates to a method and apparatus for correcting an image balance, and to correcting an image balance such that luminance is flattened for each pixel of a digital radiation detection image having a nonlinear response characteristic, the response characteristic of each section according to the response characteristic of the image sensor The present invention relates to a method and apparatus for correcting a balance of an image in which a nonlinear correction method and a linear correction method are suitably used, thereby reducing the complexity of the correction while reducing the number of subdivision sections and enabling precise balance correction.

일반적으로 환자의 암이나 기타 환부의 상태를 진단하기 위하여 X-레이 등 방사선을 환부에 조사하고 해당 영상을 획득하는 기술을 이용한다. X-레이 촬영 기술에는 일반적인 아날로그 필름에 정지 영상을 획득하는 GR(General Radiography) 방식과 이미지 판에 조사되는 X-레이를 소정 리더기로 읽어 영상을 디지털화하는 CR(Computed Radiography) 방식 등이 널리 사용되고 있다. 이외에도, 최근에는 반도체 이미지 센서를 이용해 직접 X-레이를 감지하여 영상을 디지털화하는 DR(Digital Radiography) 방식도 등장하였다. In general, in order to diagnose a patient's cancer or other lesions, X-ray radiation is used to irradiate the lesion and acquire an image. X-ray imaging technology has been widely used, such as a general radiography (GR) method for acquiring still images on a general analog film, and a computed radiography (CR) method for digitizing an image by reading an X-ray irradiated onto an image plate with a predetermined reader. . In addition, recently, a digital radiography (DR) method of digitizing an image by directly detecting an X-ray using a semiconductor image sensor has been introduced.

DR 방식에 있어서 최근에 반도체 평판 이미지 센서를 이용한 다양한 방사선 진단 기기들이 개발 및 상품화되고 있다. 이와 같은 반도체 평판 이미지 센서를 이용한 방사선 진단 방식은 그 편리함으로 인해 유방암이나 기타 환부의 상태를 진단하기 위하여 널리 보급될 전망이다.Recently, in the DR method, various radiation diagnostic devices using semiconductor flat panel image sensors have been developed and commercialized. Radiation diagnostic method using such a semiconductor flat plate image sensor is expected to be widely used for diagnosing the condition of breast cancer or other lesions due to its convenience.

이러한 각종 진단 장비에 사용되고 있는 반도체 평판 이미지 센서는 수백만 개의 화소로 이루어져 있고, 센서의 특성상 모든 화소가 조금씩 휘도 특성이 다를 수 밖에 없어 이를 보정하기 위한 전처리 과정을 거치게 된다. The semiconductor flat panel image sensor used in such various diagnostic equipment is composed of millions of pixels, and all the pixels have little brightness characteristics due to the characteristics of the sensor.

전처리 과정중 각 화소마다 조금씩 다른 휘도 특성을 보정해 영상을 평탄화시키는 이미지 밸런스 보정 방법으로서, 완전히 선형응답특성을 가지는 이미지 센서의 경우 이미지 센서에 일정한 입력을 주고 획득한 화이트 이미지의 전체 휘도 평균값을 실제 측정된 이미지의 휘도값과 곱한후 각 화소의 화이트 이미지의 휘도값으로 나눠주는 방법이 이용될 수 있다.An image balance correction method that flattens an image by correcting slightly different luminance characteristics for each pixel during the preprocessing process.In the case of an image sensor having a completely linear response characteristic, the average luminance value of the white image obtained by giving a constant input to the image sensor is actually used. A method of multiplying by the luminance value of the measured image and dividing by the luminance value of the white image of each pixel may be used.

하지만, 비선형 응답특성을 가지는 이미지 센서의 경우에, 위와 같은 선형응답특성에 대한 처리 방법으로는 임의의 입력에 대해 밸런스 영상을 얻을 수 없다. 일반적으로 비선형 응답 특성을 가지는 이미지 센서의 휘도 특성의 밸런스 보정은, 입력구간을 세분화하고 각 세분화된 구간에서는 응답특성이 선형이라는 가정하에, 구간별로 위에서 언급한 일반적 선형 보정을 수행하는 방법을 사용하거나, 다차 곡선형태의 응답 특성을 가질 경우에 구간을 세분화하여 각 구간에서는 2차 곡선을 만든 후 2차 다항식에 의해 보정 절차를 진행하는 방식을 사용하게 된다. 하지만, 이와 같은 구간별 선형 보정 방식에서는 정확한 보정을 위하여 세부 구간의 수가 많아져 구간별로 보정된 이미지를 얻는 과정이 오래 걸리고 복잡해진다. 또한, 위와 같은 다항보정방식에 있어서는, 세분화 구간의 수가 적어 보정 이미지 얻는 과정을 줄일 수 있으나 보정 이미지의 적용식이 선형보정 방식에 비해 복잡해지는 단점이 있다.However, in the case of an image sensor having a nonlinear response characteristic, a balance image cannot be obtained for an arbitrary input by the above-described processing method for the linear response characteristic. In general, the balance correction of the luminance characteristics of an image sensor having a nonlinear response characteristic uses a method of performing the above-described general linear correction for each section on the assumption that the input section is segmented and the response characteristic is linear in each segmented section. In the case of having a response characteristic in the form of a multi-order curve, a section is subdivided to make a quadratic curve in each section, and then a correction procedure is performed by a quadratic polynomial. However, in such a linear correction scheme for each section, the number of detailed sections increases for accurate correction, and the process of obtaining a corrected image for each section takes a long time and is complicated. In addition, in the above polynomial correction method, the number of subdivision sections can be reduced to reduce the process of obtaining a corrected image, but the application formula of the corrected image is more complicated than the linear correction method.

본 발명은 상술한 문제점을 해결하기 위한 것으로서, 본 발명의 목적은, 비선형 응답특성을 가지는 디지털 방사선 검출 이미지의 각 화소 마다 휘도 특성이 평탄화되도록 이미지 밸런스를 보정하기 위하여, 이미지 센서의 응답 특성에 따라 구간별 응답 특성에 맞게 비선형 보정 방식과 선형보정 방식을 적절히 혼합 사용함으로써, 세분화 구간의 수를 줄이면서 보정 복잡도를 줄이며 정밀한 밸런스 보정이 이루어질 수 있는 이미지의 밸런스 보정 방법 및 장치를 제공하는 데 있다.SUMMARY OF THE INVENTION The present invention has been made to solve the above-described problems, and an object of the present invention is to adjust an image balance such that luminance is flattened for each pixel of a digital radiation detection image having a nonlinear response characteristic, according to a response characteristic of an image sensor. The present invention provides a method and apparatus for correcting a balance of an image in which a nonlinear correction method and a linear correction method are appropriately mixed according to the response characteristics of each section, thereby reducing the complexity of the correction while reducing the number of subdivision sections and enabling precise balance correction.

먼저, 본 발명의 특징을 요약하면, 상기의 목적을 달성하기 위한 본 발명의 일면에 따른 이미지의 밸런스를 위한 보정 방법은, 입력 휘도값별 출력 휘도값을 나타내는 상기 이미지 센서의 응답 특성 데이터에 대한 2차 기울기 정보를 임계치 상하로 구분되는 복수의 구간으로 나누는 단계; 및 상기 이미지 센서로부터 수신한 화소 데이터의 상기 복수의 구간 중 해당 구간에 대응된 보정 방식으로 상기 화소 데이터를 보정하는 단계를 포함한다.First, to summarize the features of the present invention, a method for correcting the balance of the image according to an aspect of the present invention for achieving the above object, 2 for the response characteristic data of the image sensor indicating the output luminance value for each input luminance value Dividing the difference slope information into a plurality of sections divided into upper and lower thresholds; And correcting the pixel data in a correction scheme corresponding to a corresponding section among the plurality of sections of pixel data received from the image sensor.

상기 화소 데이터는 방사선 진단 장비에서 방사선에 대한 영상 데이터를 획득하는 이미지 센서로부터 검출된 각 화소의 디지털 영상 신호일 수 있다.The pixel data may be a digital image signal of each pixel detected from an image sensor that acquires image data of radiation in a radiation diagnostic equipment.

상기 2차 기울기 정보는 상기 응답 특성 데이터에 대한 기울기 정보에 대하여 재차 기울기를 계산하고 절대값을 취한 정보일 수 있다.The second slope information may be information obtained by calculating a slope again and taking an absolute value with respect to the slope information of the response characteristic data.

상기 임계치는 0일 수 있다.The threshold may be zero.

상기 보정하는 단계는, 상기 임계치 이상의 구간에 대하여 비선형 보정을 수행하는 단계; 및 상기 임계치 보다 작은 구간에 대하여 선형 보정을 수행하는 단계를 포함한다.The correcting may include: performing nonlinear correction on a section above the threshold; And performing linear correction on the section smaller than the threshold.

상기 비선형 보정은, 해당 구간에 대하여 상기 응답 특성 데이터에 대하여 피팅한 2차 이상의 다차 함수에 기초하여 상기 화소 데이터를 보정하는 것이고, 상기 선형 보정은, 해당 구간에 대하여 상기 응답 특성 데이터에 대하여 피팅한 1차 함수에 기초하여 상기 화소 데이터를 보정하는 것을 포함한다.The nonlinear correction is to correct the pixel data based on a quadratic or higher quadratic function fitted to the response characteristic data for the interval, and the linear correction is to fit the response characteristic data for the interval. Correcting the pixel data based on a linear function.

상기 비선형 보정은, 상기 다차 함수에 의하여 결정된 해당 구간에 대한 전체 화소의 휘도 평균값을 상기 화소 데이터와 곱한 후, 곱해진 값을 해당 구간에서의 화소의 평균 휘도에 대응된 휘도값으로 나누는 보정이고, 상기 선형 보정은, 상기 1차 함수에 의하여 결정된 해당 구간에 전체 화소의 휘도 평균값을 상기 화소 데이터와 곱한 후, 곱해진 값을 해당 구간에서의 화소의 평균 휘도에 대응된 휘도값으로 나누는 보정인 것을 포함한다.The nonlinear correction is a correction that multiplies the luminance average value of all pixels for the corresponding section determined by the multi-order function with the pixel data, and then divides the multiplied value by the luminance value corresponding to the average luminance of the pixels in the corresponding section, The linear correction is a correction that multiplies the luminance average value of all pixels in the interval determined by the linear function with the pixel data, and then divides the multiplied value by the luminance value corresponding to the average luminance of the pixels in the interval. Include.

그리고, 본 발명의 다른 일면에 따른 이미지의 밸런스를 위한 보정 장치는, 입력 휘도값별 출력 휘도값을 나타내는 이미지 센서의 응답 특성 데이터에 대한 2 차 기울기 정보를 임계치 상하로 구분되는 복수의 구간으로 나누는 구간 설정부; 및 상기 이미지 센서로부터 수신하는 화소 데이터의 상기 복수의 구간 중 해당 구간에 대응된 보정 방식으로 상기 화소 데이터를 보정하는 보정부를 포함한다.According to another aspect of the present invention, an apparatus for correcting an image balance includes: dividing second slope information of response characteristic data of an image sensor indicating an output luminance value for each input luminance value into a plurality of sections divided by upper and lower thresholds; Setting unit; And a correction unit configured to correct the pixel data in a correction scheme corresponding to a corresponding section among the plurality of sections of pixel data received from the image sensor.

상기 보정 장치는, 상기 응답 특성 데이터에 대한 기울기 정보를 생성하는 제1 기울기 계산부; 및 상기 제1 기울기 계산부로부터의 기울기 정보에 대한 기울기 정보를 구하고 절대값을 계산하여 상기 2차 기울기 정보를 생성하는 제2 기울기 계산부를 더 포함한다.The correction device may include a first slope calculator configured to generate slope information of the response characteristic data; And a second slope calculator which obtains slope information on the slope information from the first slope calculator and calculates an absolute value to generate the second slope information.

상기 보정부는 상기 임계치 이상의 구간에 대하여 비선형 보정을 수행하고, 상기 임계치 보다 작은 구간에 대하여 선형 보정을 수행할 수 있다.The correction unit may perform nonlinear correction on a section above the threshold and perform linear correction on a section below the threshold.

상기 보정부는, 상기 비선형 보정으로서, 해당 구간에 대하여 상기 응답 특성 데이터에 대하여 피팅한 2차 이상의 다차 함수에 기초하여 상기 화소 데이터를 보정하고, 상기 선형 보정으로서, 해당 구간에 대하여 상기 응답 특성 데이터에 대하여 피팅한 1차 함수에 기초하여 상기 화소 데이터를 보정할 수 있다.The correction unit corrects the pixel data based on a quadratic or higher quadratic function fitted to the response characteristic data for the corresponding section as the nonlinear correction, and applies the linear characteristic to the response characteristic data for the section. The pixel data may be corrected based on a first-order function fitted with respect to each other.

상기 보정부는, 상기 비선형 보정에서, 상기 다차 함수에 의하여 결정된 해당 구간에 대한 전체 화소의 휘도 평균값을 상기 화소 데이터와 곱한 후, 곱해진 값을 해당 구간에서의 화소의 평균 휘도에 대응된 휘도값으로 나누고, 상기 선형 보정에서, 상기 1차 함수에 의하여 결정된 해당 구간에 전체 화소의 휘도 평균값을 상기 화소 데이터와 곱한 후, 곱해진 값을 해당 구간에서의 화소의 평균 휘도에 대응된 휘도값으로 나눌 수 있다.The correction unit, in the nonlinear correction, multiplies the luminance average value of all pixels for the corresponding section determined by the multi-order function with the pixel data, and then multiplies the multiplied value by the luminance value corresponding to the average luminance of the pixels in the corresponding section. In the linear correction, a luminance average value of all pixels may be multiplied by the pixel data in the corresponding interval determined by the first order function, and then the multiplied value may be divided by a luminance value corresponding to the average luminance of the pixels in the corresponding interval. have.

본 발명에 따른 이미지 밸런스 보정 방법 및 장치에 따르면, 기존의 이미지 보정 방법 대신 이미지 센서의 구간별 응답 특성을 고려한 적절한 보정방식을 사용함으로써 보정의 정확도는 높이면서도 보정 이미지 획득 과정과 적용의 복잡도를 낮춰, 결과적으로 방사선 검출 이미지의 최종적인 결과물의 화질과 신뢰성을 개선하는 효과를 얻을 수 있다. 특히, 터모신테시스(tomosynthesis) 시스템과 같이 저선량을 사용하는 진단 시스템에 적용 시 정밀 보정이 더욱 중요하게 되어 기존 방법 대비 높은 효과를 거둘 수 있다. According to the image balance correction method and apparatus according to the present invention, by using an appropriate correction method in consideration of the response characteristics of each section of the image sensor instead of the conventional image correction method, while reducing the accuracy of the correction image acquisition process and application complexity while increasing the accuracy of the correction As a result, the image quality and reliability of the final result of the radiation detection image can be obtained. In particular, when applied to a diagnostic system that uses low doses, such as the thermosynthesis system, the precision correction becomes more important than the conventional method.

        본 발명과 본 발명의 동작상의 이점 및 본 발명의 실시에 의하여 달성되는 목적을 충분히 이해하기 위해서는 본 발명의 바람직한 실시 예를 예시하는 첨부 도면 및 첨부 도면에 기재된 내용을 참조하여야만 한다. In order to fully understand the present invention, the operational advantages of the present invention, and the objects achieved by the practice of the present invention, reference should be made to the accompanying drawings which illustrate preferred embodiments of the present invention and the contents described in the accompanying drawings.

이하, 첨부한 도면을 참조하여 본 발명의 바람직한 실시예를 설명함으로써, 본 발명을 상세히 설명한다.Hereinafter, exemplary embodiments of the present invention will be described in detail with reference to the accompanying drawings.

본 발명에서는 X-선 등 방사선을 이용해 환자의 암이나 기타 환부의 상태를 진단하기 위한 진단 장비에서 반도체 이미지 센서 등을 이용해 영상 데이터를 획득하여 해당 영상을 제공하는 과정에서, 이미지 센서로부터 검출된 각 화소의 디지털 영상 신호를 이미지 센서의 구간별 응답 특성을 고려한 적절한 보정방식을 사용하여 보정함으로써, 보정의 정확도는 높이면서도 보정 이미지 획득 과정과 적용의 복잡도를 낮춰, 결과적으로 방사선 검출 이미지의 최종적인 결과물의 화질과 신뢰성을 개선할 수 있도록 하였다. In the present invention, in the diagnostic equipment for diagnosing the condition of cancer or other lesions of the patient using X-rays and other radiation, the angle detected from the image sensor in the process of acquiring the image data using a semiconductor image sensor and the like and providing the corresponding image By correcting the digital image signal of the pixel using an appropriate correction method that takes into account the response characteristics of each section of the image sensor, the accuracy of the correction is increased while the complexity of the acquisition and application of the correction image is reduced, resulting in the final result of the radiation detection image. To improve the image quality and reliability.

본 발명에서, 위와 같은 이미지 센서로부터의 디지털 영상 신호를 보정하여 이미지의 밸런스(또는 평탄화)를 위한 보정 방법은, 입력 휘도값별 출력 휘도값을 나타내는 도 1과 같은 이미지 센서의 응답 특성 데이터에 대한 2차 기울기 정보(도 3 참조)를 임계치(예를 들어, 0) 상하로 구분되는 복수의 구간으로 나누는 과정과, 이와 같은 구간별 응답 특성을 고려하여, 이미지 센서로부터 수신한 화소 데이터에 대하여 위와 같이 나누어진 복수의 구간 중 해당 구간에 대응된 보정 방식으로 상기 수신된 화소 데이터를 보정하는 과정으로 이루어진다. 도 3과 같은 2차 기울기 정보는 도 1과 같은 응답 특성 데이터에 대한 1차 기울기 정보(도 2참조)에 대한 기울기를 계산함으로써 획득될 수 있다. In the present invention, the correction method for the balance (or flattening) of the image by correcting the digital image signal from the image sensor as described above, 2 for the response characteristic data of the image sensor as shown in FIG. Dividing the difference slope information (refer to FIG. 3) into a plurality of sections that are divided into upper and lower thresholds (for example, 0) and the pixel data received from the image sensor in consideration of the response characteristics of the sections as described above. Compensating the received pixel data in a correction method corresponding to the corresponding section of the plurality of divided sections. The second slope information as shown in FIG. 3 may be obtained by calculating a slope for the first slope information (see FIG. 2) with respect to the response characteristic data as shown in FIG. 1.

도 1과 같이 입력 휘도값별(예를 들어, 휘도 0~255 입력 디지털 영상 신호) 출력 휘도값(예를 들어, 휘도 0~255 출력 디지털 영상 신호)을 나타내는 이미지 센서의 비선형 응답 특성 데이터에서 볼 수 있는 바와 같이, 이를 전체적으로 선형 또는 2차 이상의 다차 함수에 따른 보정 방식으로 처리 할 경우 보정 에러가 발생하며 또한 각 방식에 따른 단점이 존재하게 된다.As shown in FIG. 1, the nonlinear response characteristic data of an image sensor indicating an output luminance value (eg, luminance 0 to 255 output digital image signal) for each input luminance value (eg, luminance 0 to 255 input digital image signal) may be viewed. As can be seen, when this is processed as a correction method according to the linear or quadratic multi-order function as a whole, a correction error occurs and there are disadvantages according to each method.

예를 들어, 도1 과 같은 형태의 응답 특성 데이터에서 구간별(A, B, C)로 서로 다른 응답 특성을 가지고 있음을 알 수 있다. 즉, 낮은 입력 휘도값을 가지는 구간(A)은 2차 곡선형태의 응답특성으로 볼 수 있으며, 중간 구간(B)은 선형 응답 특성으로 볼 수 있고, 높은 입력 휘도값을 가지는 구간(C)에서는 다시 2차 곡선 형태가 된다. 따라서, 낮은 입력 휘도값을 가지는 구간(A)은 응답 특성이 2차 곡선의 형태를 띠게 되므로, 2차 이상의 다차 함수에 기초한 비선형 보정 방식을 적용하 고, 중간 구간(B)에서는 1차 함수에 기초한 선형 보정 방식을 적용하며, 마지막 구간(C)에서는 역시 2차 이상의 다차 함수에 기초한 비선형 보정 방식을 적용함으로써, 좀더 정확한 보정을 하면서도 구간의 수를 줄이고 실 측정값에 적용 시의 연산량도 적절하도록 계산될 수 있게 된다. For example, it can be seen that the response characteristics data of FIG. 1 have different response characteristics for each section (A, B, C). That is, the section A having a low input luminance value can be seen as a response characteristic in the form of a quadratic curve, and the middle section B can be viewed as a linear response characteristic, and in the section C having a high input luminance value. It is again quadratic. Therefore, in the section A having a low input luminance value, since the response characteristic is in the form of a quadratic curve, a nonlinear correction method based on a second or more ordered quadratic function is applied. Based on the linear correction method, the last section (C) also applies the nonlinear correction method based on the quadratic or higher quadratic function, so that the number of sections is reduced and the calculation amount is appropriate when applied to the actual measured value while performing more accurate correction. Can be calculated.

이때 위와 같은 적절한 구간의 분배는 다음과 같은 방법으로 이루어 질 수 있다. 먼저, 소정 기울기 계산 수단을 이용해 도 1과 같은 응답 특성 데이터에 대한 기울기 정보(예를 들어, 각 입력 휘도값에 대한 기울기 값의 정보)(1차 기울기 정보)를 생성할 수 있다. 이와 같이 응답 특성 데이터에 대한 기울기 정보는 도 2와 같은 형태로 나타낼 수 있다. 즉, 응답 특성 데이터에 대한 기울기 정보는 A 구간에서 지속적으로 증가하는 기울기 값을 나타내다가, B 구간에서 일정한 기울기 값을 나타내며, 다시 C 구간에서 기울기 값이 감소 하는 형태로 나타난다. At this time, the distribution of the appropriate section as described above may be made in the following way. First, inclination information (for example, inclination information for each input luminance value) for the response characteristic data as shown in FIG. 1 may be generated using a predetermined inclination calculation means (first order inclination information). As such, the slope information on the response characteristic data may be represented as shown in FIG. 2. That is, the slope information on the response characteristic data indicates a slope value that is continuously increased in the section A, a constant slope value in the section B, and the slope value decreases in the section C again.

여기서 다시 소정 기울기 계산 수단을 이용해 도 2와 같은 기울기 정보에 대한 기울기 정보(예를 들어, 각 입력 휘도값에 대한 기울기 값의 정보)를 생성할 수 있다. 이때 C 구간에서는 기울기 값이 음의 값을 가지므로 절대값을 취하여 나타내면, 이때의 기울기 정보(2차 기울기 정보)는 도 3과 같은 형태로 나타낼 수 있다. 도3과 같은 기울기 정보는 응답 특성의 변화를 반영하게 되며 이렇게 얻어진 도3과 같은 기울기 정보가 일정한 임계치를 넘게 되는 포인트를 구간 구분값으로 설정할 수 있다. 또한, 원하는 정밀도에 따라 임계치를 적절히 조정 할 수 있다. 만약, 임계치를 0으로 정하게 되면 도 3과 같이 3개 구간(A, B, C)으로 나누어 질 수 있으며, A 구간(2차 기울기값이 임계치 이상의 구간)에 대하여는 응답 특성이 2차 곡 선의 형태를 띠게 되므로 2차 이상의 다차 함수에 기초한 비선형 보정 방식을 적용하고, B 구간(2차 기울기값이 임계치 보다 작은 구간)에서는 1차 함수에 기초한 선형 보정 방식을 적용하며, 마지막 C 구간(2차 기울기값이 임계치 이상의 구간)에서는 역시 2차 이상의 다차 함수에 기초한 비선형 보정 방식을 적용함으로써, 좀더 정확한 보정을 하면서도 구간의 수를 줄이고 실 측정값에 적용 시의 연산량도 적절하도록 계산될 수 있게 된다. Here again, the slope information (for example, the information of the slope value for each input luminance value) with respect to the slope information may be generated using the predetermined slope calculation means. In this case, since the inclination value has a negative value in the section C, when the absolute value is taken and displayed, the inclination information (secondary inclination information) at this time may be represented as shown in FIG. 3. The inclination information as shown in FIG. 3 reflects the change in response characteristics, and thus, the point at which the inclination information as shown in FIG. 3 exceeds a predetermined threshold value may be set as the interval division value. In addition, the threshold can be appropriately adjusted according to the desired precision. If the threshold value is set to 0, it can be divided into three sections (A, B, and C) as shown in FIG. 3, and in response to the section A (section 2 having a second slope greater than or equal to the threshold), response characteristics are in the form of a secondary curve. Since the nonlinear correction method based on the quadratic function of the quadratic or higher order is applied, the linear correction method based on the linear function is applied in the B section (the second slope is smaller than the threshold), and the last C section (the second slope). In the case where the value is greater than or equal to the threshold value, the nonlinear correction method based on the second-order or higher quadratic function can also be applied to reduce the number of intervals and to calculate the calculation amount when applied to the actual measured value while performing more accurate correction.

여기서, 비선형 보정은, 해당 구간(2차 기울기값이 임계치 이상의 구간)에 대하여 응답 특성 데이터에 대하여 소정 수치해석법을 적용해 그 데이터와 유사한 커브를 갖는2차 이상의 다차 함수로 피팅(fitting)하여 해당 다차 함수를 획득하고, 획득한 다차 함수에 기초하여 이미지 센서로부터 수신한 화소 데이터를 보정하는 것이다. Here, the nonlinear correction is performed by applying a predetermined numerical method to response characteristic data for a corresponding section (section where the second inclination value is greater than or equal to the threshold value) and fitting it to a quadratic or higher quadratic function having a curve similar to the data. Acquiring a multi-order function, and correcting the pixel data received from the image sensor based on the obtained multi-order function.

예를 들어, 화소 데이터를 비선형 보정하기 위하여는, 먼저, 위와 같이 각 화소마다 다차 함수를 결정하고, 소정 계산 수단으로 화소 데이터가 나타내는 휘도값의 해당 구간에 대한 전체 화소의 휘도 평균값을 구한다. 이와 같은 화소 데이터가 나타내는 휘도값의 해당 구간에 대한 전체 화소의 휘도 평균값은 해당 구간의 다차 함수의 평균을 이용해 구할 수 있다. 이와 같이 획득된 평균값을 수신된 화소 데이터와 곱한 후, 곱해진 값을 해당 구간에서의 화소의 평균 휘도에 대응된 휘도값으로 나누어 줌으로써 화소 데이터의 비선형 보정이 완성될 수 있다. 여기서 해당 구간에서의 화소의 평균 휘도는 전체 화소가 아닌 자신 화소의 해당 구간에서의 평균적 휘도이다.For example, in order to nonlinearly correct pixel data, first, a multi-order function is determined for each pixel as described above, and the luminance average value of all pixels for a corresponding section of the luminance value represented by the pixel data is determined by a predetermined calculation means. The luminance average value of all the pixels of the corresponding interval of the luminance value indicated by the pixel data may be obtained using the average of the polynomial functions of the corresponding interval. The non-linear correction of the pixel data may be completed by multiplying the average value thus obtained by the received pixel data and dividing the multiplied value by the luminance value corresponding to the average luminance of the pixel in the corresponding section. Here, the average luminance of the pixels in the corresponding section is not the all pixels but the average luminance of the corresponding pixels in the corresponding section.

또한, 선형 보정은, 해당 구간(2차 기울기값이 임계치 보다 작은 구간)에 대하여 응답 특성 데이터에 대하여 소정 수치해석법을 적용해 그 데이터와 유사한 커브를 갖는1차 함수로 피팅(fitting)하여 해당 1차 함수를 획득하고, 획득한 함수에 기초하여 이미지 센서로부터 수신한 화소 데이터를 보정하는 것이다. In addition, linear correction is performed by applying a predetermined numerical analysis method to response characteristic data for a corresponding section (section where the second slope value is smaller than the threshold value) and fitting it to a linear function having a curve similar to that data. The difference function is obtained, and the pixel data received from the image sensor is corrected based on the obtained function.

예를 들어, 화소 데이터를 선형 보정하기 위하여는, 먼저, 위와 같이 각 화소마다 1차 함수를 결정하고, 소정 계산 수단으로 화소 데이터가 나타내는 휘도값의 해당 구간에 대한 전체 화소의 휘도 평균값을 구한다. 이와 같은 화소 데이터가 나타내는 휘도값의 해당 구간에 대한 전체 화소의 휘도 평균값은 해당 구간의 1차 함수의 평균을 이용해 구할 수 있다. 이와 같이 획득된 평균값을 수신된 화소 데이터와 곱한 후, 곱해진 값을 해당 구간에서의 화소의 평균 휘도에 대응된 휘도값으로 나누어 줌으로써 화소 데이터의 선형 보정이 완성될 수 있다. 여기서 해당 구간에서의 화소의 평균 휘도는 전체 화소가 아닌 자신 화소의 해당 구간에서의 평균적 휘도이다.For example, in order to linearly correct pixel data, first, a linear function is determined for each pixel as described above, and the luminance average value of all pixels for a corresponding section of the luminance value represented by the pixel data is determined by a predetermined calculation means. The luminance average value of all the pixels of the corresponding interval of the luminance value represented by the pixel data may be obtained using the average of the linear function of the corresponding interval. After multiplying the obtained average value by the received pixel data, the linear correction of the pixel data may be completed by dividing the multiplied value by the luminance value corresponding to the average luminance of the pixel in the corresponding section. Here, the average luminance of the pixels in the corresponding section is not the all pixels but the average luminance of the corresponding pixels in the corresponding section.

이와 같은 선형 및 비선형 보정을 통해 도 4와 같이 방사선 검출 이미지의 각 화소 마다 휘도 특성이 평탄화되도록 이미지 밸런스가 보정될 수 있다. 도 4의 410은 보정전 전체 화소의 일정 휘도값에 대한 영상 출력의 결과이고, 도 4의 420은 같은 휘도값에 대하여 위와 같은 선형 및 비선형 보정 방식을 적용한 후의 영상 출력의 결과이다. Through such linear and nonlinear correction, the image balance may be corrected such that luminance characteristics are flattened for each pixel of the radiation detection image as shown in FIG. 4. 410 of FIG. 4 is a result of image output for a constant luminance value of all pixels before correction, and 420 of FIG. 4 is a result of image output after applying the linear and nonlinear correction schemes to the same luminance value.

이와 같은 이미지 밸런스를 위한 보정 방법은 도 5와 같은 본 발명의 일실시예에 따른 이미지 보정 장치(500)에 의하여 구현될 수 있다.The correction method for the image balance may be implemented by the image correction device 500 according to an embodiment of the present invention as shown in FIG.

도 5를 참조하면, 본 발명의 일실시예에 따른 이미지 보정 장치(500)는, 제1 기울기 계산부(510), 제2 기울기 계산부(520), 구간 설정부(530), 및 보정부(540)를 포함한다.Referring to FIG. 5, the image correcting apparatus 500 according to an embodiment of the present invention may include a first tilt calculator 510, a second tilt calculator 520, a section setting unit 530, and a corrector. 540.

제1 기울기 계산부(510)는 도 1과 같은 입력 휘도값별 출력 휘도값을 나타내는 이미지 센서의 응답 특성 데이터를 수신하여 해당 기울기 정보를 생성하여 도 2와 같은 1차 기울기 정보를 생성한다. 1차 기울기 정보는 각 입력 휘도값에 대한 응답 특성 데이터의 기울기 값에 대한 정보이다.The first slope calculator 510 receives response characteristic data of an image sensor indicating output luminance values for each input luminance value as shown in FIG. 1, generates corresponding slope information, and generates first slope information as shown in FIG. 2. The first slope information is information on a slope value of response characteristic data for each input luminance value.

제2 기울기 계산부(520)는 상기 1차 기울기 정보에 대한 기울기 정보를 구하고 절대값을 계산하여 도 3과 같은 2차 기울기 정보를 생성한다. 2차 기울기 정보는 각 입력 휘도값에 대한 1차 기울기 정보의 기울기 값에 대한 정보이다.The second slope calculator 520 obtains slope information of the first slope information, calculates an absolute value, and generates second slope information as illustrated in FIG. 3. The second slope information is information on a slope value of the first slope information for each input luminance value.

구간 설정부(530)는 위와 같이 계산된 2차 기울기 정보를 임계치(예를 들어, 0) 상하로 구분되는 복수의 구간으로 나누어 각 구간을 식별하기 위한 정보를 생성한다. 도 3은 임계치를 0으로 한 경우에 2차 기울기 정보가 3개 구간(A, B, C)으로 나뉘어지는 것을 나타낸 것으로서 임계치를 다른 값으로 설정하는 경우에 임계치 상하로 구분되는 2차 기울기 정보의 구간들은 더 많은 구간으로 나뉘어질 수도 있다. The section setting unit 530 generates information for identifying each section by dividing the second slope information calculated as above into a plurality of sections that are divided into upper and lower thresholds (for example, 0). FIG. 3 illustrates that secondary slope information is divided into three sections A, B, and C when the threshold is set to 0. When the threshold is set to another value, the second slope information divided into upper and lower thresholds is shown. The sections may be divided into more sections.

보정부(540)는 이미지 센서로부터 수신하는 화소 데이터에 대하여, 구간 설정부(530)에서 설정한 구간 식별 정보에 따른 복수의 구간 중 해당 구간에 대응된 보정 방식으로 상기 화소 데이터를 보정하게 된다. 예를 들어, 보정부(540)는 위와 같은 해당 임계치 이상의 구간에 대하여 비선형 보정을 수행할 수 있고, 해당 임계 치 보다 작은 구간에 대하여 선형 보정을 수행할 수 있다.The correction unit 540 corrects the pixel data with respect to the pixel data received from the image sensor by a correction method corresponding to the corresponding section among the plurality of sections according to the section identification information set by the section setting unit 530. For example, the correction unit 540 may perform nonlinear correction on a section above the threshold, and may perform linear correction on a section below the threshold.

예를 들어, 도 3의 경우에, 보정부(540)는 A 구간(2차 기울기값이 임계치 이상의 구간)에 대하여는 응답 특성이 2차 곡선의 형태를 띠게 되므로 2차 이상의 다차 함수에 기초한 비선형 보정 방식을 적용하고, B 구간(2차 기울기값이 임계치 보다 작은 구간)에서는 1차 함수에 기초한 선형 보정 방식을 적용하며, 마지막 C 구간(2차 기울기값이 임계치 이상의 구간)에서는 역시 2차 이상의 다차 함수에 기초한 비선형 보정 방식을 적용하여, 수신한 화소 데이터를 보정할 수 있다.For example, in the case of FIG. 3, the correction unit 540 has a nonlinear correction based on the quadratic function of the quadratic or higher because the response characteristic is in the form of a quadratic curve for the section A (the quadratic slope is greater than or equal to the threshold). In the B section (secondary slope value is smaller than the threshold), the linear correction method is applied based on the linear function, and in the last C section (secondary slope value is above the threshold), the second order or more The received pixel data may be corrected by applying a nonlinear correction method based on a function.

여기서, 비선형 보정을 위하여 보정부(540)는, 해당 구간(2차 기울기값이 임계치 이상의 구간)에 대하여 응답 특성 데이터에 대하여 소정 수치해석법을 적용해 그 데이터와 유사한 커브를 갖는2차 이상의 다차 함수로 피팅(fitting)하여 해당 다차 함수를 획득하고, 획득한 다차 함수에 기초하여 이미지 센서로부터 수신한 화소 데이터를 보정할 수 있다. Here, for nonlinear correction, the correction unit 540 applies a predetermined numerical analysis method to response characteristic data for a corresponding section (section where the second inclination value is greater than or equal to the threshold value) and has a quadratic or higher quadratic function having a curve similar to the data. The corresponding multi-order function may be obtained by fitting a function, and the pixel data received from the image sensor may be corrected based on the obtained multi-order function.

예를 들어, 화소 데이터를 비선형 보정하기 위하여, 먼저, 보정부(540)는 위와 같은 방법을 적용하여 도 6과 같이 각 화소마다 다차 함수를 결정하고, 화소 데이터가 나타내는 휘도값의 해당 구간에 대한 전체 화소의 휘도 평균값을 구할 수 있다. 이와 같은 화소 데이터가 나타내는 휘도값의 해당 구간에 대한 전체 화소의 휘도 평균값(각 화소의 커브들에 대한 평균 커브에서의 평균 휘도)은 해당 구간의 다차 함수의 평균을 이용해 구해질 수 있다. 보정부(540)는 이와 같이 획득된 평균값을 수신된 화소 데이터와 곱한 후, 곱해진 값을 해당 구간에서의 화소의 평균 휘도(해당 화소의 커브에서의 평균 휘도)에 대응된 휘도값으로 나누어 줌으로써 화소 데이터의 비선형 보정을 완성할 수 있다. 여기서 해당 구간에서의 화소의 평균 휘도는 전체 화소가 아닌 자신 화소의 해당 구간에서의 평균적 휘도이다.For example, in order to nonlinearly correct pixel data, first, the correction unit 540 determines the multi-order function for each pixel as shown in FIG. The luminance average value of all the pixels can be obtained. The luminance average value (average luminance in the average curve for the curves of each pixel) of all the pixels for the corresponding interval of the luminance value represented by such pixel data may be calculated using an average of the multi-order function of the corresponding interval. The correction unit 540 multiplies the average value thus obtained by the received pixel data and divides the multiplied value by the luminance value corresponding to the average luminance (average luminance in the curve of the pixel) of the pixel in the corresponding section. Nonlinear correction of the pixel data can be completed. Here, the average luminance of the pixels in the corresponding section is not the all pixels but the average luminance of the corresponding pixels in the corresponding section.

또한, 선형 보정을 위하여 보정부(540)는, 해당 구간(2차 기울기값이 임계치 보다 작은 구간)에 대하여 응답 특성 데이터에 대하여 소정 수치해석법을 적용해 그 데이터와 유사한 커브를 갖는1차 함수로 피팅(fitting)하여 해당 1차 함수를 획득하고, 획득한 함수에 기초하여 이미지 센서로부터 수신한 화소 데이터를 보정할 수 있다. 예를 들어, 화소 데이터를 선형 보정하기 위하여, 먼저, 보정부(540)는 위와 같은 방법을 적용하여 도 6과 같이 각 화소마다 1차 함수를 결정하고, 소정 계산 수단으로 화소 데이터가 나타내는 휘도값의 해당 구간에 대한 전체 화소의 휘도 평균값을 구할 수 있다. 이와 같은 화소 데이터가 나타내는 휘도값의 해당 구간에 대한 전체 화소의 휘도 평균값은 해당 구간의 1차 함수의 평균을 이용해 구해질 수 있다. 보정부(540)는 이와 같이 획득된 평균값을 수신된 화소 데이터와 곱한 후, 곱해진 값을 해당 구간에서의 화소의 평균 휘도(해당 화소의 커브에서의 평균 휘도)에 대응된 휘도값으로 나누어 줌으로써 화소 데이터의 선형 보정을 완성할 수 있다. 여기서도 해당 구간에서의 화소의 평균 휘도는 전체 화소가 아닌 자신 화소의 해당 구간에서의 평균적 휘도이다.In addition, for linear correction, the correction unit 540 applies a predetermined numerical method to response characteristic data for a corresponding section (a section in which the second inclination value is smaller than the threshold value) to a linear function having a curve similar to the data. The first linear function may be acquired by fitting, and the pixel data received from the image sensor may be corrected based on the obtained function. For example, in order to linearly correct pixel data, first, the correction unit 540 determines the linear function of each pixel as shown in FIG. The luminance average value of all pixels for the corresponding interval of may be obtained. The luminance average value of all pixels of the corresponding interval of the luminance value represented by the pixel data may be calculated using the average of the linear function of the corresponding interval. The correction unit 540 multiplies the average value thus obtained by the received pixel data and divides the multiplied value by the luminance value corresponding to the average luminance (average luminance in the curve of the pixel) of the pixel in the corresponding section. Linear correction of the pixel data can be completed. Here, the average luminance of the pixels in the corresponding section is not the entire pixel but the average luminance of the corresponding pixels in the corresponding section.

이와 같은 선형 및 비선형 보정을 통해, 도 4의 410과 같이 방사선 검출 이미지의 보정전 전체 화소에서의 평탄화되지 않은 영상이, 보정 후에는, 도 4의 420과 같이 전체 화소에서 평탄화되어 밸런스가 이루어진 영상으로 출력될 수 있다.Through such linear and nonlinear correction, an unflattened image of all the pixels before correction of the radiation detection image as shown in 410 of FIG. 4 is flattened and balanced in all the pixels of 420 of FIG. 4 after correction. It can be output as

이상에서와 같이 도면과 명세서에서 최적 실시예가 개시되었다. 여기서 특정 한 용어들이 사용되었으나, 이는 단지 본 발명을 설명하기 위한 목적에서 사용된 것이지 의미한정이나 특허청구범위에 기재된 본 발명의 범위를 제한하기 위하여 사용된 것은 아니다. 그러므로 본 기술 분야의 통상의 지식을 가진 자라면 이로부터 다양한 변형 및 균등한 타 실시예가 가능하다는 점을 이해할 것이다. 따라서, 본 발명의 진정한 기술적 보호 범위는 첨부된 특허청구범위의 기술적 사상에 의해 정해져야 할 것이다. As described above, optimal embodiments have been disclosed in the drawings and the specification. Although specific terms have been used herein, they are used only for the purpose of describing the present invention and are not intended to limit the scope of the present invention as defined in the claims or the claims. Therefore, those skilled in the art will understand that various modifications and equivalent other embodiments are possible from this. Therefore, the true technical protection scope of the present invention will be defined by the technical spirit of the appended claims.

도 1은 방사선 이미지 데이터의 비선형 응답 특성을 설명하기 위한 입력 신호에 대한 출력 신호를 나타내는 그래프이다.1 is a graph showing an output signal with respect to an input signal for explaining the nonlinear response characteristics of radiographic image data.

도 2는 도 1의 기울기를 나타내는 그래프이다.2 is a graph illustrating the inclination of FIG. 1.

도 3은 도 2의 기울기를 나타내는 그래프이다.3 is a graph illustrating the inclination of FIG. 2.

도 4는 본 발명의 일실시예에 따른 이미지 보정 방법에 따른 이미지 처리의 결과를 설명하기 위한 도면이다.4 is a view for explaining a result of image processing according to the image correction method according to an embodiment of the present invention.

도 5는 본 발명의 일실시예에 따른 이미지 보정 장치의 블록도이다.5 is a block diagram of an image correction device according to an embodiment of the present invention.

도 6은 본 발명의 일실시예에 따른 보정부의 선형/비선형 보정 수행에 대하여 설명하기 위한 그래프이다.6 is a graph illustrating a linear / nonlinear correction performed by a correction unit according to an embodiment of the present invention.

Claims (12)

입력 휘도값별 출력 휘도값을 나타내는 이미지 센서의 응답 특성 데이터에 대한 2차 기울기 정보를 임계치 이상의 구간과 상기 임계치 보다 작은 구간을 포함하는 복수의 구간으로 나누는 단계; 및 Dividing the second slope information on the response characteristic data of the image sensor indicating the output luminance value for each input luminance value into a plurality of sections including a section above a threshold and a section below the threshold; And 상기 이미지 센서로부터 수신한 화소 데이터의 상기 복수의 구간 중 해당 구간에 대응된 보정 방식으로 상기 화소 데이터를 보정하는 단계Correcting the pixel data in a correction scheme corresponding to a corresponding section among the plurality of sections of pixel data received from the image sensor. 를 포함하는 것을 특징으로 하는 이미지의 밸런스를 위한 보정 방법.Correction method for the balance of the image comprising a. 제1항에 있어서, The method of claim 1, 상기 화소 데이터는 방사선 진단 장비에서 방사선에 대한 영상 데이터를 획득하는 이미지 센서로부터 검출된 각 화소의 디지털 영상 신호인 것을 특징으로 하는 이미지의 밸런스를 위한 보정 방법.And the pixel data is a digital image signal of each pixel detected from an image sensor for obtaining image data of radiation in a radiation diagnostic equipment. 제1항에 있어서,The method of claim 1, 상기 2차 기울기 정보는 상기 응답 특성 데이터에 대한 기울기 정보에 대하여 재차 기울기를 계산하고 절대값을 취한 정보인 것을 특징으로 하는 이미지의 밸런스를 위한 보정 방법.And the second inclination information is information obtained by calculating an inclination again and taking an absolute value with respect to the inclination information with respect to the response characteristic data. 제1항에 있어서, 상기 임계치는 0인 것을 특징으로 하는 이미지의 밸런스를 위한 보정 방법.The method of claim 1, wherein the threshold is zero. 제1항에 있어서, 상기 보정하는 단계는,The method of claim 1, wherein the correcting step, 상기 임계치 이상의 구간에 대하여 비선형 보정을 수행하는 단계; 및Performing nonlinear correction on a section above the threshold; And 상기 임계치 보다 작은 구간에 대하여 선형 보정을 수행하는 단계Performing linear correction on a section smaller than the threshold 를 포함하는 것을 특징으로 하는 이미지의 밸런스를 위한 보정 방법.Correction method for the balance of the image comprising a. 제5항에 있어서, The method of claim 5, 상기 비선형 보정은, 해당 구간에 대하여 상기 응답 특성 데이터에 대하여 피팅한 2차 이상의 다차 함수에 기초하여 상기 화소 데이터를 보정하는 것이고,The nonlinear correction is to correct the pixel data based on a quadratic or higher quadratic function fitted to the response characteristic data for a corresponding section, 상기 선형 보정은, 해당 구간에 대하여 상기 응답 특성 데이터에 대하여 피팅한 1차 함수에 기초하여 상기 화소 데이터를 보정하는 것을 특징으로 하는 이미지의 밸런스를 위한 보정 방법.And the linear correction corrects the pixel data based on a linear function fitted to the response characteristic data for a corresponding section. 제6항에 있어서, The method of claim 6, 상기 비선형 보정은, 상기 다차 함수에 의하여 결정된 해당 구간에 대한 전체 화소의 휘도 평균값을 상기 화소 데이터와 곱한 후, 곱해진 값을 해당 구간에서의 화소의 평균 휘도에 대응된 휘도값으로 나누는 보정이고,The nonlinear correction is a correction that multiplies the luminance average value of all pixels for the corresponding section determined by the multi-order function with the pixel data, and then divides the multiplied value by the luminance value corresponding to the average luminance of the pixels in the corresponding section, 상기 선형 보정은, 상기 1차 함수에 의하여 결정된 해당 구간에 전체 화소의 휘도 평균값을 상기 화소 데이터와 곱한 후, 곱해진 값을 해당 구간에서의 화소의 평균 휘도에 대응된 휘도값으로 나누는 보정인 것을 특징으로 하는 이미지의 밸런 스를 위한 보정 방법.The linear correction is a correction that multiplies the luminance average value of all pixels in the interval determined by the linear function with the pixel data, and then divides the multiplied value by the luminance value corresponding to the average luminance of the pixels in the interval. A correction method for balancing an image characterized by. 입력 휘도값별 출력 휘도값을 나타내는 이미지 센서의 응답 특성 데이터에 대한 2차 기울기 정보를 임계치 이상의 구간과 상기 임계치 보다 작은 구간을 포함하는 복수의 구간으로 나누는 구간 설정부; 및 A section setting unit dividing second gradient information on response characteristic data of an image sensor indicating output luminance values for each input luminance value into a plurality of sections including a section above a threshold and a section below the threshold; And 상기 이미지 센서로부터 수신하는 화소 데이터의 상기 복수의 구간 중 해당 구간에 대응된 보정 방식으로 상기 화소 데이터를 보정하는 보정부A correction unit configured to correct the pixel data in a correction scheme corresponding to a corresponding section among the plurality of sections of pixel data received from the image sensor. 를 포함하는 것을 특징으로 하는 이미지의 밸런스를 위한 보정 장치.Correction apparatus for the balance of the image, characterized in that it comprises a. 제8항에 있어서,The method of claim 8, 상기 응답 특성 데이터에 대한 기울기 정보를 생성하는 제1 기울기 계산부; 및A first slope calculator configured to generate slope information on the response characteristic data; And 상기 제1 기울기 계산부로부터의 기울기 정보에 대한 기울기 정보를 구하고 절대값을 계산하여 상기 2차 기울기 정보를 생성하는 제2 기울기 계산부A second slope calculator which obtains slope information on the slope information from the first slope calculator and calculates an absolute value to generate the second slope information 를 더 포함하는 것을 특징으로 하는 이미지의 밸런스를 위한 보정 장치.Correction apparatus for the balance of the image, characterized in that it further comprises. 제8항에 있어서, 상기 보정부는,The method of claim 8, wherein the correction unit, 상기 임계치 이상의 구간에 대하여 비선형 보정을 수행하고, 상기 임계치 보다 작은 구간에 대하여 선형 보정을 수행하는 것을 특징으로 하는 이미지의 밸런스를 위한 보정 장치.And performing non-linear correction on the section above the threshold, and performing linear correction on the section below the threshold. 제10항에 있어서, 상기 보정부는,The method of claim 10, wherein the correction unit, 상기 비선형 보정으로서, 해당 구간에 대하여 상기 응답 특성 데이터에 대하여 피팅한 2차 이상의 다차 함수에 기초하여 상기 화소 데이터를 보정하고,As the nonlinear correction, the pixel data is corrected based on a quadratic or higher quadratic function fitted to the response characteristic data for the interval, 상기 선형 보정으로서, 해당 구간에 대하여 상기 응답 특성 데이터에 대하여 피팅한 1차 함수에 기초하여 상기 화소 데이터를 보정하는 것을 특징으로 하는 이미지의 밸런스를 위한 보정 장치.And correcting the pixel data based on a linear function fitted to the response characteristic data for the section as the linear correction. 제11항에 있어서, 상기 보정부는,The method of claim 11, wherein the correction unit, 상기 비선형 보정에서, 상기 다차 함수에 의하여 결정된 해당 구간에 대한 전체 화소의 휘도 평균값을 상기 화소 데이터와 곱한 후, 곱해진 값을 해당 구간에서의 화소의 평균 휘도에 대응된 휘도값으로 나누고,In the nonlinear correction, a luminance average value of all pixels for the corresponding interval determined by the multi-order function is multiplied by the pixel data, and then the multiplied value is divided by the luminance value corresponding to the average luminance of the pixels in the corresponding interval, 상기 선형 보정에서, 상기 1차 함수에 의하여 결정된 해당 구간에 전체 화소의 휘도 평균값을 상기 화소 데이터와 곱한 후, 곱해진 값을 해당 구간에서의 화소의 평균 휘도에 대응된 휘도값으로 나누는 것을 특징으로 하는 이미지의 밸런스를 위한 보정 장치.In the linear correction, a luminance average value of all pixels is multiplied by the pixel data in the corresponding interval determined by the first order function, and then the multiplied value is divided by the luminance value corresponding to the average luminance of the pixels in the corresponding interval. Correction device for the balance of the image.
KR1020090116235A 2009-11-27 2009-11-27 Method and Apparatus for Correcting Balance of Digital Radiation Detection Image having Nonlinear Response Characteristics KR101092207B1 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
KR1020090116235A KR101092207B1 (en) 2009-11-27 2009-11-27 Method and Apparatus for Correcting Balance of Digital Radiation Detection Image having Nonlinear Response Characteristics

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
KR1020090116235A KR101092207B1 (en) 2009-11-27 2009-11-27 Method and Apparatus for Correcting Balance of Digital Radiation Detection Image having Nonlinear Response Characteristics

Publications (2)

Publication Number Publication Date
KR20110059486A KR20110059486A (en) 2011-06-02
KR101092207B1 true KR101092207B1 (en) 2011-12-12

Family

ID=44394631

Family Applications (1)

Application Number Title Priority Date Filing Date
KR1020090116235A KR101092207B1 (en) 2009-11-27 2009-11-27 Method and Apparatus for Correcting Balance of Digital Radiation Detection Image having Nonlinear Response Characteristics

Country Status (1)

Country Link
KR (1) KR101092207B1 (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101347254B1 (en) * 2012-09-10 2014-01-07 연세대학교 원주산학협력단 Correction method for non-uniformity of photon counting detector

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2002236179A (en) 2001-02-06 2002-08-23 Canon Inc Radiation detection device and radiation imaging system using it

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2002236179A (en) 2001-02-06 2002-08-23 Canon Inc Radiation detection device and radiation imaging system using it

Also Published As

Publication number Publication date
KR20110059486A (en) 2011-06-02

Similar Documents

Publication Publication Date Title
CN108158597B (en) Method and device for determining original X-ray energy data and CT (computed tomography) equipment
US7156554B2 (en) Method of compensating for image faults in an x-ray image recording
US7632016B1 (en) Digital detector calibration with known exposure
US6460003B1 (en) Apparatus and method for resolution calibration of radiographic images
US9360571B2 (en) System and method for correction of vignetting effect in multi-camera flat panel x-ray detectors
CN104599234A (en) Image processing apparatus, radiation imaging system, and image processing method
US20230306657A1 (en) Noise suppression using deep convolutional networks
US8406376B2 (en) Radiographic apparatus
US10182783B2 (en) Visualization of exposure index values in digital radiography
CN115266779A (en) Flat panel detector pixel response nonlinear correction method for CT scanning
US7496221B2 (en) Radiological image diagnostic system and data processing method thereof
US20010048734A1 (en) Method for dynamic range management of digital radiographic images
US20070116348A1 (en) Adaptive image processing and display for digital and computed radiography images
US7274770B2 (en) X-ray examination apparatus and a method of controlling an output of an X-ray source of an X-ray examination apparatus
CN105374014A (en) Image correction method, image correction device, medical image generation method and medical image generation device
US8611625B2 (en) Tomographic apparatus
KR101092207B1 (en) Method and Apparatus for Correcting Balance of Digital Radiation Detection Image having Nonlinear Response Characteristics
CN106651768B (en) Image correction method and device and X-ray photographing equipment
US20090190716A1 (en) Sensitometric response mapping for radiological images
TW201529044A (en) Radiation imaging device
JP2002034961A (en) Radiographing apparatus and radiographing method
US8611498B2 (en) Image processing apparatus, image processing method, radiation imaging system, and computer-readable recording medium
JP2004195233A (en) Method and apparatus to reduce spectrum-sensitive type artifacts
JPH05130987A (en) Beam hardening correcting method
CN111096761B (en) Method, device and related equipment for correcting scattering of wedge-shaped filter

Legal Events

Date Code Title Description
A201 Request for examination
E902 Notification of reason for refusal
E701 Decision to grant or registration of patent right
GRNT Written decision to grant
FPAY Annual fee payment

Payment date: 20141204

Year of fee payment: 4

FPAY Annual fee payment

Payment date: 20161206

Year of fee payment: 6

FPAY Annual fee payment

Payment date: 20171201

Year of fee payment: 7

FPAY Annual fee payment

Payment date: 20181204

Year of fee payment: 8

FPAY Annual fee payment

Payment date: 20191204

Year of fee payment: 9