R-sig-geo'ers:

Continuing on my quest to *precisely* reproject/resize a categorical raster
where the output cell is the category which has the maximum areal cover of
that cell (this is not the same as nearest neighbor), it is looking more
and more like the solution is to first convert a class raster to a polygon
surface, reproject/resize, then re-rasterize the output.

As such, is there a function within R that does this?  I see raster has the
"rasterToPolygons" but as far as I can tell, this produces one polygon per
cell (e.g. it doesn't merge the cells if they contain the same values --
Robert, correct me if I'm wrong here) -- this could cause some memory
issues, I imagine, on some big files.  I'm looking for something like
what gdal_polygonize.py
does.

Any thoughts?

--j

-- 
Jonathan A. Greenberg, PhD
Assistant Professor
Global Environmental Analysis and Remote Sensing (GEARS) Laboratory
Department of Geography and Geographic Information Science
University of Illinois at Urbana-Champaign
607 South Mathews Avenue, MC 150
Urbana, IL 61801
Phone: 217-300-1924
http://www.geog.illinois.edu/~jgrn/
AIM: jgrn307, MSN: jgrn307@hotmail.com, Gchat: jgrn307, Skype: jgrn3007

	[[alternative HTML version deleted]]

