V.krige GSoC 2009: Difference between revisions

From GRASS-Wiki
Jump to navigation Jump to search
(renamed todo section)
No edit summary
Line 23: Line 23:
Question: another layout can be a Choicebox with the three packages and in respect of the choice the drawing of corresponding widgets. Opinions welcome.  
Question: another layout can be a Choicebox with the three packages and in respect of the choice the drawing of corresponding widgets. Opinions welcome.  


On GRASS side, I succeeded in populating a wxChoicebox with vector layers of the mapset in use, furthermore discarding maps with no points. This second filter requires some time if there are lot of maps, so I'm considering to first set up the interface and in the meanwhile the user inspects the options perform this operation. It is not dead long (few seconds on my 2GB-RAM and Intel Core2Duo laptop, using spearfish dataset), but it creates an evident lag between the call of the module and the creation of interface.
The lsit of point layers is provided by VectorSelect class (gselect module) to provide uniformity of layout with wxGUI. The filtering is done once at the start of the module, then if the user adds a suitable point layer she/he will press a Refresh button to force addition to the list. Automatic refresh on popup has been tested as too time and resource consuming.


GRASS and R integration are really smooth with Python. In next days I'll build the minimum-option kriging procedure via automap, hopefully for next report.  
GRASS and R integration are really smooth with Python. In next days I'll build the minimum-option kriging procedure via automap, hopefully for next report.  
Line 39: Line 39:
** Initial '''sanity checks''':
** Initial '''sanity checks''':
*** Check '''dependencies''' all before interface load. If something is missing, a popup window will give a clear message.
*** Check '''dependencies''' all before interface load. If something is missing, a popup window will give a clear message.
*** '''Filter suitable point vector maps''' without creating lag before interface creation -> splash screen
*** Create a '''package''' to delegate dependencies check!
*** Create a '''package''' to delegate dependencies check!
** '''Variogram fit''': in a separated window, with allowed interactive model fit?
** '''Variogram fit''': in a separated window, with allowed interactive model fit?

Revision as of 15:31, 24 June 2009

v.autokrige: Python porting and wxPython GUI addition

Anne Ghisla's Google Summer of Code 2009 project, mentored by Martin Landa and Michael Barton

Aim of the project

As GRASS presently lacks kriging capability, it is performed via an add-on, v.autokrige, that delegates analysis to R (package automap). This module is written in ksh and has the classical autogenerated GUI. The project aims to rewrite the module in Python, creating a new GUI in wxPython that allows the user to refine parameters. I therefore plan to examinate present v.autokrige code and port it into Python, possibly improving it at the same time, and add a wxPython GUI. All time left will be dedicated to addition of further functionality, in respect of the most needed.

State of the art

Writing the wxPython interface, with GRASS bindings and part of R bindings - have a look at it from svn:

Attention: The module is now called v.krige, as its functions go beyond automatic kriging.

svn co https://svn.osgeo.org/grass/grass-addons/vector/v.krige v.krige

It is quite complete, allows the user to choose among automap (default choice), gstat (advanced functions) and geoR.

These three packages are not all required; if the user already uses R for kriging, the interface will load the package(s) present in R library. If none of the packages is available a popup message will ask for installation of at least one of the packages and abort. That is, crash gracefully and with clear explanation. I know what it means to crawl among missing dependencies and I don't want my users to be puzzled by inexplicable error messages.

Question: another layout can be a Choicebox with the three packages and in respect of the choice the drawing of corresponding widgets. Opinions welcome.

The lsit of point layers is provided by VectorSelect class (gselect module) to provide uniformity of layout with wxGUI. The filtering is done once at the start of the module, then if the user adds a suitable point layer she/he will press a Refresh button to force addition to the list. Automatic refresh on popup has been tested as too time and resource consuming.

GRASS and R integration are really smooth with Python. In next days I'll build the minimum-option kriging procedure via automap, hopefully for next report. Variogram autofit (package automap) is up and running!!

Bulletin Board

  • TODOs
    • Document dependencies - continue fill up html template page
    • Set CLI (SWIG?)
    • Refactor the code where widgets are created and set into sizers, like LayerBook's __createModifyPage in dbm.py
    • Start integrating gstat functions
  • IDEAS TO DISCUSS
    • Initial sanity checks:
      • Check dependencies all before interface load. If something is missing, a popup window will give a clear message.
      • Create a package to delegate dependencies check!
    • Variogram fit: in a separated window, with allowed interactive model fit?
  • BLOCKING ISSUES
    • Automap's autoKrige() internal function fails managing projections - wait for patch/workaround

Documentation

Kriging theory from Isaaks and Srivastava's "An Introduction to Applied Geostatistics" [0]. Great book.

Users feedback and other kriging software

  • Call for users has been done on gfoss.it, grass-it and grass-users, r-sig-geo and r-sig-ecology mailing lists. Thanks to Giovanni Manghi and João Tiago from GFOSS.pt and Giovanni Allegri for user advices, and to Edzer Pebesma, Paul Hiemstra, Ebrahim Jahanshiri, Dylan Beaudette for support on R side.
  • There are (have been) some GRASS modules that perform kriging, like:
    • v.surf.krige, v.variogram: no longer developed by authors, difficult to port.
    • s.surf.krig: deprecated since GRASS 6 series.
    • scripts including R and spgrass6: are the concrete basis upon which build the module.
  • Obtained ArcGIS availability for comparison of interface and functionality, because it is AFAIK the widespread tool used for kriging. Also Isats [1] is a valuable source of ispiration for interface design.

Planned Timeline

Draft:

  • create a wxPython interface for ordinary kriging. Will follow Humane Interface rules [2].
  • integrate R package automap
  • midterm deadline: have a working module that performs ordinary kriging.
  • integrate more R functions from gstat and geoR, giving the user the choice between gstat and geoR (kriging results vary in respect of implementation)

Weekly reports are available on SoC mailing list [3]

References