Advances in OptoElectronics

Volume 2015, Article ID 496863, 5 pages

http://dx.doi.org/10.1155/2015/496863

## Block-Based Compressed Sensing for Neutron Radiation Image Using WDFB

Faculty of Information Science & Engineering, Ningbo University, Ningbo 315211, China

Received 13 December 2014; Accepted 23 March 2015

Academic Editor: Michele Norgia

Copyright © 2015 Wei Jin et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

An ideal compression method for neutron radiation image should have high compression ratio while keeping more details of the original image. Compressed sensing (CS), which can break through the restrictions of sampling theorem, is likely to offer an efficient compression scheme for the neutron radiation image. Combining wavelet transform with directional filter banks, a novel nonredundant multiscale geometry analysis transform named Wavelet Directional Filter Banks (WDFB) is constructed and applied to represent neutron radiation image sparsely. Then, the block-based CS technique is introduced and a high performance CS scheme for neutron radiation image is proposed. By performing two-step iterative shrinkage algorithm the problem of norm minimization is solved to reconstruct neutron radiation image from random measurements. The experiment results demonstrate that the scheme not only improves the quality of reconstructed image obviously but also retains more details of original image.

#### 1. Introduction

Digital neutron radiography, one type of nondestructive testing technology, is a newly developing high-technology which is the combination of a set of technologies such as nuclear science, optoelectromechanical integration, and digital image processing [1]. With the development of digital neutron imaging system, the dynamic range of neutron radiation image is improved [2]; meanwhile the amount of image data is increased rapidly. The huge amount of data occupies plenty of the storage space and network resources and sometimes even becomes the bottleneck of the whole system performance improvement. For the traditional image compression methods are hard to preserve the details which are the most valuable part of the neutron radiation images at high compression ratio, how to find a rational method to reduce the amount of data of neutron radiation image has become an increasingly urgent issue. In recent years, a new theory of compressed sensing (CS) has been put forward [3, 4], and it points out that for the “sparse” signal (or the signal’s transformed coefficients are “sparse” in a specific transform domain), using a random measurement matrix, the original signal is projected onto a low dimensional space to get a small part of the random measurements, and by solving an optimization problem the original signal can be reconstructed from this part of the measurements. Because the neutron radiation image meets the sparsity condition, by using the rational designed measurement matrix, the original image can be reconstructed from small projection measurements. With the help of this scheme, the sampling rate of neutron imaging system and the amount of neutron radiation image may be reduced effectively. In general, in order to reconstruct the ideal image from the random measurements, three issues should be concerned including sparse representation, measurement matrix, and reconstruction algorithm. Combining wavelet transform with directional filter banks [5, 6], a novel nonredundant multiscale geometry analysis transform named Wavelet Directional Filter Banks (WDFB) is constructed and applied to represent neutron radiation image sparsely. In order to ensure that the linear projection of neutron radiation image can keep the original information and reduce the time and space complexity of sparse measurement, we introduce block-based compressed sensing (BCS) technology [4] by using the blocked measurement matrix which satisfies restricted isometry property (RIP) to perform random measurement. Finally, the original neutron radiation image was reconstructed from the measurement by solving the problem of norm minimization based on two-step iterative shrinkage reconstruction algorithm [7, 8]. The proposed scheme is applied to the neutron radiation image and compared with the traditional method in terms of reconstruction speed, visual effect, and the objective evaluation of reconstructed image.

#### 2. Wavelet Directional Filter Banks (WDFB)

We know that the parts of dramatic change in the image, such as edge and contour, are the important local characteristics which often contain more information than other parts. In order to represent such directional features of neutron radiation image sparsely, we use direction filter bank (DFB) to implement multidirection decomposition to express the high dimensional singularity of image effectively. Considering that DFB can not depict image in multiresolution manner and its poor performance to represent the low frequency component of image, we use wavelet decomposition combined with directional filter banks to construct a nonredundant multiscale geometric analysis transform, and we call this transform Wavelet Directional Filter Banks (WDFB).

By using WDFB, multiscale analysis is performed via wavelet transform to capture point singularity of the image, followed by DFB to synthesize the singular points in the same direction into a coefficient. In a word, WDFB use the basic structure like contour segment to approximate the original image and the support interval of the basis function has an elongated structure with the change of directions and scales. Compared with wavelet analysis, WDFB can represent the image more sparsely. Specifically, the image is decomposed into an approximation subband and several detail subbands , , and after level wavelet decomposition; then, WDFB carries out DFB to the detail subbands of wavelet decomposition to obtain the direction information of image. In order to meet the anisotropy scaling relation, in the realization of DFB decomposition, the direction numbers were decreased in dyadic manner from fine to coarse scale. That is to say, if the direction numbers of the finest scale are , then the numbers of its adjacent coarser scale will be , and so on until the coarsest scale. In practical applications, the decomposition levels of wavelet and the maximum direction numbers of DFB are determined by the characteristics of actual image and the specific tasks.

#### 3. Block-Based Compressed Sensing

In recent years, CS technology where theoretical frameworks were established by Donoho, Candes, and so forth has been developed quickly and attracted the attention of scholars. Suppose that the signal with length can be represented sparsely by the orthogonal basis or tight frame, that is, , that is to say, there are only nonzero coefficients in domain (), we call a sparse signal in domain (in this paper, refers to performing WDFB transform for ). In that case, an measurement matrix can be constructed and the linear measurements can be obtained by

Then, the original signal can be reconstructed perfectly from measurements by solving the optimization problem based on the known :where norm counts the nonzero elements in the corresponding vector. Considering that the traditional CS is not suitable for sensing natural image in real-time mode for the sampling process requiring accessing the entire target at once and the reconstruction algorithms are generally very expensive, in [2], a block-based CS (BCS) for natural images was proposed. Motivated by the great success of block DCT coding systems, in BCS scheme, the original image is divided into small blocks and each block is sampled independently using the same measurement matrix.

BCS can reduce the size of compressed sensing measurement matrix significantly; moreover, the computational complexity of the reconstruction algorithm can be greatly reduced. Consider an original image with pixels, where and represent the rows and columns of image, respectively; we divided the image into small blocks with size of . Let represent the vectorized signal of the th block through raster scanning. Then, design an orthonormalized i.i.d. Gaussian matrix with size of as the measurement matrix for each subblock. According to the traditional CS theory, the th block can be reconstructed from random measurements . Finally, the original image can be feasibly restored by integrating accordingly. Obviously, if the measurement rate is set as , then . For there is a trade-off in the selection of block size , small requires less memory in storage and faster implementation but causes mosaic in the reconstructed image, while large offers better reconstruction performance but causes high computational and storage complexity. From empirical studies, we set block size in this paper.

#### 4. Neutron Radiation Image Reconstruction Using Two-Step Iterative Shrinkage/Threshold

How to reconstruct the original signal from linear measurements is the key to the application of CS technique. The norm optimization like (2) is both numerically unstable and nondeterministic polynomial (NP) complete problem, and it can hardly be solved in polynomial time so several alternative solution procedures have been proposed [9]. Perhaps the most prominent of these is basis pursuit (BP) which applies a convex relaxation to the problem resulting in an optimization. Although BP can be implemented effectively with linear programming, its computational complexity is often high. Other variant algorithms include conjugate gradient (CG), matching pursuit (MP), and orthogonal matching pursuit (OMP). Unfortunately, such algorithms reduce computational complexity at the cost of lower reconstruction quality.

In Daubechies et al.’s [7], optimization can be solved by iterative shrinkage/threshold (IST) algorithm. Bioucas-Dias and Figueiredo [8] further develop this algorithm and propose a two-step iterative shrinkage/threshold (TwIST) algorithm in order to improve the convergence speed of IST. Unlike IST, in the iteration of TwIST, the new estimate depends on the two previous estimates so that the optimized results can be obtained more quickly. In this paper, we use TwIST to realize BCS reconstruction for the neutron radiation image. The following is the description of the proposed scheme in detail.

*Step 1. *The original neutron radiation image is divided into small blocks with size . represent the vectorized signal of the th block through raster scanning and then performing random measurement using the measurement matrix .

*Step 2. *Estimating the initial reconstruction subblock , set iteration counter .

*Step 3. *In order to reduce blocking effect, Wiener-filtering was carried out to the reconstructed image which was integrated by each subblock and then calculating the objective function corresponds to aswhere denotes WDFB transform, the parameter can be regarded as a regularization term, and tuning this parameter can make balance between two parts of the objective function.

*Step 4. *Thresholding process , where the threshold function iswhere threshold is hard thresholding function, to set a proper threshold , we imply the universal threshold method , and is the number of the transform coefficients, and is estimated using a robust median estimator. is the iteration step and set its initial value equal to 1; then, is calculated using (3); if set to recalculate ; otherwise continue to the next step.

*Step 5. *Updating the iteration counter, , and then using two previous estimates to update the current value,where the parameters and determine the speed of convergence given by [8]Calculating the objective function , if then return to Step 3; otherwise continue to the next step.

*Step 6. *The iteration will be terminated if , where is the preset small number; else return to Step 5. is defined asThen the final reconstruction image can be obtained via integrating accordingly.

#### 5. Experimental Results and Analysis

In order to verify the effectiveness of the proposed scheme, the compressed sensing experiments were performed on the neutron radiation image of a test specimen. The digital neutron imaging system was carried out by a virtual radiography station using Monte Carlo simulation. In this simulated system, the neutron source has a cold spectrum; the diameter of the collimator is 4 cm at the narrowest part. The specimen is positioned 100 cm from the aperture. The scintillating screen is located 20 cm behind the object and the neutron radiation image is recorded by use of the ^{6}Li(n,a)-reaction. The test specimen is an aluminum cylinder which contains iron hoop, iron balls, screw, polyethylene rack, cadmium bar, and so forth. Figure 1(a) is the profile of the test specimen while Figure 1(b) is the corresponding neutron radiation image. The test specimen simulates the structure of an engine and has a good representative for the applications of neutron radiography in nondestructive testing.