Going back to my question about about cutting out sections of maps. What I am looking for is being able to project a raster map of say temperature over the globe (180, -180, 90, -90), then I bring in a vector map of say New Mexico. My projection at this point has the globe temperature data and a small vector map of New Mexico. Is there a way to cut out a piece of that map leaving new mexico with the the temperature data over it, leaving behind all the rest of the global temperature data. The directions for v.digit, I do not think quite does this operation.
Aaron
On Wednesday 18 July 2007 12:24, goldneaa@onid.orst.edu wrote:
Going back to my question about about cutting out sections of maps.
What I am looking for is being able to project a raster map of say
temperature over the globe (180, -180, 90, -90), then I bring in a
vector map of say New Mexico. My projection at this point has the
globe temperature data and a small vector map of New Mexico. Is there
a way to cut out a piece of that map leaving new mexico with the the
temperature data over it, leaving behind all the rest of the global
temperature data. The directions for v.digit, I do not think quite
does this operation.
Aaron
Decide on a common projection, I would go with something like AEA or LCC for
the state of NW, then project your temp data into that location. You can cut
out the temp data by
1. convert NM outline into raster
2. rename raster to MASK
3. make a copy of the projected temp data
4. the new copy with be clipped by the raster MASK
dylan
goldneaa@onid.orst.edu wrote:
Going back to my question about about cutting out sections of maps.
What I am looking for is being able to project a raster map of say
temperature over the globe (180, -180, 90, -90), then I bring in a
vector map of say New Mexico. My projection at this point has the
globe temperature data and a small vector map of New Mexico. Is there
a way to cut out a piece of that map leaving new mexico with the the
temperature data over it, leaving behind all the rest of the global
temperature data. The directions for v.digit, I do not think quite
does this operation.
v.to.rast ... output=MASK
r.resample ...
g.remove rast=MASK
If a raster map named "MASK" exists, any module which reads raster
maps will automatically see nulls wherever the MASK cells are zero.
--
Glynn Clements <glynn@gclements.plus.com>
goldneaa@onid.orst.edu wrote:
Is there a way to cut out a piece of that map leaving new mexico
with the the temperature data over it, leaving behind all the rest
of the global temperature data. The directions for v.digit, I do
not think quite does this operation.
See the "v.cutter" thread from few hours ago.
Maciek