Cross-sensor calibration between Ikonos and Landsat ETM+ for spectral mixture analysis

Spectral mixture analysis is an algorithm that is developed to overcome the weakness in traditional land-use/land-cover (LULC) classification where each picture element (pixel) from remote sensing is assigned to one and only one LULC type. In reality, a remotely sensed signal from a pixel is often a spectral mixture from several LULC types. Spectral mixture analysis can derive subpixel proportions for the endmembers from remotely sensed data. However, one frequently faces the problem in determining the spectral signatures for the endmembers. This study provides a cross-sensor calibration algorithm that enables us to obtain the endmember signatures from an Ikonos multispectral image for spectral mixture analysis using Landsat ETM+ images. The calibration algorithm first converts the raw digital numbers from both sensors into at-satellite reflectance. Then, the Ikonos at-satellite reflectance image is degraded to match the spatial resolution of the Landsat ETM+ image. The histograms at the same spatial resolution from the two images are matched, and the signatures from the pure pixels in the Ikonos image are used as the endmember signatures. Validation of the spectral mixture analysis indicates that the simple algorithm works effectively. The algorithm is not limited to Ikonos and Landsat sensors. It is, in general, applicable to spectral mixture analysis where a high spatial resolution sensor and a low spatial resolution sensor with similar spectral resolutions are available as long as images collected by the two sensors are close in time over the same place.     

地物空间关系匹配的高分辨率遥感影像检索

从低层视觉特征与地物空间关系特征对影像内容进行描述,建立检索模板与目标影像间的相似性直方图表达,提出一种适用于高分辨率遥感影像检索的新方法。首先,利用Quin+树将大幅面原始遥感影像分解为一系列同尺寸的序列子块;然后,分别提取各子块的低层视觉特征与地物关系特征,并以子块为基元构建候选子块的特征直方图;最后,对比检索模板与候选子块间的特征直方图相似性,实现高分辨率遥感影像的检索。使用多幅多源高分辨率遥感影像进行实验,结果表明本文方法对耕地、水系、建筑物等地类的检索精度大都维持在0.8以上,且各项检索性能指标均优于已有的两种遥感图像检索算法。     

Full year crop monitoring and separability assessment with fully-polarimetric L-band UAVSAR: A case study in the Sacramento Valley,California

Spatial and temporal information on plant and soil conditions is needed urgently for monitoring of crop productivity. Remote sensing has been considered as an effective means for crop growth monitoring due to its timely updating and complete coverage. In this paper, we explored the potential of L-band fully-polarimetric Uninhabited Aerial Vehicle Synthetic Aperture Radar (UAVSAR) data for crop monitoring and classification. The study site was located in the Sacramento Valley, in California where the cropping system is relatively diverse. Full season polarimetric signatures, as well as scattering mechanisms, for several crops, including almond, walnut, alfalfa, winter wheat, corn, sunflower, and tomato, were analyzed with linear polarizations (HH, HV, and VV) and polarimetric decomposition (Cloude–Pottier and Freeman–Durden) parameters, respectively. The separability amongst crop types was assessed across a full calendar year based on both linear polarizations and decomposition parameters. The unique structure-related polarimetric signature of each crop was provided by multitemporal UAVSAR data with a fine temporal resolution. Permanent tree crops (almond and walnut) and alfalfa demonstrated stable radar backscattering values across the growing season, whereas winter wheat and summer crops (corn, sunflower, and tomato) presented drastically different patterns, with rapid increase from the emergence stage to the peak biomass stage, followed by a significant decrease during the senescence stage. In general, the polarimetric signature was heterogeneous during June and October, while homogeneous during March-to-May and July-to-August. The scattering mechanisms depend heavily upon crop type and phenological stage. The primary scattering mechanism for tree crops was volume scattering (>40%), while surface scattering (>40%) dominated for alfalfa and winter wheat, although double-bounce scattering (>30%) was notable for alfalfa during March-to-September. Surface scattering was also dominant (>40%) for summer crops across the growing season except for sunflower and tomato during June and corn during July-to-October when volume scattering (>40%) was the primary scattering mechanism. Crops were better discriminated with decomposition parameters than with linear polarizations, and the greatest separability occurred during the peak biomass stage (July-August). All crop types were completely separable from the others when simultaneously using UAVSAR data spanning the whole growing season. The results demonstrate the feasibility of L-band SAR for crop monitoring and classification, without the need for optical data, and should serve as a guideline for future research.     

基于MRF随机场和广义混合模型的遥感图像分级聚类

有限混合模型FM的分级聚类已广泛应用于不同领域,然而,它的计算复杂度与观测数据的平方成正比,因此,在海量数据方面的应用就受到了限制。另一方面,多光谱图像数据中同时包含有空间和光谱两类信息,但大多数基于像素的多光谱图像聚类方法,仅使用了其频谱信息而忽视了空间信息。本文提出了一种新的基于广义有限混合模型GFM的分级聚类方法,该算法把MRF随机场和GFM模型结合在一起,分类数可以通过PLIC准则自动确定。算法在执行过程中,采用K均值聚类方式获得过分类图像,分级聚类从过分类图像开始,代替原来从单点类开始的方式,这样可以方便获取GFM模型成分密度的初始参数。最后,采用由Gibbs采样器生成的仿真测试图对算法的精度进行了定量评价,通过与K均值聚类和FM聚类的比较说明了本文算法的优越性,同时用荷兰Flevoland农业地区的极化SAR图像验证了本文算法的有效性。     

Establishing robustness of a spatial dataset in a tolerance‐based vector model

Spatial data are usually described through a vector model in which geometries are represented by a set of coordinates embedded into an Euclidean space. The use of a finite representation, instead of the real numbers theoretically required, causes many robustness problems which are well known in the literature. Such problems are made even worse in a distributed context, where data is exchanged between different systems and several perturbations can be introduced in the data representation. In order to discuss the robustness of a spatial dataset, two implementation models have to be distinguished: the identity and the tolerance model. The robustness of a dataset in the identity model has been widely discussed in the literature and some algorithms of the Snap Rounding (SR) family can be successfully applied in such contexts. Conversely, this problem has been less explored in the tolerance model. The aim of this article is to propose an algorithm inspired by those of the SR family for establishing or restoring the robustness of a vector dataset in the tolerance model. The main ideas are to introduce an additional operation which spreads instead of snapping geometries, in order to preserve the original relation between them, and to use a tolerance region for such an operation instead of a single snapping location. Finally, some experiments on real‐world datasets are presented, confirming how the proposed algorithm can establish the robustness of a dataset.     

矢栅一体化拓扑关系的度量描述研究

首先,借助传统的9交模型确定了矢栅一体化数据模型拓扑关系的定性表达;其次,利用重叠部分的栅格数目占各个目标的比率来确定交叠分量,同时根据栅格的统计量来描述两目标接近程度的邻近分量;最终,将定性的拓扑关系、用以度量化描述的交叠分量及邻近分量以三元组的形式来描述目标间的拓扑关系,从而更加有效地实现了多类型地理目标拓扑关系的度量化描述。     

Spectral unmixing based spatiotemporal downscaling fusion approach

Time-series remote sensing data are important in monitoring land surface dynamics. Due to technical limitations, satellite sensors have a trade-off between temporal, spatial and spectral resolutions when acquiring remote sensing images. In order to obtain remote sensing images with high spatial resolution and high temporal frequency, spatiotemporal fusion methods have been developed. In this paper, we propose a Linear Spectral Unmixing-based Spatiotemporal Data Fusion Model (LSUSDFM) for spatial and temporal data fusion. In this model, the endmember abundance of the low-resolution image pixel is calculated based on that of the high-resolution image by the spectral mixture analysis. The endmember spectrum signals of low-resolution images are then calculated continuously within an optimized moving window. Subsequently, the fused image is reconstructed according to the endmember spectrum and its corresponding abundance map. A simulated dataset and real satellite images are used to test the fusion model, and the fusion results are compared with a current spectral unmixing based downscaling fusion model (SUDFM). Our experimental work shows that, compared to the SUDFM, the proposed LSUSDFM can achieve better quality and accuracy of fused images, especially in effectively eliminating the “plaque” phenomenon in the results by the SUDFM. The LSUSDFM has great potential in generating images with both high spatial resolution and high temporal frequency, as well as increasing the number of spectral bands of the high spatial resolution data.     

基于MRF随机场的多光谱遥感影像最优化分级聚类

有限混合模型FM的分级聚类已广泛应用于不同领域,然而,由于它的计算复杂度与观测数据量平方成正比,致使在遥感影像方面应用受到了限制。另外,多光谱图像能提供空间和光谱两类信息详细的数据,但是,大多数多光谱图像聚类方法是基于像素的聚类,仅使用了其光谱信息而忽视了空间信息。本文定义一个相对混合密度函数,通过引入一个q-参数来调节各成分密度对其混合分布的贡献,提出一种广义有限混合模型GFM．设计一种新的适用于多光谱遥感影像的GFM分级聚类算法。该算法把MRF随机场和GFM模型结合在了一起,分类数通过PLIC准则自动确定。最后,利用仿真结果验证该算法的有效性,同时通过与K均值聚类、FM分级聚类以及SVMM分级聚类的比较说明本文算法的优越性。     

Development of a spatial decision support system for monitoring earthquake-induced landslides based on aerial photographs and the finite element method

The evaluation of slope stability is essential for the management of landslide hazards. The integration of spatial information and geomechanical modeling facilitates the understanding and evaluation of landslide hazards. In this study, we use a spatial decision support system (SDSS)—incorporating aerial photographic data, GIS techniques, field investigations, and finite element geomechanical modeling—to analyze the mechanisms of the Hungtaiping landslide, which was induced by the 1999 Chi-Chi earthquake. The analysis clarifies the slide mechanisms that cannot be revealed either by examining aerial photographic or underground exploration data alone. The finite element modeling calibrated using digital aerial photographic data shows that the landslide results from the deformation and slides of the thick colluvium. Surficial displacements in the twenties of meters are attributed to the slide between the colluvium and the bedrock as well as the shear deformation and slides within the colluvium. The landslide SDSS can help determine model parameters, evaluate slide mechanisms and remediation measures, and predict slope behavior for a subsequent earthquake event.     

城市三维GIS中的海量数据组织方法

针对城市三维GIS中的海量数据组织需求,本文设计并实现了基于关系数据库的三维模型数据库,主要特点是支持地形数据、三维模型几何与纹理数据的一体化存储和属性的自定义动态管理。结合城市三维GIS中的海量数据访问和管理需求,从空间索引、纹理缓存、属性存取、数据加密、数据更新等方面介绍了该架构,指出按模型类别建立分类索引,能够提高空间索引查询效率。     

区域Gamma混合模型的SAR图像分割

针对传统Gamma混合模型用于SAR图像分割时忽略像素间空间相关性,导致分割结果不连续并产生大量误分割的现象,提出了区域Gamma混合模型的SAR图像分割算法。首先对图像进行分水岭分割,得到过分割区域块,然后将其作为输入样本进行基于Gamma混合模型的聚类,在模型的参数估计过程中进一步考虑区域间的空间相关性,设计邻域因子融入到迭代过程,得到邻域加权类分布概率。该算法充分利用像素间的空间相关性,能够降低噪声对分割结果的影响。通过合成图像和真实SAR图像的实验表明,本文算法能够实现SAR图像的准确分割。     

基于网格模型的导航道路图渐进式化简方法

针对将导航数据精度降低时的道路图形综合转换为基于网格模型的线图形化简问题,提出了一种基于网格模型的导航道路数据渐进式协同图形综合方法。该方法能保证综合后的单条道路和道路网都与原始高精度数据保持空间方向关系相似。在数据处理过程中,用导航需求来约束道路的渐进式图形简化,并对道路形状和道路交叉处的空间方向关系进行维护,顾及了导航中道路的可视化效果。已用实际生产中的导航数据进行了验证,证明了本文方法是有效且实用的。     

地勘工程3维空间数据模型及其数据结构设计

3维现象随着研究领域的不同，描述空间实体的方法存在较大差异，不可能设计一种数据模型来适合所有的应用领域，应根据研究领域空间实体分布特征，设计出专用的3维空间数据模型。本文对传统的3维数据模型进行了分析，以地质勘探工程的各种3维空间对象为研究背景，提出了适合地质勘探工程的矢量与栅格集成的面向对象混合数据模型，以复杂体－体－面－线－点对象之间的逻辑关系来建立对象之间的拓扑关系，设计了部分对象的数据结构，并以一玢岩型铁矿床为例，对数据结构进行了描述。     

Animation of Reliability in Computer-generated Dot Maps and Elevation Models

Visualization of reliability in spatial data has been the subject of considerable recent research activity. Animation has been suggested as one method to achieve this, and its application to various measures associated with class-area maps (classified satellite images and soil maps) has been discussed elsewhere. Animation is achieved by randomly selecting a location and then assigning it to a different map unit (cover or soil type) according to the information on the reliability associated with the original map or map units. In this article the same basic method is extended to mapping locational reliability in dot maps and surface error in a digital elevation model (DEM). In the former case, the dots, which are randomly located in the first place, are randomly relocated so that any meaningless positional information implicit in the location of the dots is lost while the meaningful information (the number of dots within a region) remains constant. In the DEM, animation uses a random field as an error surface, based upon the root mean squared error (RMSE). The amount of error at a location is constantly changed, giving no impression that the elevation is precisely known. The ability to vary the spatial autocorrelation within the error field provides a graphic illustration that the usual RMSE is not a sufficient method for the reporting of error in spatial databases. In both examples animation of reliability is believed to make a novel, but appropriate, use of the computer in cartography.     

A 3-dimensional data model for visualizing cloverleaf junction in a city model

The research work has been seldom done about cloverleaf junction expression in a 3-dimensional city model (3DCM). The main reason is that the cloverleaf junction is often in a complex and enormous construction. Its main body is bestraddle in air, and has aerial intersections between its parts. This complex feature made cloverleaf junction quite different from buildings and terrain, therefore, it is difficult to express this kind of spatial objects in the same way as for buildings and terrain. In this paper, authors analyze spatial characteristics of cloverleaf junction, propose an all-constraint points TIN algorithm to partition cloverleaf junction road surface, and develop a method to visualize cloverleaf junction road, surface using TIN. In order to manage cloverleaf junction data efficiently, the authors also analyzed the mechanism of 3DCM data management, extended BLOB type in relational database, and combined R-tree index to manage 3D spatial data. Based on this extension, an appropriate data structure for cloverleaf junction in 3DCM is proposed.     

树模型的三维可视化研究

树的表达无论是在3DGIS，VRGIS还是在3DCM中，都是一个难点问题。与一般实体对象不同的是，树明显具有分形空间特征，因此用常规简单的建模方法很难表达一颗树。现在分形算法常用来生成树模型，具有逼真的可视效果，但所产生的数据量对于GIS环境来说太大，可视化困难；另外，用树的图像替代树模型也是一种常用的方法，但在近视点视觉效果差，同时也没有表达出树木所固有的空间特征。本文作者提出一种方法，使用LOD技术，将树的分形模型与其图像替代模型两者结合起来，初步解决树模型在GIS中的可视化问题，可达到逼真、高速的可视效果。     

空间数据处理模型误差和不确定性分析

在GIS应用中,涉及到大量的模型应用,这些模型包括了利用GIS进行空间信息处理的大部分阶段中所用到的模型。模型处理以及分析结果往往是进行下一步应用的基础,因此模型处理结果的误差和不确定性制约了实际的GIS应用。影响空间数据处理模型的误差和不确定性的因素主要包括:定位和特征信息,制图,空间分析,空间数据库以及空间数据处理模型等所具有的误差和不确定性。主要分析了空间数据处理模型误差和不确定性的表达、来源以及分析方法。     

基于遥感与碳循环过程模型估算土壤有机碳储量

土壤有机碳含量的估算是当前全球碳循环研究的热点之一,但不同学者之间的估算值差异较大。从估算方法看,主要有基于土壤剖面的直接估算法和基于生态系统碳循环过程模型的间接估算法,这两种方法各有优缺点。直接估算法由于只反映了不同土壤或植被类型的土壤有机碳含量平均值的差异,因而空间分辨率较低。而间接估算法由于参数的简化,影响了土壤有机碳估算的空间精度。作者将遥感的高时空分辨率特征、反映生态系统碳循环动态变化的过程模型、实际测量的土壤有机碳结合起来,以求提高土壤有机碳估算的空间分辨率。考虑到受温度、水分的影响,土壤呼吸与土壤有机碳含量的关系并不好,而土壤基础呼吸由于剔除了温度和水分的影响,从而使其与土壤有机碳的关系非常密切,其测定系数R^2可达0.78。采用了结合遥感和碳循环过程的CASA模型及Van＇t Hoff土壤呼吸模型,首先估算了8km分辨率的土壤基础呼吸的空间分布,在此基础上结合实测的土壤有机碳估算了8km分辨率的土壤有机碳的空间分布。     

空间数据处理模型误差和不确定性分析

在GIS应用中,涉及到大量的模型应用,这些模型包括了利用GIS进行空间信息处理的大部分阶段中所用到的模型.模型处理以及分析结果往往是进行下一步应用的基础,因此模型处理结果的误差和不确定性制约了实际的GIS应用.影响空间数据处理模型的误差和不确定性的因素主要包括:定位和特征信息,制图,空间分析,空间数据库以及空间数据处理模型等所具有的误差和不确定性.主要分析了空间数据处理模型误差和不确定性的表达、来源以及分析方法.     

Automatic interpretation of digital maps

In the past, the availability and/or the acquisition of spatial data were often the main problems of the realization of spatial applications. Meanwhile this situation has changed: on one hand, comprehensive spatial datasets already exist and on the other hand, new sensor technologies have the ability to capture fast and with high quality large amounts of spatial data. More and more responsible for the increasing accessibility of spatial data are also collaborative mapping techniques which enable users to create maps by themselves and to make them available in the internet. However, the potential of this diversity of spatial data can only hardly be utilized. Especially maps in the internet are represented very often only with graphical elements and no explicit information about the map’s scale, extension and content is available. Nevertheless, humans are able to extract this information and to interpret maps. For example, it is possible for a human to distinguish between rural and industrial areas only by looking at the objects’ geometries. Furthermore, a human can easily identify and group map objects that belong together. Also the type, scale and extension of a map can be identified under certain conditions only by looking at the objects’ geometries. All these examples can be subsumed under the term “map interpretation”. In this paper it is discussed how map interpretation can be automated and how automatic map interpretation can be used in order to support other processes. The different kinds of automatic map interpretation are discussed and two approaches are shown in detail.