Hi,
I might have found a solution to my problem, which I shared with you earlier today. I am writing about it in case someone else has a similar issue and this may be a way of solving the problem. However, I am not certain if this is the easiest or the most practical way of getting the results.
The problem was the presence of few negative value cells in ‘cubic’ reprojected DEMs (from Latlong to UTM). I ran r.mapcalc to make these negative cells NULLs. Then, I ran ‘r.fillnulls’ in that part of the DEM and patched the ‘filled’ map with the original map (the one that has NULLs). Finally, I interpolated the resultant (i.e., patched) imagery from 30m to 10m. The values in 10m DEM look accurate and free of negative value cells.
If you know another way of doing it and share with me, I will greatly appreciate.
–
BÜLENT ARIKAN
School of Human Evolution and Social Change
Arizona State University
Tempe - AZ
85287-2402