Extract Something From an amsr Object
Extract something from the
data slot of an
[[ method works for all
oce objects, i.e.
objects inheriting from
oce-class. The purpose,
as with the related replacement method,
is to insulate users from the internal details of
objects, by looking for items within the various storage
slots of the object. Items not actually stored can also be
extracted, including derived data, units of measurement,
and data-quality flags.
The method uses a two-step process to try to find the
requested information. First, a class-specific function
is used to try to access the requested information (see
“Details of the specialized ... method”).
Second, if no match is found, a general function is used
(see ‘Details of the general method’). If neither
method can locates the requested item,
NULL is returned.
## S4 method for signature 'amsr' x[[i, j, ...]]
Character string indicating the name of item to extract.
Optional additional information on the
Optional additional information (ignored).
Partial matches for
are permitted for
j is ignored for
Data within the
data slot may be found directly, e.g.
j="SSTDay" will yield sea-surface temperature in the daytime
j="SSTNight" is used to access the nighttime data. In
j="SST" yields an average of the night and day values
(using just one of these, if the other is missing). This scheme works for
all the data stored in
amsr objects, namely:
rain. In each case, the default
is to calculate values in scientific units, unless
which case the raw data are returned.
"raw" mode can be useful
in decoding the various types of missing value that are used by
as.raw(255) for land,
a missing observation,
as.raw(253) for a bad observation,
as.raw(252) for sea ice, or
as.raw(251) for missing SST
due to rain or missing water vapour due to heavy rain. Note that
something special has to be done for e.g.
because the idea is that this syntax (as opposed to specifying
"SSTDay") is a request to try to find good
data by looking at both the Day and Night measurements. The scheme
employed is quite detailed. Denote by "A" the raw value of the desired field
in the daytime pass, and by "B" the corresponding value in the
nighttime pass. If either A or B is 255, the code for land, then the
result will be 255. If A is 254 (i.e. there is no observation),
then B is returned, and the reverse holds also. Similarly, if either
A or B equals 253 (bad observation), then the other is returned.
The same is done for code 252 (ice) and code 251 (rain).
In all cases, the returned value is a matrix with
NA values inserted at locations where
the raw data equal
as.raw(251:255), as explained
Details of the general method
If the specialized method produces no matches, the following generalized
method is applied. As with the specialized method, the procedure hinges
first on the value of
First, a check is made as to whether
i names one of the standard
oce slots, and returns the slot contents if so. Thus,
x[["metadata"]] will retrieve the
i is a string ending in the
"Unit", then the characters preceding that string
are taken to be the name of an item in the data object, and a list
containing the unit is returned. This list consists of an item
unit, which is an
an item named
scale, which is a string describing the
measurement scale. If the string ends in
" unit", e.g.
x[["temperature unit"]], then just the expression is returned,
and if it ends in
" scale", then just the scale is returned.
i is a string ending in
"Flag", then the corresponding
data-quality flag is returned (or
NULL if there is no such flag).
x[["salinityFlag"]] returns a vector of salinity
x is a ctd object.
If none of the preceding conditions are met, a check is done
to see if the
metadata slot contains an item with
the provided name, and that is returned, if so. A direct match
is required for this condition.
data slot is checked to see if it contains
an item with the name indicated by
i. In this case, a partial
match will work; this is accomplished by using
If none of the above-listed conditions holds, then
NULL is returned.
Other functions that extract parts of
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17
## Not run: # Show a daytime SST image, along with an indication of whether # the NA values are from rain. library(oce) earth <- read.amsr("f34_20160102v7.2.gz") fclat <- subset(earth , 35 <= latitude & latitude <= 55) fc <- subset(fclat , -70 <= longitude & longitude <= -30) par(mfrow=c(2,1)) plot(fc, "SSTDay") rainy <- fc[["SSTDay", "raw"]] == as.raw(0xfb) lon <- fc[["longitude"]] lat <- fc[["latitude"]] asp <- 1 / cos(pi*mean(lat)/180) imagep(lon, lat, rainy, asp=asp) mtext("red: too rainy to sense SSTDay") ## End(Not run)
Want to suggest features or report bugs for rdrr.io? Use the GitHub issue tracker.