Dataset Information

Restoration of clipped seismic waveforms using projection onto convex sets method.

ABSTRACT: The seismic waveforms would be clipped when the amplitude exceeds the upper-limit dynamic range of seismometer. Clipped waveforms are typically assumed not useful and seldom used in waveform-based research. Here, we assume the clipped components of the waveform share the same frequency content with the un-clipped components. We leverage this similarity to convert clipped waveforms to true waveforms by iteratively reconstructing the frequency spectrum using the projection onto convex sets method. Using artificially clipped data we find that statistically the restoration error is ~1% and ~5% when clipped at 70% and 40% peak amplitude, respectively. We verify our method using real data recorded at co-located seismometers that have different gain controls, one set to record large amplitudes on scale and the other set to record low amplitudes on scale. Using our restoration method we recover 87 out of 93 clipped broadband records from the 2013?Mw6.6 Lushan earthquake. Estimating that we recover 20 clipped waveforms for each M5.0+ earthquake, so for the ~1,500 M5.0+ events that occur each year we could restore ~30,000 clipped waveforms each year, which would greatly enhance useable waveform data archives. These restored waveform data would also improve the azimuthal station coverage and spatial footprint.

SUBMITTER: Zhang J

PROVIDER: S-EPMC5155270 | biostudies-literature | 2016 Dec

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Restoration of clipped seismic waveforms using projection onto convex sets method.

Zhang Jinhai J Hao Jinlai J Zhao Xu X Wang Shuqin S Zhao Lianfeng L Wang Weimin W Yao Zhenxing Z

Scientific reports 20161214

The seismic waveforms would be clipped when the amplitude exceeds the upper-limit dynamic range of seismometer. Clipped waveforms are typically assumed not useful and seldom used in waveform-based research. Here, we assume the clipped components of the waveform share the same frequency content with the un-clipped components. We leverage this similarity to convert clipped waveforms to true waveforms by iteratively reconstructing the frequency spectrum using the projection onto convex sets method. ...[more]

PMID: 27966618

Similar Datasets

Project description:Four-dimensional computed tomography (4D-CT) is commonly used to account for respiratory motion of target volumes in radiotherapy to the thorax. From the 4D-CT acquisition, a maximum-intensity projection (MIP) image set can be created and used to help define the tumor motion envelope or the internal gross tumor volume (iGTV). The purpose of this study was to quantify the differences in automatically contoured target volumes for usage in the delivery of stereotactic body radiation therapy using MIP data sets generated from one of the four methods: (1) 4D-CT phase-binned (PB) based on retrospective phase calculations, (2) 4D-CT phase-corrected phase-binned (PC-PB) based on motion extrema, (3) 4D-CT amplitude-binned (AB), and (4) cine CT built from all available images.MIP image data sets using each of the four methods were generated for a cohort of 28 patients who had prior thoracic 4D-CT scans that exhibited lung tumor motion of at least 1 cm. Each MIP image set was automatically contoured on commercial radiation treatment planning system. Margins were added to the iGTV to observe differences in the final simulated planning target volumes (PTVs).For all patients, the iGTV measured on the MIP generated from the entire cine CT data set (iGTVcine) was the largest. Expressed as a percentage of iGTVcine, 4D-CT iGTV (all sorting methods) ranged from 83.8% to 99.1%, representing differences in the absolute volume ranging from 0.02 to 4.20 cm3; the largest average and range of 4D-CT iGTV measurements was from the PC-PB data set. Expressed as a percentage of PTVcine (expansions applied to iGTVeine), the 4D-CT PTV ranged from 87.6% to 99.6%, representing differences in the absolute volume ranging from 0.08 to 7.42 cm3. Regions of the measured respiratory waveform corresponding to a rapid change of phase or amplitude showed an increased susceptibility to the selection of identical images for adjacent bins. Duplicate image selection was most common in the AB implementation, followed by the PC-PB method. The authors also found that the image associated with the minimum amplitude measurement did not always correlate with the image that showed maximum tumor motion extent.The authors identified cases in which the MIP generated from a 4D-CT sorting process under-represented the iGTV by more than 10% or up to 4.2 cm3 when compared to the iGTVcine. They suggest utilization of a MIP generated from the full cine CT data set to ensure maximum inclusive tumor extent.

Dataset Information

Restoration of clipped seismic waveforms using projection onto convex sets method.

Publications

Restoration of clipped seismic waveforms using projection onto convex sets method.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets