R/feature.R

In GreKro/cloudcompare: R interface to CloudCompare command line mode

#' Compute Geometric Feature
#'
#' @param file a
#' @param type a
#' @param radius a
#' @param output_dir a
#' @param global_shift a
#' @param global_shift_type a
#' @param filter_sf a
#' @param filter_value a
#' @param c_export_fmt a
#' @param c_ext a
#' @param silent a
#' @param no_timestamp a
#'
#' @return
#' @export
#'
#' @examples
#' \dontrun{
#' #' pc_files <- list.files('C://', pattern='.las$',full.names=T)
#' CC(merge_clouds(files= pc_files))
#' }
#'
#'
#'
feature <- function(file, type='SUM_OF_EIGENVALUES', radius=0.01, output_dir='C:\\output.las',
                    global_shift = F, global_shift_type= 'AUTO', filter_sf= F, filter_value=c(0,1),
                    c_export_fmt= 'LAS', c_ext= 'las', silent= T, no_timestamp= T) {


  if(class(file) != 'character' ){
    stop(paste('Wrong file is a ', class(file),', must be a character'))

  }

  if(!type %in% c('SUM_OF_EIGENVALUES','OMNIVARIANCE','EIGENTROPY',  'ANISOTROPY', 'PLANARITY','LINEARITY','PCA1', 'PCA2', 'SURFACE_VARIATION', 'SPHERICITY', 'VERTICALITY')){
    stop(paste('Wrong type: ',type,' (',class(type),'), must be character value, one of these : SUM_OF_EIGENVALUES, OMNIVARIANCE, EIGENTROPY, ANISOTROPY, PLANARITY, LINEARITY, PCA1, PCA2, SURFACE_VARIATION, SPHERICITY, VERTICALITY',sep=''))

  }


  if(class(radius) != 'numeric' ){
    stop(paste('Wrong nop is a ', class(file),', must be a numeric'))

  }



  if(class(filter_sf) != "logical" ){
    stop('Wrong filter_sf must be a logical value - TRUE or FALSE')
  }


  if(length(filter_value) != 2){
    stop('Wrong filter_value must be vector with 2 values for minimum and maximum treshold in filter_sf tool')
  }


  if(class(global_shift) != "logical" ){
    stop('Wrong global_shift must be a logical value - TRUE or FALSE')
  }


  if(!(global_shift_type %in% c('AUTO', 'FIRST')) && (class(global_shift_type) != 'numeric' || length(global_shift_type) != 3)){
    stop('Wrong global_shift_type must be character value : "AUTO" or "FIRST" (avaliable since CC v.2.11) or numeric vector with 3 values for dimensions x,y,z ')
  }

  if(class(no_timestamp) != "logical" ){
    stop('Wrong no_timestamp must be a logical value - TRUE or FALSE')
  }

  if(!c_export_fmt %in% c('ASC','BIN', 'PLY', 'LAS', 'E57', 'VTK', 'PCD', 'SOI', 'PN', 'PV')){
    stop(paste('Wrong c_export_fmt: ',c_export_fmt,' (',class(c_export_fmt),'), must be character value, one of these : ASC, BIN, PLY, LAS, E57, VTK, PCD, SOI, PN, PV',sep=''))

  }

  if( (!c_ext %in% c('ASC','BIN', 'PLY', 'LAS', 'LAZ', 'E57', 'VTK', 'PCD', 'SOI', 'PN', 'PV') && !c_ext %in% tolower(c('ASC','BIN', 'PLY', 'LAS','LAZ', 'E57', 'VTK', 'PCD', 'SOI', 'PN', 'PV'))) ){
    stop(paste('Wrong c_ext: ',c_ext,' (',class(c_ext),'),
               must be character value, one of these :
               ASC,BIN, PLY, LAS, LAZ, E57, VTK, PCD, SOI, PN, PV
               or
               asc, bin, ply, las, laz, ,e57, ,vtk, pcd, soi, pn, pv',sep=''))

  }


  auto_save_off <- '-AUTO_SAVE OFF'

  if (silent==T){
    s1= '-SILENT'

  }else{
    s1= ''
  }


  if (no_timestamp== T) {
    s2= '-NO_TIMESTAMP'
  } else{
    s2=''
  }

  s3= paste('-C_EXPORT_FMT',c_export_fmt,'-EXT',c_ext,sep=' ')


  if (global_shift == T){

    s4= paste('-O -GLOBAL_SHIFT',global_shift_type, file,sep=' ')

  } else {

    s4= paste('-O',file,sep=' ')
  }


  s5= paste('-FEATURE',type, radius, sep=' ')

  if (filter_sf == T) {
    s6= paste('-FILTER_SF',filter_value[1], filter_value[2], sep=' ')
  } else{
    s6=''
  }


  s7= paste('-SAVE_CLOUDS' , 'FILE', output_dir, sep=' ')




  cc_syntax <- paste(s1,s2,auto_save_off,s3,s4,s5,s6,s7,sep=' ')

  return(cc_syntax)
}