JBiSE  Vol.9 No.11 , October 2016
Dark Current Compensation and Sensitivity Adjustment on Gallium Arsenide Linear Array Detector for X-Ray Imaging
ABSTRACT
For the last several years, the linear array x-ray detector for x-ray imaging with gallium arsenide direct conversion sensitive elements has been developed and tested at the In-stitute for High Energy Physics. The array consists of 16 sensitive modules. Each module has 128 gallium arsenide (GaAs) sensitive elements with 200 μm pitch. Current article describes two key program procedures of initial dark current compensation of each sensitive element in the linear array, and sensitivity adjustment for alignment of strip pattern in the raw image data. As a part of evaluation process a modular transfer function (MTF) was measured with the slanted sharp-edge object under RQA5 technique as it described in the International Electrotechnical Commission 62220-1 standard (high voltage 70 kVp, additional aluminium filter 21 mm) for images with compensated dark currents and adjusted sensitivity of detector elements. The 10% level of the calculated MTF function has spatial resolution within 2 - 3 pair of lines per mm for both vertical and horizontal orientation.

1. Introduction

The using of semiconductor materials with direct conversion of x-ray photon energy into electron charges is the one of the primary techniques in developing modern detector systems for different kinds of applications. System based on selenium and cadmium telluride (CdTe and CZT) semiconductors are being used for medical purpose and destructive testing both in current integration and counting mode. In other hand the gallium arsenide (GaAs) material doesn't used in commercial x-ray detection systems yet, but there were several examples of linear arrays [1] - [4] developed in the past decade. All of them were using different types of detector sensitive elements arrangement into array, different ASIC to form raw data for acquisition system and different acquisition system elements to form an actual image. The formation of the images in those systems requires several procedures to transform raw data into image “for processing” as DICOM 3.0 standard requires for medical imaging. In case of linear array detectors with a direct conversion there are two sources of non-uniformity: the raw sequence of signals formed by different dark current of each element without x-ray flux and sensitivity of each element to x-ray exposure. There are many different algorithms to reduce this type non-uniformity which available mainly for infrared imaging systems [5] - [8] and because infrared linear arrays have very similar non-uniformity pattern. But in case of x-ray imaging such algorithms are the proprietary methods, unavailable as software libraries. This article will focus on two algorithms of software dark current and of detector elements sensitivity compensation developed for the GaAs arrays.

2. Detector Array Structure and Adjustment

2.1. Linear Array Structure

The key unit of the linear array is a module, consisting of 128 sensitive elements of gallium arsenide chromium doped resistive type material which was designed, manufactured and assembled in Siberian Physical-Technical Institute of Tomsk State University [9] [10] . The pitch size of a single element can be 100 × 100 μm2 or 200 × 200 μm2 and the sensitive length of material vary from 600 μm for pad elements up to 4 mm for strip elements. Figure 1 shows two detectors modules of strip sensitive elements with 200 μm pixel pitch, mounted on the printed circuit board. Arrows show the direction of photon flux. The maximum array length was 400 mm assembled from 16 modules. Due to technological reasons is it impossible to make all the modules with same sensitivity

Figure 1. Simplified example of strip sensitive elements arrangement on detector module and two detector modules on an acquisition board.

to the x-ray exposure, moreover the sensitivity of pads and strips vary not only from one module to another but also inside each module. The transformation of signals from detectors into useful voltage signal is done by 128-channel ASIC chip [11] . Every channel of ASIC consists of: input circuit for matching of bias voltage of sensitive element, personal 8-bit compensation current generator and 8-bit storage memory, compensation current generator for all channels, integrator and buffer. ASIC output doesn’t have any programmable gain amplifiers to adjust the sensitivity of each element. That is the reason why we have to do a programmable software adjustment. The connection between sensitive elements and ASIC is made by ultrasound welding. The chip is converting the electron charges, formed by the x-ray flux into serial sequence of voltages.

2.2. Principle of Control and Adjustment

There are two possibilities to control the output serial voltage signals of ASIC. The first one is to adjust the average level of output signals and the second one is to adjust input current for every sensitive element to minimize the spread of dark currents of all detectors on the module. To adjust both average level of output signals value and individual dark currents the chip has several analog and digital inputs. The adjustment of the mean level is controlled manually with the potentiometer and with the software regulations of the dark currents by writing 8-bit value codes in ASIC. These two possibilities allow control a shadow signals only when there is no x-ray flux.

For each sensitive element there is a transfer characteristic function of output voltage in dependence of input 8-bit code from 0 to 255. By calculation of optimal input codes (sometimes it’s called a mask) it is possible to get output serial voltage signals as a straight line. The current of each input of ASIC includes three components, they are: Idet-current formed by sensitive element under x-ray exposure, I1-average current of all 128 sensitive elements, I2-compensation current, k-input 8-bit value. The output voltage signal for every one is proportional to input current and clock frequency.

(1)

(2)

Clock frequency is constant during the exposure and acquisition. Voltage sequence transformed by ADC into counts has a form as it showed in a Figure 2.

The ADC count has three main ranges (3), two of them are constants, and one is the working range where the dark current can be compensated.

(3)

To make all signals have the same or close to same Abase count value it’s necessary to calculate Nbase input value for each chip sensitive element input. To do that it’s necessary to form a response matrix A.

Figure 2. Transfer characteristic of the ADC count of the single sensitive element according to 8-digit input code of ASIC.

(4)

where m is number of input values in chip, and n is amount of sensitive elements. In our case the maximum m is equal to 255 and n is equal to 128. The m values could be an array of equidistant distribution with different step from 0 to 255, e.g. .

The easiest way to find the Nj = Nbase or it nearest neighbor is to calculate the index of minimum absolute difference between Abase and each element of column vector A.

(5)

where

After finding all compensation values for the array and writing them into ASICs, in result the output voltage signal will be with minimum spread. If an equidistant distribution step not equal to 1 the Figure 3 shows the spread of ADC counts from a single module.

2.3. Compensation of Different Sensitive Elements X-Ray Response

However, there is still a different sensitivity of detector elements to x-ray flux, even when all values for dark current compensation are correct and output signals without x-ray flux are a uniform. That pattern depends on many parameters such as bias voltage of detector elements, sampling clock frequency, high voltage of x-ray tube etc. To reduce such pattern we use a calibration test-object which has a shape of a ramp. The idea of sensitivity adjustment is based on a calculated response from a reference uniform signal for different intensities. Such shape allows to simulate the variation of x-ray flux intensity and get a response of each sensitive element under the same uniform x-ray flux. In our case the maximum difference in the ramp thickness is 20 mm that covers the major part of studies where x-ray high voltage is equal or less 70 kVp. Figure 4 demonstrates the ramp object scanned from top to bottom by the linear array with 5 detector modules. The data presented after dark current compensation. The vertical pattern represents the different sensitivity from one detector element to another.

For a typical x-ray tube that has a bremsstrahlung spectrum as well as characteristic lines, the mono-energetic formula of intensity attenuation like (6)

(6)

can’t be used to calculate the response of each sensitive element and compensate it, thus the response of each sensitive element is represented as a polynomial dependency (7), where d is a ramp thickness and ki are coefficients, which allows to take into account both the scattering and absorption of x-ray photons.

(7)

Figure 3. The box and whisker spreads of output amplitudes after dark current compensation for different number of equidistant steps.

Figure 5 shows, as example, a simulation of spectrum of the x-ray tube source with tungsten target, with 70 kVp anode voltage, and built-in 2.5 mm plus additional 20 mm aluminum filter. The spectrum generated with Monte-Carlo simulation program XMIM-SIM [12] [13] . The detector current response for such difference in intensity and x-ray energy is easier to describe as a polynomial expression.

To find the ki coefficients we use the data of the ramp object exposure. The data of a raw image for sensitivity compensation is B, the ramp test-object is represented as a matrix R where each row is an array’s data on the same ramp thickness. The vector Z is the values proportional to di from (7). Number of rows in R and Z is the same.

(8)

In order to compensate the amplitude difference between all sensitive elements, we calculate an adjustment value for each pixel in B.

To compensate sensitivity of all elements to the same level we use a ramp as a reference

Figure 4. Reference ramp object for sensitivity adjustment of detector linear array.

Figure 5. Simulation of an x-ray tube spectrum with different aluminum filters.

and for each raw image pixel we calculate a minimum index of the absolute difference between pixel and ramp matrix column. Column vector Y is formed for each element of matrix B.

(9)

As a result the compensated pixel will be

(10)

In order to increase a number of rows in matrix R the interpolation could be done by calculating additional data for all rows. Equations (11) and (7) allow to find ki coefficients with linear least square approach and recalculate row data for expanded Z.

(11)

The resulted image formed from pixels, calculated from equation 10 is represented as “for processing” in DICOM 3.0 format for further PACS systems.

3. Examples of Sensitivity Adjustment on Raw Images

Images described below were obtained with 3 modules which have minimum pixel intensity fluctuation due to scanner array movement and defects in input channels of ASIC which affects on output signals. Raw image example with compensated dark currents is showed in Figure 6 together with adjusted sensitivity by using the ramp object calibration and by MIRE [14] algorithm. The scan was done from top to the bottom and the exposure time was 1.4 seconds. Algorithms which use only image data itself without additional information will leave artifacts in regions with significant intensity gradients.

The acquisition parameters were:

・ clock frequency: 19 MHz (1.2 MHz per detector module)

・ raw scan-line formation time: 444 us

・ scan-line formation time after resample to 1:1 scale: 2.2 ms

・ detector voltage: −20V

・ detector operation temperature: +12.5˚C.

・ high voltage of x-ray tube: 80 kVp.

・ focal distance: 100 cm.

One of the main characteristics of the x-ray imaging system is a modular transfer

Figure 6. Example of sensitivity adjustment, from left to right: raw image, adjusted with MIRE algorithm, adjusted using ramp object.

function (MTF), which is measured by means of slanted edge test object throughout the additional 20 mm aluminium filter. High voltage of x-ray tube during MTF measurement was 70 kVp and the distance between x-ray tube and test-object was 95 cm. Figure 7 represents raw and adjusted images of slanted sharp-edge object oriented in parallel and perpendicular to array detector during the scan. Resulted MTF functions were calculated in ImageJ [15] program using two plugins SE_MTF [16] and COQ [17] [18] .

Examples of images before and after sensitivity adjustment are shown on Figure 8 from full length of linear array of 16 sensitive modules and 2048 sensitive elements. Resulted MTF is represented in Figure 9 for two rescale filters (Cubic and Lanczos) of raw image obtained with the same exposition parameters.

The 10% level of MTF function lies within the 2 - 3 lp/mm range for both vertical and horizontal orientation of the slanted edge in agreement with the spatial resolution pattern image.

Another very important evaluation parameter on development stage of a scanning array is a normalized noise power spectrum (NNPS). It is fluctuating because of different characteristics of detector modules and ASICs. NNPS is presented in Figure 10 for the resulted image with 20 mm aluminum filter under uniform exposition for vertical, horizontal and radial direction of the scan. NPS data were measured with 3 modules (384 sensitive elements). The peak in vertical direction represents periodic electronic noise in output signals from some ASICs.

The presented technique of sensitivity adjustment has several disadvantages such as need a calibration ramp or other type of uniform object. The misalignment between a ramp object and a linear array may cause additional artifacts in a final image after the adjustment procedure. In case of the large thickness of the object the brightest or darkest

Figure 7. Images of slanted edge object before and after detector array sensitivity adjustment.

Figure 8. Examples of images before and after sensitivity adjustment.

pixels in the image could be adjusted wrongly and some artifacts will still remain.

Figure 9. Modular transfer function for two convolution matrices Cubic and Lanczos (raw image rescale filters) with the same exposure and acquisition parameters.

4. Results

There are two algorithms described in article: dark current compensation algorithm to reduce signal variance on each detective module, and sensitivity adjustment to compensate sensitivity of each detection element. The represented algorithm of dark current compensation allows minimizing spread of output 12-bit ADC values up to 30 ADC counts (sigma of Gaussian distribution). The noise reduction of the raw image data by sensitivity adjustment via reference ramp object could be a key element of a program algorithm when it’s impossible to use data of the raw image itself and the resulted image quality is crucial. Different exposure parameters, such as different high voltages of x-ray tube and different anode currents can be adjusted well with the usage of several ramp objects. This procedure must be adopted initially before any other exposures with real objects or patients. The same technique could be used with number of

Figure 10. Normalized noise power spectrum for horizontal, vertical and radial directions of the scan.

images with uniform exposure of the whole region of interest. The modular transfer function and spatial resolution for GaAs linear array detector were measured with dark current compensation and sensitivity adjustment of detector elements using the developed algorithms with reasonable output. For 200 um pitch sensitive array element 10% MTF has value between 2 - 3 pair of lines per mm.

Cite this paper
Polkovnikov, M. (2016) Dark Current Compensation and Sensitivity Adjustment on Gallium Arsenide Linear Array Detector for X-Ray Imaging. Journal of Biomedical Science and Engineering, 9, 532-543. doi: 10.4236/jbise.2016.911047.
References
[1]   Ardashev, E.N., Gorokhov, S.A., Polkovnikov, M.K., Lobanov, I.S. and Vorobiev, A.P. (2011) A Portable X-Ray Apparatus with GaAs Linear Array. Nuclear Instruments and Methods in Physics Research Section A, 648, S8-S11.
http://dx.doi.org/10.1016/j.nima.2011.02.030

[2]   Dubecky, F., Perd’ochová, A., Scepko, P., et al. (2005) Digital X-Ray Portable Scanner Based on Monolithic Semi-Insulating GaAs Detectors: General Description and First “Quantum” Images. Nuclear Instruments and Methods in Physics Research Section A, 546, 118-124.
http://dx.doi.org/10.1016/j.nima.2005.03.020

[3]   Dvoryankin, V.F., Dikaev, Yu.M., Krikunov, A.I., Panova, T.M. and Telegin, A.A. (2004) Multielement X-Ray Row Detector on GaAs with Spatial Resolution of 108 μm. Nuclear Instruments and Methods in Physics Research Section A, 531, 87-88.
http://dx.doi.org/10.1016/j.nima.2004.05.097

[4]   Dvoryankin, V.F., Dvoryankina, G.G., Dikaev, Yu.M., Ermakov, M.G., Kudryashov, A.A., Petrov, A.G. and Telegin, A.A. (2012). Photovoltaic GaAs Detectors for Digital X-Ray Imaging, Mammography - Recent Advances, Dr. Nachiko Uchiyama (Ed.), InTech.
http://dx.doi.org/10.5772/32431

[5]   Tendero, Y. and Gilles, J. (2012) ADMIRE: A Locally Adaptive Single-Image, Non-Uniformity Correction and Denoising Algorithm: Application to Uncooled IR Camera. Proceedings SPIE 8353, Infrared Technology and Applications XXXVIII, 83531O.

[6]   Harris, J.G. and Chiang, Y.M. (1999) Nonuniformity Correction of Infrared Image Sequences Using the Constant-Statistics Constraint. IEEE Transactions on Image Processing, 8, 1148-1151.
http://dx.doi.org/10.1109/83.777098

[7]   Li, H.-C. and Song, H.-F. (2011) A Rel-ative Radiometric Correction Method for Linear Array Push-Broom Imagery. Proceedings SPIE 8194, International Symposium on Photoelectronic Detection and Imaging 2011: Advances in Imag-ing Detectors and Applications, 81941L.
http://dx.doi.org/10.1117/12.900195

[8]   Ma, W.H., Dong, T., Tian, H. and Ni, J.P. (2014) Line-Scan CCD Camera Calibration in 2D Coordinate Measurement. Optik - International Journal for Light and Electron Optics, 125, 4795-4798.
http://dx.doi.org/10.1016/j.ijleo.2014.04.057

[9]   Tyazhev, A.V., Budnitsky, D.L., Koretskay, O.B., Novikov, V.A., Okaevich, L.S., Potapov, A.I., Tolbanov, O.P. and Vorobiev, A.P. (2003) GaAs Radiation Imaging Detectors with an Active Layer Thickness up to 1 mm. Nuclear Instruments and Methods in Physics Research Section A, 509, 34-39.
http://dx.doi.org/10.1016/S0168-9002(03)01545-6

[10]   Veale, M.C., Bell, S.J., Duarte, D.D., French, M.J., Schneider, A., Seller, P., Wilson, M.D., Lozinskaya, A.D., Novikov, V.A., Tolbanov, O.P., Tyazhev, A. and Zarubin, A.N. (2014) Chromium Compensated Gallium Arsenide Detectors for X-Ray and γ-Ray Spectroscopic Imaging. Nuclear Instruments and Methods in Physics Research Section A, 762, 6-14.
http://dx.doi.org/10.1016/j.nima.2014.03.033

[11]   Borodin, D.V., Osipov, V.V., Shushkevich, N.A. and Lopuhin A.A. (2003) Testing the Readout Integrated Circuits. Applied Physics, 1, 154-156 (in Russian).
http://applphys.orionir.ru

[12]   XMI-MSIM: Monte Carlo Simulation of Energy-Dispersive X-Ray Fluorescence Spectrometers.
http://github.com/tschoonj/xmimsim

[13]   Schoonjans, T., Brunetti, A., Vincze, L., et al. (2011) The Xraylib Library for X-Ray-Matter Interactions. Recent Developments. Spectrochimica Acta B, 66, 776-784.
http://dx.doi.org/10.1016/j.sab.2011.09.011

[14]   Tendero, Y., Landeau, S. and Gilles, J. (2012) Non-Uniformity Correction of Infrared Images by Midway Equalization. Image Processing On Line, 2, 134-146.
http://dx.doi.org/10.5201/ipol.2012.glmt-mire

[15]   ImageJ. https://imagej.nih.gov/ij/index.html

[16]   SE_MTFplugin.
http://rsb.info.nih.gov/ij/plugins/se-mtf/index.html

[17]   Nitrosi, A., Bertolini, M., et al. (2009) Application of QC_DR Software for Acceptance Testing and Routine Quality Control of Direct Digital Radiography Systems: Initial Experiences using the Italian Association of Physicist in Medicine Quality Control Protocol. Journal of Digital Imaging, 22, 656-666.
http://dx.doi.org/10.1007/s10278-008-9150-z

[18]   Donini, B., Rivetti, S., Lanconelli, N., et al. (2014) Free Software for Performing Physical Analysis of Systems for Digital Radiography and Mammography. Medical Physics, 41, Article ID: 051903.
http://dx.doi.org/10.1118/1.4870955

 
 
Top