We’ll visualize several bathymetric datasets. Bathymetry data documents the terrain of the floors of water bodies. In our case, we’ll produce maps of the ocean floor in several parts of the Atlantic Ocean.

These datasets come in XYZ coordinate formats and must be converted to .bathy formats before they can be visualized using the wireframe function from the lattice library. The two packages we use in this exercise are marmap and lattice.

First, we look at an included dataset of the marmap package.

## 
## Attaching package: 'marmap'
## 
## The following object is masked from 'package:grDevices':
## 
##     as.raster

plot of chunk unnamed-chunk-1

Now, we look at three randomly sampled patches of ocean floor, taken from UCSD Scripps Institution of Oceanography’s global topography data extractor (link: http://topex.ucsd.edu/cgi-bin/get_data.cgi)

##      V1     V2    V3
## 1 359.5 -10.11 -5408
## 2 359.5 -10.11 -5392
## 3 359.5 -10.11 -5382
## 4 359.6 -10.11 -5386
## 5 359.6 -10.11 -5412
## 6 359.6 -10.11 -5456

plot of chunk unnamed-chunk-2

##      V1  V2 V3
## 1 300.5 -30 43
## 2 300.5 -30 43
## 3 300.5 -30 41
## 4 300.6 -30 39
## 5 300.6 -30 39
## 6 300.6 -30 39

plot of chunk unnamed-chunk-2

The last bathymetric map looks particularly interesting, so we expanded our coordinates for the area where it was sampled from and produce the map below. We can clearly see some sort of trench feature, as well as a potential underwater mountain.

##      V1  V2 V3
## 1 298.0 -28 83
## 2 298.0 -28 83
## 3 298.0 -28 81
## 4 298.1 -28 81
## 5 298.1 -28 81
## 6 298.1 -28 79

plot of chunk unnamed-chunk-3