CN101256676A - Potential field multi-direction multi-scale edge detection method - Google Patents
Potential field multi-direction multi-scale edge detection method Download PDFInfo
- Publication number
- CN101256676A CN101256676A CNA2008100066766A CN200810006676A CN101256676A CN 101256676 A CN101256676 A CN 101256676A CN A2008100066766 A CNA2008100066766 A CN A2008100066766A CN 200810006676 A CN200810006676 A CN 200810006676A CN 101256676 A CN101256676 A CN 101256676A
- Authority
- CN
- China
- Prior art keywords
- prime
- partiald
- potential field
- wavelet transformation
- theta
- 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 28
- 238000003708 edge detection Methods 0.000 title abstract 3
- 230000009466 transformation Effects 0.000 claims description 22
- 239000000284 extract Substances 0.000 claims description 6
- 238000002592 echocardiography Methods 0.000 claims description 2
- 230000002159 abnormal effect Effects 0.000 abstract 1
- 230000007547 defect Effects 0.000 abstract 1
- 238000004458 analytical method Methods 0.000 description 1
- 238000007405 data analysis Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 230000011664 signaling Effects 0.000 description 1
Images
Abstract
The invention discloses a bit field multi-direction multi-scale edge detection method, which is a technology for quickly inverting a shallow three-dimensional crustal structure based on bit field data. The invention realizes the calculation of the modulus and the argument of the wavelet transform after the coordinate system rotates by adopting the direction wavelet transform, thereby realizing the method that the modulus maximum point of the wavelet transform is adopted as the multi-scale edge of the potential field along the vertical direction connecting line of the gradient vector after the coordinate system rotates by different angles. The invention overcomes the defect that the analytic signal method, the Euler deconvolution method, the horizontal derivative method, the potential field multi-scale edge detection method and the like are insensitive to the information of the directionality, and can quickly obtain the complete and accurate abnormal boundary position information of the geophysical potential field in different depth ranges, thereby generating a more accurate shallow crust three-dimensional structure chart.
Description
Technical field
The present invention relates to a kind of method for checking margin of potential field polydirectionally and multiple dimensions, it is a technology of carrying out the three-dimensional crustal structure fast inversion of superficial part based on the potential field data, relate to fields such as wavelet analysis, Flame Image Process and geophysics, can directly apply to potential field data analysis disposal system, all kinds of geologic examination (areal geology, engineering geology, environmental hazard evaluation, ore prospecting) etc.
Background technology
Anomaly source boundary position detection method based on the potential field data mainly comprises analytic signal method, Euler deconvolution method, horizontal derivative method, potential field multi-scale morphology method etc. at present, its common problem that exists is insensitive to the information of directivity, can not obtain complete, exception boundary position accurately.How accurately, obtain the establish important meaning of different depth anomaly source boundary information fast for superficial part earth's crust three-dimensional structure.
Summary of the invention
The objective of the invention is to above-mentioned deficiency, a kind of method for checking margin of potential field polydirectionally and multiple dimensions is provided, can access the anomaly source boundary information of different directions and different scale (degree of depth), thereby set up superficial part earth's crust three-dimensional structure fast at prior art.
For realizing such purpose, the present invention has designed multi-direction multi-scale morphology method, is applied to geophysics potential field data processing, can obtain more complete, different deep exception boundary positional information accurately.
Method for checking margin of potential field polydirectionally and multiple dimensions of the present invention comprises following concrete steps:
1. potential field data pre-service: handle on the expansion limit of at first carrying out the potential field data, need convert pseudo-gravity to for magnetic field data.
2. according to the potential field data characteristics, select suitable wavelet basis.
3. with plane right-angle coordinate XOY anglec of rotation β, at new coordinate system X ' OY ' down, calculate the mould and the argument of potential field data different scale wavelet transformation, the curve that the modulus maximum point that extracts wavelet transformation connects into along the vertical direction of gradient vector is as the multi-scale edge of potential field.
4. select a plurality of anglec of rotation β to calculate (common 4,8 or 16), contrast the edge of same yardstick different directions, select the edge of its union, carry out cutting edge according to actual bit field data scope then and handle as this yardstick.
5. based on the different scale edge, extract and generate superficial part earth's crust tomograph.
Method for checking margin of potential field polydirectionally and multiple dimensions of the present invention, wherein with plane right-angle coordinate XOY anglec of rotation β, can realize with the direction wavelet transformation at the new following potential field wavelet transformation of coordinate system X ' OY ':
ψ
1(x, y) and ψ
2(x y) is equivalent to two mutually perpendicular direction small echos, and corresponding direction is 0 He
If (x y) is two-dimentional smooth function, ψ to θ
1(x, y) and ψ
2(x y) can use direction small echo ψ
0(x, y) and
Express:
With plane right-angle coordinate XOY anglec of rotation β, at the new following ψ of coordinate system X ' OY '
1(x ', y ') and ψ
2(x ', y ') can use direction small echo ψ
0+ β(x ', y ') and
Express:
At new coordinate system X ' OY ' down, for two-dimentional quadractically integrable function f (x ', y ') ∈ L
2(R
2), can be defined as at the wavelet transformation of yardstick s with two components:
Wavelet transformation can be represented with its mould and argument:
Wherein, Mod (W
sF (x ', y '), Arg (W
sF (x ', y ') represents the mould and the argument of wavelet transformation respectively.
Advantage of the present invention, the present invention has realized the mould of wavelet transformation and the calculating of argument after the coordinate system rotation by adopting the direction wavelet transformation, thereby the modulus maximum point of having realized adopting after the coordinate system rotation different angles wavelet transformation is along the vertical direction line of the gradient vector method as the multi-scale edge of potential field.The present invention has overcome parsing signalling, Euler deconvolution method, horizontal derivative method, potential field multi-scale morphology method etc. to the insensitive shortcoming of the information of directivity, can obtain geophysics potential field, accurately exception boundary positional information complete fast, thereby generate more accurate superficial part earth's crust tomograph in the different depth scope.
Description of drawings
Fig. 1 is multi-direction multi-scale morphology result and multi-scale morphology result's contrast.
Wherein, yardstick s=2
3, Fig. 1 (a) is the multi-scale morphology result, Fig. 1 (b) is multi-direction multi-scale morphology result.
Fig. 2 is area, a middle pasture superficial part earth's crust tomograph, and base map is that area, middle pasture boat magnetic Δ T schemes unusually.
Embodiment
Technical scheme for a better understanding of the present invention is further described embodiments of the present invention below in conjunction with accompanying drawing.
Area, pasture high precision boat magnetic data is carried out the multi-direction multi-scale morphology of potential field in employing the north, Yunnan, comprises following concrete implementation detail:
1) data pre-service: handle on the expansion limit of at first carrying out the potential field data, converts area, middle pasture boat magnetic data to pseudo-gravity.
2) boat magnetic data characteristics in area, pasture in the basis, the quadratic spline function of selecting tight support is as wavelet basis.
3), at new coordinate system X ' OY ' down, calculate potential field data different scale s=2 with plane right-angle coordinate XOY anglec of rotation β
jThe curve that the mould of wavelet transformation and argument, the modulus maximum point that extracts wavelet transformation connect into along the vertical direction of gradient vector is as the multi-scale edge of potential field.
4) select 16 anglec of rotation β to calculate, contrast the edge of same yardstick different directions, select the edge of its union, carry out cutting edge according to boat magnetic data actual range then and handle as this yardstick.Accompanying drawing 1 is multi-direction multi-scale morphology result and multi-scale morphology result's contrast, wherein, and yardstick s=2
3, Fig. 1 (a) is the multi-scale morphology result, and Fig. 1 (b) is multi-direction multi-scale morphology result, and the result shows that the inventive method can obtain complete, exception boundary position accurately, obviously is better than the multi-scale morphology method.
5) based on different scale s=2
jThe edge of (j=2,3,4) extracts and generates superficial part earth's crust tomograph, sees Fig. 2.
Claims (3)
1, method for checking margin of potential field polydirectionally and multiple dimensions is characterized in that:
By adopting the direction wavelet transformation to realize the mould of wavelet transformation and the calculating of argument after the coordinate system rotation, thereby the modulus maximum point of having realized adopting after the coordinate system rotation different angles wavelet transformation is along the vertical direction line of the gradient vector method as the multi-scale edge of potential field.
2, method for checking margin of potential field polydirectionally and multiple dimensions of the present invention comprises following concrete steps:
1. potential field data pre-service: handle on the expansion limit of at first carrying out the potential field data, need convert pseudo-gravity to for magnetic field data.
2. according to the potential field data characteristics, select suitable wavelet basis.
3. with plane right-angle coordinate XOY anglec of rotation β, at new coordinate system X ' OY ' down, calculate the mould and the argument of potential field data different scale wavelet transformation, the curve that the modulus maximum point that extracts wavelet transformation connects into along the vertical direction of gradient vector is as the multi-scale edge of potential field.
4. select a plurality of anglec of rotation β to calculate (common 4,8 or 16), contrast the edge of same yardstick different directions, select the edge of its union as this yardstick.
5. based on the different scale edge, extract and generate superficial part earth's crust tomograph.
3, method for checking margin of potential field polydirectionally and multiple dimensions of the present invention, wherein with plane right-angle coordinate XOY anglec of rotation β, can realize with the direction wavelet transformation at the new following potential field wavelet transformation of coordinate system X ' OY ':
ψ
1(x, y) and ψ
2(x y) is equivalent to two mutually perpendicular direction small echos, and corresponding direction is 0 He
If (x y) is two-dimentional smooth function, ψ to θ
1(x, y) and ψ
2(x y) can use direction small echo ψ
0(x, y) and
Express:
With plane right-angle coordinate XOY anglec of rotation β, at the new following ψ of coordinate system X ' OY '
1(x ', y ') and ψ
2(x ', y ') can use direction small echo ψ
0+ β(x ', y ') and
Express:
At new coordinate system X ' OY ' down, for two-dimentional quadractically integrable function f (x ', y ') ∈ L
2(R
2), can be defined as at the wavelet transformation of yardstick s with two components:
Wavelet transformation can be represented with its mould and argument:
Wherein, Mod (W
sF (x ', y '), Arg (W
sF (x ', y ') represents the mould and the argument of wavelet transformation respectively.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB2008100066766A CN100559400C (en) | 2008-01-31 | 2008-01-31 | Potential field multi-direction multi-scale edge detection method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB2008100066766A CN100559400C (en) | 2008-01-31 | 2008-01-31 | Potential field multi-direction multi-scale edge detection method |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101256676A true CN101256676A (en) | 2008-09-03 |
CN100559400C CN100559400C (en) | 2009-11-11 |
Family
ID=39891476
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CNB2008100066766A Active CN100559400C (en) | 2008-01-31 | 2008-01-31 | Potential field multi-direction multi-scale edge detection method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN100559400C (en) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103020637A (en) * | 2012-11-22 | 2013-04-03 | 北京航空航天大学 | Point cloud data segmentation method on top surface of building based on K-plane algorithm |
CN103020629A (en) * | 2012-11-30 | 2013-04-03 | 中国地质科学院矿产资源研究所 | Method for extracting and identifying bit field data construction information |
CN104658037A (en) * | 2015-03-04 | 2015-05-27 | 中国地质科学院矿产资源研究所 | Automatic extraction method for potential field structural framework |
CN104965232A (en) * | 2015-06-04 | 2015-10-07 | 中国地质科学院矿产资源研究所 | Automatic extraction method of magnetic structure grillwork in low latitude region |
WO2016138874A1 (en) * | 2015-03-04 | 2016-09-09 | 中国地质科学院矿产资源研究所 | Method for automatically extracting tectonic framework of potential field |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0561593B1 (en) * | 1992-03-17 | 1997-07-16 | Sony Corporation | Image compression apparatus |
KR100492517B1 (en) * | 2002-03-07 | 2005-06-03 | 장형수 | Manufacturing apparatus for bast fibers paper using bast fibers extracted from paper mulberry tree, manufacturing methode of the same and bast fibers paper thereby |
JP2007316983A (en) * | 2006-05-26 | 2007-12-06 | Matsushita Electric Ind Co Ltd | Image processor, image processing method and program |
-
2008
- 2008-01-31 CN CNB2008100066766A patent/CN100559400C/en active Active
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103020637A (en) * | 2012-11-22 | 2013-04-03 | 北京航空航天大学 | Point cloud data segmentation method on top surface of building based on K-plane algorithm |
CN103020637B (en) * | 2012-11-22 | 2016-01-27 | 北京航空航天大学 | A kind of buildings end face point cloud data segmentation method based on K-plane algorithm |
CN103020629A (en) * | 2012-11-30 | 2013-04-03 | 中国地质科学院矿产资源研究所 | Method for extracting and identifying bit field data construction information |
CN103020629B (en) * | 2012-11-30 | 2016-01-06 | 中国地质科学院矿产资源研究所 | Method for extracting and identifying bit field data construction information |
CN104658037A (en) * | 2015-03-04 | 2015-05-27 | 中国地质科学院矿产资源研究所 | Automatic extraction method for potential field structural framework |
WO2016138874A1 (en) * | 2015-03-04 | 2016-09-09 | 中国地质科学院矿产资源研究所 | Method for automatically extracting tectonic framework of potential field |
CN104658037B (en) * | 2015-03-04 | 2017-04-05 | 中国地质科学院矿产资源研究所 | Automatic extraction method for potential field structural framework |
RU2664488C1 (en) * | 2015-03-04 | 2018-08-17 | Инститьют Оф Минерал Рисорсиз, Чайниз Акедеми Оф Джиолоджикал Сайенсиз | Method of automatic generation of potential field data structure |
AU2016228027B2 (en) * | 2015-03-04 | 2018-11-22 | Institute of Mineral Resources, Chinese Academy of Geological Sciences | Method for automatically extracting tectonic framework of potential field |
US10884161B2 (en) | 2015-03-04 | 2021-01-05 | Institute of Mineral Resources, Chinese Academy of Geological Sciences | Method for automatically extracting structural framework from potential field data |
CN104965232A (en) * | 2015-06-04 | 2015-10-07 | 中国地质科学院矿产资源研究所 | Automatic extraction method of magnetic structure grillwork in low latitude region |
Also Published As
Publication number | Publication date |
---|---|
CN100559400C (en) | 2009-11-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Salem et al. | Analytic signal approach and its applicability in environmental magnetic investigations | |
Biondi | 3D seismic imaging | |
Boniger et al. | Subsurface utility extraction and characterization: Combining GPR symmetry and polarization attributes | |
Kim et al. | Seismic imaging of the Cocos plate subduction zone system in central Mexico | |
CN100559400C (en) | Potential field multi-direction multi-scale edge detection method | |
Key et al. | Determining the orientation of marine CSEM receivers using orthogonal Procrustes rotation analysis | |
Zhang et al. | NAV-Edge: Edge detection of potential-field sources using normalized anisotropy variance | |
CN105629299A (en) | Travel-time table and angle table acquisition method for angle domain prestack depth migration and imaging method | |
Fedi et al. | Composite continuous wavelet transform of potential fields with different choices of analyzing wavelets | |
Li et al. | Estimating the location of magnetic sources using magnetic gradient tensor data | |
CN102944905A (en) | Gravity and magnetic anomaly handling method based on directional wavelet analysis | |
CN102419457B (en) | Method for determining deep rock structural surface attitude by utilizing television image of single vertical drilling hole | |
Guo et al. | A hybrid positive-and-negative curvature approach for detection of the edges of magnetic anomalies, and its application in the South China Sea | |
Yang et al. | Mantle flow in the vicinity of the eastern edge of the Pacific‐Yakutat slab: Constraints from shear wave splitting analyses | |
Smith et al. | Constraining depth to basement for mineral exploration using microtremor: A demonstration study from remote inland Australia | |
Cheong et al. | Integrated offshore seismic survey using an unmanned wave glider | |
Cevallos | Interpreting the direction of the gravity gradient tensor eigenvectors: Their relation to curvature parameters of the gravity field | |
Zhou | Depth estimation method based on the ratio of gravity and full tensor gradient invariant | |
CN104658037B (en) | Automatic extraction method for potential field structural framework | |
Guo et al. | Global correlation imaging of magnetic total field gradients | |
CN112684498A (en) | Reservoir fracture prediction method and system based on wide-azimuth seismic data | |
Price et al. | 3D CSEM over Frigg-dealing with cultural noise | |
CN104459699B (en) | Water surface and underground target classification and identification method based on target navigation depth | |
Tomar et al. | Sub‐sample time shift and horizontal displacement measurements using phase‐correlation method in time‐lapse seismic | |
Flanagan et al. | Improvements in magnetic depth estimation: application of depth and width extent nomographs to standard depth estimation techniques |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
DD01 | Delivery of document by public notice |
Addressee: Cao Dianhua Document name: payment instructions |