Spatiotemporal filtering of regional GNSS network’s position time series with missing data using principle component analysis

The existing spatiotemporal analysis methods suppose that the involved time series are complete and have the same data interval. However missing data inevitably occur in the position time series of Global Navigation Satellite Systems networks for many reasons. In this paper, we develop a modified principal component analysis to extract the Common Mode Error (CME) from the incomplete position time series. The principle of the proposed method is that a time series can be reproduced from its principle components. The method is equivalent to the method of Dong et al. (J Geophys Res 111:3405–3421, 2006) in case of no missing data in the time series and to the extended ‘stacking’ approach under the assumption of a uniformly spatial response. The new method is first applied to extract the CME from the position time series of the Crustal Movement Observation Network of China (CMONOC) over the period of 1999–2009 where the missing data occur in all stations with the different gaps. The results show that the CMEs are significant in CMONOC. The size of the first principle components for the North, East and Up coordinates are as large as 40, 41 and 37 % of total principle components and their spatial responses are not uniform. The minimum amplitudes of the first eigenvectors are only 41, 15 and 29 % for the North, East and Up coordinate components, respectively. The extracted CMEs of our method are close to the data filling method, and the Root Mean Squared error (RMS) values computed from the differences of maximum CMEs between two methods are only 0.31, 0.52 and 1.55 mm for North, East and Up coordinates, respectively. The RMS of the position time series is greatly reduced after filtering out the CMEs. The accuracies of the reconstructed missing data using the two methods are also comparable. To further comprehensively test the efficiency of our method, the repeated experiments are then carried out by randomly deleting different percentages of data at some stations. The results show that the CMEs can be extracted with high accuracy at the non missing-data epochs. And at the missing-data epochs, the accuracy of extracted CMEs has a strong dependence on the number of stations with missing data.     

嫦娥三号导航相机测图能力分析及地形重建

基于立体相机成像模型并结合相机参数对嫦娥三号导航相机3维测图能力进行分析,利用摄影测量原理和误差传播定律对巡视器30 m范围内的DEM精度进行了理论分析,推导出导航相机立体影像获得的采样点精度公式,并绘制了DEM的平面精度图和高程精度图;同时使用多线程技术开发了基于导航相机立体影像的地形快速重建算法,利用多线程技术完成影像的特征匹配和密集匹配,并通过分块内插生成DEM。该技术应用于嫦娥三号任务中,有力地支持了嫦娥三号遥操作路径规划相关任务。     

Mapping Wetland Areas Using Landsat-Derived NDVI and LSWI: A Case Study of West Songnen Plain,Northeast China

Increasing interest in wetlands for environmental management requires an understanding of the location, spatial extent, and configuration of the resource. The National Wetlands Inventory is the most commonly used data source for this information. However, its accuracy is limited in some contexts, such as agricultural and forested wetlands. An large number of studies have mapped wetlands worldwide from the perspective of land use and land cover change. However, information on the actual wetland planting areas annually is limited, which greatly impacts ongoing research. In this case study of the West Songnen Plain, we developed a simple algorithm for the quick mapping of wetlands by utilizing their unique physical features, such as annual display of phenological land-cover change of exposed soils, shallow flooding water, and plants from multi-temporal Landsat images. Temporal variations of the Normalized Difference Vegetation Index (NDVI) and Land Surface Water Index (LSWI) derived from Landsat images in 2010 for wetlands at different growth stages were analyzed. Results show that during the ante-tillering phase, the NDVI value (above zero) is lower than the LSWI value of paddies because of flooding of shallow water; during the reproductive and ripening phases, the NDVI value is higher than the LSWI value (above zero); and during the post-harvest wetland planting phase, the NDVI value is still higher than the LSWI value, but the LSWI value is negative. Wetland areas can be detected using one or two images in the optimum time window. The algorithm based on the difference of NDVI and LSWI values derived from Landsat images was used to extract the actual wetland planting area. Validated alongside statistical data, the algorithm showed high accuracy. Therefore, this algorithm highlights the unique features of wetlands and can help in mapping the actual wetland area annually on a regional scale. Results further indicate that the new method has a classification accuracy of 92 %. In comparison, two traditional methods based on Landsat-7/ETM registered accuracy rates of only 83 % and 87 % respectively.     

改善MODIS BRDF产品热点效应的方法研究

核驱动的Ross Thick-LiSparse Reciprocal(RTLSR)双向反射函数(BRDF)模型已广泛地应用于MODIS等星载传感器的业务化产品处理中。但是,对于多年MODIS二向反射产品历史数据,如何基于RTMLSR模型发展一种简单有效快速的方法,进行热点效应的校正是一个迫切需要解决的问题。本文提出了一种简单有效的方法,不需要对观测数据重新反演,直接在现有MODIS二向反射产品的基础上进行热点校正,方便用户对历史MODIS二向反射产品的使用。该方法应用POLDER-3/BRDF数据库和部分经MODIS业务化算法筛选的反射率数据进行验证,并与RTLSR模型和RTMLSR模型的结果进行比较,结果表明:(1)该方法比现有的MODIS业务化RTLSR算法,对热点反射率有明显改善,拟合相对误差平均降低了10.12%;(2)该方法相对于RTMLSR模型在热点反射率拟合效果上差别不大,相对误差相差2.10%;(3)该方法对热点和冷点归一化的植被指数(NDHD)的估算效果有一定程度的改善,相对于RTLSR模型降低了约4.99%,与RTMLSR模型的相对误差相差1.32%,该方法对直接应用现有MODIS BRDF产品,基于热点方向反射率反演植被结构参数(如植被聚集指数)的精度提高有现实应用价值。     

A new method for indentifying the model error of adjustment system

Some theory problems affecting parameter estimation are discussed in this paper. Influence and transformation between errors of stochastic and functional models is pointed out as well. For choosing the best adjustment model, a formula, which is different from the literatures existing methods, for estimating and identifying the model error, is proposed. On the basis of the proposed formula, an effective approach of selecting the best model of adjustment system is given.     

Shadow detection and compensation for color aerial images

A method for shadow detection and compensation for color aerial images is presented. It is considered that the intensity value of each image pixel is the product of illumination function and ground object reflection, and the shadowed regions on the image are mainly caused by the short of illumination, so the information compensation for the shadowed regions should concentrate on the illumination adjustment of concerned area on the basis of the analysis of whole image. The shadow detection and compensation procedure proposed by this paper consists of four steps.     

Automatic change detection for road networks from images based on GIS

Up to now, detailed strategies and algorithms of automatic change detection for road networks based on GIS have not been discussed. This paper discusses two different strategies of automatic change detection for images with low resolution and high resolution using old GIS data, and presents a buffer detection and tracing algorithm for detecting road from low-resolution images and a new profile tracing algorithm for detecting road from high-resolution images. For feature-level change detection (FL-CD), a so-called buffer detection algorithm is proposed to detect changes of features. Some ideas and algorithms of using GIS prior information and some context information such as substructures of road in high-resolution images to assist road detection and extraction are described in detail.     

The integrated spatial databases of GeoStar

GeoStar is the registered trademark of GIS software made by WTUSM in China. By means of the GeoStar, multi-scale images, DEMs, graphics and attributes integrated in very large seamless databases can be created, and the multi-dimensional dynamic visualization and information extraction are also available. This paper describes the fundamental characteristics of such huge integrated databases, for instance, the data models, database structures and the spatial index strategies. At last, the typical applications of GeoStar for a few pilot projects like the Shanghai CyberCity and the Guangdong provincial spatial data infrastructure (SDI) are illustrated and several concluding remarks are stressed.     

Correction on effect of Earth’s oblateness in inversion of GPS occultation data

By using observed CHAMP orbit ephemeredes and MSISE-90 dry air model and regarding the earth as a sphere and an ellipsoid respectively, phase delays are simulated and the simulated data are retrieved under different schemes. The comparison between the inverted temperature profiles and the model temperature profiles shows that by inverting observed data, we will get temperature results with large errors if the effect of Earth’s oblateness is omitted. The correction method is proved to be effective because the temperature errors decreased obviously with this method.     

Geoinformatics education in China

The paper gives an overview of the current status of education in geoinformatics in China. First, the paper provides a brief introduction to the history of geoinformatics education in China and a general review of the scientific and technological development. It then presents how the development affects the education and training in China. In the paper, universities and institutes in China that can award academic degrees related to geoinformatics are summarized, and undergraduate majors are briefly introduced. Next, the paper reports the work having been done by the national expert group on Surveying and Mapping, including the revision of discipline catalog and guide for graduate education and requirements. A list of typical curricula in geoinformatics education is suggested. Activities on promoting the graduate student exchange platform are presented. Finally, a case study of geoinformatics education in Wuhan University is discussed.