[GRASS-user] High resolution raster interpolation

[kapo coulibaly]

I tried r.surf.contour but using the vector contours converted to raster, it
didn't work either. I'll try the actual lidar point and get back to you.
Maybe the contours leave wide empty spaces (it is a pretty flat terrain) and
it makes computation more difficult.

Thanks.

On Mon, Sep 20, 2010 at 11:19 AM, Markus Metz <
markus.metz.giswork at googlemail.com> wrote:

> kapo coulibaly wrote:
> > I tried v.surf.contour
>
> Please try r.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
> >
> ...
> > rows:       7246
> > cols:       8750
> > cells:      63402500
> >
> Weird, that's only 63 million cells, shouldn't cause memory troubles.
> With 8GB RAM, interpolations should be no problem (I tested with
> several 100 million cells), but depending on the settings, can take
> some time.
>
> Markus M
>
> >
> > Markus Metz 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 at lists.osgeo.org
> >> > http://lists.osgeo.org/mailman/listinfo/grass-user
> >> >
> >> >
> >
> >
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://lists.osgeo.org/pipermail/grass-user/attachments/20100920/80a56561/attachment.html