BioVeL ESW DIFF - ENM Statistical Workflow with raster difference computation

Created: 2013-07-08 15:06:00      Last updated: 2016-06-22 14:45:16

The ENM Statistical Difference Workflow (ESW DIFF) allows the computation of the extent and intensity of change in species potential distribution through computation of the differences between two raster layers using the R statistical environment (R Core Team 2013). The difference file is computed from two input files (in this case present projection and 2050 projection) coming from the Ecological Niche Modelling (ENM) Workflow (http://www.myexperiment.org/workflows/3355).

The difference between each corresponding raster cell value is computed and stored in the difference file, regardless of the input files’ geographical extent and origin. If the files have a different geographical extent and/or origin the raster-diff workflow automatically crops them to the same extent and resamples the values using the “nearest neighbour” method, resulting in a perfect cell match between the two rasters. The resulting value in the difference file ranges from the negative and positive maximum values of the input files to represent the maximum possible change in both directions.

The diff result values are classified into six positive and negative classes. The predicted differences are presented for each species as a heat map, where cells with colors from green to red indicate an increase and from green to blue a decrease of predicted potential for a species. The workflow enables the computation of the raster centerpoints, overall coverage, overall intensity and the difference in intensity or coverage between two raster layers, as well as the distance between the centerpoints of the two input files (shift vector). The centerpoint is computed as the weighted average of raster cells inside the mask, overall coverage as the percentage of raster cells with values >0, overall intensity as the sum of all valued cells divided by the number of raster cells, and finally the shift vector as the distance between the centerpoints of the two input files using the Vincenty’s formulae. Summary statistics for the raster layers are also provided: number of raster cells, mean intensity, median, CV, SD, Min value, Max value and number of cells outside the mask.

This workflow store the computed difference image to the BioVeL Geoserver and create PNG overwiews contains the input and output files and the distribution points. The difference result image could be shown in the BioSTIF web visualisation client

The workflow contains the computation methods:

- Computation of the difference layer:

difference = inputfile_B - inputfile_A

- Computation of the overall coverage and intensity overall_coverage = ((ncell(x)-cellStats(x, 'countNA')-count(x,0))/(ncell(x)-cellStats(x, 'countNA')))*100 overall_intensity = cellStats(x, sum)/(ncell(x)-cellStats(x, 'countNA'))

- Used standard R functions:

cellStats: https://r-forge.r-project.org/scm/viewvc.php/pkg/raster/R/cellStats.R?view=markup&root=raster

The PNG overview will be created with these values (corresponding to the Geoserver style) ERDASImagine file:

breakpoints = c(0,1,25,50,75,100) GeoTIFF file breakpoints = c(0,2,63,127,190,254) colors = c(rgb(1.0,1.0,1.0), rgb(0.996,0.898,0.851), rgb(0.988,0.682,0.569), rgb(0.984,0.416,0.290), rgb(0.871,0.176,0.149), rgb(0.647,0.059,0.082))

Expected as input files are those raster images as a URL reference e.g. output files from openModeller created by an ENM workflow. *A file upload into the workflow is not possible, only a reference to the input file as URL*

Supported formats are:

- ERDASImagine (.img) files with a value range 0 - 100, noData value 101 - GeoTIFF (.tif) files with a value range 0- 254, noData Value 255

Input files with a different extent will be intersected to a common extent and the current layer (inputfile 1) will be resampled to the prediction layer (inputfile 2) using the "nearest neighbor" sampling method.

The difference layer will be computed as inputfile_2 - inputfile_1

For working with the ESW in a own localhost an R instance should by installed on a own machine (http://cran.r-project.org/). For the computation the followed extensions must be installed: rgdal, raster and rserve. Rserve must be started before the computation.

--------------------

This workflow has been created by the Biodiversity Virtual e-Laboratory (BioVeL http://www.biovel.eu/) project. BioVeL is funded by the EU’s Seventh Framework Program, grant no. 283359.

Information Preview

Information Run

Run this Workflow in the Taverna Workbench...

Option 1:

Copy and paste this link into File > 'Open workflow location...'
http://www.myexperiment.org/workflows/3647/download?version=16
[ More InfoExpand ]


Information Workflow Components

Information Authors (1)
Information Titles (1)
Information Descriptions (1)
Information Dependencies (0)
Inputs (6)
Processors (2)
Beanshells (1)
Outputs (10)
Datalinks (17)
Coordinations (0)

Information Workflow Type

Taverna 2

Information Uploader

Information License

All versions of this Workflow are licensed under:

Information Version 16 (latest) (of 16)

View version:

Information Credits (1)

(People/Groups)

Information Attributions (0)

(Workflows/Files)

None

Information Tags (12)

Log in to add Tags

Information Shared with Groups (2)

Information Featured In Packs (0)

None

Log in to add to one of your Packs

Information Attributed By (1)

(Workflows/Files)

Information Favourited By (0)

No one

Information Statistics

 

Citations (0)

None


Version History

In chronological order:



Reviews Reviews (0)

No reviews yet

Be the first to review!



Comments Comments (0)

No comments yet

Log in to make a comment




Workflow Other workflows that use similar services (0)

There are no workflows in myExperiment that use similar services to this Workflow.