Rich wrote: > The approach that worked for me was described by > Helmut Kudrnovsky: > > 1) zoom to your area of interest of your larger > region. > > 2) v.in.region e.g. "v.in.region output=testarea", > so the extent of your > area of interest is defined by the vector > "testarea" > > 3) "g.region vect=testarea", by this the > region is set to your area of > interest for analysis. "This region defines the > geographic area in which > all GRASS displays and raster analyses will be > done." > > Works like a charm.
the v.in.region and g.region vect= steps are redundant, but make sure when you zoom you update the computational region, not just the display one in the GUI (d.zoom + d.mon does this automatically) Hamish _______________________________________________ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user