Write MKDE to a GRASS GIS 3D ASCII raster file.
Write the MKDE to a VTK file.
writeToGRASS(mkde.obj, fname, nodat="NA", cumprob=FALSE)
mkde.obj |
3D MKDE object created with |
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. |
This function writes a GRASS GIS ASCII raster file that can be
imported using the r3.in.ascii function.
No value is returned
Jeff A. Tracey, PhD
USGS Western Ecological Research Center, San Diego Field Station
jatracey@usgs.gov
James Sheppard, PhD
San Diego Zoo Institute for Conservation Research
jsheppard@sandiegozoo.org
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")
Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.