cloudSampleR: Clip random matched locations from multiple point cloud...

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

Clip random matched locations from multiple point cloud datasets (lidar or dap or ifsar etc.)

Description

<Delete and Replace>

Usage

cloud2xSample(
  pathClipData = "c:/fusion/clipdata.exe",
  pathOutA = "",
  pathOutB = "",
  pathLasA = "",
  pathLasB = NA,
  patternA = ".*[.]la[s|z]{1}$",
  patternB = ".*[.]la[s|z]{1}$",
  pathDTMA = NA,
  pathDTMB = NA,
  patternDTMA = ".*[.]asc$",
  patternDTMB = ".*[.]asc$",
  extentPolyA = NA,
  extentPolyB = NA,
  sampleShpA = NA,
  proj4A = NA,
  proj4B = NA,
  zMult = c(1, 0.3048, 1/0.3048)[1],
  extentSample = c(llx = NA, lly = NA, ulx = NA, uly = NA),
  nSample = 100,
  radii = list(feet = c(FtTenthAc = 37.2, FtAcre = 117.8, Ft5Acres = 263.3), meters =
    c(MtenthAc = 11.3, MAcre = 35.9, M5Acres = 80.3))[[1]],
  sampleShape = c("circle", "square", "userPoly")[1],
  sampleType = c("regular", "random", "hexagonal"),
  doSteps = c("sample", "clip", "metrics"),
  switchesClipdata = "",
  switchesCloudmetrics = "",
  procMethod = c("lidR", "FUSION"),
  nCore = 2,
  temp = "c:\\temp\\clipdata"
)

Arguments

pathClipData	where to find FUSION clipdata function - only used with procMethod = "FUSION"
pathOutA	directory to place clipped plots from project A
pathOutB	directory to place clipped plots from project B
pathLasA	directory containing project A las or laz files
pathLasB	directory containing project B las or laz files
patternA	grep pattern to match files on for project A
patternB	grep pattern to match files on for project B
pathDTMA	(Semi-functional with procMethod = "FUSION", NOT TESTED YET) DTM to use with project A to normalize point heights
pathDTMB	(Semi-functional with procMethod = "FUSION", NOT TESTED YET) DTM to use with project B to normalize point heights
extentPolyA	polygon of project A extent
extentPolyB	polygon of project B extent
sampleShpA	(Optional) use an existing shapefile with plot locations
proj4A	(Optional) Provide proj4 strings if projects A and B are in different projections but do not contain projection information (e.g. missing .prj files)
proj4B	(Optional) Provide proj4 strings if projects A and B are in different projections but do not contain projection information (e.g. missing .prj files)
zMult	specify multiplier for vertical units to go e.g. from feet to meters
extentSample	(Optional) If the overlap between projects A and B is known, you can provide it here
nSample	number of plots to sample
radii	a vector of plot radii to use in clipping plots
sampleShape	shape of clipped plots ("circle" or "square")
sampleType	type of spatial distribution of of sample points ("regular" = square grid,"random","hexagonal" = hexagonal grid)
doSteps	(NON FUNCTIONING) Select which steps should be completed in this run
switchesClipdata	optional switches to FUSION's cloudmetrics.exe - see fusion documentation for options
switchesCloudmetrics	(NON FUNCTIONING)
procMethod	Should plots be clipped using R or FUSION?
nCore	How many cores should be used for clipping (the number and speed of hdds is as important as the number of machine threads)
temp	(use with procMethod = "FUSION")Where should batch files be placed
pathCloudmetrics	(NON FUNCTIONING)

Details

<Delete and Replace>

Revision History

1.0	2019 May 07 First Debug Version
1.0	2019 May 09 Add Documentation for arguments

Value

<Delete and Replace>

Author(s)

Jacob Strunk <[e-mail]>

See Also

readLAS
spsample

Examples

		library(rgdal)
		library(rgeos)
		library(cloudSampleR)

		#A read in tiles, prepare extents of data
		poly1=rgdal::readOGR("D:/data/wadnr_hood_canal/las/hood_canal_3in_DSM_2015/manage_las","las_polys")
		poly1a=gBuffer(poly1,width = 5)
		poly2=rgdal::readOGR("D:/data/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/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 = "FUSION"
			,procMethod = "lidR"
			,radii = list( feet = c(FtTenthAc = 37.2, FtAcre = 117.8, Ft5Acres = 263.3 ))
		)

jstrunk001/RSForInvt documentation built on April 18, 2022, 11:03 p.m.