Difference between revisions of "Movement Analysis"

From stgo
Jump to: navigation, search
(Download)
 
(61 intermediate revisions by 2 users not shown)
Line 1: Line 1:
We have create a toolbox for the free GIS software [http://www.openjump.org OpenJUMP] that allows to analyse movement data, called '''OpenJUMP HoRAE''' - '''Ho'''me '''R'''ange '''A'''nalysis and '''E'''stimation. The software is and can be freely distributed under the [http://www.gnu.org/licenses/gpl-2.0.txt General Public License verison 2.0 (GPL)].
+
We have create a toolbox for the free GIS software [http://www.openjump.org OpenJUMP] that allows to analyse movement data, called '''OpenJUMP HoRAE''' - '''Ho'''me '''R'''ange '''A'''nalysis and '''E'''stimation (pronounced: ''hooray''). The software is and can be freely distributed under the [http://www.gnu.org/licenses/gpl-2.0.txt General Public License verison 2.0 (GPL)] - so you can give it to your co-worker or even modify it.
 +
 
 +
'''Update: 28. Apr. 2016''' - fixed bug that prevented correct execution of Asymptote Analysis HR method<br/>
 +
'''Update: 7. Mar. 2015''' - changed some download and document links, since our file hosting provider stopped offering private file download service.<br/>
 +
'''Update: 19. Nov. 2014''' - A version packaged with OpenJUMP 1.7.1 was released. No warranties though. Please report to stefan---geosteiniger.cl if you spot differences.<br/>
 +
'''Update: 20. Feb. 2013''' - A new version has been released including a new CSV file reader, adding the GeoEllipse method by JA Long and TA Nelson (2012), and a function to extract time attributes from date fields.
  
 
== About the Home Range Analysis Toolbox ==
 
== About the Home Range Analysis Toolbox ==
Line 10: Line 15:
 
* Line Buffer
 
* Line Buffer
 
* Local Convex Hulls (LoCoH), methods: radius r, k-1 nearest neighbors, alpha region.
 
* Local Convex Hulls (LoCoH), methods: radius r, k-1 nearest neighbors, alpha region.
 +
* Potential Path Area / Geo-Ellipse (in new version only, see above)
  
 
Additionally one can:
 
Additionally one can:
* perform an Asymptote analysis for MCP, Point-KDE, Line-Buffer  
+
* Perform an asymptote analysis for MCP, Point-KDE, Line-Buffer  
 
* Calculate core areas
 
* Calculate core areas
 
* Classify home range parts
 
* Classify home range parts
* derive the Skeleton for a home range region.
+
* Derive the skeleton for a home range region.
* analyze the daily travel
+
* Analyze the daily travel
  
The toolbox was originally developed to analyse location data from grizzly bears (ursus arctos) collected with GPS collars within the Grizzly Bear Programm of the [http://foothillsresearchinstitute.ca/ Foothills Research Institute] (FRI, Alberta, Canada). Funding was provided by Canadian Phase IV [http://www.geoide.ulaval.ca/home.aspx GEOIDE] grant (#3).
+
The toolbox was originally developed to analyse location data from grizzly bears (''ursus arctos'') collected with GPS collars within the Grizzly Bear Programm of the [http://foothillsresearchinstitute.ca/ Foothills Research Institute] (FRI, Alberta, Canada). Funding was provided by Canadian Phase IV [http://www.geoide.ulaval.ca/home.aspx GEOIDE] grant (#3).
  
 
== Download ==
 
== Download ==
The software should run on Windows and Linux systems. Currently there is a bug that does not allow saving data on MacOSX systems. However, if you write me (Stefan: ssteinig) I can send you a patched version that allows saving to file put can't save to a database. Anyway - here to the downloads:
+
The software should run on '''Windows''' and '''Linux''' systems. MacOSX users can use the same version.
* [http://www.geo.uzh.ch/~sstein/ojhorae/openjump14movantools_oct2010.zip download software] (package including OpenJUMP 1.4 with the HoRAE toolbox)
+
* [http://sourceforge.net/projects/jump-pilot/files/OpenJUMP_plugins/More%20Plugins/OpenJUMP%20HORAE/OpenJUMP-171Plus_HORAE_apr2016.zip/download download software] - version from Apr. 2016 ca. 48 MB including documentation - this version is packaged with OpenJUMP 1.7.1 instead of 1.6 as for the March 2013 version; (Note, a backup copy on another server can be downloaded from [http://146.155.17.19/ojhorae/openjump171plus_moveantools_oct2014.zip here])
* [http://www.geo.uzh.ch/~sstein/ojhorae/notes.txt install instructions] -  note, the software requires [http://cds.sun.com/is-bin/INTERSHOP.enfinity/WFS/CDS-CDS_Developer-Site/en_US/-/USD/ViewProductDetail-Start?ProductRef=jre-6u22-oth-JPR@CDS-CDS_Developer Java Runtime Environment (JRE)] installed
+
*: older well-tested versions are still available: [http://horae.geosteiniger.cl/downloads/openjump_rc16moveantools_mar2013.zip HoraeMar2013] [http://horae.geosteiniger.cl/downloads/openjump143moveantools_nov2011.zip HoraeNov2011]
* a short [http://www.geo.uzh.ch/~sstein/ojhorae/oj_horae_sstein_v1_wiki.pdf article (pdf)] that describes the functions of the toolbox (article under submission)
+
* [http://horae.geosteiniger.cl/downloads/notes.txt install instructions] -  note, the software requires [http://www.oracle.com/technetwork/java/javase/downloads/index.html Java Runtime Environment (JRE)] installed
 +
* a short [http://horae.geosteiniger.cl/downloads/oj_horae_sstein_v1_wiki.pdf article (pdf)] that describes the functions of the toolbox (this is the 1st version; the final is here [http://onlinelibrary.wiley.com/doi/10.1002/wsb.168/abstract Wildlife Society Bulletin])
  
finally some preliminary results for grizzly bear data: [http://www.geo.uzh.ch/~sstein/manuscripts/sstein_gb_homerange_estimation_final.pdf GIScience 2010 extended abstract (pdf)] and [http://www.geo.uzh.ch/~sstein/manuscripts/sstein_giscience2010_hre_eval_v2_pdfscreenprint_lowres.pdf poster (pdf)]
+
Some preliminary results for a comparison of the estimators with grizzly bear data: [http://horae.geosteiniger.cl/downloads/sstein_gb_homerange_estimation_final.pdf GIScience 2010 extended abstract (pdf)]. (The link to the [http://www.geo.uzh.ch/~sstein/manuscripts/sstein_giscience2010_hre_eval_v2_pdfscreenprint_lowres.pdf poster (pdf)] seems to be dead).
 +
 
 +
A test dataset with 14 GPS points to check out the functions of the toolbox is shipped with the distribution and can be found in the folder "''/testdata/''".
 +
 
 +
'''MacOSX''' users please read: [[MacOSX Users Short Installation Guide]] - at the end we also say what keyboard-combinations are important for context menu access.
  
 
== Documentation ==
 
== Documentation ==
 +
A revised version of the user manual for OpenJUMP HoRAE is available since Feb 2013 ( [http://horae.geosteiniger.cl/downloads/horae_documentation_v1.2final.pdf download pdf]). It is also included in the software distribution. ('''Important''': In the 2011 Version was an error in the core area calculation description; but not in the calculation.)
 +
 +
The manual describes only use of the software. For more background information on the methods we recommend to have a look at the [http://horae.geosteiniger.cl/downloads/laver_abode_user_manual_2005.pdf ABODE manual] since most functions implemented were "inspired" by the ABODE toolbox (but software calculations are differently).
 +
 +
We have also created a wiki page describing [[Problems with the HoRAE Toolbox]]. We had a start with a wiki for the manual, [[Home Range Calculation Manual]], but it needs to be revised.
 +
 +
Other documents describing new methods implemented:
 +
* Steiniger, S. (2009, unpublished): ''A vector-based algorithm to detect core, edge, patch and corridor areas and comparison with its raster-based complement.'' [http://horae.geosteiniger.cl/downloads/rastervector_lepattern_sstein_v2_forGBpaper.pdf (pdf)] --- The document explains the different types of classes such as core, edge, patch and corridor that can be derived with the function <code>MOVEAN>HRE Analysis>Extract Core, Patch, Edge</code>; and: <code>MOVEAN>HRE Analysis>Extract Corridors</code>.
 +
* Steiniger, S., Hunter, A.J.S. (2013): A scaled line-based kernel density estimator for the retrieval of utilization distributions and home ranges from GPS movement tracks. ''Ecological Informatics'' 13:1-8 ([http://dx.doi.org/10.1016/j.ecoinf.2012.10.002 Elsevier link], [http://horae.geosteiniger.cl/downloads/linekde_sstein_forweb.pdf manuscript pdf]) --- The document explains how the line-based Kernel Density Estimator (Line-KDE) works and compares it with Point-based KDE and the Brownian Bridge (BB) estimator. The function is to be found under <code>MOVEAN>HRE>Line KD>...</code>.
 +
 +
If you want to '''cite/reference''' OpenJUMP HoRAE, then we ask you to use our 2012 article in the Wildlife Society Bulletin (see also above, [http://onlinelibrary.wiley.com/doi/10.1002/wsb.168/abstract WSB link]).
 +
 +
== Loading Your Data ==
 +
Currently we have no special GPS file loaders and our input files for testing have been Shapefiles (though there is a GPS extension for OpenJUMP). If you have problems to get your data into OpenJUMP, then please write us (to Stefan). Also, if you need to create the "unique days" attribute and you don't know how - let us know, so we can write a guide. There is no guide on that yet, because we had only pre-processed data and the unique days were generated with a Matlab script and not with an office application.
 +
 +
related:
 +
* [[Converting GPX files to Shapefiles]]
 +
 +
== Support and Contact ==
 +
If you have questions on the use of OpenJUMP HoRAE you can post to our user list at Google Groups:
 +
* http://groups.google.com/d/forum/ojhorae-users .
 +
Please note, that you need to login or subscribe first to send a message.
 +
 +
However, you can contact me, Stefan Steiniger, writing an email to ssteiniger(=at=)uc.cl. Alternatively you may contact [http://www.ucalgary.ca/~ahunter/ Andrew Hunter].
 +
 +
== Source Code & Future Development & Bug Reporting ==
 +
The source code for the HoRAE toolbox (an extension to OpenJUMP) can be found on SourceForge.net. OpenJUMP has there, as Jump Pilot Project, a code repository. In that SVN repository the HoRAE code can be found under:
 +
* http://jump-pilot.svn.sourceforge.net/viewvc/jump-pilot/plug-ins/horae/
 +
How to access the SVN is described [http://sourceforge.net/scm/?type=svn&group_id=118054&source=navbar here].
 +
 +
There are several ideas for future developments (e.g. to evaluate co-occurrence for two individuals). However, currently we don't have funding for that. If you wish a specific function in the toolbox, then it should be possible to implement that when you can provide funding. However, bug fixing is free of cost and depends on (my) free time. If you find a bug, then let me (Stefan) know. As the source code is freely available - others can do feature development and bug fixing too. Please let us know if you want to contribute, need help, or even write how you use HoRAE.
 +
 +
[[HoRAE TODO]]
  
We work on a wiki page with some more online documentation. However, the [http://www.geo.uzh.ch/~sstein/ojhorae/laver_abode_user_manual_2005.pdf ABODE manual] may be a good start, since most functions implemented were "inspired" by the ABODE toolbox (but softwares calculations are differently). For now you may check our [[Home Range Calculation Manual]], which contains some info on how using our toolbox.
+
== Comparison of Results with HRT ==
  
== Contact ==
+
Thanks to Christine Richter (Göttingen, GER), we got aware of the fact that results calculated with [http://flash.lakeheadu.ca/~arodgers/hre/ HRT] (by A.R. Rodgers and J.G. Kie) and OpenJUMP cannot be directly compared.
You can contact me, Stefan Steiniger, writing an email to sstein(=at=)geo.uzh.ch. Alternatively you may contact [http://www.ucalgary.ca/~ahunter/ Andrew Hunter].
+
If you use the same bandwidth parameter ''h'' for HRT and OpenJUMP HoRAE, then results (contours) will be much more detailed than those from HRT, with the latter being much smoother. However, if you take the bandwidth parameter from HRT (e.g. 63m) and multiply it with about 2.78 (e.g. 2.78 * 63 = 175m) and use this value (i.e. 175m) for HoRAE, then the contours should be about the same. The reason for the differences and the need to multiply the values is, that in ArcGIS/HRT such multiplication of h_ref is performed. This is done, as the general h_ref formula is an estimate for a Gaussian function/kernel, but will be applied to a bi-weight Kernel (in HRT). As you can see in HoRAE and on [http://en.wikipedia.org/wiki/Uniform_kernel#Kernel_functions_in_common_use Wikipedia the kernel functions] have different maximum values and properties. Internally HoRAE  does not perform a scaling so one can directly compare density results. The multiplication with the value of 2.78 for the Gaussian-based h_ref value to be used with a bi-weight kernel will however, result in comparable density estimation with a Gaussian-based density estimation (see Silverman 1986).

Latest revision as of 10:35, 2 March 2021

We have create a toolbox for the free GIS software OpenJUMP that allows to analyse movement data, called OpenJUMP HoRAE - Home Range Analysis and Estimation (pronounced: hooray). The software is and can be freely distributed under the General Public License verison 2.0 (GPL) - so you can give it to your co-worker or even modify it.

Update: 28. Apr. 2016 - fixed bug that prevented correct execution of Asymptote Analysis HR method
Update: 7. Mar. 2015 - changed some download and document links, since our file hosting provider stopped offering private file download service.
Update: 19. Nov. 2014 - A version packaged with OpenJUMP 1.7.1 was released. No warranties though. Please report to stefan---geosteiniger.cl if you spot differences.
Update: 20. Feb. 2013 - A new version has been released including a new CSV file reader, adding the GeoEllipse method by JA Long and TA Nelson (2012), and a function to extract time attributes from date fields.

About the Home Range Analysis Toolbox

OpenJUMP HoRAE Screenshot

The toolbox contains functions for visualization of movement tracks and allows the creation and analysis of animal home ranges. Home range estimation methods that have been implemented are (see also the article below):

  • Minimum Convex Polygon (MCP),
  • Kernel Density Estimation (KDE), methods: (a) href, (b) LSCV, and (c) ad-hoc (region split)
  • Line-based KDE, scaled and un-scaled
  • Brownian Bridges
  • Line Buffer
  • Local Convex Hulls (LoCoH), methods: radius r, k-1 nearest neighbors, alpha region.
  • Potential Path Area / Geo-Ellipse (in new version only, see above)

Additionally one can:

  • Perform an asymptote analysis for MCP, Point-KDE, Line-Buffer
  • Calculate core areas
  • Classify home range parts
  • Derive the skeleton for a home range region.
  • Analyze the daily travel

The toolbox was originally developed to analyse location data from grizzly bears (ursus arctos) collected with GPS collars within the Grizzly Bear Programm of the Foothills Research Institute (FRI, Alberta, Canada). Funding was provided by Canadian Phase IV GEOIDE grant (#3).

Download

The software should run on Windows and Linux systems. MacOSX users can use the same version.

Some preliminary results for a comparison of the estimators with grizzly bear data: GIScience 2010 extended abstract (pdf). (The link to the poster (pdf) seems to be dead).

A test dataset with 14 GPS points to check out the functions of the toolbox is shipped with the distribution and can be found in the folder "/testdata/".

MacOSX users please read: MacOSX Users Short Installation Guide - at the end we also say what keyboard-combinations are important for context menu access.

Documentation

A revised version of the user manual for OpenJUMP HoRAE is available since Feb 2013 ( download pdf). It is also included in the software distribution. (Important: In the 2011 Version was an error in the core area calculation description; but not in the calculation.)

The manual describes only use of the software. For more background information on the methods we recommend to have a look at the ABODE manual since most functions implemented were "inspired" by the ABODE toolbox (but software calculations are differently).

We have also created a wiki page describing Problems with the HoRAE Toolbox. We had a start with a wiki for the manual, Home Range Calculation Manual, but it needs to be revised.

Other documents describing new methods implemented:

  • Steiniger, S. (2009, unpublished): A vector-based algorithm to detect core, edge, patch and corridor areas and comparison with its raster-based complement. (pdf) --- The document explains the different types of classes such as core, edge, patch and corridor that can be derived with the function MOVEAN>HRE Analysis>Extract Core, Patch, Edge; and: MOVEAN>HRE Analysis>Extract Corridors.
  • Steiniger, S., Hunter, A.J.S. (2013): A scaled line-based kernel density estimator for the retrieval of utilization distributions and home ranges from GPS movement tracks. Ecological Informatics 13:1-8 (Elsevier link, manuscript pdf) --- The document explains how the line-based Kernel Density Estimator (Line-KDE) works and compares it with Point-based KDE and the Brownian Bridge (BB) estimator. The function is to be found under MOVEAN>HRE>Line KD>....

If you want to cite/reference OpenJUMP HoRAE, then we ask you to use our 2012 article in the Wildlife Society Bulletin (see also above, WSB link).

Loading Your Data

Currently we have no special GPS file loaders and our input files for testing have been Shapefiles (though there is a GPS extension for OpenJUMP). If you have problems to get your data into OpenJUMP, then please write us (to Stefan). Also, if you need to create the "unique days" attribute and you don't know how - let us know, so we can write a guide. There is no guide on that yet, because we had only pre-processed data and the unique days were generated with a Matlab script and not with an office application.

related:

Support and Contact

If you have questions on the use of OpenJUMP HoRAE you can post to our user list at Google Groups:

Please note, that you need to login or subscribe first to send a message.

However, you can contact me, Stefan Steiniger, writing an email to ssteiniger(=at=)uc.cl. Alternatively you may contact Andrew Hunter.

Source Code & Future Development & Bug Reporting

The source code for the HoRAE toolbox (an extension to OpenJUMP) can be found on SourceForge.net. OpenJUMP has there, as Jump Pilot Project, a code repository. In that SVN repository the HoRAE code can be found under:

How to access the SVN is described here.

There are several ideas for future developments (e.g. to evaluate co-occurrence for two individuals). However, currently we don't have funding for that. If you wish a specific function in the toolbox, then it should be possible to implement that when you can provide funding. However, bug fixing is free of cost and depends on (my) free time. If you find a bug, then let me (Stefan) know. As the source code is freely available - others can do feature development and bug fixing too. Please let us know if you want to contribute, need help, or even write how you use HoRAE.

HoRAE TODO

Comparison of Results with HRT

Thanks to Christine Richter (Göttingen, GER), we got aware of the fact that results calculated with HRT (by A.R. Rodgers and J.G. Kie) and OpenJUMP cannot be directly compared. If you use the same bandwidth parameter h for HRT and OpenJUMP HoRAE, then results (contours) will be much more detailed than those from HRT, with the latter being much smoother. However, if you take the bandwidth parameter from HRT (e.g. 63m) and multiply it with about 2.78 (e.g. 2.78 * 63 = 175m) and use this value (i.e. 175m) for HoRAE, then the contours should be about the same. The reason for the differences and the need to multiply the values is, that in ArcGIS/HRT such multiplication of h_ref is performed. This is done, as the general h_ref formula is an estimate for a Gaussian function/kernel, but will be applied to a bi-weight Kernel (in HRT). As you can see in HoRAE and on Wikipedia the kernel functions have different maximum values and properties. Internally HoRAE does not perform a scaling so one can directly compare density results. The multiplication with the value of 2.78 for the Gaussian-based h_ref value to be used with a bi-weight kernel will however, result in comparable density estimation with a Gaussian-based density estimation (see Silverman 1986).