|
@@ -1,10 +1,15 @@
|
|
|
#!/usr/bin/python
|
|
|
|
|
|
# run within GRASS Spearfish session
|
|
|
+# run this before starting python to append module search path:
|
|
|
+# export PYTHONPATH=/usr/src/grass70/swig/python
|
|
|
+# check with "import sys; sys.path"
|
|
|
+# or:
|
|
|
+# sys.path.append("/usr/src/grass70/swig/python")
|
|
|
+# FIXME: install the grass bindings in $GISBASE/lib/ ?
|
|
|
|
|
|
import os, sys
|
|
|
-import swig.grass as grasslib
|
|
|
-import swig.vector as grassvect
|
|
|
+import grass
|
|
|
|
|
|
if not os.environ.has_key("GISBASE"):
|
|
|
print "You must be in GRASS GIS to run this program."
|
|
@@ -16,35 +21,35 @@ else:
|
|
|
input = raw_input("Vector Map Name? ")
|
|
|
|
|
|
# initialize
|
|
|
-grasslib.G_gisinit('')
|
|
|
+grass.G_gisinit('')
|
|
|
|
|
|
# find map in search path
|
|
|
-mapset = grasslib.G_find_vector2(input,'')
|
|
|
+mapset = grass.G_find_vector2(input,'')
|
|
|
|
|
|
# define map structure
|
|
|
-map = grassvect.Map_info()
|
|
|
+map = grass.Map_info()
|
|
|
|
|
|
# define open level (level 2: topology)
|
|
|
-grassvect.Vect_set_open_level (2)
|
|
|
+grass.Vect_set_open_level (2)
|
|
|
|
|
|
# open existing map
|
|
|
-grassvect.Vect_open_old(map, input, mapset)
|
|
|
+grass.Vect_open_old(map, input, mapset)
|
|
|
|
|
|
# query
|
|
|
print 'Vect map: ', input
|
|
|
-print 'Vect is 3D: ', grassvect.Vect_is_3d (map)
|
|
|
-print 'Vect DB links: ', grassvect.Vect_get_num_dblinks(map)
|
|
|
-print 'Map Scale: 1:', grassvect.Vect_get_scale(map)
|
|
|
+print 'Vect is 3D: ', grass.Vect_is_3d (map)
|
|
|
+print 'Vect DB links: ', grass.Vect_get_num_dblinks(map)
|
|
|
+print 'Map Scale: 1:', grass.Vect_get_scale(map)
|
|
|
# misleading:
|
|
|
-# print 'Number of lines:', grassvect.Vect_get_num_lines(map)
|
|
|
+# print 'Number of lines:', grass.Vect_get_num_lines(map)
|
|
|
# how to access GV_POINT?
|
|
|
-# print 'Number of points: ', grassvect.Vect_get_num_primitives(map,GV_POINT)
|
|
|
+# print 'Number of points: ', grass.Vect_get_num_primitives(map,GV_POINT)
|
|
|
# confusing:
|
|
|
-#print 'Number of lines: ', Vect_get_num_primitives(map,GV_LINE)
|
|
|
-#print 'Number of areas:', Vect_get_num_primitives(map,GV_AREA)
|
|
|
-print 'Number of areas:', grassvect.Vect_get_num_areas(map)
|
|
|
+#print 'Number of lines: ', grass.Vect_get_num_primitives(map,GV_LINE)
|
|
|
+#print 'Number of areas:', grass.Vect_get_num_primitives(map,GV_AREA)
|
|
|
+print 'Number of areas:', grass.Vect_get_num_areas(map)
|
|
|
|
|
|
# close map
|
|
|
-grassvect.Vect_close(map)
|
|
|
+grass.Vect_close(map)
|
|
|
## end of the python script
|
|
|
|