WO2020181560A1 - Procédé et système d'imagerie basés sur une fluctuation d'intensité de lumière aléatoire - Google Patents

Procédé et système d'imagerie basés sur une fluctuation d'intensité de lumière aléatoire Download PDF

Info

Publication number
WO2020181560A1
WO2020181560A1 PCT/CN2019/078183 CN2019078183W WO2020181560A1 WO 2020181560 A1 WO2020181560 A1 WO 2020181560A1 CN 2019078183 W CN2019078183 W CN 2019078183W WO 2020181560 A1 WO2020181560 A1 WO 2020181560A1
Authority
WO
WIPO (PCT)
Prior art keywords
imaging target
distribution function
light intensity
expression
imaging
Prior art date
Application number
PCT/CN2019/078183
Other languages
English (en)
Chinese (zh)
Inventor
王宇
张罗莎
王魁波
朱精果
杨光华
赵复生
韩哲
亓岩
颜博霞
韩春蕊
郭馨
陈进新
崔惠绒
罗艳
谢婉露
周翊
吴晓斌
Original Assignee
中国科学院微电子研究所
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 中国科学院微电子研究所 filed Critical 中国科学院微电子研究所
Priority to PCT/CN2019/078183 priority Critical patent/WO2020181560A1/fr
Publication of WO2020181560A1 publication Critical patent/WO2020181560A1/fr

Links

Images

Classifications

    • G—PHYSICS
    • G02—OPTICS
    • G02B—OPTICAL ELEMENTS, SYSTEMS OR APPARATUS
    • G02B27/00—Optical systems or apparatus not provided for by any of the groups G02B1/00 - G02B26/00, G02B30/00

Definitions

  • the present disclosure relates to an imaging method and system based on random light intensity fluctuations.
  • correlative imaging The imaging obtained by correlating the fluctuations of the light field intensity is called correlative imaging. Correlated imaging has become a research hotspot at home and abroad because of its theoretical advantage in breaking the diffraction limit of classical optical systems.
  • the traditional dual-arm intensity-associated imaging scheme uses the laser 1a shown in Figure 1 to illuminate the rotating ground glass 2a to generate pseudothermal light with random fluctuations in light intensity.
  • the signal light illuminates the imaging target 5a and carries the imaging target.
  • the signal light of the information is received by the single-pixel detector 6a, and the reference light is received by the area array detector 4a after being split.
  • the correlation result of the single-pixel detector 6a and the area array detector 4a is used to reconstruct the image of the imaging target.
  • the image imaging and image reconstruction speed of this system is slow. Since the working principle of the area array detector 4a used in the solution determines that each frame of image requires sufficient integration time and readout time, and subsequent analog processing circuits and data acquisition circuits will take more time.
  • the sampling frequency of area array detectors is mostly at the MHz level, and the more pixels, the longer the time required for image acquisition per frame.
  • the resolution of the imaging system is limited by the size of the pixel unit of the area array detector 4a.
  • the computational ghost imaging scheme developed on the basis of the dual-arm intensity-associated imaging scheme uses a laser 1b to illuminate the microlens array or the projection device 7b to produce pseudothermographic illumination imaging targets with known light intensity fluctuations 5b, the signal light carrying the imaging target information is received by the single-pixel detector 6b, and the total light intensity information received by the single-pixel detector 6b is correlated with the known pseudo-thermal light field information to reconstruct the image of the imaging target.
  • the system needs to introduce complex optical modules such as digital microlens arrays or projection systems to generate randomly fluctuating light field distributions.
  • the pixel unit of the digital microlens array is about 10 ⁇ m, and the resolution of the light field is relatively low after transmission. Therefore, this solution can only be used for remote sensing, imaging of buildings or daily macroscopic objects.
  • the current image reconstruction of the associated imaging method is based on the correlation between the intensity of the reference light and the signal light.
  • Various improvements are based on the reference light signal, that is, the acquisition of the random fluctuation light field.
  • Information needs to be added to a cumbersome optical module, and the resolution of the system is also easily affected by the limitation of the pixel unit size of the area array detector or the transmission of the digital microlens array, making the resolution lower.
  • An aspect of the present disclosure provides an imaging method based on random light intensity fluctuations.
  • the method includes: step S1, a light source irradiates ground glass and forms a random fluctuation light field on the surface of the ground glass, the random fluctuation light field After illuminating the imaging target, it is collected by a single-pixel detector to obtain the total light intensity; step S2, rotating the ground glass N times, and repeating step S1 to obtain N total light intensity values, where N ⁇ 1; step S3, establishing the imaging target According to the distribution function expression of the random fluctuation light field, the N total light intensity values and the preset relationship between the distribution function expression, determine the distribution function and position information of the imaging target; step S4, reconstruct an image of the imaging target according to the distribution function and position information of the imaging target.
  • the distribution function expression of the imaging target includes an exponential function, Fourier function, trigonometric function, polynomial, Gaussian function, or Weibull function, and the distribution function expression is an M-dimensional K-order function, where M ⁇ 1; K ⁇ 1.
  • the distribution function expression of the imaging target is a two-dimensional first-order Gaussian function
  • x 0 and y 0 are position information of the imaging target; a, b, and c are mathematical parameters of a two-dimensional first-order Gaussian function.
  • the preset relationship is expressed by the following expression:
  • f(x, y) is the random fluctuation light field distribution
  • t(x, y) is the distribution function of the imaging target
  • I i is the total light intensity
  • i 0 ⁇ N.
  • the determining the distribution function and position information of the imaging target according to the preset relationship between the random fluctuation light field, the N total light intensity values and the distribution function expression includes :
  • the values of mathematical parameters a, b, and c in the distribution function of the imaging target and the values of position information x 0 , y 0 are determined.
  • the imaging target size is micrometers and below.
  • the system includes a light source, ground glass, an imaging target, and a single-pixel detector; wherein the light source illuminates the ground glass and forms a random pattern on the ground glass surface.
  • a fluctuating light field the random fluctuating light field illuminates the imaging target and is collected by a single-pixel detector to obtain the total light intensity, which is used to reconstruct the image of the imaging target.
  • the imaging system further includes a data processing module configured to reconstruct an image of the imaging target according to the value of the total light intensity.
  • the data processing module is configured to obtain N total light intensity values by rotating the ground glass N times, where N ⁇ 1; establishing the distribution function expression of the imaging target, and according to the random fluctuation light
  • the predetermined relationship between the field, the N total light intensity values and the distribution function expression determines the distribution function and position information of the imaging target; and reconstructs imaging according to the distribution function and position information of the imaging target The image of the target.
  • the distribution function expression of the imaging target is a two-dimensional first-order Gaussian function
  • x 0 and y 0 are the position information of the imaging target; a, b, and c are the mathematical parameters of a two-dimensional first-order Gaussian function;
  • the preset relationship is expressed by the following expression:
  • f(x, y) is the random fluctuation light field distribution
  • t(x, y) is the distribution function of the imaging target
  • I i is the total light intensity
  • i 0 ⁇ N.
  • FIG. 1 schematically shows the structure diagram of a traditional dual-arm intensity correlation imaging system provided in the prior art
  • FIG. 2 schematically shows a structural diagram of a computational ghost imaging system provided by the prior art
  • FIG. 3 schematically shows a method flowchart of an imaging method based on random light intensity fluctuations provided by an embodiment of the present disclosure
  • FIG. 4 schematically shows a structural diagram of an imaging system based on random light intensity fluctuations provided by an embodiment of the present disclosure.
  • An embodiment of the present disclosure provides an imaging method based on random light intensity fluctuations. Referring to FIG. 3 and FIG. 4, the method includes the content of step S1 to step S4:
  • step S1 the light source 1 illuminates the ground glass 2 and forms a random fluctuating light field on the surface of the ground glass 2.
  • the random fluctuating light field illuminates the imaging target 3 and is collected by the single-pixel detector 4 to obtain the total light intensity I 0 .
  • Step S2 rotating the ground glass 2 N times, and repeating step S1 to obtain N total light intensity values, where N ⁇ 1. That is, every time the ground glass 2 is rotated once and the position of the ground glass 2 is changed, a different random fluctuation light field will be formed on the surface of the ground glass 2, so the total light intensity collected by the single-pixel detector 4 will also change accordingly. After rotating the ground glass 2 N times, N total light intensity values will be obtained, I 1 I 2 I 3 I 4 ...I N.
  • Step S3 Establish a distribution function expression of the imaging target 3, and determine the imaging according to the preset relationship among the random fluctuation light field, the N total light intensity values and the distribution function expression The distribution function and location information of target 3.
  • the distribution function expression of the imaging target 3 includes an exponential function, a Fourier function, a trigonometric function, a polynomial, a Gaussian function, or a Weibull function, and the distribution function expression is an M-dimensional K-order function, where M ⁇ 1 ; K ⁇ 1.
  • the distribution function expression of the imaging target 3 is a two-dimensional first-order Gaussian function
  • the distribution function t(x, y) expression is:
  • x 0 and y 0 are the position information of the imaging target 3; a, b, and c are mathematical parameters of a two-dimensional first-order Gaussian function.
  • f(x, y) is the random fluctuation light field distribution
  • t(x, y) is the distribution function of the imaging target 3
  • I i is the total light intensity
  • i 0 ⁇ N.
  • the distribution function and position information of the imaging target 3 are determined ,include:
  • N+1 simultaneous equations are obtained; since the ground glass is rotated to a position, a set of equations can be obtained, and when the ground glass is rotated N times, N equations can be obtained. Combining with an equation when the ground glass is not rotated for the first time, N+1 simultaneous equations are obtained.
  • the values of mathematical parameters a, b, and c in the distribution function of the imaging target 3 and the values of position information x 0 , y 0 are determined. Since all the random fluctuation light field distribution f(x, y) are generated by the same piece of ground glass, the random field generated by the ground glass during one rotation is different, but it will show certain regularity. When solving +1 simultaneous equations, there is no need to know exactly the actual value of the random fluctuation light field distribution f(x, y).
  • the values of the mathematical parameters a, b, and c in the distribution function of the imaging target 3 and the values of the position information x 0 and y 0 can be determined, and the resolution of the system can also be increased by increasing the number of measurements.
  • the size of the imaging target 3 is micrometers and below.
  • Step S4 reconstruct an image of the imaging target 3 according to the distribution function and position information of the imaging target 3.
  • the present disclosure can realize the reconstruction of the imaging target through simple and compact components such as light source, ground glass, imaging target, and single-pixel detector. Without obtaining the actual value of the random fluctuation light field distribution, only By constructing the distribution function of the imaging target, measuring multiple values of total light intensity, and according to the mathematical relationship between the distribution function and multiple values of total light intensity and random fluctuation light field distribution, the imaging target can be achieved reconstruction. It avoids the need to introduce complex optical modules such as digital microlens arrays or projection systems in the prior art, and also solves the problem in the prior art that the actual value of the fluctuating light field distribution needs to be obtained to reconstruct the imaging target.
  • a two-dimensional first-order Gaussian function t(x, y) is used to describe the distribution function of the particle.
  • x 0 and y 0 are the position information of the imaging target (3); a, b, and c are mathematical parameters of a two-dimensional first-order Gaussian function.
  • the values of mathematical parameters a, b, and c in the distribution function of the imaging target and the values of position information x 0 , y 0 are determined. Since all the random fluctuation light field distribution f(x, y) are generated by the same piece of ground glass, the random field generated by the ground glass during one rotation is different, but it will show certain regularity. When solving +1 simultaneous equations, it is not necessary to know the actual value of the random fluctuation light field distribution f(x, y).
  • the values of the mathematical parameters a, b, and c in the distribution function of the imaging target and the values of the position information x 0 , y 0 can be determined, and the resolution of the system can also be increased by increasing the number of measurements.
  • the system includes: a light source 1, ground glass 2, an imaging target 3, and a single-pixel detector 4;
  • the light source 1 illuminates the ground glass 2 and forms a random fluctuation light field on the surface of the ground glass 2.
  • the random fluctuation light field illuminates the imaging target 3 and is collected by the single-pixel detector 4 to obtain the total light intensity. Used to reconstruct the image of the imaging target 3.
  • the imaging system further includes a data processing module configured to reconstruct an image of the imaging target 3 according to the value of the total light intensity.
  • the data processing module may be a computer terminal connected to the single-pixel detector 4.
  • the data processing module is used to obtain N total light intensity values by rotating the ground glass 2 N times, where N ⁇ 1; to establish the distribution function expression of the imaging target 3, according to the random fluctuation
  • the predetermined relationship between the light field, the N total light intensity values and the distribution function expression determines the distribution function and position information of the imaging target 3; according to the distribution function and position information of the imaging target 3 , To reconstruct the image of the imaging target 3.
  • the distribution function expression of the imaging target 3 is a two-dimensional first-order Gaussian function
  • x 0 and y 0 are the position information of the imaging target (3); a, b, and c are the mathematical parameters of a two-dimensional first-order Gaussian function;
  • the preset relationship is expressed by the following expression:
  • f(x, y) is the random fluctuation light field distribution
  • t(x, y) is the distribution function of the imaging target (3)
  • I i is the total light intensity
  • i 0 ⁇ N.
  • N+1 simultaneous equations are obtained; since the ground glass is rotated to a position, a set of equations can be obtained, and when the ground glass is rotated N times, N equations can be obtained. Combining with an equation when the ground glass is not rotated for the first time, N+1 simultaneous equations are obtained.
  • the values of mathematical parameters a, b, and c in the distribution function of the imaging target 3 and the values of position information x 0 , y 0 are determined. Since all the random fluctuation light field distribution f(x, y) are generated by the same piece of ground glass, the random field generated by the ground glass during one rotation is different, but it will show certain regularity. When solving +1 simultaneous equations, there is no need to know exactly the actual value of the random fluctuation light field distribution f(x, y).
  • the values of the mathematical parameters a, b, and c in the distribution function of the imaging target 3 and the values of the position information x 0 and y 0 can be determined, and the resolution of the system can also be increased by increasing the number of measurements.

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Optics & Photonics (AREA)
  • Photometry And Measurement Of Optical Pulse Characteristics (AREA)
  • Length Measuring Devices By Optical Means (AREA)

Abstract

L'invention concerne un procédé et un système d'imagerie basés sur une fluctuation d'intensité de lumière aléatoire. Le procédé comprend les étapes suivantes : une source de lumière (1) irradie du verre rodé (2) et forme un champ lumineux fluctuant aléatoire sur une surface du verre rodé (2), le champ lumineux fluctuant aléatoire éclairant une cible d'imagerie (3) et étant ensuite collecté par un détecteur à pixel unique (4) pour obtenir une intensité lumineuse totale ; rotation du verre rodé (2) N fois pour acquérir N valeurs d'intensité lumineuse totale, avec N ≥ 1 ; établissement d'une expression de la fonction de distribution de la cible d'imagerie (3) et détermination, selon une relation prédéfinie entre le champ lumineux fluctuant aléatoire, les N valeurs d'intensité lumineuse totale et l'expression de la fonction de distribution, d'une fonction de distribution et d'informations de position de la cible d'imagerie (3) ; et rétablissement d'une image de la cible d'imagerie (3) conformément à la fonction de distribution et aux informations de position de la cible d'imagerie (3). Le rétablissement d'une cible d'imagerie est réalisé sans qu'il soit nécessaire d'acquérir une valeur de distribution réelle d'un champ lumineux fluctuant aléatoire.
PCT/CN2019/078183 2019-03-14 2019-03-14 Procédé et système d'imagerie basés sur une fluctuation d'intensité de lumière aléatoire WO2020181560A1 (fr)

Priority Applications (1)

Application Number Priority Date Filing Date Title
PCT/CN2019/078183 WO2020181560A1 (fr) 2019-03-14 2019-03-14 Procédé et système d'imagerie basés sur une fluctuation d'intensité de lumière aléatoire

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
PCT/CN2019/078183 WO2020181560A1 (fr) 2019-03-14 2019-03-14 Procédé et système d'imagerie basés sur une fluctuation d'intensité de lumière aléatoire

Publications (1)

Publication Number Publication Date
WO2020181560A1 true WO2020181560A1 (fr) 2020-09-17

Family

ID=72427776

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/CN2019/078183 WO2020181560A1 (fr) 2019-03-14 2019-03-14 Procédé et système d'imagerie basés sur une fluctuation d'intensité de lumière aléatoire

Country Status (1)

Country Link
WO (1) WO2020181560A1 (fr)

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20060055936A1 (en) * 2004-09-10 2006-03-16 The General Hospital Corporation System and method for optical coherence imaging
CN103777206A (zh) * 2014-01-26 2014-05-07 上海交通大学 一种基于偏振关联成像的单像素成像***

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20060055936A1 (en) * 2004-09-10 2006-03-16 The General Hospital Corporation System and method for optical coherence imaging
CN103777206A (zh) * 2014-01-26 2014-05-07 上海交通大学 一种基于偏振关联成像的单像素成像***

Similar Documents

Publication Publication Date Title
CN109752844B (zh) 一种基于随机光强涨落的成像方法和***
US10136078B2 (en) Random grating based compressive sensing wideband hyperspectral imaging system
EP2035861B1 (fr) Dispositif de localisation et d'imagerie de sources de rayonnement gamma ou x.
KR102254627B1 (ko) 높은 스루풋 및 저비용 높이 삼각측량 시스템 및 방법
JP2019537713A (ja) チップ欠陥検出装置及びチップ欠陥検出方法
CN111458044A (zh) 基于快照光谱成像技术的瞬态温度测量装置和方法
JP2024026464A (ja) 装置および方法
JPH10267636A (ja) 表面検査方法および表面検査装置
CN105738073B (zh) 一种在空间频率域进行像素响应函数测量的方法
CN105067489B (zh) 一种基于动态光散射技术的悬浮颗粒粒径测量装置及方法
US20180293790A1 (en) Method and Apparatus for Measuring 3Dimensional Structures
CN106767523B (zh) 一种提高相位精度的方法及装置
CN112585521A (zh) 用于捕获具有湍流衰减的显微全光图像的设备和方法
WO2020181560A1 (fr) Procédé et système d'imagerie basés sur une fluctuation d'intensité de lumière aléatoire
JPWO2015132831A1 (ja) 撮像装置
FR2965921A1 (fr) Procede de mesure de l'orientation et de la deformation elastique de grains dans des materiaux multicristallins
US20200341125A1 (en) Method of generating a time domain echo waveform and electromagnetic radiation echo waveform generation system
JP2014238299A (ja) 被検物の計測装置、算出装置、計測方法および物品の製造方法
WO2008149003A2 (fr) Dispositif d'imagerie par rayons x à source poly-chromatique
US8836944B2 (en) Spatially resolved imaging of opto-electrical property variations
JP2010197370A (ja) 光応用計測装置
US8355119B2 (en) Method for optical measurement of velocities according to the spatial frequency filtering method and sensor for optical measurement of velocities
Cowsik et al. Design and performance of a sub-nanoradian resolution autocollimating optical lever
CN110533706B (zh) 一种基于计算鬼成像的三维目标重构方法
US20220412731A1 (en) Apparatus and method for quantifying the surface flatness of three-dimensional point cloud data

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 19918780

Country of ref document: EP

Kind code of ref document: A1

NENP Non-entry into the national phase

Ref country code: DE

122 Ep: pct application non-entry in european phase

Ref document number: 19918780

Country of ref document: EP

Kind code of ref document: A1