Adrien/ACOUSTIQUE-PNR_BARONNIES-2020
Adrien/ACOUSTIQUE-PNR_BARONNIES-2020
1
0
Fork 0
Code Issues Pull Requests Projects Releases Wiki Activity

Addition of parallel calculation scripts for the indices of P1 & P2 on the cluster

Browse Source
This commit is contained in:
A-d-r-i 2020-07-28 16:54:55 +02:00
parent 2a96d6aa48
commit 1e4b2d9d75
2 changed files with 74 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

43
CLUSTER/indices-P1.r Normal file
View File

@ -0,0 +1,43 @@
library(parallel)
soundRep<-"../sound/"
ResultRep<-"../result/"
ListSound<-dir(soundRep)
#### PARALLELE ####
no_cores<-detectCores()
c2<-makeCluster(no_cores)
parLapply(c2,ListSound, function(x){
library(tuneR)
library(seewave)
source("AcouIndexAlpha.r")
soundRep<-"../sound/"
ResultRep<-"../result/"
ListSound<-dir(soundRep)
TableTotal<-NULL
a=0
b=1
for (i in 1: 6)
{
wave<-readWave(paste(soundRep,x, sep="/"),from=a,to=b,units="minutes")
Result<-AcouIndexAlpha(wave, stereo=FALSE, min_freq = 1000, max_freq = 22000,anthro_min = 300,
anthro_max = 2000,bio_min=2000,bio_max=18000, wl=512, j=5,
AcouOccupancy=FALSE,Bioac=TRUE,Hf=FALSE,Ht=FALSE,H=FALSE,ACI=TRUE,
AEI_villa=FALSE,M=FALSE,NDSI=TRUE,ADI=FALSE,NP=TRUE)
Result$Mono_left$ID<-paste(x,i, sep="_")
TableTotal<-rbind(TableTotal,Result$Mono_left)
a=a+10
b=b+10
}
write.table(TableTotal,file=paste(ResultRep,"/",x,".csv", sep=""), sep=";", row.names = F)
})

31
CLUSTER/indices-P2.r Normal file
View File

@ -0,0 +1,31 @@
library(parallel)
soundRep<-"../sound/"
ResultRep<-"../result/"
ListSound<-dir(soundRep)
#### PARALLELE ####
no_cores<-detectCores()
c2<-makeCluster(no_cores)
parLapply(c2,ListSound, function(x){
library(tuneR)
library(seewave)
source("AcouIndexAlpha.r")
soundRep<-"../sound/"
ResultRep<-"../result/"
ListSound<-dir(soundRep)
wave<-readWave(paste(soundRep,x, sep="/"),from=0,to=1,units="minutes")
Result<-AcouIndexAlpha(wave, stereo=FALSE, min_freq = 1000, max_freq = 22000,anthro_min = 300,
anthro_max = 2000,bio_min=2000,bio_max=18000, wl=512, j=5,
AcouOccupancy=FALSE,Bioac=TRUE,Hf=FALSE,Ht=FALSE,H=FALSE,ACI=TRUE,
AEI_villa=FALSE,M=FALSE,NDSI=TRUE,ADI=FALSE,NP=TRUE)$Mono_left
#write.table(TableTotal,file=paste(ResultRep,"/",x,".csv", sep=""), sep=";", row.names = F)
Result$ID<-x
save(Result, file=paste(ResultRep,"/",x,".Rdata", sep=""))
})