data-raw/MakeDatafor_abspatial.r

#rm(list=ls())   # Cleanup the R console if required
# Set up the run ----------


options("show.signif.stars"=FALSE,"stringsAsFactors"=FALSE,
        "max.print"=50000,"width"=240)

datadir <- "C:/A_Mal/Rcode/Abalone/abspatial/data-raw/"
#  listFunctions("C:/A_Mal/Rcode/Abalone/abspatial/R/spatial_utils.r")

# taspoly -----------
# now called trawldeep_scalefish
filename="tas_100k.csv"
taspoly <- read.csv(paste0(datadir,filename),header=TRUE)


save(taspoly,file=paste0(datadir,"taspoly.rda"))
# actually saved as sysdata.rda inside the R directory


#  check data ---------------
tools::checkRdaFiles(paths=datadir)
tools::resaveRdaFiles(paths=datadir,compress="auto")
tools::checkRdaFiles(paths=datadir)


# make abdat.Rdata ---------------------------------------------------------

filename <- "C:/A_Mal/Rcode/AbGeo026/geostat/G1Ha_2019_05_13.rds"
ab <- readRDS(filename)#get spatial obj
abdat <- getlatlong(ab)     # convert to data.frame + lat long
props <- properties(abdat)  # characterize each field in data.frame

repeatoid <- apply(abdat[,5:11],1,countgtzero)
abdat$occur <- NA
abdat$occur <- repeatoid
save(abdat,file=paste0(datadir,"abdat.RData"))
haddonm/abspatial documentation built on June 7, 2019, 9:54 a.m.