nders and ArcGIS combined with three-dimensional photorealistic modeling for mapping outcrops in the Slick Hills, Oklahoma

The mapping of geology is conventionally done visually in a hands-on fashion, and the data are recorded in a fi eld book or with photography. An alternative technique that combines refl ectorless laser rangefi nders or high-speed terrestrial laser scanners, global positioning system, and the Environmental Systems Research Institute (ESRI) ArcGIS software platform has been developed that is effective for mapping geology at a distance and in three dimensions. Portable handheld refl ectorless lasers are used to capture geologic features such as contacts and terrain and can be combined with digital elevation models in ArcGIS software. Fast terrestrial laser scanners capture an entire exposure at the detail and accuracy of (3D) photorealistic (virtual) models with the additional color information from image pixels. This latter method is expensive and complicated and requires signifi cant amounts of fi eld and processing effort. The laser gun approach is simple, portable, and cost effective. When integrated with ESRI ArcGIS software and a module, such as our recently developed ArcGIS extension 3DLT (laser tool), a simple yet sophisticated platform exists for mapping, visualizing, and analyzing outcrops in real time in the fi eld. The potential of laser mapping is demonstrated in the Paleozoic outcrops of a structural geology teaching site in the Slick Hills, Oklahoma. Fast laser scanning and digital photography are used to build a 3D photorealistic model of an area of the anticline. The 3DLT is used for mapping specifi c detailed features such as contacts and faults. Three-dimensional quantitative information can be extracted from the geology with these methods. A laser rangefi nder combined with 3DLT can image and display terrain and outcrop features in the fi eld, in real time. Mapping with fast scanners requires several steps in processing of the point cloud data utilizing a variety of sophisticated and expensive software, but can capture an entire outcrop, such as a mountainside. The resulting model then can be analyzed in the lab. When combined with digital photography, virtual photorealistic models derived from point clouds can be even more effectively analyzed. The most appropriate method for digitally mapping geology depends on a variety of issues, such as cost, time, complexity, portability, and the project goals.


INTRODUCTION
Common methods for digital geologic mapping have included aerial photography, global positioning system (GPS), and the utilization of conventional laser and/or optical based surveying tools. Digital geologic mapping of outcrops has also been carried out with various combinations of GPS, cameras, high-speed laser scanners, and survey instruments (Xu et al., 2001;McCaffrey et al., 2005;Thurmond et al., 2005;Bellian et al., 2005;Pringle et al., 2006;Oldow et al., 2006). Methods and technology developed in geographic information system (GIS) mapping applications, such as remote sensing, are used in geologic mapping, but for detailed mapping of exposures on steep slopes and cliffs it is necessary to map digitally and obliquely at relatively close range (often <1 km) using portable and perhaps handheld equipment. Refl ectorless laser rangefi nders (Gilbert, 1996) such as handheld guns and total stations (Lyman et al., 1997;Nielsen et al., 1999Nielsen et al., , 2000Xu et al., , 2001Zeng et al., 2004;Xu, 2000;Ramirez-Ugalde, 2002) have been an effective technology for geological mapping, operating at 1-250 points/s out to ranges of as much as 1 km and vertical and horizontal angle accuracy as precise as 0.001° (total station) and 0.01° (handheld rangefi nders). Terrestrial laser scanners (TLS) programmed to run at high speed (thousands to hundreds of thousands of points per second) generate dense point clouds (millions of points) at angle and range accuracies better than the handheld rangefi nders. These scanners have been used to map geology, creating three-dimensional (3D) models that capture entire mountainsides (Thurmond et al., 2005;Bellian et al., 2005).
Improvements have been made to these methods for digital mapping of geology by using enhancements to the laser rangefi nders (e.g., Bluetooth communications), centimeter accuracy real-time kinematic global positioning system (RTK GPS), and an ArcGIS-based mapping software system. The improvements create more effective 3D analysis and allow for real-time visualization of the results. TLS scans, when combined with digital photography, can be used to build 3D photorealistic (virtual) models to capture large areas of geology at great detail, and this workfl ow has been changed and made more effi cient with new software and techniques. How these laser mapping approaches can be or should be applied is explored, and these concepts are applied to map geology in the Slick Hills, Oklahoma.

Refl ectorless Laser Rangefi nding
Refl ectorless laser rangefi nders are designed to measure distance to low-refl ectance surfaces such as buildings, soil, rock, and vegetation. They do not require highly refl ective targets such as corner cube prisms (Lichti et al., 2002). The maximum range that can be measured is dependent upon the physical properties of the target surface (electric permittivity and magnetic permeability, which determine the absorption characteristics of the surface), geometric characteristics of the surface, such as roughness, and the angle of incidence of the laser beam to the surface (Fowles, 1975). A dark surface, such as coal or basalt, will absorb more of the laser energy than a lighter surface such as limestone. A low angle of incidence will increase the refl ection of the beam away from the laser and will spread the beam over a greater area, diminishing the strength of the returned signal. A smooth surface, such as fi ne, compacted dirt, will cause the laser beam to refl ect away from the laser rather than backscatter toward it. The negative impact on range caused by the surface smoothness increases as the angle of incidence between the laser beam and the surface decreases, and at very low angles of incidence there will be little returned signal. When specifi cations of the range of a laser rangefi nder are noted by a manufacturer, they list two ranges depending on the degree of refl ectivity of the material; 20% refl ectivity (typical for most rocks) and 80% refl ectivity (white surfaces typical for constructed objects such as buildings). For a Laser Atlanta Advantage CI rangefi nder, this would be 800 m and 1000 m, respectively. The laser beam diverges at a rate of 3 milliradians, resulting in a laser spot almost 1 m in diameter at 1000 m. These are the dominant factors affecting the maximum target range in geological applications (Amann et al., 2001).
The angular accuracy of the handheld rangefi nders with an internal fl uxgate compass is typically 0.5° in azimuth. By our own tests the angle accuracy can be improved by more than 10× by mounting the rangefi nder on an angular encoder. The use of the encoder rather than an internal compass requires that a backsight be used, as is done with a total station, in order to georeference the measurements. The vertical inclinometers in the rangefi nders typically have accuracies to 0.2°, but are limited to maximum zenith angles between 45° and 60°, and can also be improved in accuracy when using a vertical encoder by a factor of 10. If the internal compass is used for horizontal angles rather than an angular encoder, caution must be exercised to ensure that the measurements are not affected by spurious magnetic fi elds such as power lines. The data from each laser shot are stored on media compatible with the particular laser, such as internal memory and fl ash memory devices (PC card, compact fl ash card, secure digital card), or the rangefi nder can be connected directly to a portable computer.
Various software programs have been developed to support connection of the rangefi nder directly to a portable computer. The CyberMap program of the University of Texas at Dallas (UTD) (built with Environmental Systems Research Institute [ESRI] Arcview and MapObjects)  provides an interface to both a rangefi nder and a total station and maps the points as they are shot. Other programs that provide for an interface to rangefi nders are Condor Earth Technologies PenMap (Hillman, 1995), the University of California at Berkeley GeoMapper, which uses ESRI PenMap (Brimhall and Vanegas, 2001), Carlson's Tsunami (Turner, 2002), and the Geological Survey of Canada's (GSC) GanFeld, using the ESRI ArcPad (Guy and Buller, 2004). The GSC was a pioneer in developing data logging geology software, where one stands on a point and records into a program what is measured or observed (Brodaric, 2000). Far more common are nongeological application software programs, such as the Trimble Navigation GPS receiver software program TerraSync, which can be linked to a rangefi nder. However, the lack of a continuous sampling mode and the inability to view the data as points, lines, and polygons in a 3D viewer on the computer have made them less effective, especially for geological applications.
In order to capture the features of a geological outcrop or terrain, a laser rangefi nder is used to sketch the features to form a 3D image of the feature. A total station or rangefi nder in continuous sampling mode can be swiveled to sketch the terrain or the features of the geological outcrop. Single-shot measurements can be made, but the preferable mode of operation, due to operator fatigue, is a continuous sampling, or "trigger on" mode. The resolution with which the feature is sketched is determined by the rate at which the laser is swept across the feature. A faster sweep rate results in a lower resolution sketch. By sweeping the mounted rangefi nder along salient features in the terrain, the features are sketched as a 3D digital image. Xu (2000) used a handheld rangefi nder mounted on a custom-designed cradle with an angular encoder to map terrain features such as bedding layers. The rangefi nder was operated at a slow 1 shot/s rate up to the maximum rate of 250 shots/s. Features as much as 1 km away were sketched in this fashion.
Laser rangefi nders with continuous sampling have proven to be effective in mapping terrain and outcrop features, especially when mapping linear features such as bedding layers, faults, and contacts between different lithologies. The details of the characteristics of the different rangefi nder guns vary, but the dominant issues in determining whether it is acceptable for a given application are equipment cost (much less than that of a high-speed terrestrial scanner), range accuracy (10 cm), sample rate (1-250 Hz), horizontal accuracy (0.5°), vertical accuracy (0.2°) with built-in compass and inclinometer, and whether it can be adapted to angular encoders. In addition, the size, weight, and power requirements of the rangefi nder will affect the applications for which it is suitable, or preferable to high-speed terrestrial scanners.

3DLT Extension
Because GIS software platforms provide powerful tools to manage, edit, and visualize spatial data, a software extension to ESRI ArcGIS has been developed to support 3D fi eld analysis of geological features and terrain using laser rangefi nders. The program combines many of the capabilities of the UTD Cybermapping program and the GCS GanFeld program, integrating the rangefi nder data with GPS, and implementing it in the ArcGIS environment. The 3DLT extension operates the laser rangefi nder from the ArcGIS environment ( Fig. 1) and in real time displays the data in 2D and 3D views as points. The ability to display the work in progress in the 3D ArcScene viewer allows the geologist to monitor his or her work, as it is being captured, on top of a previously loaded digital elevation model (DEM) or other 3D model of the terrain or outcrop. There is the ability to assign feature names to the point attribute table while capturing the data, as well as the ability to convert the points to lines and polygons as the outcrop or terrain are being scanned. These abilities provide real-time feedback to the geologist as he or she sketches the features.
The 3DLT extension is written in the Visual Basic for Application (VBA) language and combined with ESRI ArcObjects. ArcObjects are the C++ building blocks of the ArcGIS software platform. They provide services to support GIS applications on the computer desktop. With ArcObjects, menus, tools, workfl ows, applications, and custom feature classes can be created for use by the VBA (Visual Basic for Applications) code, which implements the 3DLT extension. In real time, the 3DLT extension can create 3D point, polyline, and polygon shapefi les from the rangefi nder point data. The rangefi nder can operate in either single-shot mode or continuous-shot mode, and the data are displayed on the screen as they are captured by the rangefi nder. The 3DLT extension will accept angular data for each point from both the rangefi nder's internal compass and inclinometer or from external angular encoders. If the scanner was georeferenced with a backsight prior to starting work, the 3D features that are produced will be georeferenced at the time of capture. The point data (X, Y, and Z) can be exported as ASCII (American Standard Code for Information Interchange) fi les, and ASCII point data can be imported into the program.
The 3DLT extension was tested against the clean lines of a constructed structure, a corridor at UTD (Fig. 2). The points in Figure 2 are shown in 3D mode as they were captured. The corners, walls, and doors of the corridor were sketched with the rangefi nder and then converted to polygons and polylines ( Fig. 3): the walls are polygons and the door frames are polylines. Sketching the corridor in continuous sampling mode resulted in the straight lines of the wall intersections and doorways not being perfectly straight. This was caused by a combination of the operator's inability to hold the rangefi nder absolutely steady when tracing a straight line and the accuracy of the laser ranging and the accuracy of the angles from the vertical inclinometer (a horizontal encoder was used).

TLS
An alternative to the handheld laser rangefi nder is the high-speed TLS, or ground light detection and ranging (LIDAR). Refl ectorless ground LIDAR systems can scan outcrops and terrain from thousands of points per second to hundreds of thousand of points per second. Systems such as Cyrax (McCaffrey et al., 2005) use rotating mirrors to defl ect the laser beam in order to scan across the image fi eld. Systems such as the Riegl LPM swivel the laser up and down and rotate the laser horizontally in order to scan the image fi eld.
The scanners are larger and heavier than the handheld rangefi nders and considerably more expensive. They require more power than the handheld scanner, which involves carrying larger batteries when the scan site is not accessible by vehicle. Real-time conversion of the scanned image into a 3D georeferenced framework is not yet available. The scanned point clouds are post-processed and georeferenced by scanning georeferenced control points during the work phase and then using a six variable affi ne transformation to convert the point coordinates from the scanner coordinate system to the georeferenced coordinate system. When viewed from the perspective of the scanner, the point clouds appear almost like a photograph, especially if the points have been texture mapped with an integrated camera image, as can be done with the Riegl Z420i. The intensity of the refl ected signal for each point is saved and can be viewed in grayscale, which often highlights features in the outcrop or terrain. The point cloud can later be converted into a triangular irregular network (TIN) model, and photographs that were taken of the outcrop or terrain can be draped onto the model , creating a photorealistic model of the outcrop or terrain. These techniques have been used to map geological structures with point clouds (Bellian et al., 2005), and by others to create surface models using TIN meshes and photographic texture mapping (Thurmond et al., 2005). Few, so far, have used these techniques for mapping geological structures due to the complexity and cost of the equipment and associated specialized software.

Digitally Mapping Geological Features
In order to assess the effectiveness of the developed tools and software, geological features in the Slick Hills of Oklahoma were mapped using both the 3DLT handheld laser rangefi nder system and a high-speed terrestrial LIDAR. UTD developed software was used to convert the LIDAR point cloud into a 3D photo-realistic model. The Slick Hills area is ~24 km north of Lawton, Oklahoma, along State Highway 58 (Fig. 4). A UTD structural geology class (lead by Weldon Beachamp) combined with the UTD Cybermapping Laboratory to map part of the Kimbell anticline located in the Kimbell Ranch, which is in the southern part of Blue Creek Canyon (Fig. 5). The area includes rocks from the Kindblade and Cool Creek Formations of the Arbuckle Group that are exposed along the Kimbell anticline. Both the Kindblade and Cool Creek Formations are marine carbonate buildups deposited during the Early Ordovician Period (Rigby and Toomey, 1978). Throughout the Pennsylvanian Period, tectonics associated with the formation of the Wichita Mountains resulted in left-lateral shearing throughout the canyon that led to the development of several major faults, creating a horst and graben terrain (Donovan et al., 1986). Over time, the intense tectonic forces acting on this region began to subside, and erosion during the Permian reduced it to its present Slick Hills topography (Donovan et al., 1988). It is dominated by largescale deformation; there are several major anticlines and synclines present. The steeply dipping resistant ledges characteristic of the Slick Hills makes the anticlines and synclines obvious. Vegetation in the area is low to the ground; therefore the geology is very well exposed. The area that was digitally mapped includes rocks from the Kindblade and Cool Creek Formations

Figure 3. The conversion of three-dimensional (3D) point data to lines and polygons. The corners of the corridor were mapped (upper left) and then used to defi ne polygons by converting the points into lines (upper right) and then polygons (lower right) and then volumes (lower right).
of the Arbuckle Group exposed along the Kimbell anticline. A hillside exposure of the Kimbell anticline consists of intermittent bed-scale outcrops in tall dry grasses and small juniper trees. This was the area of detailed mapping.
A soft-copy raster image of a hand-drawn map was georeferenced, then converted into a shapefi le to be able to geographically compare in GIS the newly captured features. The features were then integrated with the vectorized georeferenced version of the original hard-copy map and an available 10 m DEM. All shapefi les were imported into GPS Track Maker (GPSTM), free software, and then used in conjunction with shapefi les, ArcGIS, and Google Earth software to form an integrated model of the area.

Digitally Mapping the Kimbell Anticline
Geologic contacts such as stratigraphy were captured as points using the 3DLT system, RTK GPS for georeference control, and a modifi ed Laser Atlanta Optics rangefi nder with a Bogen/ horizontal encoder confi guration tripod. This is not the most stable confi guration, but it is relatively lightweight and portable. The points were then converted to lines using 3DLT to defi ne stratigraphic contacts (Fig. 6) and subsequently were overlain on the previous geologic mapping and a digital terrain model (Fig. 7). The points were converted into polygons in order to defi ne a plane that allows computation of their strike and dip orientations. At this stage the various strike and dip routines available can be used, such as MicroImages TOOL SCRIPT (MicroImages, 2005), and those of Vacher (1989), and Fienen (2005), to compute the strikes and dips. The data were imported into GPS TrackMaker to import and/or export shapefi les from and to ArcGIS and Google Earth, exported as shapefi les into ArcGIS, and then imported into Google Earth and overlain on terrain models of the Slick Hills (Fig. 8). Integrating the data with Google Earth maps provides a global context and a visualization evaluation of the accuracy of the captured fi eld data points as well as lines and planes processed from the data in the lab.

TLS of the Kimbell Anticline
A TLS scan was carried out to a range of 150-200 m on the north side of the anticline with a Riegl LPM i800 scanner (Fig. 9) from the University of Idaho. At this distance a scan increment of 4 mgon (400 gons = 360°) suffi ciently captures at ~1000 points/s detail outcrops at an interval of a few centimeters. Note in the fi gures the gaps in the data, due to the angle of the scanners with relation to the terrain across drainage. At least two scan locations would have been required to completely map the area to fi ll in where shadowing creates gaps in the point data and possibly in the surfaces, but for our purpose the scans were suffi cient. Vegetation usually does not allow coherent mesh building and therefore results in holes. Four scans, over 2 h, from one site captured the area of interest (Fig. 10). While the scanner scanned the outcrop, digital photographs were taken of the hillside. After data acquisition, the raw scan data were loaded into Riegl LPMScan software and exported as an ASCII fi le with X, Y, Z and intensity data. The data were then loaded into point cloud software, InnovMetric Polyworks software. The IMAlign module was used to create a merge of the point cloud (Fig. 11). The merged point was loaded into the IMMerge module and a polygonal surface was created. The polygonal mesh (triangulated irregular network, TIN model) was edited and made ready for photo draping, which can be very time consuming if outward appearance is important. If the geologic mapping appearance is not as important, then less editing would be necessary, except in the area to be analyzed.
It has been said that the point cloud is sufficient for any 3D geologic analysis (Bellian et al., 2005) with 3D data extraction accomplished by point-cloud editing software such as Polyworks. Even a colored point cloud loses its information when it is observed at close range (Aiken et al., 2004a(Aiken et al., , 2004b, whereas a photographic image is only limited by the pixel density on a   Focus Area: The points represent the data captured using 3DLT system and the lines are hand drawn by Beauchamp Captured point data using 3DLT integrated with hand-drawn structural geology by Beauchamp displayed in ArcMap triangulated mesh or other such surface. For this reason, if surface analysis is to be performed on a model generated from a laser scan point cloud, it is best to have the photographic texture draped onto the model. The photographs highlight discontinuities in the lithology of the surface.

Photorealistic Modeling
Methods used for generating 3D photorealistic models have been developed at UTD using a series of nonlinear transformations to relate features in the real world of XYZ to the pixels on a photograph in UV space (used to defi ne x, y positions on an image). This results in a model where the photograph is accurately "draped" onto the model and is not the result of rubber sheeting methods (Xu, 2000;Xu et al., 2000, Olariu et al., 2006Thurmond et al., 2005). A new work fl ow (Fig. 12) utilizing the UTD 3DPM software was used, and resulted in working with the raw point clouds in Riegl RiScan Pro, then exporting them as point clouds in an ASCII format, and subsequently importing them into Innovmetric Polyworks software, in which the different scans from various positions are merged as if referenced to a single ref-erence scan position. In Polyworks, editing of the point cloud and then the triangulated (TIN) mesh can be carried out. The editing is done to the extent necessary to yield the desired results. This is usually the most time-consuming step, taking days to months depending on the size and problems with the data sets (e.g., a lot of vegetation will produce tangled triangulated meshes that may have to be repaired, such as fi lling the holes and cleaning up the model; the work is then exported as a model in a Wavefront fi le format [.obj]). All of these steps are performed under Polyworks. Next, at least four common points are picked between each photo and the 3D model using Polyworks and our own 3D photorealistic mapping (3DPM) software. UTD designed software solves for the transformation parameters relating the surface mesh to the photos. Finally, the photos are draped (really transformed) onto the mesh models using the UTD 3DPM software, resulting in a 3D photorealistic (virtual) model (Fig. 13). At this stage in the processing, Polyworks or other such similar 3D analysis or CADD (computer assisted design and drafting) software can be used to trace (digitize) geologic contacts on the 3D photorealistic model (Fig. 14). This is equivalent to digitizing the rocks directly with a laser gun in 3DLT mapping. In either case a contact mapped as a series of points can be used for solving three-point problems, such as strikes and dips on surfaces fi t to the tracing. 3DLT has the capability to carry out this kind of 3D analysis on the virtual model as well as with a laser on a real outcrop (Fig. 15) (e.g., MicroImages, 2005).

CONCLUSIONS
The 3D point clouds, triangulated meshes, or the photo models can and were used for extraction of 3D surfaces (Figs. 16A, 16B) so that strikes and dips can be computed. Specifi c stratigraphic contacts were mapped by 3DLT on both sides of a ridge and overlain onto Google Earth models. Care must be taken when extracting and/or digitizing these geometric shapes. The Slick Hills topography has more resistant layers that stand out in relief, or lines of narrow and parallel ledges. These types of exposures are captured with laser scanners, but raise a question. With either a rangefi nder in the fi eld or with a cursor on a 3D photorealistic mesh model in the lab, which linear feature should be digitized in order to defi ne the contact? For example, is   the digitizing done along the edge of the ledge, or the contact with the next layer? This is something a geologist must decide. The edges of the resistant layers were consistently digitized in this case. Preferably, it should be done by the geologist familiar with the geology and the features to be extracted, because such extraction is a step in the analysis and interpretation of geology.
The 3DLT approach is very fast and cost effective at an accuracy of centimeters to deci-meters, but it only maps features of interest at that time the work is performed. If more information needs to be mapped, then the site must be revisited. The scanner-based 3D photorealistic approach captures everything in an area within line of sight at higher accuracy (centimeters), but it requires much more expensive hardware and software, the procedure is much more complicated, and it takes longer to acquire and build the model (although this is improv-   Fig. 9, color coded, but seen from above using Polyworks to display the gaps in the data due to shadowing (blockage of the scanning). ing). In the Slick Hills example, both of these approaches proved effective in this steeply dipping, grassy environment. The 3DLT provided results in real time at centimeter to decimeter accuracy for immediate evaluation and discussion. 3D photorealistic modeling, after considerable post-processing in the lab, was successful at an accuracy level of a few centimeters and was compared with conventional fi eld and 3DLT mapping. Digitizing with a cursor on the 3D model produced almost the same 3D quantitative results as direct fi eld "digitizing" by laser mapping (tracing) when these points are fi t to surfaces for analysis of their orientation. These results were generated at many locations of the geologic exposures, exceeding that which could be done manually with a compass, and these digital results provide statistics of those measurements (e.g., standard deviations) that can be archived allowing future reanalysis. Importing the resultant lines and surfaces in the model also allows effective visualization of the topography when mapped onto Google Earth scenes. These systems complement each other and both belong to the geologists' toolbox for mapping and analysis.