read.oceancolor: Read satellite data from Oceancolor hdf files

In r-forge/satin: Functions for downloading, reading and displaying satellite data for oceanographic applications with R

Description

Remote sensing data downloaded from the Oceancolor web page (see details) for a user defined area is extracted from hdf files

Usage

read.oceancolor(file, lons, lats)

Arguments

file	name of the hdf file(s) to be read
lons	a vector with western- and eastern-most longitudes
lats	a vector with southern- and northern-most latitudes

Details

This function relies on the h5r package, therefore only hdf version 5 (v5) files can be read. It has been tested with files with data from Aqua Modis and SeaWiFS sensors, downloaded from the oceancolor web page (http://oceancolor.gsfc.nasa.gov/cgi/l3). These hdf files are available as v4, but they were converted to v5 with h4h5tools (http://www.hdfgroup.org/h4toh5). The user must specify the minimum and maximum latitude and longitude limits within the interval (-90, 90) and (-180, 180) respectively.

Value

An object of class “satin” (a list) with components:

longitude	a vector of longitudes of the area of interest (aoi)
latitude	a vector of latitudes of the aoi
rs.data	an array of remote sensed data
rs.name	full name of remote sensed variable
rs.units	variable measurement units
period	date or averaging period
itype	a string indicating the type of product

Author(s)

Héctor Villalobos and Eduardo González-Rodríguez

Examples

## Not run: 
 
 # read a single file
 sst <- read.oceancolor("A20101822010212.L3m_MO_SST4_4.h5", 
                        lats=c(20, 30), lons=c(-130, -105))

 # read several files residing in a directory 
 setwd("D:/data")
 fi <- list.files(pattern = glob2rx("*.h5"))
 lats <- c(20, 25)
 lons <- c(-115, -110)
 mSST <- read.oceancolor(fi, lats, lons)      

 # plotting the first processed file
 plotSatin(mSST) 

 # plotting the second processed file
 plotSatin(mSST, image = 2)

## End(Not run)

r-forge/satin documentation built on Feb. 18, 2022, 2:39 a.m.