首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 140 毫秒
1.
 Two iterative vector methods for computing geodetic coordinates (φ, h) from rectangular coordinates (x, y, z) are presented. The methods are conceptually simple, work without modification at any latitude and are easy to program. Geodetic latitude and height can be calculated to acceptable precision in one iteration over the height range from −106 to +109 m. Received: 13 December 2000 / Accepted: 13 July 2001  相似文献   

2.
Gravity field terrain effect computations by FFT   总被引：2，自引：2，他引：2  
The widespread availability of detailed gridded topographic and bathymetric data for many areas of the earth has resulted in a need for efficient terrain effect computation techniques, especially for applications in gravity field modelling. Compared to conventional integration techniques, Fourier transform methods provide extremely efficient computations due to the speed of the Fast Fourier Transform (FFT. The Fourier techniques rely on linearization and series expansions of the basically unlinear terrain effect integrals, typically involving transformation of the heights/depths and their squares. TheFFT methods will especially be suited for terrain reduction of land gravity data and satellite altimetry geoid data. In the paper the basic formulas will be outlined, and special emphasis will be put on the practial implementation, where a special coarse/detailed grid pair formulation must be used in order to minimize the unavoidable edge effects ofFFT, and the special properties ofFFT are utilized to limit the actual number of data transformations needed. Actual results are presented for gravity and geoid terrain effects in test areas of the USA, Greenland and the North Atlantic. The results are evaluated against a conventional integration program: thus, e.g., in an area of East Greenland (with terrain corrections up to10 mgal), the accuracy ofFFT-computed terrain corrections in actual gravity stations showed anr.m.s. error of0.25 mgal, using height data from a detailed photogrammetric digital terrain model. Similarly, isostatic ocean geoid effects in the Faeroe Islands region were found to be computed withr.m.s. errors around0.03 m  相似文献   

3.
《测量评论》2013,45(9)
Abstract

The following method will be found better and quicker than the usual logarithmic process in computing the co-ordinates of intersected points in minor triangulation and traverse work. Let A and B be two stations whose co-ordinates (x 1 y 1), (x 2 y 2) are known. Let P be an intersected point whose co-ordinates (x, y) we wish to determine. Let α and β be the observed angles at A and B respectively.  相似文献   

4.
IntroductionThe map is a basic form of geographic informationvisualization[1]. To provide space attributes or geo-graphic orders is the basic function of a map. Incartography, according to the different measure ofphenomenal quantitative attribute, four fo…  相似文献   

5.
The Google Earth terrain model could prove beneficial for extraction of positional data in the future. At present, only an aging independent benchmark study (Potere, D., 2008. Horizontal position accuracy of Google Earth's high-resolution imagery archive. Sensors, 8, 7973–7981) provides constraints on positional accuracy for Google Earth imagery. In this investigation, we compared virtually traced positions against high-precision (<1 m) field measurements along three stratigraphic unconformity sub-sections in the Big Bend region to determine current positional accuracy for the Google Earth terrain model. A horizontal position accuracy of 2.64 m RMSEr was determined for the Google Earth terrain model with mean offset distance being 6.95 m. A vertical position accuracy of 1.63 m RMSEz with mean offset distance of 2.66 m was also calculated for the terrain model. Results suggest data extracted from the Google Earth terrain model could plausibly be used in future studies. However, we urge caution in using Google Earth data due to limited information disclosures by developers.  相似文献   

6.
Geocentre motion signals measured by satellite geodesy and those predicted from the observed mass redistribution in the ocean, atmosphere and terrestrial waters over 1993.1–2003.0 are analysed and compared under two viewpoints: the amplitudes and phases of the seasonal components, and the spectral signature of the non-seasonal components. The geodetic signals partly match the geophysical variations in the seasonal band, with possible remaining annual and semi-annual errors in both techniques, at the millimetre level in the equatorial plane for Satellite laser ranging (SLR) and Doppler Orbitography and radiopositioning integrated on Satellite (DORIS), and at the centimetre level in T z (Z-axis translation) for DORIS. Unlike SLR, the DORIS annual signatures in all three geocentre components have strongly varying amplitudes after 1996. The amplitude of the annual geophysical signal in T y is slowly growing with time. All three geophysical fluids contribute to this effect. The magnitude of the geophysically derived long-term geocentre motion is of the same magnitude in the T x , T y and T z directions, with a 0.5–1.0 mm Allan standard deviation for the 1-year sampling time, while the geodetic values are 2 mm in the equatorial plane for both SLR and DORIS, 4 mm for SLR and 9 mm for DORIS in the T z direction. The mismatch of the geodetic signal with the geophysical one in the inter-annual band is suggested to be due partly to excessive geodetic noise and partly to underestimated geophysical signal.  相似文献   

7.
Harmonic maps     
Harmonic maps are generated as a certain class of optimal map projections. For instance, if the distortion energy over a meridian strip of the International Reference Ellipsoid is minimized, we are led to the Laplace–Beltrami vector-valued partial differential equation. Harmonic functions x(L,B), y(L,B) given as functions of ellipsoidal surface parameters of Gauss ellipsoidal longitude L and Gauss ellipsoidal latitude B, as well as x(,q), y(,q) given as functions of relative isometric longitude =L–L0 and relative isometric latitude q=Q–Q0 gauged to a vector-valued boundary condition of special symmetry are constructed. The easting and northing {x(b,),y(b,)} of the new harmonic map is then given. Distortion energy analysis of the new harmonic map is presented, as well as case studies for (1) B[–40°,+40°], L[–31°,+49°], B0= ±30°, L0=9° and (2) B[46°,56°], L{[4.5°, 7.5°]; [7.5°, 10.5°]; [10.5°,13.5°]; [13.5°,16.5°]}, B0= 51°, L0 {6°,9°,12°,15°}.  相似文献   

8.
Based on the estimating rule of the normal vector angles between two adjacent terrain units, we use the concept of terrain complexity factor to quantify the terrain complexity of DEM, and then the formula of terrain complexity factor in Raster DEM and TIN DEM is deduced theoretically. In order to make clear how the terrain complexity factor E CF and the average elevation h affect the accuracy of DEM terrain representation RMSE Et , the formula of Gauss synthetical surface is applied to simulate several real terrain surfaces, each of which has different terrain complexity. Through the statistical analysis of linear regression in simulation data, the linear equation between accuracy of DEM terrain representation RMSE Et , terrain complexity factor E CF and the average elevation h is achieved. A new method is provided to estimate the accuracy of DEM terrain representation RMSE Et with a certain terrain complexity and it gives convincing theoretical evidence for DEM production and the corresponding error research in the future.  相似文献   

9.
Light detection and ranging (lidar) technologies provide a practical solution to 3-D terrain mapping through laser ranging and scanning technologies. A lidar footprint, which is the critical parameter describing the size of laser sampling area, varies with the scanning geometry and the local topography encountered. Integrating the effects of the scanning geometry and terrain orientations, this letter analyzes the intersection geometry of a laser beam and the terrain and develops a set of rigorous generalized footprint equations on inclined terrain for both across-track and Palmer scanning systems.   相似文献   

10.
Spatial selectivity estimation is crucial to choose the cheapest execution plan for a given query in a query optimizer. This article proposes an accurate spatial selectivity estimation method based on the cumulative density (CD) histograms, which can deal with any arbitrary spatial query window. In this method, the selectivity can be estimated in original logic of the CD histogram, after the four corner values of a query window have been accurately interpolated on the continuous surface of the elevation histogram. For the interpolation of any corner points, we first identify the cells that can affect the value of point (x, y) in the CD histogram. These cells can be categorized into two classes: ones within the range from (0, 0) to (x, y) and the other overlapping the range from (0, 0) to (x, y). The values of the former class can be used directly, whereas we revise the values of any cells falling in the latter class by the number of vertices in the corresponding cell and the area ratio covered by the range from (0, 0) to (x, y). This revision makes the estimation method more accurate. The CD histograms and estimation method have been implemented in INGRES. Experiment results show that the method can accurately estimate the selectivity of arbitrary query windows and can help the optimizer choose a cheaper query plan.  相似文献   

11.
Résumé Une des techniques de détermination fine et globale du champ de gravitation terrestre U est la gradiométrie spatiale, dans laquelle on mesure à bord d'un satellite sur orbite basse certaines combinations linéaires des composantes du tenseur ∂2 U/∂xi ∂xj dans des axes {x i } liés au satellite. Un tel projet, appelé GRADIO, est actuellement à l'étude en France et pourrait aboutir à partir de 1990. Après avoir rappelé les objectifs scientifiques d'une telle mission, nous en donnons les spécifications—étayées par une série d'études analytiques; nous définissons ensuite le satellite porteur et ses caractères techniques, en insistant sur les points délicats de la faisabilité (facteurs d'échelle des micro-accéléromètres constituant l'appareil, connaissance de l'attitude...) et en présentant des idées de solution en cours d'approfondissement.
Summary Satellite gradiometry arises as one of the methods for improving our knowledge of the global Earth gravity field at high resolution: by means of micro-accelerometers on board a low orbiting spacecraft, linear combination of the gravity tensor components ∂2 U/∂xi ∂xj are measured in a satellite-fixed reference frame {x i }. Based on this technique, a project named GRADIO is presently under study in France and could fly in 1990 at the earliest. After the scientific objectives of that experiment have been reviewed, the measurement specifications are given as coming from various analytical studies. The platform and its characteristics are then defined: the critical realization problems (scale factors of the micro-accelerometers, spacecraft attitude control and restitution) are pointed out together with some ideas for their solution which are under analysis and require further study.
  相似文献   

12.
Applying a one-step integrated process, i.e. by simultaneously processing all data and determining all satellite orbits involved, a Terrestrial Reference Frame (TRF) consisting of a geometric as well as a dynamic part has been determined at the observation level using the EPOS-OC software of Deutsches GeoForschungsZentrum. The satellite systems involved comprise the Global Positioning System (GPS) as well as the twin GRACE spacecrafts. Applying a novel approach, the inherent datum defect has been overcome empirically. In order not to rely on theoretical assumptions this is done by carrying out the TRF estimation based on simulated observations and using the associated satellite orbits as background truth. The datum defect is identified here as the total of all three translations as well as the rotation about the z-axis of the ground station network leading to a rank-deficient estimation problem. To rectify this singularity, datum constraints comprising no-net translation (NNT) conditions in x, y, and z as well as a no-net rotation (NNR) condition about the z-axis are imposed. Thus minimally constrained, the TRF solution covers a time span of roughly a year with daily resolution. For the geometric part the focus is put on Helmert transformations between the a priori and the estimated sets of ground station positions, and the dynamic part is represented by gravity field coefficients of degree one and two. The results of a reference solution reveal the TRF parameters to be estimated reliably with high precision. Moreover, carrying out a comparable two-step approach using the same data and models leads to parameters and observational residuals of worse quality. A validation w.r.t. external sources shows the dynamic origin to coincide at a level of 5 mm or better in x and y, and mostly better than 15 mm in z. Comparing the derived GPS orbits to IGS final orbits as well as analysing the SLR residuals for the GRACE satellites reveals an orbit quality on the few cm level. Additional TRF test solutions demonstrate that K-Band Range-Rate observations between both GRACE spacecrafts are crucial for accurately estimating the dynamic frame’s orientation, and reveal the importance of the NNT- and NNR-conditions imposed for estimating the components of the dynamic geocenter.  相似文献   

13.
激光雷达森林参数反演研究进展   总被引：6，自引：0，他引：6  
李增元  刘清旺  庞勇 《遥感学报》2016,20(5):1138-1150
激光雷达通过发射激光能量和接收返回信号的方式,来获取高精度的森林空间结构和林下地形信息。全波形激光雷达通过记录返回信号的全部能量,得到亚米级植被垂直剖面;离散回波激光雷达记录的单个或多个回波,表示来自不同冠层的回波信号。星载激光雷达一般采用全波形或光子计数激光剖面系统,仅能获取卫星轨道下方的单波束或多波束数据,用于区域/全球范围的森林垂直结构及变化观测。机载激光雷达多采用离散回波或全波形激光扫描系统,能够获取飞行轨迹下方特定视场范围内的扫描数据,用于林分/区域范围的森林结构观测。地基激光雷达多采用离散回波激光扫描系统,获取以测站为中心的球形空间内扫描数据,用于单木/样地范围的森林结构观测。激光雷达单木因子估测方法可分为CHM单木法、NPC单木法和体元单木法3类。CHM单木法通过局部最大值识别树冠顶点,采用区域生长或图像分割算法识别树冠边界或树冠主方向,NPC单木法一般通过空间聚类或形态学算法识别单木,体元单木法在3维体元空间采用区域生长或空间聚类算法识别树冠。根据激光雷达冠层高度分布可以估测林分因子,冠层高度分布特征来自于离散点云或全波形。多时相激光雷达可用于森林生长量、生物量变化等监测,以及森林采伐、灾害等引起的结构变化监测。随着激光雷达技术的发展,它将在森林调查、生态环境建模等生产与科学研究领域中得到更为广泛的应用。  相似文献   

14.
Identifying local spatial association in flow data   总被引：5，自引：1，他引：4  
In this paper we develop a spatial association statistic for flow data by generalizing the statistic of Getis-Ord, G i (and G i *). This local measure of spatial association, G ij, is associated with each origin-destination pair. We define spatial weight matrices with different metrics in flow space. These spatial weight matrices focus on different aspects of local spatial association. We also define measures which control for generation or attraction nonstationarity. The measures are implemented to examine the spatial association of residuals from two different models. Using the permutation approach, significance bounds are computed for each statistic. In contrast to the G i statistic, the normal approximation is often appropriate, but the statistics are still correlated. Small sample properties are also briefly discussed. Received: 18 February 1998/Accepted: 29 September 1998  相似文献   

15.
ABSTRACT

The effect of terrain shadow, including the self and cast shadows, is one of the main obstacles for accurate retrieval of vegetation parameters by remote sensing in rugged terrains. A shadow- eliminated vegetation index (SEVI) was developed, which was computed from only red and near-infrared top-of-atmosphere reflectance without other heterogeneous data and topographic correction. After introduction of the conceptual model and feature analysis of conventional wavebands, the SEVI was constructed by ratio vegetation index (RVI), shadow vegetation index (SVI) and adjustment factor (f (Δ)). Then three methods were used to validate the SEVI accuracy in elimination of terrain shadow effects, including relative error analysis, correlation analysis between the cosine of solar incidence angle (cosi) and vegetation indices, and comparison analysis between SEVI and conventional vegetation indices with topographic correction. The validation results based on 532 samples showed that the SEVI relative errors for self and cast shadows were 4.32% and 1.51% respectively. The coefficient of determination between cosi and SEVI was only 0.032 and the coefficient of variation (std/mean) for SEVI was 12.59%. The results indicate that the proposed SEVI effectively eliminated the effect of terrain shadows and achieved similar or better results than conventional vegetation indices with topographic correction.  相似文献   

16.
Following a brief review of the kernel regression approach to estimating surface models of the form z=f(x,y) +ε, this article will consider the situation where f is not a continuous surface function, and in particular where the discontinuities take the form of one‐dimensional breaks in the surface, and are not specified a priori. This form of model is particularly useful when visualizing some social and economic data where very rapid changes in geographical characteristics may occur – such as crime rates or house prices. The article briefly reviews approaches to this problem and proposes a novel approach (Bilateral Kernel Regression) adapting an algorithm from the field field of image processing (Bilateral Filtering), giving example analyses of synthetic and real‐world data. Techniques for enhancing the basic algorithm are also considered.  相似文献   

17.
We used RapidEye and Moderate Resolution Imaging Spectroradiometer (MODIS)/Terra data to study terrain illumination effects on 3 vegetation indices (VIs) and 11 phenological metrics over seasonal deciduous forests in southern Brazil. We applied TIMESAT for the analysis of the Enhanced Vegetation Index (EVI) and the Normalized Difference Vegetation Index (NDVI) derived from the MOD13Q1 product to calculate phenological metrics. We related the VIs with the cosine of the incidence angle i (Cos i) and inspected percentage changes in VIs before and after topographic C-correction. The results showed that the EVI was more sensitive to seasonal changes in canopy biophysical attributes than the NDVI and Red-Edge NDVI, as indicated by analysis of non-topographically corrected RapidEye images from the summer and winter. On the other hand, the EVI was more sensitive to terrain illumination, presenting higher correlation coefficients with Cos i that decreased with reduction in the canopy background L factor. After C-correction, the RapidEye Red-Edge NDVI, NDVI, and EVI decreased 2%, 1%, and 13% over sunlit surfaces and increased up to 5%, 14%, and 89% over shaded surfaces, respectively. The EVI-related phenological metrics were also much more affected by topographic effects than the NDVI-derived metrics. From the set of 11 metrics, the 2 that described the period of lower photosynthetic activity and seasonal VI amplitude presented the largest correlation coefficients with Cos i. The results showed that terrain illumination is a factor of spectral variability in the seasonal analysis of phenological metrics, especially for VIs that are not spectrally normalized.  相似文献   

18.
During a 4-year period starting in July 1996 and using intervals ranging from 3 days to 4 years, four precise polar motion (PM) series have been compared to excitation by atmospheric angular momentum (AAM) augmented with oceanic angular momentum (OAM) data. The first three series (C03, C04 and Bulletin A) are multi-technique combinations generated by the International Earth Rotation and Reference Systems Service (IERS) and the fourth combined series (IGS00P02) is produced by the International GPS Service (IGS) using only GPS data. The IGS PM compared the best with the combined excitations of atmosphere and oceans (AAM+OAM) at all intervals, showing high overall correlation of 0.8–0.9. Even for the interval of only three days, the IGS PM gave a significant correlation of about 0.6. Moreover, during the interval of February 1999 – July 2000, which should be representative of the current precision of the IGS PM, a significant correlation (>0.4) extended to periods as short as 2.2 days and 2.5 days for the xp and yp PM components, respectively. When using the IERS Bulletin B (C04) PM and an interval of almost 6 years, starting in November 1994, the combined OAM+AAM accounted for practically all the annual, semi-annual and Chandler wobble (CW) PM signals. When only AAM was used, either the US National Centers for Environment Prediction reanalysis data, which were used throughout this study, or the Japanese Meteorological Agency data, two large and well-resolved amplitude peaks of about 0.1 mas/day, remained at the retrograde annual and CW periods.  相似文献   

19.
This paper proposes a new approach to 3D modeling on the basis of laser scanner data. It presents the attempt to automatically detect two planes and an edge in one processing algorithm. The algorithm is based on the M split estimation, which is a recently developed regression method of multi functional models. Plane equations in three dimensional space are derived on the basis of the presented method. M split estimation divides the functional model of least squares adjustment in particular into two groups. Both simulated and real data from airborne laser scanning is tested to present the results. In each case the mehod gave good results which encourages for further work on the application of this method for 3D modeling.  相似文献   

20.
This paper presents a novel approach to building roof modeling, including roof plane segmentation and roof model reconstruction, from airborne laser scanning data. Segmentation is performed by minimizing an energy function formulated as multiphase level set. The energy function is minimized when each segment corresponds to one or several roof plans of the same normal vector. With this formulation, maximum n regions are segmented at a time by applying log2n level set functions. The roof ridges or step edges are then delineated by the union of the zero level contours of the level set functions. In the final step of segmentation, coplanar and parallel roof segments are separated into individual roof segments based on their connectivity and homogeneity. To reconstruct a 3D roof model, roof structure points are determined by intersecting adjacent roof segments or line segments of building boundary and then connected based on their topological relations inferred from the segmentation result. As a global solution to the segmentation problem, the proposed approach determines multiple roof segments at the same time, which leads to topological consistency among the segment boundaries. The paper describes the principle and solution of the multiphase level set approach and demonstrates its performance and properties with two airborne laser scanning data sets.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号