mkde (version 0.4)

writeToGRASS: Write MKDE to a GRASS GIS 3D ASCII raster file.

Description

Write the MKDE to a VTK file.

Usage

writeToGRASS(mkde.obj, fname, nodat="NA", cumprob=FALSE)

Value

No value is returned

Arguments

mkde.obj

3D MKDE object created with initialize3DMKDE and density initialized with initializeDensity

fname

The patch and file name for output VTK file

nodat

A no data character string that will be written in place of no data values.

cumprob

Indicate whether to write the voxel probabilities of cumulative probabilities.

Author

Jeff A. Tracey, PhD
USGS Western Ecological Research Center, San Diego Field Station

James Sheppard, PhD
San Diego Zoo Institute for Conservation Research

Details

This function writes a GRASS GIS ASCII raster file that can be imported using the r3.in.ascii function.

Examples

library(terra)
data(condor)
condor <- condor[1:20,] # simply to make example run more quickly
mv.dat <- initializeMovementData(condor$time, condor$x, condor$y, 
z.obs=condor$z, sig2obs=25.0, sig2obs.z=81.0, t.max=65.0)

fpath <- system.file("extdata", "condordem120.RDS", package="mkde")
condordem120 <- terra::readRDS(fpath)
cell.sz <- mean(res(condordem120))
ext <- ext(condordem120)
nx <- ncol(condordem120)
ny <- nrow(condordem120)
mkde.obj <- initializeMKDE3D(ext$xmin, cell.sz, nx, ext$ymin, cell.sz,
ny, min(values(condordem120), na.rm=TRUE), cell.sz, 25)

# note: we use a raster coarse integration time step so the
# example runs faster
dens.res <- initializeDensity(mkde.obj, mv.dat, integration.step=10.0)
mkde.obj <- dens.res$mkde.obj
mv.dat <- dens.res$move.dat

# Write file
writeToGRASS(mkde.obj, "ascii3d.txt")

# Clean up files
unlink("ascii3d.txt")