Practical Experiences with Films, Cameras and Navigation Systems

The first author commenced his experiences with aerial photography in 1945, when he joined the Royal Air Force (RAF). The aerial cameras used at that time were the Williamson F24 and F52, with Ilford HP3 and FP3 panchromatic films. Subsequently the Williamson F95 camera was tested against the USAF Sonne 11 camera. During his final years in the RAF, a variety of film types was employed: colour, monochrome infra-red and false colour infra-red. Additional experience was gained when the author was employed by the International Training Centre (ITC), The Netherlands and he saw the further development of the air camera from the Wild RC8 and Zeiss (Oberkochen) RMK 15/23 to the Wild RC20, Zeiss (Jena) LMK and Zeiss (Oberkochen) RMK TOP. He predicts that the future lies with digital cameras.
The second author started working with aerial photography in the 1950s. Details are given of cameras, visual navigation sights and survey aircraft dating from that period. Overseas operations, carried out mainly by Hunting Surveys and Fairey Surveys, posed particular problems because the majority of the work there was undertaken without the aid of existing maps. The introduction of external navigation systems is described, from the Decca ship navigation system, through the GNS-200 VLF/Omega low frequency radio system to INS, the Inertial Navigation System. The author later joined the ITC, which enabled him to participate in trials of the Computer-controlled Photo Navigation System. Eventually practical experience was gained with differential GPS, particularly in countries in southeast Asia. Aerial photography currently appears to benefit less from financial investment than the photographic, aviation and survey industries. However, the satisfaction of making a significant contribution to major overseas development projects outweighs the frustration of lack of acceptance of plans for the future of the air survey industry.     

An analysis of changes in the urban landscape composition and configuration in the Sancaktepe District of Istanbul Metropolitan City,Turkey using landscape metrics and satellite data

Changes in landscape composition and configuration patterns of Sancaktepe Municipal District in the Asian side of Istanbul Metropolitan City of Turkey were analysed using landscape metrics. Class-level and landscape-level metrics were calculated from the land cover/land use data using Patch Analyst, an extension in the Arc View GIS. The land cover/land use data were derived from classified satellite images of Landsat Thematic Mapper of 2002 and 2009 for Sancaktepe District. There was evidence of increase in agglomeration process of built-up patches as indicated by the increases in mean patch size, decrease in total edge and number of patches between 2002 and 2009. The urban expansion pattern experienced overall was not fragmented but concentrated due to infilling around existing patches. Changes in Area-Weighted Mean Shape Index and Area-Weighted Patch Fractal Dimension Index indicated that the physical shapes within built-up, forest and bareland areas were relatively complex and irregular. A conclusion is made in this study that spatial metrics are useful tools to describe the urban landscape composition and configuration in its various aspects and certain decisions whether to approve a specific development in urban planning could, for example, be based on some measures of urban growth form or pattern in terms of uniformity and irregularity, attributable to the dynamic processes of agglomeration and fragmentation of land cover/land use patches caused by urban expansion.     

Separation of global time-variable gravity signals into maximally independent components

The Gravity Recovery and Climate Experiment (GRACE) products provide valuable information about total water storage variations over the whole globe. Since GRACE detects mass variations integrated over vertical columns, it is desirable to separate its total water storage anomalies into their original sources. Among the statistical approaches, the principal component analysis (PCA) method and its extensions have been frequently proposed to decompose the GRACE products into space and time components. However, these methods only search for decorrelated components that on the one hand are not always interpretable and on the other hand often contain a superposition of independent source signals. In contrast, independent component analysis (ICA) represents a technique that separates components based on assumed statistical independence using higher-order statistical information. If one assumes that independent physical processes generate statistically independent signal components added up in the GRACE observations, separating them by ICA is a reliable strategy to identify these processes. In this paper, the performance of the conventional PCA, its rotated extension and ICA are investigated when applied to the GRACE-derived total water storage variations. These analyses have been tested on both a synthetic example and on the real GRACE level-2 monthly solutions derived from GeoForschungsZentrum Potsdam (GFZ RL04) and Bonn University (ITG2010). Within the synthetic example, we can show how imposing statistical independence in the framework of ICA improves the extraction of the ‘original’ signals from a GRACE-type super-position. We are therefore confident that also for the real case the ICA algorithm, without making prior assumptions about the long-term behaviour or on the frequencies contained in the signal, improves over the performance of PCA and its rotated extension in the separation of periodical and long-term components.     

Evaluating the impact of data quantity, distribution and algorithm selection on the accuracy of 3D subsurface models using synthetic grid models of varying complexity

Testing the accuracy of 3D modelling algorithms used for geological applications is extremely difficult as model results cannot be easily validated. This paper presents a new approach to evaluate the effectiveness of common interpolation algorithms used in 3D subsurface modelling, utilizing four synthetic grids to represent subsurface environments of varying geological complexity. The four grids are modelled with Inverse Distance Weighting and Ordinary Kriging, using data extracted from the synthetic grids in different spatial distribution patterns (regular, random, clustered and sparse), and with different numbers of data points (100, 256, 676 and 1,600). Utilizing synthetic grids for this evaluation allows quantitative statistical assessment of the accuracy of both interpolation algorithms in a variety of sampling conditions. Data distribution proved to be an important factor; as in many geological situations, relatively small numbers of randomly distributed data points can generate more accurate 3D models than larger amounts of clustered data. This study provides insight for optimizing the quantity and distribution of data required to accurately and cost-effectively interpolate subsurface units of varying complexity.     

MODIS多波段数据对南极海冰变化的监测研究

以归一化雪被指数法为基础根据南极的环境特点提出了MODIS影像的0．86m与1．24m的新波段组合方法分别实现了晴空与薄云下的海冰提取大幅度地提高了南极海冰监测的效率与分辨率 结合AMSR－E微波数据进行了海冰变化研究得到2002 2010年的全南极海冰范围与净冰面积都在2月份达到最小值均值分别为3．17×106　km2与2．42×106　km2 都在9月份达到最大值均值分别为18．40×106　km2与16．60×106　km2 而全南极海冰的年际变化可分为三个阶段从2002 2004年海冰基本持平2005 2007年海冰减少2008 2010年海冰重新增长这与南极海冰20a长周期的推断相符合但9年间的全南极海冰无明显的趋势性变化     

On the use of iterative re-weighting least-squares and outlier detection for empirically modelling rates of vertical displacement

The proper identification and removal of outliers in the combination of rates of vertical displacements derived from GPS, tide gauges/satellite altimetry, and GRACE observations is presented. Outlier detection is a necessary pre-screening procedure in order to ensure reliable estimates of stochastic properties of the observations in the combined least-squares adjustment (via rescaling of covariance matrices) and to ensure that the final vertical motion model is not corrupted and/or distorted by erroneous data. Results from this study indicate that typical data snooping methods are inadequate in dealing with these heterogeneous data sets and their stochastic properties. Using simulated vertical displacement rates, it is demonstrated that a large variety of outliers (random scattered and adjacent, as well as jointly influential) can be dealt with if an iterative re-weighting least-squares adjustment is combined with a robust median estimator. Moreover, robust estimators are efficient in areas weakly constrained by the data, where even high quality observations may appear to be erroneous if their estimates are largely influenced by outliers. Four combined models for the vertical motion in the region of the Great Lakes are presented. The computed vertical displacements vary between  − 2 mm/year (subsidence) along the southern shores and 3 mm/year (uplift) along the northern shores. The derived models provide reliable empirical constraints and error bounds for postglacial rebound models in the region.     

Non-stationary covariance function modelling in 2D least-squares collocation

Standard least-squares collocation (LSC) assumes 2D stationarity and 3D isotropy, and relies on a covariance function to account for spatial dependence in the observed data. However, the assumption that the spatial dependence is constant throughout the region of interest may sometimes be violated. Assuming a stationary covariance structure can result in over-smoothing of, e.g., the gravity field in mountains and under-smoothing in great plains. We introduce the kernel convolution method from spatial statistics for non-stationary covariance structures, and demonstrate its advantage for dealing with non-stationarity in geodetic data. We then compared stationary and non- stationary covariance functions in 2D LSC to the empirical example of gravity anomaly interpolation near the Darling Fault, Western Australia, where the field is anisotropic and non-stationary. The results with non-stationary covariance functions are better than standard LSC in terms of formal errors and cross-validation against data not used in the interpolation, demonstrating that the use of non-stationary covariance functions can improve upon standard (stationary) LSC.     

On the topographic bias in geoid determination by the external gravity field

The topographic bias is defined as the error/bias committed by continuing the external gravity field inside the topographic masses by a harmonic function. We study the topographic bias given by a digital terrain model defined by a spherical template, and we show that the topographic bias is given only by the potential of an inner-zone cap, and it equals the bias of the Bouguer shell, independent of the size of the cap. Then we study the effect on the real Earth by decomposing its topography into a template, and we show also in this case that the topographic bias is that of the Bouguer shell, independent of the shape of the terrain. Finally, we show that the topographic potential of the terrain at the geoid can be determined to any precision by a Taylor expansion outside the Earth’s surface. The last statement is demonstrated by a Taylor expansion to fourth order.