<!DOCTYPE html PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN">

<html>

<head>

  <meta content="text/html;charset=ISO-8859-3" http-equiv="Content-Type">

  <title></title>

</head>

<body bgcolor="#ffffff" text="#000000">

Maciej Sieczka wrote:

<blockquote cite="mid:46D44BA4.5010901@o2.pl" type="cite">

  <pre wrap="">Luigi Ponti wrote:

  </pre>

  <br>

  <blockquote type="cite">

    <pre wrap="">What I meant is that the outer Voronoi cells may extend up to the 

region's bounding box, with consequent new Voronoi cells for 

interpolation points outside the convex hull, please see 

<a class="moz-txt-link-freetext" href="http://quartese.googlepages.com/voronoi">http://quartese.googlepages.com/voronoi</a>

    </pre>

  </blockquote>

  <pre wrap=""><!----></pre>

</blockquote>

<blockquote cite="mid:46D44BA4.5010901@o2.pl" type="cite">

  <pre wrap="">

Anyway - what could be done to overcome the trimming to convex hull

effect: you need to make sure that the convex hull for your input cells

covers your whole area of interest - ie. enlarge region enough for

important input cells not to be omitted. If cells in the input raster

are distributed so that a convex hull covering the whole area of

interest is not possible, you need to enhance your input raster.</pre>

</blockquote>

Please, have also a look at<br>

<a class="moz-txt-link-freetext" href="http://www.ems-i.com/gmshelp/Interpolation/Interpolation_Schemes/Natural_Neighbor_Interpolation.htm">http://www.ems-i.com/gmshelp/Interpolation/Interpolation_Schemes/Natural_Neighbor_Interpolation.htm</a><br>

<br>

This web page suggests an approach for extrapolation that would get a

non-null value to all the points in the region: four additional input

points are included in the NN interpolation that correspond to the

corners of a user-defined region (a bounding box that includes all

input points), by assigning them a value based on a different

interpolation method (IDW). This way, the convex hull of input points

has the same extent of the bounding box of the region. Also, I think

that there is no need of giving the -W nnbathy flag a negative value

and the resulting interpolated raster would still be inside the range

of the input data because of the way IDW works.<br>

<br>

It sounds feasible to add this functionality to r.surf.nnbathy, but I

don't know if the approach is correct from the GIS point of view (i.e.,

mixing two interpolation methods) or even more useful/better than just

using IDW for the whole region.<br>

<br>

Thanks and regards,<br>

<br>

Luigi<br>

</body>

</html>