Write a survey plan to a GPX file.

Write a survey plan to a GPX file.

WriteSurveyPlan(SurveyPlan, Format = c("GPX", "Ruter"))

Arguments

SurveyPlan

The SurveyPlan process data.

Format

Character: The format of the output file from plotting functions, defaulted to "png". Other options are defined by ggsave, including "tiff", "jpeg" and "pdf".

Value

An object of StoX data type WriteSurveyPlanData.

Examples

library(sf)
#> Linking to GEOS 3.13.0, GDAL 3.8.5, PROJ 9.5.1; sf_use_s2() is TRUE

stratumFile <- system.file(
  "testresources", 
  "strata_sandeel_2020_firstCoverage.wkt", package = "RstoxBase"
 )
stratumPolygon<- DefineStratumPolygon(
  DefinitionMethod = "ResourceFile", 
  FileName = stratumFile
)

# Harbitz zigzag survey design along each stratum:
surveyPlanZZ_Along <- DefineSurveyPlan(
  DefinitionMethod = "ZigZagRectangularEnclosure", 
  StratumPolygon = stratumPolygon, 
  SurveyTime = 200, 
  SurveySpeed = 10, 
  Seed = 1, 
  Bearing = "Along"
)
#> StoX: SurveyPlan for Stratum AlbjoernLing...
#> StoX: numIter: 1
#> StoX: numIter: 2
#> StoX: SurveyPlan for Stratum Engelsk_Klondyke_2020...
#> StoX: numIter: 1
#> StoX: numIter: 2
#> StoX: SurveyPlan for Stratum Inner_Shoal_East_2016...
#> StoX: SurveyPlan for Stratum Inner_Shoal_North_2020...
#> StoX: numIter: 1
#> StoX: numIter: 2
#> StoX: SurveyPlan for Stratum Inner_Shoal_West_2018...
#> StoX: numIter: 1
#> StoX: numIter: 2
#> StoX: SurveyPlan for Stratum Nordgyden...
#> StoX: numIter: 1
#> StoX: numIter: 2
#> StoX: SurveyPlan for Stratum Ostbanken_2020...
#> StoX: numIter: 1
#> StoX: numIter: 2
#> StoX: SurveyPlan for Stratum Outer_Shoal_2020_1...
#> StoX: SurveyPlan for Stratum Vestbanken_North_2020...
#> StoX: numIter: 1
#> StoX: numIter: 2
#> StoX: SurveyPlan for Stratum VestbankenSouthEast...
#> StoX: SurveyPlan for Stratum VestbankenSouthWest...
#> StoX: SurveyPlan for Stratum Vikingbanken...
#> StoX: numIter: 1
#> StoX: numIter: 2

# Convert the survey plan to an sf object and write this as a gpx file 
# (this is done automatically by RstoxFramework in StoX):
gpxData <- WriteSurveyPlan(surveyPlanZZ_Along)
filePath <- tempfile(fileext = ".gpx")
st_write(
     gpxData,
     dsn = filePath,
     layer = "track_points",
     driver = "GPX"
 )
#> writing: substituting ENGCRS["Undefined Cartesian SRS with unknown unit"] for missing CRS
#> Writing layer `track_points' to data source 
#>   `/var/folders/x7/ch5v91h56_zbvbd1y2f600dm0000gn/T//RtmpTYC7Zm/file4109300da0d9.gpx' using driver `GPX'
#> Writing 354 features with 4 fields and geometry type Point.
 
 # Read the data back in:
 gpxData_backin <- st_read(filePath)
#> Multiple layers are present in data source /private/var/folders/x7/ch5v91h56_zbvbd1y2f600dm0000gn/T/RtmpTYC7Zm/file4109300da0d9.gpx, reading layer `waypoints'.
#> Use `st_layers' to list all layer names and their type in a data source.
#> Set the `layer' argument in `st_read' to read a particular layer.
#> Warning: automatically selected the first layer in a data source containing more than one.
#> Reading layer `waypoints' from data source 
#>   `/private/var/folders/x7/ch5v91h56_zbvbd1y2f600dm0000gn/T/RtmpTYC7Zm/file4109300da0d9.gpx' 
#>   using driver `GPX'
#> Simple feature collection with 0 features and 23 fields
#> Bounding box:  xmin: NA ymin: NA xmax: NA ymax: NA
#> Geodetic CRS:  WGS 84
 identical(gpxData, gpxData_backin)
#> [1] FALSE