CN110210359A - Tachometric survey optimization method in space filtering velocity measuring technique - Google Patents
Tachometric survey optimization method in space filtering velocity measuring technique Download PDFInfo
- Publication number
- CN110210359A CN110210359A CN201910439522.4A CN201910439522A CN110210359A CN 110210359 A CN110210359 A CN 110210359A CN 201910439522 A CN201910439522 A CN 201910439522A CN 110210359 A CN110210359 A CN 110210359A
- Authority
- CN
- China
- Prior art keywords
- matrix
- subregion
- gray value
- measured
- picture
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 44
- 238000001914 filtration Methods 0.000 title claims abstract description 29
- 238000005457 optimization Methods 0.000 title claims abstract description 15
- 230000003595 spectral effect Effects 0.000 claims abstract description 17
- 238000005070 sampling Methods 0.000 claims abstract description 10
- 239000011159 matrix material Substances 0.000 claims description 38
- 238000004088 simulation Methods 0.000 claims description 10
- 238000012545 processing Methods 0.000 claims description 9
- 230000003321 amplification Effects 0.000 claims description 8
- 238000003199 nucleic acid amplification method Methods 0.000 claims description 8
- 230000005540 biological transmission Effects 0.000 claims description 6
- 238000012360 testing method Methods 0.000 claims description 6
- 230000008569 process Effects 0.000 claims description 5
- 238000006243 chemical reaction Methods 0.000 claims description 4
- 230000002708 enhancing effect Effects 0.000 claims description 4
- 238000003384 imaging method Methods 0.000 claims description 4
- 230000008901 benefit Effects 0.000 claims description 3
- 239000002131 composite material Substances 0.000 claims description 3
- 125000004122 cyclic group Chemical group 0.000 claims description 3
- 230000001419 dependent effect Effects 0.000 claims description 3
- 239000000203 mixture Substances 0.000 claims description 3
- 238000012634 optical imaging Methods 0.000 claims description 3
- 230000000007 visual effect Effects 0.000 claims description 3
- 238000012546 transfer Methods 0.000 claims description 2
- 238000000605 extraction Methods 0.000 abstract description 10
- 238000001228 spectrum Methods 0.000 abstract description 5
- 230000006872 improvement Effects 0.000 abstract description 3
- 238000005259 measurement Methods 0.000 abstract description 3
- 238000000205 computational method Methods 0.000 abstract description 2
- 238000007781 pre-processing Methods 0.000 abstract 1
- 238000005516 engineering process Methods 0.000 description 2
- 230000009471 action Effects 0.000 description 1
- 230000003044 adaptive effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000006835 compression Effects 0.000 description 1
- 238000007906 compression Methods 0.000 description 1
- 230000001186 cumulative effect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005286 illumination Methods 0.000 description 1
- 238000000691 measurement method Methods 0.000 description 1
- 238000004321 preservation Methods 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 238000002834 transmittance Methods 0.000 description 1
- 238000000827 velocimetry Methods 0.000 description 1
Classifications
-
- G06T5/70—
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/11—Region-based segmentation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
- G06F2218/08—Feature extraction
- G06F2218/10—Feature extraction by analysing the shape of a waveform, e.g. extracting parameters relating to peaks
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T10/00—Road transport of goods or passengers
- Y02T10/10—Internal combustion engine [ICE] based vehicles
- Y02T10/40—Engine management systems
Landscapes
- Engineering & Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Image Analysis (AREA)
- Image Processing (AREA)
Abstract
The present invention relates to tachometric survey optimization methods in a kind of space filtering velocity measuring technique to enhance contrast by image preprocessing, and convenient for the extraction of characteristic point, this has good improvement for measuring low-speed motion object in complex environment;By calculus of differences, eliminating part ambient noise bring influences;Suitable grid distance is obtained according to characteristic point size, can be extracted convenient for spectrum peak;Using Welch method, the power Spectral Estimation that available flatness is good, precision is good is convenient for peak extraction.On the basis of existing space filters Computational Method of Velocity Measurement, by the optimization of pretreatment and Power estimation algorithm to sampling picture, more accurate speed is obtained.
Description
Technical field
The present invention relates to a kind of velocity measuring technique, in particular to tachometric survey optimization side in a kind of space filtering velocity measuring technique
Method.
Background technique
Space filtering tests the speed (Spatial Filtering Velocimetry, SFV) technology with its optically and mechanically structure
The optional advantage of simple and stable, light source is more and more paid attention to.When progress space filtering tests the speed, the typical space of use
Filtering device includes rectangular spatial filter, circular space filter and Gauss weighted space filter etc., rectangular filter
It is easier to realize in practical applications, therefore the most commonly used.
Rectangular spatial filter can be simulated by face battle array CMOS camera, the picture of shooting is analyzed and processed,
To obtain rate results.Such as: " non-contact measurement device for measuring and its measurement method of mud-rock flow superficial velocity field distribution " (application is open
Number: CN108344879A) etc., directly ask autocorrelative Fourier transformation to obtain power Spectral Estimation, but obtain in practical application
The pictorial information obtained always contains much noise, and randomness is stronger, and the power Spectral Estimation error obtained using this method is larger, especially
Accurate speed can not be often measured in complex environment.
Summary of the invention
The problem of the present invention be directed to rate accuracy is affected in complex environment proposes a kind of space filtering and tests the speed
Tachometric survey optimization method in technology passes through the pretreatment to sampling picture on the basis of existing space filters Computational Method of Velocity Measurement
And the optimization of Power estimation algorithm, obtain more accurate speed.
The technical solution of the present invention is as follows: tachometric survey optimization method in a kind of space filtering velocity measuring technique, specifically include as
Lower step:
1) under LED white-light illuminating, face battle array CMOS camera acquisition frame rate is set as 1000fps, image model is adopted for gray scale
Collection carries out high speed to region to be measured using face battle array CMOS camera and continuously takes pictures 1 second, and the picture transfer of sampling is entered computer and is waited
Processing, collected picture sample;
2) piecemeal is carried out according to the requirement of user's resolution ratio to acquired image, selects subregion to be measured, it is carried out
Mini-value filtering carries out the binary conversion treatment of image followed by maximum variance between clusters;
3) simulation space-filtering operation is realized using face battle array CMOS camera:
If f (x), f (y) be respectively by camera pixel array it is collected, by ash in step 2) treated picture
Cross, longitudinal gray value function in the matrix of angle value composition;H (x), h (y) are respectively the cross of grating, longitudinal transmission function;
Space filtering theoretical formula is as follows:
s(xr)=∫ f (xr-x)h(x)dx
s(yr)=∫ f (yr-y)h(y)dy
Wherein, xr=vxt+c1,yr=vyt+c2, vxFor the lateral average speed size of subregion to be measured, vyFor sub-district to be measured
Longitudinal average speed size in domain;c1With c2For the dependent constant determined by optical imaging system amplification factor, t is the time for exposure,
R refers to r frame picture;X is the lateral coordinates of the subregion pixel to be measured;Y is the vertical of the subregion pixel to be measured
To coordinate;s(xr) and s (yr) it is that gray value exports;
To the parameter processing of space filtering:
Area to be measured is obtained using picture after enhancing contrast firstly the need of the size of the suitable raster grid spacing of determination
Characteristic point average equivalent area diameter projected is d in domain, then takes raster grid spacing p=2d;
For h (x), h (y), grating simulation is carried out here with the pixel array of face battle array CMOS camera, that is, is obtaining son
In area image, it is sampling interval according to grid distance p, is sampled;
For h (x), the 1st, p+1 ..., np+1 of r frame picture subregion image pixel matrix is classified as by " grating "
Grayscale information after filtering out, and each column extracted are recombinated, it is named as matrix Ar;For h (y), r frame picture subregion figure
As the grayscale information after " grating " filters out of the 1st, p+1 ..., np+1 behavior of picture element matrix, and each row weight that will be extracted
Group is named as matrix Br;Matrix A hereinr、BrThe i.e. visual frame makees the result of product of gray value function and transmission function;
By A obtained in selected framerThe odd column gray value of matrix is added to obtainEven column gray value is added to obtainThen differentiated signal isAll frames are handled by this method, can must export square
Battle array:R=1,2 ..., 1000;
By B obtained in selected framerThe odd-numbered line gray value of matrix is added to obtainEven number line gray value is added
It arrivesThen differentiated signal isAll frames are handled by this method, can must export square
Battle array:R=1,2 ..., 1000;
The two-dimensional matrix that gray value information after simulation grating filters out is constituted;
4) Matrix C that Len got in step 3) is 1000 is divided into K sections, every segment length is N and overlaps each other, to every section
Signal carries out windowing process using Hamming window;I-th section of signal may be expressed as:
Ci(m)=C (m+iE), m=0 ..., N-1;I=0 ..., K-1
M is that number is arranged in signal, and E is the length of overlapping;
Every section of cyclic graph power Spectral Estimation is
W (m) is Hamming window function, and U is the energy of Hamming window function, and f is time-domain signal Ci(m) frequency;
Then it is averaged, finally obtains the power Spectral Estimation of entire C, i.e.,
Frequency from obtaining peak amplitude in entire power Spectral Estimation is denoted as f as maximum frequencyxmax;
Same treatment similarly is made to D, obtains maximum frequency and is denoted as fymax;
5) it is built into scale bar in the plane coplanar with object under test, imaging amplification factor is obtained after comparison, is denoted as M;
6) it seeks the velocity component in all directions: utilizing following formula:
The lateral average speed v of subregion to be measured can be obtainedxSize and longitudinal average speed vySize;
7) velocity composite:
The velocity vector of the available subregion to be measured.
The beneficial effects of the present invention are: tachometric survey optimization method in space filtering velocity measuring technique of the present invention passes through figure
As pretreatment, enhance contrast, convenient for the extraction of characteristic point, this has very well for measuring low-speed motion object in complex environment
Improvement;By calculus of differences, eliminating part ambient noise bring influences;Using Welch method, available flatness is good,
The good power Spectral Estimation of precision is convenient for peak extraction;Suitable grid distance is obtained according to characteristic point size, power spectrum can be convenient for
Peak extraction.
Detailed description of the invention
Fig. 1 is acquired original image;
Fig. 2 is that the present invention carries out pretreated image;
Fig. 3 is the spectrogram directly obtained to the target area of original image using traditional auto-correlation algorithm;
Fig. 4 is that the present invention obtains spectrogram using the algorithm of optimization to pretreated target area.
Specific embodiment
Tachometric survey optimization method in a kind of space filtering velocity measuring technique, comprising the following steps:
Step 1 sets face battle array CMOS camera acquisition frame rate as 1000fps, schemes under the LED white-light illuminating of certain brightness
Picture mode is sampled grey, carries out high speed to region to be measured using face battle array CMOS camera and continuously takes pictures 1 second, by the picture of sampling
It is to be processed to be conveyed into computer etc., collected picture sample;
Step 2 can carry out piecemeal according to the requirement of user's resolution ratio to acquired image, select subregion to be measured,
Mini-value filtering is carried out to it, the binary conversion treatment of image is carried out followed by maximum variance between clusters;
Mini-value filtering is a kind of image processing means, first sequence surrounding pixel sum of the grayscale values center pixel gray value,
Then by center pixel gray value compared with minimum gradation value, if smaller than minimum gradation value, minimum gradation value is substituted for
Heart grey scale pixel value, can eliminate the even bring of uneven illumination using this method influences;
Maximum variance between clusters are a kind of methods that adaptive threshold value determines.It is the gamma characteristic by image, will be schemed
As being divided into two parts of background and target.Inter-class variance between background and target is bigger, illustrates to constitute the two-part of image
Difference is bigger, when partial target mistake is divided into background or part background mistake is divided into target and all two parts difference can be caused to become smaller.Therefore,
The maximum segmentation of inter-class variance is set to mean misclassification probability minimum.It can help to find suitable threshold value using this method, with
Just binary conversion treatment, the picture superposition after binaryzation, convenient for the extraction of characteristic point are carried out;
Step 3 realizes simulation space-filtering operation using face battle array CMOS camera:
If f (x), f (y) be respectively by camera pixel array it is collected, by ash in step 2 treated picture
Cross, longitudinal gray value function in the matrix of angle value composition;H (x), h (y) are respectively the cross of grating, longitudinal transmission function.Space
Filtering theory formula is as follows:
s(xr)=∫ f (xr-x)h(x)dx
s(yr)=∫ f (yr-y)h(y)dy
Wherein, xr=vxt+c1,yr=vyt+c2, vxFor the lateral average speed size of subregion to be measured, vyFor sub-district to be measured
Longitudinal average speed size in domain;c1With c2For the dependent constant determined by optical imaging system amplification factor, t is the time for exposure,
R refers to r frame picture;X is the lateral coordinates of the subregion pixel to be measured;Y is the vertical of the subregion pixel to be measured
To coordinate;s(xr) and s (yr) it is that gray value exports.
To the parameter processing of space filtering:
Area to be measured is obtained using picture after enhancing contrast firstly the need of the size of the suitable raster grid spacing of determination
All characteristic point average equivalent area diameter projecteds are d in domain, then take raster grid spacing p=2d;
For h (x), h (y), grating simulation is carried out here with the pixel array of face battle array CMOS camera, that is, is obtaining son
In area image, it is sampling interval according to grid distance p, is sampled.
For h (x), the 1st, p+1 ..., np+1 of r frame picture subregion image pixel matrix is classified as by " grating "
Grayscale information after filtering out, and each column extracted are recombinated, it is named as matrix Ar;For h (y), r frame picture subregion figure
As the grayscale information after " grating " filters out of the 1st, p+1 ..., np+1 behavior of picture element matrix, and each row weight that will be extracted
Group is named as matrix Br.Matrix A hereinr、BrThe i.e. visual frame makees gray value function and transmission function (i.e. rectangular space filtering
Device) result of product.Due to the introducing of additional light source, need to carry out difference processing to signal, to eliminate part ambient noise.
By A obtained in selected framerThe odd column gray value of matrix is added to obtainEven column gray value is added to obtainThen differentiated signal isAll frames are handled by this method, can must export square
Battle array:R=1,2 ..., 1000.
By B obtained in selected framerThe odd-numbered line gray value of matrix is added to obtainEven number line gray value is added to obtainThen differentiated signal isAll frames are handled by this method, can must export square
Battle array:R=1,2 ..., 1000.
The two-dimensional matrix that gray value information after simulation grating filters out is constituted can be by pixel by above-mentioned difference method
Matrix compression is a gray value, and which greatly enhances Computing efficiency.
Step 4, to output matrix power Power estimation:
The Matrix C that Len got in step 3 is 1000 is divided into K sections, every segment length is N and overlaps each other, and is believed every section
Number utilize Hamming window carry out windowing process.I-th section of signal may be expressed as: Ci(m)=C (m+iE), m=0 ..., N-1;i
=0 ..., K-1
M is that number is arranged in signal, and E is the length of overlapping.E=N/2 indicates that half overlaps each other, and E=N expression does not have
Overlapping.In this way, every section of cyclic graph power Spectral Estimation is
W (m) is Hamming window function, and U is the energy of Hamming window function, and f is time-domain signal Ci(m) frequency:
Then it is averaged, finally obtains the power Spectral Estimation of entire C, i.e.,
Frequency from obtaining peak amplitude in entire power Spectral Estimation is denoted as f as maximum frequencyxmax.Similarly D is made
Same treatment obtains maximum frequency and is denoted as fymax。
Step 5 utilizes the pixel dimension parameter of camera, comparison by measuring pixel size shared by feature object in picture
Feature object actual size, can calculate imaging amplification factor, be denoted as M;
Step 6 seeks the velocity component in all directions:
Utilize formula:
The lateral average speed v of subregion to be measured can be obtainedxSize and longitudinal average speed vySize;
Step 7, velocity composite:
It utilizes
The velocity vector of the available subregion to be measured.
In order to be easy to understand the technical means, the creative features, the aims and the efficiencies achieved by the present invention, lead to below
It crosses embodiment and tachometric survey optimization algorithm provided by the invention is specifically addressed in conjunction with attached drawing.
The picture of acquisition is utilized first and pre-processes enhancing contrast, in picture after treatment, is obtained in region to be measured
Characteristic point average equivalent area diameter projected is d, sets grid distance p=2d.
The size of collected 1000 frame image is 1547*768, and wherein useful area is 1300*600, is utilized
MATLAB divides it, and obtaining size is the subregion for being 135*50.
For subregion to be measured, simulation shading treatment is carried out using MATLAB according to the size of grid distance p, by transmittance section
Each row data divided carry out cumulative summation and obtain the matrix of 1*135.Even number line, odd-numbered line are separately summed, then after will add up
Even number line and odd-numbered line after being added subtract each other the gray scale difference score value for just having obtained the region in this frame image.According to phase Tongfang
Residue frame is made same treatment by method.This method is for acquiring longitudinal velocity vy, seek vxWhen, column are operated.This link is available
Gray value sequence:
Utilize the tachometric survey optimization algorithm of above-mentioned any one, which comprises the following steps:
Step 1 sets face battle array CMOS camera acquisition frame rate as 1000fps, and image model is sampled grey, continuous acquisition
Image;
Step 2 acquires image 1 second, and collected photo is 1000 frames at this time, and preservation is to be processed to computer etc., is adopted
The picture sample collected is as shown in Figure 1;
Step 3, in the plane coplanar with object under test, merging is decorated with the card of 1cm straight line as scale bar scale,
A picture is acquired using face battle array CMOS camera, for measuring imaging amplification factor M;
The size of step 4, collected 1000 frame image is 1547*768, and wherein useful area is 1300*600, benefit
Directly the 1st frame image can be pre-processed from tape function ordfilt2 () and graythresh () with MATLAB, handled
Picture afterwards is as shown in Figure 2.Then it is divided, is divided into the subregion of 135*50.Shown under sub-district domain matrix after division;
Step 5 is sampling interval according to grid distance p, is sampled in obtaining sub-district area image.Sub-district area image
The 1st, p+1 of picture element matrix ..., np+1 is classified as the grayscale information after " grating " filters out, and each column extracted are heavy
Group is named as matrix A1, by A1The odd column gray value of matrix is added to obtainEven column gray value is added to obtain
Then differentiated signal isSimilar processing is carried out to remaining frame by this method, output matrix can be obtained:Similarly by row sampling, can obtain
Step 6, using MATLAB software, carries out power Spectral Estimation, parameter setting using Welch method to it for C, D
Are as follows: window type is Hamming window, and length is 100 points, and FFT points take obtain power convenient for Fast Fourier Transform (FFT) processing at 1024 points
Power estimation figure, (as shown in figure 4, abscissa is frequency f, compared with tradition spectrum estimation method (Fig. 3 institute in obtained power Spectral Estimation figure
Showing) peak value is more obvious, easily facilitate extraction);
Step 7 obtains the frequency f at peak amplitudexmaxAnd fymax, utilize formula
Obtain the velocity component in both direction;
Step 8 is summed using vector
It can obtain the velocity vector of the subregion;
In the present embodiment tachometric survey optimization algorithm process the following steps are included:
Step S1, setting face battle array CMOS camera parameter;
Step S2, Image Acquisition;
Step S3, merging scale bar obtain amplification factor;
Step S4, picture pretreatment and sub-zone dividing;
Step S5, subregion gray value difference processing;
Step S6 carries out Welch method power Spectral Estimation to result is exported after difference;
Step S7, analyzes power spectrum chart, obtains frequency at peak amplitude, calculates the speed point on different directions
Amount;
Step S8 carries out vector summation to the speed in both direction, obtains velocity vector;
The action and effect of embodiment:
By pretreatment, enhance contrast, convenient for the extraction of characteristic point, this in complex environment for measuring low-speed motion
Object has good improvement;By calculus of differences, eliminating part ambient noise bring influences;It is available using Welch method
The power Spectral Estimation that flatness is good, precision is good is convenient for peak extraction;Suitable grid distance is obtained according to characteristic point size, it can
Convenient for spectrum peak extraction.
Claims (1)
1. tachometric survey optimization method in a kind of space filtering velocity measuring technique, which is characterized in that specifically comprise the following steps:
1) under LED white-light illuminating, face battle array CMOS camera acquisition frame rate is set as 1000fps, image model is sampled grey, benefit
Region to be measured continuously take pictures 1 second at a high speed with face battle array CMOS camera, it is to be processed that the picture transfer of sampling is entered computer etc.,
Collected picture sample;
2) piecemeal is carried out according to the requirement of user's resolution ratio to acquired image, selects subregion to be measured, minimum is carried out to it
Value filtering carries out the binary conversion treatment of image followed by maximum variance between clusters;
3) simulation space-filtering operation is realized using face battle array CMOS camera:
If f (x), f (y) be respectively by camera pixel array it is collected, by gray value in step 2) treated picture
Cross, longitudinal gray value function in the matrix of composition;H (x), h (y) are respectively the cross of grating, longitudinal transmission function;
Space filtering theoretical formula is as follows:
s(xr)=∫ f (xr-x)h(x)dx
s(yr)=∫ f (yr-y)h(y)dy
Wherein, xr=vxt+c1,yr=vyt+c2, vxFor the lateral average speed size of subregion to be measured, vyFor subregion to be measured
Longitudinal average speed size;c1With c2For the dependent constant determined by optical imaging system amplification factor, t is the time for exposure, and r refers to
Generation r frame picture;X is the lateral coordinates of the subregion pixel to be measured;Y is that the longitudinal of the subregion pixel to be measured sits
Mark;s(xr) and s (yr) it is that gray value exports;
To the parameter processing of space filtering:
It is obtained in region to be measured firstly the need of the size of the suitable raster grid spacing of determination using picture after enhancing contrast
Characteristic point average equivalent area diameter projected is d, then takes raster grid spacing p=2d;
For h (x), h (y), grating simulation is carried out here with the pixel array of face battle array CMOS camera, that is, is obtaining subregion
In image, it is sampling interval according to grid distance p, is sampled;
For h (x), the 1st, p+1 ..., np+1 of r frame picture subregion image pixel matrix is classified as to be filtered out by " grating "
Grayscale information afterwards, and each column extracted are recombinated, it is named as matrix Ar;For h (y), r frame picture subregion image slices
Grayscale information of 1st, p+1 ..., the np+1 behavior of prime matrix after " grating " filters out, and each row extracted is recombinated,
It is named as matrix Br;Matrix A hereinr、BrThe i.e. visual frame makees the result of product of gray value function and transmission function;
By A obtained in selected framerThe odd column gray value of matrix is added to obtainEven column gray value is added to obtainThen differentiated signal isAll frames are handled by this method, can must export square
Battle array:
By B obtained in selected framerThe odd-numbered line gray value of matrix is added to obtainEven number line gray value is added to obtainThen differentiated signal isAll frames are handled by this method, can must export square
Battle array:Gray value information after simulation grating filters out is constituted
Two-dimensional matrix;
4) Matrix C that Len got in step 3) is 1000 is divided into K sections, every segment length is N and overlaps each other, to every segment signal
Windowing process is carried out using Hamming window;I-th section of signal may be expressed as:
Ci(m)=C (m+iE), m=0 ..., N-1;I=0 ..., K-1
M is that number is arranged in signal, and E is the length of overlapping;
Every section of cyclic graph power Spectral Estimation is
W (m) is Hamming window function, and U is the energy of Hamming window function, and f is time-domain signal Ci(m) frequency;
Then it is averaged, finally obtains the power Spectral Estimation of entire C, i.e.,
Frequency from obtaining peak amplitude in entire power Spectral Estimation is denoted as f as maximum frequencyxmax;
Same treatment similarly is made to D, obtains maximum frequency and is denoted as fymax;
5) it is built into scale bar in the plane coplanar with object under test, imaging amplification factor is obtained after comparison, is denoted as M;
6) it seeks the velocity component in all directions: utilizing following formula:
The lateral average speed v of subregion to be measured can be obtainedxSize and longitudinal average speed vySize;
7) velocity composite:
The velocity vector of the available subregion to be measured.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910439522.4A CN110210359B (en) | 2019-05-24 | 2019-05-24 | Speed measurement optimization method in spatial filtering speed measurement technology |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910439522.4A CN110210359B (en) | 2019-05-24 | 2019-05-24 | Speed measurement optimization method in spatial filtering speed measurement technology |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110210359A true CN110210359A (en) | 2019-09-06 |
CN110210359B CN110210359B (en) | 2023-08-22 |
Family
ID=67788589
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910439522.4A Active CN110210359B (en) | 2019-05-24 | 2019-05-24 | Speed measurement optimization method in spatial filtering speed measurement technology |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110210359B (en) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110677582A (en) * | 2019-09-27 | 2020-01-10 | 中国科学院长春光学精密机械与物理研究所 | Filtering speed measurement focus detection method and system and terminal equipment |
CN111931127A (en) * | 2020-06-29 | 2020-11-13 | 电子科技大学 | Efficient implementation method for power spectral density estimation |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20130222786A1 (en) * | 2010-08-09 | 2013-08-29 | Steen Hanson | Vector velocimeter |
CN108279317A (en) * | 2018-01-23 | 2018-07-13 | 江西省智成测控技术研究所有限责任公司 | A kind of space filtering tachogenerator device and the method for improving rate accuracy |
CN109714513A (en) * | 2019-02-15 | 2019-05-03 | 江西省智成测控技术研究所有限责任公司 | Inhibit the method for velocity calculated noise in a kind of optics speed and mileage measuring instrument |
-
2019
- 2019-05-24 CN CN201910439522.4A patent/CN110210359B/en active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20130222786A1 (en) * | 2010-08-09 | 2013-08-29 | Steen Hanson | Vector velocimeter |
CN108279317A (en) * | 2018-01-23 | 2018-07-13 | 江西省智成测控技术研究所有限责任公司 | A kind of space filtering tachogenerator device and the method for improving rate accuracy |
CN109714513A (en) * | 2019-02-15 | 2019-05-03 | 江西省智成测控技术研究所有限责任公司 | Inhibit the method for velocity calculated noise in a kind of optics speed and mileage measuring instrument |
Non-Patent Citations (2)
Title |
---|
OLESEN,AS,ET AL: "Spatial filtering velocimetry for real-time measurements of speckle dynamics due to out of plane motion", APPLIED OPTICS, vol. 55, no. 55, pages 3858 - 3865 * |
侯鹏;杨晖;李然;林世昊;华云松: "复杂颗粒流速度场分布的空间滤波测量", 光学精密工程, no. 011, pages 2632 - 2638 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110677582A (en) * | 2019-09-27 | 2020-01-10 | 中国科学院长春光学精密机械与物理研究所 | Filtering speed measurement focus detection method and system and terminal equipment |
CN110677582B (en) * | 2019-09-27 | 2020-11-24 | 中国科学院长春光学精密机械与物理研究所 | Filtering speed measurement focus detection method and system and terminal equipment |
CN111931127A (en) * | 2020-06-29 | 2020-11-13 | 电子科技大学 | Efficient implementation method for power spectral density estimation |
CN111931127B (en) * | 2020-06-29 | 2023-03-14 | 电子科技大学 | Method for realizing power spectral density estimation |
Also Published As
Publication number | Publication date |
---|---|
CN110210359B (en) | 2023-08-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110120020A (en) | A kind of SAR image denoising method based on multiple dimensioned empty residual error attention network | |
CN103702015B (en) | Exposure control method for human face image acquisition system under near-infrared condition | |
CN108230367A (en) | A kind of quick method for tracking and positioning to set objective in greyscale video | |
CN101567977A (en) | Flicker detection method and device thereof | |
CN103400129A (en) | Target tracking method based on frequency domain saliency | |
CN108225537A (en) | A kind of contactless small items vibration measurement method based on high-speed photography | |
CN110210359A (en) | Tachometric survey optimization method in space filtering velocity measuring technique | |
CN111062978B (en) | Texture recognition method for spatio-temporal image flow measurement based on frequency domain filtering technology | |
CN111043988B (en) | Single stripe projection measurement method based on graphics and deep learning | |
CN103295241A (en) | Frequency domain saliency target detection method based on Gabor wavelets | |
CN1786658A (en) | Method and apparatus for measuring profile of object by double wavelength structural light | |
CN105225216A (en) | Based on the Iris preprocessing algorithm of space apart from circle mark rim detection | |
CN104835177A (en) | Star point segmentation method under interference of moonlight and FPGA device of star point segmentation method | |
CN110097634A (en) | A kind of terrible imaging method of the three-dimensional of self-adapting multi-dimension | |
CN108871197A (en) | Displacement field measurement method, device, equipment and storage medium for material surface | |
CN113139904A (en) | Image blind super-resolution method and system | |
CN105930793A (en) | Human body detection method based on SAE characteristic visual learning | |
CN110111347A (en) | Logos extracting method, device and storage medium | |
US10075702B2 (en) | Electronic device with an upscaling processor and associated methods | |
CN108507476A (en) | Displacement field measurement method, device, equipment and storage medium for material surface | |
Yang et al. | A high-efficiency acquisition method of LED multispectral images using Gray code based square wave frequency division modulation | |
CN116183226A (en) | Bearing test bed vibration displacement measurement and modal analysis algorithm based on phase | |
Yang et al. | Optimized lighting method of applying shaped-function signal for increasing the dynamic range of LED-multispectral imaging system | |
CN105631897B (en) | Based on the film nuclear magnetic resonance image sequence motion method of estimation for singly drilling signal characteristic distance and cross-correlation transformation optical flow algorithm | |
CN107038706A (en) | Infrared image confidence level estimation device and method based on adaptive mesh |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |