# this is a script to run the notebook with a grass environment available on the command line
# so it is possible to run grass commands adding the "!" prefix
if [ -n "$LD_LIBRARY_PATH" ] ; then
LD_LIBRARY_PATH="/usr/local/grass-7.1.svn/lib:$LD_LIBRARY_PATH"
else
LD_LIBRARY_PATH="/usr/local/grass-7.1.svn/lib"
fi
if [ -n "$PYTHONPATH" ] ; then
PYTHONPATH="/usr/local/grass-7.1.svn/etc/python:$PYTHONPATH"
else
PYTHONPATH="/usr/local/grass-7.1.svn/etc/python"
fi
GISBASE="/usr/local/grass-7.1.svn/"
PATH="$GISBASE/bin:$GISBASE/scripts:$PATH"
GIS_LOCK="$$"
GISRC="/home/$USER/.grass7/rc"
export LD_LIBRARY_PATH PYTHONPATH GISBASE PATH GIS_LOCK GISRC
!g.mapset location=loc_ncarolina_spm_base0.3.1 mapset=user1
WARNING: <user1> is already the current mapset
!g.list vect
boundary_region boundary_state census firestations geology geonames history_markers hospitals points_of_interest railroads roadsmajor schools streams streets zipcodes
!g.list rast
basin_50K basins elevation@PERMANENT elevation@user1 elevation_shade geology lakes landuse
!g.region rast=basin_50K -ap
projection: 99 (Lambert Conformal Conic) zone: 0 datum: nad83 ellipsoid: a=6378137 es=0.006694380022900787 north: 228500 south: 215000 west: 630000 east: 645000 nsres: 10 ewres: 10 rows: 1350 cols: 1500 cells: 2025000
from IPython.core.display import Image
import os
import sys
from grass.script import core as grass
from grass.script import gisenv
from grass.pygrass.modules.shortcuts import display as d
from grass.pygrass.modules.shortcuts import general as g
!rm -rf basin_50K.png
os.environ['GRASS_RENDER_IMMEDIATE'] = 'cairo'
os.environ['GRASS_RENDER_FILE'] = 'basin_50K.png'
os.environ['GRASS_RENDER_FILE_COMPRESSION'] = '9'
os.environ['GRASS_RENDER_WIDTH'] = '640'
os.environ['GRASS_RENDER_HEIGHT'] = '480'
os.environ['GRASS_RENDER_TRANSPARENT']='TRUE'
os.environ['GRASS_RENDER_READ_FILE']='TRUE'
os.environ['GRASS_RENDER_TRUECOLOR']='TRUE'
os.environ['GRASS_RENDER_PNG_AUTO_WRITE']='TRUE'
monitor_old = None
genv = gisenv()
if 'MONITOR' in genv:
monitor_old = genv['MONITOR']
g.gisenv(unset='MONITOR')
d.rast(map='basin_50K')
d.vect(map='streams')
#d.vect(map='p')
#d.redraw()
Image('basin_50K.png')