[GRASS-user] High resolution raster interpolation

[Markus Metz]

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
>> >
>> >
>
>