CN114721006A - 高精度单像素成像方法及系统 - Google Patents

高精度单像素成像方法及系统 Download PDF

Info

Publication number
CN114721006A
CN114721006A CN202210641951.1A CN202210641951A CN114721006A CN 114721006 A CN114721006 A CN 114721006A CN 202210641951 A CN202210641951 A CN 202210641951A CN 114721006 A CN114721006 A CN 114721006A
Authority
CN
China
Prior art keywords
target object
light
chebyshev
different orders
dimensional
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202210641951.1A
Other languages
English (en)
Other versions
CN114721006B (zh
Inventor
韩凯
来文昌
崔文达
孟琪
雷国忠
王彦
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
National University of Defense Technology
Original Assignee
National University of Defense Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by National University of Defense Technology filed Critical National University of Defense Technology
Priority to CN202210641951.1A priority Critical patent/CN114721006B/zh
Publication of CN114721006A publication Critical patent/CN114721006A/zh
Application granted granted Critical
Publication of CN114721006B publication Critical patent/CN114721006B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S17/00—Systems using the reflection or reradiation of electromagnetic waves other than radio waves, e.g. lidar systems
    • G01S17/88—Lidar systems specially adapted for specific applications
    • G01S17/89—Lidar systems specially adapted for specific applications for mapping or imaging

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Electromagnetism (AREA)
  • General Physics & Mathematics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Investigating Or Analysing Materials By Optical Means (AREA)
  • Optical Radar Systems And Details Thereof (AREA)

Abstract

本发明提供一种高精度单像素成像方法及系统，包括：根据目标物体，生成一系列不同阶数的切比雪夫多项式分布的二维光场并加载到空间光调制器上；空间光调制器利用所述一系列不同阶数的切比雪夫多项式分布的二维光场对光源发出的激光进行调制，调制后的激光照射到目标物体上，采集目标物体的反射回光；基于所述反射回光以及所述一系列不同阶数的切比雪夫多项式分布的二维光场，对目标物体图像进行重构。本发明利用切比雪夫多项式不同阶数间的正交性，可以完整重构出物体图像，与散斑光场相比，该方案可大幅减少测量次数，且成像质量较高。

Description

高精度单像素成像方法及系统
技术领域
本发明涉及光学探测成像技术领域，具体是一种采用正交矩光场的高精度单像素成像方法及系统。
背景技术
传统的光学成像广泛采用以硅为光敏材料的阵列像素探测器(如CCD和CMOS)，硅的光谱响应范围集中在可见光及其附近波段（约400-800 nm）。对于非可见光波段的成像（如红外波段、太赫兹、X光波段等），目前市场上还未见报道的具有高分辨率且低成本的阵列探测器。而制作工作在非可见波段的单像素探测器的成本低得多，因此单像素探测成像技术可有效解决非可见光波段成像成本高的难题，同时，单像素探测器具有更高的光强灵敏度、更快的响应速度，因此在医学成像、远距离目标探测等弱光探测领域具有明显的优势。
单像素成像技术使用无空间分辨能力的单点探测器（如光电倍增管）对目标进行时域上的多次采样。与传统扫描式单像素成像采用局部单点采样方式不同，现在的单像素成像技术通常指基于空间光调制的单像素成像，采用全局采样方式对物体空间信息进行编码，再通过单像素探测器对编码后的信息依次进行采集，最后通过计算对空间信息进行解码，重构物体图像。
单像素成像技术最早源于量子鬼成像技术，鬼成像技术通常采用随机图案作为空间光调制图案，最常见的方式是将激光通过旋转的毛玻璃产生随机散斑光场，通过随机散斑对目标进行强度调制，记录随机散斑场的分布和探测相应的回光强度，可重构物体图像。然而基于随机散斑场获取目标信息需要大量的测量次数才能有效重构物体图案，数据采集时间长，成像效率较低。后来压缩感知技术被应用于鬼成像技术可有效减少数据采集时间。然而，压缩感知技术的算法运行时间往往耗时较长，且成像质量取决于物体的稀疏性，难以实时获取高质量物体图像。
发明内容
针对目前单像素成像技术中数据采集时间较长、成像质量不高的问题，本发明提出了一种高精度单像素成像方法及系统。
为实现上述技术目的，本发明提出的技术方案为：
一方面，本发明提供一种高精度单像素成像方法，包括：
(S1) 根据目标物体，生成一系列不同阶数的切比雪夫多项式分布的二维光场并加载到空间光调制器上；
(S2) 空间光调制器利用所述一系列不同阶数的切比雪夫多项式分布的二维光场对光源发出的激光进行调制，调制后的激光照射到目标物体上，采集目标物体的反射回光；
(S3) 基于所述反射回光以及所述一系列不同阶数的切比雪夫多项式分布的二维光场，对目标物体图像进行重构。
进一步地，所述步骤(S1)包括：
目标物体离散化大小为M*N像素的矩阵，目标物体的实际大小为
Figure 577988DEST_PATH_IMAGE001
的矩形，其中M、N为正整数，
Figure 524210DEST_PATH_IMAGE002
、
Figure 219633DEST_PATH_IMAGE003
分别为一个像素在y，x方向的几何尺寸；
生成一系列不同阶数的切比雪夫多项式分布的二维光场Q(x, y, n, m)，其中x，y为x方向和y方向的像素点坐标，n和m分别为x方向和y方向的切比雪夫多项式的阶数，x=1,2,...,N-1；n=1,2,...,N-1；y=1,2,...,M-1；m=1,2,...,M-1；
将不同阶数(n, m)的二维光场Q(x, y, n, m)拆分为光场Q + (x, y, n, m)和光场Q - (x, y, n, m)后依次加载在空间光调制器上，其中
Figure 930101DEST_PATH_IMAGE004
，Q + (x, y, n, m)保留Q(x, y, n, m)的正数元素，其余元素用零补齐，Q - (x, y, n, m)保留Q(x, y, n, m)中的负数元素并在负数元素相应位置取其绝对值，其余元素用零补齐。
进一步地，所述步骤(S3)包括：
利用所述反射回光计算目标物体图像的切比雪夫矩；
基于目标物体图像的切比雪夫矩以及所述一系列不同阶数的切比雪夫多项式分布的二维光场，重构出目标物体图像。
进一步地，所述步骤(S3)中目标物体图像的切比雪夫矩为：
Figure 674066DEST_PATH_IMAGE005
其中
Figure 770066DEST_PATH_IMAGE006
为光场Q + (x, y, n, m)照射目标物体后产生的反射回光，
Figure 574074DEST_PATH_IMAGE007
为光场Q - (x, y, n, m)照射目标物体后产生的反射回光；
进一步地，所述步骤(S3)中重构出的目标物体图像f(x,y)为：
Figure 568575DEST_PATH_IMAGE008
。
另一方面，本发明提供一种高精度单像素成像系统，包括：
空间光调制图案生成模块，用于根据目标物体，生成一系列不同阶数的切比雪夫多项式分布的二维光场并加载到空间光调制器上；
光源，用于产生激光，并将产生的激光投影照射到空间光调制器上；
空间光调制器，利用所述一系列不同阶数的切比雪夫多项式分布的二维光场对光源发出的激光进行调制，调制后的激光照射到目标物体上；
单像素光电探测器，用于采集目标物体的反射回光；
成像模块，用于基于所述反射回光以及所述一系列不同阶数的切比雪夫多项式分布的二维光场，对目标物体图像进行重构。
作为优选实施方案，所述空间光调制图案生成模块和成像模块由计算机实现。
作为优选实施方案，所述空间光调制器可以是数字微镜阵列（DMD）或者空间光相位调制器。
作为优选实施方案，所述光源包括激光器以及准直扩束系统。
与现有技术相比，本发明的优点在于：
本发明采用具有确定数学函数表达式的满足切比雪夫多项式分布的二维结构光场替代鬼成像技术中的随机散斑光场，在此基础上重构物体图像，该方案利用切比雪夫多项式不同阶数间的正交性，可以完整重构出物体图像，与鬼成像技术中的随机散斑光场相比，该方案可大幅减少测量次数，且成像质量较高。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图示出的结构获得其他的附图。
图1是本发明一实施例的流程图；
图2是本发明一实施例的结构图；
图3为本发明一实施例采用的不同阶数的切比雪夫多项式分布的二维光场；
图4为本发明一实施例中单像素探测器探测到的不同阶数的切比雪夫矩；
图5是本发明一实施例中重建的目标物体图像；
图6是本发明一实施例中在欠采样条件下重建的目标物体图像。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚明白，下面将以附图及详细叙述来清楚说明本发明所揭示内容的精神，任何所属技术领域技术人员在了解本发明内容的实施例后，当可由本发明内容所教示的技术，加以改变及修饰，其并不脱离本发明内容的精神与范围。本发明的示意性实施例及其说明用于解释本发明，但并不作为对本发明的限定。
在一实施例中，参照图1，提供一种高精度单像素成像方法，包括：
(S1) 根据目标物体，生成一系列不同阶数的切比雪夫多项式分布的二维光场并加载到空间光调制器上；
(S2) 空间光调制器利用所述一系列不同阶数的切比雪夫多项式分布的二维光场对光源发出的激光进行调制，调制后的激光照射到目标物体上，采集目标物体的反射回光；
(S3) 基于所述反射回光以及所述一系列不同阶数的切比雪夫多项式分布的二维光场，对目标物体图像进行重构。
在一实施例中，所述步骤(S1)包括：
目标物体离散化大小为M*N像素的矩阵，目标物体的实际大小为
Figure 116231DEST_PATH_IMAGE009
的矩形，其中M、N为正整数，
Figure 83050DEST_PATH_IMAGE010
、
Figure 74271DEST_PATH_IMAGE011
分别为一个像素在y，x方向的几何尺寸；
生成一系列不同阶数的切比雪夫多项式分布的二维光场Q(x, y, n, m)，其中x，y为x方向和y方向的像素点坐标，n和m分别为x方向和y方向的切比雪夫多项式的阶数，x=1,2,...,N-1；n=1,2,...,N-1；y=1,2,...,M-1；m=1,2,...,M-1；
将不同阶数(n, m)的二维光场Q(x, y, n, m)拆分为光场Q + (x, y, n, m)和光场Q - (x, y, n, m)后依次加载在空间光调制器上，其中Q + (x, y, n, m)保留Q(x, y, n, m)的正数元素，其余元素用零补齐，Q - (x, y, n, m)保留Q(x, y, n, m)中的负数元素并在负数元素相应位置取其绝对值，其余元素用零补齐，
Figure 228172DEST_PATH_IMAGE012
。
在一实施例中，所述步骤(S3)包括：
利用所述反射回光计算目标物体图像的切比雪夫矩；
Figure 579519DEST_PATH_IMAGE013
其中
Figure 400844DEST_PATH_IMAGE014
为光场Q + (x, y, n, m)照射目标物体后产生的反射回光，
Figure 61502DEST_PATH_IMAGE015
为光场Q - (x, y, n, m)照射目标物体后产生的反射回光。
基于目标物体图像的切比雪夫矩以及所述一系列不同阶数的切比雪夫多项式分布的二维光场，利用Q + (x, y, n, m)、Q - (x, y, n, m)的正交性以及目标物体图像的切比雪夫矩重构出目标物体图像。
重构出的目标物体图像f(x,y)为：
Figure 968278DEST_PATH_IMAGE016
。
本发明的理论依据如下：
在x方向，离散的切比雪夫多项式的表达式为：
Figure 654474DEST_PATH_IMAGE017
(1)
其中x为离散坐标，n为多项式的阶数，N为离散坐标总数。
不同阶数的切比雪夫多项式满足下列正交关系：
Figure 595885DEST_PATH_IMAGE018
(2)
其中归一化参数
Figure 178176DEST_PATH_IMAGE019
满足：
Figure 322981DEST_PATH_IMAGE020
(3)
为使单像素成像系统更为简洁，简化成像算法，故采用归一化的切比雪夫多项式，表达式为：
Figure 750551DEST_PATH_IMAGE021
(4)
对于目标物体图像f(x,y)，归一化后的切比雪夫矩为：
Figure 77627DEST_PATH_IMAGE022
(5)
根据切比雪夫多项式的正交性可重构出目标物体图像f(x,y)，即目标物体图像f(x,y)满足：
Figure 830820DEST_PATH_IMAGE023
(6)
结合单像素成像系统的特性，二维光场Q(x, y, n, m)可设置为:
Figure 446609DEST_PATH_IMAGE024
(7)
由于Q(x, y, n, m)中的元素有正有负，将不同阶数(n, m)的二维光场Q(x, y, n, m)拆分为光场Q + (x, y, n, m)和光场Q - (x, y, n, m)后依次加载在空间光调制器上，其中
Figure 192717DEST_PATH_IMAGE025
，Q + (x, y, n, m)保留Q(x, y, n, m)的正数元素，其余元素用零补齐，Q - (x, y, n, m)保留Q(x, y, n, m)中的负数元素并在负数元素相应位置取其绝对值，其余元素用零补齐。
目标物体图像的切比雪夫矩T nm 可表示为：
Figure 639879DEST_PATH_IMAGE026
(8)
其中，
Figure 298393DEST_PATH_IMAGE027
为光场Q + (x, y, n, m)照射目标物体后产生的反射回光，
Figure 667058DEST_PATH_IMAGE028
为光场Q - (x, y, n, m)照射目标物体后产生的反射回光。
通过依次将不同阶数(n, m)的光场Q + (x, y, n, m)和光场Q - (x, y, n, m)依次照射到目标物体上并收集反射回光
Figure 436430DEST_PATH_IMAGE027
和
Figure 957673DEST_PATH_IMAGE028
，可以重构出目标物体图像f(x,y)，重构目标物体图像f(x,y)满足：
Figure 583826DEST_PATH_IMAGE029
(9)
参照图2，一实施例提供一种高精度单像素成像系统，包括：
空间光调制图案生成模块，用于根据目标物体，生成一系列不同阶数的切比雪夫多项式分布的二维光场并加载到空间光调制器上；
光源3，用于产生激光，并将产生的激光投影照射到空间光调制器上；
空间光调制器2，利用所述一系列不同阶数的切比雪夫多项式分布的二维光场对光源发出的激光进行调制，调制后的激光照射到目标物体4上；
单像素光电探测器5，用于采集目标物体4的反射回光；
成像模块，用于基于所述反射回光以及所述一系列不同阶数的切比雪夫多项式分布的二维光场，对目标物体图像进行重构。
其中所述空间光调制图案生成模块和成像模块的功能均由计算机1实现。
一实施例中，所述空间光调制图案生成模块，包括：
输入模块，确定目标物体，将目标物体离散化大小为M*N像素的矩阵，目标物体的实际大小为
Figure 439787DEST_PATH_IMAGE030
的矩形，其中M、N为正整数，
Figure 12851DEST_PATH_IMAGE031
、
Figure 372288DEST_PATH_IMAGE032
分别为一个像素在y，x方向的几何尺寸；
切比雪夫多项式分布光场生成模块，用于生成一系列不同阶数的切比雪夫多项式分布的二维光场Q(x, y, n, m)，其中x，y为x方向和y方向的像素点坐标，n和m分别为x方向和y方向的切比雪夫多项式的阶数，x=1,2,...,N-1；n=1,2,...,N-1；y=1,2,...,M-1；m=1,2,...,M-1；
处理模块，用于将不同阶数(n, m)的二维光场Q(x, y, n, m)拆分为光场Q + (x, y, n, m)和光场Q - (x, y, n, m)后依次加载在空间光调制器上，其中
Figure 434922DEST_PATH_IMAGE033
，Q + (x, y, n, m)保留Q(x, y, n, m)的正数元素，其余元素用零补齐，Q - (x, y, n, m)保留Q(x, y, n, m)中的负数元素并在负数元素相应位置取其绝对值，其余元素用零补齐。
一实施例中，所述成像模块包括:
切比雪夫矩计算模块，用于利用所述反射回光计算目标物体图像的切比雪夫矩；
目标重构模块，用于基于目标物体图像的切比雪夫矩以及所述一系列不同阶数的切比雪夫多项式分布的二维光场，重构出目标物体图像。
上述各模块功能的实现方法，可以采用前述各实施例中相同的方法实现，在此不再赘述。
本发明中的空间光调制器2的类型不限，可以是数字微镜阵列（DMD）或者空间光相位调制器（SLM），其中数字微镜阵列对二值化图案的刷新速率可达到22 KHz，SLM的刷新速率约为数十Hz。
一实施例采用图2所示结构，采用数字微镜阵列作为空间光调制器，主要用于加载空间光调制图案，对光源发出的光进行强度调制产生结构照明光场。光源3包括激光器以及准直扩束系统，本实施例采用的激光器工作在可见光波段，波长为532 nm，准直扩束系统由准直扩束透镜组成，对激光器发出的光进行扩束、准直，并投影照射到空间光调制器2上。本实施例中采用的目标物体4为一副自然图像（Lena人像图片）。本实施例采用的单像素光电探测器5为光电倍增管结构（Thorlabs PMT-PMM02），用于获取目标物体的反射光强度并转化为电信号输出到计算机1进行目标图像重构。
图3为本实施例采用的不同阶数的切比雪夫多项式分布的二维光场。
图4为本实施例中单像素探测器探测到的不同阶数的切比雪夫矩；对于一个大小为M*N的矩阵，x方向的阶数n取值在0至N-1之间，y方向的阶数m的取值在0至M-1之间，当满采样的情况下，总共有M*N阶二维切比雪夫矩，因此单像素光电探测器5采集到的切比雪夫矩可构成一个大小为M*N的矩阵。由于照明光场将第M*N阶切比雪夫光场Q(x, y, n, m)拆分成两幅光场Q + (x, y, n, m)和Q - (x, y, n, m)，相对应的单像素光电探测器5有两个探测值
Figure 761867DEST_PATH_IMAGE034
和
Figure 404201DEST_PATH_IMAGE035
，那么对应的第m*n阶切比雪夫矩T mn 满足公式(8)。
单像素光电探测器5采集到的
Figure 618144DEST_PATH_IMAGE034
和
Figure 851679DEST_PATH_IMAGE035
输入到计算机1中，计算机1基于公式(10)可重构出目标物体图像f(x,y)。在完整采样条件下，即采样数为M*N，重构出的目标物体图像如图5所示。切比雪夫矩也可以在欠采样的条件下重构出物体图像，当x方向的阶数n依次从0取到0.5*N，y方向的阶数m依次从0取到0.5*M时，此时采样数为0.25*M*N，基于欠采样的条件下切比雪夫光场的单像素成像方法所重构的目标物体如图6所示，此时图像仍然清晰可分辨，只是成像质量相比满采样条件下的目标物体成像结果（如图5所示）有所下降。
以上实施例的各技术特征可以进行任意的组合，为使描述简洁，未对上述实施例中的各个技术特征所有可能的组合都进行描述，然而，只要这些技术特征的组合不存在矛盾，都应当认为是本说明书记载的范围。
以上所述实施例仅表达了本申请的几种实施方式，其描述较为具体和详细，但并不能因此而理解为对发明专利范围的限制。应当指出的是，对于本领域的普通技术人员来说，在不脱离本申请构思的前提下，还可以做出若干变形和改进，这些都属于本申请的保护范围。因此，本申请专利的保护范围应以所附权利要求为准。

Claims (10)

1.一种高精度单像素成像方法，其特征在于，包括：
(S1) 根据目标物体，生成一系列不同阶数的切比雪夫多项式分布的二维光场并加载到空间光调制器上，方法包括：
目标物体离散化大小为M*N像素的矩阵，目标物体的实际大小为
Figure 201672DEST_PATH_IMAGE001
的矩形，其中M、N为正整数，
Figure 587654DEST_PATH_IMAGE002
、
Figure 632970DEST_PATH_IMAGE003
分别为一个像素在y，x方向的几何尺寸；
生成一系列不同阶数的切比雪夫多项式分布的二维光场Q(x, y, n, m)，其中x，y为x方向和y方向的像素点坐标，n和m分别为x方向和y方向的切比雪夫多项式的阶数，x=1,2,...,N-1；n=1,2,...,N-1；y=1,2,...,M-1；m=1,2,...,M-1；
将不同阶数(n, m)的二维光场Q(x, y, n, m)拆分为光场Q+(x, y, n, m)和光场Q-(x, y, n, m)后依次加载在空间光调制器上，其中Q+(x, y, n, m)保留Q(x, y, n, m)的正数元素，其余元素用零补齐，Q-(x, y, n, m)保留Q(x, y, n, m)矩阵二维调制矩阵中的负数元素并在负数元素相应位置取其绝对值，其余元素用零补齐；
(S2) 空间光调制器利用所述一系列不同阶数的切比雪夫多项式分布的二维光场对光源发出的激光进行调制，调制后的激光照射到目标物体上，采集目标物体的反射回光；
(S3) 基于所述反射回光以及所述一系列不同阶数的切比雪夫多项式分布的二维光场，对目标物体图像进行重构。
2.根据权利要求1所述的高精度单像素成像方法，其特征在于，所述步骤(S3)包括：
利用所述反射回光计算目标物体图像的切比雪夫矩；
基于目标物体图像的切比雪夫矩以及所述一系列不同阶数的切比雪夫多项式分布的的二维光场，重构出目标物体图像。
3.根据权利要求2所述的高精度单像素成像方法，其特征在于，所述目标物体图像的切比雪夫矩为：
Figure 247753DEST_PATH_IMAGE004
其中
Figure 530967DEST_PATH_IMAGE005
为光场Q + (x, y, n, m)照射目标物体后产生的反射回光，
Figure 556692DEST_PATH_IMAGE006
为光场Q - (x, y, n, m)照射目标物体后产生的反射回光。
4.根据权利要求3所述的高精度单像素成像方法，其特征在于，重构出的目标物体图像f(x,y)为：
Figure 292567DEST_PATH_IMAGE007
。
5.一种高精度单像素成像系统，其特征在于，包括：
空间光调制图案生成模块，用于根据目标物体，生成一系列不同阶数的切比雪夫多项式分布的二维光场并加载到空间光调制器上；
光源，用于产生激光，并将产生的激光投影照射到空间光调制器上；
空间光调制器，利用所述一系列不同阶数的切比雪夫多项式分布的二维光场对光源发出的激光进行调制，调制后的激光照射到目标物体上；
单像素光电探测器，用于采集目标物体的反射回光；
成像模块，用于基于所述反射回光以及所述一系列不同阶数的切比雪夫多项式分布的二维光场，对目标物体图像进行重构。
6.根据权利要求5所述的高精度单像素成像系统，其特征在于，所述空间光调制图案生成模块和成像模块由计算机实现。
7.根据权利要求5或6所述的高精度单像素成像系统，其特征在于，所述空间光调制图案生成模块，包括：
输入模块，确定目标物体，将目标物体离散化大小为M*N像素的矩阵，目标物体的实际大小为
Figure 757046DEST_PATH_IMAGE008
的矩形，其中M、N为正整数，
Figure 878455DEST_PATH_IMAGE009
、
Figure 340660DEST_PATH_IMAGE010
分别为一个像素在y，x方向的几何尺寸；
切比雪夫多项式分布光场生成模块，用于生成一系列不同阶数的切比雪夫多项式分布的二维光场Q(x, y, n, m)，其中x，y为x方向和y方向的像素点坐标，n和m分别为x方向和y方向的切比雪夫多项式的阶数，x=1,2,...,N-1；n=1,2,...,N-1；y=1,2,...,M-1；m=1,2,...,M-1；
处理模块，用于将不同阶数(n, m)的二维光场Q(x, y, n, m)拆分为光场Q + (x, y, n, m)和光场Q - (x, y, n, m)后依次加载在空间光调制器上，其中Q + (x, y, n, m)保留Q(x, y, n, m)的正数元素，其余元素用零补齐，Q - (x, y, n, m)保留Q(x, y, n, m)中的负数元素并在负数元素相应位置取其绝对值，其余元素用零补齐，
Figure 360569DEST_PATH_IMAGE011
。
8.根据权利要求7所述的高精度单像素成像系统，其特征在于，成像模块包括:
切比雪夫矩计算模块，用于利用所述反射回光计算目标物体图像的切比雪夫矩；
目标重构模块，用于基于目标物体图像的切比雪夫矩以及所述一系列不同阶数的切比雪夫多项式分布的二维光场，重构出目标物体图像。
9.根据权利要求5或6或8所述的高精度单像素成像系统，其特征在于，所述空间光调制器是数字微镜阵列或者空间光相位调制器。
10.根据权利要求5或6或8所述的高精度单像素成像系统，其特征在于，所述光源包括激光器以及准直扩束系统。
CN202210641951.1A 2022-06-08 2022-06-08 高精度单像素成像方法及系统 Active CN114721006B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210641951.1A CN114721006B (zh) 2022-06-08 2022-06-08 高精度单像素成像方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210641951.1A CN114721006B (zh) 2022-06-08 2022-06-08 高精度单像素成像方法及系统

Publications (2)

Publication Number Publication Date
CN114721006A true CN114721006A (zh) 2022-07-08
CN114721006B CN114721006B (zh) 2022-09-02

Family

ID=82233042

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210641951.1A Active CN114721006B (zh) 2022-06-08 2022-06-08 高精度单像素成像方法及系统

Country Status (1)

Country Link
CN (1) CN114721006B (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116609794A (zh) * 2023-07-21 2023-08-18 中国人民解放军国防科技大学 基于径向切比雪夫光场的单像素成像方法、装置及设备
CN116893429A (zh) * 2023-07-13 2023-10-17 中国人民解放军国防科技大学 基于圆谐傅里叶光场的单像素成像方法和目标识别方法

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20040093364A1 (en) * 2002-08-01 2004-05-13 Cheng Andrew F. Data compression using Chebyshev transform
US20110176709A1 (en) * 2010-01-21 2011-07-21 Samsung Electronics Co., Ltd. Method and apparatus for calculating a distance between an optical apparatus and an object
CN106371201A (zh) * 2016-11-03 2017-02-01 清华大学 基于计算鬼成像的傅里叶重叠关联成像系统及方法
CN109343078A (zh) * 2018-11-27 2019-02-15 北京理工大学 一种具有目标导向性切片三维鬼成像实现方法
CN110109083A (zh) * 2019-05-05 2019-08-09 上海交通大学 基于一维光相控阵的三维扫描激光雷达
CN110703276A (zh) * 2019-08-30 2020-01-17 清华大学深圳研究生院 在强散射条件下的傅里叶成像装置及方法
CN111094940A (zh) * 2017-07-19 2020-05-01 小川宏 断层图像拍摄装置
CN111650600A (zh) * 2020-05-28 2020-09-11 北京空间机电研究所 一种面向极微弱信号的双谱段激光成像装置
CN111833265A (zh) * 2020-06-15 2020-10-27 北京邮电大学 一种基于组稀疏循环调制的鬼成像图像恢复方案

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20040093364A1 (en) * 2002-08-01 2004-05-13 Cheng Andrew F. Data compression using Chebyshev transform
US20110176709A1 (en) * 2010-01-21 2011-07-21 Samsung Electronics Co., Ltd. Method and apparatus for calculating a distance between an optical apparatus and an object
CN106371201A (zh) * 2016-11-03 2017-02-01 清华大学 基于计算鬼成像的傅里叶重叠关联成像系统及方法
CN111094940A (zh) * 2017-07-19 2020-05-01 小川宏 断层图像拍摄装置
CN109343078A (zh) * 2018-11-27 2019-02-15 北京理工大学 一种具有目标导向性切片三维鬼成像实现方法
CN110109083A (zh) * 2019-05-05 2019-08-09 上海交通大学 基于一维光相控阵的三维扫描激光雷达
CN110703276A (zh) * 2019-08-30 2020-01-17 清华大学深圳研究生院 在强散射条件下的傅里叶成像装置及方法
CN111650600A (zh) * 2020-05-28 2020-09-11 北京空间机电研究所 一种面向极微弱信号的双谱段激光成像装置
CN111833265A (zh) * 2020-06-15 2020-10-27 北京邮电大学 一种基于组稀疏循环调制的鬼成像图像恢复方案

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
MEER NAZMUS SAKIB,: "Impact of Optical Transmission on Multiband OFDM Ultra-Wideband Wireless System With Fiber Distribution", 《JOURNAL OF LIGHTWAVE TECHNOLOGY》 *
左超 等: "计算光学成像：何来，何处，何去，何从？", 《红外与激光工程》 *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116893429A (zh) * 2023-07-13 2023-10-17 中国人民解放军国防科技大学 基于圆谐傅里叶光场的单像素成像方法和目标识别方法
CN116893429B (zh) * 2023-07-13 2024-01-30 中国人民解放军国防科技大学 基于圆谐傅里叶光场的单像素成像方法和目标识别方法
CN116609794A (zh) * 2023-07-21 2023-08-18 中国人民解放军国防科技大学 基于径向切比雪夫光场的单像素成像方法、装置及设备
CN116609794B (zh) * 2023-07-21 2023-09-26 中国人民解放军国防科技大学 基于径向切比雪夫光场的单像素成像方法、装置及设备

Also Published As

Publication number Publication date
CN114721006B (zh) 2022-09-02

Similar Documents

Publication Publication Date Title
CN114721006B (zh) 高精度单像素成像方法及系统
CN110425986B (zh) 基于单像素传感器的三维计算成像方法及装置
JP6415685B2 (ja) ランダム格子に基づいた圧縮センシング広帯域ハイパースペクトルイメージングシステム
US9961285B2 (en) Optical imaging method using single pixel detector
CN107306333B (zh) 一种高速单像素成像方法
Sun et al. Differential computational ghost imaging
JP2013535931A (ja) 圧縮撮像装置の画像取得時間の減少
CN110471082B (zh) 单像素激光计算成像装置及方法
US20150003750A1 (en) Reconstructing an image of a scene captured using a compressed sensing device
CN103913229B (zh) 编码模板多目标主动成像光谱系统及方法
CN114095718B (zh) 单像素成像系统及方法
Yu et al. Single-photon compressive imaging with some performance benefits over raster scanning
CN116893429B (zh) 基于圆谐傅里叶光场的单像素成像方法和目标识别方法
CN108007385B (zh) 一种大视场弹坑表面形貌成像系统及方法
CN111443358B (zh) 单点探测成像系统
Kuusela Single-pixel camera
CN111781733A (zh) 基于光波调制和相位恢复的多层复数域成像方法和装置
Suo et al. A self-synchronized high speed computational ghost imaging system: A leap towards dynamic capturing
CN111854945B (zh) 一种单像素紫外偏振成像方法及系统
Thapa et al. Less is more: compressive sensing in optics and image science
CN115272098A (zh) 一种基于数字微镜高斯矩阵二值表示的单像素成像方法
CN111640063B (zh) 基于空间频域多尺度调制与重建的压缩成像系统及方法
CN114710615B (zh) 高效单像素成像方法及系统
JP2019092088A (ja) 撮像装置
CN112866532B (zh) 一种紧凑型双光路单像素成像系统及不均匀光源校正方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant