DigitSeis v1.5 QuickStart Guide
0Downloading DigitSeis and Other Resources

DigitSeis is available at http://www.seismology.harvard.edu/research/DigitSeis test.
html. It is provided as a .zip file containing a standalone package for either Mac OS
X or Windows, and so does not require installation of MATLAB1 . In the package,
README.txt contains information for installing and uninstalling DigitSeis2 .
A sample image of an analog seismogram is available for download; it was used to produce the step-by-step examples found in this guide. A comprehensive guide is alsoavailable both within the bundles and as an individual download. Videos demonstrating usage of DigitSeis step-by-step and answers to frequently asked questions areavailable as online resources.

1Loading and Preparing Images

Upon opening DigitSeis, the following window (Figure 1) appears.

Figure 1: File loading options presented upon opening DigitSeis.

To continue work on an existing analysis, load a .mat file by selecting Load Existing
Analysis. The user can then pick up from where the loaded analysis was saved.
1
2For MATLAB source codes, please contact the Harvard Seismology Group.
Step-by-step instructions for installation can also be found in the comprehensive DigitSeis manual.

1

To upload a new image3 , select Open an Image. The image may take some timeto load. Once it has been loaded, an image with light-coloured traces over a dark
background should be visible as in Figure 2.

I(A)

(B)

(C)

(D)

(E)

(F)

II

Figure 2: Main DigitSeis window after loading a new image. At top (I), Save (A), Crop(B), Time Mark Direction (C), Measure Width (D), Measure Offset (E), and navigation(F) tools. At left (II), fields to input time mark dimensions and begin classification.

New images must be cropped to eliminate space without data. Select the Crop tool
(Figure 2, A), then click and drag to create a box encompassing the desired area.
Finish cropping the image by double clicking anywhere within the box.
Before classifying, time mark properties must be defined. Using tools D and E in
Figure 2, measure the approximate offset and width of a time mark4 (as shown inFigure 3) and input the results in the corresponding fields in the left sidebar (Figure2, II). Ensure that the Time Mark Direction icon (Figure 2, C) displays a time markin the correct orientation; click the icon to toggle between positions5 .

345DigitSeis works with all image formats supported by the MATLAB function imread, detailed athttps://www.mathworks.com/help/matlab/import export/supported-file-formats.html.
To make this process as easy as possible, zoom in on a single time mark to measure it. See thecomprehensive DigitSeis manual for tips on using navigation tools.
For images without time marks, see the comprehensive DigitSeis manual.

2

IIIIII

IVFigure 3: Measuring time mark width and offset. (I) Tool measuring only horizontal
distance for width, (II) tool measuring only vertical distance for offset, (III) toggle indicatingtime marks are below traces for the image, and (IV) fields in left sidebar for inputting time
mark dimensions.

2Classifying Objects

Once the image has been cropped and time mark dimensions have been input, pressing Classify Objects in the left sidebar (Figure 2, II) opens a window where objectsdetected by DigitSeis can be classified, as below (Figure 4). It takes some time toopen.
For the image to be digitized, the objects within it must be properly labeled. Thougha majority of objects should already be correctly classified, the tools6 displayed inFigure 4 can be used to ensure proper classification of all objects:(A) Small Region Analysis. Creates a box which can be double-clicked to open
a new window to more closely examine a region of the image7 .
(B) Remove Pixels. The first of these two tools masks individual pixels in a linedrawn by the user; the second masks an area indicated by a shape drawn bythe user.
67For in-depth information about these and other tools, see the comprehensive DigitSeis manual.
More information regarding options in the Small Region Analysis window can be found in thecomprehensive DigitSeis manual.

3

I
(A)

(B)

(C)

(D)

(E)

(F)

(G)

II

Figure 4: DigitSeis classification window. At top, icons for Small Region Analysis (A),Remove Pixels (B), Combine Objects (C), Overlap Analysis (D), Undo (E), Classify Objects(F), and Show Time Marks (G).

(C) Combinine 0bjects. Allows for selection of multiple objects to be defined
as a single object. Once all desired objects have been selected (temporarily
highlighted in yellow), right-click to complete.
(D) Overlap Analysis. Opens a new window displaying only a single selected object. Allows for closer examination and masking, or designation of overlapping
regions8 .
(E) Undo. Undoes only one, most recent, past action. Appears as a green arrow
when an action can be undone. Cannot be redone.
(F) Classify Objects. Introduces crosshairs to define objects as noise (red), timemarks (blue), or traces (white). Once objects have been defined, right click toend selection.
(G) Show Time Marks. Boxes all defined time marks with yellow for ease of error
detection.
8More information regarding options in the Overlap Analysis window can be found in the comprehensive DigitSeis manual.

4

3

Digitizing and Correcting Traces

Once all objects in the image have been correctly classified, the classification window
can be closed. This automatically reopens the main DigitSeis window.

IIIFigure 5: Main DigitSeis window after creation of digitized traces. Traces, zero lines, andgaps are visible. At left (I), the sidebar contains options for displaying zero lines, traces,and gaps, and for correcting traces. A zoomed in view of the digitized traces (II) displays
correctable gaps, indicated in yellow.

Back in the main DigitSeis window, pressing Digitize in the left sidebar produces
digitized traces and zero lines9 along the traces in the image, as in Figure 5. Any gaps
in the digitized traces are boxed in yellow if the gaps visible box, below CorrectTrace in the left sidebar, is checked.
To correct a gap in a digitized trace, first select Correct Trace, which introducescrosshairs that can be used to select the trace to be corrected. A box, outlined inred, subsequently appears; this box can be adjusted to encompass the area to becorrected. Double clicking anywhere inside the red-outlined box opens a window tocorrect the trace, as in Figure 6.

9

To troubleshoot extra zero lines or misplaced traces, see the comprehensive DigitSeis manual.

5(A)

(B)

(C)

(D)

Figure 6: Correct trace window. At top, icons for pinpointing tools (A), Show Original
Image (B), Unmask (C), and Digitize Trace (D).

The Correct Trace window allows for masking, unmasking, and reclassification ofobjects in the selected area of the trace. Gaps in data existing in the original media
should remain, but any gaps created by errors in classification or digitization can becorrected in the Correct Trace window.
Once all objects are properly masked and classified, or a point-traced curve has beencreated and classified10 , the Digitize Trace tool (Figure 6, D) can be used to create adigitized trace. That digitized trace is then applied to the final trace either with orwithout NaN values11 , after which the window can be closed.

4

TimingOnce the digitized traces are correct, press Calculate Timing to initiate calculationof the timing anchors. The time marks are indicated by yellow marks (timebars)
within traces, and as red timebars at the ends of traces, as in Figure 7.
10

DigitSeis’s pinpointing feature is frequently used for thin, high-amplitude traces that overlap
others. For more information about using this tool, see the comprehensive DigitSeis manual.
11If NaNs are included, all NaN values for the segment of the trace in the Correct Trace window
are applied to the trace.

6(A)
III

III(B)

Figure 7: At top (A), main DigitSeis window after calculating timing. Zoomed views
shown in center are (I) the right-click options for a timebar, (II) the window for setting
absolute time of selected time bar, and (III) absolute time appearing in red once set. Atbottom (B), after time has been updated, times at 15-minute intervals line up vertically.

7

An absolute time can then be assigned to one or more timebars. This is done byright-clicking a timebar and choosing Set Time, which opens a window where theabsolute time corresponding to the timebar can be set (Figure 7, A, II). Once enough12
timebars have been set, Update Time in the left sidebar becomes clickable. Whenclicked, it assigns absolute time to each timebar that is displayed for every 15 minutes(Figure 7, B).

5

Saving Results

Data can be saved from DigitSeis in two ways – as a .mat file and as a .SAC file.
An analysis can be saved as a .mat file at any time by clicking the save icon at top(Figure 2, A). When the analysis is first saved, a window opens where information
about the analysis13 can be input. Once that window is closed, and for all subsequentsaves, the standard save window for the user’s operating system opens to allow forfile naming and choice of save location.
To extract the data as a .SAC file, allsteps for digitization must first be complete. Then, Extract SAC in theleft sidebar can be clicked, which opens
a window (Figure 8) where station information can be input. Once Generate SAC Files is clicked, the user
can choose where to store the .SAC files
through the operating system’s save window.
Throughout the process of digitizing arecord, frequent saves to new .mat files
are recommended, to make backtracking
easy in case of an error.
Figure 8: Window for inputting
station information before saving.SAC files.

12

As long as no timebars are missing, only one timebar should need to be set; for more informationon other situations, see the comprehensive DigitSeis manual.
13The window collects information for the user signature. Other information about the analysis
can be accessed through the Show Information icon in the top left corner of the main DigitSeiswindow.

8

This manual is provided to complement the DigitSeis v1.5 package for digitizationof analog seismograms available on the
http://www.seismology.harvard.edu/research/DigitSeis.html
DigitSeis manuals and supplementary materials prepared by Katelyn Lee for theHarvard Seismology Group.
For more information on the Harvard Seismology Group, please visit:
http://www.seismology.harvard.edu/index.html.
For any questions, comments, or bug reports, please contact Thomas Lee, graduate
student in the Harvard Seismology Group at thomasandrewlee@g.harvard.edu.
The development of DigitSeis was partially supported by the United States
Geological Survey Earthquake Hazard Program through the following grants:
G14AP00016
G16AP00021
G17AP00007
G18AP00018
and by the Miller Institute for Basic Sciences at the University of California,
Berkeley.

9DigitSeis uses two functions that were not written by the authors nor are built-infunctions that come with MATLAB. Permissions from the functions’ authors toinclude them in the DigitSeis package were obtained prior to inclusion. The
individual functions can be obtained from sources below.

rgb2hsv fast.m
This function was written by Dr. Alexander Ihlow (Technische UniversitatIlmenau), and is available at http://www.mathworks.com/matlabcentral/
fileexchange/15985-fast-rgb2hsv/content/rgb2hsv fast.m.

writesac.m
This function was written by Dr. David Yang (Los Alamos National Laboratory).
The latest version of this function was obtained from the author, but it can also befound included in the Seismic Analysis Code (SAC) distribution, at
http://ds.iris.edu/ds/nodes/dmc/forms/sac/.
Copyright statement for this function:
Copyright, 19, The Board of Governors of the Los Alamos National Security, LLC.
This software was produced under a U. S. Government contract(DE-AC52-06NA25396) by Los Alamos National Laboratory, which is operated by
the Los Alamos National Security, LLC for the U. S. Department of Energy. The U.
S. Government is licensed to use, reproduce, and distribute this software.
Permission is granted to the public to copy and use this software without charge,
provided that this Notice and any statement of authorship are reproduced on allcopies. Neither the Government nor the LANS makes any warranty, express orimplied, or assumes any liability or responsibility for the use of this software.

10