IP Library Granted Patent US 12,579,803
Granted Patent B2
US 12,579,803 · App. 18/162,180 · Granted Mar 17, 2026

Totagraphy for superresolution imaging and signal processing of positive, real-valued images and signals

Inventor: Okan Ersoy (West Lafayette, IN)
Assignee: WAVEFRONT ANALYSIS SYSTEMS LLC
G06V10/898G06V10/25
View Patent ↗
Loading inventors, assignments & file history…
Monitor This Case
Get email alerts when status or documents change.
Order Certified Copies
Most orders are placed with the USPTO same day — all within 24 business hours.
Order via The Patent Place →
Pre-filled with this patent's details
Quick Facts
Patent No.
US 12,579,803
App. No.
18/162,180
Granted
Mar 17, 2026
Kind
B2
Abstract

A method for recovering missing spectral amplitude and phase information of an input real, positive valued image or signal possibly carried by incoherent or coherent waves.

Claims (21)

1 . A method for recovering missing spectral amplitude and phase information of an input signal, the method comprising:

providing at least one transformation unit with an input and an output;

providing at least three masks or illumination patterns, comprising a first transparent mask and a complementary pair of unipolar masks or illumination patterns, each of the at least three masks or illumination patterns configured to be disposed at the input;

separately multiplying the input with each of the at least three masks or illumination patterns to generate a modified input from each of the masks or illumination patterns;

convolving the modified inputs with a transfer function to yield outputs comprising positive real values to be detected by an output sensor, the transfer function comprising at least one invertible region and at least one non-invertible region;

passing the outputs through the at least one transformation unit, the at least one transformation unit configured to produce a spectral representation of each output corresponding to each modified input;

defining two regions of a recovered output comprising a plurality of points, divided into Region 1 in which the transfer function has known values which are invertible, and Region 2 in which the transfer function has very small, unknown or corrupted values;

computing amplitude and phase values of the points corresponding to Region 1 in each spectral representation;

setting an amplitude at each point in Region 2 of each spectral representation to be the mean amplitude computed over all points in Region 1 of the corresponding spectral representation, and setting a corresponding phase to a random value in [0,2π] for each point in Region 2 of the corresponding spectral representation to produce a set of modified spectral representations;

calculating an Inverse Fourier Transform of each of the modified spectral representations to obtain a set of candidate intermediate outputs;

averaging the candidate intermediate outputs to yield a single averaged candidate intermediate output;

providing the single averaged candidate intermediate output at an input plane and applying each of the at least three masks or illumination patterns to generate modified inputs from each of the masks or illumination patterns;

repeating the above steps, iteratively producing new candidate intermediate outputs by recalculating the values in Region 2 and replacing Region 1 with the original Region 1 until the candidate intermediate output values converge, yielding the recovered output.

2 . The method of claim 1 , wherein the complementary pair of unipolar masks or illumination patterns is a phase mask.

3 . The method of claim 2 , wherein the complementary pair of unipolar masks is replaced by a bipolar binary mask with elements equal to +1 or −1.

4 . The method of claim 1 , wherein the input signal is an image.

5 . The method of claim 4 , wherein the image is a real, positive valued image.

6 . The method of claim 1 , wherein the complementary pair of unipolar masks or illumination patterns is a set of illumination patterns.

7 . The method of claim 1 , further comprising the step of performing a Fourier transform of the outputs to produce the spectral representation of each output.

8 . The method of claim 1 , wherein the transfer function is an optical transfer function.

9 . The method of claim 1 , wherein the transfer function is a coherent transfer function.

Assignments (2)
ASSIGNMENT OF ASSIGNOR'S INTEREST Recorded Jul 17, 2023
From: GERCHBERG OPHTHALMIC DISPENSING, PLLC
To: WAVEFRONT ANALYSIS SYSTEMS LLC
Reel/Frame 064279/0558 →
ASSIGNMENT OF ASSIGNOR'S INTEREST Recorded Jul 17, 2023
From: ERSOY, OKAN
To: GERCHBERG OPHTHALMIC DISPENSING, PLLC
Reel/Frame 066239/0498 →
Continuity (2)
Provisional Application 63305302 · Feb 1, 2022
Related Publication 20230245440A1 · Aug 3, 2023
References Cited (41)
US 5050220A · Marsh · 1991 [cited by applicant]
US 5086222A · Shibuya · 1992 [cited by applicant]
US 5274716A · Mitsuoka · 1993 [cited by applicant]
US 5426521A · Chen · 1995 [cited by applicant]
US 5454047A · Chang · 1995 [cited by applicant]
US 5768242A · Juday · 1998 [cited by applicant]
US 6097856A · Hammond, Jr. · 2000 [cited by applicant]
US 6229649B1 · Woods · 2001 [cited by applicant]
US 6369932B1 · Gerchberg · 2002 [cited by applicant]
US 6545790B2 · Gerchberg · 2003 [cited by applicant]
US 6885442B1 · Nugent · 2005 [cited by applicant]
US 6906839B2 · Gerchberg · 2005 [cited by applicant]
US 8040595B2 · Gerchberg · 2011 [cited by applicant]
US 8520213B2 · Popescu · 2013 [cited by applicant]
US 8837045B2 · Popescu · 2014 [cited by applicant]
US 9052180B2 · Popescu · 2015 [cited by applicant]
US 9404857B2 · Popescu · 2016 [cited by applicant]
US 9715098B2 · Babacan · 2017 [cited by applicant]
US 10132609B2 · Popescu · 2018 [cited by applicant]
US 20020060831A1 · Gerchberg · 2002 [cited by applicant]
US 20170059845A1 · Waller · 2017 [cited by applicant]
US 20170146788A1 · Waller · 2017 [cited by applicant]
US 20180329191A1 · Cheng · 2018 [cited by examiner]
US 20190107655A1 · Waller · 2019 [cited by applicant]
US 20190227490A1 · Waller · 2019 [cited by examiner]
US 20190310374A1 · Gerchberg · 2019 [cited by applicant]
Acquire mass and volume information in real time with label-free, quantitative imaging for live cells, assays, tissues and organoids, Phioptics, https://phioptics.com, 3 pages, 2020. [cited by applicant]
Bourquard et al. “A practical inverse-problem approach to digital holographic reconstruction,” Optics Express, vol. 21, No. 3, pp. 3417-3433, Feb. 11, 2013. [cited by applicant]
Chandra et al. “PhasePack User Guide,” https://github.com/tomgoldstein/phasepack-matlab, pp. 1-20, Nov. 29, 2017. [cited by applicant]
Cheng et al. “Phase retrieval and diffractive imaging based on Babinet's principle and complementary random sampling,” Optics Express, vol. 23, No. 22, pp. 28874-28882, Nov. 2, 2015. [cited by applicant]
Eguchi et al. “Optimization of random phase diversity for adaptive optics using an LCoS spatial light modulator,” Applied Optics, vol. 58, No. 25, pp. 6834-6840, Sep. 1, 2019. [cited by applicant]
Eguchi et al. “Single-shot phase retrieval with complex diversity,” Optics Letters, vol. 44, Issue 21, pp. 5108-5112, 2019. [cited by applicant]
Gerchberg “A New Approach to Phase Retrieval of a Wave Front,” Journal of Modern Optics,vol. 49, No. 7, pp. 1185-1196, 2002. [cited by applicant]
Gerchberg “Super-Resolution Through Error Energy Reduction,” Optica ACTA, vol. 21, No. 9, pp. 709-720, 1974. [cited by applicant]
Gerchberg, et al. “A Practical Algorithm for the Determination of Phase from Image and Diffraction Plane Pictures,” Optik, vol. 35, No. 2, pp. 237-246, 1972. [cited by applicant]
Horisaki, et al. “Single-shot phase imaging with a coded aperture,” Optics Letters, vol. 39, No. 22, pp. 6466-6469, Nov. 15, 2014. [cited by applicant]
Pozzi et al. “Fast Calculation of Computer Generated Holograms for 3D Photostimulation through Compressive-Sensing Gerchberg-Saxton Algorithm,” Methods and Protocols, pp. 1-11, 2019. [cited by applicant]
Yuan et al. “Phase Retrieval via Reweighted Wirtinger Flow,” Appl. Opt., vol. 56, No. 9, pp. 1-21, Mar. 2017. [cited by applicant]
Yurtsever et al. “Sketchy Decisions: Convex Low-Rank Matrix Optimization with Optimal Storage,” pp. 1-16, Feb. 22, 2017. [cited by applicant]
Zeng et al. “Coordinate Descent Algorithms for Phase Retrieval,” IEEE, pp. 1-13, Jun. 2017. [cited by applicant]
Zheng et al. “Wide-field, High-resolution Fourier Ptychographic Microscopy,” Nature Photonics, pp. 739-745, vol. 7, No. 9, Sep. 1, 2013. [cited by applicant]