GRASS and Python
Python extensions for GRASS GIS
Python-SWIG-GRASS interface
There is a prototype GRASS-SWIG interface available (thanks to Sajith VK), find it in GRASS 6-CVS: swig/python/
Background: SWIG (Simplified Wrapper and Interface Generator) is:
- A compiler that turns ANSI C/C++ declarations into scripting language interfaces.
- Completely automated (produces a fully working Python extension module).
- Language neutral. SWIG can also target Tcl, Perl, Guile, MATLAB, etc...
- Attempts to eliminate the tedium of writing extension modules.
TODO: Implement modules support in a Python class using --interface-description and a Python-XML parser. This should be a generic class with module's name as parameter, returning back an object which describes the module (description, flags, parameters, status of not/required). See GRASS 6 wxPython interface for inspiration. Important is to auto-generate the GRASS-Python class at compile time with a Python script.
Python-GRASS add-ons
Stand-alone addons:
- Jáchym Čepický's G-ps.map, a GUI to typeset printable maps with ps.map (http://les-ejk.cz/?cat=gpsmap)
- Jáchym Čepický's v.pydigit, a GUI to v.edit (http://les-ejk.cz/?cat=vpydigit)
- Jáchym Čepický's PyWPS, GRASS-Web Processing Service (http://les-ejk.cz/?cat=pywps)
Using Grass gui.tcl in python
Here is some example code to use the grass automatically generated guis in python code. This could (should) all be bundled up and abstracted away so that the implementation can be replaced later.
import Tkinter import os # Startup (once): tk = Tkinter.Tk() tk.eval ("wm withdraw .") tk.eval ("source $env(GISBASE)/etc/gui.tcl") # Here you could do various things to change what the gui does # See gui.tcl and README.GUI # Make a gui (per dialog) # This sets up a window for the command. # This can be different to integrate with tkinter: tk.eval ('set path ".dialog$dlg"') tk.eval ('toplevel .dialog$dlg') # Load the code for this command: fd = os.popen ("d.vect --tcltk") gui = fd.read() # Run it tk.eval(gui) dlg = tk.eval('set dlg') # This is used later to get and set # Get the current command in the gui we just made: currentcommand = tk.eval ("dialog_get_command $dlg") # Set the command in the dialog we just made: tk.eval ("dialog_set_command $dlg 'd.vect map=roads'")
Links
- GRASS-SWIG interface http://mpa.itc.it/markus/grass61progman/swig/
- SWIG http://www.swig.org/
- SIP (C/C++ bindings generator) http://directory.fsf.org/all/Python-SIP.html
- Python and GDAL/OGR:
- Python bindings to PROJ:
- Other external projects
- Scientific Python: http://www.scipy.org/