Description: State Game Lands within the State of Pennsylvania, United States.
This dataset consists of LiDAR derivatives covering the State Game Lands (SGLs) of Pennsylvania. It was created with the objective of automatically detecting all relict charcoal hearths (RCHs) in these SGLs. For details on these results, please see “When Computers Dream of Charcoal: Using Deep Learning, Open Tools and Open Data to Identify Relict Charcoal Hearths in and around State Game Lands in Pennsylvania” .
Herein we detail the steps taken to assemble the raw data and process it into a usable format. The unprocessed LiDAR data used for this project was obtained through Pennsylvania Spatial Data Access (PASDA) . The process included identifying the data needed to encompass the SGLs and a 1 kilometer buffer around them, downloading the raw LiDAR files in .las format, organizing these files, and processing them into the following LiDAR derivatives: digital elevation map (DEM), hillshade, and slope analysis. The text below and (Figure 1) describe the detailed steps.
Pennsylvania’s LiDAR data is organized in two grids for the north and the south of the state, which overlap slightly [3, 4]. These grids were merged into a single shapefile via the “Merge Vector Layers” tool in QGIS. Shapefiles for these grids contain attributes such as direct links to download data for each grid tile. Next, we downloaded a shapefile containing SGL boundaries  provided by Pennsylvania and added a one-kilometer buffer around each SGL, which was given the same identifying number as the SGL it surrounded. By overlaying this buffer with the unified grid and using the “Join Attributes by Location” function, any grid tile that touched an SGL or its buffer was given that SGL’s number as an attribute. This provided a single shapefile of the grid with two central attributes for each tile, one with the SGL/buffer number and one with a download link for the LiDAR data for that tile.
Exporting the attribute table for the state’s tile grid into Microsoft Excel allowed us to manipulate the data more easily. We deleted all grids that were not associated with a SGL/buffer. This resulted in a selection of 3,925 unique tiles to download out of the 13,797 (28.4%) total tiles in the grid. One field given for each grid tile was a direct link to download the .las file, the raw LiDAR point cloud. Copying this list of URLs into a text document allowed us to use simple bash scripting and the wget command to automate the downloading of these files. Note that the list of tiles downloaded was actually over 4,000, not exactly 3,925. This is because tiles were downloaded consecutively by each individual SGL and some tiles overlapped more than one SGL.
In order to use bash scripting on a Windows computer, we enabled the Linux Subsystem and download a Linux Distro, which is described below. These instructions are current for Microsoft Windows 10 Home, Version 10.0.19041 Build 19041.
Enabling a Linux Subsystem on Windows
Enable the Linux Subsystem by following these steps:
- On a computer running Windows, open “Settings” from the Start Menu
- Click “Apps” | “Apps & features” | “Optional features”
- Click “More Windows features” under the “Related settings” submenu
- Scroll down to “Windows Subsystem for Linux” and check the box
- Click “OK”
- When prompted, restart computer
Second, download a Linux distribution through the Microsoft Store app. We used Ubuntu, published by Canonical Group Limited, and our instructions assume you are using the 20.04 LTS version . Once you have downloaded Ubuntu, open the app and follow the prompts to create a username and password.
It is important to note that files saved through this system will be saved to the Linux Subsystem. It can be somewhat difficult to locate these files if you are not familiar with the program. To find where your files will be saved, open the File Explorer and follow these steps:
- Navigate to %localappdata% in the File Explorer search bar
- Open the “Packages” folder
- Housed within should be the folder: “CanonicalGroupLimited.UbuntuonWindows_79rhkp1fndgsc” or the equivalent for your Linux distro
- Within this particular folder, navigate down through these subfolders:
- – “LocalState” - “rootfs” - “home” - Username
This Username folder is where you must store and place files for your Linux subsystem to read and edit. It is also where files you download using Linux will be stored.
Downloading .las Files with Wget
To download the .las files using this Linux subsystem, we used the wget command. This command tells the computer to download a file from a specific URL or to download multiple files from a list of URLs in a text file. For this dataset, we copied the column from Excel containing the URLs of .las files for each SGL and pasted it into a text document. There was a single URL per line— this is critical. Blank lines separated each SGL into clusters within the text document due to Excel formatting, but this did not impact the functionality of using wget.
The wget command was formatted with several modifiers. To run wget as an administrator, which is necessary to execute the command, you must use sudo which means “superuser do.” Using the modifier -i allowed us to specify a text file to be read, and using the modifier -P allowed us to specify where the downloads should be saved. The final command executed in Ubuntu is:
sudo wget -i SampleList.txt -P sample_file_folder
In this example, SampleList.txt is the name of the text file containing the URLs and sample_file_folder is the name of the folder where downloads should be saved. Both this file and this folder should be in the ‘Username’ folder discussed earlier prior to running the command.
After pressing enter to run the command, wget will download the files from the list and save them in the specified folder. The Linux app will display progress and the files will appear in the specified folder.
This process may take some time, and we left wget running overnight. Using this process, over 4,000 LAS files for a total of about 350 gigabytes of data in zip files was downloaded. The files contained the LiDAR data for every tile touching a State Game Land, or their buffer, in Pennsylvania.
Although we had downloaded all of the data, it was all stored within a single folder. Separating these files into separate folders by SGL was accomplished with a program (ProjectLambda) created by Moritz “Moe” Schiesser. Moe Schiesser’s program read an Excel sheet containing two columns – Column A: the current file location, Column B: desired file location – and copied each file to the appropriate folder. (See Table 1). Note that the file structure in Column B had been created prior to running this program.
|F:\New LAS Files\26001780PAS.las||E:\Pennsylvania Iron and Charcoal\LAS Files\73\26001780PAS.las|
Processing LiDAR Data
The .las files, which were now separated into folders by SGL, were processed through LAStools into three LiDAR derivatives: a digital elevation model (DEM), hillshade, and slope analysis . This process was automated using another program (ProjectKappa) written by Schiesser. This program read the default file names for Pennsylvania’s .las files to determine which coordinate reference system (CRS) to use (State Plane 83 PA South or State Plane 83 PA North) and saved the three LiDAR derivative outputs to the same SGL folder. All tiles for each SGL were combined into one image (for each derivative) in TIF format. For those SGLs that contained portions in both the North and South CRS, it processed the files for each CRS separately, which means some SGLs had two DEMs, hillshades, and slope analyses. In the end, the total dataset, including all .las files and LiDAR derivatives was 1.3 terabytes.
The project studied all SGLs in Pennsylvania. We added a 1km buffer outside of the SGL boundary for each file in order to ensure objects on the edges of SGLs were fully visible. All LAS tiles of which any portion fell within the SGL or buffer were included.
We wrote a program to load each file into QGIS according to the list of SGLs in order to validate each file was present and opened correctly. An example of this program loads the layers Northeast Pennsylvania: https://github.com/jeffblackadar/charcoalhearths/blob/master/qgis_3_load_rasters_north_east.py.
The files for SGLs 264 and 258 could not be processed due to problems in the source files.
(3) Dataset description
Since the dataset contains multiple objects, filenames are displayed in Zenodo.
Primary data: Digital Elevation Model, Slope, Hillshade files. Secondary data: SGL boundaries.
Format names and versions
TIFF, Shapefile, Geojson
01/01/2020 – 31/08/2020
Weston Conner, archaeologist, author, data collection and verification, GIS expert, Lehigh University. https://orcid.org/0000-0001-9906-3762
Moritz “Moe” Schiesser, Programmer.
- State Game Lands with a 1 km buffer
- Grid tiles with SGL information included
- Description for downloading LiDAR tiles (in LAS format): https://zenodo.org/record/4596905
- Program for arranging files downloaded with the above, known as ProjectLambda: https://zenodo.org/record/4572993
- Program for processing original LiDAR data (in LAS format) to yield Digital Elevation Model, Hillshade and Slope, known as ProjectKappa: https://zenodo.org/record/4573004
- Files resulting from the above for State Game Lands (in GeoTIFF format) 
- SGL 12–50: https://zenodo.org/record/4596903
- SGL 51–80: https://zenodo.org/record/4597896
- SGL 81–110: https://zenodo.org/record/4597943
- SGL 111–160: https://zenodo.org/record/4597941
- SGL 161–200: https://zenodo.org/record/4597994
- SGL 201–250: https://zenodo.org/record/4598045
- SGL 251–260: https://zenodo.org/record/4598089
- SGL 261–280: https://zenodo.org/record/4598096
- SGL 281–335: https://zenodo.org/record/4598100
(4) Reuse potential
This dataset is useful for further study of Pennsylvania’s SGLs. Generally, SGLs are relatively undeveloped today, preserving archaeological evidence such as RCHs, building foundations and roadbeds. The Shapefile and GeoTIFFs may be used for further study of the SGLs for RCHs or other types of objects. The method and programs are useful for the collection and processing of LiDAR files to study of areas of Pennsylvania outside of the SGLs. The process may also be adapted for use with LiDAR datasets for geographic areas outside of Pennsylvania.