Monday, 9 March 2015

Some minor enhancements in qProf

qProf is a QGIS plugin for the creation of topographic and geologic profiles. Some minor enhancements were added to qProf in the new, experimental releases 0.2.9 and 0.3.0.

In version 0.2.9 the following options were added:
- the calculation of the absolute slope along the profile; previously only the directional slope was calculated;
- the possibility to flip a profile horizontally (option "Reverse profile direction"), as well as to reverse the orientation of the x axis (option "Reverse x axes direction) were added.

See examples in Figs. 1-4 (below).


Fig. 1. Map of Mt. Alpi zone (Basilicata, Southern Italy), with profile in red (direction from left to right). DEM data: TINITALY.
Fig. 2. Topographic (top) and absolute slope (bottom) profiles (see map trace in Fig. 1).

Fig. 3. Same profile as in Fig. 2, but with profile direction reversed.

Fig. 4. Same profile as in previous figures, but with reversed x axes directions.

In version 0.3.0, released on 2015-03-08, the possibility of configuring the figure save parameters (width, resolution, font size, subplot configuration parameters) was added. These graphic parameters can be saved in a text file and subsequently loaded for being applied to future plots.
This option is available from Export - Figure. Fig. 5 represents the window for defining the parameters, saving or loading them, and also for saving the figure as a file in PDF, SVG or TIF formats.

Fig. 5. Window for saving a figure as a graphic file, and for configuring graphic parameters, and saving or loading them.


The two versions are experimental since the impact of reversing the profile on the geological plotting has not been fully tested. For geological plots, the suggested version is 0.2.8.

Gerrit Tombrink, PhD student at the Goettingen University (Germany), is thanked for proposing the implemented options.

Saturday, 3 January 2015

Cross-section intersections of line and polygon layers with qProf in QGis

A new release of qProf, plugin for QGis, allows to determine the intersections of line and polygon layers, representing for instance faults and geological outcrops, on a vertical cross-section.
These new functionalities complement those already available for projecting geological attitudes and lines on a cross-section.
The resulting data can be saved as shapefiles and csv.

Example of cross-section representing geological polygon layer intersections with a topographic profile.

The plugin can be installed/updated via the plugin manager of QGis. Alternatively it can be downloaded from the QGis plugin repository:
http://plugins.qgis.org/plugins/qProf/
or from the Bitbucket repository:
https://bitbucket.org/mauroalberti/qprof/downloads



Tuesday, 11 November 2014

R is progressing fast in TIOBE Index

According the TIOBE Index for November 2014, R, an open source statistical language,  progressed from position 15 in October to position 12 in November, thus heading for entering in the "top 10". It is preceded by Visual Basic and Visual Basic .NET. For instance, Python is at position 7.
The much younger Julia language is at position 126 while MATLAB is # 24.
More details:
http://www.tiobe.com/index.php/content/paperinfo/tpci/index.html


Thursday, 25 September 2014

Changes in qgSurf plugin tools

Due to persisting plugin incompatibilities, related to Matplotlib, I've created a new qgSurf version, # 0.3.4. This version is made up by just two tools: "Best fit plane" and "DEM-plane intersection".
The removed 3D surfaces simulation and deformation tools are now available in a pure Python project, "simSurf", developed and tested in Linux Mint and published at: https://github.com/mauroalberti/simSurf.

Friday, 12 September 2014

Simulation and deformation of georeferenced geological surfaces: a pure Python implementation

A pure Python project for simulation and deformation of georeferenced geological surfaces is presented at:
https://github.com/mauroalberti/simSurf
It is an updated and Python-based, independent version of the simulation and deformation modules originally presented in the qgSurf plugin (vers. 0.3.3) for Quantum GIS. It is developed and tested in Linux Mint.
Two modules, "geosurface_simulation.py" and "geosurface_deformation.py" make possible to simulate and deform georeferenced surfaces (see example in Fig. below).


Example of a sheared and rotated sinusoidal surface with geographic parameters matching those of the Mt. Alpi - Mt. Raparo Aster DEM (Lucania, Southern Italy). The view is from NE to SW. Mt. Raparo is at the right. 3D visualization created with ArcScene (ESRI).




Thursday, 28 August 2014

Tuesday, 24 June 2014

New version of qgSurf for QGIS

A new version of qgSurf, 0.3.3, has been released. It adds the possibility to use point or line layers as sources for the best-fit-plane calculation.
The plugin should work also in QGIS 2.2.0 for Windows 8 - 64 bit, even if the QGIS installer complains about a lacking dependency, Tkinter, a problem that I will investigate in the next days. If however you accept the installation, close QGIS and then reopen it, QGIS will present the plugin as installed and available. No problem of this kind has been observed in Windows Vista and Ubuntu.

Friday, 13 June 2014

Backward pathline calculation for steady vector fields in Quantum GIS

Calculating backward pathlines in steady vector fields is possible using the new version of VectorFieldCalc, a Python plugin for Quantum GIS. The pathline calculation is implemented via the Runge-Kutta-Fehlberg method (see references at the end).

For forward calculations, the user defines positive values of time steps and total times. For backward calculations, use negative values for both total time and time steps. An example for the last case, using the example data provided in the plugin folder, is below.





The result of the backward pathline calculation is presented in the figure below, with pathline points symbolised by increasingly negative total time, together with labels of the total times. The flows in this Antarctic glacier in Terra Nova Bay (data from D. Biscaro PhD thesis) are from top-left to bottom-right. Higher velocities of the flow field are in red, lower velocities in blue. Starting points are represented by red squares in the bottom-left section of the glacier.



 The plugin can be installed via the usual procedure for Quantum GIS, using the plugin manager.

For additional details (based on the previous version) you can see:
Vector Field Processing
'Vector field processings': un plugin Quantum GIS per l'analisi di campi vettoriali 2D (in Italian)


Wednesday, 28 May 2014

Released first beta version of gvSIG CE 1.0.0 (64 bit OS)

From gvSIG CE developer announcement:

The first beta version of gvSIG CE 1.0.0 has just been released.
Downloads for Linux, Mac OS X and Windows here:
https://sourceforge.net/projects/gvsigce/files/beta-1/
All provided packages are for 64 bit operating systems.

Read the full announcement at:
http://gvsigce.sourceforge.net/wiki/index.php/GvSIG_CE_1.0.0_beta_1_release_notes

Sunday, 25 May 2014

Plot of geological attitudes and traces in vertical sections along fold axes with Quantum GIS and qProf

The new released version of qProf (0.2.4), a Python plugin for Quantum GIS, allows to project geological attitudes and traces onto a vertical section along fold axes. On-the-fly reprojection of layers is supported.

The first step in using this plugin is creating a topographic profile by defining one (or more) source DEMs and a line layer that define the section, plus a line densify distance (generally of the order of the used DEM cell resolution) and an optional line order field (required only when the input geometry of the line layer is not sufficiently simple, see previous posts on the qProf plugin).
While for creating only topographic profiles the input line can be made up of many points and more than one DEM can be used, for geological projections the input line must be made up by just two points and only one DEM may be selected.
After having defined and created the profile using the tools in the "profile" tab, it is then possible to project geological attitudes or traces by using the commands in the "Project" tab.


Projection of geological attitudes

Geological attitudes will be stored in a point layer, with the measured dip direction and dip angle of the features in two distinct fields of the attribute table.
The commands are available from the "Geological attitudes" subwindow.
The projection of geological attitudes may be applied via different methodologies: (a) the nearest intersection point between the geological plane of a feature and the section plane; (b) using a common projection ("fold") axis for all features; (c) using individual axes for each attitude, with two numeric fields storing trend and plunge values.
If some data are selected, only those selected will be plotted, otherwise all layer features will be used. Labels can also be added to projected measures.
It is possible to cumulate in a single analysis session more sub-analyses, distinguishing them via different colors. Previous sub-analyses can be removed with the 'Reset geological attitudes' button.
Data can be saved as csv files or shapefiles.


A plot example, with attitudes labelled by original trend and plunge values, is in the figure below.



Projection of line traces

Line traces can represent the intersection of geological surfaces with the topography. When these surfaces are cylindrically folded, it is possible to projected the traces onto the section via a user-defined fold axis. The elevation information will be derived from the chosen DEM.
The commands are available from the "Geological traces" subwindow.
A line layer has to be chosen, with the name of the id field and a value for the line densification.
Trend and plunge of the projection axis will be defined, and a color chosen for the line plotting.
Different subanalyses, with different projection axes, can be differentiated by using different colors.
Previous analyses can be eliminated with the "Reset traces" command.
Results can be saved as text file, storing id-s-z values (where s is the horizontal distance from the section start) for each resulting point.



An example of the output is below. In this case, only geological traces are plotted. In general, both geological traces and attitudes can be plotted in the same section.



Installation

qProf can be installed from Quantum GIS via the plugin manager, or downloaded and unzipped in the Python plugin folder from:


For bug segnalation:
or write directly to alberti.m65@gmail.com (also for improvement requests or questions).
 

Wednesday, 21 May 2014

Friday, 9 May 2014

Interactive map of small particle air pollution in world cities

In the Guardian web site an interactive map, based on the data of the World Health Organization, is available:
http://www.theguardian.com/news/datablog/interactive/2014/may/08/exposure-air-pollution-city-map
North America has lower pollution levels than in European states, with south-eastern European countries having the highest pollutions for Europe. India, China, Mongolia, Pakistan, Afghanistan, Iran present the highest levels in the world, together with many Middle East cities.

Monday, 5 May 2014

Traduzione del manuale di QGis in italiano: vuoi contribuire con una donazione?

"QGIS è un software open source molto potente sviluppato interamente da volontari. Ma cosa sarebbe un software senza un'accurata documentazione? 
Una nuova versione di QGIS viene rilasciata ogni 4 mesi e di conseguenza la documentazione cambia ed è necessario aggiornare il manuale con le nuove funzionalità in modo che tutti siano in grado di sfruttarlo al massimo.

Siamo un gruppo di volontari che nel tempo libero crea, mantiene e traduce la documentazione di QGIS. Aiutaci con un piccolo contributo!"