weatherOz for BOM
Jonathan Carroll, Dean Marchiori, Paul Melloy, Mark Padgham, Hugh Parsonage, Keith Pembleton and Adam H. Sparks
Source:vignettes/weatherOz_for_BOM.Rmd
weatherOz_for_BOM.RmdUsing {weatherOz} With BOM Data
Several functions are provided by {weatherOz} to retrieve Australian Bureau of Meteorology (BOM) data. A family of functions retrieve data files from BOM or parse local files from BOM and return data frames;
-
get_precis_forecast()and `parse_precis_forecast(), which return the précis (short) forecast; -
get_ag_bulletin()and `parse_ag_bulletin(), which return the agriculture bulletin; -
get_coastal_forecast()and `parse_coastal_forecast(), which return coastal waters forecasts for each state.
Using get_precis_forecast()
This function only takes one argument, state. The
state parameter allows the user to select the forecast for
just one state or a national forecast. States or territories are
specified using the official postal codes or full name with fuzzy
matching performed via agrep().
ACT - Australian Capital Territory
NSW - New South Wales
NT - Northern Territory
QLD - Queensland
SA - South Australia
TAS - Tasmania
VIC - Victoria
WA - Western Australia
AUS - Australia, returns national forecast including all states, NT and ACT.
Results
The function, get_precis_forecast(), will return a data
frame of the weather forecast for the daily forecast for selected towns.
See Appendix 1 for a full description of the fields and values.
Example 1: Getting a Forecast
Following is an example fetching the forecast for Queensland.
library(weatherOz)
(QLD_forecast <- get_precis_forecast(state = "QLD"))
#> ----------- Australian Bureau of Meteorology (BOM) Précis Forecast. -
#> The HTML version of Short Form (Précis) Forecast for
#> QLD can be found at:
#> <http://www.bom.gov.au/qld/forecasts/state.shtml>
#> Please note information at the page
#> <http://www.bom.gov.au/catalogue/data-feeds.shtml#precis>.
#> --------------------------------------------------------------------------------
#> index product_id state town aac lat lon elev start_time_local end_time_local
#> <fctr> <char> <char> <char> <char> <num> <num> <num> <POSc> <POSc>
#> 1: 0 IDQ11295 QLD Brisbane QLD_PT001 -27.48080 153.0389 8.1 2024-04-19 08:38:18 2024-04-20
#> 2: 1 IDQ11295 QLD Brisbane QLD_PT001 -27.48080 153.0389 8.1 2024-04-20 00:00:00 2024-04-21
#> 3: 2 IDQ11295 QLD Brisbane QLD_PT001 -27.48080 153.0389 8.1 2024-04-21 00:00:00 2024-04-22
#> 4: 3 IDQ11295 QLD Brisbane QLD_PT001 -27.48080 153.0389 8.1 2024-04-22 00:00:00 2024-04-23
#> 5: 4 IDQ11295 QLD Brisbane QLD_PT001 -27.48080 153.0389 8.1 2024-04-23 00:00:00 2024-04-24
#> ---
#> 787: 2 IDQ11295 QLD Port Douglas QLD_PT254 -16.48681 145.4635 70.4 2024-04-21 00:00:00 2024-04-22
#> 788: 3 IDQ11295 QLD Port Douglas QLD_PT254 -16.48681 145.4635 70.4 2024-04-22 00:00:00 2024-04-23
#> 789: 4 IDQ11295 QLD Port Douglas QLD_PT254 -16.48681 145.4635 70.4 2024-04-23 00:00:00 2024-04-24
#> 790: 5 IDQ11295 QLD Port Douglas QLD_PT254 -16.48681 145.4635 70.4 2024-04-24 00:00:00 2024-04-25
#> 791: 6 IDQ11295 QLD Port Douglas QLD_PT254 -16.48681 145.4635 70.4 2024-04-25 00:00:00 2024-04-26
#> utc_offset start_time_utc end_time_utc minimum_temperature maximum_temperature
#> <fctr> <POSc> <POSc> <num> <num>
#> 1: 10:00 2024-04-18 22:38:18 2024-04-19 14:00:00 NA 28
#> 2: 10:00 2024-04-19 14:00:00 2024-04-20 14:00:00 17 25
#> 3: 10:00 2024-04-20 14:00:00 2024-04-21 14:00:00 17 24
#> 4: 10:00 2024-04-21 14:00:00 2024-04-22 14:00:00 19 25
#> 5: 10:00 2024-04-22 14:00:00 2024-04-23 14:00:00 19 26
#> ---
#> 787: 10:00 2024-04-20 14:00:00 2024-04-21 14:00:00 23 29
#> 788: 10:00 2024-04-21 14:00:00 2024-04-22 14:00:00 21 28
#> 789: 10:00 2024-04-22 14:00:00 2024-04-23 14:00:00 21 29
#> 790: 10:00 2024-04-23 14:00:00 2024-04-24 14:00:00 22 29
#> 791: 10:00 2024-04-24 14:00:00 2024-04-25 14:00:00 22 29
#> lower_precipitation_limit upper_precipitation_limit precis probability_of_precipitation
#> <num> <num> <char> <num>
#> 1: NA NA Becoming partly cloudy. 20
#> 2: 0 10 Rain. 80
#> 3: 2 25 Rain. 95
#> 4: 0 6 Shower or two. 70
#> 5: 0 2 Shower or two. 60
#> ---
#> 787: 0 1 Partly cloudy. 30
#> 788: NA NA Mostly sunny. 20
#> 789: 0 1 Partly cloudy. 30
#> 790: 0 1 Possible shower. 40
#> 791: 0 1 Partly cloudy. 30Using get_ag_bulletin()
get_ag_bulletin() only takes one argument,
state. The state parameter allows the user to
select the bulletin for just one state or a national forecast. States or
territories are specified using the official postal codes or full name
with fuzzy matching performed via agrep().
NSW - New South Wales
NT - Northern Territory
QLD - Queensland
SA - South Australia
TAS - Tasmania
VIC - Victoria
WA - Western Australia
AUS - Australia, returns bulletin for all states and NT.
Results
The function, get_ag_bulletin(), will return a data
frame of the agriculture bulletin for selected stations. See Appendix 3
for a full list and description of the fields and values.
Example 2: Getting an Ag Bulletin
Following is an example fetching the ag bulletin for Western Australia.
library(weatherOz)
(WA_bulletin <- get_ag_bulletin(state = "WA"))
#> ----------- Australian Bureau of Meteorology (BOM) Ag Bulletin. -----------
#> Please note information at the foot of
#> <http://www.bom.gov.au/cgi-bin/wrap_fwo.pl?IDW65176.html>,
#> the HTML version of Agricultural Observations Bulletin for
#> WA. Also see
#> <http://www.bom.gov.au/catalogue/observations/about-agricultural.shtml>.
#> ----------------------------------------------------------------------------
#> product_id state station site obs_time_local obs_time_utc time_zone r tn tx
#> <fctr> <char> <char> <fctr> <POSc> <POSc> <char> <num> <num> <num>
#> 1: IDW65176 WA Beverley 10515 2024-04-19 09:00:00 2024-04-19 01:00:00 WST 0 10.8 32.0
#> 2: IDW65176 WA Brookton 10524 2024-04-19 09:00:00 2024-04-19 01:00:00 WST 0 10.4 30.2
#> 3: IDW65176 WA Broome 03003 2024-04-19 09:00:00 2024-04-19 01:00:00 WST 0 23.8 36.0
#> 4: IDW65176 WA Carnarvon 06011 2024-04-19 09:00:00 2024-04-19 01:00:00 WST 0 17.5 36.1
#> 5: IDW65176 WA Donnybrook 09534 2024-04-19 09:00:00 2024-04-19 01:00:00 WST 0 12.5 32.8
#> ---
#> 29: IDW65176 WA Port Hedland 04032 2024-04-19 09:00:00 2024-04-19 01:00:00 WST 0 21.5 36.8
#> 30: IDW65176 WA Ravensthorpe 10633 2024-04-19 09:00:00 2024-04-19 01:00:00 WST 0 12.5 27.0
#> 31: IDW65176 WA Salmon Gums RS 12071 2024-04-19 09:00:00 2024-04-19 01:00:00 WST 0 9.4 25.7
#> 32: IDW65176 WA Wandering 10917 2024-04-19 09:00:00 2024-04-19 01:00:00 WST 0 9.3 29.9
#> 33: IDW65176 WA York 10311 2024-04-19 09:00:00 2024-04-19 01:00:00 WST 0 8.1 32.5
#> twd ev tg sn solr t5 t10 t20 t50 t1m wr
#> <num> <num> <num> <num> <num> <num> <num> <num> <num> <num> <num>
#> 1: 5.2 NA NA NA 16.9 NA NA NA NA NA NA
#> 2: 5.2 NA NA NA 16.4 NA NA NA NA NA NA
#> 3: 12.6 NA NA NA 22.2 NA NA NA NA NA NA
#> 4: 13.1 NA NA NA 19.6 NA NA NA NA NA NA
#> 5: 5.5 NA NA NA 15.9 NA NA NA NA NA NA
#> ---
#> 29: 14.0 NA NA NA 20.7 NA NA NA NA NA NA
#> 30: 4.0 NA NA NA 16.2 NA NA NA NA NA NA
#> 31: NA NA NA NA 16.5 NA NA NA NA NA NA
#> 32: NA NA NA NA 16.4 NA NA NA NA NA NA
#> 33: 5.8 NA NA NA 16.8 NA NA NA NA NA NAUsing get_coastal_forecast()
This function only takes one argument, state. The
state parameter allows the user to select the forecast for
just one state or a national forecast. States or territories are
specified using the official postal codes or full name with fuzzy
matching performed via agrep()
ACT - Australian Capital Territory
NSW - New South Wales
NT - Northern Territory
QLD - Queensland
SA - South Australia
TAS - Tasmania
VIC - Victoria
WA - Western Australia
AUS - Australia, returns national forecast including all states, NT and ACT.
Results
The function, get_coastal_forecast(), will return a data
frame of the coastal waters forecast for marine zones in each state. See
Appendix 6 for a full description of the fields and values.
Example 3: Getting a Coastal Forecast
Following is an example fetching the forecast for Queensland.
library(weatherOz)
(QLD_coastal_forecast <- get_coastal_forecast(state = "QLD"))
#> ------- Australian Bureau of Meteorology (BOM) Coastal Waters Forecast. -------
#> Please note information at the foot of
#> <http://www.bom.gov.au/cgi-bin/wrap_fwo.pl?IDQ11290.html>,
#> the HTML version of Coastal Waters Forecast for
#> QLD.
#> Also see
#> <http://www.bom.gov.au/catalogue/observations/about-coastal-observations.shtml>.
#> --------------------------------------------------------------------------------
#> index product_id type state_code dist_name pt_1_name
#> <fctr> <char> <char> <char> <char> <char>
#> 1: <NA> IDQ11290 <NA> QLD Queensland <NA>
#> 2: <NA> IDQ11290 <NA> QLD Queensland Gulf of Carpentaria <NA>
#> 3: <NA> IDQ11290 <NA> QLD North Queensland <NA>
#> 4: <NA> IDQ11290 <NA> QLD South Queensland <NA>
#> 5: 0 IDQ11290 <NA> QLD South East Gulf of Carpentaria: QLD-NT Border to Cape Keerweer <NA>
#> ---
#> 45: 1 IDQ11290 <NA> QLD Gold Coast Waters: Cape Moreton to Point Danger <NA>
#> 46: 2 IDQ11290 <NA> QLD Gold Coast Waters: Cape Moreton to Point Danger <NA>
#> 47: 0 IDQ11290 <NA> QLD Great Barrier Reef Offshore Waters <NA>
#> 48: 1 IDQ11290 <NA> QLD Great Barrier Reef Offshore Waters <NA>
#> 49: 2 IDQ11290 <NA> QLD Great Barrier Reef Offshore Waters <NA>
#> pt_2_name aac start_time_local end_time_local utc_offset start_time_utc end_time_utc
#> <char> <char> <POSc> <POSc> <fctr> <POSc> <POSc>
#> 1: <NA> QLD_FA001 2024-04-19 11:21:06 2024-04-19 11:21:06 10:00 2024-04-19 11:21:06 2024-04-19 11:21:06
#> 2: <NA> QLD_FA002 2024-04-19 05:00:00 2024-04-19 05:00:00 10:00 2024-04-19 05:00:00 2024-04-19 05:00:00
#> 3: <NA> QLD_FA003 2024-04-19 05:00:00 2024-04-19 05:00:00 10:00 2024-04-19 05:00:00 2024-04-19 05:00:00
#> 4: <NA> QLD_FA004 2024-04-19 05:00:00 2024-04-19 05:00:00 10:00 2024-04-19 05:00:00 2024-04-19 05:00:00
#> 5: <NA> QLD_MW001 2024-04-19 11:02:29 2024-04-19 11:02:29 10:00 2024-04-19 11:02:29 2024-04-19 11:02:29
#> ---
#> 45: <NA> QLD_MW014 2024-04-20 00:00:00 2024-04-20 00:00:00 10:00 2024-04-20 00:00:00 2024-04-20 00:00:00
#> 46: <NA> QLD_MW014 2024-04-21 00:00:00 2024-04-21 00:00:00 10:00 2024-04-21 00:00:00 2024-04-21 00:00:00
#> 47: <NA> QLD_MW015 2024-04-19 05:00:00 2024-04-19 05:00:00 10:00 2024-04-19 05:00:00 2024-04-19 05:00:00
#> 48: <NA> QLD_MW015 2024-04-20 00:00:00 2024-04-20 00:00:00 10:00 2024-04-20 00:00:00 2024-04-20 00:00:00
#> 49: <NA> QLD_MW015 2024-04-21 00:00:00 2024-04-21 00:00:00 10:00 2024-04-21 00:00:00 2024-04-21 00:00:00
#> forecast_seas
#> <char>
#> 1: <NA>
#> 2: <NA>
#> 3: <NA>
#> 4: <NA>
#> 5: Around 1 metre, increasing to 1 to 1.5 metres offshore.
#> ---
#> 45: Below 1 metre, increasing to 1 to 1.5 metres during the morning.
#> 46: 1 to 2 metres, increasing to 2 to 3 metres during the afternoon.
#> 47: Below 1 metre, increasing to 1 to 1.5 metres north of Creal Reef.
#> 48: 1 to 1.5 metres.
#> 49: 1.5 to 2.5 metres.
#> forecast_weather
#> <char>
#> 1: <NA>
#> 2: <NA>
#> 3: <NA>
#> 4: <NA>
#> 5: Partly cloudy. 50% chance of showers north of Kowanyama, near zero chance elsewhere. The chance of a thunderstorm.
#> ---
#> 45: Cloudy. 95% chance of rain. The chance of a thunderstorm.
#> 46: Cloudy. 90% chance of rain. The chance of a thunderstorm.
#> 47: Partly cloudy.
#> 48: Partly cloudy.
#> 49: Partly cloudy. 70% chance of showers.
#> forecast_winds
#> <char>
#> 1: <NA>
#> 2: <NA>
#> 3: <NA>
#> 4: <NA>
#> 5: Southeasterly 15 to 20 knots.
#> ---
#> 45: Southeasterly 15 to 20 knots, reaching up to 25 knots inshore in the late afternoon.
#> 46: Southeasterly 15 to 25 knots increasing to 25 to 30 knots during the day.
#> 47: Southeasterly 15 to 20 knots.
#> 48: Southeasterly 15 to 20 knots.
#> 49: Southeasterly 20 to 30 knots.
#> forecast_swell1
#> <char>
#> 1: <NA>
#> 2: <NA>
#> 3: <NA>
#> 4: <NA>
#> 5: <NA>
#> ---
#> 45: Northeast to southeasterly around 1 metre.
#> 46: Southeasterly 1 to 1.5 metres, increasing to 1.5 to 2.5 metres during the afternoon.
#> 47: Easterly 1 to 1.5 metres, decreasing to around 1 metre by early evening.
#> 48: Easterly around 1 metre.
#> 49: East to southeasterly around 1 metre, increasing to 1 to 2 metres through the Capricorn Channel during the morning.
#> forecast_swell2 forecast_caution marine_forecast tropical_system_location forecast_waves
#> <char> <char> <char> <lgcl> <lgcl>
#> 1: <NA> <NA> <NA> NA NA
#> 2: <NA> <NA> <NA> NA NA
#> 3: <NA> <NA> <NA> NA NA
#> 4: <NA> <NA> <NA> NA NA
#> 5: <NA> <NA> <NA> NA NA
#> ---
#> 45: <NA> <NA> <NA> NA NA
#> 46: <NA> <NA> <NA> NA NA
#> 47: <NA> <NA> <NA> NA NA
#> 48: <NA> <NA> <NA> NA NA
#> 49: <NA> <NA> <NA> NA NAWorking with BOM Image Files
A second family of functions retrieve information pertaining to satellite and radar imagery.
-
get_available_imagery(), which returns available satellite imagery; -
get_satelllite_imagery()which returns a [terra::SpatRaster] or [stars] object as requested; -
get_available_radar(), which returns available radar images; and -
get_radar_imagery()which returns radar images as a [magick] object.
Using {weatherOz} to retrieve BOM satellite imagery
{weatherOz} provides functionality to retrieve high-definition GeoTIFF satellite imagery provided by BOM through public FTP with the following types of imagery being available: i.) Infrared images, ii.) Visible images and iii.) Clouds/surface composite.
Working with Satellite Imagery
Valid BOM satellite Product IDs for GeoTIFF files include:
| Product ID | Description | Type | Delete time |
|---|---|---|---|
| IDE00420 | AHI cloud cover only 2km FD GEOS | Satellite | 24 |
| IDE00421 | AHI IR (Ch13) greyscale 2km FD GEOS | Satellite | 24 |
| IDE00422 | AHI VIS (Ch3) greyscale 2km FD GEOS | Satellite | 24 |
| IDE00423 | AHI IR (Ch13) Zehr 2km FD GEOS | Satellite | 24 |
| IDE00425 | AHI VIS (true colour) / IR (Ch13 greyscale) composite 1km FD GEOS | Satellite | 24 |
| IDE00426 | AHI VIS (true colour) / IR (Ch13 greyscale) composite 2km FD GEOS | Satellite | 24 |
| IDE00427 | AHI WV (Ch8) 2km FD GEOS | Satellite | 24 |
| IDE00430 | AHI cloud cover only 2km AUS equirect. | Satellite | 24 |
| IDE00431 | AHI IR (Ch13) greyscale 2km AUS equirect. | Satellite | 24 |
| IDE00432 | AHI VIS (Ch3) greyscale 2km AUS equirect. | Satellite | 24 |
| IDE00433 | AHI IR (Ch13) Zehr 2km AUS equirect. | Satellite | 24 |
| IDE00435 | AHI VIS (true colour) / IR (Ch13 greyscale) composite 1km AUS equirect. | Satellite | 24 |
| IDE00436 | AHI VIS (true colour) / IR (Ch13 greyscale) composite 2km AUS equirect. | Satellite | 24 |
| IDE00437 | AHI WV (Ch8) 2km AUS equirect. | Satellite | 24 |
| IDE00439 | AHI VIS (Ch3) greyscale 0.5km AUS equirect. | Satellite | 24 |
| Information gathered from Australian Bureau of Meteorology (BOM) | |||
Using get_available_imagery()
get_available_imagery() only takes one argument,
product_id, a BOM identifier for the imagery that you wish
to check for available imagery. Using this function will fetch a listing
of BOM GeoTIFF satellite imagery from ftp://ftp.bom.gov.au/anon/gen/gms/
to display which files are currently available for download. These files
are available at ten minute update frequency with a 24 hour delete time.
This function can be used see the most recent files available and then
specify in the _imagery() function. If no valid Product ID
is supplied, defaults to all GeoTIFF images currently available.
Example 5: Checking Available Imagery
library(weatherOz)
(avail <- get_available_imagery(product_id = "IDE00425"))
#> [1] "IDE00425.202404180410.tif" "IDE00425.202404180420.tif" "IDE00425.202404180430.tif" "IDE00425.202404180440.tif"
#> [5] "IDE00425.202404180450.tif" "IDE00425.202404180500.tif" "IDE00425.202404180510.tif" "IDE00425.202404180520.tif"
#> [9] "IDE00425.202404180530.tif" "IDE00425.202404180540.tif" "IDE00425.202404180550.tif" "IDE00425.202404180600.tif"
#> [13] "IDE00425.202404180610.tif" "IDE00425.202404180620.tif" "IDE00425.202404180630.tif" "IDE00425.202404180640.tif"
#> [17] "IDE00425.202404180650.tif" "IDE00425.202404180700.tif" "IDE00425.202404180710.tif" "IDE00425.202404180720.tif"
#> [21] "IDE00425.202404180730.tif" "IDE00425.202404180740.tif" "IDE00425.202404180750.tif" "IDE00425.202404180800.tif"
#> [25] "IDE00425.202404180810.tif" "IDE00425.202404180820.tif" "IDE00425.202404180830.tif" "IDE00425.202404180840.tif"
#> [29] "IDE00425.202404180850.tif" "IDE00425.202404180900.tif" "IDE00425.202404180910.tif" "IDE00425.202404180920.tif"
#> [33] "IDE00425.202404180930.tif" "IDE00425.202404180940.tif" "IDE00425.202404180950.tif" "IDE00425.202404181000.tif"
#> [37] "IDE00425.202404181010.tif" "IDE00425.202404181020.tif" "IDE00425.202404181030.tif" "IDE00425.202404181040.tif"
#> [41] "IDE00425.202404181050.tif" "IDE00425.202404181100.tif" "IDE00425.202404181110.tif" "IDE00425.202404181120.tif"
#> [45] "IDE00425.202404181130.tif" "IDE00425.202404181140.tif" "IDE00425.202404181150.tif" "IDE00425.202404181200.tif"
#> [49] "IDE00425.202404181210.tif" "IDE00425.202404181220.tif" "IDE00425.202404181230.tif" "IDE00425.202404181240.tif"
#> [53] "IDE00425.202404181250.tif" "IDE00425.202404181300.tif" "IDE00425.202404181310.tif" "IDE00425.202404181320.tif"
#> [57] "IDE00425.202404181330.tif" "IDE00425.202404181340.tif" "IDE00425.202404181350.tif" "IDE00425.202404181400.tif"
#> [61] "IDE00425.202404181410.tif" "IDE00425.202404181420.tif" "IDE00425.202404181430.tif" "IDE00425.202404181450.tif"
#> [65] "IDE00425.202404181500.tif" "IDE00425.202404181510.tif" "IDE00425.202404181520.tif" "IDE00425.202404181530.tif"
#> [69] "IDE00425.202404181540.tif" "IDE00425.202404181550.tif" "IDE00425.202404181600.tif" "IDE00425.202404181610.tif"
#> [73] "IDE00425.202404181620.tif" "IDE00425.202404181630.tif" "IDE00425.202404181640.tif" "IDE00425.202404181650.tif"
#> [77] "IDE00425.202404181700.tif" "IDE00425.202404181710.tif" "IDE00425.202404181720.tif" "IDE00425.202404181730.tif"
#> [81] "IDE00425.202404181740.tif" "IDE00425.202404181750.tif" "IDE00425.202404181800.tif" "IDE00425.202404181810.tif"
#> [85] "IDE00425.202404181820.tif" "IDE00425.202404181830.tif" "IDE00425.202404181840.tif" "IDE00425.202404181850.tif"
#> [89] "IDE00425.202404181900.tif" "IDE00425.202404181910.tif" "IDE00425.202404181920.tif" "IDE00425.202404181930.tif"
#> [93] "IDE00425.202404181940.tif" "IDE00425.202404181950.tif" "IDE00425.202404182000.tif" "IDE00425.202404182010.tif"
#> [97] "IDE00425.202404182020.tif" "IDE00425.202404182030.tif" "IDE00425.202404182040.tif" "IDE00425.202404182050.tif"
#> [101] "IDE00425.202404182100.tif" "IDE00425.202404182110.tif" "IDE00425.202404182120.tif" "IDE00425.202404182130.tif"
#> [105] "IDE00425.202404182140.tif" "IDE00425.202404182150.tif" "IDE00425.202404182200.tif" "IDE00425.202404182210.tif"
#> [109] "IDE00425.202404182220.tif" "IDE00425.202404182230.tif" "IDE00425.202404182240.tif" "IDE00425.202404182250.tif"
#> [113] "IDE00425.202404182300.tif" "IDE00425.202404182310.tif" "IDE00425.202404182320.tif" "IDE00425.202404182330.tif"
#> [117] "IDE00425.202404182340.tif" "IDE00425.202404182350.tif" "IDE00425.202404190000.tif" "IDE00425.202404190010.tif"
#> [121] "IDE00425.202404190020.tif" "IDE00425.202404190030.tif" "IDE00425.202404190040.tif" "IDE00425.202404190050.tif"
#> [125] "IDE00425.202404190100.tif" "IDE00425.202404190110.tif" "IDE00425.202404190120.tif" "IDE00425.202404190130.tif"
#> [129] "IDE00425.202404190140.tif" "IDE00425.202404190150.tif" "IDE00425.202404190200.tif" "IDE00425.202404190210.tif"
#> [133] "IDE00425.202404190220.tif" "IDE00425.202404190230.tif" "IDE00425.202404190250.tif" "IDE00425.202404190300.tif"
#> [137] "IDE00425.202404190310.tif" "IDE00425.202404190320.tif" "IDE00425.202404190330.tif" "IDE00425.202404190340.tif"
#> [141] "IDE00425.202404190350.tif" "IDE00425.202404190400.tif" "IDE00425.202404190410.tif" "IDE00425.202404190420.tif"
#> [145] "IDE00425.202404190430.tif" "IDE00425.202404190440.tif" "IDE00425.202404190450.tif"Using get_satellite_imagery()
get_satellite_imagery() fetches BOM satellite GeoTIFF
imagery, returning a SpatRaster object and takes two arguments. Files
are available at ten minute update frequency with a 24 hour delete time.
It is suggested to check file availability first by using
get_available_imagery(). The arguments are:
product_id, a character value of the BOM product ID to download. Alternatively, a vector of values fromget_available_imagery()may be used here. This argument is mandatory.scansa numeric value for the number of scans to download, starting with the most recent and progressing backwards, e.g.,1- the most recent single scan available ,6- the most recent hour available,12- the most recent 2 hours available, etc. Negating will return the oldest files first. Defaults to 1. This argument is optional.
Example 6: Fetching Satellite Imagery and Viewing It
library(weatherOz)
# Specify product ID and scans
i <- get_satellite_imagery(product_id = "IDE00425", scans = 1)terra::plot() has been re-exported to simplify
visualising these files while using {weatherOz}.
plot(i)
plot of chunk plot_satellite
Using {weatherOz} With BOM Radar Imagery
{weatherOz} provides functionality to retrieve the latest radar imagery provided by BOM through a public FTP server. These are the latest snapshots for each radar locations at various radar ranges e.g., 512km, 256km, 128km and 64km for some stations.
Using get_available_radar()
get_available_radar() fetches the available radar
imagery from the BOM FTP and returns a data frame for reference. This
data frame contains the product_id, which is required when using the
get_radar_imagery() function. The files available are the
latest .png files of BOM radar imagery which are typically
updated each 6-10 minutes. Only the most recent image is retrieved for
each radar location. There are usually several radar ranges available
for each radar location, such as 512km, 256km, 128km and possibly 64km.
The arguments are:
-
radar_idwhich is the BOM radar ID number; this defaults to ‘all’ which will return a data frame of all radar IDs in Australia.
Example 7: Getting Available Radar Imagery
library(weatherOz)
x <- get_available_radar()
head(x)
#> product_id LocationID range Name Longitude Latitude Radar_id Full_Name IDRnn0name IDRnn1name
#> <char> <char> <char> <fctr> <num> <num> <int> <fctr> <fctr> <fctr>
#> 1: IDR641 64 512km Adelaide 138.4689 -34.6169 64 Adelaide (Buckland Park) BuckPk BucklandPk
#> 2: IDR642 64 256km Adelaide 138.4689 -34.6169 64 Adelaide (Buckland Park) BuckPk BucklandPk
#> 3: IDR643 64 128km Adelaide 138.4689 -34.6169 64 Adelaide (Buckland Park) BuckPk BucklandPk
#> 4: IDR644 64 64km Adelaide 138.4689 -34.6169 64 Adelaide (Buckland Park) BuckPk BucklandPk
#> 5: IDR311 31 512km Albany 117.8163 -34.9418 31 Albany Albany Albany
#> 6: IDR312 31 256km Albany 117.8163 -34.9418 31 Albany Albany Albany
#> State Type Group_ Status Archive as.is
#> <fctr> <fctr> <fctr> <fctr> <fctr> <lgcl>
#> 1: SA Doppler Yes Public BuckPk TRUE
#> 2: SA Doppler Yes Public BuckPk TRUE
#> 3: SA Doppler Yes Public BuckPk TRUE
#> 4: SA Doppler Yes Public BuckPk TRUE
#> 5: WA Doppler Yes Public Albany TRUE
#> 6: WA Doppler Yes Public Albany TRUEUsing get_radar_imagery()
get_radar_imagery() fetches the latest BOM radar imagery
for a given product ID. The files available are the latest
.png files of BOM radar imagery, which are typically
updated each 6-10 minutes. Only the most recent image is retrieved for
each radar location. There are usually several radar ranges available
for each radar location, such as 512km, 256km, 128km and possibly 64km.
The only argument is:
-
product_idthe BOM product_id associated with each radar imagery file. These can be obtained from theget_available_radar()function. This value must be specified and the function will accept only one at a time.
Example 8: Fetching Radar Imagery
library(weatherOz)
y <- get_radar_imagery(product_id = "IDR032")
#> file downloaded to:/var/folders/r4/wwsd3hsn48j5gck6qv6npkpc0000gr/T//Rtmp4Yk6lY/file6ad77634cf2e.png
plot(y)
plot of chunk get_radar_imagery
References
Australian Bureau of Meteorology (BOM) Weather Data Services
Australian Bureau of Meteorology (BOM) FTP Public Products
Australian Bureau of Meteorology (BOM) Weather Data Services Agriculture Bulletins
Australian Bureau of Meteorology (BOM) Weather Data Services Observation of Rainfall
Australian Bureau of Meteorology (BOM) High-definition satellite images
Appendix 1 - Output From get_precis_forecast()
The functions, get_precis_forecast() or
`parse_precis_forecast(), will return a data frame of the 7 day short
forecast with the following fields:
- index
- Forecast index number, 0 = current day … 7 day
- product_id
- BOM Product ID from which the data are derived
- state
- State name (postal code abbreviation)
- town
- Town name for forecast location
- aac
- AMOC Area Code, e.g., WA_MW008, a unique identifier for each location
- lat
- Latitude of named location (decimal degrees)
- lon
- Longitude of named location (decimal degrees)
- elev
- Elevation of named location (metres)
- start_time_local
- Start of forecast date and time in local TZ
- end_time_local
- End of forecast date and time in local TZ
- UTC_offset
-
Hours offset from difference in hours and minutes from Coordinated
Universal Time (UTC) for
start_time_localandend_time_local - start_time_utc
- Start of forecast date and time in UTC
- end_time_utc
- End of forecast date and time in UTC
- maximum_temperature
- Maximum forecast temperature (degrees Celsius)
- minimum_temperature
- Minimum forecast temperature (degrees Celsius)
- lower_precipitation_limit
- Lower forecast precipitation limit (millimetres)
- upper_precipitation_limit
- Upper forecast precipitation limit (millimetres)
- precis
- Précis forecast (a short summary, less than 30 characters)
- probability_of_precipitation
- Probability of precipitation (percent)
Appendix 2 - Output From get_ag_bulletin()
The functions, get_ag_bulletin() or
`parse_ag_bulletin(), will return a data frame of the agriculture
bulletin with the following fields:
- product_id
- BOM Product ID from which the data are derived
- state
- State name (postal code abbreviation)
- dist
- BOM rainfall district
- name
- Full station name (some stations have been retired so “station” will be same, this is the full designation
- wmo
- World Meteorological Organization number (unique ID used worldwide)
- site
- Unique BOM identifier for each station
- station
- Station name
- obs-time-local
- Observation time
- obs-time-utc
- Observation time (time in UTC)
- time-zone
- Time zone for observation
- lat
- Latitude (decimal degrees)
- lon
- Longitude (decimal degrees)
- elev_m
- Station elevation (metres)
- bar_ht
- Bar height (metres)
- station
- BOM station name
- start
- Year data collection starts
- end
- Year data collection ends (will always be current)
- r
- Rain to 9am (millimetres). Trace will be reported as 0.01
- tn
- Minimum temperature (degrees Celsius)
- tx
- Maximum temperature (degrees Celsius)
- twd
- Wet bulb depression (degrees Celsius)
- ev
- Evaporation (millimetres)
- tg
- Terrestrial minimum temperature (degrees Celsius)
- sn
- Sunshine (hours)
- solr
- Solar Radiation MJ/sq m
- t5
- 5cm soil temperature (degrees Celsius)
- t10
- 10cm soil temperature (degrees Celsius)
- t20
- 20cm soil temperature (degrees Celsius)
- t50
- 50cm soil temperature (degrees Celsius)
- t1m
- 1m soil temperature (degrees Celsius)
- wr
- Wind run (kilometres)
Appendix 3 - Output From get_coastal_forecast()
The output of get_coastal_forecast() or
parse_coastal_forecast() will return a data frame with
coastal waters forecast values of each area within the given state with
the following fields:
- index
- Forecast index number. 0 = current day
- product_id
- BOM Product ID from which the data are derived
- type
- Forecast Region type e.g. Coastal
- state_code
- State name (postal code abbreviation)
- dist_name
- Name of forecast district
- pt_1_name
- Start of forecast district
- pt_2_name
- End of forecast district
- aac
- AMOC Area Code, e.g., WA_MW008, a unique identifier for each location
- start_time_local
- Start of forecast date and time in local TZ
- end_time_local
- End of forecast date and time in local TZ
- UTC_offset
-
Hours offset from difference in hours and minutes from Coordinated
Universal Time (UTC) for
start_time_localandend_time_local - start_time_utc
- Start of forecast date and time in UTC
- end_time_utc
- End of forecast date and time in UTC
- forecast_seas
- Forecast sea conditions
- forecast_weather
- Forecast weather summary
- forecast_winds
- Forecast winds summary
- forecast_swell1
- Forecast primary swell summary
- forecast_swell2
- Forecast seondary swell summary (not always provided)
- forecast_caution
- Forecast caution issued (not always provided)
- marine_forecast
- Additional marine forecast warning information (not always provided)