[Hyperspec-help] Plotting of stacked spectra

Cornelia Hett cornelia.hett at gmail.com
Mon Jul 18 10:16:57 CEST 2011


Hi there,

I'm a total newbie with R and hyperSpec. So I struggle with some basic
plotting functions.

Background / dataset:
-------------------------------
I am analyzing land cover changes using hyper-temporal satellite imagery: I
have a time series of images over 9 years with imagery every 16 days. So in
total 207 spectral values for one x-/y-coordinate. There's only one spectral
band.

I have around 200 calibration points in the landscape, all with land cover
type "Paddy rice", but from different years (2001 to 2009). I've created a
hyperSpec object "rice" from a text-file with the following command:

rice <- read.txt.wide (file = "pdy.txt", cols = list (x = "X-coord", y =
"Y-coord", ident = "ID_ALL", yr = "Year", cnum = "CNUM", cl = "LC-Class",
spc="EVI spectral values", .wavelength = "Day in year"), header = TRUE,
check.names = FALSE)

So my hyperSpec object has 6 other describing variables, and then the
spectra.

Problem:
-------------
In odrer to visually analyze if my data spectral data varies from year to
year I would like to create a stacked plot with grouping of the spectra per
year. How do I do this?

I've tried different things, I think the trivk i to use stacked = ... in the
correct way.
>plotspc(rice, spc.nmax = 200, stacked = rice$YEAR)
... does not work.

Thanks for your help,
Cornelia
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.r-forge.r-project.org/pipermail/hyperspec-help/attachments/20110718/b6d9e0bf/attachment.htm>


More information about the hyperspec-help mailing list