CN104793253A - Airborne electromagnetic data denoising method based on mathematical morphology - Google Patents

Airborne electromagnetic data denoising method based on mathematical morphology Download PDF

Info

Publication number
CN104793253A
CN104793253A CN201510193706.9A CN201510193706A CN104793253A CN 104793253 A CN104793253 A CN 104793253A CN 201510193706 A CN201510193706 A CN 201510193706A CN 104793253 A CN104793253 A CN 104793253A
Authority
CN
China
Prior art keywords
signal
filtering
carry
opening
structural
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
Application number
CN201510193706.9A
Other languages
Chinese (zh)
Other versions
CN104793253B (en
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.)
Jilin University
Original Assignee
Jilin University
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 Jilin University filed Critical Jilin University
Priority to CN201510193706.9A priority Critical patent/CN104793253B/en
Publication of CN104793253A publication Critical patent/CN104793253A/en
Application granted granted Critical
Publication of CN104793253B publication Critical patent/CN104793253B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Abstract

The invention relates to an airborne electromagnetic data denoising method based on mathematical morphology. The airborne electromagnetic data denoising method includes: acquiring airborne electromagnetic detection data through an experiment, and combining triangular structure elements with semicircular structure elements to perform adaptive multiscale composite morphological filtering on the airborne electromagnetic date; determining a length range of the structure elements according to a minimum value and a maximum value of adjacent peak value intervals in original signals, and determining corresponding analysis scale K according to the length range; determining a height range according to a minimum value and a maximum value of signal peak values; adopting structure element sets for composite morphological operation on the original signals, and acquiring an average value as an output result. The adaptive multiscale composite morphological filtering method overcomes the defect of randomness in conventional selection of morphological filtering structure elements; by the airborne electromagnetic data denoising method, structure element type and scale can be selected adaptively according to local features and noise characteristics of signals for filtering airborne electromagnetic signals.

Description

Based on the aviation electromagnetic data de-noising method of mathematical morphology
Technical field
The present invention relates to a kind of data de-noising method of aviation electromagnetic field, especially time domain aviation electromagnetic field, specifically a kind of aviation electromagnetic data de-noising method based on mathematical morphology.
Background technology
Airborne eleectomagnetics is a kind of is carrier with aircraft, carries out the one prospecting detection method of geophysical exploration, is mainly used to quick census metal ore, large area geologic mapping, hydrogeology, the field such as engineering geologic investigation and environmental monitoring.
Mathematical morphology is a subject being based upon in strict mathematical theoretical foundation, now be successfully applied to the engineering practice fields such as image procossing, pattern analysis, pattern-recognition, computer vision, Power Disturbance, mechanical vibration and earthquake detection, and cause extensive attention.The method computing is simple, and its fundamental operation comprises that burn into expands, opening operation and closed operation and the form opening and closing of drawing on this basis and form make and break computing.Signal antinoise method based on mathematical morphology only depends on the local feature of pending signal, the geometric properties of structural element to signal is utilized to mate or revise, retain the principal shape of echo signal, to reach restraint speckle, to extract the object of useful information and reservation details composition simultaneously.
The main bugbear existed in mathematical morphology filter method is choosing of structural element, is the key factor affecting its filter effect on choosing of structural element type.Adopt dissimilar structural element shape facilities different in echo signal can be extracted, choosing of structural element will as much as possible close to the features of shape of pending signal itself, could reach best filter effect as far as possible like this, common structural element type has linear pattern, rectangle, collar plate shape, parabolic type, triangle and other polygons to combine.
Existing shape filtering method only adopts single scale structural element mostly, single scale morphology only selects a fixing structural element to carry out morphological analysis to signal, although this method is simple and be easy to realize, but the quality of its treatment effect relies on relevant priori greatly, and priori is often difficult to obtain accurately and effectively.In addition owing to usually comprising more than a kind of noise type in the signal, and noise in the signal often neither be equally distributed.
Existing shape filtering method all can not carry out complex morphological filtering to dissimilar in the structural element signal under different scale with the noise contribution of intensity.
Summary of the invention
Object of the present invention is just for the deficiencies in the prior art, provides a kind of while restraint speckle, extraction useful information, the aviation electromagnetic data de-noising method based on mathematical morphology of stick signal minutia better.
The object of the invention is to be achieved through the following technical solutions:
Main thought of the present invention is: for the multiple interference component comprised in aviation electromagnetic data, propose the morphologic filtering of self-adapting multi-dimension, have chosen triangular structure element and semicircular structure element, the positive and negative impulsive noise in filtered signal and random noise.
Based on the aviation electromagnetic data de-noising method of mathematical morphology, comprise the following steps:
A, by time domain helicopter electromagnetic detection experiment, through data acquisition hardware circuit, carry out timing equal interval sampling, collect original aerial electromagnetic data;
B, employing triangular structure element combine with semicircular structure element, aviation electromagnetic signal are carried out self-adapting multi-dimension complex morphological filtering process, the positive and negative impulsive noise in filtering echo signal and random noise.
Step B comprises the following steps:
A, first determine to analyze scale size K accordingly with this according to the minimum value at adjacent peak interval in original signal and the length range of maximal value determination structural element;
B, again according to minimum value and the maximal value determination altitude range of signal peak;
C, then utilize little/large length correspondence little/highly determine greatly each structural element in multiscale analysis and form corresponding Multi-structure elements collection;
D, finally each structural element set pair original signal is adopted to carry out complex morphological computing and average as Output rusults.
Described to carry out complex morphological computing to original signal be first carry out erosion operation, dilation operation, opening operation and closed operation respectively to aviation electromagnetic signal, then carry out Opening-closing filter and form make and break filtering respectively:
Erosion operation: ( fΘg ) ( n ) = min m = 0,1 , . . . , M - 1 { f ( n + m ) - g ( m ) } , n = 0,1 , . . . , N + M - 2
Dilation operation: ( f ⊕ g ) ( n ) = max m = 0,1 , . . . , M - 1 { f ( n - m ) + g ( m ) } , n = 0,1 , . . . , N - M
Opening operation:
Closed operation: ( f · g ) ( n ) = [ ( f ⊕ g ) Θg ] ( n )
Wherein, Θ represents erosion operation, represent dilation operation, ο represents opening operation, and represent closed operation, f is original aerial electromagnetic waveforms data, and g is selected morphological structuring elements, m and n is that the discrete sampling of input signal is counted, N > > M;
Carry out Opening-closing filter and form make and break filtering more respectively:
Opening-closing filter:
Form make and break filtering:
Finally, complex morphological wave filter is constructed:
Self-adaptation complex morphological wave filter: z ( n ) = 1 2 Σ i = 1 N [ CFco ( f ( n ) , G i ) + CFoc ( f ( n ) , G i ) ]
Wherein, G=(g 1, g 2..., g i), represent one group of Multi-structure elements collection, N is the number of types of structural element.
Beneficial effect: based on the aviation electromagnetic data de-noising method of mathematical morphology, make it not only be extracted useful information when denoising, also remain the detailed information of aviation electromagnetic signal simultaneously well, so that further to data analysis; The present invention proposes self-adapting multi-dimension complex morphological filtering method, overcome existing shape filtering structural element and choose random problem, according to signal local feature and the adaptive choice structure element type of noise behavior and size, then can carry out filtering to aviation electromagnetic signal.
Accompanying drawing illustrates:
Accompanying drawing 1 is based on the aviation electromagnetic data de-noising method flow diagram of mathematical morphology
Accompanying drawing 2 self-adaptation complex morphological filtering method figure
Accompanying drawing 3 aviation electromagnetic raw-data map
Accompanying drawing 4 aviation electromagnetic original signal and mathematical morphology filter effect contrast figure
Accompanying drawing 5 self-adaptation complex morphological filter effect figure
Embodiment:
Below in conjunction with drawings and Examples, the present invention is described in further detail.
Based on the aviation electromagnetic data de-noising method of mathematical morphology, comprise the following steps:
A, by time domain helicopter electromagnetic detection experiment, through data acquisition hardware circuit, carry out timing equal interval sampling, collect original aerial electromagnetic data;
B, employing triangular structure element combine with semicircular structure element, aviation electromagnetic signal are carried out self-adapting multi-dimension complex morphological filtering process, the positive and negative impulsive noise in filtering echo signal and random noise.
Step B comprises the following steps:
A, first determine to analyze scale size K accordingly with this according to the minimum value at adjacent peak interval in original signal and the length range of maximal value determination structural element;
B, again according to minimum value and the maximal value determination altitude range of signal peak;
C, then utilize little/large length correspondence little/highly determine greatly each structural element in multiscale analysis and form corresponding Multi-structure elements collection;
D, finally each structural element set pair original signal is adopted to carry out complex morphological computing and average as Output rusults.
Described to carry out complex morphological computing to original signal be first carry out erosion operation, dilation operation, opening operation and closed operation respectively to aviation electromagnetic signal, then carry out Opening-closing filter and form make and break filtering respectively:
Erosion operation: ( fΘg ) ( n ) = min m = 0,1 , . . . , M - 1 { f ( n + m ) - g ( m ) } , n = 0,1 , . . . , N + M - 2
Dilation operation: ( f ⊕ g ) ( n ) = max m = 0,1 , . . . , M - 1 { f ( n - m ) + g ( m ) } , n = 0,1 , . . . , N - M
Opening operation:
Closed operation: ( f · g ) ( n ) = [ ( f ⊕ g ) Θg ] ( n )
Wherein, Θ represents erosion operation, represent dilation operation, ο represents opening operation, and represent closed operation, f is original aerial electromagnetic waveforms data, and g is selected morphological structuring elements, m and n is that the discrete sampling of input signal is counted, N > > M;
Carry out Opening-closing filter and form make and break filtering more respectively:
Opening-closing filter:
Form make and break filtering:
Finally, complex morphological wave filter is constructed:
Self-adaptation complex morphological wave filter: z ( n ) = 1 2 Σ i = 1 N [ CFco ( f ( n ) , G i ) + CFoc ( f ( n ) , G i ) ]
Wherein, G=(g 1, g 2..., g i), represent one group of Multi-structure elements collection, N is the number of types of structural element.
Embodiment 1:
(1) original aerial electromagnetic data is obtained by time domain helicopter electromagnetic survey mission test.Be specially: utilize aviation electromagnetic detection system data acquisition hardware circuit to obtain result of detection Wave data, carry out timing equal interval sampling through data acquisition hardware circuit, collect aviation electromagnetic original signal data.
(2) triangular structure element is suitable for the positive and negative impulse noise interference of filtering, and semicircular structure element is suitable for filtering random noise disturbance.Therefore consider to adopt triangle and semicircular structure element to carry out complex morphological filtering to original signal.Detailed process is as follows:
(1) first, respectively erosion operation, dilation operation, opening operation and closed operation are carried out to signal:
Erosion operation: ( fΘg ) ( n ) = min m = 0,1 , . . . , M - 1 { f ( n + m ) - g ( m ) } , n = 0,1 , . . . , N + M - 2
Dilation operation: ( f ⊕ g ) ( n ) = max m = 0,1 , . . . , M - 1 { f ( n - m ) + g ( m ) } , n = 0,1 , . . . , N - M
Opening operation:
Closed operation: ( f · g ) ( n ) = [ ( f ⊕ g ) Θg ] ( n )
(2) Opening-closing filter and form make and break filtering is carried out again respectively:
Opening-closing filter:
Form make and break filtering:
Construct a kind of self-adaptation complex morphological wave filter, thus improve filter effect.
The filtering of self-adaptation complex morphological: z ( n ) = 1 2 Σ i = 1 N [ CFco ( f ( n ) , G i ) + CFoc ( f ( n ) , G i ) ]
Wherein
(3) mathematic(al) representation of triangular structure element and semicircular structure element is as follows respectively:
(1) triangular structure element
g ( i ) = H × [ 1 - | i | L ] , ( i = - L , . . . , 0 , . . . , L )
(2) semicircular structure element
g ( i ) = H × [ 1 - ( i L ) 2 ] , ( i = - L , . . . , 0 , . . . , L )
Wherein, the length of L representative structure element, the height of H representative structure element.
(4) determination of structural element length
If original signal is X={x i| i=1,2 ..., N} (N is the number of data points of original signal), first calculates the local maximum value sequence of original signal, and has all carried out equalization process before the computation.If PE={PE i| i=1,2 ..., N pEbe the local maximum value sequence of original signal, N pEfor the number of local maximum value sequence.If NE={NE i| i=1,2 ..., N nEbe the local minimum value sequence of original signal, N nEfor the number of local minimum value sequence.Definition local maximum is spaced apart local minimum is spaced apart by the feature of the minimum interval of the local maximum obtained and triangle, semicircular structure element, the minimum value K of corresponding morphological structuring elements length dimension can be calculated lminwith maximal value K lmax.
In formula, for the operational symbol that rounds up, for downward rounding operation accords with.
The length dimension sequence K of structural element can be obtained thus lfor:
K l={K lmin,K lmin+1,K lmax-1,K lmax}
(5) determination of structural element height
Due to the amplitude of the height respective signal of structural element, therefore determine the height of structural element according to the minimizing amplitude size of signal local maximum.If maximal value and the minimum value of local maximum sequence are respectively p p maxand p p min, then the minimizing height value of local maximum and local of signal is respectively H pe=p p max-p p minand H ne=p n max-p n min.For making full use of the height local feature information of signal, the height value of the local extremum of definition signal is H e
H e=max(H pe,H ne)
In order to make height value sequence H lwith structural element length dimension K lcorresponding, definable structural element high degree of sequence H lfor:
H l={α·[H e/(K max-K min+1)+(j-1)·H e/(K max-K min+1)]}
j=1,2,...,K max-K min+1
In formula, α is height ratio coefficient, and the present embodiment gets 0.05.
(6) definition of unit structure element
For triangular structure element, the structural element selecting three data points is unit structural element B, and (K-1 dilation operation), during yardstick K=1, KB={0, 1, 0}, during yardstick K=2, KB={0,1, 2, 1,0}, the like, underscore represents the position of initial point.
(7) determination of each mesostructure element
G 1=H l(i)·K l(i)B Ti=1,2,...,K max-K min+1
G 2=H l(i)·K l(i)B Si=1,2,...,K max-K min+1
Wherein, B tand B sbe respectively unit triangle structural element and unit semicircular structure element.

Claims (3)

1., based on an aviation electromagnetic data de-noising method for mathematical morphology, it is characterized in that, comprise the following steps:
A, by time domain helicopter electromagnetic detection experiment, through data acquisition hardware circuit, carry out timing equal interval sampling, collect original aerial electromagnetic data;
B, employing triangular structure element combine with semicircular structure element, aviation electromagnetic signal are carried out self-adapting multi-dimension complex morphological filtering process, the positive and negative impulsive noise in filtering echo signal and random noise.
2., according to the aviation electromagnetic data de-noising method based on mathematical morphology according to claim 1, it is characterized in that, step B comprises the following steps:
A, first determine to analyze scale size K accordingly with this according to the minimum value at adjacent peak interval in original signal and the length range of maximal value determination structural element;
B, again according to minimum value and the maximal value determination altitude range of signal peak;
C, then utilize little/large length correspondence little/highly determine greatly each structural element in multiscale analysis and form corresponding Multi-structure elements collection;
D, finally each structural element set pair original signal is adopted to carry out complex morphological computing and average as Output rusults.
3. according to the aviation electromagnetic data de-noising method based on mathematical morphology according to claim 2, it is characterized in that, described to carry out complex morphological computing to original signal be first carry out erosion operation, dilation operation, opening operation and closed operation respectively to aviation electromagnetic signal, then carry out Opening-closing filter and form make and break filtering respectively:
Erosion operation: ( fΘg ) ( n ) = min m = 0,1 , . . . , M - 1 { f ( n + m ) - g ( m ) } n=0,1,...,N+M-2
Dilation operation: ( f ⊕ g ) ( n ) = min m = 0,1 , . . . , M - 1 { f ( n - m ) + g ( m ) } n=0,1,...,N-M
Opening operation:
Closed operation: ( f · g ) ( n ) = [ ( f ⊕ g ) Θg ] ( n )
Wherein, Θ represents erosion operation, represent dilation operation, ο represents opening operation, and represent closed operation, f is original aerial electromagnetic waveforms data, and g is selected morphological structuring elements, m and n is that the discrete sampling of input signal is counted, N > > M;
Carry out Opening-closing filter and form make and break filtering more respectively:
Opening-closing filter: Foc (f (n))=f ο gg
Form make and break filtering: Fco (f (n))=fg ο g
Finally, complex morphological wave filter is constructed:
Self-adaptation complex morphological wave filter: z ( n ) = 1 2 Σ i = 1 N [ CFco ( f ( n ) , G i ) + CFoc ( f ( n ) , G i ) ]
Wherein, G=(g 1, g 2..., g i), represent one group of Multi-structure elements collection, N is the number of types of structural element.
CN201510193706.9A 2015-04-22 2015-04-22 Aviation electromagnetic data de-noising method based on mathematical morphology Expired - Fee Related CN104793253B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510193706.9A CN104793253B (en) 2015-04-22 2015-04-22 Aviation electromagnetic data de-noising method based on mathematical morphology

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510193706.9A CN104793253B (en) 2015-04-22 2015-04-22 Aviation electromagnetic data de-noising method based on mathematical morphology

Publications (2)

Publication Number Publication Date
CN104793253A true CN104793253A (en) 2015-07-22
CN104793253B CN104793253B (en) 2018-01-26

Family

ID=53558214

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510193706.9A Expired - Fee Related CN104793253B (en) 2015-04-22 2015-04-22 Aviation electromagnetic data de-noising method based on mathematical morphology

Country Status (1)

Country Link
CN (1) CN104793253B (en)

Cited By (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105068131A (en) * 2015-08-03 2015-11-18 中国科学院电子学研究所 Aeromagnetic data leveling method
CN105207645A (en) * 2015-08-25 2015-12-30 苏州汇川技术有限公司 Filtering method based on mathematical morphology and filtering system thereof
CN105629317A (en) * 2016-04-08 2016-06-01 中国矿业大学(北京) Magnetotelluric noise suppressing method based on intersite transfer function
CN105652325A (en) * 2016-01-05 2016-06-08 吉林大学 Exponential fit-adaptive Kalman-based ground-air electromagnetic data de-noising method
CN106154344A (en) * 2016-08-01 2016-11-23 湖南文理学院 A kind of Magnetotelluric signal denoising method based on combined filter
CN106814402A (en) * 2016-12-22 2017-06-09 中石化石油工程技术服务有限公司 Transient electromagnetic signal Prestack Noise Suppression Methods
CN109044365A (en) * 2018-07-02 2018-12-21 苏州大学 The recognition methods of two dimensional motion state based on brain hemoglobin information
CN109143369A (en) * 2017-06-28 2019-01-04 中国石油化工股份有限公司 The method for efficiently filtering out seismic data random disturbances based on shape filtering
CN109712129A (en) * 2018-12-25 2019-05-03 河北工业大学 A kind of arc image processing method based on mathematical morphology
WO2019144486A1 (en) * 2018-01-24 2019-08-01 吉林大学 Method and device for suppressing noise in airborne electromagnetic data
CN110286412A (en) * 2019-07-02 2019-09-27 自然资源部第二海洋研究所 A kind of seabed natural potential observed pattern carrying the double towed bodys of transient electromagnetic and electrical interference bearing calibration
CN110412656A (en) * 2019-07-18 2019-11-05 长江大学 A kind of method and system that Magnetotelluric Data time-domain pressure is made an uproar

Non-Patent Citations (11)

* Cited by examiner, † Cited by third party
Title
CHEN SHI-HAI AND WANG EN-YUAN: "Electromagnetic Radiation Signals of Coal or Rock Denoising Based on Morphological Filter", 《PROCEDIA ENGINEERING》 *
Electromagnetic Radiation Signals of Coal or Rock Denoising Based on Morphological Filter;CHEN Shi-hai and WANG En-yuan;《Procedia Engineering》;20111231;第588-593页 *
刘俊峰,等: "基于多结构元素形态滤波的大地电磁去噪", 《物探与化探》 *
卢志才,等: "多尺度多元素形态滤波自适应降噪研究", 《河北科技大学学报》 *
基于多结构元素形态滤波的大地电磁去噪;刘俊峰,等;《物探与化探》;20140228;第38卷(第1期);第110-112页 *
多尺度多元素形态滤波自适应降噪研究;卢志才,等;《河北科技大学学报》;20110630;第32卷(第3期);第228-232页 *
数学形态滤波与大地电磁噪声压制;汤井田,等;《地球物理学报》;20120531;第55卷(第5期);第1784-1793页 *
数学形态滤波在音频大地电磁去噪中的应用;汤磊,等;《工程地球物理学报》;20130731;第10卷(第4期);第534、535页 *
汤井田,等: "数学形态滤波与大地电磁噪声压制", 《地球物理学报》 *
汤磊,等: "数学形态滤波在音频大地电磁去噪中的应用", 《工程地球物理学报》 *
章立军,等: "自适应多尺度形态学分析及其在轴承故障诊断中的应用", 《北京科技大学学报》 *

Cited By (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105068131A (en) * 2015-08-03 2015-11-18 中国科学院电子学研究所 Aeromagnetic data leveling method
CN105207645A (en) * 2015-08-25 2015-12-30 苏州汇川技术有限公司 Filtering method based on mathematical morphology and filtering system thereof
CN105207645B (en) * 2015-08-25 2018-09-07 苏州汇川技术有限公司 Filtering method based on mathematical morphology and filtering system
CN105652325B (en) * 2016-01-05 2017-09-19 吉林大学 Air-ground electromagnetic data denoising method based on exponential fitting adaptive Kalman
CN105652325A (en) * 2016-01-05 2016-06-08 吉林大学 Exponential fit-adaptive Kalman-based ground-air electromagnetic data de-noising method
CN105629317A (en) * 2016-04-08 2016-06-01 中国矿业大学(北京) Magnetotelluric noise suppressing method based on intersite transfer function
CN106154344A (en) * 2016-08-01 2016-11-23 湖南文理学院 A kind of Magnetotelluric signal denoising method based on combined filter
CN106814402A (en) * 2016-12-22 2017-06-09 中石化石油工程技术服务有限公司 Transient electromagnetic signal Prestack Noise Suppression Methods
CN106814402B (en) * 2016-12-22 2019-08-23 中石化石油工程技术服务有限公司 Transient electromagnetic signal Prestack Noise Suppression Methods
CN109143369A (en) * 2017-06-28 2019-01-04 中国石油化工股份有限公司 The method for efficiently filtering out seismic data random disturbances based on shape filtering
WO2019144486A1 (en) * 2018-01-24 2019-08-01 吉林大学 Method and device for suppressing noise in airborne electromagnetic data
CN109044365A (en) * 2018-07-02 2018-12-21 苏州大学 The recognition methods of two dimensional motion state based on brain hemoglobin information
CN109712129A (en) * 2018-12-25 2019-05-03 河北工业大学 A kind of arc image processing method based on mathematical morphology
CN110286412A (en) * 2019-07-02 2019-09-27 自然资源部第二海洋研究所 A kind of seabed natural potential observed pattern carrying the double towed bodys of transient electromagnetic and electrical interference bearing calibration
CN110412656A (en) * 2019-07-18 2019-11-05 长江大学 A kind of method and system that Magnetotelluric Data time-domain pressure is made an uproar
CN110412656B (en) * 2019-07-18 2021-05-04 长江大学 Magnetotelluric sounding data time domain noise suppression method and system

Also Published As

Publication number Publication date
CN104793253B (en) 2018-01-26

Similar Documents

Publication Publication Date Title
CN104793253A (en) Airborne electromagnetic data denoising method based on mathematical morphology
CN101727662B (en) SAR image nonlocal mean value speckle filtering method
CN106778610A (en) A kind of intra-pulse modulation recognition methods based on time-frequency image feature
CN104360332B (en) Atmospheric phase screen extraction method based on ground-based SAR (synthetic aperture radar) interference
CN103985125B (en) Complicated background SAR image naval ship tail track detection method
CN104462803B (en) A kind of autonomous type underwater robot fault identification method based on small echo approximate entropy
CN105303526A (en) Ship target detection method based on coastline data and spectral analysis
CN101964050A (en) Method for identifying modal parameter based on model order determination and signal noise elimination
Liu et al. LPI radar signal detection based on radial integration of Choi-Williams time-frequency image
CN104331583B (en) A kind of Multifractal Modeling method based on Observed sea clutter
CN106886044A (en) A kind of microseism first break pickup method based on shearing wave Yu Akaike's Information Criterion
CN103473559A (en) SAR image change detection method based on NSCT domain synthetic kernels
CN104316936B (en) A kind of comprehensive DME pulse interference suppression method
CN105699947A (en) SAR image sidelobe inhibition method
CN107632326A (en) Geophysical signal denoising method
CN105652325B (en) Air-ground electromagnetic data denoising method based on exponential fitting adaptive Kalman
Zheng et al. A quantified self-adaptive filtering method: Effective IMFs selection based on CEEMD
CN101650429A (en) Method for extracting chaotic characteristic parameter from low-signal-to-noise ratio underwater target signal
Ji et al. Noise reduction of grounded electrical source airborne transient electromagnetic data using an exponential fitting-adaptive Kalman filter
CN107329137B (en) Ground Penetrating Radar buried target localization method and system based on box spline wavelets frame
CN104076404B (en) Use multichannel to be concerned with and suppress the magnetic anomaly detection method of earth magnetism background noise
CN105960613A (en) System identification device
CN104463325A (en) Noise suppression method for polar ice-penetrating radar original data
Pace et al. Nyquist folding analog-to-information receiver: Autonomous information recovery using quadrature mirror filtering
CN101930605B (en) Synthetic Aperture Radar (SAR) image target extraction method and system based on two-dimensional mixing transform

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
EXSB Decision made by sipo to initiate substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20180126

Termination date: 20190422