Hamish wrote:
> g.region vect=NZMG_s27
> v.in.region output=cutter type=line --overwrite
> v.overlay ainput=wellington_contours atype=line \
> binput=cutter output=clipped operator=and
and I think you'll want v.in.region to use type=area, not line, otherwise
you risk only to get contour lines which cross the boundary and not those
that fall within it.
The manual [1] for v.overlay states that btype must be area. Shouldn't
v.overlay complain if you try to use a line type?
Craig
[1] http://grass.itc.it/grass64/manuals/html64_user/v.overlay.html
Well after about 27 hours on an imac it work ed, the command for v.overlay was
v.overlay ainput=S28 atype=line binput=cutter output=S28_Clipped operator=and
Unfortunatly I have about few hundred to do! If I can't make it any faster, might have to look at using the amazon compute cloud to do this in a reasonable time! Don't really want to do that though!
Cheers
Matthew
On 14/09/2008, at 6:08 AM, Craig Leat wrote:
Hamish wrote:
g.region vect=NZMG_s27
v.in.region output=cutter type=line --overwrite
v.overlay ainput=wellington_contours atype=line \
binput=cutter output=clipped operator=and
and I think you'll want v.in.region to use type=area, not line, otherwise
you risk only to get contour lines which cross the boundary and not those
that fall within it.
The manual [1] for v.overlay states that btype must be area. Shouldn't
v.overlay complain if you try to use a line type?
Craig
[1] http://grass.itc.it/grass64/manuals/html64_user/v.overlay.html
_______________________________________________
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user