_read_me.txt file for Zenodo repository Author: Billy Browning (bbrowning@nhcweb.com) 2024/10/12 Satellite imagery by PlanetScope used in the satellite validation section is not open source and not provided in the Zenodo repository. Therefore, input data unavailable for SAT_VAL scripts. However, satellite imagery metadata are available in Supplementary Material section. ________________________________________________________________ Input Data Extraction ________________________________________________________________ The following section provides instruction or how to download the CDEM and CanVec water features data in the Geogratis website, and how to download the LiDAR point cloud and DSM data from BC LiDAR portal for the Lillooet Stream Study Reach. Canadian digital elevation model (CDEM) and CanVec water features data download 1) Open the Geospatial Data Extraction webpage: https://maps.canada.ca/czs/index-en.html 2) In the 'select data to be extracted' tab select 'CanVec' to extract the water features shapefiles.. 3) In the 'Select clipping area' tab click 'Area from a shapefile' and choose a clipping extent from file such as 'geogratis_extent.zip' file that contains the shapefile. 4) In the 'Select options and submit job' tab, select 'Lakes and rivers - Hydrographic features', 'ESRI Shapefile' from 'Output format choice'. 5) Enter email address and the CanVec water features shapefiles will be sent to your inbox. 'watercourse_1.shp' provides line features of water courses and 'waterbody_2' provides polygon features of large water surfaces. 6) Similarly, to download the CDEM data, select 'Elevation' in the 'Select data to be extracted' tab. 7) In the 'Select clipping area' tab click 'Area from a shapefile' and choose a clipping extent from file such as 'geogratis_extent.zip' file that contains the shapefile. 8) In the 'Select options and submit job' tab, select 'Canadian Digital Elevation Model' in the 'Elevation product type' drop down and tick 'Digital Elevation Model (DEM). 5) Enter email address and the CDEM data will be sent to your inbox. LiDAR data from BC LiDAR portal (Point Cloud and DSM). 1) Open the BC LiDAR data download portal webpage: https://lidar.gov.bc.ca/pages/download-discovery 2) In the top left corner, in the search bar type in your area of interest - e.g. 'Pemberton'. 3) To download the LiDAR point cloud data; in the 'Layer List' tab, only select 'Point cloud Index - 1:2,500 Grid'. 4) In the 'Select' tab, ensure only 'Point cloud Index - 1:2,500 Grid' is selected, click 'Select'. 5) Use the cursor to draw a box around the tiles that you require. In the case of the Lillooet River, draw a box that is from the most northern grid box along the Lillooet River south towards the northern tip of Lillooet Lake. In some cases you may need to hold shift and create multiple boxes, to select specific regions. 6) Once the appropriate grids have been highlighted, click on the horizontal ellipsis button to the right of the 'Point cloud Index - 1:2,500 Grid' layer in the 'Select tab'and click 'View in Attribute Table. 7) This now brings up the attribute table of the selected tiles. Go through row by row clicking on the link in the 'Download Link' column to download each tile. Each tile must be downloaded individually. 8) For this project, save the downloaded .laz files in the following pathway: 'MyProjectFolder/data/LiDAR/lax_catalog/Lillooet'. 9) To download the DSM tiles, following steps 1 to 7, but ensure that you only select the 'DEM Index - 1:20,000 Grid' layer. 10) For this project, save the downloaded DSM.tif files in the following pathway: 'MyProjectFolder/data/LiDAR/lax_catalog/Lillooet'. A list of tiles used to download the LiDAR Point Cloud and DSM data are listed at the end of this _read_me.txt file. ________________________________________________________________ .R script file structure ________________________________________________________________ Create an R project in a project folder, then create folders named code, data, figures, and output for the project. The following files should be copied to the code folder. - 00_functions.R: Contains the relevant functions including creating appropriately spaces points of interest (POIs), Calculating maximum and minimum horizon angles, calculating shade and view factors, among others. - 01_data_prep.R: Prepping data inputs such as CDEM, LiDAR point cloud and LiDAR DSM. Delineating stream surface, create POIs. - 02.1_SAT_VAL_imagery_retrieval_and_shade_model.R: Satellite validation - retrieve imagery and run shade model on image timesteps. - 02.2_SAT_VAL_satellite_validation.R: analysis for satellite validation. - 02.3_SAT_VAL_results.R: results for satellite validation. - 04.2_ROTATE_run_shade_model.R: calculate horizon angles, run shade model and calculate view factors for Lillooet Study reach. - 05.1_RESULTS.R: results of shade model outputs. ________________________________________________________________ LiDAR tile list ________________________________________________________________ 'Point cloud Index - 1:2,500 Grid' TILES bc_092j036_2_4_2_xyes_8_utm10_20170601.laz bc_092j036_2_4_3_xyes_8_utm10_20170601.laz bc_092j036_2_4_4_xyes_8_utm10_20170601.laz bc_092j036_4_2_1_xyes_8_utm10_20170601.laz bc_092j036_4_2_2_xyes_8_utm10_20170601.laz bc_092j036_4_2_3_xyes_8_utm10_20170601.laz bc_092j036_4_3_2_xyes_8_utm10_20170601.laz bc_092j036_4_3_3_xyes_8_utm10_20170601.laz bc_092j036_4_3_4_xyes_8_utm10_20170601.laz bc_092j036_4_4_1_xyes_8_utm10_20170601.laz bc_092j037_1_1_3_xyes_8_utm10_20170601.laz bc_092j037_1_1_4_xyes_8_utm10_20170601.laz bc_092j037_1_3_1_xyes_8_utm10_20170601.laz bc_092j037_1_3_2_xyes_8_utm10_20170601.laz bc_092j046_1_2_4_xyes_8_utm10_20170601.laz bc_092j046_1_4_2_xyes_8_utm10_20170601.laz bc_092j046_1_4_4_xyes_8_utm10_20170601.laz bc_092j046_2_1_1_xyes_8_utm10_20170601.laz bc_092j046_2_1_2_xyes_8_utm10_20170601.laz bc_092j046_2_1_3_xyes_8_utm10_20170601.laz bc_092j046_2_3_1_xyes_8_utm10_20170601.laz bc_092j046_3_1_4_xyes_8_utm10_20170601.laz bc_092j046_3_2_1_xyes_8_utm10_20170601.laz bc_092j046_3_2_2_xyes_8_utm10_20170601.laz bc_092j046_3_2_3_xyes_8_utm10_20170601.laz bc_092j046_3_2_4_xyes_8_utm10_20170601.laz bc_092j046_3_3_2_xyes_8_utm10_20170601.laz bc_092j046_3_3_3_xyes_8_utm10_20170601.laz bc_092j046_3_3_4_xyes_8_utm10_20170601.laz bc_092j046_3_4_1_xyes_8_utm10_20170601.laz bc_092j046_3_4_3_xyes_8_utm10_20170601.laz bc_092j055_2_2_2_xyes_8_utm10_20170601.laz bc_092j055_2_2_3_xyes_8_utm10_20170601.laz bc_092j055_2_2_4_xyes_8_utm10_20170601.laz bc_092j056_1_1_1_xyes_8_utm10_20170601.laz bc_092j056_1_1_2_xyes_8_utm10_20170601.laz bc_092j056_1_1_3_xyes_8_utm10_20170601.laz 'DEM Index - 1:20,000 Grid' TILES bc_092j036_xl1m_utm10_170601.tif bc_092j037_xl1m_utm10_170601.tif bc_092j045_xl1m_utm10_170601.tif bc_092j046_xl1m_utm10_170601.tif bc_092j055_xl1m_utm10_170601.tif bc_092j056_xl1m_utm10_170601.tif