CN110400282A - A kind of high-resolution terahertz image processing method - Google Patents

A kind of high-resolution terahertz image processing method Download PDF

Info

Publication number
CN110400282A
CN110400282A CN201810371343.7A CN201810371343A CN110400282A CN 110400282 A CN110400282 A CN 110400282A CN 201810371343 A CN201810371343 A CN 201810371343A CN 110400282 A CN110400282 A CN 110400282A
Authority
CN
China
Prior art keywords
matrix
imaging
resolution
spread function
point spread
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.)
Pending
Application number
CN201810371343.7A
Other languages
Chinese (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.)
Shenyang Institute of Automation of CAS
Original Assignee
Shenyang Institute of Automation 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 Shenyang Institute of Automation of CAS filed Critical Shenyang Institute of Automation of CAS
Priority to CN201810371343.7A priority Critical patent/CN110400282A/en
Publication of CN110400282A publication Critical patent/CN110400282A/en
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING; COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T3/00Geometric image transformation in the plane of the image
    • G06T3/40Scaling the whole image or part thereof
    • G06T3/4053Super resolution, i.e. output image resolution higher than sensor resolution
    • GPHYSICS
    • G06COMPUTING; CALCULATING; COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00Image enhancement or restoration
    • G06T5/50Image enhancement or restoration by the use of more than one image, e.g. averaging, subtraction

Abstract

The present invention relates to a kind of high-resolution terahertz image processing methods, including following procedure: obtaining the amplitude and phase matrix of Imaged samples, obtain the amplitude and phase matrix of the point spread function of imaging system;By the complex matrix of the amplitude and phase matrix group compound imaging sample of Imaged samples, by the complex matrix of the amplitude and phase matrix group synthetic point spread function of point spread function;The complex matrix of complex matrix and point spread function to Imaged samples carries out non-blind Deconvolution calculation, obtains full resolution pricture.The present invention is not needed using expensive Near-Field Radar Imaging method, it can use traditional sensors and optical device realize high-resolution imaging, wave characteristic based on electromagnetic field, we introduce phase information to realize that preferable rebuild realizes inexpensive high-resolution imaging in turn, the physical limit for making imaging resolution beyond tradition detector and hot spot greatly expands the performance of imaging system.

Description

A kind of high-resolution terahertz image processing method
Technical field
The present invention relates to terahertz image process field, specifically a kind of high-resolution terahertz image processing method.
Background technique
Terahertz is generally referred to as frequency in the electromagnetic wave of 0.1-10THz frequency range, due to its unique characteristic, in recent years by To people's extensive concern.The vibration of many large biological molecules is at the frequency range in nature, this is the detection of biological nature Provide effective means;Terahertz photon energy is low, will not damage to detection object, can be realized non-destructive testing;Terahertz There is penetration capacity to many dielectric materials and apolar substance, can be used as the means of covered object detection.
Terahertz imaging is the important component of above-mentioned application, it is generally the case that improves terahertz image resolution There are two types of approach: first is that using highly sensitive detector by hardware device;Second is that utilizing figure by founding mathematical models As the method for processing improves the resolution ratio of terahertz image.The former is long the period, restriction expensive and by the prior art, The latter is at low cost, it is easy to accomplish, it has also become the important research direction of terahertz image process field.
Existing its imaging resolution of Terahertz focal plane imaging system depends on the physics of focal beam spot and sensor Size.Although near-field microscope, which can reduce spot size by the use of probe structure, improves imaging resolution, its It is at high cost, and use process middle probe is needed very close to body surface, therefore use environment is limited, it is very inconvenient, simultaneously Since the introducing of probe reduces the dynamic range of system.It is proposed that direction be based on image blur is caused to electromagnetic wave diffraction Physical process understanding, acquire blurred picture using ordinary sensors and carry out image recovery, utilize amplitude and phase information By operation of deconvoluting come reconstruction image.This method does not need to can use traditional sensing using expensive Near-Field Radar Imaging method Device and optical device realize high-resolution imaging.Wave characteristic based on electromagnetic field, we introduce phase information to realize preferably Reconstruction so that realize inexpensive high-resolution imaging, make the physical limit of imaging resolution beyond tradition detector and hot spot, The performance of imaging system is greatly expanded, this is that existing imaging system institute is inaccessiable.
Summary of the invention
In view of the deficiencies of the prior art, the present invention provides a kind of high-resolution terahertz image processing method, using amplitude and Two information completion terahertz images of phase are deconvoluted, and improve its imaging resolution, solve the diffusion of terahertz imaging system point The problem of function estimation inaccuracy.
Present invention technical solution used for the above purpose is:
A kind of high-resolution terahertz image processing method, including following procedure:
Step 1: obtaining the amplitude and phase matrix of Imaged samples, obtain the amplitude and phase of the point spread function of imaging system Bit matrix;
Step 2: by the complex matrix of the amplitude and phase matrix group compound imaging sample of Imaged samples, by point spread function Amplitude and phase matrix group synthetic point spread function complex matrix;
Step 3: the complex matrix of complex matrix and point spread function to Imaged samples carries out non-blind Deconvolution calculation, obtains To full resolution pricture.
The imaging system is transmission-type terahertz imaging system.
The complex matrix of the Imaged samples are as follows: g=A1.*exp(P1/180*π*i)
Wherein, A1For the magnitude matrix for obtaining Imaged samples, P1Represent the phasing matrix for obtaining Imaged samples.
The complex matrix of the point spread function are as follows: k=A2.*exp(P2/180*π*i)
Wherein, A2For the magnitude matrix for obtaining Imaged samples, P2For the phasing matrix for obtaining Imaged samples.
The invention has the following beneficial effects and advantage:
1. improving its imaging resolution by operation of deconvoluting present invention introduces amplitude and phase information come reconstruction image;
2. the present invention applies traditional detection systems, expensive near-field probe part is not needed, use scope at low cost is extensive.
Detailed description of the invention
Fig. 1 is flow chart of the method for the present invention;
Fig. 2 is the point spread function amplitude picture and phase picture of test;
Fig. 3 is the sample amplitude picture and phase picture of test;
Fig. 4 is the image after being deconvoluted using amplitude and phase information;
Specific embodiment
The present invention is described in further detail with reference to the accompanying drawings and embodiments.
It is as shown in Figure 1 flow chart of the method for the present invention.
A kind of high-resolution terahertz image processing method, including following procedure: sample is placed on focal plane, using too Designed Sample Scan is imaged in hertz transmission-type imaging, obtains the amplitude and phase matrix of Imaged samples, Imaged samples Complex matrix are as follows: g=A1.*exp(P1/ 180* π * i), wherein A1Represent the magnitude matrix for obtaining Imaged samples, P1It represents and obtains The phasing matrix of Imaged samples.Obtain the amplitude and phase matrix of the point spread function of transmission-type terahertz imaging system;Point expands Dissipate the complex matrix of function are as follows: k=A2.*exp(P2/ 180* π * i), wherein A2Represent the magnitude matrix for obtaining Imaged samples, P2 Represent the phasing matrix for obtaining Imaged samples.By the plural square of the amplitude and phase matrix group compound imaging sample of Imaged samples Battle array, by the complex matrix of the amplitude and phase matrix group synthetic point spread function of point spread function;To the plural square of Imaged samples The complex matrix of battle array and point spread function carries out non-blind Deconvolution calculation, obtains full resolution pricture.
Using the complex matrix of amplitude and phase matrix group synthetic point spread function and Imaged samples, in conjunction with the side deconvoluted Method, including Constraint least square algorithm filtering or Wiener filtering or Richardson-Lucy algorithm or super Laplace prior or height Full resolution pricture processing is realized in the non-blind deconvolution algorithm such as this priori or sparse prior.
Embodiment 1. constrains least square deconvolution method
In image procossing, image degradation model described in following equations is widely used, such as following formula (1):
G (x, y)=f (x, y) * k (x, y)+n (x, y) (1)
Wherein, g (x, y) is the degraded image that system testing is arrived, and f (x, y) is the clear image to be deconvoluted, and k (x, y) is The point spread function of imaging system, n (x, y) indicate the noise of system.The recovery optimality measured according to smoothness is available The frequency domain solution of the above problem, formula such as following formula (2):
Wherein, γ is adjusting parameter, and P (u, v) is the Fourier transformation of Laplace operator.K (u, v) indicates point spread function Several Fourier transformations, G (u, v) indicate the Fourier transformation of degraded image,That indicate is Fu of image after deconvoluting In leaf transformation.
The present invention uses Terahertz transmission-type imaging system, tests its correspondence image of point spread function magnitude matrix such as Fig. 1 (a) shown in, shown in test its correspondence image of point spread function phasing matrix such as Fig. 1 (b).Testing Imaged samples magnitude matrix, its is right It answers shown in image such as Fig. 2 (a), shown in test its correspondence image of Imaged samples phasing matrix such as Fig. 2 (b).
Specific step is as follows for deconvolution method: the magnitude matrix normalizing first to the point spread function and sample tested Change, gained normalization matrix tells mode according to claims with corresponding phasing matrix respectively and is combined into complex matrix, benefit The full resolution pricture being restored with Constraint least square algorithm deconvolution algorithm, gained full resolution pricture is as shown in figure 3, Fig. 3 (a) Carried out by Imaged samples amplitude imaging figure, Imaged samples phase imaging figure carried out by Fig. 3 (b).

Claims (4)

1. a kind of high-resolution terahertz image processing method, it is characterised in that: including following procedure:
Step 1: obtaining the amplitude and phase matrix of Imaged samples, obtain the amplitude and phase square of the point spread function of imaging system Battle array;
Step 2: by the complex matrix of the amplitude and phase matrix group compound imaging sample of Imaged samples, by the width of point spread function Degree and phasing matrix are combined into the complex matrix of point spread function;
Step 3: the complex matrix of complex matrix and point spread function to Imaged samples carries out non-blind Deconvolution calculation, obtains height Resolution image.
2. high-resolution terahertz image processing method according to claim 1, it is characterised in that: the imaging system is Penetrate formula terahertz imaging system.
3. high-resolution terahertz image processing method according to claim 1, it is characterised in that: the Imaged samples are answered Matrix number are as follows: g=A1.*exp(P1/180*π*i)
Wherein, A1For the magnitude matrix for obtaining Imaged samples, P1Represent the phasing matrix for obtaining Imaged samples.
4. high-resolution terahertz image processing method according to claim 1, it is characterised in that: the point spread function Complex matrix are as follows: k=A2.*exp(P2/180*π*i)
Wherein, A2For the magnitude matrix for obtaining Imaged samples, P2For the phasing matrix for obtaining Imaged samples.
CN201810371343.7A 2018-04-24 2018-04-24 A kind of high-resolution terahertz image processing method Pending CN110400282A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810371343.7A CN110400282A (en) 2018-04-24 2018-04-24 A kind of high-resolution terahertz image processing method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810371343.7A CN110400282A (en) 2018-04-24 2018-04-24 A kind of high-resolution terahertz image processing method

Publications (1)

Publication Number Publication Date
CN110400282A true CN110400282A (en) 2019-11-01

Family

ID=68320118

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810371343.7A Pending CN110400282A (en) 2018-04-24 2018-04-24 A kind of high-resolution terahertz image processing method

Country Status (1)

Country Link
CN (1) CN110400282A (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111579521A (en) * 2020-05-06 2020-08-25 中国科学院沈阳自动化研究所 Terahertz compression imaging optimization method and system based on data selection
CN112014349A (en) * 2020-11-02 2020-12-01 季华实验室 Terahertz time-domain spectral imaging restoration method and device, storage medium and terminal

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101699509A (en) * 2009-11-11 2010-04-28 耿则勋 Method for recovering atmosphere fuzzy remote image with meteorological data
CN104020456A (en) * 2014-06-18 2014-09-03 中国科学院电子学研究所 Linear array imaging radar system amplitude and phase error correction method based on multiple dominant scattering centers
US20140352437A1 (en) * 2013-05-31 2014-12-04 Seiko Epson Corporation Ultrasonic measuring device, ultrasonic image device, and method for processing ultrasonic image
CN106485182A (en) * 2016-06-27 2017-03-08 中国计量大学 A kind of fuzzy Q R code restored method based on affine transformation

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101699509A (en) * 2009-11-11 2010-04-28 耿则勋 Method for recovering atmosphere fuzzy remote image with meteorological data
US20140352437A1 (en) * 2013-05-31 2014-12-04 Seiko Epson Corporation Ultrasonic measuring device, ultrasonic image device, and method for processing ultrasonic image
CN104020456A (en) * 2014-06-18 2014-09-03 中国科学院电子学研究所 Linear array imaging radar system amplitude and phase error correction method based on multiple dominant scattering centers
CN106485182A (en) * 2016-06-27 2017-03-08 中国计量大学 A kind of fuzzy Q R code restored method based on affine transformation

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
DAN C. POPESCU ET AL: "Phantom-Based Point Spread Function Estimation for Terahertz Imaging System", 《INTERNATIONAL CONFERENCE ON ADVANCED CONCEPTS FOR INTELLIGENT VISION SYSTEMS》 *
肖嵩等: "《计算机图形学原理及应用》", 30 June 2014 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111579521A (en) * 2020-05-06 2020-08-25 中国科学院沈阳自动化研究所 Terahertz compression imaging optimization method and system based on data selection
CN112014349A (en) * 2020-11-02 2020-12-01 季华实验室 Terahertz time-domain spectral imaging restoration method and device, storage medium and terminal

Similar Documents

Publication Publication Date Title
Ercius et al. Electron tomography: a three‐dimensional analytic tool for hard and soft materials research
Pierre et al. Multiscale reconstruction for MR fingerprinting
Eriksson et al. Isotropic diffusion weighting in PGSE NMR by magic-angle spinning of the q-vector
Sang et al. Revolving scanning transmission electron microscopy: Correcting sample drift distortion without prior knowledge
Davis et al. Using a light microscope to measure motions with nanometer accuracy
Britton et al. Factors affecting the accuracy of high resolution electron backscatter diffraction when using simulated patterns
Merlet et al. Continuous diffusion signal, EAP and ODF estimation via compressive sensing in diffusion MRI
Amineh et al. Three-dimensional near-field microwave holography using reflected and transmitted signals
CN102197302B (en) Analysis method, radiation imaging apparatus using analysis method, and analysis program for executing analysis method
Poupon et al. New diffusion phantoms dedicated to the study and validation of high‐angular‐resolution diffusion imaging (HARDI) models
Mougin et al. Magnetization transfer phenomenon in the human brain at 7 T
Faris et al. Three-dimensional beam-deflection optical tomography of a supersonic jet
Jaeger et al. Evaluation of 1H NMR relaxometry for the assessment of pore‐size distribution in soil samples
Chang et al. Nonlinear GRAPPA: A kernel approach to parallel MRI reconstruction
Di Ianni et al. Compressive sensing of full wave field data for structural health monitoring applications
US10551520B2 (en) Surface relaxivity calculation using nuclear magnetic resonance (NMR) measurement, three dimensional (3D) rock model and NMR response simulation
Nguyen et al. T2prep three‐dimensional spiral imaging with efficient whole brain coverage for myelin water quantification at 1.5 tesla
Diebolder et al. Pushing the resolution limits in cryo electron tomography of biological structures
Latham et al. Image registration: enhancing and calibrating X-ray micro-CT imaging
Hild et al. Three-dimensional analysis of a compression test on stone wool
Heo et al. Accelerating chemical exchange saturation transfer (CEST) MRI by combining compressed sensing and sensitivity encoding techniques
Song et al. Least-square NUFFT methods applied to 2-D and 3-D radially encoded MR image reconstruction
CN104939828B (en) MR imaging method
CN102928796A (en) Fast-diffused magnetic resonance imaging and restoring method
Young Quantitative microscopy

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