Recent literature in cartography and geographic information science

Tissot's Indicatrix and regular grids have been used for assessing map projection accuracies. Despite their broad applicability for accuracy assessment, they have limitations in quantifying resampling errors caused by map projections. This is due to the structural uncertainty with regard to the placement and pattern of grids. It is also difficult to calculate the absolute amount of resampling error in each projection. As an alternative to traditional testing methods, the use of random points was investigated. Specifically, random point generation, resampling with spherical block search algorithms, resampling accuracy with a perfect grid, and resampling accuracy with eight projections were investigated and are discussed here. Eight global referencing methods were tested: the equal-area cylindrical, sinusoidal, Mollweide, Eckert IV, Hammer-Aitoff, interrupted Goode homolosine, integerized sinusoidal projections, and the equal area global gridding with a fixed latitudinal metric distance. The resampling accuracy with a perfect grid is about 75 percent. Results showed the sinusoidal and the integerized sinusoidal projections and equal-area global gridding to achieve the highest accuracies.     

High Resolution Köppen‐Geiger Classifications of Paleoclimate Simulations

The development and application of an algorithm to compute Köppen‐Geiger climate classifications from the Coupled Model Intercomparison Project (CMIP) and Paleo Model Intercomparison Project (PMIP) climate model simulation data is described in this study. The classification algorithm was applied to data from the PMIP III paleoclimate experiments for the Last Glacial Maximum, 21k years before present (yBP), Mid‐Holocene (6k yBP) and the Pre‐Industrial (0k yBP, control run) time slices. To infer detailed classification maps, the simulation datasets were interpolated to a higher resolution. The classification method presented is based on the application of Open Source Software, and the implementation is described with attention to detail. The source code and the exact input data sets as well as the resulting data sets are provided to enable the application of the presented approach.     

Modeling errors in upward continuation for INS gravity compensation

Accurate upward continuation of gravity anomalies supports future precision, free-inertial navigation systems, since the latter cannot by themselves sense the gravitational field and thus require appropriate gravity compensation. This compensation is in the form of horizontal gravity components. An analysis of the model errors in upward continuation using derivatives of the standard Pizzetti integral solution (spherical approximation) shows that discretization of the data and truncation of the integral are the major sources of error in the predicted horizontal components of the gravity disturbance. The irregular shape of the data boundary, even the relatively rough topography of a simulated mountainous region, has only secondary effect, except when the data resolution is very high (small discretization error). Other errors due to spherical approximation are even less important. The analysis excluded all measurement errors in the gravity anomaly data in order to quantify just the model errors. Based on a consistent gravity field/topographic surface simulation, upward continuation errors in the derivatives of the Pizzetti integral to mean altitudes of about 3,000 and 1,500 m above the mean surface ranged from less than 1 mGal (standard deviation) to less than 2 mGal (standard deviation), respectively, in the case of 2 arcmin data resolution. Least-squares collocation performs better than this, but may require significantly greater computational resources.     

Supervised SAR Image MPM Segmentation Based on Region-Based Hierarchical Model

This letter presents a novel method of supervised multiresolution segmentation for synthetic aperture radar images. The method uses a region-based half-tree hierarchical Markov random field model for multiresolution segmentation. To form the region-based multilayer model, the watershed algorithm is employed at each resolution level independently. The nodes of a quadtree in the proposed model are defined as regions instead of pixels. The relationship over scale is studied, and the region-based upward and downward maximization of posterior marginal estimations are deduced. The experimental results for the segmentation of homogeneous areas prove the region-based model much better in terms of robustness to speckle and preservation of edges compared to the pixel-based hierarchical model and the Gibbs sampler with the single-resolution model     

Geomasking sensitive health data and privacy protection: an evaluation using an E911 database

Geomasking is used to provide privacy protection for individual address information while maintaining spatial resolution for mapping purposes. Donut geomasking and other random perturbation geomasking algorithms rely on the assumption of a homogeneously distributed population to calculate displacement distances, leading to possible under-protection of individuals when this condition is not met. Using household data from 2007, we evaluated the performance of donut geomasking in Orange County, North Carolina. We calculated the estimated k-anonymity for every household based on the assumption of uniform household distribution. We then determined the actual k-anonymity by revealing household locations contained in the county E911 database. Census block groups in mixed-use areas with high population distribution heterogeneity were the most likely to have privacy protection below selected criteria. For heterogeneous populations, we suggest tripling the minimum displacement area in the donut to protect privacy with a less than 1% error rate.     

A community-based urban forest inventory using online mapping services and consumer-grade digital images

Community involvement in gathering and submitting spatially referenced data via web mapping applications has recently been gaining momentum. Urban forest inventory data analyzed by programs such as the i-Tree ECO inventory method is a good candidate for such an approach. In this research, we tested the feasibility of using spatially referenced data gathered and submitted by non-professional individuals through a web application to augment urban forest inventory data. We examined the use of close range photogrammetry solutions of images taken by consumer-grade cameras to extract quantitative metric information such as crown diameter, tree heights and trunk diameters. Several tests were performed to evaluate the accuracy of the photogrammetric solutions and to examine their use in addition to existing aerial image data to supplement or partially substitute for standard i-Tree ECO field measurements.     

Kinematic precise point positioning at remote marine platforms

Precise kinematic differential positioning using the global positioning system (GPS) at a marine platform usually requires a relatively short distance (e.g. <500 km) to a land-based reference station. As an alternative, precise point positioning (PPP) is normally considered free from this limiting requirement. However, due to the prerequisite of network-based satellite products, PPP at a remote marine platform may still be affected by its distance to the reference network. Hence, this paper investigates this scenario by configuring rings of reference stations with different radii centered on a to-be-positioned marine platform. Particularly, we applied ambiguity resolution at a single station to PPP by estimating uncalibrated phase delays (UPDs). We used three rings of reference stations centered on a vessel, with radii of roughly 900, 2,000 and 3,600 km, to determine satellite clocks and UPDs independently. For comparison, we also performed differential positioning based on a single reference station with baseline lengths of about 400, 1,700 and 2,800 km. We demonstrate that, despite the increasing ring-network radius to a few 1,000 km, the overall change in accuracy of the satellite clocks that are used at the vessel is smaller than 0.02 ns, and the RMS values of differences between the three sets of narrow-lane UPD estimates are around 0.05 cycles only. Moreover, the kinematic positioning accuracy of PPP is affected by the increasing ring-network radius, but can still achieve several centimeters after ambiguity resolution when the vessel is over a few 1,000 km away from the ring network, showing better performance than that of differential positioning. Therefore, we propose that ambiguity-fixed PPP can be used at remote marine platforms that support precise oceanographic and geophysical applications in open oceans.     

Contribution of new AAM data source to δLOD excitation

Data sets of the changes of the length of day (ΔLOD) measured by space geodetic techniques and of the atmospheric angular momentum (AAM) derived from global meteorological data by the National Meteorological Center (NMC) and the National Centers for Environmental Prediction/National Center for Atmospheric Research (NCEP/NCAR) are used to reanalyze and study the excitations of ΔLOD, and to examine and compare the contribution of each AAM series to ΔLOD. The main results are as follows. 1. The AAM reanalyzed with the assimilated global meteorological data by NCEP/NCAR are more accurate and have lower noise than the original AAM derived by NMC. The NCEP/NCAR-based AAM is more consistent with the changes of the length-of-day series. 2. The NCEP reanalysed AAM data may better explain the non-tidal LOD variations on intraseasonal to interannual time scales, especially on the quasi-biennial time scale. The atmosphere cannot, however, explain all LOD variations; other excitation sources are possible. 3. The effects of atmosphere on the estimated values of tides for ΔLOD components up to a month are weak. The preliminary results of the annual and semiannual tides can be estimated after removing the effects of atmosphere from ΔLOD. Received: 27 May 1998 / Accepted: 22 March 1999     

Progressive Transformation of a Flux Rope to an ICME

The solar wind conditions at one astronomical unit (AU) can be strongly disturbed by interplanetary coronal mass ejections (ICMEs). A subset, called magnetic clouds (MCs), is formed by twisted flux ropes that transport an important amount of magnetic flux and helicity, which is released in CMEs. At 1 AU from the Sun, the magnetic structure of MCs is generally modeled by neglecting their expansion during the spacecraft crossing. However, in some cases, MCs present a significant expansion. We present here an analysis of the huge and significantly expanding MC observed by the Wind spacecraft during 9 – 10 November 2004. This MC was embedded in an ICME. After determining an approximate orientation for the flux rope using the minimum variance method, we obtain a precise orientation of the cloud axis by relating its front and rear magnetic discontinuities using a direct method. This method takes into account the conservation of the azimuthal magnetic flux between the inbound and outbound branches and is valid for a finite impact parameter (i.e., not necessarily a small distance between the spacecraft trajectory and the cloud axis). The MC is also studied using dynamic models with isotropic expansion. We have found (6.2±1.5)×10²⁰ Mx for the axial flux and (78±18)×10²⁰ Mx for the azimuthal flux. Moreover, using the direct method, we find that the ICME is formed by a flux rope (MC) followed by an extended coherent magnetic region. These observations are interpreted by considering the existence of a previously larger flux rope, which partially reconnected with its environment in the front. We estimate that the reconnection process started close to the Sun. These findings imply that the ejected flux rope is progressively peeled by reconnection and transformed to the observed ICME (with a remnant flux rope in the front part).     

Estimation of Partial Pressure of Carbon Dioxide and Air-Sea Fluxes in Hooghly Estuary Based on In Situ and Satellite Observations

An empirical model is developed and used with remotely sensed predictors: sea surface temperature (SST) and chlorophyll-a concentration (Chl-a), to compute surface water partial pressure of carbon dioxide (pCO_2w) and air-sea fluxes of CO₂ in the Hooghly estuary and its adjacent coastal oceans. In situ observations used here were based on measurements carried out in this region during winter and summer periods in 2008. The estimated pCO_2w compares well with the in situ observations at root mean square error ±18 μatm. In winter, estimated pCO_2w ranges between 320 and 500 μatm with large values (>400 μatm) on the south-western and south-eastern flanks of the coastal domain and lower values (340–375 μatm) on the main-channel. In summer, it remained spatially uniform at 450 μatm. Extrapolation of the results over the study region based on the Moderate Imaging Specroradiometer (MODIS) measured SST and Chl-a suggests that the region is a strong source of atmospheric CO₂ during the summer with net release of 0.095 Tg C year^?1 (equivalent to mean flux of 90 molC m^?2 year^?1) and is a weak source during the winter with net release of 0.006 Tg C yr^?1 (0.5 molC m^?2 year^?1) from the geographical extent of 6000 Km² area.