Direct creation of construction drawings for Kwan-yin Statue by close-range photogrammetry

In order to provide the model point accuracy of ±2 mm, an extra high accuracy industrial control net with four surveying piers with accuracy ±0.1 mm was set up around the model. Tens of orientation marks have been placed on the different parts and measured with accuracy ±0.2 mm from the above-mentioned piers. The main four stereopairs taken with the P31 camera around the model are simultaneously processed on the BC2 analytical plotter. The accuracy of absolute orientation is better than ±0.9 mm. Finally, about 300 sections construction drawings have been directly offered, and the closerange photogrammetric data has been used for the design and construction of the statue.     

Improving Orthometric Heights Using Precise GPS Measurements

The vertical component obtained from the Global Positioning System (GPS) observations is from the ellipsoid (a mathematical surface), and therefore needs to be converted to the orthometric height, which is from the geoid (represented by the mean sea level). The common practice is to use existing bench marks (around the four corners of a project area and interpolate for the rest of the area), but in many areas bench marks may not be available, in which case an existing geoid undulation is used. Present available global geoid undulation values are not generally as detailed as needed, and in many areas they are not known better than ±1 to ±5 m, because of many limitations. This article explains the difficulties encountered in obtaining precise geoid undulation with some example computations, and proposes a technique of applying corrections to the best available global geoid undulations using detailed free-air gravity anomalies (within a 2° × 2° area) to get relative centimeter accuracy. Several test computations have been performed to decide the optimal block sizes and the effective spherical distances to compute the regional and the local effects of gravity anomalies on geoid undulations by using the Stokes integral. In one test computation a 2° × 2° area was subdivided into smaller surface elements. A difference of 37.34 ± 1.6 cm in geoid undulation was obtained over the same 2° × 2° area when 1° × 1° block sizes were replaced by a combination of 5' × 5' and 1' × 1' subdivision integration elements (block sizes).     

顾及GRACE季节影响的华北平原水储量变化反演

基于移去-恢复方法,利用GRACE卫星重力数据研究了华北平原2003年1月—2014年6月的陆地水变化,提出了一种顾及季节影响尺度因子计算方法,用于减小GRACE后处理误差。将本文方法所得结果与水文模式和降水模型结合进行比较分析,结果表明：2008年之前华北平原发生比较明显的干旱现象,陆地水和地表水分别以（-7.9±2.4）mm/a和（-7.3±2.8）mm/a的速度下降,但地下水仅以（-0.6±1.4）mm/a的速度减少;2008年之后,陆地水和地表水分别以（4.3±1.3）mm/a和（10.9±2.1）mm/a的速度上升,期间地下水的超采严重,平均以（-6.5±1.2）mm/a的速度下降,但下降速度以0.9 mm/a²的趋势减缓。研究期间内,华北平原的陆地水、地表水和地下水整体分别呈现（-2.0±0.6）、（2.9±0.7）和（-4.8±0.7）mm/a的变化趋势。最后,利用交叉小波谱分析了GRACE滞后于降水的现象。研究结论表明,降水和地下水开采是影响华北平原水储量变化的两大决定性因素。     

航测多功能综合检校场自动对中地标的设置

为满足摄影测量免像控点的直接法测图(DR)的生产市场需求,对相机进行独立检验的现状已无法满足DR测图相机的检校需求,必须建设针对DR测图设备的集成系统进行全面检校的综合场地。本文结合863项目的实施,对项目实施中所建立的地对地模拟实验场建设中的地标问题进行了探讨,根据模拟飞行高度50m的大比例尺地形测量需要,确定地标点的测量精度为±1㎜的精度,设计了适合摄影用的地标形状、大小、颜色,在此基础上设计出满足±1㎜的强制对中设备,可满足不同摄影相机的需要,特别是对作业用的对中设备的检验方法进行了全面系统检验方案设计。     

Improvement of the geoid in local areas by satellite-to-satellite tracking

Summary The concept of satellite-to-satellite tracking measuring the relative velocity of two orbiting satellites spaced some hundreds kilometers on a close orbit, provides now possibilities for the investigation of the Earth’s gravity field. In the paper only medium and short wave length effects affecting the measured relative velocity have been considered. Collocation is used in such an analysis of local geoid improvement, because this method allows to combine heterogeneous data in a consistent way. Covariance functions relevant for the particular case of a circular equatorial orbit are given. Two kinds of observation equations have been formulated. The choice of observation equation with regard to satellites configuration is discussed. It is found that it is sufficient to have a limited number of satellite-to-satellite observations in a 7^o×7^o area around the estimation point with distances between profiles of about 1^o.5 and between the two satellites forming the pair of 200+350 km; the altitude of satellite-to-satellite observations should be as low as possible. The accuracy of the geoid determination strongly depends on the degree and order of the reference field used. An accuracy of about ±1 m can be achieved with an assumed reference field of (40,40). The influence of measuring errors is discussed and it is shown that only satellite-to-satellite observations with accuracy better then 0.1 mm/sec will give an improvement of the geoid. Finally, some results on the combination of low-low satellite-to-satellite tracking and terrestrial gravity data are given. The proposed method seems to be especially interesting for unsurveyed areas. Furthermore, it has the practical advantage that only a local coverage data is needed.     

The ERS-2 altimetric bias and gravity field enhancement using dual crossovers between ERS and TOPEX/Poseidon

 Dual satellite crossovers (DXO) between the two European Remote Sensing satellites ERS-1 and ERS-2 and TOPEX/Poseidon are used to (1) refine the Earth's gravity field and (2) extend the study of the ERS-2 altimetric range stability to cover the first four years of its operation. The enhanced gravity field model, AGM-98, is validated by several methodologies and will be shown to provide, in particular, low geographically correlated orbital error for ERS-2. For the ERS-2 altimetric range study, TOPEX/Poseidon is first calibrated through comparison against in situ tide gauge data. A time series of the ERS-2 altimeter bias has been recovered along with other geophysical correction terms using tables for bias jumps in the range measurements at the single point target response (SPTR) events. On utilising the original version of the SPTR tables the overall bias drift is seen to be 2.6±1.0 mm/yr with an RMS of fit of 12.2 mm but with discontinuities at the centimetre level at the SPTR events. On utilising the recently released revised tables, SPTR2000, the drift is better defined at 2.4±0.6 mm/yr with the RMS of fit reduced to 3.7 mm. Investigations identify the sea-state bias as a source of error with corrections affecting the overall drift by close to 1.2 mm/yr. Received: 25 May 2000 / Accepted: 24 January 2001     

STEREOMICROSCOPE PHOTOGRAMMETRY FOR THE MEASUREMENT OF SMALL OBJECTS

A stereoscopic microscope has been used in combination with the body of a 35 mm camera to obtain stereophotographs, at scales up to 6:1, of objects only a few millimetres in size. Subsequent stereocomparator observations have been used to deduce height differences on the objects 10 a precision, for the mean of five sets of readings, of ±0.01 mm for closely spaced points. Using more widely separated points, the precision was estimated to be ± 0.05 mm. Depths were computed using × parallaxes after a number of parameters were fixed by calibration. The large convergence angle of the microscope axes was a crucial parameter; rotations about the camera axis were found to be significant variables.     

Plate Motion of India and Interseismic Strain in the Nepal Himalaya from GPS and DORIS Measurements

We analyse geodetically estimated deformation across the Nepal Himalaya in order to determine the geodetic rate of shortening between Southern Tibet and India, previously proposed to range from 12 to 21 mm yr^?1. The dataset includes spirit-levelling data along a road going from the Indian to the Tibetan border across Central Nepal, data from the DORIS station on Everest, which has been analysed since 1993, GPS campaign measurements from surveys carried on between 1995 and 2001, as well as data from continuous GPS stations along a transect at the logitude of Kathmandu operated continuously since 1997. The GPS data were processed in International Terrestrial Reference Frame 2000 (ITRF2000), together with the data from 20 International GNSS Service (IGS) stations and then combined using quasi- observation combination analysis (QOCA). Finally, spatially complementary velocities at stations in Southern Tibet, initially determined in ITRF97, were expressed in ITRF2000. After analysing previous studies by different authors, we determined the pole of rotation of the Indian tectonic plate to be located in ITRF2000 at 51.409±1.560° N and ?10.915±5.556°E, with an angular velocity of 0.483±0.015°. Myr^?1. Internal deformation of India is found to be small, corresponding to less than about 2 mm yr^?1 of baseline change between Southern India and the Himalayan piedmont. Based on an elastic dislocation model of interseismic strain and taking into account the uncertainty on India plate motion, the mean convergence rate across Central and Eastern Nepal is estimated to 19±2.5 mm yr^?1, (at the 67% confidence level). The main himalayan thrust (MHT) fault was found to be locked from the surface to a depth of about 20 km over a width of about 115 km. In these regions, the model parameters are well constrained, thanks to the long and continuous time-series from the permanent GPS as well as DORIS data. Further west, a convergence rate of 13.4±5 mm yr^?1, as well as a fault zone, locked over 150 km, are proposed. The slight discrepancy between the geologically estimated deformation rate of 21±1.5 mm yr^?1 and the 19±2.5 mm yr^?1 geodetic rate in Central and Eastern Nepal, as well as the lower geodetic rate in Western Nepal compared to Eastern Nepal, places bounds on possible temporal variations of the pattern and rate of strain in the period between large earthquakes in this region.     

GPS relative positioning at a precision level of one part per billion

Six days of data from the GIG'91 experiment have been analyzed with a fiducial strategy. The results obtained with our orbital software GEOSAT, show daily horizontal and length repeatabilities at the level of 1 part in 10⁹ plus 2 mm for baseline lengths up to 4000 km. A direct comparison with results from the GIPSY software shows, with some exceptions, mean differences at the sub-cm level. After transformation to ITRF'90 the rms of the coordinate differences is 14.8 mm. Studies of orbital predictions and comparisons with external high precision orbits indicate a mean orbit precision and accuracy of around 35 cm in each cartesian coordinate. Correlations between the GEOSAT and GIPSY solutions indicate some common model deficiencies.     

Determination of geoidal heights and deflections of the vertical for the hellenic area using heterogeneous data

Mean gravity anomalies, deflections of the vertical, and a geopotential model complete to degree and order180 are combined in order to determine geoidal heights in the area bounded by [34°≦ϕ≤42°, 18°≦λ≦28°]. Moreover, employing point gravity anomalies simultaneously with the above data, an attempt is made to predict deflections of the vertical in the same area. The method used in the computations is least squares collocation. Using empirical covariance functions for the data, the suitable errors for the different sources of observations, and the optimum cap radius around each point of evaluation, an accuracy better than±0.60m for geoidal heights and±1″.5 for deflections of the vertical is obtained taking into account existing systematic effects. This accuracy refers to the comparison between observed and predicted values.     

VLBI data analysis with the GEOSAT software

The GEOSAT software for satellite geodetic and geodynamics applications has lately been extended to Very Long Baseline Interferometry (VLBI) applications. Geodetic VLBI data have been analyzed using this software and station coordinates and velocities, source coordinates, and earth orientation parameters have been derived in a global mode using data from the Extended Research and Development Experiment (ERDE) and the Research and Development series (R&D) within the NASA's CDP.It is demonstrated using the ERDE data, that the program is capable of calculating station coordinates with an accuracy of a few mm in the local horizontal plane and about 7 mm in the vertical direction. Analysis of the R&D dataset yields errors 2 to 3 times higher in all directions.     

ITRF2008: an improved solution of the international terrestrial reference frame

ITRF2008 is a refined version of the International Terrestrial Reference Frame based on reprocessed solutions of the four space geodetic techniques: VLBI, SLR, GPS and DORIS, spanning 29, 26, 12.5 and 16?years of observations, respectively. The input data used in its elaboration are time series (weekly from satellite techniques and 24-h session-wise from VLBI) of station positions and daily Earth Orientation Parameters (EOPs). The ITRF2008 origin is defined in such a way that it has zero translations and translation rates with respect to the mean Earth center of mass, averaged by the SLR time series. Its scale is defined by nullifying the scale factor and its rate with respect to the mean of VLBI and SLR long-term solutions as obtained by stacking their respective time series. The scale agreement between these two technique solutions is estimated to be 1.05 ± 0.13 ppb at epoch 2005.0 and 0.049 ± 0.010?ppb/yr. The ITRF2008 orientation (at epoch 2005.0) and its rate are aligned to the ITRF2005 using 179 stations of high geodetic quality. An estimate of the origin components from ITRF2008 to ITRF2005 (both origins are defined by SLR) indicates differences at epoch 2005.0, namely: ?0.5, ?0.9 and ?4.7?mm along X, Y and Z-axis, respectively. The translation rate differences between the two frames are zero for Y and Z, while we observe an X-translation rate of 0.3?mm/yr. The estimated formal errors of these parameters are 0.2?mm and 0.2?mm/yr, respectively. The high level of origin agreement between ITRF2008 and ITRF2005 is an indication of an imprecise ITRF2000 origin that exhibits a Z-translation drift of 1.8?mm/yr with respect to ITRF2005. An evaluation of the ITRF2008 origin accuracy based on the level of its agreement with ITRF2005 is believed to be at the level of 1?cm over the time-span of the SLR observations. Considering the level of scale consistency between VLBI and SLR, the ITRF2008 scale accuracy is evaluated to be at the level of 1.2?ppb (8?mm at the equator) over the common time-span of the observations of both techniques. Although the performance of the ITRF2008 is demonstrated to be higher than ITRF2005, future ITRF improvement resides in improving the consistency between local ties in co-location sites and space geodesy estimates.     

CBERS-02B卫星遥感影像的区域网平差

针对中巴资源一号卫星(CBERS-02B)卫星遥感影像姿态角误差较大的特点,提出了利用区域网平差方法提高其对地目标定位精度的策略和具体计算过程。首先对参与平差的每景影像选取4个地面控制点进行影像姿态角常差检校,然后采用与地形无关方案解求各自的RPC参数,最后选取带仿射变换项的有理函数模型(RFM)进行多重覆盖影像的区域网平差。对两个地区的0级CBERS-02B单条CCD立体影像对的区域网平差试验表明,对地目标定位的平面和高程精度均达到了±3个像元的水平,且高程精度明显优于平面精度。相对于常规的卫星遥感影像区域网平差方法,平面和高程精度均有明显提升,几乎达到国外同等高分辨率卫星遥感影像的几何定位精度。这说明中国卫星遥感影像亦具有较好的几何定位潜力,在区域网平差之前进行系统误差预改正是必要和可行的。     

高速铁路CPⅢ平面网复测若干技术问题探讨

高速铁路CPⅢ平面网需要在不同施工阶段进行复测,复测过程中的主要技术困难有大跨度连续梁段埋设的CPⅢ点随梁体伸缩变化存在多值性,无法满足随时使用的坐标准确性小于等于±3mm。采用线性回归模型对大跨度连续梁段的CPⅢ点平面坐标进行预测,基本解决了多值性的问题,同时提出在复测阶段选用合格的cPⅢ点代替CPII对其平面网进行约束,使复测后的成果尽量与原测成果相匹配。该方法可增强CPⅢ点使用的准确性,为现场施工提高效率。     

Accuracy assessment of the GPS-TEC calibration constants by means of a simulation technique

During the last 2 decades, Global Positioning System (GPS) measurements have become a very important data-source for ionospheric studies. However, it is not a direct and easy task to obtain accurate ionospheric information from these measurements because it is necessary to perform a careful estimation of the calibration constants affecting the GPS observations, the so-called differential code biases (DCBs). In this paper, the most common approximations used in several GPS calibration methods, e.g. the La Plata Ionospheric Model (LPIM), are applied to a set of specially computed synthetic slant Total Electron Content datasets to assess the accuracy of the DCB estimation in a global scale scenario. These synthetic datasets were generated using a modified version of the NeQuick model, and have two important features: they show a realistic temporal and spatial behavior and all a-priori DCBs are set to zero by construction. Then, after the application of the calibration method the deviations from zero of the estimated DCBs are direct indicators of the accuracy of the method. To evaluate the effect of the solar activity radiation level the analysis was performed for years 2001 (high solar activity) and 2006 (low solar activity). To take into account seasonal changes of the ionosphere behavior, the analysis was repeated for three consecutive days close to each equinox and solstice of every year. Then, a data package comprising 24 days from approximately 200 IGS permanent stations was processed. In order to avoid unwanted geomagnetic storms effects, the selected days correspond to periods of quiet geomagnetic conditions. The most important results of this work are: i) the estimated DCBs can be affected by errors around ±8 TECu for high solar activity and ±3 TECu for low solar activity; and ii) DCB errors present a systematic behavior depending on the modip coordinate, that is more evident for the positive modip region.     

新一代全球海底地形模型BAT_WHU2020

本文利用由多源卫星测高资料计算的新版全球重力异常Grav_Alti_WHU,联合船测水深资料,构建了全球75°S—70°N范围的1′×1′海底地形模型BAT_WHU2020。以船测水深、现有模型和多波束测深数据为参考,对模型精度进行了分析评价。结果表明,在中国海域及邻区(104°E—160°E,0°N—50°N),本文模型与船测水深之差值的标准差约70 m,与SIO V19.1模型精度相当,优于ETOPO1、DTU10、GEBCO_08等模型,较此前发布的BAT_VGG模型精度提高了约30%,说明本文模型构建方法可靠、数据处理准确、精度较高。在全球范围内,BAT_WHU2020模型与船测水深之差值的标准差为50～65 m,差值在±200 m范围内的比率超过95%,与SIO V19.1模型精度相当,优于ETOPO1、DTU10、GEBCO_08等模型,较BAT_VGG模型精度提高了27%～36%。以SIO V19.1模型为参考,模型之差的标准差为90～110 m,约90%格网点差值在±200 m以内,约95%格网点差值在±300 m以内,两者一致性良好。最后,讨论了地壳均衡、Parker公式高次项等对成果精度的影响,模型的真实空间分辨率,以及以多波束测深为参考的模型精度问题。分析认为,BAT_WHU2020模型空间分辨率为10～18 km,在马里亚纳海沟、麦夸里海岭地区相对精度为5%～6%。     

联合卫星重力、卫星测高和海洋资料研究全球海平面变化

利用GRACE、卫星测高和海洋实测温盐数据,探讨了2003~2012年间全球海平面、比容海平面和海水质量等的变化特征,并讨论了南极冰盖和格陵兰冰盖消融对全球海平面变化的影响。全球海平面整体呈上升趋势,上升速度为2.72±0.07 mm/a,且存在显著的空间分布特征。全球海平面、比容海平面和海水质量等的变化还具有显著的季节性特征,其中全球海平面变化的年周期振幅为4.6±0.3 mm。使用经验正交函数分析(EOF)得到全球海平面和比容海平面的季节性变化在南北半球存在显著的差异,但海水质量季节性变化不存在这种差异。南极冰盖和格陵兰冰盖的消融速率分别为-75.7±12.3 Gt/a和-124.1±2.9 Gt/a,对海平面的长期趋势项贡献分别为0.21±0.03 mm/a和0.34±0.01 mm/a,仅占全球海水质量增加速度1.80±0.10 mm/a的12%和19%,总计占31%,因此,两极冰盖质量消融并不是2003-2012年间海水质量增加的最主要因素。     

Stereo-Photogrammetric Mapping of Tooth Replicas Incorporating Texture

Commercial digital photogrammetric software has been applied to convergent stereoscopic photography of human tooth replicas prepared to exhibit optical texture resulting in successful generation of 3D coordinate data. Tooth replicas were imaged using a semi-metric 35 mm camera and f = 100 mm macro lens on extension bellows. Model precision was within acceptable limits of 12  μ m or better for manual target matching and 21  μ m or better for automatic image matching. Further improvement in optical texture is required to achieve automatic image matching precision comparable to that of manual target matching. Small errors in interior orientation parameters attributed to instability in the bellows as well as small errors in the relative orientation resulted in some systematic errors. The use of a fixed camera lens system is expected to reduce these errors. When combined with commercially available, moderately priced, digital SLR cameras this brings 3D model generation closer to everyday clinical dental practice.     

Land cover and land use classification performance of machine learning algorithms in a boreal landscape using Sentinel-2 data

ABSTRACT

In recent years, the data science and remote sensing communities have started to align due to user-friendly programming tools, access to high-end consumer computing power, and the availability of free satellite data. In particular, publicly available data from the European Space Agency’s Sentinel missions have been used in various remote sensing applications. However, there is a lack of studies that utilize these data to assess the performance of machine learning algorithms in complex boreal landscapes. In this article, I compare the classification performance of four non-parametric algorithms: support vector machines (SVM), random forests (RF), extreme gradient boosting (Xgboost), and deep learning (DL). The study area chosen is a complex mixed-use landscape in south-central Sweden with eight land-cover and land-use (LCLU) classes. The satellite imagery used for the classification were multi-temporal scenes from Sentinel-2 covering spring, summer, autumn and winter conditions. Using stratified random sampling, each LCLU class was allocated 1477 samples, which were divided into training (70%) and evaluation (30%) subsets. Accuracy was assessed through metrics derived from an error matrix, but primarily overall accuracy was used in allocating algorithm hierarchy. A two-proportion Z-test was used to compare the proportions of correctly classified pixels of the algorithms and a McNemar’s chi-square test was used to compare class-wise predictions. The results show that the highest overall accuracy was produced by support vector machines (0.758 ± 0.017), closely followed by extreme gradient boosting (0.751 ± 0.017), random forests (0.739 ± 0.018), and finally deep learning (0.733 ± 0.0023). The Z-test comparison of classifiers showed that a third of algorithm pairings were statistically different. On a class-wise basis, McNemar’s test results showed that 62% of class-wise predictions were significant from one another at the 5% level or less. Variable importance metrics show that nearly half of the top twenty Sentinel-2 bands belonged to the red edge (25%) and shortwave infrared (23%) portions of the electromagnetic spectrum, and were dominated by scenes from spring (38%) and summer (40%). The results are discussed within the scope of recent studies involving machine learning and Sentinel-2 data and key knowledge gaps identified. The article concludes with recommendations for future research.