本发明属于光学波段显微成像系统及成像方法。
背景技术：
申请人于2005年申请的cn200510018994.0发明专利“光学波段的近场显微镜”，利用循环s矩阵构成的编码模板进行编码，并基于阿达玛变换解码实现光学成像，
所述循环s矩阵所构成的编码模板为制作在透明基板上的金属薄膜上的小孔阵列，阵列形式由循环s矩阵所决定，循环s矩阵可以是3、15、35、63、1023阶循环s矩阵，小孔形状可以是圆形、c形、三角形、栅栏状等异型，孔径可以等于波长、小于波长和大于波长。
基于阿达玛变换实现光学成像技术的基本特征，利用循环s矩阵编码、解码成像时的理论信噪比改善量为
本申请文本中，利用平面xy直角坐标系，xy平面中，x轴为水平轴，向右为正方向，y轴为垂直轴，向上为正方向。
技术实现要素：
本发明提供一种循环s矩阵成像系统及其成像方法，解决现有编码成像方法基于阿达玛变换实现光学成像所存在的信噪比改善量小的问题，以得到更好的图像质量。
本发明所提供的一种循环s矩阵成像系统，包括光源、透明载物片、编码模板、二维微动平台、采样窗口、光学镜头、单元探测器、ad转换器，其特征在于：
所述光源和透明载物片位于编码模板一侧，光学镜头和单元探测器依次位于编码模板另一侧，采样窗口位于光学镜头和编码模板5之间，编码模板能够在二维微动平台驱动下进行移动，单元探测器的敏感面位于光学镜头焦点上，单元探测器接收到的光学信号经过ad转换器转换为数字信号送到计算机10进行处理，以便得到重构图像；
所述编码模板由镀覆在透明基板上的金属薄膜构成，金属薄膜上具有0和1构成的m行n列小孔阵列，其中0为不透光小孔，1为透光小孔；
所述m行n列小孔阵列由r阶循环s矩阵构成，
所述采样窗口为p行q列的矩形窗口，p×q＝r，
当r＝3、7、15、35时，m＝p，n＝(p×q)+(q-1)；例如：3阶循环s矩阵，则采样窗口m＝1、p＝3，n＝3+2、n＝3+2；又如15阶循环s矩阵，可折叠为m×p＝3×5或者m×p＝5×3，若取m＝5、p＝3，则n＝15+2；
当r＝15、35、63、255、511、1023时，所述m行n列小孔阵列中，m＝2p-1，n＝2q-1。
所述r阶循环s矩阵，
r＝3时，首行为101，采样窗口尺寸p×q＝1×3或p×q＝3×1；
r＝7时，首行为1110100，p×q＝1×7或p×q＝7×1；
r＝15时，首行为000100110101111，p×q＝3×5或p×q＝5×3；
r＝35时，首行为00100110101000010011101111100011101，p×q＝7×5或p×q＝5×7：
r＝63时，首行为：
00000100001100010100111101000111001
0010110111011001101010111111，
p×q＝9×7或p×q＝7×9；
r＝225时，首行为：
p×q＝9×25或p×q＝25×9；
以下字符对应关系：以八进制0、1、2、3、4、5、6、7分别表示二进制000、001、010、011、100、101、110、111；
r＝511时，首行为：
p×q＝7×73或p×q＝73×7；
r＝1023时，首行为：
p×q＝33×31或者p×q＝31×33。
本发明所提供的一种循环s矩阵成像方法，包括设置步骤和采样步骤，其特征在于：
(1)设置步骤：
将光源、透明载物片、样品置于编码模板一侧，光学镜头和单元探测器依次置于编码模板另一侧，单元探测器位于光学镜头焦点上，采样窗口位于光学镜头和编码模板之间，待测目标置于光源和编码模板之间，单元探测器通过光学镜头对准编码模板，待测目标、采样窗口和探测器固定；
当编码模板上的小孔为亚波长尺寸时，待测目标与编码模板的距离为近场距离；
(2)采样步骤：
编码模板沿x方向负向移动，移动步长为小孔阵列两相邻元素间距，每移动一次，单元探测器对采样窗口下的编码模板的相应二维区域采样，形成一维行向量光信号，送往ad转换器，所述ad转换器将一维行向量光信号转换为数字信号送到计算机。
计算机可以再借助压缩感知的正交匹配追踪算法进行图像重构。
所述的循环s矩阵成像方法，其进一步特征在于：
所述采样步骤中，采样窗口内二维原始图像对应区域像素按照矩阵方式从第一行起逐行矩阵元素首尾顺序连接构成一维图像行向量
编码模板按步长沿x方向负向每移动一次，编码模板上对应采样窗口内的r阶循环s矩阵从第一行起逐行矩阵元素首尾顺序连接构成一维测量行向量αi，探测器对采样窗口下的编码模板采样，得到测量值
所以对于具有r＝p×q阶循环s矩阵的编码模板，共计沿x方向负向移动r-1次，包含初始位置的测量总共形成r个一维行向量光信号，每个一维行向量光信号经ad转换器转换为数字信号送到计算机，在计算机内可应用正交匹配追踪算法重构出m行p列的目标图像。
当r阶循环s矩阵中，r＝15、35、63、255、511、1023时，所述m行n列小孔阵列中，m＝2p-1，n＝2q-1。可以采用下述成像方法：
本发明所提供的又一种循环s矩阵成像方法，包括设置步骤和采样步骤：
(1)设置步骤：
将光源、、透明载物片、样品置于编码模板一侧，光学镜头和单元探测器依次置于编码模板另一侧，单元探测器位于光学镜头焦点上，采样窗口位于光学镜头和编码模板之间，待测目标置于光源和编码模板之间，单元探测器通过光学镜头对准编码模板，待测目标、采样窗口和单元探测器固定；
当编码模板上的小孔为亚波长尺寸时，待测目标与编码模板的距离为近场距离；
(2)采样步骤：
编码模板沿x方向负向移动，移动步长为小孔阵列两相邻元素间距，每移动一次，单元探测器对采样窗口下的编码模板的相应二维区域采样，形成一维行向量光信号，送往ad转换器，所述ad转换器将一维行向量光信号转换为数字信号送到计算机；当编码模板沿x方向负向移动到采样窗口右边缘时，编码模板复位；
编码模板沿y方向正向移动，移动步长为小孔阵列两相邻元素间距，然后再沿x方向负向移动，移动步长为小孔阵列两相邻元素间距，每移动一次，单元探测器对采样窗口下的编码模板的相应二维区域采样，形成一维行向量光信号，送往ad转换器，所述ad转换器将一维行向量光信号转换为数字信号送到计算机；直至编码模板沿y方向正向移动到采样窗口下边缘为止；
计算机可以再借助压缩感知的正交匹配追踪算法进行图像重构。
所述采样步骤中，采样窗口内二维原始图像对应区域像素按照矩阵方式从第一行起逐行矩阵元素首尾顺序连接构成一维图像行向量
编码模板按步长沿x方向负向以及沿y方向正向每移动一次，编码模板上对应采样窗口内的r阶循环s矩阵从第一行起逐行矩阵元素首尾顺序连接构成一维测量行向量αi，探测器对采样窗口下的编码模板采样，得到测量值
对于具有r＝p×q阶循环s矩阵的编码模板，共计沿x方向负向和沿y方向正向移动r-1次，包含初始位置的测量总共形成r个一维行向量光信号，每个一维行向量光信号经ad转换器转换为数字信号送到计算机，在计算机内可应用正交匹配追踪算法重构出m行p列的目标图像。
所述ad转换器送来的序列数字信号在计算机内可以借助通用的压缩感知方法中的正交匹配追踪算法(omp)进行目标图像的重构，这种方法的核心思想是通过多次迭代、逐步选取的方式从采样窗口(测量矩阵)中获取用来表示列原子的行向量。每次迭代中都要保证选择的原子与当前残余(冗余)向量有最大的相关性，都能满足最佳原子匹配的准则，再完成从观测向量中减去相关部分的步骤来更新残差向量，直到迭代次数达到稀疏度k，然后停止迭代。
所述单元探测器为与待探测的波段相适应的单元探测器。
采样时，单元探测器和目标固定，单元探测器对应于如图3所示的采样窗口内编码模板中的一个适当折叠(例如接近正方形的折叠)的二维区域采样，下一次采样前编码模板需要相应移动(例如图3中所示5×3的小区)。
本发明通过矩阵标准乘法的方式来获取测量值y＝a×x，a为测量矩阵，x为图像(目标)矩阵。所以，为方便计算，把二维原始图像x展开成一维向量信号
本发明利用单个单元探测器，借助编码模板、通过采样窗口对目标进行采样，所得到的r个一维行向量光信号便于再借用压缩感知方法进行处理重构图像，获得了比阿达玛变换光学成像更高的信噪比改善量，最终得到的图像质量更好，适用于从紫外波段到太赫兹波段的整个光学波段。
附图说明
图1为本发明的成像系统结构示意图；
图2为采样窗口采样方式示意图；
图3为采样所形成的15阶循环s矩阵示意图；
图4(a)左、右依次为3×1部分原始图像、3阶重构图像；
图4(b)左、右依次为5×3部分原始图像、15阶重构图像；
图4(c)左、右依次为7×5部分原始图像、35阶重构图像；
图4(d)左、右依次为9×7部分原始图像、63阶重构图像；
图4(e)左、右依次为33×31部分原始图像、1023阶重构图像。
具体实施方式
以下结合实施例和附图对本发明进一步说明。
如图1所示，本发明的循环s矩阵成像系统，包括光源1、编码模板5、二维微动平台4、采样窗口6、光学镜头7、单元探测器8、ad转换器9；
所述光源1位于编码模板5一侧，光学镜头7和单元探测器8依次位于编码模板另一侧，采样窗口6位于光学镜头7和编码模板5之间，编码模板能够在二维微动平台4驱动下进行移动，单元探测器的敏感面位于光学镜头焦点上，单元探测器接收到的光学信号经过ad转换器9转换为数字信号送到计算机10进行处理，以便得到重构图像；
所述编码模板由镀覆在透明基板上的金属薄膜构成，金属薄膜上具有0和1构成的m行n列小孔阵列，其中0为不透光小孔，1为透光小孔；
所述m行n列小孔阵列由r阶循环s矩阵构成，r＝3、15、35、63、255、511、1023，可折叠为二维阵列；
所述采样窗口为m行p列的矩形窗口，m×p＝r，n＝(m×p)+(p-1)。
采样时，待测目标3放置在承载板2上，置于光源1和编码模板5之间，单元探测器8通过光学镜头7对准编码模板5，待测目标3、采样窗口6和单元探测器8固定。
如图2、图3所示，以15阶循环s矩阵为例，测量矩阵a每一行的行向量①、②、③…分别对应图3的采样窗口内编码模板上对应二维区域①、②、③…，
如行向量①、实际就是编码模板上对应采样窗口内的第一个方框逐行顺序连接起来的一行，采集开始时，如图3所示，单元探测器对应采样窗口内的编码模板二维区域①，即一个5行3列的矩阵采集一次图像，然后目标图像和单元探测器固定不动，编码模板5沿x方向负向移动一位，单元探测器对应采样窗口内的编码模板二维区域②，即又一个5行3列的矩阵再采集第二次图像，按照这种模式，即可完成从左往右的循环移动编码采集过程，得到15个值，作为重构算法所需测量值y。
实施例：3阶、15阶、35阶、63阶和1023阶可折叠循环s矩阵，分别按照1×3；3×5；5×7；7×9；31×33折叠作成采样窗口，对女性头像图采样并且经正交匹配追踪法(omp)方法重构的相应区域的图像并给出其重构特性，
表1.几个不同阶数的s矩阵在正交匹配追踪法(omp)下的重构特性
表1中，mse为均方误差，re为相对误差，psnr为峰值信噪比，mdb为信噪比改善量。
在上述折叠方式下，五种阶数的循环s矩阵的重构效果分别比对应的经典阿达玛变换成像方法的理论信噪比改善量提高了多倍。图4(a)～图4(e)各图的右侧图则为相应的重构图像，几乎和原图相差不大。
表1中，分别取女性头像图中间3×1、5×3、7×5、9×7、33×31的部分区域，然后利用相应阶数的循环s矩阵，为便于和经典的阿达玛变换成像方法的理论结果比较，这里给出在100％采样率下，通过omp重构算法计算得到的结果。
各阶循环s矩阵设计的单探测器条件下超分辨成像模板都有很好的成像效果。重构图像的效果比较传统阿达玛变换光学成像的理论信噪比改善量大幅提高。这些矩阵也为实现单探测器条件下的压缩感知超衍射分辨宽视场光学成像提供了重要的参考。
如果要得到大面积的目标图像，可以采用分区采样重构再合成的办法实现。