航空重力测量中FIR低通滤波器的设计与应用

航空重力测量数据含有许多噪声，需用高质量的数字滤波器加以平滑。有限冲激响应（ＦＩＲ）滤波器的优点是具有精确线性相位而且总是稳定的。本文利用契比雪夫逼近方法，设计了适用于航空重力测量的ＦＩＲ低通滤波器。为检测其性能，文中建立了一个模拟加速度模型。检测结果表明：滤波器性能可靠，适用于航空重力测量数据处理。     

国外物探仪器新进展

本文介绍了国外航空、地面及井中物探仪器的最新进展，和一些新型仪器的性能。     

Modeling Mediterranean forest structure using airborne laser scanning data

The conservation of biological diversity is recognized as a fundamental component of sustainable development, and forests contribute greatly to its preservation. Structural complexity increases the potential biological diversity of a forest by creating multiple niches that can host a wide variety of species. To facilitate greater understanding of the contributions of forest structure to forest biological diversity, we modeled relationships between 14 forest structure variables and airborne laser scanning (ALS) data for two Italian study areas representing two common Mediterranean forests, conifer plantations and coppice oaks subjected to irregular intervals of unplanned and non-standard silvicultural interventions. The objectives were twofold: (i) to compare model prediction accuracies when using two types of ALS metrics, echo-based metrics and canopy height model (CHM)-based metrics, and (ii) to construct inferences in the form of confidence intervals for large area structural complexity parameters.Our results showed that the effects of the two study areas on accuracies were greater than the effects of the two types of ALS metrics. In particular, accuracies were less for the more complex study area in terms of species composition and forest structure. However, accuracies achieved using the echo-based metrics were only slightly greater than when using the CHM-based metrics, thus demonstrating that both options yield reliable and comparable results. Accuracies were greatest for dominant height (Hd) (R² = 0.91; RMSE% = 8.2%) and mean height weighted by basal area (R² = 0.83; RMSE% = 10.5%) when using the echo-based metrics, 99th percentile of the echo height distribution and interquantile distance. For the forested area, the generalized regression (GREG) estimate of mean Hd was similar to the simple random sampling (SRS) estimate, 15.5 m for GREG and 16.2 m SRS. Further, the GREG estimator with standard error of 0.10 m was considerable more precise than the SRS estimator with standard error of 0.69 m.     

Effects of topographic–isostatic masses on gravitational functionals at the Earth’s surface and at airborne and satellite altitudes

Topographic–isostatic masses represent an important source of gravity field information, especially in the high-frequency band, even if the detailed mass-density distribution inside the topographic masses is unknown. If this information is used within a remove-restore procedure, then the instability problems in downward continuation of gravity observations from aircraft or satellite altitudes can be reduced. In this article, integral formulae are derived for determination of gravitational effects of topographic–isostatic masses on the first- and second-order derivatives of the gravitational potential for three topographic–isostatic models. The application of these formulas is useful for airborne gravimetry/gradiometry and satellite gravity gradiometry. The formulas are presented in spherical approximation by separating the 3D integration in an analytical integration in the radial direction and 2D integration over the mean sphere. Therefore, spherical volume elements can be considered as being approximated by mass-lines located at the centre of the discretization compartments (the mass of the tesseroid is condensed mathematically along its vertical axis). The errors of this approximation are investigated for the second-order derivatives of the topographic–isostatic gravitational potential in the vicinity of the Earth’s surface. The formulas are then applied to various scenarios of airborne gravimetry/gradiometry and satellite gradiometry. The components of the gravitational vector at aircraft altitudes of 4 and 10 km have been determined, as well as the gravitational tensor components at a satellite altitude of 250 km envisaged for the forthcoming GOCE (gravity field and steady-state ocean-circulation explorer) mission. The numerical computations are based on digital elevation models with a 5-arc-minute resolution for satellite gravity gradiometry and 1-arc-minute resolution for airborne gravity/gradiometry.     

Classifying airborne radiometry data with Agglomerative Hierarchical Clustering: A tool for geological mapping in context of rainforest (French Guiana)

In highly weathered environments, it is crucial that geological maps provide information concerning both the regolith and the bedrock, for societal needs, such as land-use, mineral or water resources management. Often, geologists are facing the challenge of upgrading existing maps, as relevant information concerning weathering processes and pedogenesis is currently missing. In rugged areas in particular, where access to the field is difficult, ground observations are sparsely available, and need therefore to be complemented using methods based on remotely sensed data.     

Sea Surface Determination Using Long-Range Kinematic GPS Positioning and Laser Airborne Depth Sounder Techniques

Precise long-range kinematic GPS positioning requires the use of carrier phase measurements, the data processing of which suffers from the technical challenges of ambiguity resolution and cycle slip repair. In this paper, the combination of an ambiguity recovery technique and a linear bias correction method has been used to overcome such problems. An experiment was conducted to test the utility of this technique to determine aircraft height to high accuracy, over very long baselines (of the order of one thousand kilometres), in support of the Laser Airborne Depth Sounder (LADS). From a comparison of four independently derived trajectories, this airborne GPS kinematic positioning experiment has confirmed that the sea surface can be determined to centimetre accuracy. The sea surface profiles thus obtained can be used to correct the errors introduced by long period ocean swells.     

航空重力傅里叶基追踪低通滤波方法研究

航空重力测量受到各种各样的高频噪声干扰,因此,低通滤波是提取重力信号的重要环节,其关键在于设计性能优越的低通滤波器.目前航空重力测量中常用FIR（Finite Impulse Response）低通滤波方法存在明显的滤波边缘效应,导致不得不舍弃边缘部分数据.针对这一问题,本文引入一种可以有效抑制边缘效应的新方法——傅里叶基追踪低通滤波方法（Fourier Basis Pursuit Low Pass Filter,FBPLPF）.该方法通过基追踪准则,选择全局优化,采用凸优化中的内点算法,将低频信号挤压在低频上,实现低频信号与高频信号的有效分离,能够有效减少有限时间序列造成的谱污染和谱泄漏.最后利用仿真实验和实测数据对该方法进行了验证,均方根误差（RMS）东西测线为0.7×10^-5 m·s^-2,南北测线为1.4×10^-5 m·s^-2,与FIR低通滤波方法舍弃边缘数据后统计的均方根误差相当.表明该方法可以不舍弃或者舍弃少量边缘数据,提高航空重力数据的利用率.     

Research Note: First results of a low‐altitude unmanned aircraft system gamma survey by comparison with the terrestrial and aerial gamma survey data

In order to increase the efficiency of geological survey in difficult landscape‐morphological conditions, the authors have developed low‐altitude gamma technology surveys based on a multi‐rotor unmanned aircraft system, also known informally as Drones. These results have been compared against those from standard terrestrial and aerial gamma surveys. The successful results obtained at the reference site show that the developed methodology enables compilation of high‐quality data in difficult landscape conditions. These results yield the same level of information content as standard terrestrial gamma surveys. The exact unmanned aircraft system altitude over the earth must be maintained for carrying out high‐quality surveys. It was found that a CsI(Tl) detector with sizes 80 × 80 mm (vol. 0.4 dm³) enables obtaining of high‐quality radiometric data at altitudes of up to 70 m, at a speed of at least 20 km/h. The authors believe that low altitude ‘quasi‐terrestrial’ gamma surveys with a terrain drape may in the future replace terrestrial gamma surveys since they can offer better effectiveness at a lower cost.     

Slopes of an airborne electromagnetic resistivity model interpolated jointly with borehole data for 3D geological modelling

We investigate a novel way to introduce resistivity models deriving from airborne electromagnetic surveys into regional geological modelling. Standard geometrical geological modelling can be strengthened using geophysical data. Here, we propose to extract information contained in a resistivity model in the form of local slopes that constrain the modelling of geological interfaces. The proposed method is illustrated on an airborne electromagnetic survey conducted in the region of Courtenay in France. First, a resistivity contrast corresponding to the clay/chalk interface was interpreted confronting the electromagnetic soundings to boreholes. Slopes were then sampled on this geophysical model and jointly interpolated with the clay/chalk interface documented in boreholes using an implicit 3D potential‐field method. In order to evaluate this new joint geophysical–geological model, its accuracy was compared with that of both pure geological and pure geophysical models for various borehole configurations. The proposed joint modelling yields the most accurate clay/chalk interface whatever the number and location of boreholes taken into account for modelling and validation. Compared with standard geological modelling, the approach introduces in between boreholes geometrical information derived from geophysical results. Compared with conventional resistivity interpretation of the geophysical model, it reduces drift effects and honours the boreholes. The method therefore improves what is commonly obtained with geological or geophysical data separately, making it very attractive for robust 3D geological modelling of the subsurface.