#2671: r.surf.idw won't make interpolation
---------------------------------------+-------------------------
Reporter: kuszinger | Owner: grass-dev@…
Type: defect | Status: new
Priority: normal | Milestone: 7.1.0
Component: Raster | Version: svn-trunk
Keywords: interpolation, e.surf.idw | CPU: x86-64
Platform: Linux |
---------------------------------------+-------------------------
However '''r.surf.idw''' won't make the interpolation for my dataset, I
tested on NC sample dataset where it was ok. ??
I was not able to reproduce it with any other dataset but mine. I've
created a cut of my data and accessible in the next comment.
Projection: UTM 48N / WGS84 (on the area of Cambodia)
r.surf.idw runs for the whole coverage. No error message displayed. When
finished, the target raster map is the very same as the input was.
Command run:
''r.surf.idw --overwrite input=debugfrom@lidar output=debugto npoints=7
''
Purpose: generate a continuous surface from a very sparse set of points on
the map.
#2671: r.surf.idw won't make interpolation
------------------------+---------------------------------------
Reporter: kuszinger | Owner: grass-dev@…
Type: defect | Status: new
Priority: normal | Milestone: 7.1.0
Component: Raster | Version: svn-trunk
Resolution: | Keywords: interpolation, e.surf.idw
CPU: x86-64 | Platform: Linux
------------------------+---------------------------------------
Comment (by kuszinger):
Test data is accessible here: [http://www.giscom.hu/debugfrom.pack\].
4.4 MB and I tested r.surf.idw against it with the command above. Output
was the same as input.
#2671: r.surf.idw won't make interpolation
------------------------+---------------------------------------
Reporter: kuszinger | Owner: grass-dev@…
Type: defect | Status: new
Priority: normal | Milestone: 7.1.0
Component: Raster | Version: svn-trunk
Resolution: | Keywords: interpolation, e.surf.idw
CPU: x86-64 | Platform: Linux
------------------------+---------------------------------------
Comment (by annakrat):
Replying to [comment:1 kuszinger]:
>
> Test data is accessible here: [http://www.giscom.hu/debugfrom.pack\].
> 4.4 MB and I tested r.surf.idw against it with the command above. Output
was the same as input.
>
>
Confirmed. Does it work in 6.4? What example did you use with NC dataset?
I was not able to make it work.
#2671: r.surf.idw won't make interpolation
------------------------+---------------------------------------
Reporter: kuszinger | Owner: grass-dev@…
Type: defect | Status: new
Priority: normal | Milestone: 7.1.0
Component: Raster | Version: svn-trunk
Resolution: | Keywords: interpolation, r.surf.idw
CPU: x86-64 | Platform: Linux
------------------------+---------------------------------------
Comment (by kuszinger):
'''NC''': used '''elevation''', the whole, set g.region res=5 so the
target dataset is bigger, it ran well with good results. Anyway, the
source dataset is not sparse. I''' rerun on the NC elevation and making it
sparse with random nulls, ok?'''
'''6.4: yes, worked like charm''', less than 2GB memory usage (same
machine same everything).
One more extra: while the above debug region is:
#2671: r.surf.idw won't make interpolation
------------------------+---------------------------------------
Reporter: kuszinger | Owner: grass-dev@…
Type: defect | Status: new
Priority: normal | Milestone: 7.1.0
Component: Raster | Version: svn-trunk
Resolution: | Keywords: interpolation, r.surf.idw
CPU: x86-64 | Platform: Linux
------------------------+---------------------------------------
Comment (by neteler):
The indicated test is now performed as follows:
{{{
# GRASS 7.1.svn (nc_spm_08_grass7)
g.region raster=elevation@PERMANENT -p
r.random -z --overwrite input=elevation@PERMANENT npoints=100000
raster=elev_sparse
r.surf.idw --overwrite input=elev_sparse output=elev_sparse_idw npoints=7
ERROR: This module currently only works for integer (CELL) maps
}}}
So, no result should be expected unless #2672 is solved.
Replying to [comment:6 annakrat]:
> I think I fixed it in r65214. Null values were not identified correctly.
I still don't understand why 64 version actually works.
Backported in r65261.
> By the way, r.surf.idw works only with integer rasters. But this is for
a different ticket.