Digital photogrammetry project with very high-resolution stereo pairs acquired by DigitalGlobe, Inc. satellite Worldview-2

Transcription

1 White PAPER Greater area of the City of La Paz, Bolivia Digital photogrammetry project with very high-resolution stereo pairs acquired by DigitalGlobe, Inc. satellite Worldview-2 By: Engineers Nelson Mattie, Alexis Piña, Ramón López, and Leonel Perozo Summary A digital photogrammetry study was conducted on an area of 750 Has in the city of La Paz, Bolivia. The project s aim was to perform a survey with level curves equidistant at 1m and to render in 3D elements such as roads, hydrography, vegetation, and buildings. Also a digital terrain model (DTM) on a one (1) square meter resolution grid was produced. All of the foregoing is generated by means of photogrammetric processing of 50 cm resolution stereo pairs from the WorldView-2 satellite. Accuracy results of coordinates (X, Y) for all the area was less than 38 cm RMSE and elevation accuracy (Z) for all the area was 29 cm RMSE and determined by 5 main geodesic survey points for triangulation, 50 geodesic control points measured in the field by means of GPS and electronic level, and 35 checkpoints. These geodesic control points were used for aerial triangulation and altimetric leveling of stereo pairs. The difference between values (X, Y, Z) from the 50 points measured in the field and the points measured in the stereo block or model was 2.8 cm RMSE. By means of these state-of-the-art digital sensors, and by employing advanced mathematical models (software), the precision obtained by our satellite s photogrammetric products is similar to that obtained through conventional topographic processes. The applications for our products range from highway planning to management of agricultural production, feasibility studies for hydroelectric stations, urban planning, environmental analysis, land registries, oil surveying, or road infrastructure projects among others. Location These 750 Has are located within two work polygons, North and South in the greater area of the City of La Paz, Bolivia (Figures 1 and 2). Figure 1. North Polygon. The round yellow points indicate the location of geodesic control points 1

2 Figure 2. South polygon. Yellow points indicate the location of geodesic control points Geodesy: Based on the SIRGAS points laid over the territory of Bolivia in 1995, the Military Geographical Institute (IGM for its acronym in Spanish) has created a national geodesic reference framework network constituted by a continuous GPS, a semicontinuous GPS network, and the passive GPS network that currently operates with 125 points. The Military Geographical Institute of Bolivia, alongside the Andes Project central have installed eight (8) GPS stations for data collection, which have been used as reference for the Bolivian State s MARGEN framework network. The stations were used to process our data For differential GPS measurements all necessary information was collected, which included: General information regarding the INGA point (point belonging to the MARGEN framework net work of the State of Bolivia), which a permanent station for continuous GPS data collection by the Military Geographical Institute of Bolivia. Design of technical specifications and geodesic parameters. For the execution of static GPS measurements for the Control Points, the following specifications were used: Capture interval: 5 seconds. Cut-off: 15 º Measurement time: Depends on the length of base lines (higher after 1.5 hours of stay). 2

3 Figure 3. Data from the main GPS measurements of geodesic control points. For the height of the different vertices that constitute the geodesic network, the national geodesic system of the Plurinational State of Bolivia was taken as reference by means of its main BM or points of known height. The reference for height is the average sea level as reference, this means that such height is orthometric and is located at a source point named Arica on which all heights from this system are based. A limitation in the State of Bolivia is the lack of a GEOIDAL reference framework, which would allow the calculation of Geoidal heights, through a model approved and certified by Bolivia s IGM, without the need to obtain them by geometric leveling. Figure 4. Geometric leveling of geodesic control points.. 3

4 Capture of stereo pairs Only a week after the satellite was programmed, Digital Globe Inc. managed to successfully complete a stereoscopic pair with satellite Worldview-2 on September 1, (For further reference please see www. digitalglobe.com) The first stage of the photogrammetric project based on stereoscopic pairs taken by the WV-2 satellite consisted on performing aero triangulation supported by the ground geodesic control points measured in the area of study with L1/L2 frequency GPS equipment to guarantee ground controls as part of accuracy adjustments and assurance. In this case 5 main geodesic vertices were used, as well as 50 control points and 35 checkpoints. Triangulation and adjustment were performed on the photogrammetric blocks; homologous points were calculated between images to produce epipolar imagery and then proceed to automatically and manually extract topographic data by means of algorithms developed for this purpose. Afterwards, digital surface models (DSM) were obtained as well as the digital terrain models (DTM), level curves were manually edited. The precision results were: Average Residual Error X, Y (rmsex, rmsey)= m; Average Residual Error Z (RmseZ)=0.290m. Figure 5. Stereo pair and location of the geodesic control points. The overlapping percentage between the images of the stereo pair was 99,30%. 4

5 Figure 6. Measurement of checkpoints and homologous points on the stereo pair. Figure 7. Topographic drawings with level curves at 1m 5

6 Figure 8. Measurement of checkpoints and homologous points on the stereo pair. Figure 9. X,Y, Z point cloud generated from WV-2 stereo pairs used to generate the digital terrain model. Approximately 52.5 million points were generated in LAS format. Figure 11. X,Y,Z dot cloud generated from the WV-2 stereo pairs used to generate the digital terrain model. Approximately 52.5 million points were generated in LAS format. 6

7 COMPARATIVE CHART OF Z VALUES (ORTHOMETRIC HEIGHT) BETWEEN THE MAIN 25 GROUND CONTROL POINTS MEASURED WITH TOPOGRAPHIC EQUIPMENT IN THE FIELD AND THE SAME POINTS MEASURED IN THE STEREO MODEL POINT NORTH (m) EAST (m) HEIGHT (H) Ground Height (H) Stereo Residual TL TL TL TL TL TL TL TL TL TL TL TL TL TL TL TL TL TL TL TL TL TL TL TL TL Avg: pt/

8 Digital terrain model for the North polygon Digital terrain model for the South polygon 8

9 Figure 8. 3D rendering of roads, vegetation, construction, and hydrography Figure 9. 3D rendering of roads, vegetation, construction, and hydrography 9