最大后验概率重建算法在发射CT中的应用

最大似然期望最大(Maximum Likelihood-Expectation Maximization, MLEM)重建算法已在发射CT的二维重建中得到了广泛的研究和应用,但MLEM有两点不足:1)一定迭代次数以后噪声随着迭代增多而增多.2)重建需要太长时间.而最大后验概率重建算法(maximum a posteriori, MAP)根据图像的先验知识分布,则可以有效的抑制噪声,后来逐渐发展起来的多种加速迭代方法结合MAP重建算法可以缩短重建时间.本文对近年来MAP重建算法的发展做了一个回顾,并对各类算法了做了比较和讨论.     

Responses of phreatophyte transpiration to falling water table in hyper-arid and arid regions,Northwest China

Quantitative assessment of the impact of groundwater depletion on phreatophytes in (hyper-) arid regions is key to sustainable groundwater management. However, a parsimonious model for predicting the response of phreatophytes to a decrease of the water table is lacking. A variable saturated flow model, HYDRUS-1D, was used to numerically assess the influences of depth to the water table (DWT) and mean annual precipitation (MAP) on transpiration of groundwater-dependent vegetation in (hyper-) arid regions of northwest China. An exponential relationship is found for the normalized transpiration (a ratio of transpiration at a certain DWT to transpiration at 1 m depth, T_a*) with increasing DWT, while a positive linear relationship is identified between T_a* and annual precipitation. Sensitivity analysis shows that the model is insensitive to parameters, such as saturated soil hydraulic conductivity and water stress parameters, indicated by an insignificant variation (less than 20% in most cases) under ± 50% changes of these parameters. Based on these two relationships, a universal model has been developed to predict the response of phreatophyte transpiration to groundwater drawdown for (hyper-) arid regions using MAP only. The estimated T_a* from the model is reasonable by comparing with published measured values.© 2021 China Geology Editorial Office.     

Multiparametric Cartographic Visualisation of Glacier Rheology

Abstract

In cooperation between remote sensing experts and cartographers interested in glaciology, new types of maps showing the glacier dynamics have been developed. The maps make use of the original phase gradient approach to glacier rheology modelling based on repeat-pass ERTS SAR interferograms. Careful map design and, in particular, colour assignment allow the visualisation of the glacier dynamics in its locally changing velocity with an estimated accuracy of approximately 2.0 cm per day. Two map derivates – a differential interferogram showing the glacier velocity and another product displaying the glacier strain rate – have been designed. Moreover, maps displaying the glacier marginal changes within the space of four years have been generated. The strain rate maps evidence that spots with high values frequently correspond with crevasse-prone areas which are even detectable under thick layers of snow. In this sense, the latter visualisations can be seen as maps of crevasse danger zones. The Svartisen in Norway and the Hintereis Glacier in Austria served as testbeds for the development of these different types of maps which are at the scales of 1:25 000, 1:50 000 and 1:100 000.     

Drainage Typification Based on Dendritic Decomposition

Abstract

Map generalisation is an abstraction process that seeks to transform the representation of cartographic objects from the original version into a coarser one. The characteristics of cartographic objects and the arrangement of map features have to be observed and preserved in a generalisation process. A method is developed for typifying drainages while preserving their structural characteristics, i.e.presenting the drainages with reduced number of rivers under the constraint of preserving the original structure in terms of the type and distribution of the rivers. We apply Töpfer's radical law to calculate the amount of the rivers to be retained on the generalised map. The drainages share the amount of retained rivers in proportion to the number of their tributaries. In each of the drainages, the shared amount is divided among the rivers based on the dendritic decomposition of the drainage. We implement and test the method in Java Environment. Results from case studies show that the method effectively preserves the original structures of the drainages on the generalised maps.     

'Be Map Aware: Cutting-edge…Not Close to the Edge'

Abstract

An assessment is made of the accuracy of the map and written survey, which is essentially a list of land-holders and the areas held by each, and of their mutual consistency. The geometry of the map is compared quantitatively with that of the first edition (1852) Ordnance Survey 6-inch (1 : 10 560) map. The Saxton map is found to have a scale of 1 : 5140, leading to the conclusion that the perch used in the scale bar represents 7¹/₉ yards or thereabouts. The areas quoted in the written survey are consistent with acres based on this value. The positional accuracy of points on the map is found to be ～20 m, and areas of typical plots determined from the map are uncertain by between about 4 and 12%. The relationship between these two facts can be understood through a detailed study of the statistics of the errors in separations of points on the map. The written survey is found to omit a substantial number of plots of land from its listings. The map thus appears to be more reliable overall than the written survey.     

Map use: reading,analysis, interpretation

This paper presents a method of combining text and icon label placement in a real-time computing environment. The method computes label configurations based on placement properties, cartographic disturbance, and label overlap. The process is divided into four phases. In the first phase, candidate positions of the text labels are chosen. In the second phase, the same is done for the icon labels. The choice of candidate positions is based on cartographic preference and cartographic disturbance. The removal of overlap between labels is solved, in the third phase, by means of a combinatorial optimization technique (simulated annealing). When there are label pairs in conflict that could not be resolved, the fourth and final step is executed to remove one label in the pair. The success of the proposed method lies in the ability to effectively reduce the search space for the combinatorial optimization. A number of strategies for reducing search space have been evaluated in a case study. The results show that a good search-space-reduction strategy will lead to acceptable solutions for text and icon labeling within a limited processing time.     

Establishing Classification and Hierarchy in Populated Place Labeling for Multiscale Mapping for The National Map

Current standards for federal mapping call for use of the Geographic Names Information System (GNIS) point layer for placement of United States populated place labels. However, this point layer contains limited classification information and hierarchy information, resulting in problems of map quality for database-driven, multi-scale, reference mapping, such as maps served by The National Map Viewer from USGS. Database-driven mapping often relies simply on what labels fit best in the map frame. Our research investigates alternative sources for labeling populated places, including polygons defined by the U.S. Census Bureau, such as incorporated place, census designated place (CDP), and economic place. Within each of these polygon layers we investigate relevant attributes from the decennial and economic censuses, such as population for incorporated places and CDPs, and the number of employees for economic places. The data selected are available for the entire country to serve national mapping requirements. This combination of data allows a more refined classification of populated places on maps that better represents relative importance. Visual importance on maps through scale should derive from more than simply residential population, but also economic importance, though comparison is made to this simpler case. We differentiate a fourth category of GNIS populated place points, essentially “neighborhoods” and related features—which are not incorporated places, CDPs, nor economic places. Populated places in this fourth class do not have federally defined boundaries, necessitating an alternative method for determining hierarchy in label presentation through scale.     

Generalization Services on the Web—Classification and an Initial Prototype Implementation

Much progress has been made in the field of web-based cartography through standards developed by the Open Geospatial Consortium (OGC). While automated access and presentation of cartographic data have been defined, the services for automated generalization are yet to be standardized. This paper aims to show advantages of applying the service concept to generalization and suggests several classification schemas of generalization services at different levels of granularity. A detailed explanation of a real implemented Generalization Service is provided. We show how software developers can make their generalization functionality available as a service and how these services can be accessed dynamically. For the implementation, the open source Java Unified Mapping Platform (JUMP) was extended to work as a framework for generalization. Generalization services could be used in different application scenarios, for instance as a middleware component extending a web map service with adaptive zooming or as stand-alone services supporting the production of topographic maps by national mapping agencies. They may also allow the development of a common research platform, where researchers would have access to a common generalization framework.     

Distinctive Recent Maps

Successful implementation of an intelligent system for automated map generalization requires formalization of cartographic principles that are, in many cases, only intuitively understood. Formalizing these principles requires acquisition and re-expression in the form of semantic nets, frames, production rules, or similar formalization methods. The various techniques for cartographic knowledge acquisition have been discussed on a theoretical basis; however, little empirical research has been conducted. This paper reports on empirical acquisition of cartographic knowledge by reverse engineering; that is, on trying to recapitulate decisions made on published documents or maps. The work is based on a computer-assisted multi-scale inventory of the Austrian National Topographic Map Series. Queries of the relational database, within which inventory data are stored, lead to the formulation of prototype production rules for modifying map symbols during automatic scale changes. Components of map generalization expressed in such rules include the selection behavior of settlement, transportation, and hydrographic objects, and the degree of simplification of settlement domains and building clusters. The acquired cartographic knowledge reveals quantitative relations between map elements and the changes in these relations that occur with scale transition. These insights can guide subsequent knowledge refinement using other acquisition methods. This paper provides, in addition, a conceptual framework by which other topographic map series may be compared at multiple scales.     

World Maps With Natural Boundaries

The graticule of meridians and parallels is a largely artificial type of map boundary that can detract from the display of irregular features such as oceans. Such natural boundaries as shorelines may be used instead as the boundary of world maps. The principle of natural boundaries has been applied to several examples of equal-area or conformal world ocean maps with single or multiple lobes. By careful selection of the poles and centers, these maps can show both oceans and continents in their entirety on a single map.