<div dir="ltr">I have tested my solution and it works. <div>Can I do the push on SVN? My question is for Martin Landa who is mentor for the authors of the command and did the last commit of nnbathy.py command?<div><div><br></div><div><div><br></div><div>I found also two more problem in the v.surf.nnbathy command. </div><div><br></div><div>1) The command crash. It is probaby due to the v.out.ascii -r command (the -r option do not works correctly, at least on GRASS 7.4). It is possible to remove the -r flag in this command and the v.in.nnbathy works correctly, even if the output is different. Perhaps is a better output, because it use also the points outside the computational region for the interpolation.</div><div><br></div><div>2) An other problem is due to the command to write output ascii file for interpolation --> fout.write(parts[0]+' '+parts[1]+' '+parts[3]) is valid only for 2D geometry while for 3D geometry the correct command is fout.write(parts[0]+' '+parts[1]+' '+parts[4])</div><div><br></div><div>I think we need an if condition to check if the geometry is 2 or 3D. I have no idea of a possible python command.. I will do some search.. Have you any suggestion?</div><div><br></div><div>R</div><div> </div></div></div></div></div><div class="gmail_extra"><br><div class="gmail_quote">2018-03-23 18:05 GMT+01:00 Roberto Marzocchi <span dir="ltr"><<a href="mailto:roberto.marzocchi@gmail.com" target="_blank">roberto.marzocchi@gmail.com</a>></span>:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Dear all, <div><br></div><div>I am using r.surf.nnbathy (on GRASS 7.4)</div><div><br></div><div>Using a projected location (meters) I obtain a result</div><div><br></div><div>Using lat/lon location I see the following error: </div><div><div><br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><br>Traceback (most recent call last):<br> File<br>"/home/roberto/.grass7/addons/<wbr>scripts/r.surf.nnbathy", line<br>80, in <module><br> main()<br> File<br>"/home/roberto/.grass7/addons/<wbr>scripts/r.surf.nnbathy", line<br>74, in main<br> obj = Nnbathy_raster(options)<br> File "/home/roberto/.grass7/addons/<wbr>etc/v.surf.nnbathy/nnba<br>thy.py", line 107, in __init__<br> Nnbathy.__init__(self, options)<br> File "/home/roberto/.grass7/addons/<wbr>etc/v.surf.nnbathy/nnba<br>thy.py", line 13, in __init__<br> self.region()<br> File "/home/roberto/.grass7/addons/<wbr>etc/v.surf.nnbathy/nnba<br>thy.py", line 20, in region<br> reg_N = float(kv['north'])<br>ValueError: invalid literal for float(): 47:00:29.111155N</blockquote><div><br></div><div><br></div><div><br></div><div>I think the error is due to the <a href="https://trac.osgeo.org/grass/browser/grass-addons/grass7/vector/v.surf.nnbathy/nnbathy.py" target="_blank">https://trac.osgeo.org/<wbr>grass/browser/grass-addons/<wbr>grass7/vector/v.surf.nnbathy/<wbr>nnbathy.py</a> file and in particular to the region() function</div><div><br></div><div><div><div><span style="white-space:pre-wrap"> </span> </div></div></div></div><blockquote style="margin:0 0 0 40px;border:none;padding:0px"><div><div><div><div>reg = grass.read_command("g.region", flags='p')</div></div></div></div><div><div><div><div> kv = grass.parse_key_val(reg, sep=':')</div></div></div></div><div><div><div><div>reg_N = float(kv['north'])</div></div></div></div><div><div><div><div>reg_W = float(kv['west'])</div></div></div></div><div><div><div><div>reg_S = float(kv['south'])</div></div></div></div><div><div><div><div>reg_E = float(kv['east'])</div></div></div></div><div><div><div><div>nsres = float(kv['nsres'])</div></div></div></div><div><div><div><div> ewres = float(kv['ewres'])</div></div></div></div></blockquote><div><div><div><br></div><div>perhaps this can be a solution:</div><div><br></div></div></div><blockquote style="margin:0 0 0 40px;border:none;padding:0px"><div><div><div>in_region = gscript.region()</div></div></div><div><div><div>reg_E=float(in_region['e'])</div></div></div><div><div><div>reg_W=float(in_region['w'])</div></div></div><div><div><div>reg_N=float(in_region['s'])</div></div></div><div><div><div>reg_S=float(in_region['s'])</div></div></div><div><div>nsres = float(in_region['nsres'])</div></div><div><div>ewres = float(in_region['ewres'])</div></div></blockquote><div><div><br></div><div>If you agree I will do some test and I can update the add-on</div><div><br></div><div>Any suggestions? </div><div><br></div><div>Thanks in advanced for your help</div><span class="HOEnZb"><font color="#888888"><div>R <br></div></font></span></div></div>
</blockquote></div><br></div>