create_tracking_history: Wrapper function to strt the JPDAF algorithm

In ick003/vesselett: A suite of target tracking algorithm, with a specialization for marine radar

Description

Wrapper function to strt the JPDAF algorithm

Usage

create_tracking_history(
  n = 2,
  nscans = 30,
  delta_time = 1,
  starting_position = c(0, 0),
  starting_speed = c(0, 0),
  measure_rate = 1,
  target_type = "point",
  with_radar = TRUE,
  prob_missed = 0.1,
  prob_clutter = 0.1,
  clutter_intensity = 0,
  scan_range = 200,
  turning_rate = 0.1,
  prob_ext = 0.9,
  kinematic_parameters
)

Arguments

n	integer: number of targets
nscans	integer: number of scans by the radar
delta_time	scalar: time between two scans
starting_position	vector: location for the initial position. If multiple targets, locations are randomly chosen around the vector
starting_speed	vector: speed for the target birth.
measure_rate	scalar: Average number of measure generated by one target for a given scan.
target_type	string: at the moment, only point or extended.
with_radar	logical: create a trajectory for the radar
prob_missed	numeric: Probability of target not being detected, between 0 and 1
prob_clutter	numeric: Probability of false target being detected, between 0 and 1
clutter_intensity	integer: Number of potential false targets
scan_range	scalar: Extent of the radar scan, at each time.
turning_rate	scale:
prob_ext	numeric: Probability of target vanishing
kinematic_parameters	list:

Value

If type is point, data.frame with time, utm north and east (in meters) and unique target identifier per scan. If extended, returns values for speed, turning rate, size and shape.

Examples

## Not run: 
library(ggplot2)
library(ggforce)
library(tidyverse)

## Example 1
x <- create_tracking_history()
ggplot(x$obs_target, aes(x = east, y = north, col = target_id)) +
geom_point(size = 0.5) +
geom_path(data = x$obs_target %>% dplyr::filter(.data$target_id != "0"))

## Example 2
 x <- create_tracking_history(n = 1, target_type = "extended", measure_rate = 150)$obs_target
ggplot(data = x) +
geom_point(aes(x = east, y = north, col = time), size = 1) +
theme_bw()


## End(Not run)

ick003/vesselett documentation built on July 20, 2020, 9:08 p.m.