Hi list
I have a vector, which are remotely sensed chlorophyll values for a certain
pixel in 11 years, then i have a flag or label vector (t_max) that consists
of 0 and 1, which tells me where i have the annual peak (one per year, then
eleven 1 and all the rest are 0).
I'm interested in extracting the date in which the bloom beggins. For that
matter I'm using a threshold of 5% above the median of the whole series.
I need to go over the chlorophyll data vector (cla) and every time I find
the yearly maximum value (where t_max == 1), i need to go backwards in cla
as many time steps as needed untill I cross the threshold (5% over the
median of cl)... once i crossed that threshold, i need to be sure that at
least for 2 time steps, the value of cla keeps going down (it's below the
threshold)... once that condition is met, i want to get the first position
where cla is higher than that threshold...
that would be the phytoplankton bloom starting date...
This is more or less what i was trying, but it does not work and i'm
stuck... can you please help me??
cla<- ts of 506 values (every 46 values, i have one year)
t_max<- vector of 0 & 1; 1 where the yearly cla max occurs (11 ones & all
rest 0)
threshold = 0.05*median(cla)+median(cla)
for ( t in 1:length(t_max) ) {
if ( t_max[t] == 0 )
{
next;
} else {
if ( (cla[t-1] < threshold) && (cla[t-2)] < threshold) )
{
return which(cla[t])
}
}
}
Thanks a lot in advance!! :)
Vero
[[alternative HTML version deleted]]