一种适用于极区的横轴墨卡托海图

极区位于高纬度地区,极区航海图投影采用墨卡托投影时,图上极区附近变形较大,在极点处变形达到无穷大;采用极球面投影时,图上恒向线不是直线;心射切面投影不是正形投影,不能直接在图上量取方位和距离;兰伯特正形投影图上,经度是由角度表示,而不是真的经度;采用横轴墨卡托投影时,几乎所有的经线以及纬线在图上为曲线。以上常用的海图投影不利于船舶在极区安全航行,也不利于航海人员进行海图作业。针对此问题,将构建的伪经纬线网与经纬线网进行转换,通过Matlab得到基于伪经纬线网的横轴墨卡托投影图,此图具有类似于墨卡托海图的特点,图上极区变形失真小,图上伪经纬线的特点可方便极区海图作业。     

1949年以来兰州城市资本密度空间变化及其机制

构建了计划经济体制和中国转型期城市资本密度空间分布理论模型,并以兰州市为例进行了实证.利用高精度卫星影象和大规模实地调查方法获取数据,采用GIS手段和数理统计方法,利用建筑高度指标,分区、分行业对兰州城市资本密度空间变化进行了深入研究.同时,采用资本存量增量(或增长率)/建成区面积增量(或增长率)的指标衡量资本角度的城市紧凑化程度,研究了城市紧凑化变动过程.结果表明:①中国计划经济体制下的城市资本空间密度为缓慢波动下降或均衡分布规律,甚至出现上升趋势.转轨期(中国)城市资本空间密度由中心到外围的相应规律则居于两者之间,存在向市场经济空间模型靠近的趋势;②中国城市计划经济体制下的用地空间由中心到外围呈现商务→住宅→工业→农业的模式,而转轨期依然基本保留了此特征;③转轨期的建筑高度提高速度存在加速趋势,住宅、商业、办公等类建筑高度明显高于计划时期,其空间分布模型也更接近市场经济体制模型;④从资本视角审视,中国城市将更加紧凑化,且随时间的推移,城市资本密度空间变化曲线大致存在较为明显的"雁行波动上升式"规律.     

基于地图方里网的全球剖分系统

针对现有剖分模型的不足,提出了一种基于地图方里网的全球剖分系统,有效避免了传统经纬度格网模型在高纬度地区的形状退化和正多面体格网模型的面片形状不规则问题。制定了相应的编码,实现了面片编码与传统地理坐标之间的转换和邻接关系的计算,最后对方里网在地球椭球面上的变形规律进行了研究。研究结果表明:基于地图方里网的剖分系统不仅具有科学的数据组织形式,而且面片单元的变形面积小、变形规律稳定。     

A compromise aspect-adaptive cylindrical projection for world maps

There are two problems with current cylindrical projections for world maps. First, existing cylindrical map projections have a static height-to-width aspect ratio and do not automatically adjust their aspect ratio in order to optimally use available canvas space. Second, many of the commonly used cylindrical compromise projections show areas and shapes at higher latitudes with considerable distortion. This article introduces a new compromise cylindrical map projection that adjusts the distribution of parallels to the aspect ratio of a canvas. The goal of designing this projection was to show land masses at central latitudes with a visually balanced appearance similar to how they appear on a globe. The projection was constructed using a visual design procedure where a series of graphically optimized projections was defined for a select number of aspect ratios. The visually designed projections were approximated by polynomial expressions that define a cylindrical projection for any height-to-width ratio between 0.3:1 and 1:1. The resulting equations for converting spherical to Cartesian coordinates require a small number of coefficients and are fast to execute. The presented aspect-adaptive cylindrical projection is well suited for digital maps embedded in web pages with responsive web design, as well as GIS applications where the size of the map canvas is unknown a priori. We highlight the projection with a height-to-width ratio of 0.6:1, which we call the Compact Miller projection because it is inspired by the Miller Cylindrical projection. Unlike the Miller Cylindrical projection, the Compact Miller projection has a smaller height-to-width ratio and shows the world with less areal distortion at higher latitudes. A user study with 448 participants verified that the Compact Miller – together with the Plate Carrée projection – is the most preferred cylindrical compromise projection.     

A new pseudocylindrical equal-area projection for adaptive composite map projections

The recently introduced adaptive composite map projection technique changes the projection to the geographic area shown on a map. It is meant as a replacement for the commonly used web Mercator projection, which grossly distorts areas when representing the entire world. The original equal-area version of the adaptive composite map projection technique uses the Lambert azimuthal projection for regional maps and three alternative projections for world maps. Adaptive composite map projections can include a variety of other equal-area projections when the transformation between the Lambert azimuthal and the world projections uses Wagner’s method. To select the most suitable pseudocylindrical projection, the distortion characteristics of a pseudocylindrical projection family are analyzed, and a user study among experts in the area of map projections is carried out. Based on the results of the distortion analysis and the user study, a new pseudocylindrical projection is recommended for extending adaptive composite map projections. The new projection is equal-area throughout the transformation to the Lambert azimuthal projection and has better distortion characteristics then small-scale projections currently included in the adaptive composite map projection technique.     

GIS支持下的东川区泥石流危险度区划

在建立东川区泥石流信息系统的基础上,根据泥石流形成原因,初步选择面积、相对切割度、主沟沟床比降、相对高差、主沟长度、圆状率作为评价指标,利用东川区21条典型泥石流沟的数据,通过主成分分析,选择出面积、圆状率、相对高差三个指标,由主成分载荷得到危险度评价公式,在ARC／INFO和ARCVIEW软件的支持下,把东川区划分为60个评价单元,计算出各个单元内的危险度总分值,依据分值划分为四个级别,重度危险区、危险区、轻度危险区和基本无危险区,最后根据泥石流危险度区划原则,得到东川区泥石流危险度区划图。     

Computation and visualisation of the accuracy of old maps using differential distortion analysis

The accuracy of old maps can hold interesting historical information, and is therefore studied using distortion analysis methods. These methods start from a set of ground control points that are identified both on the old map and on a modern reference map or globe, and conclude with techniques that compute and visualise distortion. Such techniques have advanced over the years, but leave room for improvement, as the current ones result in approximate values and a coarse spatial resolution. We propose a more elegant and more accurate way to compute distortion of old maps by translating the technique of differential distortion analysis, used in map projection theory, to the setting where an old map and a reference map are directly compared. This enables the application of various useful distortion metrics to the study of old maps, such as the area scale factor, the maximum angular distortion and the Tissot indicatrices. As such a technique is always embedded in a full distortion analysis method we start by putting forward an optimal analysis method for a general-purpose study, which then serves as the foundation for the development of our technique. Thereto, we discuss the structure of distortion analysis methods and the various options available for every step of the process, including the different settings in which the old map can be compared to its modern counterpart, the techniques that can be used to interpolate between both, and the techniques available to compute and visualise the distortion. We conclude by applying our general-purpose method, including the differential distortion analysis technique, to an example map also used in other literature.     

Graphical design of world map projections

The design of new map projections has up until now required mathematical and cartographic expertise that has limited this activity to a small group of specialists. This article introduces the background mathematics for a software-based method that enables cartographers to easily design new small-scale world map projections. The software is usable even by those without mathematical expertise. A new projection is designed interactively in an iterative process that allows the designer to graphically and numerically assess the graticule, the representation of the continents, and the distortion properties of the new projection. The method has been implemented in Flex Projector, a free and open-source application enabling users to quickly create new map projections and modify existing projections. We also introduce new tools that help evaluate the distortion properties of projections, namely a configurable acceptance index to assess areal and angular distortion, a derived acceptance visualization, and interactive profiles through the distortion space of a projection. To illustrate the proposed method, a new projection, the Cropped Ginzburg VIII projection, is presented.     

A quantitative assessment of vegetable farming on vacant lots in an urban fringe area in Metro Manila: Can it sustain long-term local vegetable demand?

The aim of this study was to quantitatively assess vegetable production and consumption in an urban fringe area in Metro Manila, the Philippines. We conducted field investigations into vegetable production at a farm within a subdivision, monitored customers for 1 week at the farm's vegetable stand, interviewed customers of the stand, and conducted household monitoring and interviews to identify purchasing behavior. We also conducted an object-based WorldView-2 image analysis of a larger area to identify current and potential future vegetable farming plots and a spatial analysis using GIS to visualize the vegetable production/consumption ratio. We found that about one-third of current vegetable demand could be supplied by greater promotion of vegetable farming in vacant lots and by setting up more vegetable stands at spatial intervals less than 1-km. Subdivision regulations also influenced vegetable farming because of the level of social acceptance of caretakers of vacant lots. The current planning and institutional framework does not consider social contexts in land use and planning. We suggested that a new framework is necessary; it should consider balancing the formation of social capital and formalize planning and institutional procedures to support sustainable local vegetable production.     

Density inhomogeneities of the upper mantle of the Central Balkans

A 3-D density model was created for the Central Balkans area down to a depth of 670  km on the basis of seismic (both artificial sources and earthquakes) and gravity data. This model is based on density columns constructed for the main geological units of the study region. The densities for these columns were obtained using the density variation method. This method makes it possible to extrapolate the density distribution from the well-studied uppermost layers to the deeper levels of the Earth. The constructed 3-D density model was interpreted in relation to the available data on the heat flow and the seismicity of the region. The subdivision of the region by the Maritza fault into two parts—the southern part including the Rhodope massif and the northern part including the structures of Alpine activation of Srednogorie and the Balkanides—was confirmed. The upraised position of the 400  km boundary within the upper mantle, which was established from the density modelling, is assumed to be a sign of development of recent geodynamical processes over the Srednogorie block. From the viewpoint of seismicity prediction, a finding of mantle inhomogeneities orthogonal to the Maritza suture is of great importance.     

Novel shape indices for vector landscape pattern analysis

The formation of an anisotropic landscape is influenced by natural and/or human processes, which can then be inferred on the basis of geometric indices. In this study, two minimal bounding rectangles in consideration of the principles of mechanics (i.e. minimal width bounding (MWB) box and moment bounding (MB) box) were introduced. Based on these boxes, four novel shape indices, namely MBLW (the length-to-width ratio of MB box), PAMBA (area ratio between patch and MB box), PPMBP (perimeter ratio between patch and MB box) and ODI (orientation difference index between MB and MWB boxes), were introduced to capture multiple aspects of landscape features including patch elongation, patch compactness, patch roughness and patch symmetry. Landscape pattern was, thus, quantified by considering both patch directionality and patch shape simultaneously, which is especially suitable for anisotropic landscape analysis. The effectiveness of the new indices were tested with real landscape data consisting of three kinds of saline soil patches (i.e. the elongated shaped slightly saline soil class, the circular or half-moon shaped moderately saline soil, and the large and complex severely saline soil patches). The resulting classification was found to be more accurate and robust than that based on traditional shape complexity indices.     

西南丘陵区土地整理的预期景观效应分析——以重庆市江北区五宝镇大树村为例

在GIS技术的支持下,选取了不同坡度早地的斑面积占项目区实施块总面积、斑块平均面积、斑块平均周长、斑块个数、斑块平均密度、斑块总面积的比例等指标分析了土地整理前后斑块的变化;运用廊道长度、宽度、密度、连通性、环度、各类型廊道长度占项目区廊道总长度的比例等指标分析了整理前后廊道的变化特征。研究结果表明：土地整理后,斑块的平均面积、平均周长增大,斑块个数减少和斑块平均密度变小,形状更加规范化;耕地景观破碎度降低,廊道布局更加合理;各类型廊道结构更加完善,功能更加齐全,各类型廊道的连通性和环度明显提高。     

Storage-efficient method for generating contours focusing on roundness

A storage-efficient contour generation method, focusing on planar contours, is developed. Given cartographic elevations on a rectangular lattice, a continuous bivariate function, z = f(x, y), is determined by interpolating the elevation values. Then, we focus on a contour determined by z = constant. The contour curve is partitioned into multiple sections, each of which is exactly or approximately round. Three curvature types are introduced to evaluate the roundness of each section. The area and perimeter of the contour are computed by one-dimensional line integration using Green’s theorem. If the contour is open, it is divided into two curves starting from the same initial point, with the control points advancing in opposite directions. Two types of numerical experiments are performed to validate the effectiveness of the proposed method. One experiment uses an analytically defined elevation function and investigates the number of control points and computation time for a resulting computation error. The second experiment uses actual digital elevation model data of an isolated island in Japan and compares the proposed method with existing ones. Because the algorithm does not require lattice subdivision and the number of control points is drastically reduced, the proposed method is storage efficient.     

Quasi indicatrix approach for distortion visualization and analysis for map projections

Tissot’s indicatrix or ellipse of distortion is a diagram that is the projection of an infinitesimal circle on the original surface. It is normally an ellipse of which elongation depends on the amount of distortion caused by map projection. It provides a medium for analyzing existing projections and developing new ones. The ellipse can be scaled and depicted on the map for visualization purposes. This paper presents an alternative approach, in which the projection of a finite small circle on the sphere is used. Its projection is normally an ellipse that can be very close to Tissot’s indicatrix, and is called quasi indicatrix, here. Its parameters can be derived from the forward projection equations without using partial derivatives. Therefore, it is a useful and practical approach from a programmer’s point of view. The quasi indicatrix approach is also numerically tested on Aitoff–Hammer projection with a set of points. The indicatrix parameters obtained by using this approach deviate 0.5% from the ground truths at most, being the average less than 0.2%.     

Equal arc ratio projection and a new spherical triangle quadtree model

A projection family, named Equal Arc Ratio Projection (EARP), is proposed in this paper, which has two kinds of geometric cases. Coordinates in the projections are determined by the ratios between four special arcs. The projection for spherical octahedron is derived, which is available for Dutton's Quaternary Triangular Mesh (QTM), and both projections of Goodchild's and Otoo's for QTM are found to be two special implementations of EARP. The icosahedral projection and thus the corresponding spherical triangle quadtree partition model are also developed. The partition model is compared with several other global grid models.     

Influence of repeated freeze-thaw on dynamic modulus and damping ratio properties of silty sand

Under repeated freezing and thawing in deep seasonal frozen regions, the stability and strength of the soil are imposed in the form of large uneven settlement, instability and strength reduction, which...     

Gauss Krüger projection for areas of wide longitudinal extent

This paper describes a new algorithm for calculating all equations (direct and inverse transformation, convergence of meridians, linear distortion, calculus of surfaces, and arc‐to‐chord correction) for the Gauss–Krüger projection. Instead of using different equations for each problem, all the calculi are based on the equations used to obtain the direct transformation. These equations are also more accurate than previous ones and can be extended to an arbitrary width. This paper also explains how the Gauss–Krüger projection may satisfy the needs and requirements of civil engineers as well as Geographical Information Systems (GIS) users, particularly if extended beyond its usual 4° longitudinal range, and therefore can be used as a global reference system for GIS.     

张家界旅游企业空间成长对城镇空间形态的影响

运用核密度估计法和标准差椭圆分析法,分析张家界旅游企业空间成长特征及演化过程;选取紧凑度指标,结合城镇空间均衡分布椭圆,分析张家界城镇空间形态演化过程;通过比较不同阶段旅游企业椭圆和城镇空间椭圆的重合面积,剖析旅游企业空间成长对城镇空间形态的影响程度。研究发现：① 张家界旅游企业空间成长过程和特征为：由出现期的“一核多极点”逐渐演变为生成期的“核心-边缘”,最后在发展期形成一个高集聚中心和三个次中心联动的“组团扩展型”格局;② 张家界建成区面积呈现增长态势,城镇空间不断扩展,建成区用地紧凑度逐渐下降;③ 不同阶段,张家界旅游企业空间成长对城镇空间形态的影响程度有一定的差异性,整体呈现出逐渐增强的趋势。研究结论对于重新认识旅游对城镇化建设的作用具有重要的现实意义。     

喀斯特洼地的分形特性研究*

以贵州高原中部普定马官地区和云南高原边缘西畴峰丛-洼地区为试验区,研究和探讨了喀斯特洼地发育过程中的分形特性。     

SRG剖分编码与地理坐标的转换算法研究

提出一种基于SRG剖分方法的剖分编码与经纬度坐标之间的转换算法,不涉及任何投影变换,计算过程只需加、减、乘、除简单算术运算,计算速度快,且SRG剖分编码本身就具有固定的方向性,利于邻近搜索。算法中采用坐标系辅助区分一些难以区分的菱形块,降低了转化误码率,进而提高了转换的精度。