USACE Bare Earth LiDAR and derived projects for Fort Indiantown Gap near Lebanon, Pennsylvania


Identification_Information:
Citation:
Citation_Information:
Originator:Pennsylvania Department of Military Veterans Affairs
Publication_Date:2017
Title:
USACE Bare Earth LiDAR and derived projects for Fort Indiantown Gap near Lebanon, Pennsylvania
Online_Linkage: http://www.pasda.psu.edu
Description:
Abstract:
The National Guard Bureau (NGB) required high accruacy classified LiDAR data in combination with raster digital elevation models and hydrographic breaklines. For this effort, Continental Mapping Consultants (Continental) will collect and process high accuracy classified LiDAR data in .LAS format as well as a combination of raster digital elevation models.
Purpose:
The National Guard Bureau (NGB) requires the collection and processing of high accuracy classified LiDAR data in .LAS format as well as a combination of raster digital elevation models and additional hydrographic breaklines. The data is to be acquired during the Spring 2017 timeframe, during leaf-off conditions. The acquired LiDAR data will be used for various planning, design, research and mapping purposes. The NGB requires this data collection for Fort Indiantown Gap near Lebanon, Pennsylvania.
Time_Period_of_Content:
Time_Period_Information:
Single_Date/Time:
Calendar_Date:2017
Currentness_Reference:
ground condition
Status:
Progress:Complete
Maintenance_and_Update_Frequency:Unknown
Spatial_Domain:
Bounding_Coordinates:
West_Bounding_Coordinate:-76.755715
East_Bounding_Coordinate:-76512210
North_Bounding_Coordinate:40.493479
South_Bounding_Coordinate:40.384648
Keywords:
Theme:
Theme_Keyword_Thesaurus:None
Theme_Keyword:Elevation Data
Theme_Keyword:LiDAR
Theme_Keyword:DTM
Theme_Keyword:Land Surface
Theme_Keyword:Bare Earth LAS
Theme_Keyword:DEM
Theme_Keyword:LAS
Theme_Keyword:Ground Control
Theme_Keyword:Point Cloud
Theme_Keyword:GeoTIFF
Theme_Keyword:ERDAS Imagine
Theme_Keyword:Tile Index
Theme_Keyword:US Army Corps of Engineers
Theme_Keyword:USACE
Theme_Keyword:IMG
Place:
Place_Keyword_Thesaurus:None
Place_Keyword:Fort Indiantown Gap
Place_Keyword:Lebanon
Place_Keyword:Lebanon County
Place_Keyword:Pennsylvania
Place_Keyword:USA
Place_Keyword:PA
Access_Constraints:None
Use_Constraints:
None
Point_of_Contact:
Contact_Information:
Contact_Organization_Primary:
Contact_Organization:Continental Mapping Consultants
Contact_Address:
Address_Type:mailing and physical address
Address:
121 S Bristol Street
City:Sun Prairie
State_or_Province:WI
Postal_Code:53590
Country:USA
Contact_Voice_Telephone:888-815-3327
Data_Set_Credit:
Data credit goes to Continental Mapping Consultants who was contracted by US Army Corps of Engineers St. Loius.
Security_Information:
Security_Classification_System:N/A
Security_Classification:Unclassified
Security_Handling_Description:N/A
Native_Data_Set_Environment:
Environment as of Metadata Creation: Microsoft Windows 7 Professional; Esri ArcGIS 10.2.2; Microstation V8i Select Series 3 (version 8.11.09.459); TerraScan (version 016.001)
Back to Top
Data_Quality_Information:
Attribute_Accuracy:
Attribute_Accuracy_Report:
N/A
Logical_Consistency_Report:
Survey data have been confirmed to be in proper units, coordinate systems and format.  The terrain data have been confirmed as complete LAS format data files.
Completeness_Report:
Survey data have been checked for completeness, points have been collected in correct vegetation locations, and distributed throughout the AOI.  The terrain data have been checked for completeness against AOI polygons.
Positional_Accuracy:
Horizontal_Positional_Accuracy:
Horizontal_Positional_Accuracy_Report:
N/A
Vertical_Positional_Accuracy:
Vertical_Positional_Accuracy_Report:
USACE EM 1110-1-1000, 30, APR 2015 Requirements
Vertical Accuracy Class: 10 cm
RMSE Non-Vegetated: 10 cm
NVA at 95% Confidence Level: 19.6 cm
VVA at 95th Percentile: 30 cm

Non-Vegetated Vertical Accuracy
Vertical RMSE: 0.028

Vegetated Vertical Accuracy
Vertical RMSE: 0.041
Quantitative_Vertical_Positional_Accuracy_Assessment:
Vertical_Positional_Accuracy_Value:0.035 centimeters
Vertical_Positional_Accuracy_Explanation:
Non-Vegetated Vertical Accuracy (NVA) equal to the 95th Percentile confidence level calculated in open terrain.  Reported in centimeters.
Quantitative_Vertical_Positional_Accuracy_Assessment:
Vertical_Positional_Accuracy_Value:0.035 centimeters
Vertical_Positional_Accuracy_Explanation:
It passes the Vegetated Vertical Accuracy (VVA) requirements at 95th percentile.  Reported in centimeters.
Lineage:
Process_Step:
Process_Description:
Continental Mapping Consultants was tasked with performing a LiDAR ground control survey in support of LiDAR collection of Fort Indiantown Gap (FIG) in central Pennsylvania. The Global Positioning System (GPS) was used to establish the ground control network, and a total station was also used to survey the woods VVA checkpoints. The map in figure 1 shows the location of the ten ground control points (GCP, shown as red dots) and the twenty-five checkpoints (CP, shown as green dots). Also shown are two control stations which provided a tie to the NSRS (red cross inside white square).

The National Spatial Reference System (NSRS) was used to provide control for the network. An existing ground station of the National Spatial Reference System (NSRS) was occupied during the survey to provide a direct tie to the NSRS. In addition, a private CORS which is rigorously tied to the NSRS was also included.  The horizontal datum was NAD83 (2011) epoch 2010.0, and the vertical
datum was NAVD88, realized through the use of the GEOID12B model.

The field survey was done by using the KeyNetVirtual Reference Station Network. In the VRS method, a GPS receiver (a R10 dual frequency receivers was used) equipped with a data connection is used as a rover to occupy each survey point. Data is received via the cellular network for a “virtual” reference station nearby. The data is processed in real time to compute a position and elevation for the unknown station. The rover is positioned from the nearest physical CORS point and the vector data is stored in a data collector as a vector from this position, along with statistical data for the solution.
Process_Date:2017
Process_Step:
Process_Description:
The aerial LiDAR acquisition was completed by Keystone Aerial Surveys.  The sensor used was an Optech Galaxy.  There were 25 project flight lines, and 1 cross flights collected.  The LiDAR acquisition was completed in one (1) aircraft lift occurring on March 29, 2017.  

Altitude: 1800 meters
Aircraft Ground Speed: 113 knots average
Pulse Rate: 300 kHz
Scan Rate: 72 Hz
Full Field of View: 35 degrees
Full Swath Width: 680 meters
Swath Overlap: 30%
Average Point Density: 4.179 pts/m^2
Process_Date:2017
Process_Step:
Process_Description:
The Continental team utilized Optech's LMS to initially process the data and convert into LAS format.  TerraMatch was then used to refine the calibration of the LiDAR dataset.
Process_Date:2017
Process_Step:
Process_Description:
The trajectory files and point cloud swaths are imported into GeoCue to perform project setup.  This project set up phase sets the project parameters, tiling scheme, and is the platform for initial macro runs. After import, checkpoints are run against the point cloud to verify the accuracy of the data prior to classification. The detailed description of this process is below in 4.0 Accuracy Assessment. After verifying the accuracy, the processing continues. Multiple macros are run through TerraScan to flag overlap, and to classify the ground. Due to differing terrain, this step may take multiple iterations. Once the analyst has verified the results with the ground macro, the ground classification QC begins. During the QC phase, analysts are reclassifying the point cloud in areas where the macro was not able to, or were misclassified.
Multiple macros are run on the dataset after the ground classification is complete including the building and water macros. The same rigorous quality steps are performed each classification.  The bare earth lidar points that were within 1 meter of the water were classified to class 10. After the analysts have completed the QC process in TerraScan, raster files were produced into 32-bit floating GeoTiffs and Erdas Imagine IMG files using Global Mapper. These files were created using only the ground and building classes.
Process_Date:2017
Process_Step:
Process_Description:
The calibrated and controlled lidar files were processed using automatic point classification routines in proprietary software. These routines operate against the entire collection (all swaths, all lifts), eliminating character differences between files.  The trajectory files and point cloud swaths were imported into GeoCue to perform project setup.  This project set up phase set the project parameters, tiling scheme, and was the platform for initial macro runs.  After import, checkpoints were run against the point cloud to verify the accuracy of the data prior to classification.  After verifying the NVA accuracy, the processing continued.  Multiple macros were run through TerraScan to flag overlap, ground, low points, etc.  Points below ground surface that were identified as low points were classified to class 7. A bare earth ground surface was derived from the unclassed points and put to class 2 using a suitable macro function for the project's terrain type. Points above the treeline that were not identified as a feature were classified to class 18. Ground points inside of water features were classifed to class 9 to represent water and ground points outside of hydro features but within 1 meter of hydro breaklines were classified to class 10. All remaining points were classified to class 1.  Final lidar LAS delivery classes for the fully classed LAS tiles consistent with ASPRS LAS classes to be compliant with USGS LiDAR Guidelines and Base Specifications v13 consist of: Class 1 - Unclassified; Class 2 - Ground; Class 7 – Low Noise; Class 8 – Model Key; Class 9 - Water; Class 10 - Ignored Ground (including 2-meter buffer around water breaklines); Class 12 - Overlap points exceeding scan angles of +/- 20 degrees; and Class 18 - High Noise.  Once the analyst was comfortable with the ground macro results, the ground classification QC began using TerraScan.  During the QC phase, analysts reclassified the point cloud in areas where the macro was not able to, or created misclassifications.  The same rigorous quality steps were performed on each classification.

Data were then distributed as virtual tiles to experienced lidar analysts for localized automatic classification, manual editing, and peer-based QC checks. Supervisory QC monitoring of work in progress and completed editing ensured consistency of classification character and adherence to project requirements across the entire project.

Upon completion of point classification, an automated process was executed to turn any points to class 1 in the Fully Classified LAS files and to delete any point for the Bare Earth and Key Point LAS Files that fell outside of the provided project buffered boundary.

Breaklines were digitized at water elevation for any bodies of water over the entire project area including streams greater than 25ft in nominal width, water bodies greater than 0.5 acres in area, and islands greater than 0.5 acre. A macro was then run to classify points that lay at nominal water elevation to class 9 that fell within bodies of water. Concurrently a 1-meter buffer zone around water polygons was derived from the ground class and put to class 10. Water bodies and streams with flow were hydro-enforced using the Ponds, Islands, and Double Line Drain shapefiles to demonstrate the removal of unnatural surface artifacts in both ponds and streams and to show downward flow for streams.

A DEM data set was generated for both the Bare Earth LAS set and Model Key Point LAS set as 32-bit Erdas Imagine .IMG files at a resolution of 1.0 m (meter) supported by the Hydro Breaklines.

1-foot contours were generated in LP360 using the hydro breaklines and Bare Earth data.  
Process_Date:2017
Process_Step:
Process_Description:
Continental utilized various software packages and techniques to verify the accuracy of the data. Utilizing QCoherent’s LP360, Continental ran a survey to las check, followed by seamline analysis. The survey to las check calculates the deviation between the survey point elevation and the point cloud elevation and export a Non-Vegetated Vertical Accuracy (NVA) report. The second check, calculates the deviation between the seamlines of the point cloud swaths. The third check, the Vegetated Vertical Accuracy (VVA) testing occurred after the ground classification has been completed. The VVA testing utilized both the points that were labeled VVA and NVA if those points were either located in or under vegetation. Other software like Terra Solid and Global Mapper were utilized to verify the results of LP360.

Continental appreciates the time-critical nature of this product. The field survey and aerial survey teams were deployed at the first opportunity based on availability of acceptable weather conditions and base personnel (for coordination). The area of interest contains subareas of dense vegetation which present fewer bare ground returns, higher variability, and potentially less accuracy than typical vegetated areas.  Per the table below the survey accuracy results meet industry standards for both NVA and VVA based on the ground survey control points collected in December.

Once all of the deliverables have been produced and verified, the data was moved to the Quality office for final review. The Quality Office verifies that the correct procedures were followed, tests the data, and verifies that all of the deliverables in the SOW are finished.
Process_Date:2017
Back to Top
Spatial_Data_Organization_Information:
Direct_Spatial_Reference_Method:Vector
Point_and_Vector_Object_Information:
Back to Top
Spatial_Reference_Information:
Horizontal_Coordinate_System_Definition:
Planar:
Planar_Coordinate_Information:
Planar_Coordinate_Encoding_Method:coordinate pair
Coordinate_Representation:
Abscissa_Resolution:0.500000
Ordinate_Resolution:0.500000
Planar_Distance_Units:meters
Geodetic_Model:
Horizontal_Datum_Name:D_NAD_1983_2011
Ellipsoid_Name:WGS_1984
Semi-major_Axis:6378137.00
Denominator_of_Flattening_Ratio:298.257224
Vertical_Coordinate_System_Definition:
Altitude_System_Definition:
Altitude_Datum_Name:North American Vertical Datum of 1988
Altitude_Resolution:0.01
Altitude_Distance_Units:meters
Altitude_Encoding_Method:Explicit elevation coordinate included with horizontal coordinates
Back to Top
Entity_and_Attribute_Information:
Detailed_Description:
Entity_Type:
Entity_Type_Label:NGB_FIG_2017_Lidar_Boundary_UTM18N_WGS84
Attribute:
Attribute_Label:FID
Attribute_Definition:
Internal feature number.
Attribute_Definition_Source:
Esri
Attribute_Domain_Values:
Unrepresentable_Domain:
Sequential unique whole numbers that are automatically generated.
Attribute:
Attribute_Label:Shape
Attribute_Definition:
Feature geometry.
Attribute_Definition_Source:
Esri
Attribute_Domain_Values:
Unrepresentable_Domain:
Coordinates defining the features.
Attribute:
Attribute_Label:Id
Attribute:
Attribute_Label:Proj_Desc
Attribute:
Attribute_Label:ProjCoord
Back to Top
Distribution_Information:
Distributor:
Contact_Information:
Contact_Organization_Primary:
Contact_Organization:Pennsylvania Spatial Data Access (PASDA)
Contact_Address:
Address_Type:mailing address
Address:
115 Land and Water building
City:University Park
State_or_Province:Pennsylvania
Postal_Code:16802
Country:United States
Contact_Voice_Telephone:(814) 865 - 8792
Contact_Electronic_Mail_Address:pasda@psu.edu
Resource_Description:N/A
Distribution_Liability:
The USER shall indemnify, save harmless, and, if requested, defend those parties involved with the development and distribution of this data, their officers, agents, and employees from and against any suits, claims, or actions for injury, death, or property damage arising out of the use of or any defect in the FILES or any accompanying documentation. Those parties involved with the development and distribution excluded any and all implied warranties, including warranties or merchantability and fitness for a particular purpose and makes no warranty or representation, either express or implied, with respect to the FILES or accompanying documentation, including its quality, performance, merchantability, or fitness for a particular purpose. The FILES and documentation are provided "as is" and the USER assumes the entire risk as to its quality and performance. Those parties involved with the development and distribution of this data will not be liable for any direct, indirect, special, incidental, or consequential damages arising out of the use or inability to use the FILES or any accompanying documentation.
Standard_Order_Process:
Non-digital_Form:
N/A
Custom_Order_Process:
N/A
Technical_Prerequisites:
N/A
Back to Top
Metadata_Reference_Information:
Metadata_Date:2017
Metadata_Contact:
Contact_Information:
Contact_Organization_Primary:
Contact_Organization:Continental Mapping Consultants
Contact_Address:
Address_Type:mailing and physical address
Address:
121 S Bristol Street
City:Sun Prairie
State_or_Province:WI
Postal_Code:53590
Country:USA
Contact_Voice_Telephone:888-815-3327
Metadata_Standard_Name:FGDC Content Standards for Digital Geospatial Metadata
Metadata_Standard_Version:FGDC-STD-001-1998
Metadata_Time_Convention:local time
Metadata_Access_Constraints:None
Metadata_Use_Constraints:
None
Metadata_Security_Information:
Metadata_Security_Classification_System:None
Metadata_Security_Classification:Unclassified
Metadata_Security_Handling_Description:
None
Metadata_Extensions:
Online_Linkage: None
Profile_Name:None
Back to Top