Behandling af rå accelerometer data

NB: Vi har flyttet vores supportsider til vores SENS hjemmeside. Indholdet her vil ikke længere blive opdateret. For de nyeste opdateringer, henviser vi til følgende side.

Baggrund

Denne artikel indeholder eksempler på den indledende loading af rå accelerometer data til RStudio eller Python. 

Begge eksempler behandler en .bin fil.

For at eksportere en .bin fil, se denne artikel.

RStudio

Dette er et eksempel på et script i RStudio, der importerer og organiserer accelerometer data.

#######
#Import af data
#######

fileBlock <- readRaw("22_1.HEX", width=24, signed=TRUE)
data<-blockValue(fileBlock)
data2<-as.data.frame(data)
df2 <- as.data.frame(matrix(data2[,1], byrow=TRUE, ncol = 25))
df2$zero<-"0x"

# Organise time variable

df2$time <- paste(df2$zero,df2$V1,df2$V2,df2$V3,df2$V4,df2$V5,df2$V6,df2$V7,df2$V8,df2$V9,df2$V10,df2$V11,df2$V12, sep="")
df2$time<-as.numeric(df2$time)
options(digits.secs=6)
df2$time2<-as.POSIXct(df2$time/1000, origin="1970-01-01", tz="GMT")

# Organise X axis data
df2$x<-paste(df2$V13,df2$V14,df2$V15,df2$V16, sep="")
df2$x2 <- strtoi(df2$x, base=16)
df2$x3 <- ifelse(bitwAnd(df2$x2, 0x8000) > 0, -0xFFFF-1, 0) + df2$x2
df2$x3<-df2$x3*0.0078125

# Organise Y axis data
df2$y<-paste(df2$V17,df2$V18,df2$V19,df2$V20, sep="")
df2$y2 <- strtoi(df2$y, base=16)
df2$y3 <- ifelse(bitwAnd(df2$y2, 0x8000) > 0, -0xFFFF-1, 0) + df2$y2
df2$y3<-df2$y3*0.0078125

# Organise Z axis data
df2$z<-paste(df2$V21,df2$V22,df2$V23,df2$V24, sep="")
df2$z2 <- strtoi(df2$z, base=16)
df2$z3 <- ifelse(bitwAnd(df2$z2, 0x8000) > 0, -0xFFFF-1, 0) + df2$z2
df2$z3<-df2$z3*0.0078125
###
finaldata<-data.frame(time=df2$time2, x=df2$x3, y=df2$y3, z=df2$z3)
data<-finaldata

Python

Denne Python function læser en .bin fil og gemmer det i et NumPy array med 4 kolonner. Kolonne 0 er timestamp, kolonne 1-3 er x, y, z. X, y og z skaleres med 0.0078125 G.