总结《Video rate spectral imaging using a coded aperture snapshot spectral imager》

改进

本模型¹是在single disperser CASSI system（SD-CASSI）² 的基础上提出的改进的单色散元件的CASSI系统，区别主要如下图：

图1（a）SD-CASSI from [2]


图1（b）the direct view CASSI from [1]

1. (a)中的CASSI在色散元件前后均添加了Relay Optics（中继镜），而(b)则只在色散元件后添加一组F/8的中继镜组。
   
   
2. （b）将色散元件更换为一个双阿米西棱镜（a double Amici prism），如图2，以希望获得一个直视配置，则有效地消除了变形失真（a direct view configuration that effectively removes the anamorphic distortion）。所以，从直观上来说，（a）的光路传播明显有一个曲折，而（b）则没有，所以从结构上会更加便捷。
   
   

图2 射线束通过双阿米奇棱镜，并在不同的波长上以不同的方式散射。

接下来将主要以the direct view CASSI为主模型，介绍其构成与数学模型。

光学元件

该成像仪主要包括

1. 物镜(Objective Lens)，物镜将景物成像到编码光圈的平面上。
2. 孔径编码（Coded Aperture），所有实验中使用的代码都是随机的256×248元素二进制图案，如图3所示，最小的代码特征是2×2 CCD像素（每边19.8微米）。这限制了重建数据立方体的空间大小，最大为256×248个空间元素
   
   
   

图3 孔径编码的模板

3. 一个带通滤波器（a bandpass filter）, 将系统的光谱范围限制在450-650纳米
4. 中继镜（A relay lens），中继镜头(Edmund Optics库存部件C45762)将图像从编码光圈的平面传输到CCD上。
5. 色散元件（a dispersive element）用的是一个双阿米西棱镜（a double Amici prism），如图2所示，其中一共包含了三个棱镜，利用棱镜间的多次折射，使得中心波长为550nm的光线通过棱镜是不会发生偏移，与相邻的波长所对应的光线，分别散开再中心轴的两侧，形成一个直视性棱镜，这样可以确保成像仪中的所有元件都可以放置在同一条直线上，系统组件之间的对齐和矫正会变得较之前更容易。同时，考虑到由棱镜到CCD上的不同波段的光线的距离都较为相似，可以极小化基于波长产生的畸变。
6. 一个单色电荷耦合器件（a monochrome charge-coupled device (CCD)），是可以将图像转化为电信号的半导体元件。本系统采用的CCD 为一8位（8-bit）相机，像素为9.9×9.9，在500nm处的响应最强，450-650nm之间的相对响应大于0.7。

数学模型

CASSI成像利用到了压缩感知的原理，这里也需要对其数学原理进行一定的梳理。
　　根据多个博士毕业论文³⁴⁵，the direct view CASSI属于基于空间维的编码形式。（即先通过孔径编码模板后色散，与DD-CASSI 不同。DD由于在通过编码孔径模板之前已经经过一次色散，在采样后光谱维信息会有所变化，所以不同）注意，在数学模型，出于简化模型和更直观的需要，忽视了在光穿过仪器时，光的点扩散方程的在传播过程用的影响。
　　光线通过物镜成像后直接到达编码模板（由于没有预先经过色散，这里编码的仅是二维空间信息，不影响光谱维信息），编码后的光线直接穿过一个450-650nm的带通滤波器，紧接着通过中继镜，到达双阿米西棱镜，，色散棱镜会将经过空间编码的光谱信息沿着光谱维进行色散。光线穿过色散元件，到达CCD传感器上。
　　设物体光谱信息为$f_0(x^{\prime \prime },y^{\prime \prime },\lambda )$ ，其中$x,y$ 表示数据的空间信息，$\lambda$则表示物体的光谱信息。通过孔径编码后的光谱强度分布/光谱信息/功率谱密度（the power spectral density）（这里的几个术语着实有点分不清，）为：
　　$$f_1(x^{\prime },y^{\prime },\lambda )=T(x^{\prime },y^{\prime })f_0(x^{\prime \prime },y^{\prime \prime },\lambda )\text{\hspace{1em}(1.1)}$$
其中$T(x^{\prime },y^{\prime })$为编码孔径模板，即数字化表示的图3。论文⁶中描述，在这里可看成数据立方体 f 0 ( x