Learn R Programming

oce (version 0.9-19)

sectionSmooth: Smooth a Section

Description

Smooth a section in the lateral (alpha version that may change).

Usage

sectionSmooth(section, method = c("spline", "barnes"), debug = getOption("oceDebug"), ...)

Arguments

section
A section object containing the section to be smoothed. For method="spline", the pressure levels must match for each station in the section.
method
Specifies the method to use; see ‘Details’.
debug
A flag that turns on debugging. Set to 1 to get a moderate amount of debugging information, or to 2 to get more.
...
Optional extra arguments, passed to either smooth.spline or interpBarnes.

Value

An object of section-class that ordered in some way.

Details

This function should be used with caution, as should any operation that changes data. Although smoothing may be desirable to produce aesthetically-pleasing plots, it can also introduce artifacts that can lead to erroneous conclusions. The prudent analyst starts by comparing plots of the raw data with plots of the smoothed data. For method="spline", the section is smoothed using smooth.spline on individual pressure levels, with any parameters listed in parameters being passed to that function. If df is not present in parameters, then this function sets it to the number of stations divided by 5. Smoothing is done separately for temperature, salinity, and sigma-theta.

For the (much slower) method="barnes" method, smoothing is done across both horizontal and vertical coordinates, using interpBarnes. Any arguments in ... being passed to that function; see ‘Examples’.

See Also

Other things related to section data: [[,section-method, [[<-,section-method, as.section, handleFlags,section-method, plot,section-method, read.section, section-class, sectionAddStation, sectionGrid, sectionSort, section, subset,section-method, summary,section-method

Examples

Run this code
library(oce)
data(section)
gs <- subset(section, 109<=stationId&stationId<=129)
gsg <- sectionGrid(gs, p=seq(0, 5000, 150))
gss1 <- sectionSmooth(gsg, "spline", df=16)
plot(gss1)
gss2 <- sectionSmooth(gsg, "barnes", xr=24, yr=100)
plot(gss2)

Run the code above in your browser using DataLab