<?xml version="1.0" encoding="UTF-8"?><!DOCTYPE article  PUBLIC "-//NLM//DTD Journal Publishing DTD v3.0 20080202//EN" "http://dtd.nlm.nih.gov/publishing/3.0/journalpublishing3.dtd"><article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" dtd-version="3.0" xml:lang="en" article-type="research article"><front><journal-meta><journal-id journal-id-type="publisher-id">JAMP</journal-id><journal-title-group><journal-title>Journal of Applied Mathematics and Physics</journal-title></journal-title-group><issn pub-type="epub">2327-4352</issn><publisher><publisher-name>Scientific Research Publishing</publisher-name></publisher></journal-meta><article-meta><article-id pub-id-type="doi">10.4236/jamp.2015.32039</article-id><article-id pub-id-type="publisher-id">JAMP-53709</article-id><article-categories><subj-group subj-group-type="heading"><subject>Articles</subject></subj-group><subj-group subj-group-type="Discipline-v2"><subject>Physics&amp;Mathematics</subject></subj-group></article-categories><title-group><article-title>
 
 
  Spectral Imagers with Linear Detector Imager Systems Based on Spectrum Compressed
 
</article-title></title-group><contrib-group><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Xiaoming</surname><given-names>Zhong</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref><xref ref-type="corresp" rid="cor1"><sup>*</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Huang</surname><given-names>Li</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Huang</surname><given-names>Li</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib></contrib-group><aff id="aff1"><addr-line>Beijing Institute of Space Mechanics &amp;amp; Electricity, Beijing, China</addr-line></aff><author-notes><corresp id="cor1">* E-mail:<email>xiaoming.zhong@163.com(XZ)</email>;</corresp></author-notes><pub-date pub-type="epub"><day>30</day><month>01</month><year>2015</year></pub-date><volume>03</volume><issue>02</issue><fpage>267</fpage><lpage>271</lpage><history><date date-type="received"><day>14</day>	<month>December</month>	<year>2014</year></date></history><permissions><copyright-statement>&#169; Copyright  2014 by authors and Scientific Research Publishing Inc. </copyright-statement><copyright-year>2014</copyright-year><license><license-p>This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/</license-p></license></permissions><abstract><p>
 
 
   Traditional spectral imagers require 2-dimensional detectors. We present a new method to implement spectral imagers with linear detector imager systems based on spectrum compressed. Using 1-dimension detectors instead of 2-dimension detectors to get 3-dimensional data cubes, the spectral imagers could get both the spectral information and the spatial information of each ground object. By the method of characteristics decoupling, we make high precision reconstruction of compressed data. Theoretical analysis and simulations show that it not only ensures the imaging quality but also reduces the dimension of the detectors and complexity of imaging system greatly. 
 
</p></abstract><kwd-group><kwd>Linear Detector Imager</kwd><kwd> Spectrum Compressed</kwd><kwd> Characteristics Decoupling</kwd></kwd-group></article-meta></front><body><sec id="s1"><title>1. Introduction</title><p>In most push-broom imaging spectrometers, they employ split structures. This will cause the spatial resolution decreasing. We present a new approach to implement computational spectral imaging, bring forward research on the mechanism of spectral imagers with linear detector imager systems on spectrum compressed, present a detailed theoretical analysis and give a simulation on this algorithm.</p></sec><sec id="s2"><title>2. The Traditional Push-Broom Imaging Spectrometer’s Principle</title><p>Push-broom imaging principle of scanning imaging spectrometer is shown as <xref ref-type="fig" rid="fig1">Figure 1</xref> shows: through a slit in the direction parallel to the slit dimension and imaging in the vertical slit direction to obtain a two-dimensional target imaging [<xref ref-type="bibr" rid="scirp.53709-ref1">1</xref>].</p></sec><sec id="s3"><title>3. Spectral Imaging with Linear Detector Imager System Model</title><sec id="s3_1"><title>3.1. System Introduction</title><p>Spectral imager with linear detector imager systems are based on the traditional dispersive spectral imaging</p><fig id="fig1"  position="float"><label><xref ref-type="fig" rid="fig1">Figure 1</xref></label><caption><title>Traditional push-broom imaging spectrometer</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/53709x3.png"/></fig><p>technology. The code template has a special mathematical code form. The spatial information and spectral of the original target can be modulated with the code template, using the theory of compressive sensing. We can reconstruct one-dimensional special information and one-dimensional spectral information from the map of aliasing imaging. We can obtain the two-dimensional special information with the push-broom imaging model.</p><p>As <xref ref-type="fig" rid="fig2">Figure 2</xref> shows, the system consists of (1) an objective lens, (2) a coded aperture, (3) relay lens, (4) dispersion prism, and (5) a monochrome linear charge-coupled device (CCD) detector. The objective lens images the scene on to the plane of the coded aperture.</p></sec><sec id="s3_2"><title>3.2. Code Process</title><p>The spectral imager with linear detector imager system realizes a single shot compressive spectral imaging system. As <xref ref-type="fig" rid="fig3">Figure 3</xref> shows, it encoded both 1D spatial and 1D spectral information of objects through an aperture code projection that is captured after it propagates through a dispersive element, a linear detector then collects all light passing through the aperture and the dispersive element. It is important to emphasize that the code</p><p>aperture pattern remains fix in the sampling process. Suppose that the scene or object is represented by <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/53709x4.png" xlink:type="simple"/></inline-formula> [<xref ref-type="bibr" rid="scirp.53709-ref2">2</xref>], where <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/53709x5.png" xlink:type="simple"/></inline-formula> is the wavelength and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/53709x6.png" xlink:type="simple"/></inline-formula> correspond to the spatial position, in discrete form it is denoted as<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/53709x7.png" xlink:type="simple"/></inline-formula>. Suppose that the code pattern is <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/53709x8.png" xlink:type="simple"/></inline-formula> [<xref ref-type="bibr" rid="scirp.53709-ref3">3</xref>].<sup> </sup>The term <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/53709x9.png" xlink:type="simple"/></inline-formula> takes into account all possible noise source. Then the</p><p>signal in front of the array detector can be expressed by (1) [<xref ref-type="bibr" rid="scirp.53709-ref4">4</xref>].</p><p>Example</p><disp-formula id="scirp.53709-formula548"><label>(1)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/53709x10.png"  xlink:type="simple"/></disp-formula></sec><sec id="s3_3"><title>3.3. The Reconstruction Algorithm</title><p>TwIST (Two-step Iterative Shrinkage Thresholding) is an algorithm framework introduced by Bioucas-Dias and Figueiredo [<xref ref-type="bibr" rid="scirp.53709-ref5">5</xref>]. In the context of image reconstruction for a CASSI system, in this paper, we reconstruct the image</p><fig id="fig2"  position="float"><label><xref ref-type="fig" rid="fig2">Figure 2</xref></label><caption><title> Spectral imagers system composition</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/53709x11.png"/></fig><fig id="fig3"  position="float"><label><xref ref-type="fig" rid="fig3">Figure 3</xref></label><caption><title> Spectral imagers on spectrum compressed</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/53709x12.png"/></fig><p>of spectral imagers with linear detector imager system with TwIST. TwIST describes a solution as the solution to the generalized non-linear, unconstrained minimization problem:</p><disp-formula id="scirp.53709-formula549"><label>(2)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/53709x13.png"  xlink:type="simple"/></disp-formula><p>Here the linear operator <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/53709x14.png" xlink:type="simple"/></inline-formula> represents the system forward model. Reconstruction of <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/53709x15.png" xlink:type="simple"/></inline-formula> is attained by solving the optimization problem. The first term minimizes the difference between the model and the measurement<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/53709x16.png" xlink:type="simple"/></inline-formula>. The variable tau &gt; 0 controls the level of sparsity attained in the reconstruction. The sparser the source<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/53709x17.png" xlink:type="simple"/></inline-formula>, the better the performance of the reconstruction algorithm. In this work, the CS reconstruction algorithm was used to solve the Equation (2). The above procedure tries to recover the overall data cube with only one measurement and hence it often yields a low SNR output performance.</p><p>Before recovering spectral data, spectra data should be reprocessed to be normalized. Normalized data can be directly used algorithm for recovering. Restoration of the accuracy of the results depend on the end of the iteration condition.</p><disp-formula id="scirp.53709-formula550"><label>(3)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/53709x18.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.53709-formula551"><label>(4)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/53709x19.png"  xlink:type="simple"/></disp-formula><p>Before restoration, we need to set two parameters, the balance to the noise function and the end of iterative threshold tol.</p><p>The termination function in defined as</p><disp-formula id="scirp.53709-formula552"><label>(5)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/53709x20.png"  xlink:type="simple"/></disp-formula><p>The termination function is tol. We stop iterativing if<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/53709x21.png" xlink:type="simple"/></inline-formula>.</p><p>The process for the algorithm is as following:</p><p>(1) obtain the parameters<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/53709x22.png" xlink:type="simple"/></inline-formula>, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/53709x22.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/53709x23.png" xlink:type="simple"/></inline-formula>,<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/53709x22.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/53709x23.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/53709x24.png" xlink:type="simple"/></inline-formula>.</p><p>(2) imaging initialization, enter the initial value<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/53709x25.png" xlink:type="simple"/></inline-formula>, setting calculation tol, the corresponding objective function<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/53709x25.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/53709x26.png" xlink:type="simple"/></inline-formula>;</p><p>(3) deal with <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/53709x27.png" xlink:type="simple"/></inline-formula> for removing the noise, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/53709x27.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/53709x28.png" xlink:type="simple"/></inline-formula>, use the noise function to get a new<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/53709x27.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/53709x28.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/53709x29.png" xlink:type="simple"/></inline-formula>, numerical noise function is selected , minimizing the total variational method to obtain the corresponding objective function, calculate<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/53709x27.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/53709x28.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/53709x29.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/53709x30.png" xlink:type="simple"/></inline-formula>;</p><p>(4) compare <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/53709x31.png" xlink:type="simple"/></inline-formula> with<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/53709x31.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/53709x32.png" xlink:type="simple"/></inline-formula>, if that the calculated result is correct, this time need to be adjusted to balance parameters in noise function;</p><p>(5) calculate <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/53709x33.png" xlink:type="simple"/></inline-formula> with <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/53709x33.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/53709x34.png" xlink:type="simple"/></inline-formula> and<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/53709x33.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/53709x34.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/53709x35.png" xlink:type="simple"/></inline-formula>, repeat step 1 - 3, then update <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/53709x33.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/53709x34.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/53709x35.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/53709x36.png" xlink:type="simple"/></inline-formula> with<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/53709x33.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/53709x34.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/53709x35.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/53709x36.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/53709x37.png" xlink:type="simple"/></inline-formula>, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/53709x33.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/53709x34.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/53709x35.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/53709x36.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/53709x37.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/53709x38.png" xlink:type="simple"/></inline-formula>with<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/53709x33.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/53709x34.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/53709x35.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/53709x36.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/53709x37.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/53709x38.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/53709x39.png" xlink:type="simple"/></inline-formula>, the iteration continues;<sub> </sub></p><p>(6) the iteration does not stop until the function is less than tol, if<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/53709x40.png" xlink:type="simple"/></inline-formula>, the iteration stops, then we can obtain the data cube reconstruction.</p></sec></sec><sec id="s4"><title>4. Simulation</title><p>The test uses MATLAB to carry on the simulation and reconstruction. The process is as following.</p><p>Firstly, from the realistic point of view, we choose online downloading remote sensing image as a simulation of ground targets. As <xref ref-type="fig" rid="fig4">Figure 4</xref> shows, this multispectral database has 33 bands from 400 nm to 720 nm wavelength with a 10 nm step. The image size is 480 &#215; 752 pixels. In the push-broom imaging simulation process, simulation design parameters as follows: the spectral range is 400 - 720 nm; spectrum channel number is 33, detector pixels is 280 &#215; 1, Intercepting 256 &#215; 280 pixels format as the push-broom scene spectrometer remote sensing image above.</p><p>Secondly, simulate smashing on the CCD. Different from CASSI, ours system coded model is 1-dimen- sional. Be sure that spectral dispersion direction and detectors arrangement direction is consistent. Linear array detector will get aliasing information, and with the push-broom imaging process, the entire scene information will be obtained. <xref ref-type="fig" rid="fig5">Figure 5</xref> and <xref ref-type="fig" rid="fig6">Figure 6</xref> show the coding model and the smashing result on the CCD.</p><p>Finally, with the prior information of the arithmetic coding model and algorithms, we can recover the spectral data cube from the aliasing image data. In our simulation, we choose TwIST for recovering for its accuracy and speed are all well. <xref ref-type="fig" rid="fig7">Figure 7</xref> shows the result of image reconstruction.</p><fig id="fig4"  position="float"><label><xref ref-type="fig" rid="fig4">Figure 4</xref></label><caption><title> Original scene</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/53709x41.png"/></fig><fig id="fig5"  position="float"><label><xref ref-type="fig" rid="fig5">Figure 5</xref></label><caption><title> Coding model</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/53709x42.png"/></fig><fig id="fig6"  position="float"><label><xref ref-type="fig" rid="fig6">Figure 6</xref></label><caption><title> Coding imaging</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/53709x43.png"/></fig><fig id="fig7"  position="float"><label><xref ref-type="fig" rid="fig7">Figure 7</xref></label><caption><title> Reconstructing result</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/53709x44.png"/></fig></sec><sec id="s5"><title>5. Result</title><p>From above, we can see that spectral imagers with linear detector imager system on spectrum compressed can make use of a 1-dimension detector instead of a 2-dimension detector to get a 3-dimensional data cube, which includes both the spectral information and the spatial information of each ground object.</p><p>It can decrease the dimensions of the detectors and complexity of imaging system greatly. This work has important and far-reaching theoretical significance and extensive applied values to improve the technique levels of digital remote sensing, and it can satisfy exigent requirements for remote sensing spectral images internally.</p></sec><sec id="s6"><title>Acknowledgements</title><p>This work was supported by the National Natural Science Foundation of China (No. 11204014).</p></sec><sec id="s7"><title>Cite this paper</title><p>Xiaoming Zhong,Huang Li,Huang Li, (2015) Spectral Imagers with Linear Detector Imager Systems Based on Spectrum Compressed. Journal of Applied Mathematics and Physics,03,267-271. doi: 10.4236/jamp.2015.32039</p></sec></body><back><ref-list><title>References</title><ref id="scirp.53709-ref1"><label>1</label><mixed-citation publication-type="other" xlink:type="simple">Bioucas-Dias, J.M. and Figueiredo, M.A. (2007) A New TWIST: Two-Step Iterative Shrinkage, Thresholding Algorithms for Image Restoration. IEEE Trans Image Process, 16, 2992-3004. 
http://dx.doi.org/10.1109/TIP.2007.909319</mixed-citation></ref><ref id="scirp.53709-ref2"><label>2</label><mixed-citation publication-type="other" xlink:type="simple">Wagadarikar, A.A., Pitsianis, N.P. and Sun, X.B. (2008) Spectral Image Estimation for Coded Aperture Snapshot Spectral Imagers. Proc SPIE, 7076, 7076-7076-02.</mixed-citation></ref><ref id="scirp.53709-ref3"><label>3</label><mixed-citation publication-type="other" xlink:type="simple">Arguello, H., Ye, P. and Arce, G.R. (2010) Spectral Aperture Code Design for Multi-Shot Compressive Spectral Imaging. OSA Optics &amp; Photonics Congress, 978-1-55752-887-2.</mixed-citation></ref><ref id="scirp.53709-ref4"><label>4</label><mixed-citation publication-type="other" xlink:type="simple">Arguello, H. and Arce, G.R. (2011) Code Aperture Agile Spectral Imaging (CAASI). OSA Optics and Photonics Congress, 978-1-55752-914-5.</mixed-citation></ref><ref id="scirp.53709-ref5"><label>5</label><mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Shen</surname><given-names> Z. </given-names></name>,<etal>et al</etal>. (<year>2002</year>)<article-title>The Principle of the Spaceborne Hyperspectrum Imager</article-title><source> Spacecraft Recovery&amp; Remote Sensing</source><volume> 23</volume>,<fpage> 28</fpage>-<lpage>34</lpage>.<pub-id pub-id-type="doi"></pub-id></mixed-citation></ref></ref-list></back></article>