Package index • gmsp

Time-series conversion

End-to-end workflows that produce consistent AT (acceleration), VT (velocity) and DT (displacement) time series from a single input type. Each combines regularisation, STFT parameter selection, anti-alias resampling, edge tapering and spectral-domain integration/differentiation.

AT2TS(): Construct acceleration, velocity, and displacement histories
VT2TS(): Construct acceleration, velocity, and displacement from velocity
DT2TS(): Construct acceleration, velocity, and displacement from displacement
TSL2TSW(): Reshape canonical long time series to wide columns
TSW2TSL(): Reshape wide time-series columns to canonical long rows
normalizeTS(): Normalize a long time-series table to unit amplitude per channel.

Decomposition and spectra

Empirical / variational mode decomposition with mode-based filtering, and elastic SDOF response spectra including D50 and D100 horizontal spectra through TSL2PS.

TS2IMF(): Decompose one time series into intrinsic mode functions.
TSL2PS(): Convert canonical long time series to response spectra.
PSL2PSW(): Convert long response spectra to wide form.
PSW2PSL(): Convert wide response spectra to long form.

Intensity measures

TSL2IM(): Compute intensity measures from a canonical long time-series table.
getIntensity(): Compute intensity measures from a long time-series table.
IML2IMW(): Cast canonical long intensity measures to wide columns.
TSL2AI(): Normalize cumulative Arias-intensity prefixes.
TSL2CAV(): Normalize cumulative absolute-velocity prefixes.
TSL2CAV5(): Normalize gmsp's 0.05-g thresholded CAV prefixes.

Spectral matching

In-memory algorithm core of the spectral-matching workflow: per-record modal shaping against a target band, suite-level record scaling, and compatibility diagnostics. The batch counterpart is runMatch().

fitModalFactor(): Fit the modal shaping factors of one record
fitRecordFactor(): Fit the suite record-scaling factors
matchCompatibility(): Compatibility diagnostics of a matched suite

Provider-format parsers

Read one record file from a single provider format and return a long (t, OCID, s) table. Used by parseRecord() and reusable on their own.

readAT2(): Read a PEER NGA-West2 AT2 acceleration record.
readV2(): Read a CESMD V2 acceleration record (multi-channel V2 or single-channel V2c).
readV2A(): Read a NWZ V2A acceleration record (3D-BLOCK, 1 file = 3 components).
readAC(): Read a 3D-COL acceleration record (ACA, ACB, LIS).
readTR(): Read a TRA/TRZ/TRB/TRC acceleration record (GSC and SGC families).
readTwoCol(): Read a 2-column whitespace-delimited ASCII record.

Record extraction pipeline

Parse one record across components, classify channels by direction (H1/H2/UP), align lengths, scale to canonical units (mm), and write raw/AT.<RecordID>.csv plus its sidecar JSON.

parseRecord(): Parse one record (event x station x owner) into a LONG time-series table.
mapComponents(): Map provider components to canonical processed components.
alignComponents(): Equalize NP across components of a parsed record.
extractRecord(): Extract one record to raw/<KIND>.<RecordID>.csv + <KIND>.<RecordID>.json.
archiveRawOwner(): Compress raw.owner/ to raw.owner.tar.gz and delete the directory.

Indexing

Build per-owner raw file, record, and intensity index tables from the per-record CSV/JSON pairs produced by extractRecord().

buildRawFileTable(): Build the per-owner RawFileTable CSV (provider file inventory, post-archive safe).
buildRawRecordTable(): Build the canonical RawRecordTable for one or more owners.
buildRawIntensityTable(): Build the canonical RawIntensityTable for one or more owners (WIDE).
getRawIntensities(): Compute KIND-derivable intensities for one raw record.

Selection and IO

selectRecords(): Select records from a record metadata table.
writeSelection(): Write a selection table to selection/<name>.csv / .json.
readAT(): Read acceleration records via readTS() with kind = "AT".
readTS(): Read parsed time-series records into the shape AT2TS() / VT2TS() / DT2TS() expect.
readVT(): Read velocity records via readTS() with kind = "VT".
readDT(): Read displacement records via readTS() with kind = "DT".
readISEE(): Read a Micromate ISEE blasting record.

Project runner wrappers

Optional batch runners for file-based record sets. Use direct helpers for in-memory work; use these wrappers when a project JSON declares the input and output contract.

runSelect() runGMSP() runTrim() runIMF() runPSA() runProcess() runStage0() runMatch() runPlot() runExport(): Run a gmsp batch pipeline

Audit

auditSite(): Audit site / station information in a record metadata table.
auditDistances(): Audit distances in a record metadata table.
auditParsers(): Audit parsers: dry-run parseRecord on every record of an owner.