[GRASS-user] r.fillnulls error

GRASS GRASS user
1

HI:
When I use the r.fillnulls command to resample a map:
r.fillnulls input=aa output=ab tension=40 smooth=0.1
It show me that: taking too long to find points for interpolation–please change the region to area where your points are. Continuing calculations…
What does it mean?

2

On Mon, Mar 1, 2010 at 8:13 AM, maven apache <apachemaven0@gmail.com> wrote:

HI:
When I use the r.fillnulls command to resample a map:
r.fillnulls input=aa output=ab tension=40 smooth=0.1
It show me that: taking too long to find points for interpolation--please
change the region to area where your points are. Continuing calculations...
What does it mean?

It is a friendly warning to tell you that you *may* be out of area of interest.
But in the case of r.fillnulls it is likely that you have large NULL areas,
so you can ignore the message. Be sure to validate the result of course (as
usual).

Markus

3

maven apache wrote:

When I use the r.fillnulls command to resample a map:
r.fillnulls input=aa output=ab tension=40 smooth=0.1
It show me that: taking too long to find points for
interpolation--please change the region to area where your
points are. Continuing calculations...

What does it mean?

It means you should read the fabulous "NOTES" section of the r.fillnulls
man page. :wink:

(the warning is harmless in this context, ignore it)

Hamish

ps::devs- use r.to.vect -zb in the r.fillnulls script? v.surf.rst doesn't
need them and it could save a lot of memory etc if lidar data with many
holes is thrown at it.

4

2010/3/1 Hamish <hamish_b@yahoo.com>

maven apache wrote:

When I use the r.fillnulls command to resample a map:
r.fillnulls input=aa output=ab tension=40 smooth=0.1
It show me that: taking too long to find points for
interpolation–please change the region to area where your
points are. Continuing calculations…

What does it mean?

It means you should read the fabulous “NOTES” section of the r.fillnulls
man page. :wink:

(the warning is harmless in this context, ignore it)

Hamish

ps::devs- use r.to.vect -zb in the r.fillnulls script? v.surf.rst doesn’t
need them and it could save a lot of memory etc if lidar data with many
holes is thrown at it.

Oh, thanks, I see it.
However what perplexing me these days is to call grass in command line model in windows platform,forward a message:

HI:
I want to run grass in a command line model, however when I write a bat,it works just for 6.3 version , and some command can not work for example I can run g.version ,rather than r.fillnulls.
The grass variable bat is :

@echo off

rem Environment variables****

rem Set GRASS Installation Directory Variable
set GRASSDIR=F:\anzhuang\GRASS

rem set GISBASE=D:\GrassDataBase

set GISRC=C:\Documents and Settings\Administrator.
grassrc6

rem Set Path to utilities (libraries and bynaries) used by GRASS

set PATH=%GRASSDIR%\msys\bin;%PATH%
set PATH=%GRASSDIR%\extrabin;%GRASSDIR%\extralib;%GRASSDIR%\extrabin;%PATH%
set PATH=%GRASSDIR%\tcl-tk\bin;%GRASSDIR%\sqlite\bin;%GRASSDIR%\gpsbabel;%PATH%
set PATH=%PATH%;%GRASSDIR%\bin;%GRASSDIR%\lib;%GRASSDIR%\scripts;%GRASSDIR%\etc;%GRASSDIR%\sqlite\lib;%GRASSDIR%\sqlite\bin;

rem Path to the proj files (notably the epsg projection list)
set GRASS_PROJSHARE=%GRASSDIR%\proj

And I write a work bat :

r.in.gdal input=D:\data\hdf5\new\OMI-Aura_L3-OMTO3e_2009m0729_v003-2009m0731t020127.tif output=aa location=al
g.gisenv set=LOCATION_NAME=al
g.region rast=aa
r.null map=aa setnull=-1.267e+30–1.268e+30
r.fillnulls input=aa output=ab tension=40 smooth=0.1
r.out.gdal input=ab format=GTiff output=D:/OMI-Aura_L3-OMTO3e_2009m0729_v003-2009m0731t020127.tif

It does not work,.
I wonder why?

In fact what I am trying to do is to call grass in my web application(like the web processing service),now I want to handle the hdf5 data,the following is my idea:

  1. use gdalinfo to get the information of the he5,(the gdal model in grass 6.3 does not support this command, the 6.4 support)
  2. change the subdataset to geotiff, use gdal_translate
  3. import the tiff to grass and create a new location according the tif
  4. log in to the new location use g.gisenv
  5. set some sepical cell data to null by using r.null…
  6. Interpolation the null data using r.fillnulls.
  7. export the map to geotiff.

So anyone can give some advise?

5

Dear GRASS users,

I’m having some problems trying to find a way to snap/move all points from one layer to the nearest node in another layer.

I found that with v.distance, I can find distances to all nodes in the desired layer (the information is in the output window, although I haven’t figured out yet how to manipulate this).

Also, v.edit has a move and a snap option, although I believe the module expects some input like coordinates, not the name of another layer.

So, I’m guessing that the solution will probably involve working with these two funcions, but I haven’t figured out how

Perhaps I’m just looking in the wrong place and there is a function in GRASS that performs this “snap to another layer” in one go.

I’m sorry if this is kind of a beginner question…
Thanks for any help.
António

6

I've done it using v.distance. It can output the x and y coords of the
nearest feature (check upload option). So then I just created a new
point vector with the coords given by v.distance

Cheers
Daniel

2010/3/1 "António M. Rodrigues" <amrodrigues@fcsh.unl.pt>:

Dear GRASS users,
I'm having some problems trying to find a way to snap/move all points from
one layer to the nearest node in another layer.
I found that with v.distance, I can find distances to all nodes in the
desired layer (the information is in the output window, although I haven't
figured out yet how to manipulate this).
Also, v.edit has a move and a snap option, although I believe the module
expects some input like coordinates, not the name of another layer.
So, I'm guessing that the solution will probably involve working with these
two funcions, but I haven't figured out how
Perhaps I'm just looking in the wrong place and there is a function in GRASS
that performs this "snap to another layer" in one go.
I'm sorry if this is kind of a beginner question...
Thanks for any help.
António

_______________________________________________
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

7

Yes, that works.

thanks.

On Mar 1, 2010, at 11:59 AM, Daniel Victoria wrote:

I've done it using v.distance. It can output the x and y coords of the
nearest feature (check upload option). So then I just created a new
point vector with the coords given by v.distance

Cheers
Daniel

2010/3/1 "António M. Rodrigues" <amrodrigues@fcsh.unl.pt>:

Dear GRASS users,
I'm having some problems trying to find a way to snap/move all points from
one layer to the nearest node in another layer.
I found that with v.distance, I can find distances to all nodes in the
desired layer (the information is in the output window, although I haven't
figured out yet how to manipulate this).
Also, v.edit has a move and a snap option, although I believe the module
expects some input like coordinates, not the name of another layer.
So, I'm guessing that the solution will probably involve working with these
two funcions, but I haven't figured out how
Perhaps I'm just looking in the wrong place and there is a function in GRASS
that performs this "snap to another layer" in one go.
I'm sorry if this is kind of a beginner question...
Thanks for any help.
António

_______________________________________________
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

8

2010/3/1 maven apache <apachemaven0@gmail.com>

2010/3/1 Hamish <hamish_b@yahoo.com>

maven apache wrote:

When I use the r.fillnulls command to resample a map:
r.fillnulls input=aa output=ab tension=40 smooth=0.1
It show me that: taking too long to find points for
interpolation–please change the region to area where your
points are. Continuing calculations…

What does it mean?

It means you should read the fabulous “NOTES” section of the r.fillnulls
man page. :wink:

(the warning is harmless in this context, ignore it)

Hamish

ps::devs- use r.to.vect -zb in the r.fillnulls script? v.surf.rst doesn’t
need them and it could save a lot of memory etc if lidar data with many
holes is thrown at it.

Oh, thanks, I see it.
However what perplexing me these days is to call grass in command line model in windows platform,forward a message:

HI:
I want to run grass in a command line model, however when I write a bat,it works just for 6.3 version , and some command can not work for example I can run g.version ,rather than r.fillnulls.
The grass variable bat is :

@echo off

rem Environment variables****

rem Set GRASS Installation Directory Variable
set GRASSDIR=F:\anzhuang\GRASS

rem set GISBASE=D:\GrassDataBase

set GISRC=C:\Documents and Settings\Administrator.
grassrc6

rem Set Path to utilities (libraries and bynaries) used by GRASS

set PATH=%GRASSDIR%\msys\bin;%PATH%
set PATH=%GRASSDIR%\extrabin;%GRASSDIR%\extralib;%GRASSDIR%\extrabin;%PATH%
set PATH=%GRASSDIR%\tcl-tk\bin;%GRASSDIR%\sqlite\bin;%GRASSDIR%\gpsbabel;%PATH%
set PATH=%PATH%;%GRASSDIR%\bin;%GRASSDIR%\lib;%GRASSDIR%\scripts;%GRASSDIR%\etc;%GRASSDIR%\sqlite\lib;%GRASSDIR%\sqlite\bin;

rem Path to the proj files (notably the epsg projection list)
set GRASS_PROJSHARE=%GRASSDIR%\proj

And I write a work bat :

r.in.gdal input=D:\data\hdf5\new\OMI-Aura_L3-OMTO3e_2009m0729_v003-2009m0731t020127.tif output=aa location=al
g.gisenv set=LOCATION_NAME=al
g.region rast=aa
r.null map=aa setnull=-1.267e+30–1.268e+30

r.fillnulls input=aa output=ab tension=40 smooth=0.1

r.out.gdal input=ab format=GTiff output=D:/OMI-Aura_L3-OMTO3e_2009m0729_v003-2009m0731t020127.tif

It does not work,.
I wonder why?

In fact what I am trying to do is to call grass in my web application(like the web processing service),now I want to handle the hdf5 data,the following is my idea:

  1. use gdalinfo to get the information of the he5,(the gdal model in grass 6.3 does not support this command, the 6.4 support)
  2. change the subdataset to geotiff, use gdal_translate
  3. import the tiff to grass and create a new location according the tif
  4. log in to the new location use g.gisenv
  5. set some sepical cell data to null by using r.null…
  6. Interpolation the null data using r.fillnulls.
  7. export the map to geotiff.

So anyone can give some advise?

I am in hurry, anyone can help me?

9

Dear users,

I'm using r.cost to calculate time distances based on friction calculated using a road network and some other information.

From what I understant, r.cost only outputs "accessbility" surfaces. The problem is I also need the identity of the nearest target (catchment areas or allocation surfaces or some other way you whish to call it).

How could I achieve this?

Thanks in advance.
António