共查询到20条相似文献,搜索用时 31 毫秒
1.
2.
Factorial Kriging (FK) is a data- dependent spatial filtering method that can be used to remove both independent and correlated
noise on geological images as well as to enhance lineaments for subsequent geological interpretation. The spatial variability
of signal, noise, and lineaments, characterized by a variogram model, have been used explicitly in calculating FK filter coefficients
that are equivalent to the kriging weighting coefficients. This is in contrast to the conventional spatial filtering method
by predefined, data-independent filters, such as Gaussian and Sobel filters. The geostatistically optimal FK filter coefficients,
however, do not guarantee an optimal filtering effect, if filter geometry (size and shape) are not properly selected. The
selection of filter geometry has been investigated by examining the sensitivity of the FK filter coefficients to changes in
filter size as well as variogram characteristics, such as nugget effect, type, range of influence, and anisotropy. The efficiency
of data-dependent FK filtering relative to data-independent spatial filters has been evaluated through simulated stochastic
images by two examples. In the first example, both FK and data-independent filters are used to remove white noise in simulated
images. FK filtering results in a less blurring effect than the data-independent fillers, even for a filter size as large
as 9 × 9. In the second example, FK and data-independent filters are compared relative to the extraction of lineaments and
components showing anisotropic variability. It was determined that square windows of the filter mask are effective only for
removing Isotropie components or white noise. A nonsquare windows must be used if anisotropic components are to be filtered
out. FK filtering for lineament enhancement is shown to be resistant to image noise, whereas data-independent filters are
sensitive to the presence of noise. We also have applied the FK filtering to the GLORIA side-scan sonar image from the Gulf
of Mexico, illustrating that FK is superior to the data-independent filters in removing noise and enhancing lineaments. The
case study also demonstrate that variogram analysis and FK filtering can be used for large images if a spectral analysis and
optimal filter design in the frequency domain is prohibitive because of a large memory requirement. 相似文献
3.
Denis Marcotte 《Mathematical Geology》1995,27(6):749-762
Conditional simulation with data subject to measurement error has received little attention in the geostatistical literature. The treatment of measurement error in simulation must be different from its treatment in estimation. Two approaches are examined: pre- and post-simulation filtering of data measurement error. The pre-simulation filtering is shown to be inefficient. The post-simulation filtering performs best. It is done by factorial kriging and a modified version of factorial kriging which ensures predetermined theoretical variance for the filtered data. It also is shown that the theoretical variogram of the filtered data reproduces the underlying variogram (i.e., without noise) almost perfectly. A simulation with a high level of correlated noise is used for validation and comparison. The post-simulation filtered values show an experimental variogram in agreement with the previously identified underlying variogram. Moreover, the filtered image compares well with the true image. The theoretical variogram corresponding to the post-simulation filter can be computed beforehand. Thus, the size of the simulation grid and of the filter neighborhood can be adjusted to ensure good reproduction of the underlying variogram. 相似文献
4.
On Visualization for Assessing Kriging Outcomes 总被引:7,自引:0,他引:7
James R. Carr 《Mathematical Geology》2002,34(4):421-433
Extant opinion about kriging is that all weights should be positive. Visualizations rendered by converting kriged grids to digital images are presented to show that negative weights may be beneficial to some spatial problems. In particular, variogram models with zero-valued nuggets, already well known to minimize smoothing through kriging, result in a visual resolution substantially superior to that from kriging with a variogram model having a nonzero nugget value in application to satellite acquired data. Negative weights are more likely when using variogram models with zero-valued nuggets, but resultant visualizations often show a smoother transition between extreme data values. This is true even when a variogram model having a nugget value of zero is not optimum with respect to mean square error, as is demonstrated using a nitrate data set. An analogy to digital image processing is used to suggest that the influence of negative weights in kriging is similar to a high-boost kernel. 相似文献
5.
Kriging Regionalized Positive Variables Revisited: Sample Space and Scale Considerations 总被引:1,自引:0,他引:1
Frequently, regionalized positive variables are treated by preliminarily applying a logarithm, and kriging estimates are back-transformed
using classical formulae for the expectation of a lognormal random variable. This practice has several problems (lack of robustness,
non-optimal confidence intervals, etc.), particularly when estimating block averages. Therefore, many practitioners take exponentials
of the kriging estimates, although the final estimations are deemed as non-optimal. Another approach arises when the nature
of the sample space and the scale of the data are considered. Since these concepts can be suitably captured by an Euclidean
space structure, we may define an optimal kriging estimator for positive variables, with all properties analogous to those
of linear geostatistical techniques, even for the estimation of block averages. In this particular case, no assumption on
preservation of lognormality is needed. From a practical point of view, the proposed method coincides with the median estimator
and offers theoretical ground to this extended practice. Thus, existing software and routines remain fully applicable. 相似文献
6.
针对数字全景钻孔摄像系统获取的实测图像,提出了一种基于颜色特征的数字式钻孔图像溶隙结构识别方法。利用岩层中的典型结构,如土质层、溶隙在颜色上与普通岩石具有较大差异性的特点,首先建立了一个自适应HSV颜色空间溶隙结构检测模型,利用该模型获取溶隙结构的二值化图像;对该二值化图像进行滤波处理;然后从处理后的二值化图像进行分区像素密度统计来确定土质层或溶隙区域的深度、面积及方位角等信息,从而实现数字式钻孔图像中溶隙结构的自动识别。通过对大量数字式钻孔图像进行试验并与对应的钻孔雷达图像进行结果对比表明,其方法能对全孔图像的溶隙和土质层进行快速、准确地自动化检测与定位,为钻孔图像岩体结构的自动识别与工程应用提供了一种新的可靠方法。 相似文献
7.
Antonio M. Montufo 《Geoarchaeology》1997,12(1):71-85
The aim of this research was to assess the potential use of satellite imagery and digital image processing for the detection and surveying of ancient land-use patterns. The rural land-use patterns in the southeastern region of Mallorca (Spain) were investigated in order to locate traces of the Roman field-division system. A satellite scene was analyzed using ground control for field observations to determine the rural patterns observed in the contemporary landscape and to identify the possible remains of centuriated systems. Digital image processing involved image rectification, georeferencing, and spatial enhancement using high-pass and directional filtering. While no incontrovertible evidence of centuriation was located, satellite data proved to be of use in surveying medium-scale rural patterns. © 1997 John Wiley & Sons, Inc. 相似文献
8.
Estimation or simulation? That is the question 总被引:1,自引:0,他引:1
Jorge Kazuo Yamamoto 《Computational Geosciences》2008,12(4):573-591
The issue of smoothing in kriging has been addressed either by estimation or simulation. The solution via estimation calls
for postprocessing kriging estimates in order to correct the smoothing effect. Stochastic simulation provides equiprobable
images presenting no smoothing and reproducing the covariance model. Consequently, these images reproduce both the sample
histogram and the sample semivariogram. However, there is still a problem, which is the lack of local accuracy of simulated
images. In this paper, a postprocessing algorithm for correcting the smoothing effect of ordinary kriging estimates is compared
with sequential Gaussian simulation realizations. Based on samples drawn from exhaustive data sets, the postprocessing algorithm
is shown to be superior to any individual simulation realization yet, at the expense of providing one deterministic estimate
of the random function. 相似文献
9.
Optimally locating additional drill holes in three dimensions using grade and simulated annealing 总被引:2,自引:0,他引:2
Saeed Soltani Mohammadi Ardeshir Hezarkhani A. Erhan Tercan 《Journal of the Geological Society of India》2012,80(5):700-706
Average kriging variance is a standard tool used in optimization of the location of additional drill holes. However, this tool cannot distinguish between areas with different priorities. This limitation could be eliminated by using weighted average kriging variance. This paper extends the problem of optimal location to three dimensional cases, use grade as a weight and search optimum locations by simulated annealing. Weighted average kriging variance is used as objective function. The method is applied to a copper deposit. Results have shown that weighting of the estimation variance with ??grade?? is effective only when the difference among the grades estimated for different blocks is considerable. 相似文献
10.
Image analysis of three-dimensional microtomographic image data has become an integral component of pore scale investigations
of multiphase flow through porous media. This study focuses on the validation of image analysis algorithms for identifying
phases and estimating porosity, saturation, solid surface area, and interfacial area between fluid phases from gray-scale
X-ray microtomographic image data. The data used in this study consisted of (1) a two-phase high precision bead pack from
which porosity and solid surface area estimates were obtained and (2) three-phase cylindrical capillary tubes of three different
radii, each containing an air–water interface, from which interfacial area was estimated. The image analysis algorithm employed
here combines an anisotropic diffusion filter to remove noise from the original gray-scale image data, a k-means cluster analysis
to obtain segmented data, and the construction of isosurfaces to estimate solid surface area and interfacial area. Our method
was compared with laboratory measurements, as well as estimates obtained from a number of other image analysis algorithms
presented in the literature. Porosity estimates for the two-phase bead pack were within 1.5% error of laboratory measurements
and agreed well with estimates obtained using an indicator kriging segmentation algorithm. Additionally, our method estimated
the solid surface area of the high precision beads within 10% of the laboratory measurements, whereas solid surface area estimates
obtained from voxel counting and two-point correlation functions overestimated the surface area by 20–40%. Interfacial area
estimates for the air–water menisci contained within the capillary tubes were obtained using our image analysis algorithm,
and using other image analysis algorithms, including voxel counting, two-point correlation functions, and the porous media
marching cubes. Our image analysis algorithm, and other algorithms based on marching cubes, resulted in errors ranging from
1% to 20% of the analytical interfacial area estimates, whereas voxel counting and two-point correlation functions overestimated
the analytical interfacial area by 20–40%. In addition, the sensitivity of the image analysis algorithms on the resolution
of the microtomographic image data was investigated, and the results indicated that there was little or no improvement in
the comparison with laboratory estimates for the resolutions and conditions tested. 相似文献
11.
The postprocessing algorithm introduced by Yao for imposing the spectral amplitudes of a target covariance model is shown to be efficient in correcting the smoothing effect of estimation maps, whether obtained by kriging or any other interpolation technique. As opposed to stochastic simulation, Yao's algorithm yields a unique map starting from an original, typically smooth, estimation map. Most importantly it is shown that reproduction of a covariance/semivariogram model (global accuracy) is necessarily obtained at the cost of local accuracy reduction and increase in conditional bias. When working on one location at a time, kriging remains the most accurate (in the least squared error sense) estimator. However, kriging estimates should only be listed, not mapped, since they do not reflect the correct (target) spatial autocorrelation. This mismatch in spatial autocorrelation can be corrected via stochastic simulation, or can be imposed a posteriori via Yao's algorithm. 相似文献
12.
A combinatorial filtering method for magnetotelluric data series with strong interference 总被引:1,自引:0,他引:1
Cai Jian-hua 《Arabian Journal of Geosciences》2016,9(13):628
In highly industrialized areas, magnetotelluric (MT)-induced variations are contaminated by strong manmade noise signals. A method is described as an alternative approach for noise removal, based on a combination of empirical mode decomposition (EMD) with independent component analysis (ICA). The filtering procedure takes advantage of the fact that data are analyzed through different scale levels, which requires a minimum of human intervention and leaves good data sections unchanged. Principle and steps of method are discussed, and de-noising results are evaluated by some parameters. After the filtering stage, data is processed in the frequency domain to yield two sets of reliable MT transfer functions and the result was compared with that of the EMD-Wavelet method. Simulated signal and measured MT data series are processed. The results show that this procedure can lead to greatly improved apparent resistivity and phase curves after processing. Point defects are filtered out to eliminate their deleterious influence, which yields reliable estimates of the MT transfer functions. The EMD-ICA method provides a new method for the de-noising of MT data series under the condition of low SNR. 相似文献
13.
The Gejiu tin-copper polymetallic ore field, located at the westernmost end of the Cathaysia Block in South China, is one of the largest tin polymetallic ore fields in the world. It is associated with magmatic-hydrothermal ore-forming processes triggered by deeply seated structures and concealed granites. The Bi-dimensional empirical mode decomposition (BEMD) was used to decompose gravity data covering an area including the Gejiu tin-copper polymetallic ore field in the west and the Bozhushan silver-lead-zinc polymetallic ore field in the east, which yielded three two-dimensional intrinsic mode function (BIMFs) images and one residue (Res(m, n)) image that depict four layers of geological architectures at different wavelengths within the study area. The high-pass filtered gravity component image (BIMF1) is interpreted to depict the shallow geological architecture, which indicates that the skarn alteration and tin-copper mineralization with positive gravity anomalies are distributed around the granites characterized by negative gravity anomalies. The band-pass filtered gravity component image (BIMF2) is interpreted to depict the middle-shallow geological architecture, which indicates that the outcropped granites in the western Gejiu ore field bounded by the Gejiu fault may extend to the eastern Gejiu ore field to form one integrated granitic complex in the subsurface. The other band-pass filtered gravity component image (BIMF3) is interpreted to depict the middle-lower geological architecture, which indicates that there may be an EW-trending granite zone with negative gravity anomaly at middle-lower depth connecting the Gejiu tin-copper polymetallic field to the Bozhushan silver-lead-zinc polymetallic field, displaying the existence of an EW-trending regional polymetallic ore-forming zone within the study area. The low-pass filtered gravity component image (Res(m, n)) is interpreted to depict the lowest geological architecture within the study area, which reflects that there may be a pair of the NW-trending uplift zone of the mantle and/or the basement with positive gravity anomaly and the depression zone of the mantle and/or the basement with negative gravity anomaly. The Gejiu tin-copper polymetallic deposits and the related granites are located at the transitional zone between the uplift and depression, whereas the Bozhushan silver-lead-zinc polymetallic deposits and the related granites are situated within the depression zone, which implies that the diversity of the Late Yanshanian granites and the related polymetallic deposits in the study area may be controlled by the complexity of the crust-mantle interaction at depth. 相似文献
14.
This study compares kriging and maximum entropy estimators for spatial estimation and monitoring network design. For second-order
stationary random fields (a subset of Gaussian fields) the estimators and their associated interpolation error variances are
identical. Simple lognormal kriging differs from the lognormal maximum entropy estimator, however, in both mathematical formulation
and estimation error variances. Two numerical examples are described that compare the two estimators. Simple lognormal kriging
yields systematically higher estimates and smoother interpolation surfaces compared to those produced by the lognormal maximum
entropy estimator. The second empirical comparison applies kriging and entropy-based models to the problem of optimizing groundwater
monitoring network design, using six alternative objective functions. The maximum entropy-based sampling design approach is
shown to be the more computationally efficient of the two. 相似文献
15.
In Particle Image Velocimetry the analysis of the multi-exposed photographic plate can be carried out either by means of the Young's fringe analysis or by a direct technique. With this second approach the image of each interrogation area, enlarged by a microscope, is recorded by means of a digital image acquisition system. After some preliminary filtering operations the digitized image can be analyzed with numerical tools. In this paper a review of the digital techniques used in the analysis of these images is carried out. Both twodimensional and onedimensional techniques are analyzed and the problems involved in their use are discussed. 相似文献
16.
基于地形分区IDW的复垦排土场土壤养分空间分布 总被引:1,自引:0,他引:1
准确掌握土壤养分的空间分布,对于评价复垦排土场土壤质量,并制定相应的复垦措施具有重要意义,但是排土场的复杂地形导致一般插值方法的预测效果较差。以内蒙古锡林浩特盟胜利一号露天矿北排土场为例,采用基于地形分区的反距离权重法对土壤养分(速效氮、速效钾、有效磷和有机质)进行空间插值并与常规反距离权重法和克里金插值法进行比较。结果表明,使用分区反距离权重法插值得到的交叉验证精度均高于克里金插值法与反距离权重法,除此之外,有效磷的R2值为0.22,速效氮的R2值为0.37,速效钾以及有机质的R2值分别为0.34与0.14,同样高于克里金插值法与反距离权重法。插值结果显示,排土场速效氮含量整体由北至南逐渐减少;有效磷含量仅在中心地区稍高;速效钾在研究区东北部含量偏高;土壤有机质在中心地区含量相对偏低。研究成果为露天矿排土场土壤养分评价提供借鉴。 相似文献
17.
Ali Akbar Daya 《Journal of the Geological Society of India》2014,83(5):567-576
Most significant iron ore deposits in Iran are located in Central Iran Zone. These deposits belong to the Bafq mining district. The Bafq mining district is located in the Early Cambrian Kashmar-Kerman volcanic arc of Central Iran. Linear estimation of regionalized variables (for example by inverse distance weighting or ordinary Kriging) results in relatively high estimation variances, i.e. the estimates have very low precision. Assessment of project economics (or other critical decision making) based on linear estimation is therefore risky. Non-linear estimation methods like disjunctive kriging perform better and the lower estimation variance allows less risky economic decision-making. Another advantage of disjunctive kriging is that it allows estimation of functions of the primary variable, which here is the grade (Fe %) of the ore. In particular it permits estimation of indicator functions defined using thresholds on the primary variable. This paper is devoted to application of disjunctive kriging method in Choghart North Anomaly iron ore deposit in Central Iran, Yazd province, Iran. In this study, the Fe concentration of Choghart North Anomaly iron ore deposit was modelled and estimated. The exploration data consists of borehole samples measuring the Fe concentration. A Gaussian isofactorial model is fitted to these data and disjunctive kriging was used to estimate the regionalized variable (Fe %) at unsampled locations and to assess the probabilities that the actual concentrations exceed a threshold value at a given location. Consequently a three dimensional model of probability of exceeding a threshold value and the estimated value are provided by disjunctive kriging to divide the ore into an economic and uneconomic part on the basis of estimation of indicator functions using thresholds grades defined on point support. The tools and concepts are complemented by a set of computer programs that are applied to the case study. The study showed that disjunctive kriging can be applied successfully for modeling the grade of an ore deposit. Results showed that the correlation between the estimated value and real value at locations close to each other is 81.9%. 相似文献
18.
本文在微型计算机上实现了对研究区布格重力异常的数字滤波处理,并用假彩色图象的数字干涉技术,详细研究了地壳的现代构造格局和形变踪迹。追踪描述了由代表着地体边界断裂的东西向大型走滑断层的拉伸作用而形成的NS向正断层。识别出由NE向、NW向构造构成的走滑断裂系统。走汾构造的活动是使青藏高原南北向缩短和地震活动的主要因素之一。突出深部异常的滤波图象显示出,麻江一比如的NE向走滑断裂至少切穿了上地壳。同时也指示出深部异常的最低值位于羌塘地体,从地壳均衡角度分析,可能意味着唐古拉山有很深的山根。 相似文献
19.
Masoomeh Delbari 《Arabian Journal of Geosciences》2014,7(10):4221-4233
This paper presents the incorporation of a digital elevation model into the spatial prediction of water table elevation in Mazandaran province (Iran) using a range of interpolation techniques. The multivariate methods used are: linear regression (LR), cokriging (COK), kriging with an external drift (KED) and regression kriging (RK). The analysis is performed on 3 years (1987, 1997 and 2007) of water table elevation data from about 260 monitoring wells. Prediction performances of the different algorithms are compared with two univariate techniques, i.e. inverse distance weighting and ordinary kriging (OK), through cross validation and examination of the consistency of the generated maps with the natural phenomena. Significantly smaller prediction errors are obtained for four multivariate algorithms but, in particular, KED and RK outperform LR and COK for 3 years. The results show the potential for using elevation for a more precise mapping of water table elevation. 相似文献
20.
以表层带岩溶洞穴Velika Pasica溶洞内两个滴水点为研究对象,利用数字滤波法、最小平滑值法,以及两者结合的最小平滑值滤波法,对岩溶表层带水流进行基流分割分析。结果表明:对于特殊的岩溶表层带含水层,递归数字滤波法在滤波参数值(α)为0.925和0.95时,经过正反正三次滤波处理后得到比较好的基流分割效果;经过最小平滑值法处理后再进行一次数字滤波处理,得到结果与在α=0.95时经过三次滤波得到的结果相近。表层带基流分割处理的研究为其调蓄机制的研究提供科学依据。 相似文献