writeToVTK: Write MKDE to an ASCII Visualization Tool Kit (VTK)...

In mkde: 2D and 3D Movement-Based Kernel Density Estimates (MKDEs)

Write MKDE to an ASCII Visualization Tool Kit (VTK) structured grid file.

Description

Write the MKDE to a VTK file.

Usage

writeToVTK(mkde.obj, fname, description="3D MKDE", 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
description	A character string with a brief description that will be placed in the header of the VTK file
cumprob	Indicate whether to write the voxel probabilities of cumulative probabilities.

Details

This function writes a VTK structured grid file for a 3D MKDE that can be used in 3D visualization tools such as MayaVI or ParaView

Value

No value is returned

Author(s)

Jeff A. Tracey, PhD
USGS Western Ecological Research Center, San Diego Field Station
[e-mail]
James Sheppard, PhD
San Diego Zoo Institute for Conservation Research
[e-mail]

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)
# next two lines reduce extent of 2D space to speed execution of example
tmp <- ext(c(range(condor$x) + c(-100, 100), range(condor$y) + c(-100, 100)))
condordem120 <- crop(condordem120, tmp)

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

# Save as VTK file
writeToVTK(mkde.obj, "condor_3dMKDE.vtk",
description="Example California condor 3D MKDE")

# Clean up files
unlink("condor_3dMKDE.vtk")

mkde documentation built on April 3, 2025, 11:44 p.m.