diff --git a/vignettes/buoy_vignette.Rmd b/vignettes/buoy_vignette.Rmd index 0786d1a4..d4e5da1d 100644 --- a/vignettes/buoy_vignette.Rmd +++ b/vignettes/buoy_vignette.Rmd @@ -40,13 +40,13 @@ head(res) #> 4 41006 #> 5 41008 #> 6 41009 -#> url -#> 1 http://dods.ndbc.noaa.gov/thredds/catalog/data/cwind/41001/catalog.html -#> 2 http://dods.ndbc.noaa.gov/thredds/catalog/data/cwind/41002/catalog.html -#> 3 http://dods.ndbc.noaa.gov/thredds/catalog/data/cwind/41004/catalog.html -#> 4 http://dods.ndbc.noaa.gov/thredds/catalog/data/cwind/41006/catalog.html -#> 5 http://dods.ndbc.noaa.gov/thredds/catalog/data/cwind/41008/catalog.html -#> 6 http://dods.ndbc.noaa.gov/thredds/catalog/data/cwind/41009/catalog.html +#> url +#> 1 https://dods.ndbc.noaa.gov/thredds/catalog/data/cwind/41001/catalog.html +#> 2 https://dods.ndbc.noaa.gov/thredds/catalog/data/cwind/41002/catalog.html +#> 3 https://dods.ndbc.noaa.gov/thredds/catalog/data/cwind/41004/catalog.html +#> 4 https://dods.ndbc.noaa.gov/thredds/catalog/data/cwind/41006/catalog.html +#> 5 https://dods.ndbc.noaa.gov/thredds/catalog/data/cwind/41008/catalog.html +#> 6 https://dods.ndbc.noaa.gov/thredds/catalog/data/cwind/41009/catalog.html ``` Or browse them on the web @@ -116,24 +116,24 @@ Including year and datatype ```r -buoy(dataset = 'cwind', buoyid = 41001, year = 2008, datatype = "cc") +buoy(dataset = 'cwind', buoyid = 45005, year = 2008, datatype = "c") ``` ``` -#> Dimensions (rows/cols): [1585 X 5] +#> Dimensions (rows/cols): [29688 X 5] #> 2 variables: [wind_dir, wind_spd] #> #> time lat lon wind_dir wind_spd -#> 1 2008-05-28T16:00:00Z 34.704 -72.734 230 8.6 -#> 2 2008-05-28T16:10:00Z 34.704 -72.734 230 8.7 -#> 3 2008-05-28T16:20:00Z 34.704 -72.734 229 8.5 -#> 4 2008-05-28T16:30:00Z 34.704 -72.734 231 8.8 -#> 5 2008-05-28T16:40:00Z 34.704 -72.734 236 8.5 -#> 6 2008-05-28T16:50:00Z 34.704 -72.734 235 8.9 -#> 7 2008-05-28T17:00:00Z 34.704 -72.734 233 8.2 -#> 8 2008-05-28T17:10:00Z 34.704 -72.734 233 8.2 -#> 9 2008-05-28T17:20:00Z 34.704 -72.734 231 8.3 -#> 10 2008-05-28T17:30:00Z 34.704 -72.734 232 7.8 +#> 1 2008-04-29T09:00:00Z 41.677 -82.398 10 9.0 +#> 2 2008-04-29T09:10:00Z 41.677 -82.398 8 9.0 +#> 3 2008-04-29T09:20:00Z 41.677 -82.398 5 9.3 +#> 4 2008-04-29T09:30:00Z 41.677 -82.398 13 9.5 +#> 5 2008-04-29T09:40:00Z 41.677 -82.398 14 9.4 +#> 6 2008-04-29T09:50:00Z 41.677 -82.398 12 9.4 +#> 7 2008-04-29T14:00:00Z 41.677 -82.398 341 6.5 +#> 8 2008-04-29T14:10:00Z 41.677 -82.398 332 6.8 +#> 9 2008-04-29T14:20:00Z 41.677 -82.398 335 6.4 +#> 10 2008-04-29T14:30:00Z 41.677 -82.398 332 6.5 #> .. ... ... ... ... ... ``` @@ -141,23 +141,23 @@ Including just datatype ```r -buoy(dataset = 'cwind', buoyid = 41001, datatype = "cc") +buoy(dataset = 'cwind', buoyid = 45005, datatype = "c") ``` ``` -#> Dimensions (rows/cols): [1585 X 5] +#> Dimensions (rows/cols): [26784 X 5] #> 2 variables: [wind_dir, wind_spd] #> -#> time lat lon wind_dir wind_spd -#> 1 2008-05-28T16:00:00Z 34.704 -72.734 230 8.6 -#> 2 2008-05-28T16:10:00Z 34.704 -72.734 230 8.7 -#> 3 2008-05-28T16:20:00Z 34.704 -72.734 229 8.5 -#> 4 2008-05-28T16:30:00Z 34.704 -72.734 231 8.8 -#> 5 2008-05-28T16:40:00Z 34.704 -72.734 236 8.5 -#> 6 2008-05-28T16:50:00Z 34.704 -72.734 235 8.9 -#> 7 2008-05-28T17:00:00Z 34.704 -72.734 233 8.2 -#> 8 2008-05-28T17:10:00Z 34.704 -72.734 233 8.2 -#> 9 2008-05-28T17:20:00Z 34.704 -72.734 231 8.3 -#> 10 2008-05-28T17:30:00Z 34.704 -72.734 232 7.8 -#> .. ... ... ... ... ... +#> time lat lon wind_dir wind_spd +#> 1 1996-05-15T23:00:00Z 41.68 -82.4 337 2.2 +#> 2 1996-05-15T23:10:00Z 41.68 -82.4 282 1.0 +#> 3 1996-05-15T23:20:00Z 41.68 -82.4 282 2.2 +#> 4 1996-05-15T23:30:00Z 41.68 -82.4 258 2.6 +#> 5 1996-05-15T23:40:00Z 41.68 -82.4 254 3.0 +#> 6 1996-05-15T23:50:00Z 41.68 -82.4 252 2.7 +#> 7 1996-05-16T00:00:00Z 41.68 -82.4 240 2.1 +#> 8 1996-05-16T00:10:00Z 41.68 -82.4 246 2.4 +#> 9 1996-05-16T00:20:00Z 41.68 -82.4 251 2.7 +#> 10 1996-05-16T00:30:00Z 41.68 -82.4 253 2.9 +#> .. ... ... ... ... ... ``` diff --git a/vignettes/figure/seaice1-1.png b/vignettes/figure/seaice1-1.png index b49317bd..eb5e3f21 100644 Binary files a/vignettes/figure/seaice1-1.png and b/vignettes/figure/seaice1-1.png differ diff --git a/vignettes/figure/seaice2-1.png b/vignettes/figure/seaice2-1.png index 5b8d68c6..16057301 100644 Binary files a/vignettes/figure/seaice2-1.png and b/vignettes/figure/seaice2-1.png differ diff --git a/vignettes/figure/six-1.png b/vignettes/figure/six-1.png index 062d2287..5aae1769 100644 Binary files a/vignettes/figure/six-1.png and b/vignettes/figure/six-1.png differ diff --git a/vignettes/figure/unnamed-chunk-10-1.png b/vignettes/figure/unnamed-chunk-10-1.png index 28542681..3e26e3c6 100644 Binary files a/vignettes/figure/unnamed-chunk-10-1.png and b/vignettes/figure/unnamed-chunk-10-1.png differ diff --git a/vignettes/figure/unnamed-chunk-11-1.png b/vignettes/figure/unnamed-chunk-11-1.png index acc8232f..893799e1 100644 Binary files a/vignettes/figure/unnamed-chunk-11-1.png and b/vignettes/figure/unnamed-chunk-11-1.png differ diff --git a/vignettes/figure/unnamed-chunk-3-1.png b/vignettes/figure/unnamed-chunk-3-1.png index 21b18db9..7a5a7703 100644 Binary files a/vignettes/figure/unnamed-chunk-3-1.png and b/vignettes/figure/unnamed-chunk-3-1.png differ diff --git a/vignettes/figure/unnamed-chunk-4-1.png b/vignettes/figure/unnamed-chunk-4-1.png index 166c4796..af04bcf9 100644 Binary files a/vignettes/figure/unnamed-chunk-4-1.png and b/vignettes/figure/unnamed-chunk-4-1.png differ diff --git a/vignettes/figure/unnamed-chunk-5-1.png b/vignettes/figure/unnamed-chunk-5-1.png index 387da4e6..aabb4dae 100644 Binary files a/vignettes/figure/unnamed-chunk-5-1.png and b/vignettes/figure/unnamed-chunk-5-1.png differ diff --git a/vignettes/figure/unnamed-chunk-7-1.png b/vignettes/figure/unnamed-chunk-7-1.png index 387da4e6..aabb4dae 100644 Binary files a/vignettes/figure/unnamed-chunk-7-1.png and b/vignettes/figure/unnamed-chunk-7-1.png differ diff --git a/vignettes/figure/unnamed-chunk-8-1.png b/vignettes/figure/unnamed-chunk-8-1.png index 197d1fc8..565dd82a 100644 Binary files a/vignettes/figure/unnamed-chunk-8-1.png and b/vignettes/figure/unnamed-chunk-8-1.png differ diff --git a/vignettes/homr_vignette.Rmd b/vignettes/homr_vignette.Rmd index edb0c3f6..3b1a809c 100644 --- a/vignettes/homr_vignette.Rmd +++ b/vignettes/homr_vignette.Rmd @@ -183,19 +183,24 @@ names(res) ``` ``` -#> [1] "10001871" "10100161" "10100162" "10100164" "10100166" "20004155" -#> [7] "20004158" "20004160" "20004162" "20004163" "20004167" "20004168" -#> [13] "20004171" "20004176" "20004178" "20004179" "20004180" "20004182" -#> [19] "20004184" "20004185" "30001464" "30001561" "30001831" "30017384" -#> [25] "30020917" "30021161" "30021998" "30022674" "30026770" "30027455" -#> [31] "30032423" "30032685" "30034222" "30039554" "30043742" "30046662" -#> [37] "30046814" "30051475" "30057217" "30063570" "30064900" "30065901" -#> [43] "30067636" "30069663" "30075067" "30077378" "30077857" "30077923" -#> [49] "30077988" "30079088" "30079240" "30082430" "30084216" "30084262" -#> [55] "30084537" "30084796" "30094582" "30094639" "30094664" "30094670" -#> [61] "30094683" "30094730" "30094806" "30094830" "30094917" "30094931" -#> [67] "30094936" "30101120" "30101186" "30101294" "30101341" "30101363" -#> [73] "30101370" "30102005" +#> [1] "10001871" "10100161" "10100162" "10100164" "10100166" "20004155" +#> [7] "20004158" "20004160" "20004162" "20004163" "20004167" "20004168" +#> [13] "20004171" "20004176" "20004178" "20004179" "20004180" "20004182" +#> [19] "20004184" "20004185" "30001464" "30001561" "30001831" "30017384" +#> [25] "30020917" "30021161" "30021998" "30022674" "30026770" "30027455" +#> [31] "30032423" "30032685" "30034222" "30039554" "30043742" "30046662" +#> [37] "30046814" "30051475" "30057217" "30063570" "30064900" "30065901" +#> [43] "30067636" "30069663" "30075067" "30077378" "30077857" "30077923" +#> [49] "30077988" "30079088" "30079240" "30082430" "30084216" "30084262" +#> [55] "30084537" "30084796" "30094582" "30094639" "30094664" "30094670" +#> [61] "30094683" "30094730" "30094806" "30094830" "30094917" "30094931" +#> [67] "30094936" "30101120" "30101186" "30101294" "30101341" "30101363" +#> [73] "30101370" "30102005" "30106954" "30107175" "30107251" "30107345" +#> [79] "30107377" "30107384" "30107435" "30107641" "30107692" "30107701" +#> [85] "30107718" "30107886" "30107891" "30107899" "30107936" "30108008" +#> [91] "30108119" "30108282" "30108333" "30108338" "30108394" "30108839" +#> [97] "30108887" "30108929" "30108935" "30108946" "30109183" "30110071" +#> [103] "30111663" "30111793" "30112317" "30112507" "30112571" "30113052" ``` By country @@ -297,20 +302,20 @@ head( ldply(res, "[[", "head") ) ``` ``` -#> .id preferredName latitude_dec longitude_dec -#> 1 30077883 WEAVERVILLE 4.2 N, NC 35.7579 -82.5618 -#> 2 30029796 FAIRVIEW 1.2 S, NC 35.5058 -82.4051 -#> 3 30083542 BLACK MOUNTAIN 0.8 N, NC 35.6263 -82.3297 -#> 4 20013839 LEICESTER, NC 35.65 -82.7 -#> 5 20013817 SWANNANOA 2 E, NC 35.6 -82.36667 -#> 6 20014046 SWANNANOA 2 SSE, NC 35.57333 -82.385 +#> .id preferredName latitude_dec longitude_dec +#> 1 30063513 ASHEVILLE 5.6 NNW, NC 35.6534 -82.5709 +#> 2 30014782 WEAVERVILLE 5E, NC 35.71028 -82.50389 +#> 3 20013851 WEAVERVILLE, NC 35.7 -82.56667 +#> 4 20013789 GARREN CREEK, NC 35.51667 -82.33333 +#> 5 20013801 BILTMORE 2 SE, NC 35.56833 -82.545 +#> 6 30031683 ASHEVILLE 4.0 NNE, NC 35.6224 -82.5096 #> por.beginDate por.endDate precision #> 1 Unknown Present -#> 2 Unknown Present -#> 3 Unknown Present -#> 4 1934-11-05T00:00:00.000 1962-03-31T00:00:00.000 DDMM -#> 5 1931-03-01T00:00:00.000 1987-08-26T00:00:00.000 DDMM -#> 6 1984-01-01T00:00:00.000 2008-03-31T00:00:00.000 DDMMSS +#> 2 1965-12-01T00:00:00.000 1976-01-31T00:00:00.000 DDMMSS +#> 3 1946-03-30T00:00:00.000 1992-10-01T00:00:00.000 DDMM +#> 4 1936-09-25T00:00:00.000 1962-03-31T00:00:00.000 DDMM +#> 5 1963-10-01T00:00:00.000 2007-11-14T00:00:00.000 DDMMSS +#> 6 Unknown Present ``` ## Get header information only @@ -322,20 +327,20 @@ head( ldply(res, "[[", "head") ) ``` ``` -#> .id preferredName latitude_dec longitude_dec precision -#> 1 20004159 GEORGETOWN 5 SW, DE 38.63333 -75.45 DDMM -#> 2 30067636 LEWES 0.8 SE, DE 38.7731 -75.1385 -#> 3 30027455 DELMAR 2.8 NE, DE 38.4863 -75.5332 -#> 4 30075067 PRIME HOOK DELAWARE, DE 38.8333 -75.3333 -#> 5 20004158 REHOBOTH BEACH LBS, DE 38.61667 -75.06667 DDMM -#> 6 30039554 GREENWOOD 2.9 SE, DE 38.7731 -75.5616 -#> por.beginDate por.endDate -#> 1 1946-09-14T00:00:00.000 1997-06-01T00:00:00.000 -#> 2 Unknown Present -#> 3 Unknown Present -#> 4 Unknown Present -#> 5 Unknown Present -#> 6 Unknown Present +#> .id preferredName latitude_dec longitude_dec +#> 1 30094670 HARRINGTON 3.3 SSW, DE 38.8821 -75.6052 +#> 2 10100166 WILMINGTON PORTER RES, DE 39.7739 -75.5414 +#> 3 30001464 BEAR 2 SW, DE 39.5917 -75.7325 +#> 4 30101363 LEWES 1.5 SSW, DE 38.7588 -75.1579 +#> 5 20004172 ADAMSVILLE, DE 38.81667 -75.7 +#> 6 20004165 BRIDGEVILLE STATE POLICE, DE 38.77444 -75.6075 +#> por.beginDate por.endDate precision +#> 1 Unknown Present +#> 2 1912-07-12T00:00:00.000 Present DDddddd +#> 3 2003-02-01T00:00:00.000 2013-04-02T00:00:00.000 DDMMSS +#> 4 Unknown Present +#> 5 1967-09-01T00:00:00.000 1969-12-31T00:00:00.000 DDMM +#> 6 1975-08-01T00:00:00.000 2003-09-01T00:00:00.000 DDMMSS ``` ## Data definitions diff --git a/vignettes/ncdc_vignette.Rmd b/vignettes/ncdc_vignette.Rmd index b93dfbf5..da0c362f 100644 --- a/vignettes/ncdc_vignette.Rmd +++ b/vignettes/ncdc_vignette.Rmd @@ -31,7 +31,7 @@ ncdc_stations(datasetid='GHCND', locationid='FIPS:12017', stationid='GHCND:USC00 #> #> $data #> elevation mindate maxdate latitude name -#> 1 12.2 1899-02-01 2015-06-16 28.8029 INVERNESS 3 SE, FL US +#> 1 12.2 1899-02-01 2018-11-25 28.8029 INVERNESS 3 SE, FL US #> datacoverage id elevationUnit longitude #> 1 1 GHCND:USC00084289 METERS -82.3126 #> @@ -45,32 +45,20 @@ ncdc_stations(datasetid='GHCND', locationid='FIPS:12017', stationid='GHCND:USC00 ```r out <- ncdc(datasetid='NORMAL_DLY', datatypeid='dly-tmax-normal', startdate = '2010-05-01', enddate = '2010-05-10') out$data -#> date datatype station value fl_c -#> 1 2010-05-01T00:00:00 DLY-TMAX-NORMAL GHCND:AQW00061705 869 C -#> 2 2010-05-01T00:00:00 DLY-TMAX-NORMAL GHCND:CAW00064757 607 Q -#> 3 2010-05-01T00:00:00 DLY-TMAX-NORMAL GHCND:CQC00914080 840 R -#> 4 2010-05-01T00:00:00 DLY-TMAX-NORMAL GHCND:CQC00914801 858 R -#> 5 2010-05-01T00:00:00 DLY-TMAX-NORMAL GHCND:FMC00914395 876 P -#> 6 2010-05-01T00:00:00 DLY-TMAX-NORMAL GHCND:FMC00914419 885 P -#> 7 2010-05-01T00:00:00 DLY-TMAX-NORMAL GHCND:FMC00914446 885 P -#> 8 2010-05-01T00:00:00 DLY-TMAX-NORMAL GHCND:FMC00914482 868 R -#> 9 2010-05-01T00:00:00 DLY-TMAX-NORMAL GHCND:FMC00914720 899 R -#> 10 2010-05-01T00:00:00 DLY-TMAX-NORMAL GHCND:FMC00914761 897 P -#> 11 2010-05-01T00:00:00 DLY-TMAX-NORMAL GHCND:FMC00914831 870 P -#> 12 2010-05-01T00:00:00 DLY-TMAX-NORMAL GHCND:FMC00914892 883 P -#> 13 2010-05-01T00:00:00 DLY-TMAX-NORMAL GHCND:FMC00914898 875 P -#> 14 2010-05-01T00:00:00 DLY-TMAX-NORMAL GHCND:FMC00914911 885 P -#> 15 2010-05-01T00:00:00 DLY-TMAX-NORMAL GHCND:FMW00040308 888 S -#> 16 2010-05-01T00:00:00 DLY-TMAX-NORMAL GHCND:FMW00040504 879 C -#> 17 2010-05-01T00:00:00 DLY-TMAX-NORMAL GHCND:FMW00040505 867 S -#> 18 2010-05-01T00:00:00 DLY-TMAX-NORMAL GHCND:GQC00914025 852 P -#> 19 2010-05-01T00:00:00 DLY-TMAX-NORMAL GHCND:GQW00041415 877 C -#> 20 2010-05-01T00:00:00 DLY-TMAX-NORMAL GHCND:JQW00021603 852 P -#> 21 2010-05-01T00:00:00 DLY-TMAX-NORMAL GHCND:PSC00914519 883 P -#> 22 2010-05-01T00:00:00 DLY-TMAX-NORMAL GHCND:PSC00914712 840 P -#> 23 2010-05-01T00:00:00 DLY-TMAX-NORMAL GHCND:PSW00040309 879 S -#> 24 2010-05-01T00:00:00 DLY-TMAX-NORMAL GHCND:RMW00040604 867 S -#> 25 2010-05-01T00:00:00 DLY-TMAX-NORMAL GHCND:RMW00040710 863 C +#> # A tibble: 25 x 5 +#> date datatype station value fl_c +#> +#> 1 2010-05-01T00:00:00 DLY-TMAX-NORMAL GHCND:AQW00061705 869 C +#> 2 2010-05-01T00:00:00 DLY-TMAX-NORMAL GHCND:CAW00064757 607 Q +#> 3 2010-05-01T00:00:00 DLY-TMAX-NORMAL GHCND:CQC00914080 840 R +#> 4 2010-05-01T00:00:00 DLY-TMAX-NORMAL GHCND:CQC00914801 858 R +#> 5 2010-05-01T00:00:00 DLY-TMAX-NORMAL GHCND:FMC00914395 876 P +#> 6 2010-05-01T00:00:00 DLY-TMAX-NORMAL GHCND:FMC00914419 885 P +#> 7 2010-05-01T00:00:00 DLY-TMAX-NORMAL GHCND:FMC00914446 885 P +#> 8 2010-05-01T00:00:00 DLY-TMAX-NORMAL GHCND:FMC00914482 868 R +#> 9 2010-05-01T00:00:00 DLY-TMAX-NORMAL GHCND:FMC00914720 899 R +#> 10 2010-05-01T00:00:00 DLY-TMAX-NORMAL GHCND:FMC00914761 897 P +#> # … with 15 more rows ``` ## Plot data, super simple, but it's a start @@ -81,7 +69,7 @@ out <- ncdc(datasetid='NORMAL_DLY', stationid='GHCND:USW00014895', datatypeid='d ncdc_plot(out) ``` -![plot of chunk six](figure/six-1.png) +![plot of chunk six](figure/six-1.png) ## More on plotting @@ -157,34 +145,24 @@ Then combine with a call to `ncdc_combine` ```r df <- ncdc_combine(out1, out2) head(df[[1]]); tail(df[[1]]) -#> date datatype station value fl_m fl_q fl_so -#> 1 2010-03-01T00:00:00 PRCP GHCND:USW00014895 0 T 0 -#> 2 2010-03-02T00:00:00 PRCP GHCND:USW00014895 0 T 0 -#> 3 2010-03-03T00:00:00 PRCP GHCND:USW00014895 0 T 0 -#> 4 2010-03-04T00:00:00 PRCP GHCND:USW00014895 0 0 -#> 5 2010-03-05T00:00:00 PRCP GHCND:USW00014895 0 0 -#> 6 2010-03-06T00:00:00 PRCP GHCND:USW00014895 0 0 -#> fl_t -#> 1 2400 -#> 2 2400 -#> 3 2400 -#> 4 2400 -#> 5 2400 -#> 6 2400 -#> date datatype station value fl_m fl_q fl_so -#> 148 2010-10-26T00:00:00 PRCP GHCND:USW00014895 221 0 -#> 149 2010-10-27T00:00:00 PRCP GHCND:USW00014895 0 0 -#> 150 2010-10-28T00:00:00 PRCP GHCND:USW00014895 0 T 0 -#> 151 2010-10-29T00:00:00 PRCP GHCND:USW00014895 0 T 0 -#> 152 2010-10-30T00:00:00 PRCP GHCND:USW00014895 0 0 -#> 153 2010-10-31T00:00:00 PRCP GHCND:USW00014895 0 0 -#> fl_t -#> 148 2400 -#> 149 2400 -#> 150 2400 -#> 151 2400 -#> 152 2400 -#> 153 2400 +#> # A tibble: 6 x 8 +#> date datatype station value fl_m fl_q fl_so fl_t +#> +#> 1 2010-03-01T00:00:… PRCP GHCND:USW00014… 0 T "" 0 2400 +#> 2 2010-03-02T00:00:… PRCP GHCND:USW00014… 0 T "" 0 2400 +#> 3 2010-03-03T00:00:… PRCP GHCND:USW00014… 0 T "" 0 2400 +#> 4 2010-03-04T00:00:… PRCP GHCND:USW00014… 0 "" "" 0 2400 +#> 5 2010-03-05T00:00:… PRCP GHCND:USW00014… 0 "" "" 0 2400 +#> 6 2010-03-06T00:00:… PRCP GHCND:USW00014… 0 "" "" 0 2400 +#> # A tibble: 6 x 8 +#> date datatype station value fl_m fl_q fl_so fl_t +#> +#> 1 2010-10-26T00:00:… PRCP GHCND:USW00014… 221 "" "" 0 2400 +#> 2 2010-10-27T00:00:… PRCP GHCND:USW00014… 0 "" "" 0 2400 +#> 3 2010-10-28T00:00:… PRCP GHCND:USW00014… 0 T "" 0 2400 +#> 4 2010-10-29T00:00:… PRCP GHCND:USW00014… 0 T "" 0 2400 +#> 5 2010-10-30T00:00:… PRCP GHCND:USW00014… 0 "" "" 0 2400 +#> 6 2010-10-31T00:00:… PRCP GHCND:USW00014… 0 "" "" 0 2400 ``` Then plot - the default passing in the combined plot plots the data together. In this case it looks kind of weird since a straight line combines two distant dates. diff --git a/vignettes/ncdc_workflow.Rmd b/vignettes/ncdc_workflow.Rmd index 312d1179..81b83201 100644 --- a/vignettes/ncdc_workflow.Rmd +++ b/vignettes/ncdc_workflow.Rmd @@ -98,8 +98,8 @@ ncdc_datasets(stationid = id) #> #> $data #> uid mindate maxdate name -#> 1 gov.noaa.ncdc:C00861 1763-01-01 2017-05-01 Daily Summaries -#> 2 gov.noaa.ncdc:C00946 1763-01-01 2017-04-01 Global Summary of the Month +#> 1 gov.noaa.ncdc:C00861 1763-01-01 2018-11-25 Daily Summaries +#> 2 gov.noaa.ncdc:C00946 1763-01-01 2018-10-01 Global Summary of the Month #> datacoverage id #> 1 1 GHCND #> 2 1 GSOM @@ -119,8 +119,8 @@ ncdc_datatypes(datasetid = "GHCND", stationid = id) #> #> $data #> mindate maxdate name datacoverage id -#> 1 1781-01-01 2017-04-30 Precipitation 1 PRCP -#> 2 1840-05-01 2017-04-30 Snowfall 1 SNOW +#> 1 1781-01-01 2018-11-25 Precipitation 1 PRCP +#> 2 1840-05-01 2018-11-25 Snowfall 1 SNOW #> #> attr(,"class") #> [1] "ncdc_datatypes" @@ -153,52 +153,20 @@ ncdc(datasetid = "GHCND", stationid = id, datatypeid = "PRCP", startdate = "2012 #> #> #> $data -#> date datatype station value fl_m fl_q fl_so -#> 1 2012-10-01T00:00:00 PRCP GHCND:US1FLCT0002 13 N -#> 2 2012-10-02T00:00:00 PRCP GHCND:US1FLCT0002 3 N -#> 3 2012-10-03T00:00:00 PRCP GHCND:US1FLCT0002 15 N -#> 4 2012-10-04T00:00:00 PRCP GHCND:US1FLCT0002 142 N -#> 5 2012-10-05T00:00:00 PRCP GHCND:US1FLCT0002 244 N -#> 6 2012-10-06T00:00:00 PRCP GHCND:US1FLCT0002 655 N -#> 7 2012-10-07T00:00:00 PRCP GHCND:US1FLCT0002 23 N -#> 8 2012-10-08T00:00:00 PRCP GHCND:US1FLCT0002 0 N -#> 9 2012-10-09T00:00:00 PRCP GHCND:US1FLCT0002 3 N -#> 10 2012-10-10T00:00:00 PRCP GHCND:US1FLCT0002 0 N -#> 11 2012-10-11T00:00:00 PRCP GHCND:US1FLCT0002 0 N -#> 12 2012-10-12T00:00:00 PRCP GHCND:US1FLCT0002 0 N -#> 13 2012-10-13T00:00:00 PRCP GHCND:US1FLCT0002 0 N -#> 14 2012-10-14T00:00:00 PRCP GHCND:US1FLCT0002 0 N -#> 15 2012-10-15T00:00:00 PRCP GHCND:US1FLCT0002 0 N -#> 16 2012-10-16T00:00:00 PRCP GHCND:US1FLCT0002 0 N -#> 17 2012-10-17T00:00:00 PRCP GHCND:US1FLCT0002 0 N -#> 18 2012-10-18T00:00:00 PRCP GHCND:US1FLCT0002 0 T N -#> 19 2012-10-19T00:00:00 PRCP GHCND:US1FLCT0002 0 N -#> 20 2012-10-20T00:00:00 PRCP GHCND:US1FLCT0002 0 N -#> 21 2012-10-21T00:00:00 PRCP GHCND:US1FLCT0002 0 N -#> 22 2012-10-24T00:00:00 PRCP GHCND:US1FLCT0002 0 N -#> fl_t -#> 1 -#> 2 -#> 3 -#> 4 -#> 5 -#> 6 -#> 7 -#> 8 -#> 9 -#> 10 -#> 11 -#> 12 -#> 13 -#> 14 -#> 15 -#> 16 -#> 17 -#> 18 -#> 19 -#> 20 -#> 21 -#> 22 +#> # A tibble: 22 x 8 +#> date datatype station value fl_m fl_q fl_so fl_t +#> +#> 1 2012-10-01T00:00… PRCP GHCND:US1FLCT0… 13 "" "" N "" +#> 2 2012-10-02T00:00… PRCP GHCND:US1FLCT0… 3 "" "" N "" +#> 3 2012-10-03T00:00… PRCP GHCND:US1FLCT0… 15 "" "" N "" +#> 4 2012-10-04T00:00… PRCP GHCND:US1FLCT0… 142 "" "" N "" +#> 5 2012-10-05T00:00… PRCP GHCND:US1FLCT0… 244 "" "" N "" +#> 6 2012-10-06T00:00… PRCP GHCND:US1FLCT0… 655 "" "" N "" +#> 7 2012-10-07T00:00… PRCP GHCND:US1FLCT0… 23 "" "" N "" +#> 8 2012-10-08T00:00… PRCP GHCND:US1FLCT0… 0 "" "" N "" +#> 9 2012-10-09T00:00… PRCP GHCND:US1FLCT0… 3 "" "" N "" +#> 10 2012-10-10T00:00… PRCP GHCND:US1FLCT0… 0 "" "" N "" +#> # … with 12 more rows #> #> attr(,"class") #> [1] "ncdc_data" diff --git a/vignettes/rnoaa_ropenaq.Rmd b/vignettes/rnoaa_ropenaq.Rmd index cce2fd56..2b6610d8 100644 --- a/vignettes/rnoaa_ropenaq.Rmd +++ b/vignettes/rnoaa_ropenaq.Rmd @@ -4,6 +4,8 @@ %\VignetteEncoding{UTF-8} --> + + # Complementing air quality data with weather data using rnoaa __Author__ Maëlle Salmon @@ -149,6 +151,33 @@ all_monitors_clean <- meteo_pull_monitors(monitors, date_max = as.character(Sys.Date())) %>% rename(day = date, location = id) +``` + +``` +#> file path: /Users/sckott/Library/Caches/rnoaa/ghcnd/US1INMR0134.dly +``` + +``` +#> file last updated: 2018-12-03 11:12:04 +``` + +``` +#> file min/max dates: 2015-05-01 / 2018-11-30 +``` + +``` +#> file path: /Users/sckott/Library/Caches/rnoaa/ghcnd/US1INMR0122.dly +``` + +``` +#> file last updated: 2018-12-03 11:12:06 +``` + +``` +#> file min/max dates: 2012-11-01 / 2018-11-30 +``` + +```r all_monitors_clean %>% head() %>% knitr::kable() ``` @@ -156,12 +185,12 @@ all_monitors_clean %>% head() %>% knitr::kable() |location |day | prcp| snow| |:-----------|:----------|----:|----:| -|US1INMR0134 |2017-10-21 | NA| NA| -|US1INMR0134 |2017-10-22 | NA| NA| -|US1INMR0134 |2017-10-23 | NA| NA| -|US1INMR0134 |2017-10-24 | NA| NA| -|US1INMR0134 |2017-10-25 | NA| NA| -|US1INMR0134 |2017-10-26 | NA| NA| +|US1INMR0134 |2018-11-03 | NA| NA| +|US1INMR0134 |2018-11-04 | NA| NA| +|US1INMR0134 |2018-11-05 | NA| NA| +|US1INMR0134 |2018-11-06 | 76| NA| +|US1INMR0134 |2018-11-07 | NA| NA| +|US1INMR0134 |2018-11-08 | NA| NA| Here we notice some values are not available. Therefore, we might need to go back to weather stations searching with, for instance, a larger radius. In this case let's say we're ok with the result of the search. @@ -205,4 +234,4 @@ ggplot(data_plot) + facet_grid(parameter ~ ., scales = "free_y") ``` -![plot of chunk unnamed-chunk-9](figure/unnamed-chunk-9-1.png) +![plot of chunk unnamed-chunk-10](figure/unnamed-chunk-10-1.png) diff --git a/vignettes/storms_vignette.Rmd b/vignettes/storms_vignette.Rmd index bd731632..65811d0e 100644 --- a/vignettes/storms_vignette.Rmd +++ b/vignettes/storms_vignette.Rmd @@ -85,53 +85,53 @@ storm_data(basin='WP') ``` ``` -#> # A tibble: 39,584 × 195 -#> serial_num season num basin sub_basin name -#> -#> 1 1884177N17124 1884 1 WP MM NOT NAMED -#> 2 1884177N17124 1884 1 WP MM NOT NAMED -#> 3 1884177N17124 1884 1 WP MM NOT NAMED -#> 4 1884177N17124 1884 1 WP MM NOT NAMED -#> 5 1884177N17124 1884 1 WP MM NOT NAMED -#> 6 1884177N17124 1884 1 WP MM NOT NAMED -#> 7 1884177N17124 1884 1 WP MM NOT NAMED -#> 8 1884177N17124 1884 1 WP MM NOT NAMED -#> 9 1884177N17124 1884 1 WP MM NOT NAMED -#> 10 1884177N17124 1884 1 WP MM NOT NAMED -#> # ... with 39,574 more rows, and 189 more variables: iso_time , -#> # nature , latitude , longitude , wind.wmo. , -#> # pres.wmo. , center , wind.wmo..percentile , -#> # pres.wmo..percentile , track_type , -#> # latitude_for_mapping , longitude_for_mapping , -#> # current.basin , hurdat_atl_lat , hurdat_atl_lon , -#> # hurdat_atl_grade , hurdat_atl_wind , hurdat_atl_pres , -#> # td9636_lat , td9636_lon , td9636_grade , -#> # td9636_wind , td9636_pres , reunion_lat , -#> # reunion_lon , reunion_grade , reunion_wind , -#> # reunion_pres , atcf_lat , atcf_lon , atcf_grade , -#> # atcf_wind , atcf_pres , ds824_sh_lat , -#> # ds824_sh_lon , ds824_sh_grade , ds824_sh_wind , -#> # ds824_sh_pres , ds824_ni_lat , ds824_ni_lon , -#> # ds824_ni_grade , ds824_ni_wind , ds824_ni_pres , -#> # bom_lat , bom_lon , bom_grade , bom_wind , -#> # bom_pres , ds824_au_lat , ds824_au_lon , -#> # ds824_au_grade , ds824_au_wind , ds824_au_pres , -#> # jtwc_sh_lat , jtwc_sh_lon , jtwc_sh_grade , -#> # jtwc_sh_wind , jtwc_sh_pres , jtwc_wp_lat , -#> # jtwc_wp_lon , jtwc_wp_grade , jtwc_wp_wind , -#> # jtwc_wp_pres , td9635_lat , td9635_lon , -#> # td9635_grade , td9635_wind , td9635_pres , -#> # ds824_wp_lat , ds824_wp_lon , ds824_wp_grade , -#> # ds824_wp_wind , ds824_wp_pres , jtwc_io_lat , -#> # jtwc_io_lon , jtwc_io_grade , jtwc_io_wind , -#> # jtwc_io_pres , cma_lat , cma_lon , cma_grade , -#> # cma_wind , cma_pres , hurdat_epa_lat , +#> # A tibble: 99,028 x 200 +#> serial_num season num basin sub_basin name iso_time nature latitude +#> +#> 1 1884177N1… 1884 1 " WP" " MM" NOT … 1884-06… " NR" -999 +#> 2 1884177N1… 1884 1 " WP" " MM" NOT … 1884-06… " NR" -999 +#> 3 1884177N1… 1884 1 " WP" " MM" NOT … 1884-06… " NR" -999 +#> 4 1884177N1… 1884 1 " WP" " MM" NOT … 1884-06… " NR" -999 +#> 5 1884177N1… 1884 1 " WP" " MM" NOT … 1884-06… " NR" -999 +#> 6 1884177N1… 1884 1 " WP" " MM" NOT … 1884-06… " NR" -999 +#> 7 1884177N1… 1884 1 " WP" " MM" NOT … 1884-06… " NR" -999 +#> 8 1884177N1… 1884 1 " WP" " MM" NOT … 1884-06… " NR" -999 +#> 9 1884177N1… 1884 1 " WP" " MM" NOT … 1884-06… " NR" -999 +#> 10 1884177N1… 1884 1 " WP" " MM" NOT … 1884-06… " NR" -999 +#> # … with 99,018 more rows, and 191 more variables: longitude , +#> # wind.wmo. , pres.wmo. , center , +#> # wind.wmo..percentile , pres.wmo..percentile , +#> # track_type , latitude_for_mapping , +#> # longitude_for_mapping , current.basin , +#> # hurdat_atl_lat , hurdat_atl_lon , hurdat_atl_grade , +#> # hurdat_atl_wind , hurdat_atl_pres , td9636_lat , +#> # td9636_lon , td9636_grade , td9636_wind , +#> # td9636_pres , reunion_lat , reunion_lon , +#> # reunion_grade , reunion_wind , reunion_pres , +#> # atcf_lat , atcf_lon , atcf_grade , atcf_wind , +#> # atcf_pres , mlc_natl_lat , mlc_natl_lon , +#> # mlc_natl_grade , mlc_natl_wind , mlc_natl_pres , +#> # ds824_sh_lat , ds824_sh_lon , ds824_sh_grade , +#> # ds824_sh_wind , ds824_sh_pres , ds824_ni_lat , +#> # ds824_ni_lon , ds824_ni_grade , ds824_ni_wind , +#> # ds824_ni_pres , bom_lat , bom_lon , bom_grade , +#> # bom_wind , bom_pres , ds824_au_lat , +#> # ds824_au_lon , ds824_au_grade , ds824_au_wind , +#> # ds824_au_pres , jtwc_sh_lat , jtwc_sh_lon , +#> # jtwc_sh_grade , jtwc_sh_wind , jtwc_sh_pres , +#> # jtwc_wp_lat , jtwc_wp_lon , jtwc_wp_grade , +#> # jtwc_wp_wind , jtwc_wp_pres , td9635_lat , +#> # td9635_lon , td9635_grade , td9635_wind , +#> # td9635_pres , ds824_wp_lat , ds824_wp_lon , +#> # ds824_wp_grade , ds824_wp_wind , ds824_wp_pres , +#> # jtwc_io_lat , jtwc_io_lon , jtwc_io_grade , +#> # jtwc_io_wind , jtwc_io_pres , cma_lat , cma_lon , +#> # cma_grade , cma_wind , cma_pres , hurdat_epa_lat , #> # hurdat_epa_lon , hurdat_epa_grade , hurdat_epa_wind , #> # hurdat_epa_pres , jtwc_ep_lat , jtwc_ep_lon , #> # jtwc_ep_grade , jtwc_ep_wind , jtwc_ep_pres , #> # ds824_ep_lat , ds824_ep_lon , ds824_ep_grade , -#> # ds824_ep_wind , ds824_ep_pres , jtwc_cp_lat , -#> # jtwc_cp_lon , ... +#> # ds824_ep_wind , ds824_ep_pres , … ``` Or by storm id, see `storm_meta('storm_names')` for storm ids and their names. @@ -142,26 +142,27 @@ storm_data(storm='1970143N19091') ``` ``` -#> # A tibble: 5 × 195 -#> serial_num season num basin sub_basin name iso_time -#> -#> 1 1970143N19091 1970 2 NI BB 02B 1970-05-22 18:00:00 -#> 2 1970143N19091 1970 2 NI BB 02B 1970-05-23 00:00:00 -#> 3 1970143N19091 1970 2 NI BB 02B 1970-05-23 06:00:00 -#> 4 1970143N19091 1970 2 NI BB 02B 1970-05-23 12:00:00 -#> 5 1970143N19091 1970 2 NI BB 02B 1970-05-23 18:00:00 -#> # ... with 188 more variables: nature , latitude , -#> # longitude , wind.wmo. , pres.wmo. , center , -#> # wind.wmo..percentile , pres.wmo..percentile , -#> # track_type , latitude_for_mapping , -#> # longitude_for_mapping , current.basin , -#> # hurdat_atl_lat , hurdat_atl_lon , hurdat_atl_grade , -#> # hurdat_atl_wind , hurdat_atl_pres , td9636_lat , -#> # td9636_lon , td9636_grade , td9636_wind , -#> # td9636_pres , reunion_lat , reunion_lon , -#> # reunion_grade , reunion_wind , reunion_pres , -#> # atcf_lat , atcf_lon , atcf_grade , atcf_wind , -#> # atcf_pres , ds824_sh_lat , ds824_sh_lon , +#> # A tibble: 5 x 200 +#> serial_num season num basin sub_basin name iso_time nature latitude +#> +#> 1 1970143N1… 1970 2 " NI" " BB" 02B 1970-05… " TS" -999 +#> 2 1970143N1… 1970 2 " NI" " BB" 02B 1970-05… " TS" -999 +#> 3 1970143N1… 1970 2 " NI" " BB" 02B 1970-05… " TS" -999 +#> 4 1970143N1… 1970 2 " NI" " BB" 02B 1970-05… " TS" -999 +#> 5 1970143N1… 1970 2 " NI" " BB" 02B 1970-05… " NR" -999 +#> # … with 191 more variables: longitude , wind.wmo. , +#> # pres.wmo. , center , wind.wmo..percentile , +#> # pres.wmo..percentile , track_type , +#> # latitude_for_mapping , longitude_for_mapping , +#> # current.basin , hurdat_atl_lat , hurdat_atl_lon , +#> # hurdat_atl_grade , hurdat_atl_wind , hurdat_atl_pres , +#> # td9636_lat , td9636_lon , td9636_grade , +#> # td9636_wind , td9636_pres , reunion_lat , +#> # reunion_lon , reunion_grade , reunion_wind , +#> # reunion_pres , atcf_lat , atcf_lon , atcf_grade , +#> # atcf_wind , atcf_pres , mlc_natl_lat , +#> # mlc_natl_lon , mlc_natl_grade , mlc_natl_wind , +#> # mlc_natl_pres , ds824_sh_lat , ds824_sh_lon , #> # ds824_sh_grade , ds824_sh_wind , ds824_sh_pres , #> # ds824_ni_lat , ds824_ni_lon , ds824_ni_grade , #> # ds824_ni_wind , ds824_ni_pres , bom_lat , @@ -182,8 +183,7 @@ storm_data(storm='1970143N19091') #> # jtwc_ep_lat , jtwc_ep_lon , jtwc_ep_grade , #> # jtwc_ep_wind , jtwc_ep_pres , ds824_ep_lat , #> # ds824_ep_lon , ds824_ep_grade , ds824_ep_wind , -#> # ds824_ep_pres , jtwc_cp_lat , jtwc_cp_lon , -#> # jtwc_cp_grade , ... +#> # ds824_ep_pres , … ``` Or by year, from 1842 to 2014. @@ -194,53 +194,53 @@ storm_data(year=1941) ``` ``` -#> # A tibble: 1,766 × 195 -#> serial_num season num basin sub_basin name -#> -#> 1 1940215S18149 1941 1 SP EA NOT NAMED -#> 2 1940215S18149 1941 1 SP EA NOT NAMED -#> 3 1940215S18149 1941 1 SP EA NOT NAMED -#> 4 1940215S18149 1941 1 SP EA NOT NAMED -#> 5 1940215S18149 1941 1 SP EA NOT NAMED -#> 6 1940215S18149 1941 1 SP EA NOT NAMED -#> 7 1940215S18149 1941 1 SP EA NOT NAMED -#> 8 1940215S18149 1941 1 SP EA NOT NAMED -#> 9 1940215S18149 1941 1 SP EA NOT NAMED -#> 10 1940215S18149 1941 1 SP EA NOT NAMED -#> # ... with 1,756 more rows, and 189 more variables: iso_time , -#> # nature , latitude , longitude , wind.wmo. , -#> # pres.wmo. , center , wind.wmo..percentile , -#> # pres.wmo..percentile , track_type , -#> # latitude_for_mapping , longitude_for_mapping , -#> # current.basin , hurdat_atl_lat , hurdat_atl_lon , -#> # hurdat_atl_grade , hurdat_atl_wind , hurdat_atl_pres , -#> # td9636_lat , td9636_lon , td9636_grade , -#> # td9636_wind , td9636_pres , reunion_lat , -#> # reunion_lon , reunion_grade , reunion_wind , -#> # reunion_pres , atcf_lat , atcf_lon , atcf_grade , -#> # atcf_wind , atcf_pres , ds824_sh_lat , -#> # ds824_sh_lon , ds824_sh_grade , ds824_sh_wind , -#> # ds824_sh_pres , ds824_ni_lat , ds824_ni_lon , -#> # ds824_ni_grade , ds824_ni_wind , ds824_ni_pres , -#> # bom_lat , bom_lon , bom_grade , bom_wind , -#> # bom_pres , ds824_au_lat , ds824_au_lon , -#> # ds824_au_grade , ds824_au_wind , ds824_au_pres , -#> # jtwc_sh_lat , jtwc_sh_lon , jtwc_sh_grade , -#> # jtwc_sh_wind , jtwc_sh_pres , jtwc_wp_lat , -#> # jtwc_wp_lon , jtwc_wp_grade , jtwc_wp_wind , -#> # jtwc_wp_pres , td9635_lat , td9635_lon , -#> # td9635_grade , td9635_wind , td9635_pres , -#> # ds824_wp_lat , ds824_wp_lon , ds824_wp_grade , -#> # ds824_wp_wind , ds824_wp_pres , jtwc_io_lat , -#> # jtwc_io_lon , jtwc_io_grade , jtwc_io_wind , -#> # jtwc_io_pres , cma_lat , cma_lon , cma_grade , -#> # cma_wind , cma_pres , hurdat_epa_lat , +#> # A tibble: 1,766 x 200 +#> serial_num season num basin sub_basin name iso_time nature latitude +#> +#> 1 1940215S1… 1941 1 " SP" " EA" NOT … 1940-08… " NR" -999 +#> 2 1940215S1… 1941 1 " SP" " EA" NOT … 1940-08… " NR" -999 +#> 3 1940215S1… 1941 1 " SP" " EA" NOT … 1940-08… " NR" -999 +#> 4 1940215S1… 1941 1 " SP" " EA" NOT … 1940-08… " NR" -999 +#> 5 1940215S1… 1941 1 " SP" " EA" NOT … 1940-08… " NR" -999 +#> 6 1940215S1… 1941 1 " SP" " EA" NOT … 1940-08… " NR" -999 +#> 7 1940215S1… 1941 1 " SP" " EA" NOT … 1940-08… " NR" -999 +#> 8 1940215S1… 1941 1 " SP" " EA" NOT … 1940-08… " NR" -999 +#> 9 1940215S1… 1941 1 " SP" " EA" NOT … 1940-08… " NR" -999 +#> 10 1940215S1… 1941 1 " SP" " EA" NOT … 1940-08… " NR" -999 +#> # … with 1,756 more rows, and 191 more variables: longitude , +#> # wind.wmo. , pres.wmo. , center , +#> # wind.wmo..percentile , pres.wmo..percentile , +#> # track_type , latitude_for_mapping , +#> # longitude_for_mapping , current.basin , +#> # hurdat_atl_lat , hurdat_atl_lon , hurdat_atl_grade , +#> # hurdat_atl_wind , hurdat_atl_pres , td9636_lat , +#> # td9636_lon , td9636_grade , td9636_wind , +#> # td9636_pres , reunion_lat , reunion_lon , +#> # reunion_grade , reunion_wind , reunion_pres , +#> # atcf_lat , atcf_lon , atcf_grade , atcf_wind , +#> # atcf_pres , mlc_natl_lat , mlc_natl_lon , +#> # mlc_natl_grade , mlc_natl_wind , mlc_natl_pres , +#> # ds824_sh_lat , ds824_sh_lon , ds824_sh_grade , +#> # ds824_sh_wind , ds824_sh_pres , ds824_ni_lat , +#> # ds824_ni_lon , ds824_ni_grade , ds824_ni_wind , +#> # ds824_ni_pres , bom_lat , bom_lon , bom_grade , +#> # bom_wind , bom_pres , ds824_au_lat , +#> # ds824_au_lon , ds824_au_grade , ds824_au_wind , +#> # ds824_au_pres , jtwc_sh_lat , jtwc_sh_lon , +#> # jtwc_sh_grade , jtwc_sh_wind , jtwc_sh_pres , +#> # jtwc_wp_lat , jtwc_wp_lon , jtwc_wp_grade , +#> # jtwc_wp_wind , jtwc_wp_pres , td9635_lat , +#> # td9635_lon , td9635_grade , td9635_wind , +#> # td9635_pres , ds824_wp_lat , ds824_wp_lon , +#> # ds824_wp_grade , ds824_wp_wind , ds824_wp_pres , +#> # jtwc_io_lat , jtwc_io_lon , jtwc_io_grade , +#> # jtwc_io_wind , jtwc_io_pres , cma_lat , cma_lon , +#> # cma_grade , cma_wind , cma_pres , hurdat_epa_lat , #> # hurdat_epa_lon , hurdat_epa_grade , hurdat_epa_wind , #> # hurdat_epa_pres , jtwc_ep_lat , jtwc_ep_lon , #> # jtwc_ep_grade , jtwc_ep_wind , jtwc_ep_pres , #> # ds824_ep_lat , ds824_ep_lon , ds824_ep_grade , -#> # ds824_ep_wind , ds824_ep_pres , jtwc_cp_lat , -#> # jtwc_cp_lon , ... +#> # ds824_ep_wind , ds824_ep_pres , … ``` ## shp files @@ -254,7 +254,7 @@ The result from `storm_shp()` is the path, and some metadata, for the shp files ``` #> -#> Path: /Users/sacmac/Library/Caches/rnoaa/storms/year/Year.2013.ibtracs_all_points.v03r06.shp +#> Path: /Users/sckott/Library/Caches/rnoaa/storms/year/Year.2013.ibtracs_all_points.v03r10.shp #> Basin: #> Storm: #> Year: 2013 diff --git a/vignettes/swdi_vignette.Rmd b/vignettes/swdi_vignette.Rmd index a78ccc48..c3067dd1 100644 --- a/vignettes/swdi_vignette.Rmd +++ b/vignettes/swdi_vignette.Rmd @@ -41,10 +41,10 @@ library('plyr') swdi(dataset='nx3tvs', startdate='20060505', enddate='20060506') #> $meta #> $meta$totalCount -#> numeric(0) +#> [1] 25 #> #> $meta$totalTimeInSeconds -#> [1] 0.105 +#> [1] 0.01 #> #> #> $data @@ -141,10 +141,10 @@ out <- swdi(dataset='warn', startdate='20060506', enddate='20060507', id=533623) list(out$meta, head(out$data), head(out$shape)) #> [[1]] #> [[1]]$totalCount -#> numeric(0) +#> [1] 25 #> #> [[1]]$totalTimeInSeconds -#> [1] 3.465 +#> [1] 0.014 #> #> #> [[2]] @@ -180,7 +180,7 @@ list(out$meta, head(out$data), head(out$shape)) swdi(dataset='plsr', startdate='20060505', enddate='20060510', bbox=c(-91,30,-90,31)) #> $meta #> $meta$totalCount -#> numeric(0) +#> [1] 5 #> #> $meta$totalTimeInSeconds #> [1] 0 @@ -219,10 +219,10 @@ swdi(dataset='plsr', startdate='20060505', enddate='20060510', bbox=c(-91,30,-90 swdi(dataset='nx3tvs', startdate='20060506', enddate='20060507', tile=c(-102.12,32.62)) #> $meta #> $meta$totalCount -#> numeric(0) +#> [1] 5 #> #> $meta$totalTimeInSeconds -#> [1] 0.001 +#> [1] 0 #> #> #> $data @@ -266,7 +266,7 @@ Get daily count nx3tvs features on .1 degree grid centered at `latitude = 32.7` swdi(dataset='nx3tvs', startdate='20060505', enddate='20090516', stat='tilesum:-102.0,32.7') #> $meta #> $meta$totalCount -#> numeric(0) +#> [1] 5 #> #> $meta$totalTimeInSeconds #> [1] 0.001 @@ -274,11 +274,11 @@ swdi(dataset='nx3tvs', startdate='20060505', enddate='20090516', stat='tilesum:- #> #> $data #> day centerlat centerlon fcount -#> 1 2007-09-07 32.7 -102 1 -#> 2 2008-05-27 32.7 -102 4 -#> 3 2009-04-11 32.7 -102 1 -#> 4 2007-03-29 32.7 -102 2 -#> 5 2008-06-20 32.7 -102 2 +#> 1 2007-03-29 32.7 -102 2 +#> 2 2008-06-20 32.7 -102 2 +#> 3 2007-09-07 32.7 -102 1 +#> 4 2008-05-27 32.7 -102 4 +#> 5 2009-04-11 32.7 -102 1 #> #> $shape #> shape