CN108564090B - 一种基于信号加权的快速傅里叶单像素成像方法 - Google Patents

一种基于信号加权的快速傅里叶单像素成像方法 Download PDF

Info

Publication number
CN108564090B
CN108564090B CN201810274897.5A CN201810274897A CN108564090B CN 108564090 B CN108564090 B CN 108564090B CN 201810274897 A CN201810274897 A CN 201810274897A CN 108564090 B CN108564090 B CN 108564090B
Authority
CN
China
Prior art keywords
pattern
gray
fourier
patterns
pixel
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.)
Active
Application number
CN201810274897.5A
Other languages
English (en)
Other versions
CN108564090A (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.)
Hefei Institutes of Physical Science of CAS
Original Assignee
Hefei Institutes of Physical Science of CAS
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 Hefei Institutes of Physical Science of CAS filed Critical Hefei Institutes of Physical Science of CAS
Priority to CN201810274897.5A priority Critical patent/CN108564090B/zh
Publication of CN108564090A publication Critical patent/CN108564090A/zh
Application granted granted Critical
Publication of CN108564090B publication Critical patent/CN108564090B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V10/00—Arrangements for image or video recognition or understanding
    • G06V10/20—Image preprocessing
    • G06V10/28—Quantising the image, e.g. histogram thresholding for discrimination between background and foreground patterns
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/14—Fourier, Walsh or analogous domain transformations, e.g. Laplace, Hilbert, Karhunen-Loeve, transforms
    • G06F17/141—Discrete Fourier transforms
    • G06F17/142—Fast Fourier transforms, e.g. using a Cooley-Tukey type algorithm

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Mathematical Physics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Mathematical Analysis (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • Discrete Mathematics (AREA)
  • Multimedia (AREA)
  • Algebra (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Image Analysis (AREA)
  • Image Processing (AREA)

Abstract

本发明公开了一种基于信号加权的快速傅里叶单像素成像方法，通过将灰度傅里叶基图案分解为一系列相应的二值化图案，再根据这些二值化图案利用数字微镜器件对光源进行调制，单像素探测器依次同步测量调制光与目标物体产生相互作用信号；通过对单像素探测器测量的信号进行加权求和来替代相应灰度傅里叶图案与物体相互作用的强度信号，以此获得相应的傅里叶变换谱系数，最后对所有或者部分谱系数构成的变换谱进行逆傅里叶变换来获得目标物体的图像。本发明大幅减少投影时间，可以在不牺牲成像空间分辨率的前提下提高了傅里叶单像素成像效率。

Description

一种基于信号加权的快速傅里叶单像素成像方法
技术领域
本发明涉及光学成像技术领域，特别涉及一种基于信号加权的快速傅里叶单像素成像方法。
背景技术
单像素成像技术是一种在成像机理上和传统成像技术有着本质区别的技术，因为其具有一些独特的优点，受到人们的关注。单像素成像技术使用没有空间分辨能力的探测器(如光电倍增管)进行时间上的多次采样，实现对目标物体成像。单像素成像技术最早源于利用量子纠缠效应的鬼成像技术，后来发展出利用热光的单像素鬼成像技术。由于单像素成像技术使用仅有一个像素的探测器来进行物理信息的采样，需要使用大量在空间上进行调制的不同光场，按照时间先后顺序投射到目标物体，以获取空间分辨。是一种以牺牲时间分辨率来换取空间分辨率的成像方案。许多科研人员提出不同的方法来提高成像效率和成像质量。[Sun M J,Meng L T,Edgar M P,et al.A Russian Dolls ordering of theHadamard basis for compressive single-pixel imaging[J].Sci Rep,2017,7(1):3464]提出使用Hadamard基图案来进行单像素成像，提高成像效率。[Zhang Z,Ma X,ZhongJ.Single-pixel imaging by means of Fourier spectrum acquisition[J].NatureCommunications,2015,6:6225]提出采用傅立叶基图案来提高成像质量。采用Hadamard基图案和傅立叶基图案的成像系统比采用随机图案的成像系统无论在成像效率还是成像质量上都有巨大改进。[Zhang Z,Wang X,Zheng G,et al.Hadamard single-pixel imagingversus Fourier single-pixel imaging.[J].Optics Express,2017,25(16):19619-19639]对采用Hadamard基图案和傅里叶基图案进行了比较，指出傅立叶基图案效率更高。在单像素成像系统中数字微镜器件是广泛使用的高速空间光调制器件，目前数字微镜器件对光场的二值化调制频率可到20KHz以上，但是调制8位灰度级的图案时，频率大约为250Hz，主要原因是将灰度图案分解为8幅二值化的图案，然后将这8幅二值化图案按照时间权重进行投影，此方法增加了大量的投影时间。傅里叶单像素成像系统中采用灰度余弦空间结构光场，在使用数字微镜器件进行实际应用时，投影频率较低，严重制约了傅里叶单像素成像系统的效率。[Zhang Z,Wang X,Zheng G,et al.Fast Fourier single-pixelimaging via binary illumination[J].Scientific Reports,2017,7(1):12029]提出通过对灰度级的傅里叶基图案进行插值一定的倍数后再进行二值化，通过数字微镜器件实现高速成像，但是此方法牺牲了傅里叶单像素成像的空间分辨率。本专利提出一种在不牺牲成像空间分辨率情况下，有效提高傅里叶单像素成像系统效率的方法。
发明内容
本发明的目的是为了解决现有技术的不足，针对使用高速数字微镜器件的傅里叶单像素成像系统，提出一种基于信号加权的快速傅里叶单像素成像方法，来提高成像效率，其技术核心采用具有确定的数学函数解析表达的灰度傅里叶空间结构光场，将灰度调制图案分解为二值化图案，二值化图案按照相同的时间照明目标物体，将对应探测器端的光强信息进行权重相加，进而获得目标物体的傅立叶谱系数，最后应用傅立叶逆变换来重构物体的像。本发明可以在不牺牲成像空间分辨率情况下，降低成像时间，提高成像效率。
本发明的技术方案如下：
一种基于信号加权的快速傅里叶单像素成像方法，包括以下步骤：
(1)计算机产生一系列灰度傅里叶基图案；
(2)将灰度傅里叶基图案按照像素值的正负分解为相应的灰度图案对；然后将分解后的图案对按照对应的量化权重依次分解成一系列二值化图案，把这些二值化图案按顺序加载到数字微镜器件上，数字微镜器件对光源进行调制后与目标物体产生相互作用信号，单像素探测器同步依次测量这些相互作用信号；
(3)将单像素探测器测量到的强度信号按照灰度图案分解成二值化图案的权重系数相加，作为对应灰度傅里叶基图案与物体相互作用信号，从而得到相应频率下目标物体的傅里叶变换谱系数；
(4)对所有频率或者部分频率的谱系数构成的傅里叶变换谱进行逆傅里叶变换，获得目标物体的图像。
进一步地，在于步骤(1)中：灰度傅里叶基图案按照公式Pφ(x,y；u,v)＝[(2R-1)*cos(2πxu+2πyv+φ)]生成，傅里叶基图案大小为N x M像素，[·]为取整操作，(x,y)为二维空间坐标，(u,v)为x,y方向上的频率，参数R为灰度图案的量化位数，φ＝0或者φ＝π/2。
进一步地，将大小为N x M图案Pφ(x,y；u,v)按照空间坐标(x,y)像素值的正负分解为相应的灰度图案对
Figure BDA0001613418380000031
分解的原则是当Pφ(x,y；u,v)空间坐标(x,y)像素值为正时，分解产生的图案像素值不变，否则图案像素值为0，产生大小为N x M的图案
Figure BDA0001613418380000032
下一步，当空间坐标(x,y)像素值为负时，图案像素值为原图案像素值的绝对值，否则图案像素值为0，产生大小为N x M的图案
Figure BDA0001613418380000033
进一步地，将大小为N x M的灰度图案
Figure BDA0001613418380000034
和
Figure BDA0001613418380000035
按十进制转二进制的原则分别分解成R张大小为N x M的二值化图案。
进一步地，灰度图案
Figure BDA0001613418380000036
和
Figure BDA0001613418380000037
分别产生的R张大小为N x M的二值化图案调制光源后与目标物体相互作用，作用信号依次被单像素探测器接收，将探测信号按照灰度图案分解的次序分别按照2的a次方的权重进行相加(对于量化位数为R时，a的变化范围是从0到R-1)，作为灰度图案
Figure BDA0001613418380000038
和
Figure BDA0001613418380000039
投影照明时的等效强度值，将这两个等效的强度值做差，作为该频率下物体的傅里叶谱系数的实部。
进一步地，将灰度图案
Figure BDA00016134183800000310
和
Figure BDA00016134183800000311
分别产生的R张二值化图案调制光源后与目标物体相互作用，作用信号依次被单像素探测器接收，将探测信号按照灰度图案分解的次序分别按照2的a次方的权重进行相加(对于量化位数为R时，a的变化范围是从0到R-1)，作为灰度图案
Figure BDA00016134183800000312
和
Figure BDA00016134183800000313
投影照明时的等效强度值，将这两个等效的强度做差，作为该频率下物体的傅里叶谱系数的虚部。
本发明的理论依据如下：
当采用傅里叶基图案进行单像素成像时，频率为(un,vm)的傅里叶调制光与物体O(x,y)相互作用，单像素探测器接收的光强为Iφ(un,vm)，其可以表示为下式：
Figure BDA0001613418380000041
以6位灰度级为例，即R＝6，使用数字微镜器件调制灰度傅里叶基图案Pφ(x,y；un,vm)的效果可以用6幅二值化图案的时间权重投影来等效，即：
Figure BDA0001613418380000042
其中T为投影单位时间，Bφ,i(x,y；un,vm)为二值化图案。
将(2)式带入(1)式可得：
Figure BDA0001613418380000043
其中，Iφ,1(un,vm)、…、Iφ,6(un,vm)分别为6幅二值化图案Bφ,1(x,y；un,vm)、…、Bφ,6(x,y；un,vm)与物体相互作用的信号。由于二值化散斑Bφ,1(x,y；un,vm)、…、Bφ,6(x,y；un,vm)照明投影的时间T是相同的，单像素探测器依次接收二值化图案调制光与物体的相互作用信号，并将相应的强度信号按照2的a次方的权重进行相加(对于量化位数为6时，a的变化范围是从0到5)，即可等效为6位灰度图案的探测强度值。在投影频率为20kHz时，本发明等效的一张6位灰度傅里叶图案的投影时间为：50us x 6＝300us，即等效的投影频率高达3.3kHz，而传统的方法投影6位灰度图案的最高频率～300Hz。
本发明的有益效果是：
本发明大幅减少投影时间，可以在不牺牲成像空间分辨率的前提下提高了傅里叶单像素成像效率。
附图说明
图1为频率为(1/128,1/128)时产生的灰度傅里叶基图案。
图2为图1灰度傅里叶基图案分解产生的图案对。
图3为图2分解产生的二值化图案。
图4为装置实施示意图。
图5为部分傅里叶谱及其复原的结果。
图6为全部傅里叶谱及其复原的结果。
具体实施方式
图4为实施装置示意图，以6位灰度级的傅里叶基图案，即R＝6，傅里叶基图案以及复原图像的像素为128x 128为例进行实施说明。
一种基于信号加权的快速傅里叶单像素成像方法，包括以下步骤：
由计算机1按照公式Pφ＝0(x,y；u,v)＝[(2R-1)*cos(2πxu+2πyv+φ)]和Pφ＝π/2(x,y；u,v)＝[(2R-1)*cos(2πxu+2πyv+φ)]分别生成一系列不同频率(u,v)的128x 128像素的傅里叶基图案，其中(x,y)为128x 128的二维空间坐标。图1为频率为(1/128,1/128)时生成的灰度傅里叶基图案。
将大小为128x 128图案Pφ＝0(x,y；u,v)和Pφ＝π/2(x,y；u,v)分别按照空间坐标(x,y)像素值的正负分解为相应的灰度级图案对，即当Pφ＝0(x,y；u,v)和Pφ＝π/2(x,y；u,v)空间坐标(x,y)像素值为正时，图案像素保持不变，否则像素为0，分别生成128x 128的灰度图案
Figure BDA0001613418380000051
和
Figure BDA0001613418380000052
当Pφ＝0(x,y；u,v)和Pφ＝π/2(x,y；u,v)像素值为负时，图案像素为原图案像素的绝对值，否则像素为0，分别生成128x 128的灰度图案
Figure BDA0001613418380000053
和
Figure BDA0001613418380000054
图2为频率(1/128,1/128)灰度傅里叶基图案分别分解成的图案对
Figure BDA0001613418380000055
与
Figure BDA0001613418380000056
以及
Figure BDA0001613418380000057
与
Figure BDA0001613418380000058
将大小为128x 128的
Figure BDA0001613418380000059
Figure BDA00016134183800000510
分别分解(原理为十进制转二进制)为6张大小为128x 128的二值化图案，图3为(1/128,1/128)时分解后的二值化图案。
将分解后的二值化图案加载到数字微镜器件3上(TI Discovery V7100)，对光源2(LED白光源)进行调制经透镜4投影到目标物体5上。单像素探测器6(Thorlabs PMT-PMM02)依次接收二值化图案与物体相互作用的强度信号，数据采集卡7(NI USB-6211)实现对强度信号的数字化。
下面以un＝1/128,vm＝1/128为例进行说明。
Figure BDA0001613418380000061
分解的6张二值化图案照明物体5时，将采集卡7数字化的强度信号按照公式
Figure BDA0001613418380000062
进行权重求和，其中
Figure BDA0001613418380000063
分别为对应二值化图案
Figure BDA0001613418380000064
与物体相互作用强度信号值；
Figure BDA0001613418380000065
分解的6张二值化图案照明物体5时，将采集卡7数字化的强度信号按照公式
Figure BDA0001613418380000066
进行权重求和，其中
Figure BDA0001613418380000067
分别为对应二值化图案
Figure BDA0001613418380000068
与物体相互作用的强度信号值。因此可得该频率下目标物体傅里叶谱系数的实部
Figure BDA0001613418380000069
Figure BDA00016134183800000610
分解的6张二值化图案照明物体时，将采集卡7数字化的强度信号按照公式
Figure BDA00016134183800000611
进行权重求和，其中
Figure BDA00016134183800000612
分别为对应二值化图案
Figure BDA00016134183800000613
与物体相互作用的强度信号值；
Figure BDA00016134183800000614
分解的6张二值化图案照明物体时，将采集卡7数字化的强度信号按照公式
Figure BDA00016134183800000615
进行权重求和，其中
Figure BDA00016134183800000616
分别为对应二值化图案
Figure BDA00016134183800000617
与物体相互作用强度信号值。因此可得该频率下目标物
体傅里叶谱系数的虚部
Figure BDA00016134183800000618
从而得到该频率下目标物体傅里叶谱系数
D(1/128,1/128)＝Iφ＝0(1/128,1/128)-j·Iφ＝π/2(1/128,1/128)。
依次改变频率un(变化范围0/128～127/128)、vm(变化范围0/128～127/128)，获得所有频率下傅里叶谱系数，最后利用全部或者部分傅里叶谱系数，进行逆傅里叶变换即可实现对目标物体成像。将获得的傅里叶谱系数按照频率从低到高进行排列，图5为低25％空间频率的谱分布及其逆傅里叶变换复原的结果。图6使用全部空间频率谱分布及其逆变换复原的结果。

Claims (1)

1.一种基于信号加权的快速傅里叶单像素成像方法，其特征在于包括以下步骤：
(1)计算机产生一系列灰度傅里叶基图案；
(2)将灰度傅里叶基图案按照像素值的正负分解为相应的灰度图案对；然后将分解后的灰度图案对按照对应的量化权重依次分解成一系列二值化图案，把这些二值化图案按顺序加载到数字微镜器件上，数字微镜器件对光源进行调制后与目标物体产生相互作用信号，单像素探测器同步依次测量上述相互作用信号；
(3)将单像素探测器测量到的强度信号按照灰度图案分解成二值化图案的权重系数相加，作为对应灰度傅里叶基图案与物体相互作用信号，从而得到相应频率下目标物体的傅里叶变换谱系数；
(4)对所有频率或者部分频率的谱系数构成的傅里叶变换谱进行逆傅里叶变换，获得目标物体的图像；
所述步骤(1)中：灰度傅里叶基图案按照公式Pφ(x,y；u,v)＝[(2R-1)*cos(2πxu+2πyv+φ)]生成，傅里叶基图案大小为N x M像素，[·]为取整操作，(x,y)为二维空间坐标，(u,v)为x,y方向上的频率，参数R为灰度图案的量化位数，φ＝0或者φ＝π/2；
将所述大小为N x M图案Pφ(x,y；u,v)按照空间坐标(x,y)像素值的正负分解为相应的灰度图案对
Figure FDA0003067081310000011
分解的原则是当Pφ(x,y；u,v)空间坐标(x,y)像素值为正时，分解产生的图案像素值不变，否则图案像素值为0，产生大小为N x M的图案
Figure FDA0003067081310000012
当空间坐标(x,y)像素值为负时，图案像素值为原图案像素值的绝对值，否则图案像素值为0，产生大小为N x M的图案
Figure FDA0003067081310000013
将所述大小为N x M的灰度图案
Figure FDA0003067081310000014
和
Figure FDA0003067081310000015
按十进制转二进制的原则分别分解成R张大小为N x M的二值化图案；
将所述灰度图案
Figure FDA0003067081310000016
和
Figure FDA0003067081310000017
分别产生的R张大小为N x M的二值化图案调制光源后与目标物体相互作用，作用信号依次被单像素探测器接收，将探测信号按照灰度图案分解的次序分别按照2的a次方的权重进行相加，对于量化位数为R时，a的变化范围是从0到R-1，作为灰度图案
Figure FDA0003067081310000021
和
Figure FDA0003067081310000022
投影照明时的等效强度值，将这两个等效的强度值做差，作为该频率下物体的傅里叶谱系数的实部；
将所述灰度图案
Figure FDA0003067081310000023
和
Figure FDA0003067081310000024
分别产生的R张二值化图案调制光源后与目标物体相互作用，作用信号依次被单像素探测器接收，将探测信号按照灰度图案分解的次序分别按照2的a次方的权重进行相加，对于量化位数为R时，a的变化范围是从0到R-1，作为灰度图案
Figure FDA0003067081310000025
和
Figure FDA0003067081310000026
投影照明时的等效强度值，将这两个等效的强度做差，作为该频率下物体的傅里叶谱系数的虚部。
CN201810274897.5A 2018-03-30 2018-03-30 一种基于信号加权的快速傅里叶单像素成像方法 Active CN108564090B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810274897.5A CN108564090B (zh) 2018-03-30 2018-03-30 一种基于信号加权的快速傅里叶单像素成像方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810274897.5A CN108564090B (zh) 2018-03-30 2018-03-30 一种基于信号加权的快速傅里叶单像素成像方法

Publications (2)

Publication Number Publication Date
CN108564090A CN108564090A (zh) 2018-09-21
CN108564090B true CN108564090B (zh) 2021-07-09

Family

ID=63533615

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810274897.5A Active CN108564090B (zh) 2018-03-30 2018-03-30 一种基于信号加权的快速傅里叶单像素成像方法

Country Status (1)

Country Link
CN (1) CN108564090B (zh)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110425986B (zh) * 2019-07-17 2020-10-16 北京理工大学 基于单像素传感器的三维计算成像方法及装置
CN110427951B (zh) * 2019-07-18 2021-09-28 北京理工大学 一种基于非均匀采样策略的傅里叶单像素成像方法
CN110809102B (zh) * 2019-10-11 2020-10-30 北京理工大学 基于二值调制的成像加速方法及装置
CN111707413B (zh) * 2020-05-15 2021-11-26 中国科学院合肥物质科学研究院 一种基于单像素探测器的质心探测方法
CN112153254B (zh) * 2020-08-31 2022-02-25 合肥工业大学 一种基于基图的两步相移单像素成像方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7945098B2 (en) * 2007-05-17 2011-05-17 The United States Of America As Represented By The Secretary Of Agriculture Method for characterizing the density and cross-section morphology of trees
CN103024900A (zh) * 2013-01-07 2013-04-03 北京理工大学 无线层析成像网络的分数阶傅立叶多径伪像修正定位方法
CN105785327A (zh) * 2016-01-19 2016-07-20 西安电子科技大学 频率分集阵列合成孔径雷达高分辨宽测绘带成像方法
CN106772310A (zh) * 2016-12-30 2017-05-31 南京理工大学 一种基于计算鬼成像的目标鉴定装置及其方法
CN106842880A (zh) * 2017-03-27 2017-06-13 深圳市美誉镜界光电科技有限公司 全息图像生成方法、处理器及全息图像显示装置、设备

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7945098B2 (en) * 2007-05-17 2011-05-17 The United States Of America As Represented By The Secretary Of Agriculture Method for characterizing the density and cross-section morphology of trees
CN103024900A (zh) * 2013-01-07 2013-04-03 北京理工大学 无线层析成像网络的分数阶傅立叶多径伪像修正定位方法
CN105785327A (zh) * 2016-01-19 2016-07-20 西安电子科技大学 频率分集阵列合成孔径雷达高分辨宽测绘带成像方法
CN106772310A (zh) * 2016-12-30 2017-05-31 南京理工大学 一种基于计算鬼成像的目标鉴定装置及其方法
CN106842880A (zh) * 2017-03-27 2017-06-13 深圳市美誉镜界光电科技有限公司 全息图像生成方法、处理器及全息图像显示装置、设备

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
Accurate and fast 3D surface measurement with temporal-spatial binary encoding structured illumination;JIANGPING ZHU 等;《OPTICS EXPRESS》;20161212;第24卷(第25期);1-12 *
Enhancing resolution of single-pixel imaging system;Dongfeng Shi 等;《Optical review》;20150819;第22卷(第5期);802-808 *
Fast Fourier single-pixel imaging via binary illumination;Zibang Zhang 等;《SCIENTIFIC REPORTS》;20170920;1-9 *
基于压缩感知的单像素视频采样关键技术研究及实现;曾潇;《中国优秀硕士学位论文全文数据库 信息科技辑》;20151215(第12期);I138-877 *
磁共振图像处理中部分傅里叶重建算法的比较;何汶静 等;《重庆医学》;20160720;第45卷(第20期);2804-2809 *

Also Published As

Publication number Publication date
CN108564090A (zh) 2018-09-21

Similar Documents

Publication Publication Date Title
CN108564090B (zh) 一种基于信号加权的快速傅里叶单像素成像方法
CN107306333B (zh) 一种高速单像素成像方法
CN108895985B (zh) 一种基于单像素探测器的物体定位方法
Yu et al. Compressive microscopic imaging with “positive–negative” light modulation
Chen et al. Discrete cosine single-pixel microscopic compressive imaging via fast binary modulation
Xiao et al. Direct single-step measurement of Hadamard spectrum using single-pixel optical detection
Nie et al. Noise-robust computational ghost imaging with pink noise speckle patterns
CN113709325B (zh) 基于Hadamard频域变换矩阵阈值滤波的单像素成像方法
CN113114882B (zh) 一种高采样效率的傅里叶单像素成像方法
Guo et al. Compressed ghost edge imaging
Liang et al. Fast Fourier single-pixel imaging based on Sierra–Lite dithering algorithm
Li et al. Fourier transform profilometry using single-pixel detection based on two-dimensional discrete cosine transform
Yuan et al. Adaptive and dynamic ordering of illumination patterns with an image dictionary in single-pixel imaging
Saqueb et al. Compressive terahertz imaging using a single-bit sensor
Suo et al. A self-synchronized high speed computational ghost imaging system: A leap towards dynamic capturing
Wu et al. Development of a DMD-based compressive sampling hyperspectral imaging (CS-HSI) system
CN112153254B (zh) 一种基于基图的两步相移单像素成像方法
CN109804229B (zh) 电磁波相位振幅生成装置、电磁波相位振幅生成方法以及存储有电磁波相位振幅生成程序的非临时记录介质
CN111854945B (zh) 一种单像素紫外偏振成像方法及系统
Xue et al. Improving the signal-to-noise ratio of complementary compressive imaging with a threshold
CN111640063A (zh) 基于空间频域多尺度调制与重建的压缩成像系统及方法
US11029149B2 (en) Multipath mitigation for time of flight system
Sun et al. High-quality reconstruction of single-pixel imaging using discrete W transform
KR102258774B1 (ko) 적외선 영상 측정 장치 및 방법
Liu et al. Compressive imaging based on multi-scale modulation and reconstruction in spatial frequency domain

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