mkde (version 0.1)

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)

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.

Value

No value is returned

Details

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

Examples

library(raster)
data(condor)
condor <- condor[1:10,] # 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)

data(condordem120)
cell.sz <- mean(res(condordem120))
ext <- extent(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")