#toCWT: (aux)
##NOTE: renvoie une matrice 3D
-toCWT <- function(X, sw= 0, tw= 0, swabs= 0, nvoice= 12, noctave= 5, s0= 2, w0= 2*pi,
- lt= 24, dt= 0.5, spectra = FALSE, smooth = TRUE, scaled = FALSE, scalevector)
+toCWT <- function(X, sw=0, tw=0, swabs=0, nvoice=12, noctave=5, s0=2, w0=2*pi,
+ lt=24, dt=0.5, spectra=FALSE, smooth=TRUE, scaled=FALSE, scalevector)
{
noctave <- adjust.noctave(lt, dt, s0, tw, noctave)
if(missing(scalevector))
resArray
}
+#from sowas
+adjust.noctave <- function(N,dt,s0,tw,noctave)
+{
+ if (tw>0)
+ {
+ dumno <- as.integer((log(N*dt)-log(2*tw*s0))/log(2))
+ if (dumno<noctave)
+ {
+ cat("# noctave adjusted to time smoothing window \n")
+ noctave <- dumno
+ }
+ }
+ noctave
+}
+
#from sowas
cwt.ts <- function(ts,s0,noctave=5,nvoice=10,w0=2*pi)
{
}
#NOTE: vect2mat = as.matrix ?! (dans aux.R)
-vect2mat <- function(vect)
+vect2mat <- function(vect, delta, lscvect)
{
vect <- as.vector(vect)
matrix(vect[-(1:2)], delta, lscvect)
smCWT <- function(CWT, sw= 0, tw= 0, swabs= 0, nvoice= 12, noctave= 2, s0= 2, w0= 2*pi,
lt= 24, dt= 0.5, scalevector )
{
-#noctave <- adjust.noctave(lt, dt, s0, tw, noctave)
-#scalevector <- 2^(0:(noctave * nvoice) / nvoice) * s0
+ #noctave <- adjust.noctave(lt, dt, s0, tw, noctave)
+ #scalevector <- 2^(0:(noctave * nvoice) / nvoice) * s0
wsp <- Mod(CWT)
smwsp <- smooth.matrix(wsp, swabs)
smsmwsp <- smooth.time(smwsp, tw, dt, scalevector)
Xcwt2[i,] <- c(delta, lscvect, Xcwt4[,,i] / max(Mod(Xcwt4[,,i])) )
#rm(conso, Xcwt4); gc()
+ #Benjamin: FIX is this OK ?
+ lscvect = dim(Xcwt4)[2]
+
## _.b WER^2 distances ########
Xwer_dist <- matrix(0.0, n, n)
for(i in 1:(n - 1))
{
- mat1 <- vect2mat(Xcwt2[i,])
+#browser()
+##ERROR là sans FIX lscvect :: delta lscvect --> taille ??!
+ mat1 <- vect2mat(Xcwt2[i,], delta, lscvect)
- for(j in (i + 1):n)
+ for(j in (i + 1):n)
{
- mat2 <- vect2mat(Xcwt2[j,])
+ mat2 <- vect2mat(Xcwt2[j,], delta, lscvect)
num <- Mod(mat1 * Conj(mat2))
WX <- Mod(mat1 * Conj(mat1))
WY <- Mod(mat2 * Conj(mat2))
}
}
diag(Xwer_dist) <- numeric(n)
- Wwer_dist
+ Xwer_dist
}