Hi List
I have a raster and a vector layer but would want to limit the
analysis to the bounding box of the vector layer - how do I accomplish
this.
Is there a vector module to create a vector bounding box (v.*).
Thanks
You should be able to set the region to the vector layer and then your
set. You plan to do the analysis on the raster, correct?
John
On 2/23/09, Intengu Technologies <sindile.bidla@gmail.com> wrote:
Hi List
I have a raster and a vector layer but would want to limit the
analysis to the bounding box of the vector layer - how do I accomplish
this.
Is there a vector module to create a vector bounding box (v.*).
Thanks
--
--
John C. Tull, Ph.D.
Conservation Director
Nevada Wilderness Project
8550 White Fir Street
Reno, NV 89523
775.746.7851 (office)
775.224.2947 (mobile)
www.wildnevada.org
Yes analysis on the raster - how do I do that John as I have a region
defined already, which is larger than the extent of the vector.
On 2/24/09, John C. Tull <john.tull@wildnevada.org> wrote:
You should be able to set the region to the vector layer and then your
set. You plan to do the analysis on the raster, correct?
John
On 2/23/09, Intengu Technologies <sindile.bidla@gmail.com> wrote:
Hi List
I have a raster and a vector layer but would want to limit the
analysis to the bounding box of the vector layer - how do I accomplish
this.
Is there a vector module to create a vector bounding box (v.*).
Thanks
--
--
John C. Tull, Ph.D.
Conservation Director
Nevada Wilderness Project
8550 White Fir Street
Reno, NV 89523
775.746.7851 (office)
775.224.2947 (mobile)
www.wildnevada.org
Yes analysis on the raster - how do I do that John as I have a region
defined already, which is larger than the extent of the vector.
# set to raster map to obtain raster's resolution
g.region rast=raster_map
# set bounds to match vector map and expand bounds to current resolution.
g.region vect=vector_map -ap
I am not sure if -a works properly if you do not also use res=. So maybe
set res= by hand in the second step to match what r.info says, and make sure
it has aligned cleanly.
maybe this way:
g.region vect=vector_map
g.region align=raster_map