diff --git a/dev/CODE_OF_CONDUCT.html b/dev/CODE_OF_CONDUCT.html index c1734b97..93dbb72f 100644 --- a/dev/CODE_OF_CONDUCT.html +++ b/dev/CODE_OF_CONDUCT.html @@ -17,7 +17,7 @@
diff --git a/dev/CONTRIBUTING.html b/dev/CONTRIBUTING.html index 4298dd01..16cf6da1 100644 --- a/dev/CONTRIBUTING.html +++ b/dev/CONTRIBUTING.html @@ -17,7 +17,7 @@ diff --git a/dev/LICENSE-text.html b/dev/LICENSE-text.html index 8656e707..a7164f5a 100644 --- a/dev/LICENSE-text.html +++ b/dev/LICENSE-text.html @@ -17,7 +17,7 @@ diff --git a/dev/articles/bioRad.html b/dev/articles/bioRad.html index 6c90108c..fb7bc228 100644 --- a/dev/articles/bioRad.html +++ b/dev/articles/bioRad.html @@ -35,7 +35,7 @@ diff --git a/dev/articles/index.html b/dev/articles/index.html index 4fe1c8fe..a023899e 100644 --- a/dev/articles/index.html +++ b/dev/articles/index.html @@ -17,7 +17,7 @@ diff --git a/dev/articles/rad_aero_22.html b/dev/articles/rad_aero_22.html index 910a7948..08d584ef 100644 --- a/dev/articles/rad_aero_22.html +++ b/dev/articles/rad_aero_22.html @@ -35,7 +35,7 @@ diff --git a/dev/articles/range_correction.html b/dev/articles/range_correction.html index f20df94d..1def2785 100644 --- a/dev/articles/range_correction.html +++ b/dev/articles/range_correction.html @@ -35,7 +35,7 @@ @@ -111,10 +111,8 @@
options("sp_evolution_status" = 2L)
library(bioRad)
-## Welcome to bioRad version 0.7.3.9000
-## Attempting to load MistNet from:/home/runner/work/_temp/Library/vol2birdR/lib
-## MistNet successfully initialized.
-## using vol2birdR version 1.0.3 (MistNet installed)
+## Welcome to bioRad version 0.8.1.9000
+## using vol2birdR version 1.0.4 (MistNet installed)
+# define a range grid from 0 to 200 km: my_range <- seq(0, 200000, 1000) # plot the beam height for each range, for a 0.5 degree elevation beam: plot(my_range, beam_height(my_range, elev = .5), xlab = "range [m]", ylab = "beam height [m]")
The beam width also broadens with range:
-diff --git a/dev/reference/beam_profile_overlap.html b/dev/reference/beam_profile_overlap.html index c69fe091..7b196444 100644 --- a/dev/reference/beam_profile_overlap.html +++ b/dev/reference/beam_profile_overlap.html @@ -19,7 +19,7 @@+@@ -149,17 +147,17 @@# plot the beam width, for a beam opening angle of 1 degree (typical for most weather radars): plot(my_range, beam_width(my_range, beam_angle = 1), xlab = "range [m]", ylab = "beam width [m]")
2 Beam propagation
+-# plot the beam profile, for a 0.5 degree elevation beam at 50 km distance from the radar: plot(beam_profile(height = 0:4000, 50000, 0.5), 0:4000, xlab = "normalized radiated energy", ylab = "height [m]", main = "beam elevation: 0.5 deg, distance=50km")
diff --git a/dev/reference/beam_profile.html b/dev/reference/beam_profile.html index 98a58ab8..0dbef7fb 100644 --- a/dev/reference/beam_profile.html +++ b/dev/reference/beam_profile.html @@ -20,7 +20,7 @@+# plot the beam profile, for a 2 degree elevation beam at 50 km distance from the radar: plot(beam_profile(height = 0:4000, 50000, 2), 0:4000, xlab = "normalized radiated energy", ylab = "height [m]", main = "beam elevation: 2 deg, distance=50km")
We can also calculate the normalized radiation profile for the two beams combined:
-@@ -129,8 +129,8 @@+@@ -169,7 +167,7 @@# plot the combined beam profile for a 0.5 and 2.0 degree elevation beam at 50 km distance from the radar: plot(beam_profile(height = 0:4000, 50000, c(0.5, 2)), 0:4000, xlab = "normalized radiated energy", ylab = "height [m]", main = "beam elevations: 0.5,2 deg, distance=50km")
3 Vertical profiles and polar volum
Let us now assume we have calculated a vertical profile (of birds) for a certain polar volume file. Let’s load an example
-diff --git a/dev/reference/beam_height.html b/dev/reference/beam_height.html index fcd72c83..f1d4a810 100644 --- a/dev/reference/beam_height.html +++ b/dev/reference/beam_height.html @@ -18,7 +18,7 @@+@@ -122,8 +122,8 @@# let's load an example polar volume: pvolfile <- system.file("extdata", "volume.h5", package = "bioRad") example_pvol <- read_pvolfile(file = pvolfile) @@ -185,13 +183,13 @@
3 Vertical profiles and polar volum ## generated by: vol2bird 0.3.17
Let’s plot the vertical profile, for the quantity eta (the linear reflectivity):
-diff --git a/dev/reference/beam_distance.html b/dev/reference/beam_distance.html index 02900731..0fbc3f3d 100644 --- a/dev/reference/beam_distance.html +++ b/dev/reference/beam_distance.html @@ -18,7 +18,7 @@+plot(example_vp, quantity = "eta")
Note that eta is directly related to reflectivity factor (DBZH), i.e. a reflectivity factor of 5 dBZ amounts to the following eta at a radar wavelength of 5.3 cm:
-@@ -253,7 +251,7 @@+dbz_to_eta(5, wavelength = 5.3)
## [1] 1140.588
4 Spatial estima
5 Example: estimating a VID image
Let’s first make a PPI plot of the lowest uncorrected scan:
-diff --git a/dev/reference/attribute_table.html b/dev/reference/attribute_table.html index 2c08103e..5a3d2911 100644 --- a/dev/reference/attribute_table.html +++ b/dev/reference/attribute_table.html @@ -17,7 +17,7 @@+# extract the first scan from the polar volume: my_scan <- example_pvol$scans[[1]] # project it as a PPI on the ground: @@ -262,7 +260,7 @@
5 Example: estimating a VID imageplot(my_ppi)
Now let’s calculate the range-corrected PPI
-diff --git a/dev/reference/as.vpts.html b/dev/reference/as.vpts.html index a0f368ce..775c54f6 100644 --- a/dev/reference/as.vpts.html +++ b/dev/reference/as.vpts.html @@ -17,7 +17,7 @@+@@ -272,18 +270,20 @@# let's use a 500 metre spatial grid (res), and restrict to 100x100 km area my_corrected_ppi <- integrate_to_ppi(example_pvol, example_vp, res = 500, xlim = c(-100000, 100000), ylim = c(-100000, 100000)) my_corrected_ppi
5 Example: estimating a VID image## dims: 400 x 400 pixels
The range corrected PPI has four parameters: VIR, VID, R, overlap. Let’s plot the adjustment factor R:
-+# plot the adjustment factor R: plot(my_corrected_ppi, param = "R")
Let’s also plot the vertically integrated reflectivity:
-+plot(my_corrected_ppi, param = "VIR")
Or plot the vertically integrated density on a map:
-++if (all(sapply(c("ggspatial","prettymapr", "rosm"), requireNamespace, quietly = TRUE))){ + map(my_corrected_ppi, map=bm, param = "VIR", alpha = .5) +}bm <- "osm" -map(my_corrected_ppi, map=bm, param = "VIR", alpha = .5)
## Warning in CPL_crs_from_input(x): GDAL Message 1: +init=epsg:XXXX syntax is ## deprecated. It might return a CRS with a non-EPSG compliant axis order.
@@ -306,7 +306,7 @@## Zoom: 6
6 Overlap between ra
We can calculate overlap between emitted radiation and the biological profile as follows:
-+@@ -168,8 +168,10 @@# calculate overlap between vertical profile of birds # and the vertical radiation profile emitted by the radar: bpo <- beam_profile_overlap(example_vp, get_elevation_angles(example_pvol), seq(0, 100000, 1000), quantity = "eta") @@ -326,7 +326,7 @@
6 Overlap between ra
The range-corrected PPI also contains a field
-overlap
, which contains the same metric but then calculated on the grid of the PPI:diff --git a/dev/index.html b/dev/index.html index a1771e55..d107a0aa 100644 --- a/dev/index.html +++ b/dev/index.html @@ -37,7 +37,7 @@+plot(my_corrected_ppi, param = "overlap")
The overlap metric might be used in the future to show only certain diff --git a/dev/articles/range_correction_files/figure-html/unnamed-chunk-13-1.png b/dev/articles/range_correction_files/figure-html/unnamed-chunk-13-1.png index ef360eaa..2fcf18f4 100644 Binary files a/dev/articles/range_correction_files/figure-html/unnamed-chunk-13-1.png and b/dev/articles/range_correction_files/figure-html/unnamed-chunk-13-1.png differ diff --git a/dev/authors.html b/dev/authors.html index 1d36e8d8..7d81720d 100644 --- a/dev/authors.html +++ b/dev/authors.html @@ -17,7 +17,7 @@
Install bioRadThen load the package with:
+#> Welcome to bioRad version 0.8.1 +#> Attempting to load MistNet from:/Users/amd427/Library/R/x86_64/4.3/library/vol2birdR/lib +#> MistNet successfully initialized. +#> using vol2birdR version 1.0.3 (MistNet installed)library(bioRad) -#> Welcome to bioRad version 0.7.3 -#> using vol2birdR version 1.0.1 (MistNet installed)
@@ -69,21 +69,37 @@(optional) Enable MistNet diff --git a/dev/news/index.html b/dev/news/index.html index 53e27a97..37462700 100644 --- a/dev/news/index.html +++ b/dev/news/index.html @@ -17,7 +17,7 @@
Changelog
-@@ -195,6 +195,7 @@bioRad 0.7.3.9000
-
- speed up
-integrate_to_ppi()
andproject_as_ppi()
by using nativesf
functions (#669)-Bugfixes
-
- -
support for tidyverse select method for polar volume and polar scan objects (#668,#460)
- +
Updated the default refractive index value used in conversion of linear reflectivity (eta) to logarithmic reflectivity (dBZ) (#642). The effect is a 7% increase in animal densities in output of functions
integrate_to_ppi()
andread_cajun()
only.bioRad 0.8.12024-07-30
+-+bioRad 0.7.32023-10-20
+bioRad 0.8.02024-07-27
-New features
-
- Replaced the
+ggmap
package withggspatial
for map visualizations. This change was made asggmap
no longer provides reliable open-source basemaps without the necessity to register for an API key (#638).New features
+
- +
ENRAM VPTS data exchange format added in package data as
vpts_schema.rda
- +
VPTS files are now able to be validated with
validate_vpts()
which uses the schema to check for min/max constraint violations for specific fields, regex and datetime formatting- +
speed up
integrate_to_ppi()
andproject_as_ppi()
by using nativesf
functions (#669)support for tidyverse select method for polar volume and polar scan objects (#668,#460)
++Bugfixes
+
- +
Corrected the default refractive index value used in conversion of linear reflectivity (eta) to logarithmic reflectivity (dBZ). The effect is a 7% increase in animal densities in output of functions
integrate_to_ppi()
andread_cajun()
only (#642).- +
Fixed the handling of empty numeric vectors when plotting clutter maps (#655)
- +
Fixed warning when reading VPTS csv containing multiple values in
lat
,lon
, ‘rcs
’ orsd_vvp_threshold
(#651)- +
Fixed enabling/disabling of
single_pol
flag incalculate_vp()
(#646)- +
Updated the s3 source bucket of
download_vpfiles()
to https://aloftdata.s3-eu-west-1.amazonaws.com (#648)- +
Fixed an error in the idw method of
composite_ppi()
that emerged with the evolution of dependency package sp / deprecation of rgdal (#666)+diff --git a/dev/reference/apply_mistnet.html b/dev/reference/apply_mistnet.html index 8c4fa983..653da534 100644 --- a/dev/reference/apply_mistnet.html +++ b/dev/reference/apply_mistnet.html @@ -18,7 +18,7 @@bioRad 0.7.32023-10-20
+
- Replaced the
+ggmap
package withggspatial
for map visualizations. This change was made asggmap
no longer provides reliable open-source basemaps without the necessity to register for an API key (#638).diff --git a/dev/reference/add_expected_eta_to_scan.html b/dev/reference/add_expected_eta_to_scan.html index 717cb0e7..2d1a6920 100644 --- a/dev/reference/add_expected_eta_to_scan.html +++ b/dev/reference/add_expected_eta_to_scan.html @@ -17,7 +17,7 @@Bugfixes
Corrected incorrect mapping of to
dbz_all
data column in VPTS CSV format to correspondingDBZH
data column in bioRad vpts object (#634).- diff --git a/dev/pkgdown.yml b/dev/pkgdown.yml index 34185e9e..c7a40f0e 100644 --- a/dev/pkgdown.yml +++ b/dev/pkgdown.yml @@ -5,7 +5,7 @@ articles: bioRad: bioRad.html rad_aero_22: rad_aero_22.html range_correction: range_correction.html -last_built: 2024-07-20T17:37Z +last_built: 2024-08-27T16:18Z urls: reference: https://adriaandokter.com/bioRad/reference article: https://adriaandokter.com/bioRad/articles diff --git a/dev/reference/Ops.html b/dev/reference/Ops.html index c9050249..d35d8e6e 100644 --- a/dev/reference/Ops.html +++ b/dev/reference/Ops.html @@ -17,7 +17,7 @@
Improved the polar volume downloads by
download_pvolfiles()
: now skips a day if there are issues with fetching instead of aborting the entire download (#636)See also
Examples
diff --git a/dev/reference/as.data.frame.vp.html b/dev/reference/as.data.frame.vp.html index 95fa5996..a8950b02 100644 --- a/dev/reference/as.data.frame.vp.html +++ b/dev/reference/as.data.frame.vp.html @@ -20,7 +20,7 @@# \donttest{ # make sure you have installed the MistNet libraries and model, using: +if (requireNamespace("vol2birdR", quietly = TRUE)){ if(!vol2birdR::mistnet_exists()){ vol2birdR::install_mistnet() vol2birdR::install_mistnet_model() @@ -209,52 +210,50 @@
Examples
# Calculate MistNet segmentation mistnet_pvol <- apply_mistnet(tempfile) -#> Filename = /tmp/RtmpcO7UT4/KBGM_example1a1d708d9899, callid = KBGM -#> Reading RSL polar volume with nominal time 20191001-000542, source: RAD:KBGM,PLC:BINGHAMTON,state:NY,radar_name:KBGM -#> Running vol2birdSetUp -#> Warning: radial velocities will be dealiased... -#> Warning: disabling single-polarization precipitation filter for S-band data, continuing in DUAL polarization mode -#> Warning: using MistNet, disabling other segmentation methods -#> Running segmentScansUsingMistnet. -#> Warning: Requested elevation scan at 1.500000 degrees but selected scan at 1.318360 degrees -#> Warning: Requested elevation scan at 3.500000 degrees but selected scan at 3.120115 degrees -#> Warning: Requested elevation scan at 4.500000 degrees but selected scan at 3.999025 degrees -#> Warning: Ignoring scan(s) not used as MistNet input: 2 4 8 9 ... -#> Running MistNet...done # Print summary info for the segmented elevation scan at the 0.5 degree, # verify new parameters BIOLOGY, WEATHER, BACKGROUND and CELL have been added scan <- get_scan(mistnet_pvol, 0.5) scan -#> Polar scan (class scan) -#> -#> parameters: DBZH VRADH RHOHV WRADH WEATHER BACKGROUND CELL PHIDP BIOLOGY ZDR -#> elevation angle: 0.483395 deg -#> dims: 1201 bins x 720 rays # Project the scan as a ppi ppi <- project_as_ppi(scan, range_max = 100000) # Plot the reflectivity parameter plot(ppi, param = "DBZH") - # Plot the MistNet class probability [0-1] for weather plot(ppi, param = "WEATHER") - # Plot the MistNet class probability [0-1] for biology plot(ppi, param = "BIOLOGY") - # Plot the final segmentation result, with values >1 indicating # areas classified as weather, and value 1 pixels that fall within an # additional 5 km fringe around weather areas plot(ppi, param = "CELL") - # Remove file file.remove(tempfile) +} +#> Starting installation of MistNet... +#> Checking write permissions... +#> Installing MistNet libraries... +#> Initializing MistNet... +#> MistNet initialized successfully. +#> MistNet installation complete. +#> Filename = /tmp/RtmpkybTUi/KBGM_example5ec2616d29a8, callid = KBGM +#> Reading RSL polar volume with nominal time 20191001-000542, source: RAD:KBGM,PLC:BINGHAMTON,state:NY,radar_name:KBGM +#> Running vol2birdSetUp +#> Warning: radial velocities will be dealiased... +#> Warning: disabling single-polarization precipitation filter for S-band data, continuing in DUAL polarization mode +#> Warning: using MistNet, disabling other segmentation methods +#> Running segmentScansUsingMistnet. +#> Warning: Requested elevation scan at 1.500000 degrees but selected scan at 1.318360 degrees +#> Warning: Requested elevation scan at 3.500000 degrees but selected scan at 3.120115 degrees +#> Warning: Requested elevation scan at 4.500000 degrees but selected scan at 3.999025 degrees +#> Warning: Ignoring scan(s) not used as MistNet input: 2 4 8 9 ... +#> Running MistNet...done #> [1] TRUE # }See also
Other beam_functions:
beam_height()
, -beam_profile_overlap()
,beam_profile()
, +beam_profile_overlap()
,beam_range()
,beam_width()
,gaussian_beam_profile()
See also
Other beam_functions:
beam_distance()
, -beam_profile_overlap()
,beam_profile()
, +beam_profile_overlap()
,beam_range()
,beam_width()
,gaussian_beam_profile()
Other beam_functions:
beam_distance()
,
beam_height()
,
-beam_profile_overlap()
,
beam_profile()
,
+beam_profile_overlap()
,
beam_width()
,
gaussian_beam_profile()
Other beam_functions:
beam_distance()
,
beam_height()
,
-beam_profile_overlap()
,
beam_profile()
,
+beam_profile_overlap()
,
beam_range()
,
gaussian_beam_profile()
# Locate and read the polar volume example file
pvolfile_source <- system.file("extdata", "volume.h5", package = "bioRad")
-# Copy the file to a home directory with read/write permissions
+# Copy the file to a temporary directory with read/write permissions
pvolfile <- paste0(tempdir(),"/volume.h5")
file.copy(pvolfile_source, pvolfile)
#> [1] TRUE
# Calculate the profile
+if (requireNamespace("vol2birdR", quietly = TRUE)) {
vp <- calculate_vp(pvolfile)
-#> Running vol2birdSetUp
-#> Warning: radial velocities will be dealiased...
# Get summary info
vp
-#> Vertical profile (class vp)
-#>
-#> radar: seang
-#> source: WMO:02606,RAD:SE50,PLC:Angelholm,NOD:seang,ORG:82,CTY:643,CMT:Swedish radar
-#> nominal time: 2015-10-18 18:00:00
-#> generated by: vol2bird 0.6.0.9201
# Clean up
file.remove(pvolfile)
+
+}
+#> Running vol2birdSetUp
+#> Warning: radial velocities will be dealiased...
#> [1] TRUE
+
# create temporary directory
+if (requireNamespace("aws.s3", quietly = TRUE)) {
temp_dir <- paste0(tempdir(),"/bioRad_tmp_files")
dir.create(temp_dir)
download_pvolfiles(
@@ -137,12 +138,13 @@ Examples
directory = temp_dir,
overwrite = TRUE
)
+# Clean up
+unlink(temp_dir, recursive = TRUE)
+}
#> Downloading data from noaa-nexrad-level2 for radar KBRO spanning over 1 days
#>
#> Downloading pvol for 2002/10/01/KBRO/
#>
|
| | 0%
|
|======================================================================| 100%
-# Clean up
-unlink(temp_dir, recursive = TRUE)
Other beam_functions:
beam_distance()
,
beam_height()
,
-beam_profile_overlap()
,
beam_profile()
,
+beam_profile_overlap()
,
beam_range()
,
beam_width()
Tidyverse methods for bioRad objects
list_vpts_aloft(radars = "bejab")
-#> Warning: Radar data found between 2018-10 and 2024-07 but not every date has radar data
+ # \donttest{
+if (requireNamespace("aws.s3", quietly = TRUE)) {
+list_vpts_aloft(radars = "bejab", date_min='2018-10-01', date_max = '2018-12-31')
+}# }
#> https://aloftdata.s3-eu-west-1.amazonaws.com/baltrad/monthly/bejab/2018/bejab_vpts_201810.csv.gz
#> https://aloftdata.s3-eu-west-1.amazonaws.com/baltrad/monthly/bejab/2018/bejab_vpts_201811.csv.gz
#> https://aloftdata.s3-eu-west-1.amazonaws.com/baltrad/monthly/bejab/2018/bejab_vpts_201812.csv.gz
-#> https://aloftdata.s3-eu-west-1.amazonaws.com/baltrad/monthly/bejab/2019/bejab_vpts_201901.csv.gz
-#> https://aloftdata.s3-eu-west-1.amazonaws.com/baltrad/monthly/bejab/2023/bejab_vpts_202302.csv.gz
-#> https://aloftdata.s3-eu-west-1.amazonaws.com/baltrad/monthly/bejab/2023/bejab_vpts_202303.csv.gz
-#> https://aloftdata.s3-eu-west-1.amazonaws.com/baltrad/monthly/bejab/2023/bejab_vpts_202304.csv.gz
-#> https://aloftdata.s3-eu-west-1.amazonaws.com/baltrad/monthly/bejab/2023/bejab_vpts_202305.csv.gz
-#> https://aloftdata.s3-eu-west-1.amazonaws.com/baltrad/monthly/bejab/2023/bejab_vpts_202306.csv.gz
-#> https://aloftdata.s3-eu-west-1.amazonaws.com/baltrad/monthly/bejab/2023/bejab_vpts_202307.csv.gz
-#> https://aloftdata.s3-eu-west-1.amazonaws.com/baltrad/monthly/bejab/2023/bejab_vpts_202308.csv.gz
-#> https://aloftdata.s3-eu-west-1.amazonaws.com/baltrad/monthly/bejab/2023/bejab_vpts_202309.csv.gz
-#> https://aloftdata.s3-eu-west-1.amazonaws.com/baltrad/monthly/bejab/2023/bejab_vpts_202310.csv.gz
-#> https://aloftdata.s3-eu-west-1.amazonaws.com/baltrad/monthly/bejab/2023/bejab_vpts_202311.csv.gz
-#> https://aloftdata.s3-eu-west-1.amazonaws.com/baltrad/monthly/bejab/2023/bejab_vpts_202312.csv.gz
-#> https://aloftdata.s3-eu-west-1.amazonaws.com/baltrad/monthly/bejab/2024/bejab_vpts_202401.csv.gz
-#> https://aloftdata.s3-eu-west-1.amazonaws.com/baltrad/monthly/bejab/2024/bejab_vpts_202402.csv.gz
-#> https://aloftdata.s3-eu-west-1.amazonaws.com/baltrad/monthly/bejab/2024/bejab_vpts_202403.csv.gz
-#> https://aloftdata.s3-eu-west-1.amazonaws.com/baltrad/monthly/bejab/2024/bejab_vpts_202404.csv.gz
-#> https://aloftdata.s3-eu-west-1.amazonaws.com/baltrad/monthly/bejab/2024/bejab_vpts_202405.csv.gz
-#> https://aloftdata.s3-eu-west-1.amazonaws.com/baltrad/monthly/bejab/2024/bejab_vpts_202406.csv.gz
-#> https://aloftdata.s3-eu-west-1.amazonaws.com/baltrad/monthly/bejab/2024/bejab_vpts_202407.csv.gz
# Project a scan as a ppi
ppi <- project_as_ppi(example_scan)
# \donttest{
+if (all(sapply(c("ggspatial","prettymapr", "rosm"), requireNamespace, quietly = TRUE))) {
# Choose a basemap
basemap <- rosm::osm.types()[1]
# Map the radial velocity of the ppi onto the basemap
map(ppi, map = basemap, param = "VRADH")
-#> Zoom: 8
-#> Fetching 4 missing tiles
-#>
|
| | 0%
|
|================== | 25%
|
|=================================== | 50%
|
|==================================================== | 75%
|
|======================================================================| 100%
-#> ...complete!
-
# Extend the plotting range of velocities, from -50 to 50 m/s
map(ppi, map = basemap, param = "VRADH", zlim = c(-50, 50))
-#> Zoom: 8
-
# Map the reflectivity
map(ppi, map = basemap, param = "DBZH")
-#> Zoom: 8
-
# Change the color palette to Viridis colors
map(ppi, map = basemap, param = "DBZH", palette = viridis::viridis(100), zlim=c(-10,10))
-#> Zoom: 8
-
# Give the data more transparency
map(ppi, map = basemap, param = "DBZH", alpha = 0.3)
-#> Zoom: 8
-
# Change the appearance of the symbol indicating the radar location
map(ppi, map = basemap, radar_size = 5, radar_color = "blue")
-#> Zoom: 8
-
# Crop the map
map(ppi, map = basemap, xlim = c(12.4, 13.2), ylim = c(56, 56.5))
+}
#> Zoom: 8
-
+#> Fetching 4 missing tiles
+#>
|
| | 0%
|
|================== | 25%
|
|=================================== | 50%
|
|==================================================== | 75%
|
|======================================================================| 100%
+#> ...complete!
+
# }
Some functions require package aws.s3 +This helper function allows to skip a test if aws.s3 is not available +Inspired by https://testthat.r-lib.org/articles/skipping.html#helpers.
+skip_if_no_aws.s3()
Invisibly returns TRUE if aws.s3 is available, otherwise skips the test with a message "Package aws.s3 not installed".
+Function map depends on several spatial dependencies (ggspatial, prettymapr, rosm). +This helper function allows to skip a test if these dependencies are not available +Inspired by https://testthat.r-lib.org/articles/skipping.html#helpers.
+skip_if_no_mapping()
Invisibly returns TRUE if dependencies available, otherwise skips the test with a message "map() dependencies (ggspatial, prettymapr, rosm) not installed".
+dplyr select method require package tidyselect +This helper function allows to skip a test if tidyselect is not available +Inspired by https://testthat.r-lib.org/articles/skipping.html#helpers.
+skip_if_no_tidyselect()
Invisibly returns TRUE if tidyselect is available, otherwise skips the test with a message "Package tidyselect not installed".
+# S3 method for class 'scan'
-select(.data, ...)
select.scan(.data, ...)
A dataset containing the schema definitions for VPTS data validation.
+vpts_schema
An object of class list
of length 3.
A list object containing the data model for the standard VPTS exchange format.
+# Reload vpts_schema from package (e.g. in case it was altered)
+data(vpts_schema)
+