/
orsifronts.Rmd
68 lines (51 loc) · 2.05 KB
/
orsifronts.Rmd
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
---
title: "Southern Ocean Fronts for R"
author: "Michael Sumner"
date: "`r Sys.Date()`"
output: rmarkdown::html_vignette
vignette: >
%\VignetteIndexEntry{Southern Ocean Fronts for R}
%\VignetteEngine{knitr::rmarkdown}
\usepackage[utf8]{inputenc}
---
```{r setup, include=FALSE}
knitr::opts_chunk$set(collapse = TRUE, comment = "#>", fig.width = 7, fig.height = 7)
```
Orsi fronts for R
---------------------------
The "Orsi fronts" are a commonly used climatological map of Southern Ocean frontal distributions, once published on a
Southern Ocean Atlas Database Page after Orsi et al. (1995).
The data is provided in this package as a built-in data set. The object `orsifronts` is a 'SpatialLinesDataFrame' object, as defined in the 'sp' package. (A 'SpatialLinesDataFrame' is a "GIS vector" layer where complex line geometry objects are linked to simple attribute data in a table form, and it behaves like a 'data.frame').
```{r}
library(orsifronts)
cols <- hcl(seq(0, 240, length = nrow(orsifronts)), c = 50)
plot(orsifronts, col = cols, lwd = 2)
degAxis(1)
degAxis(2)
box()
```
There are `r nrow(orsifronts)` line objects, with some simple names and codes.
```{r}
as.data.frame(orsifronts)
plot(orsifronts, xlim = c(60, 180), col = cols, asp = 1/cos(55 * pi / 180), lwd = 4)
legend("topleft", sprintf("%s (%s)", orsifronts$front, orsifronts$name),
col = cols, lwd = 4, cex = 0.8, bty = "n")
degAxis(1)
degAxis(2)
box()
```
Park/Durand fronts for R
---------------------------
Similarly for the more recent Park/Durand fronts derived from altimetry remote sensing.
```{r}
as.data.frame(parkfronts)
pcols <- hcl(seq(0, 240, length = nrow(parkfronts)), c = 50)
plot(parkfronts,col = pcols, asp = 1/cos(55 * pi / 180), lwd = 4)
legend("topleft", sprintf("%s (%s)", parkfronts$front, parkfronts$name),
col = cols, lwd = 4, cex = 0.8, bty = "n")
degAxis(1)
degAxis(2)
box()
```
# References
Orsi, A. H., T. Whitworth III and W. D. Nowlin, Jr. (1995). On the meridional extent and fronts of the Antarctic Circumpolar Current, Deep-Sea Research I, 42, 641-673.