test/julia/test-findpeaks.R

In kongdd/PhenoAsync: Extract Remote Sensing Vegetation Phenology

julia_source("I:/Research/phenology/phenofit.R/inst/julia/findpeaks.jl")

{
    x <- seq(0, 1, len = 1024)
    pos <- c(0.1, 0.13, 0.15, 0.23, 0.25, 0.40, 0.44, 0.65, 0.76, 0.78, 0.81)
    hgt <- c(4, 5, 3, 4, 5, 4.2, 2.1, 4.3, 3.1, 5.1, 4.2)
    wdt <- c(0.005, 0.005, 0.006, 0.01, 0.01, 0.03, 0.01, 0.01, 0.005, 0.008, 0.005)

    pSignal <- numeric(length(x))
    for (i in seq(along=pos)) {
        pSignal <- pSignal + hgt[i]/(1 + abs((x - pos[i])/wdt[i]))^4
    }
    pracma::findpeaks(pSignal, npeaks=3, threshold=4, sortstr=TRUE)
}

system.time({
    for (i in 1:1e3) {
        # r = pracma::findpeaks(pSignal, npeaks=3, threshold=4, sortstr=FALSE)
        r = julia_call("findpeaks", pSignal, npeaks=3L, threshold=4, sortstr=TRUE)
    }
})