I tried v.surf.contour and after about 24 hrs of computation the progress bar was still nowhere to be seen. I'll try to acquire the points from the county. But i converted the contours to points and tried the whole series of v.surf (rst, bspline, idw ...etc) with the same result as i stated previously. But with 50 ft resolution I can complete the computation. Here is the output for g.region to answer your question about nrows and ncolumns
north: 808710.758065 south: 736239.936599 west: 377975.676116 east: 465460.787903 nsres: 10.00149344 ewres: 9.99829849 rows: 7246 cols: 8750 cells: 63402500 Cheers On Mon, Sep 20, 2010 at 10:37 AM, Markus Metz < markus.metz.gisw...@googlemail.com> wrote: > kapo coulibaly wrote: > > I'm trying to interpolate a raster from a shapefile (points or contours). > > The shapefile is a 2ft contour lidar-derived elevation and I need a 10 ft > > resolution DEM. > > Extracting contours discards a lot of information, results are more > detailed if a raster surface is interpolated straight from the LiDAR > points. Try any of the v.surf.* modules with the LiDAR points, not the > contours. > > > Obviously even on a dual core dual processor 64 bit > > workstation with 8 GB of ram I couldn't complete the computation. > > How many cells (rows, columns) are in the current region? Did you use > r.surf.contour? The -s flag might help to reduce memory requirements a > bit for very large datasets. > > Markus M > > > Is there > > an automated way to slice the domain in smaller pieces, carry out the > > computation and piece them back together afterward? I'm thinking of > writing > > a script but i was wondering if it has already been done. > > > > Thanks > > > > _______________________________________________ > > grass-user mailing list > > grass-user@lists.osgeo.org > > http://lists.osgeo.org/mailman/listinfo/grass-user > > > > >
_______________________________________________ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user