Nothing Special   »   [go: up one dir, main page]
Next Article in Journal
Noise Removal in the Developing Process of Digital Negatives
Next Article in Special Issue
A Median-Ratio Scene-Based Non-Uniformity Correction Method for Airborne Infrared Point Target Detection System
Previous Article in Journal
Enrichment-Free Rapid Detection of Phthalates in Chinese Liquor with Electrochemical Impedance Spectroscopy
You seem to have javascript disabled. Please note that many of the page functionalities won't work as expected without javascript enabled.
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Sensors
Volume 20
Issue 3
10.3390/s20030900
sensors-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

DMD Mask Construction to Suppress Blocky Structural Artifacts for Medium Wave Infrared Focal Plane Array-Based Compressive Imaging

by
Zimu Wu
and
Xia Wang
*
Key Laboratory of Optoelectronic Imaging Technology and System, Ministry of Education, School of Optoelectronics, Beijing Institute of Technology, Beijing 100081, China
*
Author to whom correspondence should be addressed.
Sensors 2020, 20(3), 900; https://doi.org/10.3390/s20030900
Submission received: 30 December 2019 / Revised: 3 February 2020 / Accepted: 5 February 2020 / Published: 7 February 2020
(This article belongs to the Special Issue Infrared Sensors and Technologies: Recent Advances)
Browse Figures
Figure 1
<p>(<b>a</b>) Schematic diagram of our MWIR FPA CI system. The radiance of imaged scene is divided into two parts: one part to medium wave infrared (MWIR) sensor, and the other to cooling board, (<b>b</b>) Photograph of the actual system.</p> "> Figure 2
<p>Aperture interference in MWIR FPA CI system with the direct application of Hadamard matrix. The patterns of DMD mask for each block and the corresponding low-resolution MWIR images are shown in (<b>a</b>–<b>p</b>), the color-mapped MWIR images are used to represent the distributions of light and dark areas. In the patterns, white square stands for “1” and grey square for “0”.</p> "> Figure 3
<p>The high-resolution medium wave infrared (MWIR) images are reconstructed at (<b>a</b>) 0.125, (<b>b</b>) 0.25, (<b>c</b>) 0.375 and (<b>d</b>) 0.5. A small area of each image is zoomed in to highlight the blocky structural artefacts (BSA).</p> "> Figure 4
<p>The process of DMD mask generation. For the purpose of simplification, the block size <math display="inline"> <semantics> <mi>B</mi> </semantics> </math> is set to 4. In the patterns, white square stands for “1” and grey square for “0”.</p> "> Figure 5
<p>Imaging process of FPA CI. The images on DMD are the same in (<b>a</b>–<b>c</b>), but different DMD masks are used in the three imaging processes. The images on the sensor correspond to the low-resolution images, and the image on DMD is the high-resolution image, which is expected to be reconstructed.</p> "> Figure 6
<p>The numbers of the patterns for all blocks in DMD masks. The patterns of all blocks are generated from a 16-order Hadamard matrix. In (<b>a</b>–<b>p</b>), the block patterns are shown in the first row, and the numbers are depicted below the corresponding block patterns.</p> "> Figure 7
<p>The low-resolution MWIR images obtained from the actual MWIR FPA CI system. In (<b>a</b>–<b>p</b>), each subfigure corresponds to one proposed DMD mask. In the first row, the patterns of the big blocks are depicted on the left, and the captured MWIR images are shown on the right. The color-mapped MWIR images in the second row are used to represent the distributions of light and dark areas. The imaged scene is a uniform blackbody.</p> "> Figure 8
<p>The high-resolution MWIR images are reconstructed at (<b>a</b>) 0.125, (<b>b</b>) 0.25, (<b>c</b>) 0.375 and (<b>d</b>) 0.5. A small area of each image is zoomed in to highlight the BSA.</p> "> Figure 9
<p>The contrast of BRMSE between the two DMD mask design strategies.</p> "> Figure 10
<p>The high-resolution MWIR images are reconstructed with (<b>a</b>–<b>d</b>) the proposed DMD masks and (<b>e</b>–<b>h</b>) the original DMD masks. In each row, from left to right, the compression ratios are 0.125, 0.25, 0.375 and 0.5 respectively.</p> ">
Versions Notes

Abstract

:
With medium wave infrared (MWIR) focal plane array-based (FPA) compressive imaging (CI), high-resolution images can be obtained with a low-resolution MWIR sensor. However, restricted by the size of digital micro-mirror devices (DMD), aperture interference is inevitable. According to the system model of FPA CI, aperture interference aggravates the blocky structural artifacts (BSA) in the reconstructed images, which reduces the image quality. In this paper, we propose a novel DMD mask design strategy, which can effectively suppress BSA and maximize the reconstruction efficiency. Compared with random binary codes, the storage space and computation cost can be significantly reduced. Based on the actual MWIR FPA CI system, we demonstrate the proposed DMD masks can effectively suppress the BSA in the reconstructed images. In addition, a new evaluation index, blocky root mean square error, is proposed to indicate the BSA in FPA CI.

1. Introduction

Medium wave infrared (MWIR), of which the spectrum region is between 3 and 5 μm, has many imaging advantages beyond visible spectrum, such as night time observation and penetrating fog imaging [1]. Depending on the thermal radiance emitted and reflected by the objects, MWIR imaging systems can create images [2]. Targets are usually highlighted in MWIR images, as a result of temperature, emissivity and reflectivity difference between targets and background. Unfortunately, the MWIR sensors in the imaging systems are typically of low-resolution [3]. In order to improve performance of target recognition and surveillance, MWIR sensors of high-resolution are desired, but they are usually greater than tens of thousands dollars [4]. Hence, in spite of the immense application potential, high-resolution MWIR sensors are beyond the reach of scholars and researchers in some research fields which could get the most benefit from the high-resolution MWIR sensors.
MWIR compressive imaging (CI) is an alternative method of obtaining high-resolution MWIR images with a low-resolution MWIR sensor. The CI theory states that high-fidelity images can be reconstructed, which have a larger number of pixels than that physically exist in the imaging sensor [5,6,7]. An extreme example of CI is a single pixel camera (SPC) [8,9,10], which can reconstruct images with several hundred thousand pixels with just one photodetector. The trade-off of SPC, of course, is that a plenty of compressed samples of the imaging scene have to be obtained. Each compressed sample corresponds to a linear projection of the whole imaging scene onto one intensity value from a set of known functions. Although SPC has demonstrated potential for several applications, including remote sensing [11], 3D imaging [12], and microscopy [13,14], it is essentially a highly sequential imaging system. The compressed sampling process can be speeded up by using a focal plane array-based (FPA) sensor instead of the photodetector in SPC. Although the FPA CI systems have been implemented in visible [15,16] and short-wave infrared (SWIR) wavebands [17], comparatively few systems are implemented in MWIR. To our knowledge, only two experimental studies in the area of MWIR FPA CI are reported: Mahalanobis of Lockheed Martin built the first MWIR FPA CI system [18], which is a “merged-type”, and we built a “separate-type” system [19].
In spite of the acceleration of imaging speed, FPA CI has brought new challenges to the quality of reconstructed high-resolution images. As observed by Chen [17], blocky structural artifacts (BSA) are present in the reconstructed images, because of slight numerical fluctuations among the compressed samples for the same block. For MWIR FPA CI, due to physical limitation of the size of digital micro-mirror devices (DMD), aperture interference phenomena would also be observed [19], which aggravate the BSA and declines the image quality. An elaborate design for a DMD mask could alleviate the aperture interference in the low-resolution MWIR images obtained by the FPA sensor. In the study by Mahalanobis of Lockheed Martin [18], random binary codes are used for DMD masks. This strategy of DMD mask design can avoid the serious BSA, however, at a high cost in terms of storage resources and computing time. In addition, for some blocks, the mask patterns may be the same in two or more DMD masks, which reduces the efficiency of image compressed sampling and image reconstruction. An alternative common strategy is using Hadamard matrix for DMD mask design [20,21,22,23]. Because of aperture interference, the direct application of Hadamard matrix for DMD masks leads to obvious variations of light and dark areas distribution across the low-resolution images. As a result, the quality of the reconstructed high-resolution MWIR images will decline with an increase in the number of the low-resolution images.
In our previous work [19], the Hadamard matrix is directly applied to the DMD mask design, where all blocks have the same pattern in each DMD mask. Due to the aperture interference phenomenon, BSA has been observed in reconstructed high-resolution MWIR images. In this work, we further explore the aperture interference phenomenon in the MWIR FPA CI system, and focus our attention on DMD masks design to suppress the BSA. According to the system model of FPA CI, the reasons why BSA exist in the FPA CI are discussed. Particularly, the relationship between aperture interference and BSA is highlighted for MWIR FPA CI. As a result of these observations, based on the Hadamard matrix, we have designed customized patterns of DMD masks. Compared with a simple application of Hadamard matrix, the proposed patterns can suppress BSA, at the cost of a little more storage space and computation cost. On the other side, storage space and computation cost can be significantly reduced by the proposed patterns than random binary codes. We have modified the smoothed projected Landweber (SPL) reconstruction algorithm proposed by Mun [24], making it compatible with the proposed DMD mask. In addition, aiming at the BSA in FPA CI, we have proposed a new evaluation index, blocky root mean square error (BRMSE), to indicate the reconstructed image quality.
The rest of this paper is organized as follows. In Section 2, we describe the actual MWIR FPA CI system and present the aperture interference phenomenon when the DMD masks are constructed directly from Hadamard matrix. Section 3 proposes the customized patterns of DMD masks, modifies the SPL algorithm and describes BRMSE in detail. In Section 4, based on experimental results obtained from actual MWIR FPA CI system, we show the proposed patterns of DMD masks can suppress aperture interference, and BRMSE objectively assesses BSA. Some conclusions are given in Section 5.

2. Aperture Interference in MWIR FPA CI

A cursory overview of our MWIR FPA CI system is necessary to define the issues and terms for the subsequent discussion, which has been described in detail in [19]. The schematic diagram and actual photo of our MWIR FPA CI system are depicted in Figure 1. The imaged scene is projected onto the DMD by the imaging lens. DMD (DLP9500, Texas Instruments, Dallas, TX, USA) is an array of numerous tiny mirrors and the resolution is 1920 × 1080. Each tiny mirror can be controlled individually, rotating to either +12° or −12° around the normal vector of DMD. The DMD mask patterns are implemented on the DMD loaded by control circuit. As each tiny mirror has two distinct rotation angles, the DMD reflect one part of the MWIR irradiation to the relay lens and the other part to the cooling board. The cooling board is placed at the symmetrical position to the imaging lens, along the optical axis of relay lens. The cooling board absorbs the radiant heat, reducing the environmental impact. Through relay lens, the MWIR irradiation of the imaged scene, the part reflected by DMD, is projected onto the MWIR sensor. Then, the irradiation is captured by a cooled MWIR sensor (LEO MW, Sofradir, Palaiseau, France), of which the resolution is 640 × 512, and the corresponding low-resolution images are generated. In the actual laboratory setup, the effective size of DMD and MWIR sensor is 1280 × 1024 and 320 × 256, respectively, realizing a strict correspondence between a block of 4 × 4 micro-mirrors in DMD and a single detector of the MWIR sensor. Finally, low-resolution MWIR images are sent to a computer for high-resolution image reconstruction. In our experiments, a plane blackbody (CDS 100-4, EOI, Goleta, CA, USA) is used as an MWIR radiation source, of which the infrared radiation intensity can be thought as approximately parallel and uniform. A transmission type of United States Air Force (USAF) resolution test chart can be placed in front of the blackbody, acting as the testing criteria for MWIR FPA CI.
In the MWIR FPA CI system, due to the material of the imaging lens and relay lens, the imaging waveband is 3.7~4.7 μm. Compared with the size of the micro-mirror, which is 10.8 μm × 10.8 μm, the MWIR wavelength is in the same order of magnitude. According to the practical rule-of-thumb [25], Fraunhofer diffraction will occur at an aperture with the greatest width b , if
R > b 2 / λ ,
where λ is the imaging wavelength. R is the smaller of the two distances, where one is from radiation source to the aperture, and the other is from the aperture to the imaging plane. Therefore, for MWIR FPA CI, the aperture interference is inevitable.
If the DMD masks are directly constructed from 16-order Hadamard matrix [19], namely, the pattern of each block is identical in each DMD mask and the patterns of all the masks are derived from the rows in a Hadamard matrix by fulling elements in column-first order, the low-resolution MWIR images have different light and dark distributions when different DMD masks are applied into the system. If a micro-mirror is placed at +12°, it is opened and indicated as “1”; otherwise, the micro-mirror is placed at −12°, it is closed and indicated as “0”. Using the blackbody as imaged scene, the low-resolution MWIR images have distinct gray-level distributions, as depicted in Figure 2. To highlight the gray level distributions of these images, different colors are used in the colormap. The gray level ranges are (150, 240) for Figure 2a and (70, 150) for Figure 2b–p. It is obvious that different DMD mask patterns lead to distinct distributions of light and dark areas in the MWIR images. It is noteworthy that the gray level distribution in each image is caused by two reasons: ( cos ) 4 vignetting [26] and aperture interference. For the ( cos ) 4 vignetting, as the “opened” micro-mirrors are not in parallel with the DMD plane, the brightest region is not located at the center of the MWIR image. The aperture interference changes the distributions of pixel gray values across the low-resolution MWIR images.
If the new coordinate systems, X 1 O 1 Y 1 and X O Y , are established for DMD plane and detector plane, respectively, according to Fraunhofer diffraction [27], the radiance obtained by the MWIR sensor at the point ( x , y ) is
E ( x , y ) = C f exp [ i k ( f + x 2 + y 2 2 f ) ] ∬ Σ E ˜ ( x 1 , y 1 ) exp [ − i k f ( x x 1 + y y 1 ) ] d x 1 d y 1 ,
where C = 1 / i λ , k = 2 π / λ . f is the focal length of relay lens. E ˜ ( x 1 , y 1 ) is the radiance modulated by the DMD at ( x 1 , y 1 ) . The radiance of ( x 1 , y 1 ) is equal to the imaging scene when the micro-mirror is ‘1’, or zero when the micro-mirror is ‘0’.
From Equation (2), we can clearly conclude that the aperture interference phenomenon in low-resolution MWIR images is closely related to the DMD masks. In particular, when the opened micro-mirrors are aligned into parallel lines, such as 2b, 2c, 2d, 2e, 2i and 2m in Figure 2, multi-slit interference occurs.
For FPA CI, the gray values of all pixels in the reconstructed high-resolution image are computed from the low-resolution MWIR images. A FPA CI system can be viewed as an array of SPCs [28,29,30], thus the system model can be expressed in a matrix notation as
y i = Φ i x i ,
where, for the ith detector of MWIR FPA sensor, y i ∈ R m is the set of m compressed samples, x i ∈ R n represents the high-resolution image block, and Φ i ∈ R m × n is the measurement matrix for the block.
In the practical computation, y i is consisted of the ith pixel of all low-resolution MWIR images, and each row of Φ i is corresponding to the ith block pattern of the DMD mask. For example, if the DMD masks of Figure 2a–d are used for image reconstruction, Equation (3) can be rewritten as
[ y i ( 1 ) y i ( 2 ) y i ( 3 ) y i ( 4 ) ] = [ 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 1 0 0 1 1 0 0 1 1 0 0 1 1 0 0 1 0 0 1 1 0 0 1 1 0 0 1 1 0 0 1 ] [ x i , 1 x i , 2 ⋮ x i , 16 ] ,
where y i = [ y i ( 1 ) , y i ( 2 ) , y i ( 3 ) , y i ( 4 ) ] T and x i = [ x i , 1 , x i , 2 , ⋯ , x i , 16 ] T . In this work, to make the symbols clear and understandable, subscript is used to indicate the pixel position and superscript is for the number of the low-resolution MWIR image. The elements of measurement matrix and x i are consisted of DMD mask patterns and the ith high-resolution MWIR image block in column-first order, respectively.
Whatever image reconstruction algorithm is selected, the essence is solving the underdetermined equations, but the detailed process is different. If y i , 1 = 2 × y i , 2 = 2 × y i , 3 = 2 × y i , 4 is not satisfied, BSA will occur. For MWIR FPA CI, aperture interference brings a great gray value difference at the same pixel position when the number of opened micro-mirrors are identical, which can be observed in Figure 2b–p.
If high-resolution images are reconstructed from these low-resolution MWIR images, in which aperture interference phenomenon is obvious, the image quality is seriously influenced by BSA. In Figure 3, the reconstructed high-resolution images with compression ratios of 0.125, 0.25, 0.375 and 0.5 are shown, and the same areas are zoomed in to highlight the BSA.

3. Image Reconstruction and BRMSE

In traditional block-based compressed sensing algorithm, all blocks have the same pattern to simplify computation [24,31]. In our previous work [19], we applied this strategy to our MWIR FPA CI system, however, this led to serious aperture interference. Random binary codes for DMD masks could avoid the pattern repetition for all blocks [18], but tremendously increase storage burden and computing time. Moreover, in the strategy of random binary code, the patterns for the same block may be identical, which reduces the efficiency of compressive sampling as well. Thus, in this work, we also designed the DMD masks with the Hadamard matrix, but extended the circular pattern area, which slightly increased storage and computation expense, but remarkably suppressed the BSA in the reconstructed high-resolution images.
In this section, we first describe the generation process of the proposed DMD masks, which can suppress BSA in MWIR FPA CI. According to the characters of the proposed DMD masks, SPL [24] is also revised. Then, we discuss the storage space and computation cost of the proposed DMD mask design strategy. Depending on the theoretical analysis, the proposed DMD masks are proved to have dominant advantages over using random binary code. At last, based on the generation causes of the BSA in FPA CI, BRMSE is proposed to quantitatively evaluate the image quality of reconstructed high-resolution images.

3.1. DMD Masks

According to the basic theory of wave optics [32], similar DMD mask patterns lead to the similarity of aperture interference in the low-resolution MWIR images, which can also be observed in Figure 2. In Figure 2c,d, for the whole mask, the patterns are two-row horizontal stripes, but with one-row displacement. This trend can also be observed in the pairs of g-h, i-m, k-l and o-p. Based on this observation, we can conclude that, if the DMD masks for MWIR FPA CI have a similar pattern, the captured low-resolution MWIR images would have a similar gray level distribution.
Therefore, the philosophy of proposed DMD mask design strategy is: binding up a group of adjacent blocks to form a big block. For one DMD mask, all big blocks have the same pattern, but the blocks in each big block are distinct. In any pair of DMD masks, the patterns of the blocks are the same, but with the positional difference in big blocks. As a result, the proposed DMD mask patterns exhibit the similarity in the whole, but are diverse within the local area.
Assume that in a MWIR FPA CI system, the block size is 2 × 2 micro-mirrors, and a big block consists of 2 × 2 blocks. The Hadamard matrix is used to construct the block pattern, which has two advantages in MWIR CI. Firstly, the Hadamard matrix is full rank, maximizing the compressive sampling efficiency. Secondly, in each block, half of the micro-mirrors are “open”, implicitly improving the signal-to-noise ratio of the MWIR image. The generation process of DMD masks involves eight steps: (1) Generate a Hadamard matrix with the order equaling to the big block size B, which is 4 for this system. (2) Revise the Hadamard matrix to applicable for DMD, turning all “−1” to “0”. Using the complement vector of the second instead of the original one, all blocks have the same number of opened micro-mirrors, which may reduce the impact of pixel crosstalk [33]. (3) Divide the revised Hadamard matrix into separate row vectors. (4) For each row vector, a small pattern for each block is generated by fulling the elements of the row in column-first order. (5) Number the small patterns, from 1 to B. (6) In a big block of the size B × B, number the blocks in a clockwise direction. (7) In big blocks, replace the numbered block with the small patterns to generate the big patterns. (8) Expand the big patterns to coverage the entire DMD area and the DMD masks are obtained. This generation process of DMD masks is illustrated in Figure 4.

3.2. Revised SPL Algorithm

In the work of Mun [24], the measurement matrix is orthonormal and identical for all blocks. However, for our MWIR FPA CI system and proposed DMD masks, the measurement matrices for the big block are non-normalized and discrepant. Thus, the SPL algorithm need to be revised for the MWIR FPA CI system, and the revised SPL algorithm are summarized in Algorithm 1.
Algorithm 1 In the revised SPL algorithm, transform domain Ψ is chosen as dual-tree discrete wavelet transform (DDWT).
1: Input: compressed samples Y , measurement matrix Φ , Ψ
2: Output: X
3: Initialize: X ( 0 ) = B ( Y , Φ ) , i = 0
4: while | D ( i ) − D ( i − 1 ) | ≥ 10 − 4 do
5:    i = i + 1
6:    X ⌢ ( i ) = Wiener ( X ( i ) )
7:    X ⌢ ⌢ ( i ) = B ⌢ ( X ⌢ ( i ) , Φ )
8:    X ⌣ ⌣ ( i ) = Ψ ⋅ X ⌢ ⌢ ( i )
9:    X ⌣ ( i ) = Threshold ( X ⌣ ⌣ ( i ) , λ )
10:    X ¯ ( i ) = Ψ − 1 ⋅ X ⌣ ( i )
11:    X ( i + 1 ) = B ¯ ( X ¯ ( i ) , Φ )
12:    D ( i ) = 1 N ‖ X ( i + 1 ) − X ⌢ ⌢ ( i ) ‖ 2
13: end while
In the revised SPL algorithm, X ( 0 ) , X ⌢ ⌢ j ( i ) and X j ( i + 1 ) should be treated specially. For a block, its measurement matrix is determined by the position in the big block, so the functions of B ( ⋅ , ⋅ ) , B ⌢ ( ⋅ , ⋅ ) and B ¯ ( ⋅ , ⋅ ) are closely related to block number j . Algorithm 2 shows the definition B ( ⋅ , ⋅ ) , B ⌢ ( ⋅ , ⋅ ) and B ¯ ( ⋅ , ⋅ ) .
Algorithm 2 In the revised SPL algorithm, B ( ⋅ , ⋅ ) , B ⌢ ( ⋅ , ⋅ ) and B ¯ ( ⋅ , ⋅ ) are all block-wise, thus they are bounded together to clarify.
1: Require: Y = [ y 1 , y 2 , ⋯ , y N ] and X = [ x 1 , x 2 , ⋯ , x N ] are block-wise
2: Require: measurement matrix Φ for the first block in big block
3: for block j = 1 to N do
4:    i = b l o c k p o s i t i o n ( j ) − 1
5:    Φ j = circshiftup ( Φ , i )
6:    B ( Y , Φ ) ⇒ x j ( 0 ) = Φ j T · y j
7:    B ⌢ ( X ⌢ ( i ) , Φ ) ⇒ x ⌢ ⌢ j ( i ) = x ⌢ j ( i ) + Φ j T ⋅ ( Φ j ⋅ Φ j T ) − 1 ⋅ ( y j − Φ j ⋅ x ⌢ j ( i ) )
8:    B ¯ ( X ¯ ( i ) , Φ ) ⇒ x j ( i + 1 ) = x ¯ j ( i ) + Φ j T ⋅ ( Φ j ⋅ Φ j T ) − 1 ⋅ ( y j − Φ j ⋅ x ¯ j ( i ) )
9: end for
Where the symbol ⇒ represents the function definitions for B ( ⋅ , ⋅ ) , B ⌢ ( ⋅ , ⋅ ) and B ¯ ( ⋅ , ⋅ ) in the jth block. b l o c k p o s i t i o n ( ⋅ ) returns the position of the jth block in a big block. In the proposed DMD masks, blocks in a big block are numbered in row-first order, and all big blocks are done in the same manner. circshiftup ( Φ , i ) circularly shifts the rows of Φ by i rows upward.
In addition, ( Φ ⋅ Φ T ) − 1 is added into the calculation of B ⌢ ( ⋅ , ⋅ ) and B ¯ ( ⋅ , ⋅ ) , to eliminate the influence of the non-normalization of measurement matrix for the FPA CI system. If ( Φ ⋅ Φ T ) − 1 is not added, all of the image information would be lost.

3.3. Advantages of the Proposed DMD Masks

To illustrate the difference of storage space and computation cost among the original DMD masks, the proposed ones and the random binary codes in FPA CI, we use the system models of the three to analyze the computing costs. The imaging process of FPA CI is shown in Figure 5. The image on DMD is the high-resolution image expected to be reconstructed, and the image on sensor is the low-resolution image. In FPA CI, the high-resolution image is computed from multiple times of low-resolution image sampling. To simplify the analysis, the resolutions of DMD and FPA sensor are set to 4 × 4 and 2 × 2 respectively. Figure 5 shows the twice low-resolution image sampling process with different strategies of DMD mask design.
In the original strategy of DMD mask design [19], as shown in Figure 5a, all blocks have the same pattern. According to Equation (3), for all pairs of DMD block and detector, the system model can be expressed as
{ [ y 1 , 1 ( 1 ) y 1 , 1 ( 2 ) ] = [ 1 1 1 1 1 1 1 0 ] [ x 1 , 1 x 2 , 1 x 1 , 2 x 2 , 2 ] [ y 2 , 1 ( 1 ) y 2 , 1 ( 2 ) ] = [ 1 1 1 1 1 1 1 0 ] [ x 3 , 1 x 4 , 1 x 3 , 2 x 4 , 2 ] [ y 1 , 2 ( 1 ) y 1 , 2 ( 2 ) ] = [ 1 1 1 1 1 1 1 0 ] [ x 1 , 3 x 2 , 3 x 1 , 4 x 2 , 4 ] [ y 2 , 2 ( 1 ) y 2 , 2 ( 2 ) ] = [ 1 1 1 1 1 1 1 0 ] [ x 3 , 3 x 4 , 3 x 3 , 4 x 4 , 4 ] ,
As the measurement matrices are same, Equation (5) can be integrated into a single expression, such as
[ y 1 , 1 ( 1 ) y 1 , 1 ( 2 ) y 2 , 1 ( 1 ) y 2 , 1 ( 2 ) y 1 , 2 ( 1 ) y 1 , 2 ( 2 ) y 2 , 2 ( 1 ) y 2 , 2 ( 2 ) ] = [ 1 1 1 1 1 1 1 0 ] [ x 1 , 1 x 2 , 1 x 1 , 2 x 2 , 2 x 3 , 1 x 4 , 1 x 3 , 2 x 4 , 2 x 1 , 3 x 2 , 3 x 1 , 4 x 2 , 4 x 3 , 3 x 4 , 3 x 3 , 4 x 4 , 4 ] ,
Analogously, in the proposed strategy of DMD mask design, a big block consists of 1 × 2 blocks, and all big blocks have the same pattern, which is drawn in yellow dashed lines, as shown in Figure 5b. The system model can be expressed as
{ [ y ¯ 1 , 1 ( 1 ) y ¯ 1 , 1 ( 2 ) y ¯ 2 , 1 ( 1 ) y ¯ 2 , 1 ( 2 ) ] = [ 1 0 1 0 1 0 0 1 ] [ x 1 , 1 x 2 , 1 x 1 , 2 x 2 , 2 x 3 , 1 x 4 , 1 x 3 , 2 x 4 , 2 ] [ y ¯ 1 , 2 ( 1 ) y ¯ 1 , 2 ( 2 ) y ¯ 2 , 2 ( 1 ) y ¯ 2 , 2 ( 2 ) ] = [ 1 0 0 1 1 0 1 0 ] [ x 1 , 3 x 2 , 3 x 1 , 4 x 2 , 4 x 3 , 3 x 4 , 3 x 3 , 4 x 4 , 4 ] ,
In Equation (7), two expressions are used for the system model, because the big block consists of two blocks. For the blocks in one big block, the measurement matrices are distinct, so they should be treated separately. However, for the blocks at the same position of all big blocks, the measurement matrices are the same, and they can be integrated into one expression. This behavior is also reflected in Equation (7).
In the strategy of random binary codes, as the measurement matrices are distinct for all blocks, as shown in Figure 5c, all pairs of DMD blocks and detectors must be treated separately, and the system model can be expressed as
{ [ y ˜ 1 , 1 ( 1 ) y ˜ 1 , 1 ( 2 ) ] = [ 1 0 0 0 0 0 0 0 ] [ x 1 , 1 x 2 , 1 x 1 , 2 x 2 , 2 ] [ y ˜ 2 , 1 ( 1 ) y ˜ 2 , 1 ( 2 ) ] = [ 1 0 1 1 0 1 0 0 ] [ x 3 , 1 x 4 , 1 x 3 , 2 x 4 , 2 ] [ y ˜ 1 , 2 ( 1 ) y ˜ 1 , 2 ( 2 ) ] = [ 0 0 1 1 0 1 0 1 ] [ x 1 , 3 x 2 , 3 x 1 , 4 x 2 , 4 ] [ y ˜ 2 , 2 ( 1 ) y ˜ 2 , 2 ( 2 ) ] = [ 1 1 1 1 1 1 1 1 ] [ x 3 , 3 x 4 , 3 x 3 , 4 x 4 , 4 ] ,
It is noteworthy that, for the last expression of Equation (8), the rows of the measurement matrix are the same. As shown in Figure 5c, in the twice low-resolution image sampling, the DMD masks are identical for one block, which is drawn in red dashed lines. Theoretically, y ˜ 2 , 2 ( 1 ) is equal to y ˜ 2 , 2 ( 2 ) , which means only one sampling is efficient for this block. The repeating pattern for some blocks, which is unavoidable for random binary codes, reduces the efficiency of image compressed sampling and image reconstruction.
In the FPA CI, the computing costs of image reconstruction is closely related to the number of matrix computation. In parallel computation, one thread is required for one matrix computation, and too many threads would dramatically increase computation burden, which is also a nightmare for software developers [34]. For the original strategy, only one matrix computation is needed. For the proposed strategy, the number of matrix computation equals to the big block size. However, for the random binary codes, the number of matrix computation equals to the number of detectors on the FPA sensor.
For larger-area FPA sensor, the gaps increases for the number of matrix computation among the three strategies. For our MWIR FPA CI system, the block size is 4 × 4 micro-mirrors, and a big block consists of 4 × 4 blocks. Thus, the number of matrix computation is 16 for the proposed strategy; that is 81,920 for random binary codes. In parallel computation, the number of threads is 1 for the original strategy, 16 for the proposed strategy, and 81,920 for random binary codes. In today’s multicore processors, 16 threads can be easily obtained, while 81,920 is far beyond reach. It is evident, compared with random binary codes, that the proposed strategy has an overwhelming advantage for computation cost. In addition, less expressions in the system model means less measurement matrices needed be stored in the computer memory; thus, storage space can also be saved in the proposed strategy to a greater extent than in using random binary codes. Consequently, based on the consideration of the algorithm complexity and real-time property of the MWIR FPA CI system, random binary codes may only be a good idea on paper.

3.4. BRMSE

In FPA CI, the BSA appear in the reconstructed high-resolution images, which reduces the image quality. Because of slight fluctuations in gray-level values of the same pixel among the low-resolution images, BSA is inevitable [17]. Aperture Interference, observed in MWIR, would aggravate the BSA. In the field of image sparse representation, the quality of the reconstructed image is generally evaluated by peak signal-to-noise ratio (PSNR). However, in FPA CI, the original high-resolution image is unobtainable, so PSNR cannot be used.
From the system perspective, each block can be treated independently, so BSA can be evaluated locally in every block. BRMSE is the mean value of the root mean square errors (RMSE) for all blocks, which can be expressed as:
B R M S E = 1 N ∑ i = 1 N R M S E ( b l o c k i )
where N is the number of blocks in one reconstructed high-resolution image, R M S E ( ⋅ ) calculates the RMSE of each block, which is defined as
R M S E ( b l o c k ) = 1 n ∑ i = 1 n ( g i − g ¯ ) 2
where R M S E ( b l o c k ) is the RMSE of the block, b l o c k , which has n pixels. g i is the gray value of the ith pixel and g ¯ the average gray value of b l o c k .
In order to evaluate BSA exactly, a uniform imaging scene—as in those of Figure 2 and Figure 3—is preferred for the calculation of BRMSE. Obviously, reconstructed images with less BSA have smaller BRMSE.

4. Experiment and Discussion

For our MWIR FPA CI system, the block size is 4 × 4 micro-mirrors, and a big block consists of 4 × 4 blocks. According to the proposed strategy of DMD mask design, a 16-order Hadamard matrix is used to generate patterns for all the blocks, and the numbers of the patterns are depicted in Figure 6. Big blocks are consisted of these blocks, and a set of low-resolution MWIR images are obtained from the actual MWIR FPA CI system, as shown in Figure 7. The big block patterns are shown on the left in the first row, of which the size is 4 × 4 blocks. To highlight the gray level distributions of these images, color-mapped images are placed in the second row, and the gray level ranges are (85, 140). With the blackbody used as imaged scene, the images have almost identical light and dark distributions, which is quite different from that in Figure 2. Therefore, compared with a simple application of the Hadamard matrix, our proposed DMD masks can effectively suppress BSA in MWIR FPA CI, visually.
If high-resolution images are reconstructed from these low-resolution MWIR images, in which aperture interference phenomenon is not obvious, BSA can be alleviated. In Figure 8, the reconstructed high-resolution images with compression ratios of 0.125, 0.25, 0.375 and 0.5 are shown, and the same areas are zoomed in to highlight the BSA.
To objectively describe the effect of the proposed DMD masks, we use BRMSE to evaluate the BSA of the reconstructed high-resolution images. With the two strategies of DMD mask design, one is the original [19] and the other is the proposed, high-resolution MWIR images, using the blackbody as imaged scene, reconstructed with different amounts of low-resolution images, from 1 to 16. For our MWIR FPA CI system, 16 is the maximum sampling number, equal to the number of micro-mirrors in one block. The contrast of BRMSE for the two DMD mask design strategies is shown in Figure 9. It is obvious that the proposed DMD masks can effectively reduce the BRMSE of the reconstructed high-resolution images, which means the BSA caused by aperture interference in MWIR FPA CI is effectively suppressed.
Using a transmission type of USAF resolution test chart as imaged scene, low-resolution MWIR images are obtained with the two strategies of DMD mask design. Then, the reconstructed high-resolution images are computed with compression ratios of 0.125, 0.25, 0.375 and 0.5. To highlight the BSA in the reconstructed images, some areas are zoomed in, which are shown in Figure 10.
In the first row, images are reconstructed using the low-resolution MWIR images obtained with the proposed DMD masks, at compression ratios of 0.125, 0.25, 0.375 and 0.5. In the second row, images are reconstructed using the images obtained with the original DMD masks. In the flat zoomed-in zones, which are highlighted by yellow borders, the proposed DMD masks are obviously superior to the original ones. For the textured zoomed-in zones, which are highlighted by red borders, the image quality of the original DMD masks is significantly reduced in Figure 10e, but looks sharper in f–h than a–d, which were obtained using the proposed DMD masks. The reason for this phenomenon is, for the original DMD masks, due to the BSA, the gray values are unevenly assigned to every two adjacent row, making one row bright and the other dark. In viewing the overall image, the image may appear shaper, but the local features of the image are seriously damaged, which is a great obstacle for traditional digital image processing algorithms [35]. However, for the proposed DMD masks, the local features of the image are better preserved, and the line pairs are not cut into pieces as with the original DMD masks. This improvement is helpful for automatic target recognition [36]. Although the BSA cannot be completely eliminated, image qualities are improved by the proposed DMD masks in comparison to the original ones.

5. Conclusions

In MWIR FPA CI, aperture interference phenomena are inevitable, which aggravate the BSA in the reconstructed high-resolution images. The quality of the reconstructed images is getting worse as the compression ratio increases. As a result, the aperture interference reduces the practicality of the MWIR FPA CI system. In this paper, we analyze the relationship between the aperture interference and BSA, and propose a novel strategy of DMD mask design, which can effectively suppress the BSA. Compared with a simple application of the Hadamard matrix, the proposed DMD masks that can alleviate the BSA of the high-resolution images reconstructed with MWIR, at the cost of a little more storage space and computation cost. Compared with random binary codes, the proposed DMD masks can significantly reduce storage space and computation cost, and maximize the reconstruction efficiency, as the patterns for each block are completely distinct between the DMD masks. To objectively evaluate the BSA in FPA CI, we propose a new evaluation index, BRMSE. With the actual MWIR FPA CI system, we demonstrate that the proposed DMD masks can significantly reduce the BRMSE, which means the BSA is effectively suppressed. The proposed DMD mask design strategy improves the practicability of MWIR FPA CI.
In FPA CI, aperture interference is closely related to the imaging wavelength. As another available waveband, long wave infrared (LWIR) is also important for infrared imaging and has a high research value. In future work, we will investigate the aperture interference in LWIR FPA CI.

Author Contributions

Z.W. conceived and designed the research; Z.W. and X.W. performed the research. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Coffey, V.C. Seeing in the Dark: Defense Applications of IR imaging. Opt. Photonics News 2011, 4, 26–31. [Google Scholar] [CrossRef]
  2. Ring, E.; Ammer, K. Infrared thermal imaging in medicine. Physiol. Meas. 2012, 3, 33–46. [Google Scholar] [CrossRef] [PubMed]
  3. Millikan, B.; Dutta, A.; Sun, Q.; Foroosh, H. Fast detection of compressively sensed IR targets using stochastically trained least squares and compressed quadratic correlation filters. IEEE Trans. Aerosp. Electron. Syst. 2017, 53, 2449–2461. [Google Scholar] [CrossRef]
  4. Huang, X.; Netravali, R.; Man, H.; Lawrence, V. Multi-sensor fusion of infrared and electro-optic signals for high resolution night images. Sensors 2012, 8, 10326–10338. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Donoho, D.L. Compressed sensing. IEEE Trans. Inf. Theory. 2006, 52, 1289–1306. [Google Scholar] [CrossRef]
  6. Wakin, M.B.; Laska, J.N.; Duarte, M.F.; Baron, D.; Sarvotham, S.; Takhar, D.; Kelly, K.F.; Baraniuk, R.G. An architecture for compressive imaging. In Proceedings of the International Conference on Image Processing, Atlanta, GA, USA, 8–11 October 2006; pp. 1273–1276. [Google Scholar]
  7. Baraniuk, R.G. Compressive sensing. IEEE Signal Process Mag. 2007, 4, 118–121. [Google Scholar] [CrossRef]
  8. Duarte, M.F.; Davenport, M.A.; Takhar, D.; Laska, J.N.; Sun, T.; Kelly, K.F.; Baraniuk, R.G. Single-pixel imaging via compressive sampling. IEEE Signal Process Mag. 2008, 2, 83–91. [Google Scholar] [CrossRef] [Green Version]
  9. Radwell, N.; Mitchell, K.J.; Gibson, G.M.; Edgar, M.P.; Bowman, R.; Padgett, M.J. Single-pixel infrared and visible microscope. Optica 2014, 5, 285–289. [Google Scholar] [CrossRef]
  10. Gibson, G.M.; Sun, B.; Edgar, M.P.; Phillips, D.B.; Hempler, N.; Maker, G.T.; Malcolm, G.P.; Padgett, M.J. Real-time imaging of methane gas leaks using a single-pixel camera. Opt. Express 2017, 4, 2998–3005. [Google Scholar] [CrossRef]
  11. Willett, R.M.; Marcia, R.F.; Nichols, J.M. Compressed sensing for practical optical imaging systems: a tutorial. Opt. Eng. 2011, 7, 72601. [Google Scholar]
  12. Sun, B.; Edgar, M.P.; Bowman, R.; Vittert, L.E.; Welsh, S.; Bowman, A.; Padgett, M.J. 3D computational imaging with single-pixel detectors. Science 2013, 340, 844–847. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Wu, Y.; Ye, P.; Mirza, I.O.; Arce, G.R.; Prather, D.W. Experimental demonstration of an optical-sectioning compressive sensing microscope (CSM). Opt. Express 2010, 24, 24565–24578. [Google Scholar] [CrossRef] [PubMed]
  14. Studer, V.; Bobin, J.; Chahid, M.; Mousavi, H.S.; Candes, E.; Dahan, M. Compressive fluorescence microscopy for biological and hyperspectral imaging. In Proceedings of the Imaging Systems and Applications 2012, Monterey, CA, USA, 24–28 June 2012; pp. 1679–1687. [Google Scholar]
  15. Ke, J.; Lam, E.Y. Object reconstruction in block-based compressive imaging. Opt. Express 2012, 20, 22102–22117. [Google Scholar] [CrossRef] [PubMed]
  16. Shepard, R.H.; Fernandez-Cull, C.; Raskar, R.; Shi, B.; Barsi, C.; Zhao, H. Optical design and characterization of an advanced computational imaging system. In Proceedings of the Optics and Photonics for Information Processing VIII, San Diego, CA, USA, 17–21 August 2014. [Google Scholar]
  17. Chen, H.; Asif, M.S.; Sankaranarayanan, A.C.; Veeraraghavan, A. FPA-CS: Focal Plane Array-based Compressive Imaging in Short-wave Infrared. In Proceedings of the IEEE Conference on Computer Vision and Pattern Recognition, Boston, MA, USA, 7–12 June 2015. [Google Scholar]
  18. Mahalanobis, A.; Shilling, R.; Murphy, R.; Muise, R. Recent results of medium wave infrared compressive sensing. Appl. Opt. 2014, 34, 8060–8070. [Google Scholar] [CrossRef] [PubMed]
  19. Wu, Z.; Wang, X. Focal plane array-based compressive imaging in medium wave infrared: Modeling, implementation, and challenges. Appl. Opt. 2019, 31, 8433–8441. [Google Scholar] [CrossRef] [PubMed]
  20. Huang, G.; Jiang, H.; Matthews, K.; Wilford, P. Lensless imaging by compressive sensing. In Proceedings of the IEEE International Conference on Image Processing, Melbourne, Australia, 15–18 September 2013. [Google Scholar]
  21. Shishkin, S.L. Fast and robust compressive sensing method using mixed Hadamard sensing matrix. IEEE J. Emerg. Sel. Top. Circ. Syst. 2012, 3, 353–361. [Google Scholar] [CrossRef]
  22. Cai, Z.; Zhao, H.; Jia, M.; Wang, G.; Shen, J. An improved Hadamard measurement matrix based on Walsh code for compressive sensing. In Proceedings of the 9th International Conference on Information, Communications & Signal Processing, Tainan, Taiwan, 10–13 December 2013. [Google Scholar]
  23. Salari, S.; Chan, F.; Chan, Y.; Read, W. TDOA estimation with compressive sensing measurements and hadamard Matrix. IEEE Trans. Aerosp. Electron. Syst. 2018, 6, 3137–3142. [Google Scholar] [CrossRef]
  24. Mun, S.; Fowler, J.E. Block compressed sensing of images using directional transforms. In Proceedings of the International Conference on Image Processing, Cairo, Egypt, 7–10 November 2009. [Google Scholar]
  25. Hecht, E. Optics, 5th ed.; Pearson Education Press: Essex, UK, 2017. [Google Scholar]
  26. Zheng, Y.; Lin, S.; Kambhamettu, C.; Yu, J.; Kang, S.B. Single-image vignetting correction. IEEE Trans. Pattern Anal. 2008, 12, 2243–2256. [Google Scholar]
  27. Pedrotti, F.L.; Pedrotti, L.M.; Pedrotti, L.S. Introduction to Optics; Cambridge University Press: Cambridge, UK, 2017. [Google Scholar]
  28. Sun, M.; Chen, W.; Liu, T.; Li, L. Image retrieval in spatial and temporal domains with a quadrant detector. IEEE Photonics J. 2017, 5, 1–6. [Google Scholar] [CrossRef]
  29. Sun, M.; Wang, H.; Huang, J. Improving the performance of computational ghost imaging by using a quadrant detector and digital micro-scanning. Sci. Rep. 2019, 9, 4105. [Google Scholar] [CrossRef] [Green Version]
  30. Herman, M.A.; Tidman, J.; Hewitt, D.; Weston, T.; McMackin, L. A higher-speed compressive sensing camera through multi-diode design. In Proceedings of the Compressive Sensing II, Baltimore, MD, USA, 29 April–3 May 2013. [Google Scholar]
  31. Unde, A.S.; Deepthi, P.P. Block compressive sensing: Individual and joint reconstruction of correlated images. J. Vis. Commun. Image Represent. 2017, 44, 187–197. [Google Scholar] [CrossRef]
  32. Gupta, S.D.; Ghosh, N.; Banerjee, A. Wave Optics: Basic Concepts and Contemporary Trends; CRC Press: Boca Raton, FL, USA, 2015. [Google Scholar]
  33. Dumas, J.P.; Lodhi, M.A.; Bajwa, W.U.; Pierce, M.C. Computational imaging with a highly parallel image-plane-coded architecture: Challenges and solutions. Opt. Express 2016, 6, 6145–6155. [Google Scholar] [CrossRef] [PubMed]
  34. Patterson, D.A.; Hennessy, J.L. Computer Organization and Design MIPS Edition: The Hardware/Software Interface; Morgan kaufmann: Cambridge, UK, 2013. [Google Scholar]
  35. Szeliski, R. Computer Vision: Agorithms and Applications; Springer Science & Business Media: Berlin, Germany, 2010. [Google Scholar]
  36. Wu, Z.; Wang, X.; Jiao, Y.; Zhu, Y.; Zhou, J. Guidance performance evaluation method for infrared imaging guided missile based on extended object-oriented Petri net. Optik 2019, 185, 88–96. [Google Scholar] [CrossRef]
Sensors 20 00900 g001 550
Figure 1. (a) Schematic diagram of our MWIR FPA CI system. The radiance of imaged scene is divided into two parts: one part to medium wave infrared (MWIR) sensor, and the other to cooling board, (b) Photograph of the actual system.
Figure 1. (a) Schematic diagram of our MWIR FPA CI system. The radiance of imaged scene is divided into two parts: one part to medium wave infrared (MWIR) sensor, and the other to cooling board, (b) Photograph of the actual system.
Sensors 20 00900 g001
Sensors 20 00900 g002 550
Figure 2. Aperture interference in MWIR FPA CI system with the direct application of Hadamard matrix. The patterns of DMD mask for each block and the corresponding low-resolution MWIR images are shown in (a–p), the color-mapped MWIR images are used to represent the distributions of light and dark areas. In the patterns, white square stands for “1” and grey square for “0”.
Figure 2. Aperture interference in MWIR FPA CI system with the direct application of Hadamard matrix. The patterns of DMD mask for each block and the corresponding low-resolution MWIR images are shown in (a–p), the color-mapped MWIR images are used to represent the distributions of light and dark areas. In the patterns, white square stands for “1” and grey square for “0”.
Sensors 20 00900 g002
Sensors 20 00900 g003 550
Figure 3. The high-resolution medium wave infrared (MWIR) images are reconstructed at (a) 0.125, (b) 0.25, (c) 0.375 and (d) 0.5. A small area of each image is zoomed in to highlight the blocky structural artefacts (BSA).
Figure 3. The high-resolution medium wave infrared (MWIR) images are reconstructed at (a) 0.125, (b) 0.25, (c) 0.375 and (d) 0.5. A small area of each image is zoomed in to highlight the blocky structural artefacts (BSA).
Sensors 20 00900 g003
Sensors 20 00900 g004 550
Figure 4. The process of DMD mask generation. For the purpose of simplification, the block size B is set to 4. In the patterns, white square stands for “1” and grey square for “0”.
Figure 4. The process of DMD mask generation. For the purpose of simplification, the block size B is set to 4. In the patterns, white square stands for “1” and grey square for “0”.
Sensors 20 00900 g004
Sensors 20 00900 g005 550
Figure 5. Imaging process of FPA CI. The images on DMD are the same in (a–c), but different DMD masks are used in the three imaging processes. The images on the sensor correspond to the low-resolution images, and the image on DMD is the high-resolution image, which is expected to be reconstructed.
Figure 5. Imaging process of FPA CI. The images on DMD are the same in (a–c), but different DMD masks are used in the three imaging processes. The images on the sensor correspond to the low-resolution images, and the image on DMD is the high-resolution image, which is expected to be reconstructed.
Sensors 20 00900 g005
Sensors 20 00900 g006 550
Figure 6. The numbers of the patterns for all blocks in DMD masks. The patterns of all blocks are generated from a 16-order Hadamard matrix. In (a–p), the block patterns are shown in the first row, and the numbers are depicted below the corresponding block patterns.
Figure 6. The numbers of the patterns for all blocks in DMD masks. The patterns of all blocks are generated from a 16-order Hadamard matrix. In (a–p), the block patterns are shown in the first row, and the numbers are depicted below the corresponding block patterns.
Sensors 20 00900 g006
Sensors 20 00900 g007 550
Figure 7. The low-resolution MWIR images obtained from the actual MWIR FPA CI system. In (a–p), each subfigure corresponds to one proposed DMD mask. In the first row, the patterns of the big blocks are depicted on the left, and the captured MWIR images are shown on the right. The color-mapped MWIR images in the second row are used to represent the distributions of light and dark areas. The imaged scene is a uniform blackbody.
Figure 7. The low-resolution MWIR images obtained from the actual MWIR FPA CI system. In (a–p), each subfigure corresponds to one proposed DMD mask. In the first row, the patterns of the big blocks are depicted on the left, and the captured MWIR images are shown on the right. The color-mapped MWIR images in the second row are used to represent the distributions of light and dark areas. The imaged scene is a uniform blackbody.
Sensors 20 00900 g007
Sensors 20 00900 g008 550
Figure 8. The high-resolution MWIR images are reconstructed at (a) 0.125, (b) 0.25, (c) 0.375 and (d) 0.5. A small area of each image is zoomed in to highlight the BSA.
Figure 8. The high-resolution MWIR images are reconstructed at (a) 0.125, (b) 0.25, (c) 0.375 and (d) 0.5. A small area of each image is zoomed in to highlight the BSA.
Sensors 20 00900 g008
Sensors 20 00900 g009 550
Figure 9. The contrast of BRMSE between the two DMD mask design strategies.
Figure 9. The contrast of BRMSE between the two DMD mask design strategies.
Sensors 20 00900 g009
Sensors 20 00900 g010 550
Figure 10. The high-resolution MWIR images are reconstructed with (a–d) the proposed DMD masks and (e–h) the original DMD masks. In each row, from left to right, the compression ratios are 0.125, 0.25, 0.375 and 0.5 respectively.
Figure 10. The high-resolution MWIR images are reconstructed with (a–d) the proposed DMD masks and (e–h) the original DMD masks. In each row, from left to right, the compression ratios are 0.125, 0.25, 0.375 and 0.5 respectively.
Sensors 20 00900 g010

Share and Cite

MDPI and ACS Style

Wu, Z.; Wang, X. DMD Mask Construction to Suppress Blocky Structural Artifacts for Medium Wave Infrared Focal Plane Array-Based Compressive Imaging. Sensors 2020, 20, 900. https://doi.org/10.3390/s20030900

AMA Style

Wu Z, Wang X. DMD Mask Construction to Suppress Blocky Structural Artifacts for Medium Wave Infrared Focal Plane Array-Based Compressive Imaging. Sensors. 2020; 20(3):900. https://doi.org/10.3390/s20030900

Chicago/Turabian Style

Wu, Zimu, and Xia Wang. 2020. "DMD Mask Construction to Suppress Blocky Structural Artifacts for Medium Wave Infrared Focal Plane Array-Based Compressive Imaging" Sensors 20, no. 3: 900. https://doi.org/10.3390/s20030900

APA Style

Wu, Z., & Wang, X. (2020). DMD Mask Construction to Suppress Blocky Structural Artifacts for Medium Wave Infrared Focal Plane Array-Based Compressive Imaging. Sensors, 20(3), 900. https://doi.org/10.3390/s20030900

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop