自适应Radon单像素成像

自适应Radon单像素成像 Adaptive radon single-pixel imaging method 1 2 first-author 王蔚松（1995-），山东威海人，硕士，2018年于中国海洋上大学获得学士学位，现为中科院长春光学精密机械与物理研究所硕士研究生，主要从事单像素成像方面的研究。E-mail： wangweisong258@126.com http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19367935&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19367943&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19367939&type= 22.01332855 32.00400925 王蔚松 （1995-），山东威海人，硕士，2018年于中国海洋上大学获得学士学位，现为中科院长春光学精密机械与物理研究所硕士研究生，主要从事单像素成像方面的研究。E-mail： wangweisong258@126.com wangweisong258@126.com 1 1 1 1 1 corresp E-mail： optlab@ciomp.ac.cn E-mail： optlab@ciomp.ac.cn optlab@ciomp.ac.cn 中国科学院 长春光学精密机械与物理研究所 光学系统先进制造技术重点 实验室，吉林 长春 130033 中国科学院 长春光学精密机械与物理研究所 光学系统先进制造技术重点 实验室，吉林 长春 130033 Key Laboratory of Optical System Advanced Manufacturing Technology， Changchun Institute of Optics， Fine Mechanics and Physics， Chinese Academy of Sciences， Changchun 130033， China Key Laboratory of Optical System Advanced Manufacturing Technology， Changchun Institute of Optics， Fine Mechanics and Physics， Chinese Academy of Sciences， Changchun 130033， China 中国科学院大学，北京 100049 中国科学院大学，北京 100049 University of Chinese Academy of Sciences， Beijing 100049， China University of Chinese Academy of Sciences， Beijing 100049， China 使用单像素探测器实现成像需要大量采样。对于目标区域仅占场景一部分的情况时，我们提出了自适应Radon单像素成像方法，能够使用单像素探测器实现目标区域的定位和成像。本文对该方法的目标定位方式、编码采样算法、重建算法等进行研究，以减少单像素成像的采样数量。基于Radon变换的基本原理，使用图像在水平和垂直方向的投影信息，以获取场景中目标区域的大小和位置。建立自适应Radon-Hadamard单像素成像模型，仅对目标区域进行单像素采样，然后使用滤波反投影技术重建目标区域。研究结果表明：所提出的自适应Radon单像素成像方法能够实现对场景中目标区域的成像，采样数量远低于重建图像的分辨率，重建图像的结构相似性系数大于95%，有效的提高了单像素成像方法的成像效率。 Single-pixel imaging requires a large amount of sampling. In this study， an adaptive Radon single-pixel imaging method is proposed for the target region that only occupies a part of the scene. This method uses single-pixel detectors to position and image the target region. We used the target positioning method， coding sampling， and reconstruction algorithms to reduce the number of single-pixel imaging samples. Based on the fundamental principle of Radon transformation， the projection information of an image in horizontal and vertical directions was used to obtain the size and position of the target region in the scene. This method established the adaptive Radon-Hadamard single-pixel imaging model. Only single-pixel sampling was performed on the target region and filtered back-projection technology was used to reconstruct the target region. The results show that the proposed adaptive Radon single-pixel imaging method can achieve imaging of the target region in a scene. The number of samples was much lower than the resolution of the reconstructed image and the structural similarity index of the reconstructed image was greater than 95%， which effectively improved the imaging efficiency of single-pixel imaging method. 单像素成像 Radon变换 目标区域 采样数量 single-pixel imaging radon transform target region sampling number 1 引　言 1 单像素成像（Single Pixel Imaging， SPI）作为新兴的计算成像技术，其利用结构光调制，仅使用一个不具备空间分辨能力的单像素探测器实现成像，能够感知微弱光强的变化，具有很高的探测效率 ［ 1 1 ］ ，随着硬件性能的提升，且价格相较于面阵探测器更便宜，SPI在不可见光波段如红外、X、荧光、太赫兹波段具有优势。SPI能够在时间维度上编码空间信息 ［ 2 2 ］ ，记录一维的光强信号，因此具有极高的时间分辨率，可应用于时间域的成像 ［ 3 3 ］ 或发展超快成像领域 ［ 4 4 - 5 5 4-5 ］ 。 单像素成像是通过牺牲时间换取空间的成像方法，对场景需要进行大量的采样。采样的结构光基底分为非正交基底和正交基底，通过参考光（结构光）与信号光之间的关联运算进行成像 ［ 6 6 – 8 8 6–8 ］ 。使用随机散斑或者二进制图案作为非正交基底对光场进行编码，需要远大于重建图像分辨率的采样次数才能够恢复清晰的图像 ［ 9 9 - 10 10 9-10 ］ ，比如Sun使用了106种随机图案进行了百万次测量来重建分辨率为256×192的图像 ［ 11 11 ］ 。即使使用空间光调制器件加快了调制速率，采样数量和成像时间仍是限制单像素成像技术应用的重要因素。压缩感知技术指出，低于奈奎斯特采样频率的信号也能够恢复原始信号 ［ 12 12 ］ ，基于压缩感知算法的单像素成像算法极大的降低采样次数，并且能够恢复低噪声和高质量的图像 ［ 13 13 - 16 16 13-16 ］ 。但随之带来较长的重建时间依然限制了成像效率。 此时一些具有正交性的感知矩阵进入研究人员的视野。根据变换基底的不同，产生了基于小波变换 ［ 17 17 - 18 18 17-18 ］ 、离散余弦变换 ［ 19 19 ］ 、傅里叶变换 ［ 20 20 - 22 22 20-22 ］ 、哈达玛变换 ［ 23 23 - 25 25 23-25 ］ 等单像素成像技术，在数学上促进了各成像领域之间的融合。单采样的次数等价于图像分辨率大小，因此仍表现为低帧速率成像。为降低采样次数，Bian提出的eSPI方法，提取并恢复频域占比最大的区域，将采样数量降低了两个数量级 ［ 26 26 ］ ，Shi等人在稀疏傅里叶单像素成像的基础上结合了压缩感知技术，在降低采样数量的同时提高了成像的质量 ［ 27 27 ］ 。实际上，场景中的有效信息仅占场景的一部分，Sun等人使用小波变换提取目标中的感兴趣区域，并使用哈达玛单像素成像方法对感兴趣区域进行高分辨率重建 ［ 17 17 ］ ，其属于多分辨率成像的范畴；Jiang等人通过傅里叶切片定理定位一个目标区域，并在目标区域自适应的进行傅里叶单像素成像，降低了采样图案的数量 ［ 28 28 ］ ，但只能够对单一目标成像；Shi给出了使用Hadamard变换实现目标边缘定位的方案，具有极高的检测速度 ［ 29 29 ］ ，但需要背景信息作为先验，并且没有实现目标区域的成像。 Radon变换建立了物体内部元素与投影数据之间的联系，Harald最早将Radon变换用于单像素成像 ［ 30 30 ］ ，但使用旋转的机械结构和标线圆盘提取Radon变换系数会产生较大的误差，且标线圆盘灵活性差；Shi通过空间光调制器产生的结构光对整个场景进行调制，实现了Radon单像素成像（Radon Single-Pixel Imaging， RSI） ［ 31 31 ］ ，需要等同于图像分辨率的采样次数才能够实现高质量成像。本文在Radon单像素成像的基础上，提出了自适应Radon单像素成像方法（Adaptive Radon Single-Pixel Imaging，ARSI），使用Radon投影变换定位目标区域的基础上，匹配的基底图案大小，对目标区域进行单像素采样和重建，提高了单像素成像的成像效率。 2 基本原理 1 2.1　基于ARSI的目标区域定位 2 如 图1 图1 （a）所示，场景的分辨率为256×256，场景中的有效信息仅集中在矩形框中，我们称之为目标区域。使用传统的单像素成像方法，在不结合压缩感知技术的前提下，至少需要256×256次采样，才能够清楚地重建整幅场景图片，以保证目标区域成像质量。实际上，背景中基本包含很少量的信息，在采样过程中可以被忽略，因此仅对目标区域进行采样与重建能够大量减少单像素成像的采样数量。ARSI方法的关键问题分为使用目标区域的定位与重建。 Radon变换可被描述为：给定图像Radon变换是图像沿给定角度的线积分 ［ 32 32 ］ ，也被形容为图像在给定方向上的投影。Radon变换的数学表达式为： <math display="block" id="M1"><mi>L</mi><mo>:</mo><mi>s</mi><mo>=</mo><mi>x</mi><mi mathvariant="normal">c</mi><mi mathvariant="normal">o</mi><mi mathvariant="normal">s</mi><mi>θ</mi><mo>+</mo><mi>y</mi><mi mathvariant="normal">s</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">n</mi><mi>θ</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19367969&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19367962&type= 32.85066605 3.72533321 ， <math display="block" id="M2"><mtable columnalign="left"><mtr><mtd><mi>r</mi><mo stretchy="false">(</mo><mi>s</mi><mo>,</mo><mi>θ</mi><mo stretchy="false">)</mo><mo>=</mo><msub><mo>∫</mo><mrow><mi>L</mi></mrow></msub><mrow><mi>f</mi><mo stretchy="false">(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo stretchy="false">)</mo><mi mathvariant="normal">d</mi><mi>L</mi></mrow><mo>=</mo></mtd></mtr><mtr><mtd><msubsup><mo>∫</mo><mrow><mo>-</mo><mi mathvariant="normal">∞</mi></mrow><mrow><mi mathvariant="normal">∞</mi></mrow></msubsup><mrow><msubsup><mo>∫</mo><mrow><mo>-</mo><mi mathvariant="normal">∞</mi></mrow><mrow><mi mathvariant="normal">∞</mi></mrow></msubsup><mrow><mi>f</mi><mo stretchy="false">(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo stretchy="false">)</mo><mi>δ</mi><mo stretchy="false">(</mo><mi>x</mi><mi mathvariant="normal">c</mi><mi mathvariant="normal">o</mi><mi mathvariant="normal">s</mi><mi>θ</mi><mo>+</mo><mi>y</mi><mi mathvariant="normal">s</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">n</mi><mi>θ</mi><mo>-</mo><mi>s</mi><mo stretchy="false">)</mo><mi mathvariant="normal">d</mi><mi>x</mi><mi mathvariant="normal">d</mi><mi>y</mi></mrow></mrow></mtd></mtr></mtable></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19367978&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19367974&type= 67.47933960 15.49400043 ， L 是 xOy 面上的直线， s 是原点到 L 的距离， <math id="M3"><mi>θ</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19368010&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19368004&type= 1.77800000 2.79399991 是 L 的法线与 x 轴正向的交角。 s 和 <math id="M4"><mi>θ</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19368010&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19368004&type= 1.77800000 2.79399991 是直线 L 的位置参数，给定 s 和 <math id="M5"><mi>θ</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19368010&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19368004&type= 1.77800000 2.79399991 时 L 被唯一确定。 <math id="M6"><mi>r</mi><mo stretchy="false">(</mo><mi>s</mi><mo>,</mo><mi>θ</mi><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19368414&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19368411&type= 10.32933331 3.55599999 的集合被称为 <math id="M7"><mi>f</mi><mo stretchy="false">(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19368432&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19368430&type= 11.43000031 3.80999994 的Radon变换，其中项 δ 是狄拉克分布函数。 根据Radon投影变换定理，图像在不同方向的投影积分曲线表现出不同的分布特征，这是ARSI定位目标区域的依据，因此我们分别提取图像在 x 方向和 y 方向的投影曲线，此时 <math id="M8"><mi>θ</mi><mo>=</mo><mn mathvariant="normal">0</mn></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19368285&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19368281&type= 8.55133343 2.79399991 和 <math id="M9"><mi>θ</mi><mo>=</mo><mi mathvariant="normal">π</mi><mo>/</mo><mn mathvariant="normal">2</mn></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19368291&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19368288&type= 12.27666664 3.30200005 ， <math display="block" id="M10"><mi>r</mi><msub><mrow><mo stretchy="false">(</mo><mi>s</mi><mo stretchy="false">)</mo></mrow><mrow><mi>θ</mi><mo>=</mo><mn mathvariant="normal">0</mn></mrow></msub><mo>=</mo><msubsup><mo>∫</mo><mrow><mo>-</mo><mi mathvariant="normal">∞</mi></mrow><mrow><mi mathvariant="normal">∞</mi></mrow></msubsup><mrow><msubsup><mo>∫</mo><mrow><mo>-</mo><mi mathvariant="normal">∞</mi></mrow><mrow><mi mathvariant="normal">∞</mi></mrow></msubsup><mrow><mi>f</mi><mo stretchy="false">(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo stretchy="false">)</mo><mi>δ</mi><mo stretchy="false">(</mo><mi>x</mi><mo>-</mo><mi>s</mi><mo stretchy="false">)</mo><mi mathvariant="normal">d</mi><mi>x</mi><mi mathvariant="normal">d</mi><mi>y</mi></mrow></mrow><mo>=</mo><msubsup><mo>∫</mo><mrow><mo>-</mo><mi mathvariant="normal">∞</mi></mrow><mrow><mi mathvariant="normal">∞</mi></mrow></msubsup><mrow><mi>f</mi><mo stretchy="false">(</mo><mi>s</mi><mo>,</mo><mi>y</mi><mo stretchy="false">)</mo><mi mathvariant="normal">d</mi><mi>y</mi></mrow></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19368059&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19368054&type= 91.44000244 6.94266701 ， <math display="block" id="M11"><mtable columnalign="left"><mtr><mtd><mi>r</mi><msub><mrow><mo stretchy="false">(</mo><mi>s</mi><mo stretchy="false">)</mo></mrow><mrow><mi>θ</mi><mo>=</mo><mi mathvariant="normal">π</mi><mo>/</mo><mn mathvariant="normal">2</mn></mrow></msub><mo>=</mo><msubsup><mo>∫</mo><mrow><mo>-</mo><mi mathvariant="normal">∞</mi></mrow><mrow><mi mathvariant="normal">∞</mi></mrow></msubsup><mrow><msubsup><mo>∫</mo><mrow><mo>-</mo><mi mathvariant="normal">∞</mi></mrow><mrow><mi mathvariant="normal">∞</mi></mrow></msubsup><mrow><mi>f</mi><mo stretchy="false">(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo stretchy="false">)</mo><mi>δ</mi><mo stretchy="false">(</mo><mi>y</mi><mo>-</mo><mi>s</mi><mo stretchy="false">)</mo><mi mathvariant="normal">d</mi><mi>x</mi><mi mathvariant="normal">d</mi><mi>y</mi></mrow></mrow><mo>=</mo></mtd></mtr><mtr><mtd><msubsup><mo>∫</mo><mrow><mo>-</mo><mi mathvariant="normal">∞</mi></mrow><mrow><mi mathvariant="normal">∞</mi></mrow></msubsup><mrow><mi>f</mi><mo stretchy="false">(</mo><mi>x</mi><mo>,</mo><mi>s</mi><mo stretchy="false">)</mo><mi mathvariant="normal">d</mi><mi>x</mi></mrow><mo>,</mo></mtd></mtr></mtable></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19368067&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19368064&type= 70.95066071 15.66333294 图1 图1 （b）为图像在 x 方向的投影分布 <math id="M12"><mi>r</mi><msub><mrow><mo stretchy="false">(</mo><mi>s</mi><mo stretchy="false">)</mo></mrow><mrow><mi>θ</mi><mo>=</mo><mn mathvariant="normal">0</mn></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19368080&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19368075&type= 11.93799973 4.48733330 ， 图1 图1 （c）为图像在 y 方向的投影分布 <math id="M13"><mi>r</mi><msub><mrow><mo stretchy="false">(</mo><mi>s</mi><mo stretchy="false">)</mo></mrow><mrow><mi>θ</mi><mo>=</mo><mi>π</mi><mo>/</mo><mn mathvariant="normal">2</mn></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19368090&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19368085&type= 13.88533401 4.48733330 。目标区域的能量分布与背景具有显著的差别，通过边缘检测 ［ 33 33 ］ 能够得到目标区域的边界点 <math id="M14"><mo stretchy="false">(</mo><msub><mrow><mi>x</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub><mo>,</mo><msub><mrow><mi>y</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19368100&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19368095&type= 11.09133339 4.74133301 <math id="M15"><mo stretchy="false">(</mo><msub><mrow><mi>x</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub><mo>,</mo><msub><mrow><mi>y</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19368110&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19368105&type= 11.09133339 4.74133301 <math id="M16"><mo stretchy="false">(</mo><msub><mrow><mi>x</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub><mo>,</mo><msub><mrow><mi>y</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19368120&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19368114&type= 11.09133339 4.74133301 <math id="M17"><mo stretchy="false">(</mo><msub><mrow><mi>x</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub><mo>,</mo><msub><mrow><mi>y</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19368129&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19368125&type= 11.09133339 4.74133301 ，目标区域的中心坐标为 <math id="M18"><mfenced separators="|"><mrow><mfrac><mrow><msub><mrow><mi>x</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub><mo>+</mo><msub><mrow><mi>x</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub></mrow><mrow><mn mathvariant="normal">2</mn></mrow></mfrac><mo>,</mo><mfrac><mrow><msub><mrow><mi>y</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub><mo>+</mo><msub><mrow><mi>y</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub></mrow><mrow><mn mathvariant="normal">2</mn></mrow></mfrac></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19368140&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19368135&type= 28.10933495 8.97466660 ，记为 <math id="M19"><mi>M</mi><mo stretchy="false">(</mo><msub><mrow><mi>x</mi></mrow><mrow><mi>r</mi></mrow></msub><mo>,</mo><msub><mrow><mi>y</mi></mrow><mrow><mi>r</mi></mrow></msub><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19368331&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19368328&type= 14.90133286 4.74133301 。 图1 图1 （d）为场景的Radon投影变换，红色曲线为四个边缘点的Radon投影结果，每一个点的投影轨迹都表现为正弦分布，并且目标区域Radon投影变换后，其能量分布集中在四条正弦曲线限制的范围内。 场景图像 f （ x ， y ）的大小为 N × N ，实现单像素成像的Radon投影变换基底模式为： <math display="block" id="M20"><mi>P</mi><mo stretchy="false">(</mo><mi>s</mi><mo>,</mo><mi>θ</mi><mo stretchy="false">)</mo><mo>=</mo><mfenced open="{" close="" separators="|"><mrow><mtable columnalign="left"><mtr><mtd><mn mathvariant="normal">1</mn><mo>,</mo><mi>x</mi><mi mathvariant="normal">c</mi><mi mathvariant="normal">o</mi><mi mathvariant="normal">s</mi><mi>θ</mi><mo>+</mo><mi>y</mi><mi mathvariant="normal">s</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">n</mi><mi>θ</mi><mo>-</mo><mi>s</mi><mo>=</mo><mn mathvariant="normal">0</mn></mtd></mtr><mtr><mtd><mn mathvariant="normal">0</mn><mo>,</mo><mi mathvariant="normal">e</mi><mi mathvariant="normal">s</mi><mi mathvariant="normal">l</mi><mi mathvariant="normal">e</mi></mtd></mtr></mtable></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19368159&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19368154&type= 56.72666931 9.56733322 . 其作用为提取图像在给定直线 <math id="M21"><mi>L</mi><mo stretchy="false">(</mo><mi>s</mi><mo>,</mo><mi>θ</mi><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19368170&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19368165&type= 10.92199993 3.55599999 上的所有像素点，遍历所有的 s 和 θ 即可获得图像Radon投影变换。而 <math id="M22"><mi>P</mi><mo stretchy="false">(</mo><mi>s</mi><mo>,</mo><mi>θ</mi><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19368180&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19368175&type= 11.17599964 3.55599999 在空间光调制器上表现为“亮”直线，编码方式虽然简单，但是每次仅获取一条直线上的像素点，SLM的运用效率极低，且获取的光强值极小。本文结合Hadamard变换，使用统计学原理获取图像在不同方向上的投影能量 ［ 31 31 ， 29 29 ］ 。Radon-Hadamard（R-H）基底为： <math display="block" id="M23"><mi>R</mi><msup><mrow><mi>H</mi></mrow><mrow><mo>+</mo></mrow></msup><msub><mrow><mo stretchy="false">(</mo><mi>s</mi><mo stretchy="false">)</mo></mrow><mrow><mi>θ</mi></mrow></msub><mo>=</mo><mfrac><mrow><mn mathvariant="normal">1</mn></mrow><mrow><mn mathvariant="normal">2</mn></mrow></mfrac><mo stretchy="false">[</mo><mn mathvariant="normal">1</mn><mo>+</mo><msup><mrow><mi>H</mi></mrow><mrow><mo>-</mo><mn mathvariant="normal">1</mn></mrow></msup><mfenced open="{" close="}" separators="|"><mrow><mi>δ</mi><mo stretchy="false">(</mo><mi>x</mi><mo>-</mo><mi>s</mi><mo stretchy="false">)</mo></mrow></mfenced><msub><mrow><mi>T</mi></mrow><mrow><mi>θ</mi></mrow></msub><mo stretchy="false">]</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19368190&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19368185&type= 61.97600174 7.45066643 ， <math display="block" id="M24"><mi>R</mi><msup><mrow><mi>H</mi></mrow><mrow><mo>-</mo></mrow></msup><msub><mrow><mo stretchy="false">(</mo><mi>s</mi><mo stretchy="false">)</mo></mrow><mrow><mi>θ</mi></mrow></msub><mo>=</mo><mfrac><mrow><mn mathvariant="normal">1</mn></mrow><mrow><mn mathvariant="normal">2</mn></mrow></mfrac><mo stretchy="false">[</mo><mn mathvariant="normal">1</mn><mo>-</mo><msup><mrow><mi>H</mi></mrow><mrow><mo>-</mo><mn mathvariant="normal">1</mn></mrow></msup><mfenced open="{" close="}" separators="|"><mrow><mi>δ</mi><mo stretchy="false">(</mo><mi>x</mi><mo>-</mo><mi>s</mi><mo stretchy="false">)</mo></mrow></mfenced><msub><mrow><mi>T</mi></mrow><mrow><mi>θ</mi></mrow></msub><mo stretchy="false">]</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19368201&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19368196&type= 61.97600174 7.45066643 ， 其中， <math id="M25"><msup><mrow><mi>H</mi></mrow><mrow><mo>-</mo><mn mathvariant="normal">1</mn></mrow></msup><mo stretchy="false">{</mo><mo stretchy="false">}</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19368210&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19368205&type= 9.99066734 3.80999994 表示哈达玛变换， <math id="M26"><msub><mrow><mi>T</mi></mrow><mrow><mi>θ</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19368221&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19368215&type= 3.72533321 3.72533321 表示旋转因子。Radon-Hadamard模式基底图案如 图2 图2 所示， 图2 图2 （a）（b）分别为获得图像在 x 和 y 方向的投影模式图案，选取 i =3阶Radon-Hadamard示例以简要显示。 使用Radon-Hadamard基底对场景图片的采样获得的测量值为： <math display="block" id="M27"><msup><mrow><mi>D</mi></mrow><mrow><mo>+</mo></mrow></msup><mo stretchy="false">(</mo><mi>s</mi><mo>,</mo><mi>θ</mi><mo stretchy="false">)</mo><mo>=</mo><mstyle displaystyle="true"><munderover><mo largeop="true">∑</mo><mrow><mi>x</mi><mo>=</mo><mn mathvariant="normal">1</mn></mrow><mrow><mi>N</mi></mrow></munderover></mstyle><mrow><mstyle displaystyle="true"><munderover><mo largeop="true">∑</mo><mrow><mi>y</mi><mo>=</mo><mn mathvariant="normal">1</mn></mrow><mrow><mi>N</mi></mrow></munderover></mstyle><mrow><mi>f</mi><mo stretchy="false">(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo stretchy="false">)</mo><mi>R</mi><msup><mrow><mi>H</mi></mrow><mrow><mo>+</mo></mrow></msup><msub><mrow><mo stretchy="false">(</mo><mi>s</mi><mo stretchy="false">)</mo></mrow><mrow><mi>θ</mi></mrow></msub></mrow></mrow></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19368243&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19368238&type= 55.96466827 8.89000034 ， <math display="block" id="M28"><msup><mrow><mi>D</mi></mrow><mrow><mo>-</mo></mrow></msup><mo stretchy="false">(</mo><mi>s</mi><mo>,</mo><mi>θ</mi><mo stretchy="false">)</mo><mo>=</mo><mstyle displaystyle="true"><munderover><mo largeop="true">∑</mo><mrow><mi>x</mi><mo>=</mo><mn mathvariant="normal">1</mn></mrow><mrow><mi>N</mi></mrow></munderover></mstyle><mrow><mstyle displaystyle="true"><munderover><mo largeop="true">∑</mo><mrow><mi>y</mi><mo>=</mo><mn mathvariant="normal">1</mn></mrow><mrow><mi>N</mi></mrow></munderover></mstyle><mrow><mi>f</mi><mo stretchy="false">(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo stretchy="false">)</mo><mi>R</mi><msup><mrow><mi>H</mi></mrow><mrow><mo>-</mo></mrow></msup><msub><mrow><mo stretchy="false">(</mo><mi>s</mi><mo stretchy="false">)</mo></mrow><mrow><mi>θ</mi></mrow></msub></mrow></mrow></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19368253&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19368248&type= 55.96466827 8.89000034 ， <math id="M29"><msup><mrow><mi>D</mi></mrow><mrow><mo>+</mo></mrow></msup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19368261&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19368257&type= 4.99533367 3.04800010 和 <math id="M30"><msup><mrow><mi>D</mi></mrow><mrow><mo>-</mo></mrow></msup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19368268&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19368265&type= 4.99533367 3.04800010 分别为单像素探测器获得的测量值，差分的测量方式能够有效的抑制噪声 ［ 24 24 ］ ，通过R-H模式获取的Radon投影变换域的系数为： <math display="block" id="M31"><mi>r</mi><mo stretchy="false">(</mo><mi>s</mi><mo>,</mo><mi>θ</mi><mo stretchy="false">)</mo><mo>=</mo><msup><mrow><mi>H</mi></mrow><mrow><mo>-</mo><mn mathvariant="normal">1</mn></mrow></msup><mo stretchy="false">{</mo><msup><mrow><mi>D</mi></mrow><mrow><mo>+</mo></mrow></msup><mo stretchy="false">(</mo><mi>s</mi><mo>,</mo><mi>θ</mi><mo stretchy="false">)</mo><mo>-</mo><msup><mrow><mi>D</mi></mrow><mrow><mo>-</mo></mrow></msup><mo stretchy="false">(</mo><mi>s</mi><mo>,</mo><mi>θ</mi><mo stretchy="false">)</mo><mo stretchy="false">}</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19368277&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19368273&type= 58.08132935 3.80999994 . 当 <math id="M32"><mi>θ</mi><mo>=</mo><mn mathvariant="normal">0</mn></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19368285&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19368281&type= 8.55133343 2.79399991 和 <math id="M33"><mi>θ</mi><mo>=</mo><mi mathvariant="normal">π</mi><mo>/</mo><mn mathvariant="normal">2</mn></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19368291&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19368288&type= 12.27666664 3.30200005 则能够获得 图1 图1 （b）和 图1 图1 （c）对应的投影曲线。考虑到Radon变换的周期特性，全采样需要的角度范围是 <math id="M34"><mo stretchy="false">[</mo><mn mathvariant="normal">0</mn><mo>,</mo><mi mathvariant="normal">π</mi><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19368300&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19368296&type= 8.89000034 3.55599999 ，因此对整个场景进行全采样所需的模式图案数量为 <math id="M35"><msub><mrow><mi>N</mi></mrow><mrow><mi>p</mi></mrow></msub><mo>=</mo><mn mathvariant="normal">2</mn><mo>×</mo><mi>m</mi><mo>×</mo><mi>N</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19368308&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19368303&type= 26.16200066 3.80999994 ， m 为投影方向的数量（ <math id="M36"><mi>m</mi><mo>=</mo><mi mathvariant="normal">π</mi><mo>/</mo><mi mathvariant="normal">Δ</mi><mi>θ</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19368317&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19368313&type= 15.66333294 3.30200005 ， <math id="M37"><mi mathvariant="normal">Δ</mi><mi>θ</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19368325&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19368321&type= 4.23333359 2.79399991 为角度的采样间隔）。 如 图3 图3 （a）所示，ARSI通过投影定位出目标区域，仅对目标区域进行单像素采样和成像，能够很大程度降低测量次数。如 图3 图3 （b）所示，我们计算目标区域的中心坐标为 <math id="M38"><mi>M</mi><mo stretchy="false">(</mo><msub><mrow><mi>x</mi></mrow><mrow><mi>r</mi></mrow></msub><mo>,</mo><msub><mrow><mi>y</mi></mrow><mrow><mi>r</mi></mrow></msub><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19368331&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19368328&type= 14.90133286 4.74133301 ，选取与目标区域大小接近的Radon-Hadamard模式作为采样图案，新的R-H模式大小为 <math id="M39"><mi>n</mi><mo>×</mo><mi>n</mi><mo stretchy="false">(</mo><mi>n</mi><mo>=</mo><msup><mrow><mn mathvariant="normal">2</mn></mrow><mrow><mi>i</mi></mrow></msup><mo>,</mo><mi>n</mi><mo>≤</mo><mi>N</mi><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19368337&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19368334&type= 33.02000046 3.89466691 ，小于场景分辨率 N × N 。更新后的目标区域大于原始目标区域，降低了对初始目标区域定位的精度要求。此时所单像素成像需要的采样数量为 <math id="M40"><msub><mrow><mi>N</mi></mrow><mrow><mi>p</mi></mrow></msub><mo>=</mo><mn mathvariant="normal">2</mn><mo>×</mo><mi>m</mi><mo>×</mo><mi>n</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19368344&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19368341&type= 25.14599991 3.80999994 。R-H模式变为： <math display="block" id="M41"><mi>R</mi><mi>H</mi><msub><mrow><mo stretchy="false">(</mo><mi>s</mi><mo stretchy="false">)</mo></mrow><mrow><mi>θ</mi></mrow></msub><mo>→</mo><mi>R</mi><mi>H</mi><msub><mrow><mo stretchy="false">(</mo><mi>s</mi><mo>_</mo><mi>r</mi><mo stretchy="false">)</mo></mrow><mrow><mi>θ</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19368350&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19368347&type= 33.86666870 4.48733330 . 目标区域的Radon变换域为： <math display="block" id="M42"><mi>r</mi><mo stretchy="false">(</mo><mi>s</mi><mo>,</mo><mi>θ</mi><mo stretchy="false">)</mo><mo>→</mo><mi>r</mi><mo stretchy="false">(</mo><mi>s</mi><mo>_</mo><mi>r</mi><mo>,</mo><mi>θ</mi><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19368365&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19368363&type= 28.70199966 3.55599999 ， 其中 <math id="M43"><mi>s</mi><mo>_</mo><mi>r</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19368371&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19368368&type= 4.48733330 3.30200005 为目标区域的投影坐标，更新后的目标区域大小为 n×n ，中心坐标不变，此时目标区域的范围变为： <math display="block" id="M44"><mi mathvariant="normal">Ω</mi><mo stretchy="false">(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo stretchy="false">)</mo><mo>,</mo><mtext> </mtext><mfenced open="{" close="" separators="|"><mrow><mtable columnalign="left"><mtr><mtd><msub><mrow><mi>x</mi></mrow><mrow><mi>r</mi></mrow></msub><mo>-</mo><mi>n</mi><mo>/</mo><mn mathvariant="normal">2</mn><mo>≤</mo><mi>x</mi><mo>≤</mo><msub><mrow><mi>x</mi></mrow><mrow><mi>r</mi></mrow></msub><mo>+</mo><mi>n</mi><mo>/</mo><mn mathvariant="normal">2</mn><mo>-</mo><mn mathvariant="normal">1</mn></mtd></mtr><mtr><mtd><msub><mrow><mi>y</mi></mrow><mrow><mi>r</mi></mrow></msub><mo>-</mo><mi>n</mi><mo>/</mo><mn mathvariant="normal">2</mn><mo>≤</mo><mi>y</mi><mo>≤</mo><msub><mrow><mi>y</mi></mrow><mrow><mi>r</mi></mrow></msub><mo>+</mo><mi>n</mi><mo>/</mo><mn mathvariant="normal">2</mn><mo>-</mo><mn mathvariant="normal">1</mn></mtd></mtr></mtable></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19368378&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19368374&type= 60.45199966 10.49866676 . 图4 图4 为ARSI对目标区域 <math id="M45"><mi mathvariant="normal">Ω</mi><mo stretchy="false">(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19368386&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19368383&type= 11.76866722 3.80999994 的采样过程，通过单像素成像方法，使用 图4 图4 （c）中不同的R-H模式图案调制目标区域，能够获取如 图4 图4 （d）所示的Radon变换域。 我们定义目标区域和场景区域的采样压缩比为： <math display="block" id="M46"><mi>C</mi><mo>=</mo><mfrac><mrow><mi>t</mi><mo>×</mo><mi>n</mi></mrow><mrow><mi>N</mi></mrow></mfrac><mo>×</mo><mn mathvariant="normal">100</mn><mi mathvariant="normal">%</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19368402&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19368399&type= 30.98800087 7.61999989 ， 其中， t 为目标个数。 2.2　图像重建 2 滤波反投影算法用于重建图像，其多用于计算机断层成像领域，能够将图像信息从Radon变换域转为空间域，其表达式为： <math display="block" id="M47"><mtable columnalign="left"><mtr><mtd><mi>f</mi><mo stretchy="false">(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo stretchy="false">)</mo><mo>=</mo><msubsup><mo>∫</mo><mrow><mn mathvariant="normal">0</mn></mrow><mrow><mi>π</mi></mrow></msubsup><mrow><mi>r</mi><mo stretchy="false">(</mo><mi>s</mi><mo>,</mo><mi>θ</mi><mo stretchy="false">)</mo><mi mathvariant="normal">*</mi><mi>h</mi><mo stretchy="false">(</mo><mi>s</mi><mo stretchy="false">)</mo><mi mathvariant="normal">d</mi><mi>θ</mi></mrow><mo>=</mo></mtd></mtr><mtr><mtd><msubsup><mo>∫</mo><mrow><mn mathvariant="normal">0</mn></mrow><mrow><mi>π</mi></mrow></msubsup><mrow><mi>r</mi><mo stretchy="false">(</mo><mi>x</mi><mi mathvariant="normal">c</mi><mi mathvariant="normal">o</mi><mi mathvariant="normal">s</mi><mi>θ</mi><mo>+</mo><mi>y</mi><mi mathvariant="normal">s</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">n</mi><mi>θ</mi><mo stretchy="false">)</mo><mi mathvariant="normal">*</mi><mi>h</mi><mo stretchy="false">(</mo><mi>s</mi><mo stretchy="false">)</mo><mi mathvariant="normal">d</mi><mi>θ</mi></mrow></mtd></mtr></mtable></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19368408&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19368405&type= 48.42933273 15.66333294 ， <math id="M48"><mi>r</mi><mo stretchy="false">(</mo><mi>s</mi><mo>,</mo><mi>θ</mi><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19368414&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19368411&type= 10.32933331 3.55599999 为单像素方法获取的Radon变换域， <math id="M49"><mi>h</mi><mo stretchy="false">(</mo><mi>s</mi><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19368427&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19368423&type= 6.85799980 3.55599999 为滤波器函数，当不存在 <math id="M50"><mi>h</mi><mo stretchy="false">(</mo><mi>s</mi><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19368427&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19368423&type= 6.85799980 3.55599999 时， 公式（15） 公式（15） 为直接反投影的过程， <math id="M51"><mi>f</mi><mo stretchy="false">(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19368432&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19368430&type= 11.43000031 3.80999994 中会存在星状伪影，造成图像质量的下降。本文采用Hamming滤波函数，避免了大量的正弦、余弦计算，所重建图像轮廓清楚。目标区域的单像素成像结果为 <math id="M52"><msub><mrow><mi>f</mi></mrow><mrow><mi>r</mi></mrow></msub><mo stretchy="false">(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19368439&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19368436&type= 11.85333347 4.74133301 ： <math display="block" id="M53"><msub><mrow><mi>f</mi></mrow><mrow><mi>r</mi></mrow></msub><mo stretchy="false">(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo stretchy="false">)</mo><mo>=</mo><msubsup><mo>∫</mo><mrow><mn mathvariant="normal">0</mn></mrow><mrow><mi>π</mi></mrow></msubsup><mrow><mi>r</mi><mo stretchy="false">(</mo><mi>s</mi><mo>_</mo><mi>r</mi><mo>,</mo><mi>θ</mi><mo stretchy="false">)</mo><mi mathvariant="normal">*</mi><mi mathvariant="normal">h</mi><mo stretchy="false">(</mo><mi>s</mi><mo>_</mo><mi>r</mi><mo stretchy="false">)</mo><mi mathvariant="normal">d</mi><mi>θ</mi></mrow></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19368444&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19368442&type= 50.79999924 6.94266701 . 为获得还原目标在场景中的位置，将目标区域与背景融合，使最终的重建场景与原场景具有相同的分辨率（N×N），坐标转换公式为： <math display="block" id="M54"><mi>f</mi><mo stretchy="false">(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo stretchy="false">)</mo><mo>=</mo><mfenced open="{" close="" separators="|"><mrow><mtable columnalign="left"><mtr><mtd><msub><mrow><mi>f</mi></mrow><mrow><mi>r</mi></mrow></msub><mo>,</mo><mo stretchy="false">(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo stretchy="false">)</mo><mo>∈</mo><mi mathvariant="normal">Ω</mi></mtd></mtr><mtr><mtd><mn mathvariant="normal">0</mn><mo>,</mo><mi mathvariant="normal">e</mi><mi mathvariant="normal">l</mi><mi mathvariant="normal">s</mi><mi mathvariant="normal">e</mi></mtd></mtr></mtable></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19368451&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19368447&type= 36.49133301 9.65200043 . 3 成像分析 1 图5 图5 为ARSI方法， 图5 图5 （a）为ARSI方法的组成结构，主要包括单像素成像硬件和ARSI算法； 图5 图5 （b）为实现ARSI的光路结构，光路类型为前调制型单像素成像光路 ［ 20 20 ］ ，选用数字投影仪器（Digital Projector， DLP）作为SLM用于产生编码图案，DLP生成时间变化的编码R-H模式序列用于照明场景，生成速率5 000 Hz，反射光经过焦距为125 mm的聚光透镜，光电二极管作为单像素探测器，数据采集设备（DAQ）用于数据采集和信号转换，计算机用于控制和处理数据，物体距离透镜0.5 m。 3.1　单目标ARSI 2 我们继续使用“小牛”图验证ARSI算法，重建场景分辨率为256×256，如 图6 图6 所示，根据投影定位算法， 图6 图6 （a）为计算出目标区域的边界，其中心坐标为M（129，168），选取64×64的R-H模式图案，因此更新后的目标区域如 图6 图6 （b）所示，采样压缩比为25%。 R-H模式基底采样时，设置角度间隔为1°，因此 m =180。通过一系列不同的R-H模式图案对目标区域进行采样，以获取Radon变换域系数，再通过滤波反投影方法重建目标区域图像，通过图像融合最终获得重建后的场景，如 图7 图7 所示，添加结构相似性系数以客观评价图像重建结果 ［ 35 35 ］ ，以研究重建目标区域与原始目标区域之间的差别。 ARSI的采样次数为23 040，RSI的采样次数为92 160，ARSI的采样次数为RSI的1/4，成像效率提升4倍，但最终的成像质量相同，证明了ARSI在降低采样数量的同时，还能够保证成像的质量。实验对真实“小鸭子”场景进行ARSI，重建场景分辨率为256×256。根据ARSI定位，128×128的Radon-Hadamard的模板对目标区域进行采样，采样压缩比为 C =50%，角度间隔为1°，实验结果如 图8 图8 所示，对于大目标区域的重建依旧表现出较高的成像质量，成像效率是全采样条件下2倍。 实验还设置了相同采样次数的条件下，RSI和ARSI成像质量的对比，分别设置 n =32，64，128，256大小的R-H模式图案执行RSI和ARSI，实验结果如 图8 图8 所示。当R-H模板仅为32×32时，其小于目标区域的大小，矩形区域放大后显示，目标区域大部分特征得以恢复但边缘信息丢失，ARSI重建的目标区域不完整，此为R-H模式大小不匹配的重建结果。此时RSI重建质量非常差，无法识别目标区域特征信息。随着R-H模式大小匹配目标区域，R-H模式分辨率小于图像分辨率时，RSI无法实现全采样，重建图像的质量明显低于ARSI，实验证明了当压缩率小于100%时，ARSI有效提高了单像素成像的成像效率。 3.2　多目标ARSI 2 在已知具有三个目标区域的情况下，使用ARSI寻找目标区域的位置坐标。如 图10 图10 所示，设置了三种不同位置分布的场景，重建场景分辨率为256×256，分别为：（a）三个“小牛”的横纵坐标互不混叠，在 x 方向和 y 方向上的投影具有三个明显的特征峰分布，因此能够准确定位出三个目标区域的位置；在 图10 图10 （b）中，“小牛2”和“小牛3”横坐标重叠，其投影的特征峰发生混叠，三个“小牛”在 x 方向上仅表现出两个特征峰，但在 y 方向上的投影互不影响，明显区分其坐标位置，因此依然能够定位三个目标位置； 图10 图10 （c）中，“小牛1”和“小牛2”纵坐标重叠，“小牛2”和“小牛3”横坐标重叠，在 x 和 y 方向上都表现出两个特征峰，尤其是peak-2其投影的能量值大于peak-1，且peak-2的宽度大于peak-1，这表示peak-2为两个目标合成的特征峰，但无法区分二者的位置，我们添加了 y = x 方向的投影分布，在该方向上能够准确区分三个特征峰，借助三个方向的投影曲线，能够精准定位三个目标区域的中心坐标。 图10 图10 （c）定位所需的R-H基底的采样次数为1 536次，使用SLM调制花费极少的时间。 图10 图10 （c）为多目标区域的重建结果，R-H模板大小为64×64，分别对三个区域进行成像，然后融合背景信息还原场景，全采样所需的采样数量为92 160，ARSI采样次数为 <math id="M55"><mn mathvariant="normal">3</mn><mo>×</mo><mn mathvariant="normal">2</mn><mo>×</mo><mn mathvariant="normal">64</mn><mo>×</mo><mn mathvariant="normal">180</mn><mo>=</mo><mn mathvariant="normal">69</mn><mtext> </mtext><mn mathvariant="normal">120</mn></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19368511&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19368507&type= 41.74066544 2.79399991 ，采样压缩比 C =75%。该实验验证了ARSI具有对场景内多目标图像定位和重建的能力，其所需的采样数量少于全采样条件下的采样数量。 3.3　局限性 2 本文所提出的方法也具有局限性，实验场景的背景信息过于简单，当目标区域位于复杂的场景中时，其定位难度将增加；Hadamard变换的对象只能为2 i 大小，限制了目标区域范围，但如果不考虑探测效率，使用高灵敏度的单像素探测器，可以仅通过线编码获取Radon变换域；对于多个目标区域的定位，目标区域之间距离不能过小，否则会造成解纠缠困难，产生定位误差。ARSI在多目标区域定位中具有潜在能力，通过更改Radon-Hadamard模式尺寸使其应用于场景的多分辨率成像，未来结合深度学习提升其在复杂目标和场景下的应用能力。 4 结　论 1 本文提出使用Radon投影变换实现目标区域定位和成像的技术。在有效信息仅占场景一部分的情况下，使用水平和垂直的Radon-Hadamard模式作为采样的基底图案，定位目标区域的中心位置和计算其大小范围。选用匹配分辨率的基底图案对目标区域进行重建，ARSI所需的采样数量远小于全采样的采样数量，且能够保证重建图像质量，提高了单像素成像的成像效率。另外尝试了多目标实验，通过增加投影方向，ARSI能够定位并成像多个目标区域，可应用于多目标探测领域。 参考文献 HADFIELD R H . Single-photon detectors for optical quantum information applications ［J］. Nature Photonics ， 2009 ， 3 （ 12 ）： 696 - 705 . Single-photon detectors for optical quantum information applications GODA K ， TSIA K K ， JALALI B . Serial time-encoded amplified imaging for real-time observation of fast dynamic phenomena ［J］. Nature ， 2009 ， 458 （ 7242 ）： 1145 - 1149 . Serial time-encoded amplified imaging for real-time observation of fast dynamic phenomena RYCZKOWSKI P ， BARBIER M ， FRIBERG A T ， et al . Ghost imaging in the time domain ［J］. Nature Photonics ， 2016 ， 10 （ 3 ）： 167 - 170 . Ghost imaging in the time domain OTA S ， HORISAKI R ， KAWAMURA Y ， et al . Ghost cytometry ［J］. Science ， 2018 ， 360 （ 6394 ）： 1246 - 1251 . Ghost cytometry XU Z H ， CHEN W ， PENUELAS J ， et al . 1000 fps computational ghost imaging using LED-based structured illumination ［J］. Optics Express ， 2018 ， 26 （ 3 ）： 2427 - 2434 . 1000 fps computational ghost imaging using LED-based structured illumination PITTMAN T B ， SHIH Y H ， STREKALOV D V ， et al . Optical imaging by means of two-photon quantum entanglement ［J］. Physical Review A ， 1995 ， 52 （ 5 ）： R3429 - R3432 . Optical imaging by means of two-photon quantum entanglement BENNINK R S ， BENTLEY S J ， BOYD R W . “ Two-Photon” coincidence imaging with a classical source ［J］. Physical Review Letters ， 2002 ， 89 （ 11 ）： 113601 . Two-Photon” coincidence imaging with a classical source SUN B ， EDGAR M P ， BOWMAN R ， et al . Differential computational ghost imaging ［C］. Imaging and Applied Optics. Arlington ， Virginia. Washington， D.C.： OSA ， 2013 ： CTu1C.4 . Differential computational ghost imaging SUN B Q ， WELSH S S ， EDGAR M P ， et al . Normalized ghost imaging ［J］. Optics Express ， 2012 ， 20 （ 15 ）： 16892 - 16901 . Normalized ghost imaging WELSH S S ， EDGAR M P ， BOWMAN R ， et al . Fast full-color computational imaging with single-pixel detectors ［J］. Optics Express ， 2013 ， 21 （ 20 ）： 23068 - 23074 . Fast full-color computational imaging with single-pixel detectors SUN B ， EDGAR M P ， BOWMAN R ， et al . 3D computational imaging with single-pixel detectors ［J］. Science ， 2013 ， 340 （ 6134 ）： 844 - 847 . 3D computational imaging with single-pixel detectors DONOHO D L . Compressed sensing ［J］. IEEE Transactions on Information Theory ， 2006 ， 52 （ 4 ）： 1289 - 1306 . Compressed sensing KATZ O ， BROMBERG Y ， SILBERBERG Y . Compressive ghost imaging ［J］. Applied Physics Letters ， 2009 ， 95 （ 13 ）： 131110 . Compressive ghost imaging TAKHAR D ， LASKA J N ， WAKIN M B ， et al . A new compressive imaging camera architecture using optical-domain compression ［C］. Proc SPIE 6065 ， Computational Imaging IV ， 2006 ， 6065 ： 606509 . A new compressive imaging camera architecture using optical-domain compression GAO Z R ， DING L X ， XIONG Q M ， et al . Image compressive sensing reconstruction based on z-score standardized group sparse representation ［J］. IEEE Access ， 2019 ， 7 ： 90640 - 90651 . Image compressive sensing reconstruction based on z-score standardized group sparse representation 彭进业 ， 金浩强 ， 石剑虹 ， 等 . 高速单像素相机数据采集系统 ［J］. 光学 精密工程 ， 2014 ， 22 （ 4 ）： 837 - 843 . 高速单像素相机数据采集系统 PENG J Y ， JIN H Q ， SHI J H ， et al . Data acquisition system for high speed single-pixel camera ［J］. Opt. Precision Eng. ， 2014 ， 22 （ 4 ）： 837 - 843 . （in Chinese） Data acquisition system for high speed single-pixel camera PHILLIPS D B ， SUN M J ， TAYLOR J M ， et al . Adaptive foveated single-pixel imaging with dynamic supersampling ［J］. Science Advances ， 2017 ， 3 （ 4 ）： e1601782 . Adaptive foveated single-pixel imaging with dynamic supersampling ROUSSET F ， DUCROS N ， FARINA A ， et al . Adaptive basis scan by wavelet prediction for single-pixel imaging ［J］. IEEE Transactions on Computational Imaging ， 2017 ， 3 （ 1 ）： 36 - 46 . Adaptive basis scan by wavelet prediction for single-pixel imaging LIU B L ， YANG Z H ， LIU X ， et al . Coloured computational imaging with single-pixel detectors based on a 2D discrete cosine transform ［J］. Journal of Modern Optics ， 2017 ， 64 （ 3 ）： 259 - 264 . Coloured computational imaging with single-pixel detectors based on a 2D discrete cosine transform 张子邦 ， 陆天傲 ， 彭军政 ， 等 . 傅里叶单像素成像技术与应用 ［J］. 红外与激光工程 ， 2019 ， 48 （ 6 ）： 22 - 40 . 傅里叶单像素成像技术与应用 ZHANG Z B ， LU T A ， PENG J ZH ， et al . Fourier single-pixel imaging techniques and applications ［J］. Infrared and Laser Engineering ， 2019 ， 48 （ 6 ）： 22 - 40 . （in Chinese） Fourier single-pixel imaging techniques and applications ZHANG Z B ， MA X ， ZHONG J G . Single-pixel imaging by means of Fourier spectrum acquisition ［J］. Nature Communications ， 2015 ， 6 ： 6225 . Single-pixel imaging by means of Fourier spectrum acquisition ZHANG Z B ， WANG X Y ， ZHENG G A ， et al . Fast Fourier single-pixel imaging via binary illumination ［J］. Scientific Reports ， 2017 ， 7 ： 12029 . Fast Fourier single-pixel imaging via binary illumination EDGAR M P ， GIBSON G M ， BOWMAN R W ， et al . Simultaneous real-time visible and infrared video with single-pixel detectors ［J］. Scientific Reports ， 2015 ， 5 ： 10669 . Simultaneous real-time visible and infrared video with single-pixel detectors ZHANG Z B ， WANG X Y ， ZHENG G A ， et al . Hadamard single-pixel imaging versus Fourier single-pixel imaging ［J］. Optics Express ， 2017 ， 25 （ 16 ）： 19619 - 19639 . Hadamard single-pixel imaging versus Fourier single-pixel imaging SUN M J ， EDGAR M P ， GIBSON G M ， et al . Single-pixel three-dimensional imaging with time-based depth resolution ［J］. Nature Communications ， 2016 ， 7 ： 12010 . Single-pixel three-dimensional imaging with time-based depth resolution BIAN L H ， SUO J L ， HU X M ， et al . Efficient single pixel imaging in Fourier space ［J］. Journal of Optics ， 2016 ， 18 （ 8 ）： 085704 . Efficient single pixel imaging in Fourier space WENWEN M ， DONGFENG S ， JIAN H ， et al . Sparse Fourier single-pixel imaging ［J］. Optics Express ， 2019 ， 27 （ 22 ）： 31490 - 31503 . Sparse Fourier single-pixel imaging JIANG H Z ， ZHU S G ， ZHAO H J ， et al . Adaptive regional single-pixel imaging based on the Fourier slice theorem ［J］. Optics Express ， 2017 ， 25 （ 13 ）： 15118 - 15130 . Adaptive regional single-pixel imaging based on the Fourier slice theorem SHI D F ， YIN K X ， HUANG J ， et al . Fast tracking of moving objects using single-pixel imaging ［J］. Optics Communications ， 2019 ， 440 ： 155 - 162 . Fast tracking of moving objects using single-pixel imaging HOVLAND H . Tomographic scanning imager ［J］. Optics Express ， 2009 ， 17 （ 14 ）： 11371 - 11387 . Tomographic scanning imager SHI D F ， HUANG J ， MENG W W ， et al . Radon single-pixel imaging with projective sampling ［J］. Optics Express ， 2019 ， 27 （ 10 ）： 14594 . Radon single-pixel imaging with projective sampling Radon J . über die bestimmung von funktionen durch ihre integralwerte längs gewisser mannigfaltigkeiten ［J］. Classic papers in modern diagnostic radiology ， 2005 ， 5 ： 21 . über die bestimmung von funktionen durch ihre integralwerte längs gewisser mannigfaltigkeiten 吴晓波 . 图像边缘特征分析 ［J］. 光学 精密工程 ， 1999 ， 7 （ 1 ）： 59 - 63 图像边缘特征分析 WU X B . Analyzing of the characteristic of image edge ［J］. Opt. Precision Eng. ， 1999 ， 7 （ 1 ）： 59 - 63 . （in Chinese） Analyzing of the characteristic of image edge 王冬云 ， 唐楚 ， 鄂世举 ， 等 . 基于导向滤波Retinex和自适应Canny的图像边缘检测 ［J］. 光学 精密工程 ， 2021 ， 29 （ 2 ）： 443 - 451 . 基于导向滤波Retinex和自适应Canny的图像边缘检测 WANG D Y ， TANG CH ， E SH J ， et al . Image edge detection based on guided filter Retinex and adaptive Canny ［J］. Opt. Precision Eng. ， 2021 ， 29 （ 2 ）： 443 - 451 . （in Chinese） Image edge detection based on guided filter Retinex and adaptive Canny HORÉ A ， ZIOU D . Image quality metrics： PSNR vs . SSIM［C］. 2010 20th International Conference on Pattern Recognition . 2326，2010 ， Istanbul， Turkey . IEEE ， 2010 ： 2366 - 2369 . Image quality metrics： PSNR vs 10.37188/OPE.20212908.1976 TP391 TB81 A 1004-924X（2021）08-1976-09 王蔚松,吴洪波,王灵杰等.自适应Radon单像素成像[J].光学精密工程,2021,29(08):1976-1984. 王蔚松,吴洪波,王灵杰等.自适应Radon单像素成像[J].光学精密工程,2021,29(08):1976-1984. WANG Wei-song,WU Hong-bo,WANG Ling-jie,et al.Adaptive radon single-pixel imaging method[J].Optics and Precision Engineering,2021,29(08):1976-1984. WANG Wei-song,WU Hong-bo,WANG Ling-jie,et al.Adaptive radon single-pixel imaging method[J].Optics and Precision Engineering,2021,29(08):1976-1984. 国家自然科学基金资助项目（No.61675100），国家自然科学基金（No.62005271） 2021-03-03 2021-03-23 2021-08-15 2021-09-01 光学精密工程 8 29