Vector Overlapping Areas: Difference between revisions

From GRASS-Wiki
Jump to navigation Jump to search
No edit summary
No edit summary
Line 9: Line 9:
   # Export it to a raster layer:
   # Export it to a raster layer:
   v.to.rast input=tmpFire layer=2 output=mf use=attr column=num
   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'


[[Category: FAQ]]
[[Category: FAQ]]

Revision as of 16:07, 7 July 2008

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'