On Wed, Apr 02, 2008 at 12:51:35PM -0700, we recorded a bogon-computron collision of the <[EMAIL PROTECTED]> flavor, containing: > > > Tom Russo wrote: > > > > ... > > v.to.rast input=o48112d2_bounds output=o48112d2_bounds use=val value=1 > > (this one generates a raster that's 1 inside the neatline, null > > outside) > > echo > > "o48112d2_withoutcollar=if(!isnull(o48112d2_bounds),o48112d2_withcollar,null())" > > | r.mapcalc > > (this populates the _withoutcollar map only where the mask is > > non-null) > > (you could also do it with > > r.mask input=o48112d2_bounds > > echo "o48112d2_withoutcollar=o48112d2_withcollar"| r.mapcalc > > r.mask -r > > ) > > ... > > > > Hi > > Another way is to use gdal_rasterize. See point e) under the Mosaicing > section of this tutorial by Frank Warmerdam: > > http://trac.osgeo.org/gdal/wiki/UserDocs/RasterProcTutorial
Neat. The same 24kgrid shapefile could be the source of the vector clipping region for that technique, too. Thanks for pointing that out, I hardly ever use gdal_rasterize and didn't know it could do that. -- Tom Russo KM5VY SAR502 DM64ux http://www.swcp.com/~russo/ Tijeras, NM QRPL#1592 K2#398 SOC#236 AHTB#1 http://kevan.org/brain.cgi?DDTNM "And, isn't sanity really just a one-trick pony anyway? I mean all you get is one trick, rational thinking, but when you're good and crazy, oooh, oooh, oooh, the sky is the limit!" --- The Tick _______________________________________________ grass-user mailing list [email protected] http://lists.osgeo.org/mailman/listinfo/grass-user
