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

<html>

<head>

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

</head>

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

Maciej Sieczka wrote:

<blockquote cite="mid:46D92701.6070203@o2.pl" type="cite">

  <pre wrap="">Luigi Ponti wrote:

  </pre>

  <blockquote type="cite">

    <pre wrap="">Maciej Sieczka wrote:

    </pre>

    <blockquote type="cite">

      <pre wrap="">Luigi Ponti wrote:

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.

      </pre>

    </blockquote>

  </blockquote>

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

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

It is possible to add such a functionality to r.surf.nnbathy, but I'd

rather leave it to the user to enhance his input raster the way he

prefers to, if necessary. Different interpolation methods or additional

data sources can be used to add the 4 corner points in question.

Cheers

Maciek

  </pre>

</blockquote>

<br>

Thanks Maciek,<br>

<br>

Any pointer to getting the four corners? Maybe they are already there,

i.e. the code that sets the working region for nnbathy?<br>

<blockquote># set the working region for nnbathy (it's cell-center

oriented)<br>

nn_n=`echo $n | awk -v res="$nsres" '{printf "%.8f",$1-res/2.0}'`<br>

nn_s=`echo $s | awk -v res="$nsres" '{printf "%.8f",$1+res/2.0}'`<br>

nn_w=`echo $w | awk -v res="$ewres" '{printf "%.8f",$1+res/2.0}'`<br>

nn_e=`echo $e | awk -v res="$ewres" '{printf "%.8f",$1-res/2.0}'`<br>

</blockquote>

Though I don't entirely understand what it does, it seems to be getting

the coordinates of the center of cells located at the 4 corners of the

region.<br>

<br>

Thanks again and regards,<br>

<br>

Luigi<br>

<br>

</body>

</html>