Hi, > use r.to.vect to map every cell to a polygon, project it to the lambert > map, do the calculation on the polygon, reproject it back to the you dont need to project it to lambert:
just run v.to.db -p option=area map=... then you get the size for every lines cell. To get the raster-map: r.mapcalc "ns_col=if(col()==1 , row() , null())" Achim _______________________________________________ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user