CN104990514A - Data processing apparatus and method for dynamic Fourier transform profilometry - Google Patents
Data processing apparatus and method for dynamic Fourier transform profilometry Download PDFInfo
- Publication number
- CN104990514A CN104990514A CN201510399889.XA CN201510399889A CN104990514A CN 104990514 A CN104990514 A CN 104990514A CN 201510399889 A CN201510399889 A CN 201510399889A CN 104990514 A CN104990514 A CN 104990514A
- Authority
- CN
- China
- Prior art keywords
- fourier transform
- module
- phase
- data
- 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.)
- Pending
Links
Abstract
The invention relates to a data processing apparatus and method for dynamic Fourier transform profilometry. The data processing apparatus comprises an image receiving unit, a frequency domain filtering processing unit, a data buffer unit, a phase-height calculation unit, and a communication unit. The image receiving unit connected with imaging equipment is used for receiving collected raster image data; the frequency domain filtering processing unit is used for carrying out Fourier transform on a raster image, completing frequency domain filtering, and extracting a fundamental component only containing to-be-measured object height information; the data buffer unit is used for carrying out data buffer and address conversion in frequency domain filtering as well as data buffer of external communication; the phase-height calculation unit is used for solving image phase distribution and calculating height distribution based on the phase-height mapping; and the communication unit is used for data transmission and communication with an external CPU. According to the invention, the computer hardware technology is applied to data processing; and the apparatus and method have advantages of fast processing speed, high integration degree, small size, low cost, and good support of real-time measurement.
Description
Technical field
The present invention discloses a kind of dynamic Fourier transform profilometry data processing method, relate to optical three-dimensional measurement field, be related specifically to and computer hardware technology is applied to dynamic Fourier transform profilometry Data processing, to improve data processing speed, support to measure in real time.
Background technology
In scientific research and commercial production, often need to carry out tri-dimensional facial type measurement to the object of observation of the dynamic processes such as such as High Rotation Speed, fluid mechanics, sport dynamics, ballistics, detonation process, impact process, to describe and to analyze the change of testee configuration of surface in dynamic process, extract the physical parameter such as structure, deformation, stress relevant to testee.In optics initiatively three-dimensional sensing measuring technique, Fourier transform profilometry (Fourier transform profilometry, be called for short FTP) only need a frame deforming stripe figure just can recover the tri-dimensional facial type of object, the tri-dimensional facial type being particularly suitable for dynamic process is measured.The light field that sinusoidal grating or Ronchi grating produce is projected to surface of moving object to be measured by dynamic process FTP, imaging system shooting is by a series of deformed grating figure of measured object dignity shape modulation, then utilize computing machine to carry out the process such as two-dimensional Fourier transform, spectral filtering, inverse Fourier transform and phase calculation to deformed grating figure, obtain the tri-dimensional facial type distribution of tested moving object through phase height mapping.The data handling procedure of dynamic process FTP comprises two-dimensional Fourier transform, spectral filtering, two-dimentional inverse Fourier transform, phase calculation, phase unwrapping and phase height mapping etc., and data calculated amount is large.At present, classic method adopts computing machine pure software mode to process usually, and data processing speed is slow, cannot realize the process real-time of data.Especially for the raster image of high-res and high frame per second, because gathered data volume is large, be ensure data processing speed, need use high performance computing machine, equipment is huge and cost is high.
Summary of the invention
The defect existed in the technology of data processing is carried out for existing employing pure software, the invention provides a kind of dynamic Fourier transform profilometry data processing method, computer hardware technology is applied to Data processing, possess that processing speed is fast, integrated level is high, volume is little, and cost is low, support the advantage measured in real time.
A kind of dynamic Fourier transform profilometry data processing equipment and method, comprise image receiving unit, frequency domain filtering processing unit, data buffer unit, phase height computing unit, communication unit.
Image receiving unit is used for being connected with imaging device, receives the raster image data collected;
Frequency domain filtering processing unit is used for carrying out two-dimensional Fourier transform to raster image, completes frequency domain filtering, extracts the fundamental component only comprising determinand elevation information;
Data buffer unit is used for data buffering in frequency domain filtering and address mapping, and with the data buffering of PERCOM peripheral communication;
Phase height computing unit is used for solving image phase distribution, calculates the distribution of determinand height according to phase height mapping;
Communication unit is used for transmitting with the data of outer CPU and communicating.
Described image receiving unit comprises camera interface module; Frequency domain filtering processing unit comprises two-dimensional Fourier transform module, two-dimensional frequency filtration module, two-dimentional inverse Fourier transform module; Data buffer unit comprises intermediate data buffer module, data output buffer module; Described phase height computing unit comprises phase calculation module, phase unwrapping module, height mapping module; Described communication unit comprises high-speed communication module.
Camera interface module is connected with imaging device; Two-dimensional Fourier transform module is connected with two-dimensional frequency filtration module, and two-dimensional frequency filtration module is connected with intermediate data buffer module, two-dimensional frequency filtration module and two-dimentional inverse Fourier transform model calling; Two dimension inverse Fourier transform module and phase calculation model calling, phase calculation module and phase unwrapping model calling, phase unwrapping module and height mapping model calling; Height mapping module and data output buffer model calling; Data output buffer module and high-speed communication model calling; High-speed communication module is connected with outer CPU.
Described camera interface module, for receiving the raster image data collected from external imaging device.
Described two-dimensional Fourier transform module, for calculating the two-dimentional Fourier spectrum of raster image;
Two-dimensional Fourier transform process comprises: adopt N number of Streaming-FFT to form parallel Pipeline and carry out one-dimensional Fourier transform by row to image, obtain one-dimensional Fourier transform result by row; Adopting N number of Streaming-FFT to form parallel Pipeline again carries out by row one-dimensional Fourier transform to described one-dimensional Fourier transform result by row, obtains the two-dimentional Fourier spectrum of image;
Two-dimensional frequency filtration module, for extracting the fundamental component only comprising determinand elevation information in Fourier spectrum;
Two-dimensional frequency filtering comprises: according to the filter factor in the type computed image of selected wave filter corresponding to each pixel, be multiplied by the filter coefficient of described correspondence with the frequency spectrum of image, obtain filtered fundamental component;
Two dimension inverse Fourier transform module, is transformed into spatial domain for Fourier spectrum filtering extracted, obtains two-dimentional spatial domain picture;
Two dimension inverse Fourier transform process comprises: adopt N number of Streaming-IFFT to form parallel Pipeline and carry out one dimension inverse Fourier transform by row to the frequency spectrum that filtering is extracted, obtain one dimension inverse Fourier transform result by row; Adopt N number of Streaming-IFFT to form parallel Pipeline again to carry out, by the one dimension inverse Fourier transform of row, obtaining two-dimentional spatial domain picture to one dimension inverse Fourier transform result by row.
Described intermediate data buffer module, for the data buffering in two-dimensional frequency filtering and address mapping;
Described data output buffer module, for the data buffering with PERCOM peripheral communication.
Described phase calculation module, for the spatial domain picture information exported frequency domain filtering unit.Be tied to the transformational relation of polar coordinate system according to Cartesian plane coordinate, calculate corresponding original phase hum pattern, namely block PHASE DISTRIBUTION;
Described phase unwrapping module, for launching there is the phase range blocked in described original phase hum pattern, obtains continuous print phase information figure;
Described height mapping module, for the distribution of described continuous phase, distributes according to phase height mapping relation computed altitude.
Described high-speed communication module, for transmitting with the data of outer CPU and communicating.
A kind of dynamic Fourier transform profilometry data processing method, is applied to dynamic Fourier transform profilometry Data processing by computer hardware technology.Hardware development plate for data processing can be any one in FPGA, DSP, MCU.
The present invention discloses a kind of dynamic Fourier transform profilometry data processing equipment and method, and technique effect is as follows:
1), relative to traditional employing computing machine pure software data processing method, computer hardware technology is applied to dynamic Fourier transform profilometry Data processing by the present invention, adopt the data-interface of the calculating of hardware full parellel and high bandwidth, improve data processing speed, the Quick Measurement of dynamic 3 D face type can be realized or measure in real time.
2), relative to traditional data handling procedure, total data processing procedure focuses on a hardware development plate and completes by the present invention, and level of integrated system improves, and volume reduces, and can be used for portable type measuring.
3), the view data large to the data volume such as high-res and high frame rate be when processing, and no longer needs high-performance computer, reduce energy consumption and measure cost.
4), dynamic Fourier transform profilometry total data processing procedure is focused on a hardware development plate complete.Computer hardware technology is applied to dynamic Fourier transform profilometry Data processing, adopts hardware to replace software to complete data processing, improve data processing speed, reduce equipment cost.
Accompanying drawing explanation
Fig. 1 is dynamic Fourier transform profilometry light path schematic diagram.
Fig. 2 is apparatus module connection layout of the present invention.
Fig. 3 is the aluminium matter thin flat plate system and device schematic diagram of measuring vibrations of the present invention.
Fig. 4 is method flow diagram of the present invention.
Fig. 5 is the aluminium matter thin flat plate that vibrates in the embodiment of the present invention 3 frames not deforming stripe figure in the same time when t=0.75s, 0.85s, 1s.
The Fourier spectrum distribution that when Fig. 6 is t=1s in the embodiment of the present invention, deforming stripe figure is corresponding.
Fig. 7 is wave filter set in the embodiment of the present invention.
Fig. 8 is to the fundamental component extracted after the Fourier spectrum filtering of deforming stripe figure during t=1s in the embodiment of the present invention.
Fig. 9 be in the embodiment of the present invention t=1s moment block PHASE DISTRIBUTION.
Figure 10 is the continuous phase distribution in t=1s moment in the embodiment of the present invention.
The aluminium matter thin flat plate tri-dimensional facial type grid distribution plan that when Figure 11 is t=0.75s, 0.85s, 1s of restoration and reconstruction in the embodiment of the present invention, 3 frames do not vibrate in the same time.
Embodiment
Below in conjunction with accompanying drawing, principle of work and embodiment, the present invention is described in further detail.
The Fourier transform profilometry measuring system index path of dynamic process is similar to common static Fourier transform profilometry index path, and its projecting light path as shown in Figure 1.6 is optical projection systems, E'
pe
pbe its optical axis, 7 is imaging systems, E'
ce
cits optical axis, two optical axises intersect at the some O in reference planes 8, and 9 is moving object to be measured, and t representative is not in the same time.Imaging system is l to the distance of testee
0, the distance between imaging system and optical projection system is d, and imaging system generally adopts ccd video camera.
By a sine (or Luo Qi) optical grating projection in reference planes 8, record its light distribution by CCD; And then project to object under test 9 surface, when object is in dynamic changing process, with a series of deforming stripe image of CCD fast recording, the 3 d shape distribution of dynamic object therefrom just can be recovered.
When adopting sinusoidal grating projection, the light distribution in reference planes is:
g
0(x,y)＝a(x,y)+b(x,y)cos[2πf
0x+φ
0(x,y)] (1)
In formula, a (x, y) represents background light intensity, and b (x, y) represents reference planes non-uniform reflection rate, f
0for the spatial frequency of grating image, φ
0(x, y) is the initial phase distribution of grating image.
The intensity distributions of the deforming stripe figure making a video recording obtained fast by CCD is:
g(x,y,t)＝a(x,y,t)+b(x,y,t)cos[2πf
0x+φ(x,y,t)]t＝t
1,t
2,…T (2)
Wherein t represents t frame, corresponding with the time, and T is total shooting frame number of whole dynamic process, a (x, y, t) background distributions when representing t frame, b (x, y, body surface reflectivity heterogeneous when t) representing t frame, φ (x, y, t) be the phase place of modulating by object height, comprise the elevation information of object.
In order to eliminate the impact of body surface non-uniform reflection rate, obtain separately the value of φ (x, y, t), order:
According to Euler's formula, (2) formula can be rewritten as:
g(x,y,t)＝a(x,y,t)+q(x,y,t)exp(j2πf
0x)+q
*(x,y,t)exp(-j2πf
0x) (4)
In formula, * represents conjugate operation, carries out two-dimensional Fourier transform, obtain the Fourier spectrum of striped to (4) formula:
G(f
x,f
y,t)＝A(f
x,f
y,t)+Q(f
x-f
0,f
y,t)+Q
*(f
x+f
0,f
y,t) (5)
Wherein G (f
x, f
y,, t), A (f
x, f
y, t), Q (f
x, f
y, t) be the two-dimensional Fourier transform frequency spectrum of g (x, y, t), a (x, y, t), q (x, y, t) respectively.By an applicable two-dimentional bandpass filter (as rectangular window, Hanning window or hamming window), filtered out by the fundamental component of wherein, then carry out inverse Fourier transform to the fundamental component separated, answering of obtaining is distributed as:
Similar conversion and filtering process are done to light distribution (1) formula in reference planes, can obtain:
Phase place change △ φ (x, y, t) caused by object height can be obtained by lower column operations:
In formula, Im{} represents the imaginary part of getting plural number.
Utilize the geometric relationship in Fig. 1 just can calculate height distribution h (x, y, t) of object:
h(x,y,t)＝l
0△φ(x,y,t)/[△φ(x,y,t)-2πf
0d]≈-l
0△φ(x,y,t)/2πf
0d (9)
Phase value is obtained by arctan function (arctangent) computing, and the phase value calculated is truncated in the main value scope (-π, π) of inverse trigonometric function, for obtaining continuous phase distribution, must carry out phase unwrapping.In the three dimension profile measurement of dynamic process, the deforming stripe intensity distributions that CCD obtains is not only x, the function of y, and the function being time t, is a three dimensional Phase field blocked by calculating the △ φ (x, y, t) obtained.Not only to carry out in x y two dimensional surface during the expansion of three dimensional Phase, but also should carry out on time shaft t direction, to ensure the continuity of adjacent moment phase face.
In order to complete above-mentioned measuring method, the invention provides a kind of dynamic Fourier transform profilometry data processing equipment, computer hardware technology is applied to Data processing, realize Quick Measurement or measure in real time, apparatus module connection layout as shown in Figure 2.
Image receiving unit 1, for: be connected with the imaging device such as camera, receive the raster image data collected.
Frequency domain filtering processing unit 2, for: Fourier transform is carried out to raster image, completes frequency domain filtering, extract the fundamental component only comprising determinand elevation information.
Data buffer unit 3, for the data buffering in frequency domain filtering and address mapping, and with the data buffering of PERCOM peripheral communication.
Phase height computing unit 4, for: solve image phase distribution, calculate the distribution of determinand height according to phase height mapping.
Communication unit 5, for: transmit with the data of outer CPU and communicate.
Described image receiving unit 1 comprises camera interface module 1a; Frequency domain filtering processing unit 2 comprises two-dimensional Fourier transform module 2a, two-dimentional inverse Fourier transform module 2b, two-dimensional frequency filtration module 2c; Data buffer unit 3 comprises intermediate data buffer module 3a, data output buffer module 3b; Described phase height computing unit 4 comprises phase calculation module 4a, phase unwrapping module 4b, height mapping module 4c; Described communication unit 5 comprises high-speed communication module 5a.
Camera interface module 1a is connected with imaging devices such as external camera; Two-dimensional Fourier transform module 2a is connected with two-dimensional frequency filtration module 2c, and two-dimensional frequency filtration module 2c is connected with intermediate data buffer module 3a, and two-dimensional frequency filtration module 2c is connected with two-dimentional inverse Fourier transform module 2b; Two dimension inverse Fourier transform module 2b is connected with phase calculation module 4a, and phase calculation module 4a is connected with phase unwrapping module 4b, and phase unwrapping module 4b is connected with height mapping module 4c; Height mapping module 4c is connected with data output buffer module 3b; Data output buffer module 3b is connected with high-speed communication module 5a; High-speed communication module 5a is connected with outer CPU.
Camera interface module 1a, for: the raster image data collected is received from imaging devices such as external camera.
Two-dimensional Fourier transform module 2a, for: the two-dimentional Fourier spectrum calculating raster image.
Two-dimensional Fourier transform process comprises: adopt N number of Streaming-FFT to form parallel Pipeline and carry out one-dimensional Fourier transform by row to image, obtain one-dimensional Fourier transform result by row; Adopting N number of Streaming-FFT to form parallel Pipeline again carries out by row one-dimensional Fourier transform to described one-dimensional Fourier transform result by row, obtains the two-dimentional Fourier spectrum of image.
Two-dimensional frequency filtration module 2c, for: extract the fundamental component only comprising elevation information in Fourier spectrum.
Two-dimensional frequency filtering comprises: according to the filter factor in the type computed image of selected wave filter corresponding to each pixel; The filter coefficient of described correspondence is multiplied with the frequency spectrum of image, obtains filtered frequency spectrum.
Two dimension inverse Fourier transform module 2b, for: image spectrum filtering extracted is transformed into spatial domain, obtains two-dimentional spatial domain picture.
Two dimension inverse Fourier transform process comprises: adopt N number of Streaming-IFFT to form parallel Pipeline and carry out one dimension inverse Fourier transform by row to image, obtain one dimension inverse Fourier transform result by row; Adopt N number of Streaming-IFFT to form parallel Pipeline again to carry out, by the one dimension inverse Fourier transform of row, obtaining two-dimentional spatial domain picture to one dimension inverse Fourier transform result by row.
Intermediate data buffer module 3a, for: the data buffering in two-dimensional frequency filtering and address mapping.
Data output buffer module 3b, for: with the data buffering of PERCOM peripheral communication.
Phase calculation module 4a, for: the spatial domain picture information exported frequency domain filtering unit, is tied to the transformational relation of polar coordinate system according to Cartesian plane coordinate, calculate corresponding original phase hum pattern.
Phase unwrapping module 4b, for: launching there is the phase range blocked in described original phase hum pattern, obtaining continuous print phase information figure.
Height mapping module 4c, for: to the distribution of described continuous phase, distribute according to Phase-height mapping relation computed altitude.
High-speed communication module 5a, for: transmit with the data of outer CPU and communicate.
Adopt described method and the aluminium matter thin flat plate of device to one piece of vibration to carry out tri-dimensional facial type measurement below, Fig. 3 is measuring system device schematic diagram, and 10 is projection arrangements, 11 is CCD imaging systems, 12 is vibration aluminium matter thin flat plates to be measured, and 13 is the hardware development plates for data processing, and 14 is PCs.Be FPGA development board for the hardware of data processing in this example, it is the SP605 development board of Xilinx company, adopts Xilinx Spartan-6XC6SL45T FPGA.
The long 280mm of tested aluminium matter thin flat plate, wide 150mm; Imaging system CCD is LumeneraCameraLu050M, and lens focus is 12mm; The filming frequency of CCD is 200fps.Optical projection system emergent pupil is about 150mm to the distance of imaging system entrance pupil, and optical projection system emergent pupil is about 800mm to the distance of plane, knocks the back side of thin flat plate while starting to record image with iron hammer, the vibration processes of shooting thin flat plate.Whole process have taken altogether 1.5 seconds, obtains 300 frame dynamic images altogether.Through preliminary images analysis, select the 50 two field picture intercepting parts wherein characterizing plate vibration process and process, rebuild the vibration processes of thin flat plate.
Idiographic flow step is as shown in Figure 4:
S1: initialization is carried out to system.
The structural parameters of recording measuring device and the inside and outside portion parameter of imaging system, arrange the parameter such as parameter and phase unwrapping reference point location of raster image processing region, filtering, these parameters be sent to FPGA.
S2: receive raster image data from imaging devices such as cameras.
Fig. 5 is wherein 3 frames not grating deforming stripe figure in the same time, corresponds respectively to t=0.75s, 0.85s, 1s.
S3: the Pipeline utilizing 2N Streaming-FFT composition to walk abreast carries out two-dimensional Fourier transform to the raster image received, obtains Fourier spectrum distribution.
The two-dimentional Fourier spectrum distribution of raster image when Fig. 6 is t=1s.
S4: carry out frequency domain filtering to the two-dimentional Fourier spectrum of image, extracts the fundamental component only comprising elevation information in frequency spectrum.
Filtering comprises: according to the filter factor in the filter type computed image selected corresponding to each pixel; The filter coefficient of described correspondence is multiplied with the frequency spectrum of image, obtains filtered frequency spectrum.
Fig. 7 is set wave filter, and Fig. 8 is the fundamental component to extracting after raster image Fourier spectrum filtering during t=1s.
S5: the Pipeline utilizing 2N Streaming-IFFT composition to walk abreast carries out two-dimentional inverse Fourier transform to Fourier's fundamental component that filtering extracts, and obtains two-dimentional spatial domain picture.
S6: carry out phase calculation to the spatial domain picture that inverse Fourier transform obtains, what obtain image blocks PHASE DISTRIBUTION.
Fig. 9 be the t=1s moment block PHASE DISTRIBUTION.
S7: launch blocking phase place, obtains the distribution of continuous print phase information.
Figure 10 is the continuous phase distribution in t=1s moment.
S8: according to phase height mapping relation, obtains the distributed intelligence of determinand height from continuous phase information.
S9: will height distributed intelligence input outer CPU.
S10: display three-dimensional height information.
Figure 11 is t=0.75s, 0.85s, 1s shown by external display tri-aluminium matter thin flat plate tri-dimensional facial type grid distribution plan of not vibrating in the same time of frame.
Dynamic Fourier transform profilometry total data processing procedure focuses on a hardware development plate and completes by the present invention.Computer hardware technology is applied to dynamic Fourier transform profilometry Data processing, adopts hardware to replace software to complete data processing, improve data processing speed, reduce equipment cost.Hardware development plate used is any one in FPGA, DSP, MCU.
Claims (10)
1. a dynamic Fourier transform profilometry data processing equipment, comprise image receiving unit (1), frequency domain filtering processing unit (2), data buffer unit (3), phase height computing unit (4), communication unit (5), it is characterized in that
Image receiving unit (1), for being connected with imaging device, receives the raster image data collected;
Frequency domain filtering processing unit (2), for carrying out Fourier transform to raster image, completes frequency domain filtering, extracts the fundamental component only comprising determinand elevation information;
Data buffer unit (3) for the data buffering in frequency domain filtering and address mapping, and with the data buffering of PERCOM peripheral communication;
Phase height computing unit (4), for solving image phase distribution, distributes according to phase height mapping computed altitude;
Communication unit (5) is for transmitting with the data of outer CPU and communicating.
2. a kind of dynamic Fourier transform profilometry data processing equipment according to claim 1, it is characterized in that, described image receiving unit (1) comprises camera interface module (1a); Frequency domain filtering processing unit (2) comprises two-dimensional Fourier transform module (2a), two-dimentional inverse Fourier transform module (2b), two-dimensional frequency filtration module (2c); Data buffer unit (3) comprises intermediate data buffer module (3a), data output buffer module (3b); Described phase height computing unit (4) comprises phase calculation module (4a), phase unwrapping module (4b), height mapping module (4c); Described communication unit (5) comprises high-speed communication module (5a).
Camera interface module (1a) is connected with imaging device; Two-dimensional Fourier transform module (2a) is connected with two-dimensional frequency filtration module (2c), two-dimensional frequency filtration module (2c) is connected with intermediate data buffer module (3a), and two-dimensional frequency filtration module (2c) is connected with two-dimentional inverse Fourier transform module (2b); Two dimension inverse Fourier transform module (2b) is connected with phase calculation module (4a), phase calculation module (4a) is connected with phase unwrapping module (4b), and phase unwrapping module (4b) is connected with height mapping module (4c); Height mapping module (4c) is connected with data output buffer module (3b); Data output buffer module (3b) is connected with high-speed communication module (5a); High-speed communication module (5a) is connected with outer CPU.
3. a kind of dynamic Fourier transform profilometry data processing equipment according to claim 2, is characterized in that, described camera interface module (1a), for receiving the raster image data collected from external imaging device.
4. a kind of dynamic Fourier transform profilometry data processing equipment according to claim 2, is characterized in that, described two-dimensional Fourier transform module (2a), for calculating the two-dimentional Fourier spectrum of raster image;
Two-dimensional Fourier transform process comprises: adopt N number of Streaming-FFT to form parallel Pipeline and carry out one-dimensional Fourier transform by row to image, obtain one-dimensional Fourier transform result by row; Adopting N number of Streaming-FFT to form parallel Pipeline again carries out by row one-dimensional Fourier transform to described one-dimensional Fourier transform result by row, obtains the two-dimentional Fourier spectrum of image;
Two-dimensional frequency filtration module (2c), for extracting the fundamental component only comprising elevation information in Fourier spectrum;
Two-dimensional frequency filtering comprises: according to the filter factor in the type computed image of selected wave filter corresponding to each pixel; The filter coefficient of described correspondence is multiplied with the frequency spectrum of image, obtains filtered frequency spectrum;
Two dimension inverse Fourier transform module (2b), is transformed into spatial domain for image spectrum filtering extracted, obtains two-dimentional spatial domain picture;
Two dimension inverse Fourier transform process comprises: adopt N number of Streaming-IFFT to form parallel Pipeline and carry out one dimension inverse Fourier transform by row to image, obtain one dimension inverse Fourier transform result by row; Adopt N number of Streaming-IFFT to form parallel Pipeline again to carry out, by the one dimension inverse Fourier transform of row, obtaining two-dimentional spatial domain picture to one dimension inverse Fourier transform result by row.
5. a kind of dynamic Fourier transform profilometry data processing equipment according to claim 2, it is characterized in that, described intermediate data buffer module (3a), for the data buffering in two-dimensional frequency filtering and address mapping; Data output buffer module (3b), for the data buffering with PERCOM peripheral communication.
6. a kind of dynamic Fourier transform profilometry data processing equipment according to claim 2, it is characterized in that, described phase calculation module (4a), for the spatial domain picture information exported frequency domain filtering unit, be tied to the transformational relation of polar coordinate system according to Cartesian plane coordinate, calculate corresponding original phase hum pattern;
Phase unwrapping module (4b), for launching there is the phase range blocked in described original phase hum pattern, obtains continuous print phase information figure;
Height mapping module (4c), for the distribution of described continuous phase, distributes according to Phase-height mapping relation computed altitude;
Described high-speed communication module (5a), for transmitting with the data of outer CPU and communicating.
7. a dynamic Fourier transform profilometry data processing method, is characterized in that, with the Pipeline that 2N Streaming-FFT composition is parallel, completes the two-dimensional Fourier transform of raster image.
8. a dynamic Fourier transform profilometry data processing method, is characterized in that, with the Pipeline that 2N Streaming-IFFT composition is parallel, completes the two-dimentional inverse Fourier transform of raster image.
9. a dynamic Fourier transform profilometry data processing method, is characterized in that, computer hardware technology is applied to dynamic Fourier transform profilometry Data processing.
10. a kind of dynamic Fourier transform profilometry data processing method according to claim 9, is characterized in that, the hardware development plate for data processing can be any one in FPGA, DSP, MCU.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510399889.XA CN104990514A (en) | 2015-07-09 | 2015-07-09 | Data processing apparatus and method for dynamic Fourier transform profilometry |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510399889.XA CN104990514A (en) | 2015-07-09 | 2015-07-09 | Data processing apparatus and method for dynamic Fourier transform profilometry |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN104990514A true CN104990514A (en) | 2015-10-21 |
Family
ID=54302357
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510399889.XA Pending CN104990514A (en) | 2015-07-09 | 2015-07-09 | Data processing apparatus and method for dynamic Fourier transform profilometry |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104990514A (en) |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105488753A (en) * | 2015-11-27 | 2016-04-13 | 武汉精测电子技术股份有限公司 | Method and device for carrying out two-dimensional Fourier transform and inverse transform on image |
| CN109105357A (en) * | 2018-07-25 | 2019-01-01 | 高雪亚 | Windowsill mosquitocide is sprayed driving mechanism |
| CN110243308A (en) * | 2019-05-31 | 2019-09-17 | 湖北大学 | Real-time phase based on zone location measures technology of profiling |
| CN111879257A (en) * | 2020-07-21 | 2020-11-03 | 南昌航空大学 | High dynamic range real-time three-dimensional measurement method based on Fourier transform profilometry |
| CN113029042A (en) * | 2021-05-25 | 2021-06-25 | 四川大学 | Dynamic measuring device and method for surface morphology of high-temperature molten metal |
| CN113358191A (en) * | 2021-06-04 | 2021-09-07 | 四川大学 | Global flood water level real-time monitoring method based on stripe projection structured light |
| CN114838910A (en) * | 2022-04-24 | 2022-08-02 | 四川大学 | Landslide surge measuring method and simulation measuring device |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20050169497A1 (en) * | 2003-11-19 | 2005-08-04 | Gustafson Ammon E. | Optimized digital watermarking functions for streaming data |
| CN1884993A (en) * | 2005-06-24 | 2006-12-27 | 通用汽车环球科技运作公司 | DSP based rough road detection system |
| CN101174258A (en) * | 2006-11-01 | 2008-05-07 | 冲电气工业株式会社 | Two-dimensional fast fourier transform calculation method and apparatus |
| CN103338368A (en) * | 2013-05-15 | 2013-10-02 | 武汉精测电子技术股份有限公司 | FPGA (field programmable gate array)-based JPEG (joint photographic experts group) parallel decoding device and decoding method |
-
2015
- 2015-07-09 CN CN201510399889.XA patent/CN104990514A/en active Pending
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20050169497A1 (en) * | 2003-11-19 | 2005-08-04 | Gustafson Ammon E. | Optimized digital watermarking functions for streaming data |
| CN1884993A (en) * | 2005-06-24 | 2006-12-27 | 通用汽车环球科技运作公司 | DSP based rough road detection system |
| CN101174258A (en) * | 2006-11-01 | 2008-05-07 | 冲电气工业株式会社 | Two-dimensional fast fourier transform calculation method and apparatus |
| CN103338368A (en) * | 2013-05-15 | 2013-10-02 | 武汉精测电子技术股份有限公司 | FPGA (field programmable gate array)-based JPEG (joint photographic experts group) parallel decoding device and decoding method |
Non-Patent Citations (5)
| Title |
|---|
| 刘艳萍等: "《EDA技术及应用教程》", 30 August 2012, 北京航空航天大学出版社 * |
| 孙丹阳: "高性能数字图像频域滤波系统研究", 《中国优秀硕士学位论文全文数据库信息科技辑》 * |
| 肖焱山: "动态过程中破裂表面的三维重建问题", 《中国优秀博硕士学位论文全文数据库 (硕士)基础科学辑》 * |
| 蔡文贵等: "《CCD技术及应用》", 30 November 1992, 电子工业出版社 * |
| 黄爱民等: "《数字图像处理与分析基础》", 31 August 2005, 中国水利水电出版社 * |
Cited By (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105488753A (en) * | 2015-11-27 | 2016-04-13 | 武汉精测电子技术股份有限公司 | Method and device for carrying out two-dimensional Fourier transform and inverse transform on image |
| CN105488753B (en) * | 2015-11-27 | 2018-12-28 | 武汉精测电子集团股份有限公司 | A kind of pair of image carries out the method and device of two-dimension fourier transform or inverse transformation |
| CN109105357A (en) * | 2018-07-25 | 2019-01-01 | 高雪亚 | Windowsill mosquitocide is sprayed driving mechanism |
| CN110243308A (en) * | 2019-05-31 | 2019-09-17 | 湖北大学 | Real-time phase based on zone location measures technology of profiling |
| CN111879257A (en) * | 2020-07-21 | 2020-11-03 | 南昌航空大学 | High dynamic range real-time three-dimensional measurement method based on Fourier transform profilometry |
| CN113029042A (en) * | 2021-05-25 | 2021-06-25 | 四川大学 | Dynamic measuring device and method for surface morphology of high-temperature molten metal |
| CN113029042B (en) * | 2021-05-25 | 2021-08-03 | 四川大学 | Dynamic measuring device and method for surface morphology of high-temperature molten metal |
| CN113358191A (en) * | 2021-06-04 | 2021-09-07 | 四川大学 | Global flood water level real-time monitoring method based on stripe projection structured light |
| CN114838910A (en) * | 2022-04-24 | 2022-08-02 | 四川大学 | Landslide surge measuring method and simulation measuring device |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104990514A (en) | Data processing apparatus and method for dynamic Fourier transform profilometry | |
| CN110599540B (en) | Real-time three-dimensional human body shape and posture reconstruction method and device under multi-viewpoint camera | |
| CN104048744B (en) | A kind of contactless real-time online vibration measurement method based on image | |
| Chen et al. | Modal identification of simple structures with high-speed video using motion magnification | |
| US8848035B2 (en) | Device for generating three dimensional surface models of moving objects | |
| CN104006765B (en) | Single width carrier frequency interference fringe phase extraction method and detecting device | |
| Gorjup et al. | Frequency domain triangulation for full-field 3D operating-deflection-shape identification | |
| Gorjup et al. | Still-camera multiview Spectral Optical Flow Imaging for 3D operating-deflection-shape identification | |
| WO2018198702A1 (en) | Stress measurement device, stress measurement system, and stress measurement method | |
| CN110108348A (en) | Thin-wall part micro breadth oscillation measurement method and system based on motion amplification optical flow tracking | |
| CN109708842B (en) | Camera lens point spread function measuring method based on single-pixel imaging | |
| Javh et al. | Measuring full-field displacement spectral components using photographs taken with a DSLR camera via an analogue Fourier integral | |
| CN110360954B (en) | Surface shape measuring method and system based on space coordinate calculation | |
| CN109931878A (en) | A kind of building curtain wall seismic deformation monitoring method based on digital speckle label | |
| WO2015197026A1 (en) | Method, apparatus and terminal for acquiring sign data of target object | |
| Li et al. | Application of integrated binocular stereo vision measurement and wireless sensor system in athlete displacement test | |
| CN106885622A (en) | A kind of big visual field multiple spot three-dimensional vibrating measuring method | |
| CN104123725B (en) | A kind of computational methods of single line array camera homography matrix H | |
| CN107990846A (en) | Master based on single frames structure light passively combines depth information acquisition method | |
| CN105579809A (en) | Measurement method, measurement device, measurement program, and computer-readable recording medium recording measurement program | |
| CN102054272A (en) | Displacement vector detection method based on joint fractional Fourier transform and device | |
| CN204154287U (en) | A kind of three-dimensional measuring apparatus of face wrinkle | |
| Qiu et al. | A shape reconstruction and visualization method for a flexible hinged plate using binocular vision | |
| Tognola et al. | A fast and reliable system for 3D surface acquisition and reconstruction | |
| Cosco et al. | Towards phase-based defect detection: A feasibility study in vibrating panels |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| RJ01 | Rejection of invention patent application after publication | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20151021 |