<br><br><div class="gmail_quote">2010/3/1 maven apache <span dir="ltr"><<a href="mailto:apachemaven0@gmail.com">apachemaven0@gmail.com</a>></span><br><blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
<br><br><div class="gmail_quote">2010/3/1 Hamish <span dir="ltr"><<a href="mailto:hamish_b@yahoo.com" target="_blank">hamish_b@yahoo.com</a>></span><div><div></div><div class="h5"><br><blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
<div><div></div><div>maven apache wrote:<br>
> When I use the r.fillnulls command to resample a map:<br>
> r.fillnulls input=aa output=ab tension=40 smooth=0.1<br>
> It show me that: taking too long to find points for<br>
> interpolation--please change the region to area where your<br>
> points are. Continuing calculations...<br>
><br>
> What does it mean?<br>
<br>
<br>
</div></div>It means you should read the fabulous "NOTES" section of the r.fillnulls<br>
man page. ;-)<br>
<br>
(the warning is harmless in this context, ignore it)<br>
<font color="#888888"><br>
<br>
Hamish<br>
</font><br>
<br>
ps::devs- use r.to.vect -zb in the r.fillnulls script? v.surf.rst doesn't<br>
need them and it could save a lot of memory etc if lidar data with many<br>
holes is thrown at it.<br>
<br>
<br>
<br></blockquote></div></div><div>Oh, thanks, I see it.<br>However what perplexing me these days is to call grass in command line model in windows platform,forward a message:<br><br><br>HI:<br>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.<br>The grass
variable bat is :<br>
-------<br>@echo off<br><br>rem *******Environment variables***********<br><br>rem
Set GRASS Installation Directory Variable<br>set
GRASSDIR=F:\anzhuang\GRASS<br><br><br>rem set GISBASE=D:\GrassDataBase<br><br>set
GISRC=C:\Documents and Settings\Administrator\.<div>grassrc6<br>
<br>rem Set Path to utilities (libraries and bynaries) used by GRASS<br><br>set
PATH=%GRASSDIR%\msys\bin;%PATH%<br>set PATH=%GRASSDIR%\extrabin;%GRASSDIR%\extralib;%GRASSDIR%\extrabin;%PATH%<br>set
PATH=%GRASSDIR%\tcl-tk\bin;%GRASSDIR%\sqlite\bin;%GRASSDIR%\gpsbabel;%PATH%<br>
set PATH=%PATH%;%GRASSDIR%\bin;%GRASSDIR%\lib;%GRASSDIR%\scripts;%GRASSDIR%\etc;%GRASSDIR%\sqlite\lib;%GRASSDIR%\sqlite\bin;<br><br><br>rem
Path to the proj files (notably the epsg projection list)<br>set
GRASS_PROJSHARE=%GRASSDIR%\proj<br>
----------<br><br>And I write a work bat :<br>---------<br>r.in.gdal
input=D:\data\hdf5\new\OMI-Aura_L3-OMTO3e_2009m0729_v003-2009m0731t020127.tif
output=aa location=al<br>g.gisenv set=LOCATION_NAME=al<br>g.region
rast=aa<br>
r.null map=aa setnull=-1.267e+30--1.268e+30<div class="im"><br>r.fillnulls input=aa
output=ab tension=40 smooth=0.1<br></div>r.out.gdal input=ab format=GTiff
output=D:/OMI-Aura_L3-OMTO3e_2009m0729_v003-2009m0731t020127.tif<br>---------<br>It
does not work,.<br>
I wonder why?<br><br>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:<br>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)<br>
2) change the subdataset to geotiff, use gdal_translate<br>3) import the
tiff to grass and create a new location according the tif<br>4) log in
to the new location use g.gisenv<br>5) set some sepical cell data to
null by using r.null....<br>
6) Interpolation the null data using r.fillnulls.<br>7) export the map
to geotiff.<br><br>So anyone can give some advise?</div> </div></div></blockquote><div>I am in hurry, anyone can help me? <br></div></div><br>