NETVOLC (Version 1.1) : User Manual
NETVOLC (Version 1.1) : User Manual
1]
http://cediac.uncu.edu.ar/software
User Manual
CONICET
National Scientific and Technical
Research Council
Argentina
July 2014
Table of Contents
1. INTRODUCTION .......................................................................................................... 3
1.1 Related Documents ...................................................................................................................... 3
1.2 Contacts ....................................................................................................................................... 3
3. SETUP ............................................................................................................................ 6
3.1 Download ..................................................................................................................................... 6
3.2 Requirements ............................................................................................................................... 6
3.3 External function installation ....................................................................................................... 6
3.4 NETVOLC installation ................................................................................................................ 7
3.5 NETVOLC parameters ................................................................................................................ 7
4. USAGE ......................................................................................................................... 10
2
1. Introduction
NETVOLC is an IDL command-line application that computes the outline or
boundary of volcanic edifices. It uses a Minimum Cost Flow algorithm jointly with
several products derived from a Digital Elevation Model (DEM). An in depth
description of the NETVOLC algorithm and its limitations can be found in the
documentation indicated in the Related Documents section.
This document is intended for those who want to install and use the NETVOLC
application. It contains step by step instructions on how to run it.
1.2 Contacts
Leonardo D. Euillades
Instituto CEDIAC – Facultad de Ingeniería
Universidad Nacional de Cuyo
Phone number: +54 261 4135000 ext 2108
E-mail: [email protected]
Pablo Grosse
CONICET & Fundación Miguel Lillo
Phone number: +54 381 4239723 ext 106
E-mail: [email protected]
Pablo A. Euillades
Instituto CEDIAC – Facultad de Ingeniería
Universidad Nacional de Cuyo
Phone number: +54 261 4135000 ext 2108
E-mail: [email protected]
For better assistance, please provide us with a detailed description of the error.
3
2. Overview and Functional Description
NETVOLC is composed of several IDL programs that are compiled using the IDL
version 7.0.1. As mentioned in the Introduction section, NETVOLC also uses
external functions and programs. Regarding external functions, NETVOLC
employs several functions from the ENVI software for extracting and/or
generating DEM-derived products. Concerning external programs, NETVOLC
uses two:
1. A DEM of the whole study area (i.e. a DEM with different volcanoes where the
user wants to apply the NETVOLC algorithm). Input DEM must be in the ENVI
format and have an associated ENVI header file.
In this mode, NETVOLC creates a singular folder for each volcano present in the
list and saves outputs in this folder. The volcanoes list must have a particular
format. It consists of a text file where each row corresponds to a volcano to be
processed, and each row contains the following parameters separated by blank
spaces:
Where:
Volcano: user-defined name/id associated to the volcano (the folder created by
NETVOLC containing outputs will have this name).
Latitude: Latitude (or Northing) of a central point within the volcano (e.g. summit, center of
crater) in decimal degrees (or meters). Latitude or Northing will depend on the geographic
projection of the supplied DEM.
Longitude: Longitude (or Easting) of a central point within the volcano (e.g. summit,
center of crater) in decimal degrees (or meters). Longitude or Easting will depend on the
geographic projection of the supplied DEM.
winSize_x: Dimension (x) in number of pixels of the window to be extracted from the input
DEM.
4
winSize_y: Dimension (y) in number of pixels of the window to be extracted from the input
DEM.
Note that the extracted window should cover the volcano of interest completely, as it
defines the maximum possible extension for the resulting outline, but should not be too
large, as processing time directly depends on window size; preferentially, inclusion of
neighboring volcanoes/landforms should be avoided. The window will be centered on
the central point defined above by Latitude (or Northing) and Longitude (or Easting).
5
3. Setup
3.1 Download
The NETVOLC package (netvolc_v1.1.zip) can be download from:
http://cediac.uncu.edu.ar/software/
3.2 Requirements
Linux
6
Note: executable will be installed in a folder named bin.Linux inside mcf-
1.3 folder. Created files must be: libmcf.a, mcf and mcflight
6. Copy compiled files to a bin folder in your Linux distribution (i.e.
/usr/local/bin) in order to be accessible to anyone.
Windows
NOTE: In version 1.1 mcf-1.3.tgz has been updated in order to obtain greater
compatibility with 64-bits systems. It is not necessary to install mcf-1.3.tgz again if
it has been installed previously.
7
The following are the parameters associated to NETVOLC. They are present in
the file named input_param.txt, which must be located in the folder where
NETVOLC will be executed. The user should set up/modify the parameters as
desired before running NETVOLC; see Euillades et al. (2013) for details.
NOTE: Don’t use Windows Notepad for modifying the input-param.txt file.
Instead use Notepad++, Microsoft WordPad o any other text processor.
[BATCH]
dem=Mariana_UTM55N_dem ;DEM file name
volcanoesList=volcanoes.txt ;Volcanoes file name
[ENVIRONMENT]
appMCF=mcflight ;MCF solver
appGREP=grep ;grep command
[NETVOLC]
net_plot_triang=y ;flag indicating if network triangulation plot is desired
(y/n). Default is y.
polar_plot_triang=y ;flag indicating if polar coordinate network
triangulation plot is desired (y/n). Default is y.
perc_threshold=100 ;percentage of arcs to be considered based on its
cost. Default is 100.
perc_threshold_arcs=99 ;percentage of arcs to be considered based on its
length. Default is 99.
bin_size=1 ;bin size for estimating histogram (see idl histogram
documentation). Default is 1.
ang_threshold=5 ;border to be considered in the polar coordinate
system for linking fictitious nodes to the network (in
degrees). It means that nodes with an azimuth polar
coordinate value within (−180; -180+ang_threshold)
degrees will be linked to the supply node. Nodes with
an azimuth polar coordinate value within (180-
ang_threshold; 180) degrees will be linked to the
demand node. Default is 5.
smooth_width=2 ;smoothing factor of the final solution: default is 2.
min_threshold_separation=0 ;minimum radial distance (in pixel units) from the
central point of the solution. Default is 0.
max_threshold_separation=200 ;maximum radial distance (in pixel units) from the
central point of the solution. Default is equal or
greater to window size extracted from the DEM.
use_dem=n ;flag indicating if height information must be
considered (y/n). Default is n.
min_dem_height_filter=2000 ;minimum height of the solution (if use_dem flag is
set to y).
max_dem_height_filter=4300 ;maximum height of the solution (if use_dem flag is
set to y).
use_slope=n ;flag indicating if slope information must be
8
considered (y/n). Default is n.
slope_factor=1. ; factor establishing the influence of slope on the
solution (when slope_factor flag is set to y). Default is
1.
iterations=10 ;number of iterations to be performed. Default is 10.
overflow=1.e+07 ;number considered an overflow. It is used for
filtering arcs with an associated cost greater than this
value. Default is 1.e+07.
stdev_factor=1 ;number of standard deviations to be considered
when statistics are calculated. Default is 1.
cancel_strategy=a ; It is used for reducing the network in every iteration.
a: indicates network must be reduced cancelling arcs.
n: indicates network must be reduced cancelling
nodes.
Default is a.
[ENVI]
KERNEL=3 ;kernel to be used for estimating DEM derived
products. Default is 3.
SHADED_RELIEF_AZIMUTH=45 ;azimuthal angle to be used for calculating shaded
relief product. Default is 45.
SHADED_RELIEF_ELEVATION=45 ;elevation angle to be used for calculating shaded
relief product. Default is 45.
[ENVIRONMENT]
appMCF= c:\cygwin\bin\mcflight ;MCF solver
appGREP= c:\cygwin\bin\grep ;grep command
9
4. Usage
Once installed, we recommend testing NETVOLC with the supplied test data
package in order to verify the correct installation and functionality. The simple
procedure indicated next for processing the test data can be used as a usage
guide for the final user.
In order to process a different area than the example, the following steps must be
reproduced:
10
1. Make available a DEM of your area of interest and place it inside the netvolc
folder.
2. Create a volcanoes file (volcanoes.txt) as indicated in the previous section.
Note: every line in this folder starting with a # character will be skipped.
3. Set up the input_param.txt as desired (see above).
4. Follow steps 4-8 of the test data procedure.
NOTE: if a failure occurs during NETVOLC execution (i.e. mcflight is not found),
NETVOLC will be stopped with an error. In order to restart the NETVOLC
execution, the following steps must be carried out just after failures have been
identified and solved:
11