Vector Overlapping Areas

From GRASS-Wiki
Revision as of 16:07, 7 July 2008 by ⚠️Rkrug (talk | contribs)
Jump to navigation Jump to search
The printable version is no longer supported and may have rendering errors. Please update your browser bookmarks and please use the default browser print function instead.

Q: How can I determine the number of overlapping features per location after import of SHAPE file?

A: You can run this (example) - note that it won't work with the DBF driver but requires SQLite, PostgreSQL, MySQL, or... driver:

  v.in.ogr dsn=./shp/ layer=fireHistory output=fireHistory
  v.category input=fireHistory output=tmpFire type=centroid
  v.db.addtable map=tmpFire table=tmpFire_2 layer=2 columns='cat int'
   
  # Export it to a raster layer:
  v.to.rast input=tmpFire layer=2 output=mf use=attr column=num
  # another solution (tested under SQLite)
  v.overlay ainput=tmpFire atype=area alayer=1 binput=tmpFire btype=area blayer=2 output=mf operator=or olayer=1,1,1 --overwrite
  v.db.update map=mf layer=1 column=b_cat value=1 where='b_cat isnull'