Excel VBA在测量数据处理中的应用

Excel是office办公软件中的一个组件,而VBA是office应用程序共享的一种通用的自动化语言,可以使常用的office应用程序自动化,并可以创建自定义的解决方案。本文对Excel2003的数据处理功能进行分析,论述通过Excel2003的VBA编辑器创建"宏"功能实现测量数据处理的优点;并结合道路工程中桩坐标及其放样数据计算实例,给出了通过VBA编程进行测量数据处理的方法。     

史前连续文化聚落的决策树分类挖掘研究——以郑州-洛阳地区为例

利用C4.5决策树分类算法和GIS组件Super Map Objects,开发实现了基于连续文化序列的空间数据挖掘系统,运用此系统对郑州-洛阳地区史前四个连续文化时期的聚落进行了数据挖掘,提取了聚落遗址面积的分类规则,并根据这些分类规则和聚落的空间分布情况对四个文化时期的I级聚落进行对比分析,得出了I级聚落的空间分布特点,即I级聚落中的一些特大面积的聚落朝着中心聚落的方向发展。     

基于LBV变换的遥感影像多步骤分类法研究

面对与日俱增的遥感数据,如何快速准确地进行影像分类成为遥感领域亟待解决的问题。本文在前人的基础上提出了一种基于LBV变换的遥感影像多步骤分类方法。L:地物的总辐射水平,集中反映了裸地的信息;B:地物的可见-红外光辐射平衡,是地面水分状况和水体存在的一个良好指标;V:地物辐射随波段变化的方向和速度向量,能集中反映地面植被状况。多步骤分类法遵循由易到难的原则,能克服类别之间的互相影响,从而提高分类精度。文章结合LBV变换和多步骤分类法的优点对影像进行了分类,结果表明该方法简单易行,且能达到良好的分类精度。     

遥感图像线性影纹理解专家系统设计与实现

线性影纹信息是遥感图像中的一类重要信息,线性影纹信息自动提取是遥感图像智能解译的重要研究领域。笔者采用基于特征对象的专家系统技术来完成线性影纹信息的简单类别提取,设计并实现了线性影纹理解专家系统;阐明了系统结构,解译规则获取、表达,基于消息的不确定性推理机及解释机的设计;最后利用系统做了线性影纹类别提取实验,结果表明基于特征对象的专家系统用于线性影纹类别自动提取是切实可行的,并具有较高的提取精度。     

用MapGis实现区域化探数据的空间分析

用VB程序,把区域化探数据点位信息转化为MapGis明码格式的点文件,再与化探数据连接,完成空间数据库的建立。笔者以河北省1∶20万区域化探数据为例,建立了省级区域化探空间数据库,利用该数据库在MapGis上进行空间分析———点对地质区的相交分析、条件查询、缓冲区分析等,实现了基于空间和地质条件约束的数据查询,解决了与成矿有关的地质问题。     

Using spectral information from the NIR water absorption features for the retrieval of canopy water content

Canopy water content (CWC) is important for mapping and monitoring the condition of the terrestrial ecosystem. Spectral information related to the water absorption features at 970 nm and 1200 nm offers possibilities for deriving information on CWC. In this study, we compare the use of derivative spectra, spectral indices and continuum removal techniques for these regions. Hyperspectral reflectance data representing a range of canopies were simulated using the combined PROSPECT + SAILH model. Best results in estimating CWC were obtained by using spectral derivatives at the slopes of the 970 nm and 1200 nm water absorption features. Real data from two different test sites were analysed. Spectral information at both test sites was obtained with an ASD FieldSpec spectrometer, whereas at the second site HyMap airborne imaging spectrometer data were also acquired. Best results were obtained for the derivative spectra. In order to avoid the potential influence of atmospheric water vapour absorption bands the derivative of the reflectance on the right slope of the canopy water absorption feature at 970 nm can best be used for estimating CWC.     

Application of Spectral Decomposition to Detection of Fracture-Cavity Carbonate Reservoir Beds in the Tahe Oiifieid, Tarim Basin, NW China

Ordovician fracture-cavity carbonate reservoir beds are the major type of producing formations in the Tahe oilfield, Tarim Basin. The seismic responses of these beds clearly changes depending on the different distance of the fracture-cavity reservoir bed from the top of the section. The seismic reflection becomes weak or is absent when the fracture-cavity reservoir beds are less than 20 ms below the top Ordovician. The effect on top Ordovician reflection became weaker with deeper burial of fracture-cavity reservoir beds but the developed deep fracture-cavity reservoir beds caused stronger reflection in the interior of the Ordovician. This interior reflection can be divided into strong long-axis, irregular and bead string reflections, and was present 80 ms below the top Ordovician. Aimed at understanding reflection characteristics, the spectral decomposition technique, which uses frequency to ＂tune-in＂ bed thickness, was used to predict Ordovician fracture-cavity carbonate formations in the Tahe oilfield. Through finely adjusting the processing parameters of spectral decomposition, it was found that the slice at 30 Hz of the tuned data cube can best represent reservoir bed development. Two large N-S-trending strong reflection belts in the mid-western part of the study area along wells TK440- TK427-TK417B and in the eastern part along wells TK404-TK409 were observed distinctly on the 30 Hz slice and 4-D time-frequency data cube carving. A small N-S trending reflection belt in the southern part along wells T403-TK446B was also clearly identified. The predicted reservoir bed development area coincides with the fracture-cavities connection area confirmed by drilling pressure testing results. Deep karst cavities occur basically in three reservoir bed-development belts identified by the Ordovician interior strong reflection. Spectral decomposition proved to be a useful technique in identifying fracture-cavity reservoir beds.     

Environmental Investigation and Evaluation of Land Subsidence in the Datong Coalfield Based on InSAR Technology

Heavy mining of Jurassic and Carboniferous horizontal coal seams in the Datong coalfield has seriously affected the local geological environment, which is mainly manifested by such geohazards as soil avalanches, landslides, mudflows, surface subsidence, surface cracks, surface solid waste accumulation and surface deformation. More seriously, coal mining causes groundwater to leak. Overpumping of groundwater has resulted in substantial land subsidence of the urban area in Datong City. Based on the previous geo-environmental investigation in the work area, the authors used radar remote sensing techniques such as InSAR （synthetic aperture radar interferometry） and D-InSAR （differential synthetic aperture radar interferometry）, supplemented by the optical remote sensing method, for geo-environmental investigation to ascertain the geo-environmental background of the Datong Jurassic and Carboniferous-Permian coalfield and evaluate the effects of the geohazards, thus providing a basis for the geo-environmental protection, geohazard control and prevention, land improvement and optimization of the human environment. In this study 8 cog-nominal ERS-1/2 SAR data frames during 1992 to 2003 were used, but the following processing was made： （1） the multitemporal SAR magnitude images were used to interpret the geological structure, vegetation, microgeomorphology and drainage system; （2） the multi-temporal InSAR coherent images were used to make a classification of surface features and evaluate the coherence change due to coal mining; and （3） the multi-temporal cog-nominal SAR images were used to complete D-InSAR processing to remove the information of differential deformation areas （sites）. In the end, a ten-year time series of differential interferograms were obtained using the multi-temporal cog-nominal SAR images. In the tests, 84 deformed areas （sites） were obtained, belonging to those in 1993-1996, 1996-1997, 1997-1998, 1998-2001, 1998-2002 and 2001-2003 respectively. Of the 84 areas, 44 are m     

The GeoForschungsZentrum Potsdam/Groupe de Recherche de Gèodésie Spatiale satellite-only and combined gravity field models: EIGEN-GL04S1 and EIGEN-GL04C

The recent improvements in the Gravity Recovery And Climate Experiment (GRACE) tracking data processing at GeoForschungsZentrum Potsdam (GFZ) and Groupe de Recherche de Géodésie Spatiale (GRGS) Toulouse, the availability of newer surface gravity data sets in the Arctic, Antarctica and North-America, and the availability of a new mean sea surface height model from altimetry processing at GFZ gave rise to the generation of two new global gravity field models. The first, EIGEN-GL04S1, a satellite-only model complete to degree and order 150 in terms of spherical harmonics, was derived by combination of the latest GFZ Potsdam GRACE-only (EIGEN-GRACE04S) and GRGS Toulouse GRACE/LAGEOS (EIGEN-GL04S) mean field solutions. The second, EIGEN-GL04S1 was combined with surface gravity data from altimetry over the oceans and gravimetry over the continents to derive a new high-resolution global gravity field model called EIGEN-GL04C. This model is complete to degree and order 360 and thus resolves geoid and gravity anomalies at half- wavelengths of 55 km at the equator. A degree-dependent combination method has been applied in order to preserve the high accuracy from the GRACE satellite data in the lower frequency band of the geopotential and to form a smooth transition to the high-frequency information coming from the surface data. Compared to pre-CHAMP global high-resolution models, the accuracy was improved at a spatial resolution of 200 km (half-wavelength) by one order of magnitude to 3 cm in terms of geoid heights. The accuracy of this model (i.e. the commission error) at its full spatial resolution is estimated to be 15 cm. The model shows a reduced artificial meridional striping and an increased correlation of EIGEN-GL04C-derived geostrophic meridional currents with World Ocean Atlas 2001 (WOA01) data. These improvements have led to select EIGEN-GL04C for JASON-1 satellite altimeter data reprocessing. Electronic Supplementary Material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

A data-driven approach to local gravity field modelling using spherical radial basis functions

We propose a methodology for local gravity field modelling from gravity data using spherical radial basis functions. The methodology comprises two steps: in step 1, gravity data (gravity anomalies and/or gravity disturbances) are used to estimate the disturbing potential using least-squares techniques. The latter is represented as a linear combination of spherical radial basis functions (SRBFs). A data-adaptive strategy is used to select the optimal number, location, and depths of the SRBFs using generalized cross validation. Variance component estimation is used to determine the optimal regularization parameter and to properly weight the different data sets. In the second step, the gravimetric height anomalies are combined with observed differences between global positioning system (GPS) ellipsoidal heights and normal heights. The data combination is written as the solution of a Cauchy boundary-value problem for the Laplace equation. This allows removal of the non-uniqueness of the problem of local gravity field modelling from terrestrial gravity data. At the same time, existing systematic distortions in the gravimetric and geometric height anomalies are also absorbed into the combination. The approach is used to compute a height reference surface for the Netherlands. The solution is compared with NLGEO2004, the official Dutch height reference surface, which has been computed using the same data but a Stokes-based approach with kernel modification and a geometric six-parameter “corrector surface” to fit the gravimetric solution to the GPS-levelling points. A direct comparison of both height reference surfaces shows an RMS difference of 0.6 cm; the maximum difference is 2.1 cm. A test at independent GPS-levelling control points, confirms that our solution is in no way inferior to NLGEO2004.