inst/examples/example_cloud2xSample.R

In jstrunk001/RSForInvt: Remote Sensing FORest INVenTory related analysis tasks, especially FUSION, lidR, NVEL

#install packages (only run 1 time, or when package is updated)
#installing dependencies takes a little while for some reason
if(F){
	#this approach requires that you first install rtools.exe
	devtools::install_github("jstrunk001/CloudSampleR")

	#alternatively, you can download the github repository, zip it, and then use R to load package from zip file
}

#load required packages
	library(rgdal)
	library(rgeos)
	library(cloudSampleR)

#Optional - build a single .vrt raster from a list of many DTM rasters (e.g. .tif and .img)
# this simplifies managing rasters with the raster package / lidR package
if(F){

	tifs = list.files("D:\\data\\wa\\wadnr_hood_canal\\hood_canal_dtm",pattern="[.]tif",full.names=T)
	gdalUtils::gdalbuildvrt(tifs,"D:\\data\\wadnr_hood_canal\\hood_canal_dtm\\dtm_demo.vrt")

}

#A read in tiles, prepare extents of data
	poly1=rgdal::readOGR("D:/data/wa/wadnr_hood_canal/las/hood_canal_3in_DSM_2015/manage_las","las_polys")
	poly1a=gBuffer(poly1,width = 5)
	poly2=rgdal::readOGR("D:/data/wa/wadnr_hood_canal/las/hood_canal_6in_DSM_2015/manage_las","las_polys")
	poly2a=gBuffer(poly2,width = 5)

#B clip plots and send to output folders
	cloud2xSample(
		pathClipData = "c:/fusion/clipdata.exe"
		,pathOutA = "d:/temp/hood_canal_test/clip3in/"
		,pathOutB = "d:/temp/hood_canal_test/clip6in/"
		,pathLasA = "D:/data/wa/wadnr_hood_canal/las/hood_canal_3in_DSM_2015/"
		,pathLasB = "D:/data/wa/wadnr_hood_canal/las/hood_canal_6in_DSM_2015/"
		,extentPolyA = poly1a
		,extentPolyB = poly2a
		,nCore = 3
		,nSample = 150
		#,procMethod = "FUSION"
		,procMethod = "lidR"
		,radii = list( feet = c(FtTenthAc = 37.2, FtAcre = 117.8, Ft5Acres = 263.3 ))
		#,radii = list( feet = c(FtTenthAc = 37.2 ))
	)

	#single buffer and bring in existing points - process with R
	cloud2xSample(
		pathClipData = "c:/fusion/clipdata.exe"
		,pathOutA = "d:/temp/hood_canal_test/clip3in/"
		,pathOutB = "d:/temp/hood_canal_test/clip6in/"
		,pathLasA = "D:/data/wadnr_hood_canal/las/hood_canal_3in_DSM_2015/"
		,pathLasB = "D:/data/wadnr_hood_canal/las/hood_canal_6in_DSM_2015/"
		,extentPolyA = poly1a
		,extentPolyB = poly2a
		,nCore = 3
		,nSample = 150
		,procMethod = "lidR"
		,radii = list( feet = c(FtTenthAc = 37.2 ))
		,sampleShpA = "D:\\temp\\hood_canal_test\\clip3in\\shapefiles\\2019Jun10165355_SamplePoints.shp"
	)

#some potential features to add:
	# reproject clips
	# compute cloudmetrics
	# merge cloudmetrics

#C compute cloudmetrics for samples and merge ?