On the spherical and spheroidal harmonic expansion of the gravitational potential of the topographic masses

This paper is devoted to the spherical and spheroidal harmonic expansion of the gravitational potential of the topographic masses in the most rigorous way. Such an expansion can be used to compute gravimetric topographic effects for geodetic and geophysical applications. It can also be used to augment a global gravity model to a much higher resolution of the gravitational potential of the topography. A formulation for a spherical harmonic expansion is developed without the spherical approximation. Then, formulas for the spheroidal harmonic expansion are derived. For the latter, Legendre’s functions of the first and second kinds with imaginary variable are expanded in Laurent series. They are then scaled into two real power series of the second eccentricity of the reference ellipsoid. Using these series, formulas for computing the spheroidal harmonic coefficients are reduced to surface harmonic analysis. Two numerical examples are presented. The first is a spherical harmonic expansion to degree and order 2700 by taking advantage of existing software. It demonstrates that rigorous spherical harmonic expansion is possible, but the computed potential on the geoid shows noticeable error pattern at Polar Regions due to the downward continuation from the bounding sphere to the geoid. The second numerical example is the spheroidal expansion to degree and order 180 for the exterior space. The power series of the second eccentricity of the reference ellipsoid is truncated at the eighth order leading to omission errors of 25 nm (RMS) for land areas, with extreme values around 0.5 mm to geoid height. The results show that the ellipsoidal correction is 1.65 m (RMS) over land areas, with maximum value of 13.19 m in the Andes. It shows also that the correction resembles the topography closely, implying that the ellipsoidal correction is rich in all frequencies of the gravity field and not only long wavelength as it is commonly assumed.     

Spherical harmonic modelling to ultra-high degree of Bouguer and isostatic anomalies

The availability of high-resolution global digital elevation data sets has raised a growing interest in the feasibility of obtaining their spherical harmonic representation at matching resolution, and from there in the modelling of induced gravity perturbations. We have therefore estimated spherical Bouguer and Airy isostatic anomalies whose spherical harmonic models are derived from the Earth’s topography harmonic expansion. These spherical anomalies differ from the classical planar ones and may be used in the context of new applications. We succeeded in meeting a number of challenges to build spherical harmonic models with no theoretical limitation on the resolution. A specific algorithm was developed to enable the computation of associated Legendre functions to any degree and order. It was successfully tested up to degree 32,400. All analyses and syntheses were performed, in 64 bits arithmetic and with semi-empirical control of the significant terms to prevent from calculus underflows and overflows, according to IEEE limitations, also in preserving the speed of a specific regular grid processing scheme. Finally, the continuation from the reference ellipsoid’s surface to the Earth’s surface was performed by high-order Taylor expansion with all grids of required partial derivatives being computed in parallel. The main application was the production of a 1′ × 1′ equiangular global Bouguer anomaly grid which was computed by spherical harmonic analysis of the Earth’s topography–bathymetry ETOPO1 data set up to degree and order 10,800, taking into account the precise boundaries and densities of major lakes and inner seas, with their own altitude, polar caps with bedrock information, and land areas below sea level. The harmonic coefficients for each entity were derived by analyzing the corresponding ETOPO1 part, and free surface data when required, at one arc minute resolution. The following approximations were made: the land, ocean and ice cap gravity spherical harmonic coefficients were computed up to the third degree of the altitude, and the harmonics of the other, smaller parts up to the second degree. Their sum constitutes what we call ETOPG1, the Earth’s TOPography derived Gravity model at 1′ resolution (half-wavelength). The EGM2008 gravity field model and ETOPG1 were then used to rigorously compute 1′ × 1′ point values of surface gravity anomalies and disturbances, respectively, worldwide, at the real Earth’s surface, i.e. at the lower limit of the atmosphere. The disturbance grid is the most interesting product of this study and can be used in various contexts. The surface gravity anomaly grid is an accurate product associated with EGM2008 and ETOPO1, but its gravity information contents are those of EGM2008. Our method was validated by comparison with a direct numerical integration approach applied to a test area in Morocco–South of Spain (Kuhn, private communication 2011) and the agreement was satisfactory. Finally isostatic corrections according to the Airy model, but in spherical geometry, with harmonic coefficients derived from the sets of the ETOPO1 different parts, were computed with a uniform depth of compensation of 30?km. The new world Bouguer and isostatic gravity maps and grids here produced will be made available through the Commission for the Geological Map of the World. Since gravity values are those of the EGM2008 model, geophysical interpretation from these products should not be done for spatial scales below 5 arc minutes (half-wavelength).     

增强型空间像元分解时空遥感影像融合算法

高空间、高时间分辨率的遥感影像对地表与大气环境的实时精细监测具有重要作用,但单一卫星传感器获取的遥感影像存在空间与时间分辨率相互制约的问题,时空融合技术发展成为了低成本、高效生成满足不同应用需求的高时空分辨率遥感影像的有效手段.近年来,国内外学者提出了大量的时空融合算法,但对于复杂的地物类型变化的空间细节修复仍存在挑战...     

Generation of high spatial and temporal resolution NDVI and its application in crop biomass estimation

Abstract

While data like HJ-1 CCD images have advantageous spatial characteristics for describing crop properties, the temporal resolution of the data is rather low, which can be easily made worse by cloud contamination. In contrast, although Moderate Resolution Imaging Spectroradiometer (MODIS) can only achieve a spatial resolution of 250 m in its normalised difference vegetation index (NDVI) product, it has a high temporal resolution, covering the Earth up to multiple times per day. To combine the high spatial resolution and high temporal resolution of different data sources, a new method (Spatial and Temporal Adaptive Vegetation index Fusion Model [STAVFM]) for blending NDVI of different spatial and temporal resolutions to produce high spatial–temporal resolution NDVI datasets was developed based on Spatial and Temporal Adaptive Reflectance Fusion Model (STARFM). STAVFM defines a time window according to the temporal variation of crops, takes crop phenophase into consideration and improves the temporal weighting algorithm. The result showed that the new method can combine the temporal information of MODIS NDVI and spatial difference information of HJ-1 CCD NDVI to generate an NDVI dataset with both high spatial and high temporal resolution. An application of the generated NDVI dataset in crop biomass estimation was provided. An average absolute error of 17.2% was achieved. The estimated winter wheat biomass correlated well with observed biomass (R ² of 0.876). We conclude that the new dataset will improve the application of crop biomass estimation by describing the crop biomass accumulation in detail. There is potential to apply the approach in many other studies, including crop production estimation, crop growth monitoring and agricultural ecosystem carbon cycle research, which will contribute to the implementation of Digital Earth by describing land surface processes in detail.     

土地利用景观格局核心指数提取:以中国广州市为例

景观格局是研究景观功能和动态的基础。景观指数分析在土地利用/土地覆盖、生态系统服务、森林监控、城市蔓生以及生物多样性等方面应用广泛。现有许多景观指数之间存在不同程度的相关性,不仅导致信息冗余,还将影响结果解译的清晰性和准确性。同时,已有研究主要针对栅格数据,针对矢量数据的景观指数分析研究较少。为解决上述问题,以广州市土地利用景观格局为例,在矢量数据格式下计算了44个常用的景观指数,并运用相关分析和因子分析相结合的多元统计分析法,提取了6个核心景观指数,这些指数描述了广州市土地利用景观格局的6个维度:景观破碎度、平均斑块面积、平均形状复杂度、空间分散度、形状复杂度差异和土地类型的相似性。     

高分辨率遥感影像的景观格局分析

以景观生态学理论为研究基础,以高分辨率遥感影像等为数据源、利用RS和GIS技术,将梅州城市规划区的土地利用类型分为耕地、林地、草地、城镇用地及公路用地、湖泊水面、河流水面、裸地、滩涂等9类,并从景观的基本构成、景观水平的破碎化与多样性、不同景观斑块的空间格局特征等进行了研究分析;并提出对研究区合理规划土地开发利用的建议...     

中国陆地区域陆表温度业务化遥感反演算法及产品运行系统

地表温度反演的裂窗算法已成功应用于NOAA系列卫星热红外遥感数据。目前，裂窗算法中应用较为广泛的一种是Becker等人于1990年提出的局地裂窗算法，主要是通过辐射传输模型模拟不同地表条件和大气状况下，地表温度和发射率对红外辐射亮温的影响，从而发展出一个利用AVHRR4，5通道亮温数据反演地表温度的线性模型。在晴空无云和地表比辐射率能精确估算的情况下，Becker算法反演地表温度的精度在1K以内。Becker算法用Lowtran程序模拟计算地表辐射量，且模型中参数主要针对NOAA-9传感器特性得到。本文在Becker算法的基础上，针对NOAA-16／17传感器热红外通道光谱响应函数特性，利用最新的、计算光谱分辨率更高的MODTRAN程序模拟不同大气状况下，不同地表温度和发射率对NOAAAVHRR4，5通道辐射亮温响应特性的影响，改进Becker算法中模型参数，使之能适用于NOAA-16／17热红外数据。同时，本文利用植被指数NDVI，在中国陆地区域lkm分辨率最新地表分类数据的基础上，得到模型中需要的地表比辐射率参数，将改进的模型应用于1km分辨率NOAA17数据，得到了旬合成中国陆地区域范围地表温度，通过地面气象台站实测数据对比验证．取得了较好的结果。     

Comparison of the Landsat Surface Reflectance Climate Data Record (CDR) and manually atmospherically corrected data in a semi-arid European study area

This study contributes to the quality assessment of atmospherically corrected Landsat surface reflectance data that are routinely generated by the Landsat Ecosystem Disturbance Adaptive Processing System (LEDAPS). This dataset, named Landsat Surface Reflectance Climate Data Record (Landsat CDR), is available at global scale and offers unprecedented opportunities to land monitoring and management services that require atmospherically corrected Earth observation (EO) data. Our assessment is based on the comparison of the Landsat CDR data against a set of Landsat and DEIMOS-1 images processed to a high degree of accuracy using an industry-standard atmospheric correction algorithm (ATCOR-2). The software package has been used for many years and its correction procedures can be considered consolidated and well-established. The dataset of Landsat and DEIMOS-1 images was acquired over a semi-arid agricultural area located in Lower Austria and was independently corrected by using a manual fine-tuning of ATCOR-2 parameters to reach the highest possible accuracy. Results show a very good correspondence of the surface reflectance in each of the six reflective spectral channels as well as for the NDVI (Normalized Difference Vegetation Index). An additional comparison against a NDVI time series from MODIS revealed also a good correspondence. Coefficients of determination (R²) between the two multi-year and multi-seasonal Landsat/DEIMOS datasets range between 0.91 (blue band) and 0.98 (nIR, SWIR-1 and SWIR-2). The results obtained for our semi-arid test site in Austria confirm previous findings and suggest that automatic atmospheric procedures, such as the one implemented by LEDAPS are accurate enough to be used in land monitoring services that require consistent multi-temporal surface reflectance data.     

Using pan-sharpened high resolution satellite data to improve impervious surfaces estimation

Impervious surface is an important environmental and socio-economic indicator for numerous urban studies. While a large number of researches have been conducted to estimate the area and distribution of impervious surface from satellite data, the accuracy for impervious surface estimation (ISE) is insufficient due to high diversity of urban land cover types. This study evaluated the use of panchromatic (PAN) data in very high resolution satellite image for improving the accuracy of ISE by various pan-sharpening approaches, with a further comprehensive analysis of its scale effects. Three benchmark pan-sharpening approaches, Gram-Schmidt (GS), PANSHARP and principal component analysis (PCA) were applied to WorldView-2 in three spots of Hong Kong. The on-screen digitization were carried out based on Google Map and the results were viewed as referenced impervious surfaces. The referenced impervious surfaces and the ISE results were then re-scaled to various spatial resolutions to obtain the percentage of impervious surfaces. The correlation coefficient (CC) and root mean square error (RMSE) were adopted as the quantitative indicator to assess the accuracy. The accuracy differences between three research areas were further illustrated by the average local variance (ALV) which was used for landscape pattern analysis. The experimental results suggested that 1) three research regions have various landscape patterns; 2) ISE accuracy extracted from pan-sharpened data was better than ISE from original multispectral (MS) data; and 3) this improvement has a noticeable scale effects with various resolutions. The improvement was reduced slightly as the resolution became coarser.     

互联网+地表覆盖验证及应用

为有效地组织大范围地表覆盖数据验证,需要发动和组织来自不同地域、不同机构的专家,但往往缺少网络化环境下的协同验证方法和工具。针对这一问题,将互联网+、服务计算技术等与地表覆盖验证业务相融合,提出了互联网+地表覆盖验证服务模型,构建了面向互联网+的地表覆盖验证业务流程和方法,研制了基于景观形状指数（landscape shape index,LSI）抽样算法的在线验证系统GLCVal（global land cover validation）,为大范围、高分辨率地表覆盖数据验证提供了新思路、新方法。利用所提出的方法与系统,与GEO（Group of Earth Observation）、UN-GGIM（the United Nations Initiative on Global Geospatial Information Management）等组织合作开展了GlobeLand30的全球验证,有效支持了跨区域专家的主动参与、多角色协同,促进了验证资源联通和共享,提高了验证效率。     

Simulation of landscape spatial layout evolution in rural-urban fringe areas: a case study of Ganjingzi District

In recent years, the rapid expansion of urban spaces has accelerated the mutual evolution of landscape types. Analyzing and simulating spatio-temporal dynamic features of urban landscape can help to reveal its driving mechanisms and facilitate reasonable planning of urban land resources. The purpose of this study was to design a hybrid cellular automata model to simulate dynamic change in urban landscapes. The model consists of four parts: a geospatial partition, a Markov chain (MC), a multi-layer perceptron artificial neural network (MLP-ANN), and cellular automata (CA). This study employed multivariate land use data for the period 2000–2015 to conduct spatial clustering for the Ganjingzi District and to simulate landscape status evolution via a divisional composite cellular automaton model. During the period of 2000–2015, construction land and forest land areas in Ganjingzi District increased by 19.43% and 15.19%, respectively, whereas farmland, garden lands, and other land areas decreased by 43.42%, 52.14%, and 75.97%, respectively. Land use conversion potentials in different sub-regions show different characteristics in space. The overall land-change prediction accuracy for the subarea-composite model is 3% higher than that of the non-partitioned model, and misses are reduced by 3.1%. Therefore, by integrating geospatial zoning and the MLP-ANN hybrid method, the land type conversion rules of different zonings can be obtained, allowing for more effective simulations of future urban land use change. The hybrid cellular automata model developed here will provide a reference for urban planning and policy formulation.     

一种保持矢量数据几何形状的数字水印算法

针对目前矢量地理空间数据数字水印技术的需求,从矢量地理空间数据特点分析出发,运用矢量地理空间数据线状、面状要素具有的冗余数据的特性,首先在矢量地理数据中增加冗余数据点,然后通过控制冗余点在线状要素上的移动来嵌入水印信息,实现了高保真的矢量地理空间数据数字水印算法,最后对嵌入水印后数据的可用性以及算法的鲁棒性进行了试验分析。     

Seasonal global mean sea level change from satellite altimeter, GRACE, and geophysical models

We estimate seasonal global mean sea level changes using different data resources, including sea level anomalies from satellite radar altimetry, ocean temperature and salinity from the World Ocean Atlas 2001, time-variable gravity observations from the Gravity Recovery and Climate Experiment (GRACE) mission, and terrestrial water storage and atmospheric water vapor changes from the NASA global land data assimilation system and National Centers for Environmental Prediction reanalysis atmospheric model. The results from all estimates are consistent in amplitude and phase at the annual period, in some cases with remarkably good agreement. The results provide a good measure of average annual variation of water stored within atmospheric, land, and ocean reservoirs. We examine how varied treatments of degree-2 and degree-1 spherical harmonics from GRACE, laser ranging, and Earth rotation variations affect GRACE mean sea level change estimates. We also show that correcting the standard equilibrium ocean pole tide correction for mass conservation is needed when using satellite altimeter data in global mean sea level studies. These encouraging results indicate that is reasonable to consider estimating longer-term time series of water storage in these reservoirs, as a way of tracking climate change.     

Global land mapping of satellite-observed CO2 total columns using spatio-temporal geostatistics

ABSTRACT

This study presents an approach for generating a global land mapping dataset of the satellite measurements of CO₂ total column (XCO₂) using spatio-temporal geostatistics, which makes full use of the joint spatial and temporal dependencies between observations. The mapping approach considers the latitude-zonal seasonal cycles and spatio-temporal correlation structure of XCO₂, and obtains global land maps of XCO₂, with a spatial grid resolution of 1° latitude by 1° longitude and temporal resolution of 3 days. We evaluate the accuracy and uncertainty of the mapping dataset in the following three ways: (1) in cross-validation, the mapping approach results in a high correlation coefficient of 0.94 between the predictions and observations, (2) in comparison with ground truth provided by the Total Carbon Column Observing Network (TCCON), the predicted XCO₂ time series and those from TCCON sites are in good agreement, with an overall bias of 0.01?ppm and a standard deviation of the difference of 1.22?ppm and (3) in comparison with model simulations, the spatio-temporal variability of XCO₂ between the mapping dataset and simulations from the CT2013 and GEOS-Chem are generally consistent. The generated mapping XCO₂ data in this study provides a new global geospatial dataset in global understanding of greenhouse gases dynamics and global warming.     

数字地球建设中的地理空间Metadata研究

为解决数字地球建设中现有海量地理空间信息的共建共享,有效管理和利用,结合国际地理空间Metadata标准规范,对数字地球建设中所需的地理空间Metadata的定义,作用,分类,内容,实现模型与实例及其应用进行详细研究和设计,提出用于数字地球建设的地理空间Metadata网络管理系统解决方案,并设计实现了原型系统,为数字地球建设中现有海量地理空间信息的有效管理和共建共享提供了理论基础和切实可行的技术支撑。对国家空间信息基础设施建设和数字地球建设具有重要的推动作用和深远的社会意义。     

Rapid,high-resolution detection of environmental change over continental scales from satellite data – the Earth Observation Data Cube

The effort and cost required to convert satellite Earth Observation (EO) data into meaningful geophysical variables has prevented the systematic analysis of all available observations. To overcome these problems, we utilise an integrated High Performance Computing and Data environment to rapidly process, restructure and analyse the Australian Landsat data archive. In this approach, the EO data are assigned to a common grid framework that spans the full geospatial and temporal extent of the observations – the EO Data Cube. This approach is pixel-based and incorporates geometric and spectral calibration and quality assurance of each Earth surface reflectance measurement. We demonstrate the utility of the approach with rapid time-series mapping of surface water across the entire Australian continent using 27 years of continuous, 25?m resolution observations. Our preliminary analysis of the Landsat archive shows how the EO Data Cube can effectively liberate high-resolution EO data from their complex sensor-specific data structures and revolutionise our ability to measure environmental change.     

FY-3D MERSI-II地表温度遥感反演与验证

地表温度是是决定地表辐射能量收支的重要变量,在岩石圈、水圈、生物圈和大气圈的能量平衡和水量平衡研究中起着重要作用。利用热红外遥感技术可实现区域和全球尺度地表温度的快速获取,其受到了研究者的广泛关注。目前,FY-3D是国内光谱分辨率最高的对地观测卫星,极大的提高了对地观测能力,其搭载的中分辨率光谱成像仪（MERSI-II）经过大幅升级改进,性能有了显著提升,热红外数据的空间分辨率达到了250 m。本文使用大气辐射传输模型MODTRAN 5模拟了MERSI-II传感器热红外通道星上观测数据。在此基础上,构建了通用劈窗地表温度反演模型,结合ASTER GED全球地表发射率产品以及MERSI-II自身大气水汽反演算法,发展了地表温度遥感反演方法。最后,利用2019-08内蒙古乌海沙漠地区及美国SURFRAD多个站点的实测地表温度数据对本文提出的方法进行了验证。研究结果表明,相较地表实测数据,构建的劈窗算法反演的地表温度RMSE在1.6—2.6 K,反演精度达到了预期目标,还具有较高的空间分辨率,可以用于业务化的地表温度的反演,同时也说明其辐射定标精度有了一定保证,有效满足了区域和全球尺度地表温度遥感监测应用需求。     

Google Earth和World Wind比较研究

Google Earth和World Wind是目前最具代表性的两款基于网络的三维地理信息浏览器,为空间信息的共享发布提供了新的解决思路和技术手段。本文首先深入分析了这两款软件的技术特点,并在此基础上作了比较研究。     

多源地理空间矢量数据关联模型设计

多源地理空间矢量数据之间存在着隐含的关联关系,这些关联关系往往隐式存在,难以直观展示,也难以与空间数据映射交互展示,更无法进行查询分析,获取所需信息。针对这种情况,本文以多源地理空间矢量数据及统计数据为研究对象,首先定义了多源地理空间矢量数据关联的概念及分类,然后以此为基础设计了多源地理空间矢量数据关联模型,并将其分为3个子模型：基于自适应四叉树编码的空间关联子模型、基于几何匹配的空间关联子模型及基于语义匹配的空间关联子模型。该模型定义了多源地理空间矢量数据之间的关联方式,为关联关系的构建奠定了理论基础。     

基于O-QTM的球面VORONOI图的生成算法

提出了基于“Q－QTM”（Octahedral Quaternary Triangular Mesh）剖分的球面Voronoi图的格网生成算法；首先介绍了球面的QTM格网划发和编码方法，并根据地址码进行邻近球面三角形的探索；然后，参照数学形态学原理，重新定义了球面三角网的膨胀操作和膨胀算子，利用球面实体的递归膨胀来生成球面Voronoi图。应用VC^ 语言在OpenGL3维平台上开发了相应的实验程序，实验结果表明：利用此算法可生成球面上任意实体的Voronoi图，且生成点、弧和曲面Voronoi图的时间复杂度是一样的；而其误差受球面距离的影响较小，主要与球面实体的位置有关。最后给出了本文研究的结论及进一步的工作。