<!DOCTYPE html PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN">
<html>
<head>
  <DEFANGED_meta content="text/html;charset=ISO-8859-1" http-equiv="Content-Type">
  <title></title>
</head>
<body bgcolor="#ffffff" text="#000000">
Hi Paul,<br>
I try to better explain me.<br>
In the grass version I have (grass-6.0.0beta2) when I interpolate a
vector point file by using v.surf.idw (with or wothout the -n flags)
the points that are not in the current region are not considered.<br>
What I want is to consider also the points outside the region to
generate the interpolated map.<br>
<br>
I understand that it is not possible to consider them (poin out of
region) by using the index point because they doesn't fall in any cell
(otherwise segmentation fault happens!) but with the old method (not
index cell) it is possible to consider also point out of region.<br>
<br>
What I simply modify in the code is to add a -a flag that causes the
use of old method (not index cell) with all the points in surface
computation. All the rest is still the same.....<br>
<br>
I'm I in wrong?<br>
Let me know.<br>
Ciao,<br>
Maxi<br>
<br>
<br>
<br>
Paul Kelly wrote:<br>
<blockquote
 cite="midPine.LNX.4.60.0503021010130.32644@agrippa.ukshells.co.uk"
 type="cite">On Wed, 2 Mar 2005, Massimiliano Cannata wrote:
  <br>
  <br>
  <blockquote type="cite">So, I've applied the change and things seems
to work now.
    <br>
In attachment i post the modified main.c file.
    <br>
Hope it can be usefull.....
    <br>
    <br>
Maxi
    <br>
    <br>
Massimiliano Cannata wrote:
    <br>
    <br>
    <blockquote type="cite">I see that in function new point (line 372)
the author exclude points outside the region due to segmentation fault
cased from new cell index solution.
      <br>
    </blockquote>
  </blockquote>
  <br>
Not sure what you mean here---in the new method the points are indexed
according to which cell they fall into. It is not possible to index
points that are outside the current region in this way (because they
don't fall inside any cells) and so they are not included.
  <br>
  <br>
The -n flag doesn't index the points at all; it just uses all the
points in the vector/sites file.
  <br>
  <br>
  <blockquote type="cite">
    <blockquote type="cite">My idea is to set a flag&nbsp; (-a interpolate
trough all the points) that automatically select the old interpolation
method to avoid any segmentation fault.......
      <br>
    </blockquote>
  </blockquote>
  <br>
Please explain what you mean about the segmentation fault.
  <br>
  <br>
  <blockquote type="cite">
    <blockquote type="cite">so it should looks:
      <br>
-n means use old interpolation method
      <br>
-a means use old interpolation method over all the sites
      <br>
    </blockquote>
  </blockquote>
  <br>
Yes but this is actually what -n does.
  <br>
  <br>
  <blockquote type="cite">
    <blockquote type="cite">else use new method over points in region
      <br>
      <br>
Impression and comments are welcome....
      <br>
Maxi
      <br>
    </blockquote>
  </blockquote>
  <br>
I'm not sure I understand what the improvement does. If you could
explain better what the problem is with the existing version, it might
help.
  <br>
  <br>
Paul
  <br>
  <br>
  <br>
</blockquote>
<br>
<div class="moz-signature">-- <br>
<DEFANGED_meta http-equiv="content-type" content="text/html; charset=ISO-8859-1">
<title>Signature</title>
<span DEFANGED_STYLE="color: rgb(102, 102, 102);">_____________________________________________________
</span><br DEFANGED_STYLE="color: rgb(102, 102, 102);">
<br DEFANGED_STYLE="color: rgb(102, 102, 102);">
<span DEFANGED_STYLE="color: rgb(102, 102, 102);">Ing. Massimiliano Cannata
</span><br DEFANGED_STYLE="color: rgb(102, 102, 102);">
<br DEFANGED_STYLE="color: rgb(102, 102, 102);">
<span DEFANGED_STYLE="color: rgb(102, 102, 102);">SUPSI - Istituto Scienze della
Terra
</span><br DEFANGED_STYLE="color: rgb(102, 102, 102);">
<span DEFANGED_STYLE="color: rgb(102, 102, 102);">C.P. 72 - CH-6952 Canobbio
(Ticino, Switzerland)
</span><br DEFANGED_STYLE="color: rgb(102, 102, 102);">
<span DEFANGED_STYLE="color: rgb(102, 102, 102);">Tel +41 (0)91 935 12 25 - Fax
+41 (0)91 935 12 09
</span><br DEFANGED_STYLE="color: rgb(102, 102, 102);">
<a class="moz-txt-link-freetext"
 href="mailto:massimiliano.cannata@supsi.ch"
 DEFANGED_STYLE="color: rgb(102, 102, 102);">mailto:massimiliano.cannata@supsi.ch</a><span
 style="color: rgb(102, 102, 102);">
</span><br DEFANGED_STYLE="color: rgb(102, 102, 102);">
<a class="moz-txt-link-freetext" href="http://www.ist.supsi.ch"
 DEFANGED_STYLE="color: rgb(102, 102, 102);">http://www.ist.supsi.ch</a><span
 style="color: rgb(102, 102, 102);">
</span><br DEFANGED_STYLE="color: rgb(102, 102, 102);">
<br>
</div>
</body>
</html>