Published

2015-09-01

High-dynamic range compressive spectral imaging by grayscale coded aperture adaptive filtering

Reconstrucción de imágenes espectrales compresivas con alto rango dinámico por medio de filtrado adaptativo en escala de grises

DOI:

https://doi.org/10.15446/ing.investig.v35n3.49868

Keywords:

Compressive sensing, hyperspectral imaging, saturation, dynamic range, coded aperture, optical imaging. (en)
Muestreo compresivo, imágenes hiperespectrales, saturación, rango dinámico, apertura codificada, imágenes ópticas. (es)

Downloads

Authors

  • Nelson Eduardo Diaz M,Sc Student at Universidad Industrial de Santander, Colombia.
  • Hoover Fabian Rueda Chacon Ph.D student at University of Delaware, Newark DE, USA.
  • Henry Arguello Fuentes Associated professor at Universidad Industrial de Santander, Colombia.

The coded aperture snapshot spectral imaging system (CASSI) is an imaging architecture which senses the three dimensional informa-tion of a scene with two dimensional (2D) focal plane array (FPA) coded projection measurements. A reconstruction algorithm takes advantage of the compressive measurements sparsity to recover the underlying 3D data cube. Traditionally, CASSI uses block-un-block coded apertures (BCA) to spatially modulate the light. In CASSI the quality of the reconstructed images depends on the design of these coded apertures and the FPA dynamic range. This work presents a new CASSI architecture based on grayscaled coded apertu-res (GCA) which reduce the FPA saturation and increase the dynamic range of the reconstructed images. The set of GCA is calculated in a real-time adaptive manner exploiting the information from the FPA compressive measurements. Extensive simulations show the attained improvement in the quality of the reconstructed images when GCA are employed.  In addition, a comparison between traditional coded apertures and GCA is realized with respect to noise tolerance.

El sistema de adquisición de imágenes espectrales basado en una apertura codificada de única captura (CASSI) es una arquitectura para muestrear imágenes espectrales que captura información de una escena en tres dimensiones mediante medidas proyectadas y codificadas del sensor en dos dimensiones (2D). Para recuperar el cubo de datos se emplea un algoritmo de reconstrucción que apro-vecha la esparcidad de las medidas compresivas. Tradicionalmente, CASSI usa aperturas codificadas binarias para modular espacial-mente la luz.  En CASSI la calidad de las imágenes reconstruidas depende del diseño de estas aperturas codificadas y de la saturación del sensor. Este trabajo presenta una nueva arquitectura CASSI  basada en aperturas codificadas de escala de grises (ACG), la cual reduce la saturación en el sensor e incrementa el rango dinámico de las imágenes reconstruidas.  El conjunto de ACG se calcula en tiempo real de forma adaptativa explotando la información contenida en las medidas compresivas del sensor. La mejora obtenida en la calidad de las imágenes reconstruidas se muestra a partir de diversas simulaciones en las que se utilizan las ACG. Adicionalmente, se lleva a cabo una comparación entre las aperturas codificadas tradicionales y las ACG con respecto a la tolerancia del ruido.
ART_6

DOI: https://doi.org/10.15446/ing.investig.v35n3.49868

High-dynamic range compressive spectral imaging
by grayscale coded aperture adaptive filtering

Reconstrucción de imágenes espectrales compresivas con alto rango dinámico por medio de filtrado adaptativo en escala de grises

N. Diaz1, H. Rueda2, and H. Arguello3


1 Nelson E. Diaz Diaz: BSc. in Computer Sciences and Master Student in Computer Sciences at Universidad Industrial de Santander, Colombia. Affiliation: Master Student at Universidad Industrial de Santander, Colombia. E-mail: nelson.diaz@correo.uis.edu.co

2 Hoover Rueda: BSc in Computer Sciences and MSc at Universidad Industrial de Santander, Colombia. Ph.D student at University of Delaware, Newark DE, USA. Affiliation: Ph.D student at University of Delaware, Newark DE, USA. E-mail: rueda@udel.edu

3 Henry Arguello: F. BSc in Electrical engineer and MSc in Electrical power at Universidad Industrial de Santander, Colombia. Ph.D in Electrical and Computer engineering from the University of Delaware, USA. Affiliation: Associated professor at Universidad Industrial de Santander, Colombia. E-mail: henarfu@uis.edu.co

How to cite: Diaz, N., Rueda, H., & Arguello, H. (2015). High-dynamic range compressive spectral imaging by grayscale coded aperture adaptive filtering. Ingeniería e Investigación, 35(3), 53-60. DOI: https://doi.org/10.15446/ing.investig.v35n3.49868


ABSTRACT

The coded aperture snapshot spectral imaging system (CASSI) is an imaging architecture which senses the three dimensional information of a scene with two dimensional (2D) focal plane array (FPA) coded projection measurements. A reconstruction algorithm takes advantage of the compressive measurements sparsity to recover the underlying 3D data cube. Traditionally, CASSI uses block-unblock coded apertures (BCA) to spatially modulate the light. In CASSI the quality of the reconstructed images depends on the design of these coded apertures and the FPA dynamic range. This work presents a new CASSI architecture based on grayscaled coded apertures (GCA) which reduce the FPA saturation and increase the dynamic range of the reconstructed images. The set of GCA is calculated in a real-time adaptive manner exploiting the information from the FPA compressive measurements. Extensive simulations show the attained improvement in the quality of the reconstructed images when GCA are employed. In addition, a comparison between traditional coded apertures and GCA is realized with respect to noise tolerance.

Keywords: Compressive sensing, hyperspectra imaging, saturation, dynamic range, coded aperture, optical imaging.


RESUMEN

El sistema de adquisición de imágenes espectrales basado en una apertura codificada de única captura (CASSI) es una arquitectura para muestrear imágenes espectrales que captura información de una escena en tres dimensiones mediante medidas proyectadas y codificadas del sensor en dos dimensiones (2D). Para recuperar el cubo de datos se emplea un algoritmo de reconstrucción que aprovecha la esparcidad de las medidas compresivas. Tradicionalmente, CASSI usa aperturas codificadas binarias para modular espacialmente la luz. En CASSI la calidad de las imágenes reconstruidas depende del diseño de estas aperturas codificadas y de la saturación del sensor. Este trabajo presenta una nueva arquitectura CASSI basada en aperturas codificadas de escala de grises (ACG), la cual reduce la saturación en el sensor e incrementa el rango dinámico de las imágenes reconstruidas. El conjunto de ACG se calcula en tiempo real de forma adaptativa explotando la información contenida en las medidas compresivas del sensor. La mejora obtenida en la calidad de las imágenes reconstruidas se muestra a partir de diversas simulaciones en las que se utilizan las ACG. Adicionalmente, se lleva a cabo una comparación entre las aperturas codificadas tradicionales y las ACG con respecto a la tolerancia del ruido.

Palabras clave: Muestreo compresivo, imágenes hiperespectrales, saturación, rango dinámico, apertura codificada, imágenes ópticas.


Received: March 27th 2015

Accepted: July 29th 2015


Introduction

The coded aperture snapshot spectral imaging system (CASSI) is an imaging architecture which senses the three dimensional information of a scene with a single two dimensional (2D) coded random projection measurement set (Wagadarikar, et al. 2008). The CASSI optical architecture comprises five optical elements: an objective lens is used to form an image of a scene in the plane of the coded aperture; a coded aperture modulates the spatial information over the complete wavelength range; a relay lens transmits the coded light field onto a dispersive element that disperses the light before it impinges on the focal plane array (FPA). Given a set of compressive measurements, compressive sensing theory (CS) (Candes 2006, Donoho 2006, Baraniuk 2007) is used to reconstruct the underlying data cube of size N × N × L from just N (N + L 1) measurements, where N is the spatial dimension and L is the spectral depth of the data cube. In CASSI, the quality of the reconstructed images relies on the design of the set of 2D coded apertures, whose entries block or transmit the light from the scene. These coded apertures are called block-unblock coded apertures (BCA).

A single shot in CASSI may not provide sufficient number of compressive measurements. A recent modification in CASSI allows multi-shot sensing procedures which increase the number of compressive measurements, (Arguello and Arce 2013, Arguello and Arce 2014, Rueda and Arguello 2013, Rueda, et al, 2014). In this modification, each shot uses a distinct coded aperture that remains fixed during the integration time of the detector. The quality of the reconstructed images improves in CASSI in proportion to the number of compressive measurements (Arguello and Arce, 2011). Each CASSI shot adds simultaneously N(N + L - 1) compressive measurements. Therefore, the total number of available measurements with K shots is KN(N + L 1).

The FPA’s dynamic range is the ratio between the smallest and the largest value detected by the sensor. The largest value is known as saturation level. Saturation occurs when compressive measurements exceed the saturation level, i.e., the quantizer’s dynamic range. In that case, the measurements take the value of the saturation level (Laska, et al, 2011). In CASSI and multi-shot CASSI, each saturated pixel in the sensor induces errors in the reconstructed image. Typically, CASSI system employs CCD or CMOS sensor, both affected by saturation depending on their dynamic range.

In order to avoid saturation the coded apertures can be used to modulate the incoming light in the system. Traditionally, BCA are implemented using piezo systems (Kittle, et al, 2010) or a digital-micromirror-device (DMD) (Wu, et al, 2011), to vary the coding pattern in each snapshot. The disadvantage of BCA is that their entries can take just two binary values {0,1}; that is, the modulation of the intensity to binary entries. In order to improve intensity modulation, we propose the use of grayscale-adaptive coded apertures (GCA) which can be also implemented using a DMD. The new design takes advantage of the fast switching time of the micro-mirrors which enables the use of a pulse-width modulation technique for the production of grayscale values. GCA can be used to yield a modulation of the intensity and to increase the dynamic range of the reconstructions.

This work extends the compressive capabilities of CASSI by replacing the traditional BCA by a set of GCA. GCA-CASSI multi-shot is motivated by the possibility of reducing saturation levels through modulation of the amplitude of the incoming scene.

Figure 1a shows the sketch of the proposed architecture. The adaptive process is accomplished through a PC real-time model (Feedback), in which the PC receives the measurement Gi from the FPA and remembers the previous compressive measurements Gi+1, Gi+2, Gi+3 and the previous coded aperture patterns. Using the adaptive filter (AF), the micro-mirrors are updated in the DMD, and the AF attenuates the coded aperture between snapshots, yielding an adaptively generated GCA . Figure 1b shows the detail of the grayscale-adaptive coded aperture where the attenuated pixel represents an oscillating DMD. In Figure 1b the pixels from the DMD are selected with 25 %, 50 %, and 75 % of duty cycle; the higher the duty cycle, the greater transmittance of the pixel. In the following sections we introduce the block-unblock CASSI optical model and present the grayscale-adaptive-based model. Simulations are performed to evaluate the improvements attained by the proposed grayscale-adaptive coded apertures.

System model

The CASSI system uses a traditional block-unblock coded aperture (BCA) T(x, y) to modulate the spatio-spectral density f0(x, y, λ), where (x, y) are the spatial dimensions and λ represents the spectral dimension. The resulting coded field is then dispersed by the dispersive element, resulting in,

(1)

(a)

(b)

Figure 1. (a) Sketch of the grayscale-adaptive coded aperture-based CASSI system (GCA-CASSI). A PC real-time approach is used to provide feedback in order to update the coded aperture pixels involved in saturation, (b) Illustration of a GCA where three grayscale pixels are represented, each pixel have a distinct duty cycle, 25 %, 50 %, and 75 %. The higher the duty cycle, the more transmittance of the pixel.

where T (x´,y´) is the transmission function representing the coded aperture, h(x S (λ), y y´) is the optical impulse response of the system and S (λ ) is the dispersion induced by the prism. The compressive measurements across the FPA are obtained by the integration of the field f2 (x, y, λ) over the detector range sensitivity. Therefore, the spectral density in front of the detector is given by . When the signal measurement is replaced in g (x, y) the resulting equation is:

(2)

When the optical impulse response of the system is assumed linear and ideal, the spectral density being integrated is given by:

(3)

Representing the detector as a spatially pixelated array with pixel size Δd, the measurements in the presence of noise ω, at the (n, m)th position can be represented as:

(4)

where m, n represent the (m, n)th pixel in the FPA sensor, and rect(.) is the rectangular function. Similarly, the coded aperture T(x, y) can be represented as a spatially pixelated array. Assuming the coded aperture pixel size as Δt and represents a binary value, (0) block and (1) unblock, the coded aperture can be expressed as,

(5)

The coded aperture in Equation (5) can be replaced in Equation (4). Thus, the detector measurement becomes

(6)

Representing the spatio-spectral source density being integrated on the detector in discrete form as Fn,m,k such that indexes the x - axis, the y - axis, and the wavelength, Equation (6) can be succinctly expressed as

(7)

where Gn,m is the intensity at the (m,n)th position of the detector with dimensions N(N + L 1), the spectral data cube F has size N × N × L, Tn,m is the (m,n)th value in the coded aperture and ω represents the noise of the system.

Grayscale CASSI system model

In this paper the BCA is replaced with a GCA, which modulates the source along the spatial coordinates. The CASSI system architecture with GCA (GCA-CASSI) is illustrated in Figure 1a, where the traditional BCA is replaced by the GCA depicted in Figure 1b. The coding is now performed by the GCA represented bywhich is applied to the spatio-spectral density source f0(x, y, λ), resulting in the coded field f1(x, y, λ). This coded field differs from the one achieved with the BC, given that a particular element of the GCA attenuates the wavelengths instead of blocking or unblocking the complete spectrum at a given spatial location. The entries of the GCA vary in the range {0 ( l 1)}, being l the number of grayscale levels of the DMD. In this way, Equation (7) can be rewritten as,

(8)

FPA saturation in CASSI

Saturation occurs when the measurements exceed the dynamic range of the sensor quantizer. The quantizer has finite dynamic range due to two reasons: the first is related to physical limitations that allow a finite range voltage to be correctly converted to bits; and the second is that only a finite number of bits are available to represent each value. Quantization with saturation is referred to as finite-range quantization (Laska, et al, 2011). The errors imposed by finite-range quantization are unbounded. CS recovery techniques only provide guarantees for noise that is bounded, or bounded with high probability (Laska, et al, 2011). Dealing with saturation is important in CASSI because it reduces the attainable reconstruction quality. Figure 2 shows examples of compressive measurements, generated using Equation (7), three distinct percentages of saturated pixels (0 %, 5 % and 10 %, respectively), and their corresponding attained reconstructions using 4 shots. Notice that the higher the saturation percentage, the lower the quality of the reconstructed image.

In order to understand how the saturation in CASSI occurs, the qth slice of the data cube in Figure 3 is analyzed. It shows a representation of the physical phenomenon in CASSI, where t is the coded aperture of the qth row, is the qth the coded and dispersed slice, and g is the FPA measurement of the qth slice.

Figure 4 shows a top-view of the integration process of the FPA sensor. The pixel g6 is analyzed in order to understand the saturation in the sensor’s pixels. The qth slice of the data cube F is represented by the matrix F. Each fn,m element is pictorially represented as a small cube. The 6th compressive measurement in the qth row is given by, where each is integrated in the sensor and is given by , k = 0,...,L1. The intensity value in g6 only depends on the coded aperture pixels from t1 to t6. The saturation at g6 can be avoided by replacing the BCA by a GCA.

Figure 2. Compressive measurements for 3 levels of saturation and their respective reconstructions. First row: Compressive measurements with 0 %, 5 %, and 10 % FPA saturation. Second row: The corresponding reconstructions for 0 %, 5 %, and 10 % saturation using 4 shots.

Figure 3. Sketch of the spectral data flow in GCA-CASSI. The qth slice of the data cube with six spectral components is coded by the qth row of the GCA falta and dispersed by the prism. The sensor captures the intensity in g integrating the coded light.

Figure 4. The qth slice of the data cube F is represented by the matrix F. The elements are pictorially represented as a small cube. The intensity value in g6 only depends on the coded aperture pixels from t1 to t6.g6 saturation can be avoided by using a GCA.

Adaptive estimation of the GCA

In order to reduce the saturation level from the compressive measurements, an adaptive filter (AF) is designed to adaptively penalize the entries of the coded apertures, so the input source is attenuated before it is integrated by the detector. Consequently, these coded apertures will exhibit non-binary values, thus generating GCA. Formally, let Vi represent a weight matrix whose entries are a measure of how many times the ith coded aperture pixel affects a saturated pixel in the sensor. In particular, the entries of Vi can be written as

(9)

where u[.] is the unit step function. is the pixel from the ith compressive measurement. and K is the number of total snapshots. The number of shots can be expressed in terms of compressive ratio. The latter can be defined as , where K represents the number of shots, and N and L the spatial and spectral dimensions. Cr can be seen as the ratio between the reconstructed measurements and the number of pixels in the reconstructed data cube; s = 2(b-1) represents the saturation level of the sensor, which depends on the number of bits b of the sensor. Notice that Equation (9) can be easily calculated in a PC real time approach (Feedback). Based on the weight matrix, a penalization function is generated by assuming that the attenuation in a pixel of the coded aperture is inversely proportional to the number of saturated pixels in the FPA. That is, the penalization function can be seen as the attenuation matrix W i whose entries are given by,

(10)

where V0 is assumed to be an one-value matrix. Notice that the attenuation matrix W takes into account the previous weighted matrices as in a memory AF approach. Given the K randomly generated coded apertures T1,...,Tk, the GCA are generated according to,

(11)

where is the Hadamard product between the matrices A and B. Notice that , that is, the first GCA remains as the original, since the AF needs feedback. After the second snapshot some measurements still remain saturated. The effect of the AF improves the quality of reconstructed images when the remaining saturated measurements are discarded. The measurements that continue to be saturated in each snapshot are discarded, reducing the number of measurements and enhancing the quality of reconstructed images.

The adaptive estimation of GCA-CASSI improves the sensing. This gain in the detector is denominated “high dynamic range”. More specifically, high dynamic range is the characteristic of the GCA-CASSI that allows sensing areas with low and high irradiance. The low irradiance areas are correctly measured taking advantages of the inhomogeneous illumination. In addition, high irradiance areas are adequately measured due to the GCA attenuation after each snapshot. For that reason, the GCA-CASSI is a system with high dynamic range due to its ability to adapt to inhomogeneous illumination, sensing adequately the compressive measurements.

Spectral image reconstruction

CS theory establishes that a spectral signalor its vector representationis S sparse in some basis Ψ, such that f = Ψθ can be approximated by a linear combination of S vectors of Ψ with . Then, f can be reconstructed from d random projections with high probability when . In CASSI the compressive measurements can be represented in matrix form such that H is a matrix of size whose structure is determined by the coded apertures and the dispersive element. Similarly, multi-shot CASSI is represented as , where represents the effect of the coded aperture (Arguello and Arce, 2011). The set of K compressive measurements with a distinct coded aperture is then assembled as . The CASSI projections can be represented alternatively as , where the matrix is the sensing matrix. The reconstructed data cube is obtained as where , θ is a S sparse representation of f on the basis Ψ, and t is a constant of regularization. Let be the measurement vector with no saturated entries whose length is where . The matrix is created by adaptively selecting the non saturated rows of the matrix H. The reconstructed data cube is obtained as (Laska, et al, 2011).

Simulation results

In this section the CASSI with grayscale-adaptive coded apertures is compared against the traditional CASSI with BCA. A set of compressive measurements is simulated using the models in Equation (7) and Equation (11). The measurements were constructed using a test spectral database obtained using a wide-band Xenon lamp as the light source, and a visible monochromator, which spans the spectral range between 450nm and 650nm. The image intensity was captured using a CCD camera exhibiting 256 × 256 pixels. The simulations are performed in a desktop architecture with an Intel i7-4770 3.4 Ghz processor, 32 GB of RAM memory and using Matlab R2012b. The test data cube F with 256 × 256 pixels of spatial resolution and L = 16 spectral bands is shown in Figure 5, and Table 1 shows the bandwidth and the central wavelength of each spectral band.

The BCA entries are realizations of a Bernoulli random variable such that the transmittance of each pattern is 25 %. The entries of the GCA are random realizations of block, unblock and attenuation elements, such that the transmittance in the first shot is 25 % and it is updated in the following shots by the AF, following Equation (11). The number of saturated pixels in the measurements is varied from 0 % to 10 %. The coded apertures are designed to have 256 × 256 pixels of spatial resolution. The CS reconstruction is carried out using the algorithm Gradient Projection for Sparse Reconstruction (GPSR) (Figueiredo, et al, 2007). The basis representation Ψ is set as the Kronecker product of two basis where Ψ1 playing the role of spatial sparsifier is the 2D-Wavelet Symmlet 8 basis, and Ψ2 being the spectral sparsifier is the 1D-DCT basis. Figure 6 shows four snapshots using the BCA and the GCA. The silhouette of the compressive saturated measurements can be observed in the grayscale-adaptive coded aperture after the first snapshot. The resulting silhouette occurs when the weighted matrix attenuates the pixels in the coded aperture, which are responsible for saturated values in the compressive measurements.

Figure 5. Spectral datacube used for simulations.

Table 1. Bandwidth and central wavelength.

Bandwidth

Central wavelength

448–454

451

455-460

458

461–468

464

469-475

472

476–485

480

486-494

490

495–505

500

506-516

511

517-529

523

530–540

535

541-557

549

558–572

565

573–594

584

595-615

605

616–639

628

640-662

651


Figure 6. Comparison between BCA and GCA for four shots. The BCA and the GCA and the corresponding compressive measurements are shown.

Figure 7 shows the average reconstruction peak signal to noise ratio (PSNR) as a function of the percentage of saturation. The grayscale adaptive and the BCA are compared for two, four, six, and eight snapshots. In both cases noiseless and with noise tests were included. Noiseless simulations are important because they show the system in ideal conditions. In contrast, noise simulation illustrated the system in a real scenario, in which the measurements were affected with gaussian noise with SNR = 10 dB. GCA outperforms in up to 11 dB the BCA when the number of saturated measurements is increased.

Figure 7. Average reconstructions PSNR as a function of the percentage of saturation. The GCA are compared with BCA. The PSNR is measured for GCA and BCA for percentages of saturation between 0 % and 10 %. Clearly, the GCA improves the quality of the reconstructed images compared with the BCA.

Figure 8 shows the average PSNR of the reconstructed data cubes as a function of the number of snapshots. The GCA is compared with BCA for distinct number of snapshots and variation of FPA saturation. The improvement of the proposed optical architecture can be quantified by averaging the PSNR of the reconstructed datacubes. GCA outperform in up to 11 dBs the BCA when the number of snapshots is increased.

Figure 8. Average PSNR of the reconstructed data cubes as a function of measurement snapshots. The block-unblock-based and the grayscale adaptive-based CASSI Imagers are compared at distinct percentages of saturation. Top-left: Saturated measurements 1 % and 4 % with noise. Top-right: Saturated measurements 7 % and 10 % with noise. Bottom-left: Saturated measurements 1 % and 4 % noiseless. Bottom-right: Saturated measurements 7 % and 10 % noiseless. The use of GCA improves in quality of reconstructed images when the number of snapshots is higher compared with BCA.

Figure 9 shows a comparison between reconstructed images using GCA and BCA where the number of snapshots varies between K = {1,3,6,8}. The BCA and GCA reconstructed images are compared including white noise in the measurements where SNR = 10 dB. All reconstructions are obtained from FPA measurements with 10 % of saturated pixels. For 1 snapshot there are no significant differences in PSNR between the BCA and the GCA, since both architectures use BCA for the first snapshot. Clearly, when more than 1 snapshot is captured, the GCA outperforms the BCA in up to 11 dB in the quality of the reconstructed images. The higher the number of snapshots, the higher the quality of reconstructed images up to 30 dB.

Figure 10 shows a comparison between GCA and BCA where 4 snapshots were captured. This time, the FPA measurements were saturated with 0 %, 3 %, 6 % and 10 % of saturated pixels. The GCA outperforms the BCA in up to 11 dB in the quality of reconstructed images. The GCA tolerates increments in the percentage of saturation better than BCA. When the percentage of saturation is 0 % the quality of the reconstructed image is similar between GCA and BCA due to the fact that AF is inactive.

The resulting reconstructed data cubes curves are compared against their respective ground truth curves from the database. In addition, to evaluate the spectral reconstruction performance, 2 spatial points were randomly chosen and the spectral signatures are plotted in Figure 11 where 8 snapshots were captured and the number of saturated pixels in the FPA measurements is 5 % and 10 %. According to Figure 11, GCA spectral signatures are more similar to the original signature than those obtained with the BCA.

Figure 9. Comparison between reconstructed images using GCA and BCA for K = {1,3,6,8} snapshots. The first and third rows are reconstructed images from the FPA noiseless measurements. The second and fourth rows are reconstructed images from FPA measurements with noise, SNR = 10 dB. All reconstructions were obtained from FPA measurements with 10 % of saturated pixels. Clearly, when more than 1 snapshot is captured, the performance of GCA is higher than the BCA.

Figure 10. Comparison between GCA and BCA for K = 4 snapshots. First row: BCA noiseless. Second row: BCA with noise. Third row: GCA noiseless. Fourth row: GCA with noise. All reconstructions were obtained from FPA measurements with 0 %, 3 %, 6 % and 10 % of saturated pixels. When the percentage of saturation is increased, the GCA presents a higher performance.

Figure 11. Spectral signatures at four selected data points, K = 8 snapshots and the number of FPA saturated pixels is 5 % and 10 %. Top-left: P1 noiseless in the FPA measurements. Top-right: P2 noiseless in the FPA measurements. Bottom-left: P1 with noise in the FPA measurements. Bottom-right: P2 with noise in the FPA measurements. GCA spectral signature is more similar to the original signature than BCA.

Conclusions

GCA has been introduced in compressive spectral imaging system CASSI to replace the traditional BCA. The proposed architecture allows to attenuate the effects of the saturation of the FPA sensors. The saturation in CASSI depends on high irradiance that enters into the system overcoming the dynamic range of the sensor. The adaptive grayscale coded aperture deals with high irradiance attenuating it. The proposed technique has advantages in sensing the high irradiance but also the low irradiances, which means that the GCA-CASSI has a high dynamic range compared with BCA-CASSI. The designed GCA outperforms the BCA in up to 11 dB in the quality of the reconstructed images. The proposed adaptive filter can be implemented in PC real-time approach.

Acknowledgments

This work was supported by the Vicerrectoría de Investigación y Extensión of the Universidad Industrial de Santander, under the grant No. 1802.

References

Arguello, H. and Arce, G. R. (2011). Code aperture optimization for spectrally agile compressive imaging, J. Opt. Soc. Am. A, Vol. 28(11), 2400-2413, DOI: 10.1364/josaa.28.002400

Arguello, H., Arce, G. R. (2014). Colored coded aperture design by concentration of measurement in compressive spectral imaging. Image Processing, IEEE Transactions on, 23(4). 1896-1908. DOI: 10.1109/TIP.2014.2310125

Arguello, H., Arce, G. R. (2013). Rank minimization coded aperture design for spectrally selective compressive imaging. Image Processing, IEEE Transactions on, 22 (3), 941-954. DOI: 10.1109/TIP.2012.2222899

Baraniuk, R. (2007). Compressive sensing. IEEE Signal Processing Magazine, 24(4), 118-121. DOI: 10.1109/MSP.2007.4286571

Candes, E. J. (2007). Compressive sampling, Proceedings of the International Congress of Mathematics, 3, 1433-1452. DOI: 10.4171/022-3/69

Donoho, D. Compressed sensing. IEEE Transactions on Information Theory, 52 (4), 1289-1306. DOI: 10.1109/TIT.2006.871582

Figueiredo, M.A.T., Nowak, R.D., and Wright, S.J. (2007). Gradient Projection for Sparse Reconstruction: Application to Compressed Sensing and Other Inverse Problems, IEEE Journal of Selected Topics in Signal Processing, 4(1), 586-597. DOI: 10.1109/JSTSP.2007.910281

Kittle, D., Choi K., Wagadarikar A. and Brady, D. J. (2010). Multiframe image estimation for coded aperture snapshot spectral imagers, Appl. Opt., 49 (36), 6824-6833. DOI: 10.1364/AO.49.006824

Laska, J. N., Boufounos, P. T., Davenport, M. A., and Baraniuk, Richard. G. (2011). Democracy in action: Quantization, saturation, and compressive sensing, 31(3), 429-443.

Rueda, H, Arguello, H. (2013). Spatial super-resolution in coded aperture-based optical compressive hyperspectral imaging systems, Revista Facultad de Ingeniería Universidad de Antioquia, (67), 7-18.

Rueda, H, Parada, A, Arguello, H. (2014). Spectral resolution enhancement of hyperspectral imagery by a multiple-aperture compressive optical imaging system, Revista Ingeniería e investigación Universidad Nacional de Colombia, 34 (3), 50-55.

Wagadarikar, A. A., John, R., Willet, R., and Brady, D. (2008). Single disperser design for coded aperture snapshot spectral imaging, Appl. Opt., 47(10), B44-B51. DOI: 10.1364/AO.47.000B44

Wu, Y., Mirza, I. O., Arce, G. R. and Prather, D. W. (2011). Development of a digital-micromirror-device-based multishot snapshot spectral imaging system, Opt., Lett., 36(14), 2692-2694. DOI: 10.1364/OL.36.00269

How to Cite

APA

Diaz, N. E., Rueda Chacon, H. F. and Arguello Fuentes, H. (2015). High-dynamic range compressive spectral imaging by grayscale coded aperture adaptive filtering. Ingeniería e Investigación, 35(3), 53–60. https://doi.org/10.15446/ing.investig.v35n3.49868

ACM

[1]
Diaz, N.E., Rueda Chacon, H.F. and Arguello Fuentes, H. 2015. High-dynamic range compressive spectral imaging by grayscale coded aperture adaptive filtering. Ingeniería e Investigación. 35, 3 (Sep. 2015), 53–60. DOI:https://doi.org/10.15446/ing.investig.v35n3.49868.

ACS

(1)
Diaz, N. E.; Rueda Chacon, H. F.; Arguello Fuentes, H. High-dynamic range compressive spectral imaging by grayscale coded aperture adaptive filtering. Ing. Inv. 2015, 35, 53-60.

ABNT

DIAZ, N. E.; RUEDA CHACON, H. F.; ARGUELLO FUENTES, H. High-dynamic range compressive spectral imaging by grayscale coded aperture adaptive filtering. Ingeniería e Investigación, [S. l.], v. 35, n. 3, p. 53–60, 2015. DOI: 10.15446/ing.investig.v35n3.49868. Disponível em: https://revistas.unal.edu.co/index.php/ingeinv/article/view/49868. Acesso em: 14 nov. 2024.

Chicago

Diaz, Nelson Eduardo, Hoover Fabian Rueda Chacon, and Henry Arguello Fuentes. 2015. “High-dynamic range compressive spectral imaging by grayscale coded aperture adaptive filtering”. Ingeniería E Investigación 35 (3):53-60. https://doi.org/10.15446/ing.investig.v35n3.49868.

Harvard

Diaz, N. E., Rueda Chacon, H. F. and Arguello Fuentes, H. (2015) “High-dynamic range compressive spectral imaging by grayscale coded aperture adaptive filtering”, Ingeniería e Investigación, 35(3), pp. 53–60. doi: 10.15446/ing.investig.v35n3.49868.

IEEE

[1]
N. E. Diaz, H. F. Rueda Chacon, and H. Arguello Fuentes, “High-dynamic range compressive spectral imaging by grayscale coded aperture adaptive filtering”, Ing. Inv., vol. 35, no. 3, pp. 53–60, Sep. 2015.

MLA

Diaz, N. E., H. F. Rueda Chacon, and H. Arguello Fuentes. “High-dynamic range compressive spectral imaging by grayscale coded aperture adaptive filtering”. Ingeniería e Investigación, vol. 35, no. 3, Sept. 2015, pp. 53-60, doi:10.15446/ing.investig.v35n3.49868.

Turabian

Diaz, Nelson Eduardo, Hoover Fabian Rueda Chacon, and Henry Arguello Fuentes. “High-dynamic range compressive spectral imaging by grayscale coded aperture adaptive filtering”. Ingeniería e Investigación 35, no. 3 (September 1, 2015): 53–60. Accessed November 14, 2024. https://revistas.unal.edu.co/index.php/ingeinv/article/view/49868.

Vancouver

1.
Diaz NE, Rueda Chacon HF, Arguello Fuentes H. High-dynamic range compressive spectral imaging by grayscale coded aperture adaptive filtering. Ing. Inv. [Internet]. 2015 Sep. 1 [cited 2024 Nov. 14];35(3):53-60. Available from: https://revistas.unal.edu.co/index.php/ingeinv/article/view/49868

Download Citation

CrossRef Cited-by

CrossRef citations9

1. Jorge Bacca, Tatiana Gelvez-Barrera, Henry Arguello. (2021). Deep Coded Aperture Design: An End-to-End Approach for Computational Imaging Tasks. IEEE Transactions on Computational Imaging, 7, p.1148. https://doi.org/10.1109/TCI.2021.3122285.

2. Yingming Lai, Jinyang Liang. (2024). Coded Optical Imaging. , p.15. https://doi.org/10.1007/978-3-031-39062-3_2.

3. Xiaohan Li, Joel A. Greenberg, Michael E. Gehm. (2019). Single-shot multispectral imaging through a thin scatterer. Optica, 6(7), p.864. https://doi.org/10.1364/OPTICA.6.000864.

4. Nelson Diaz, Juan Ramirez, Esteban Vera, Henry Arguello. (2021). Adaptive Multisensor Acquisition via Spatial Contextual Information for Compressive Spectral Image Classification. IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing, 14, p.9254. https://doi.org/10.1109/JSTARS.2021.3111508.

5. Esteban Vera, Felipe Guzman, Nelson Diaz. (2024). Coded Optical Imaging. , p.499. https://doi.org/10.1007/978-3-031-39062-3_27.

6. David P. Bannon, Nelson Diaz, Hoover Rueda, Henry Arguello. (2016). Adaptive uniform grayscale coded aperture design for high dynamic range compressive spectral imaging. Hyperspectral Imaging Sensors: Innovative Applications and Sensor Standards 2016. Hyperspectral Imaging Sensors: Innovative Applications and Sensor Standards 2016. 9860, p.98600A. https://doi.org/10.1117/12.2224378.

7. Samuel Eduardo Pinilla, Héctor Miguel Vargas García, Henry Arguello Fuentes. (2016). Probability of correct reconstruction in compressive spectral imaging. Ingeniería e Investigación, 36(2), p.68. https://doi.org/10.15446/ing.investig.v36n2.56426.

8. Nelson Diaz, Carlos Hinojosa, Henry Arguello. (2019). Adaptive grayscale compressive spectral imaging using optimal blue noise coding patterns. Optics & Laser Technology, 117, p.147. https://doi.org/10.1016/j.optlastec.2019.03.038.

9. Nelson Diaz, Hoover Rueda, Henry Arguello. (2018). Adaptive filter design via a gradient thresholding algorithm for compressive spectral imaging. Applied Optics, 57(17), p.4890. https://doi.org/10.1364/AO.57.004890.

Dimensions

PlumX

Article abstract page views

941

Downloads

Download data is not yet available.

Most read articles by the same author(s)