if (FALSE) {
dir.create("images")
## Create image sequence
traj <- simulTrajec(path="images",
nframes=30,nIndividuals=20,domain='square',
h=0.01,rho=0.9,movingNoise=TRUE,
parsMoving = list(density=20, duration=10, size=1,
speed = 10, colRange = c(0,1)),
sizes=runif(20,0.004,0.006))
## Load images
dir <- "images"
allFullImages <- loadImages (dirPictures=dir,nImages=1:30)
stillBack <- createBackground(allFullImages,method="mean")
allImages <- subtractBackground(stillBack)
partIden <- identifyParticles(allImages,threshold=-0.1,
pixelRange=c(3,400))
# select the nframes with the most identified particles
nframes <- 3
frames <- order(tapply(partIden$patchID,partIden$frame,length),
decreasing=TRUE)[1:nframes]
mId <- manuallySelect(particles=partIden,frame=frames)
}
Run the code above in your browser using DataLab