This project proposes new maps for background navigation in SeisComP applications like scolv. Improvements from basic SeisComP3 map package are:

• better global resolution using full scale ETOPO1 data
• possibility to zoom on specific targets at higher resolution, using SRTM data
• possibility to change a lot of mapping parameters (land/sea colormap, lighting, contrast, etc...)

Example of level 0 zoom (ETOPO1 data).

Example of level 3 zoom (ETOPO1 data, 1.8km resolution).

Example of level 5 zoom (SRTM3 data, 90m resolution).

Example of level 8 zoom (SRTM1 data, 30m resolution).

Download the file at seiscomp3-ipgp-maps.tgz (566 Mb).

This tar archive contains 370 tiles named to be used as background maps in SeisComP3. Includes zoom level 8 for few selected areas (Antilles, Indonesia, Réunion, ...).

To install it, just do (seiscomp must be installed and configured):

cd $SEISCOMP_ROOT/share/maps
tar zxf seiscomp3-ipgp-maps.tgz

Matlab/Octave source code to produce the maps.

• Octave 4.2.0 (or newer) or Matlab 2011b (or newer)
• some functions from mapping-matlab toolbox are necessary, in particular readhgt, dem and ibil
• ETOPO1 data, available at NGDC/NOAA: download the file etopo1_bed_g_i2.zip
• SRTM3 and SRTM1 data will be automatically downloaded by the code (needs internet connection)

Downloads codes and data and place it in a directory.

The code is a single function without argument. Some variable must be ajusted to proper local values, editing the source mksc3maps.m:

X.etopo = 'data/etopo1_bed_g_i2'; % ETOPO1 base filename (.bin and .hdr) 
X.psrtm3 = 'data/SRTM3'; % directory to write SRTM3 downloaded files
X.psrtm1 = 'data/SRTM1'; % directory to write SRTM1 downloaded files

It is mandatory to have two separated directories for SRTM1 and SRTM3 since they use the same filename. Once .hgt files are written, they won't be downloaded again from internet if the code is run again.

Default behavior of the code will make only level 4 zoom tiles using ETOPO1. You can change this

maxlevel = 4;

• 4 is default with tiles of 22.5 x 11.25 ° using ETOPO1 with bathymetry
• 6 means tiles of about 6 x 3 ° using SRTM3 (land only)
• 8 means tiles of about 1.4 x 0.7 ° using SRTM1 (land only)

Note that if maxlevel 4 is 341 total tiles, level 6 corresponds to about 1800 tiles and level 8 will make near 29000 tiles (not recommanded). Since SeisComP3 allows multiscale tiling, it is much preferred to define specific targets where you want high-resolution maps. The code will produce level 5 to 8 zoom on focused tiles defined by targets with coordinates longitude,latitude, e.g.:

targets = [ ...
   -61.6636, 16.0444;  % Soufrière, Guadeloupe
   -61.1685, 14.8113;  % Mt. Pelée, Martinique
    55.7141,-21.2486;  % Piton de la Fournaise, Réunion
    43.3588,-11.7514;  % Khartala volcano, Comores    
];

Note that targets cannot be outside latitude 60S-60N (no SRTM data).

Graphic options can be changed by editing the code script:

seacolor = [linspace(51,144)',linspace(79,161)',linspace(122,178)']/255;
landcolor = [linspace(193,230)',linspace(194,230)',linspace(159,230)']/255;

Those two variables are 100x3 matrix of RGB values to define submarine and land colormaps, respectively. These colormaps attempt to reproduce the original SeisComP3 colors. Comment these lines and use seacolor.m and landcolor.m from mapping-matlab toolbox to try colorful maps, or use any other colormap your want.

X.optdem = {'noplot','latlon','zlim',[-1e4,1e4],'landcolor',landcolor,'seacolor',seacolor,'lakezmin',0,'interp'};

This cell vector contain options for the main dem.m function that produces the lighting relief.

• 'noplot' is to avoid figure plot of each maps on the screen during the process;
• 'latlon' is mandatory to get correct z-ratio lighting;
• 'zlim' fixes the limits (min and max in meter) of colormap, must be set to avoid color discrepencies between tiles;
• 'landcolor' sets the colormap to use for land (z > 0);
• 'seacolor' sets the colormap to use for sea (z <= 0);
• 'lakezmin' fills any flat areas with sea color, for elevations > 0 (see level 8 zoom screenshot for example);
• 'interp' interpolates gaps (no values) in the data (useful for SRTM in some high relief areas)

Lighting options are default (45°W direction, linear contrast). See dem.m documentation to change them and for possible additional arguments.

From a terminal, run:

octave --eval 'mksc3maps'

will produce all world*.png images in the current directory maps/. Images are created only if not already exist. To rebuild any image, just delete it before running the script.

Image files have to be copied to $SEISCOMP_ROOT/share/maps/ to take effect in SeisComP3.

François Beauducel, IRD/IPGP, beaudu, beauducel@ipgp.fr

Ali A. Fahmi, IRD Jean-Marie Saurel, IPGP

Type "doc mksc3maps" or see the code mksc3maps.m for details.

SeisComP3 background maps, SeisComP3 documentation