New keywords

1.

edges. To use like:

  edges: "file1","file2",...;

Given boundary files will be read in and the boundaries will be used during neighborhood search. You can give both polylines and polygons.

2.

point-in-polygon. To use like:

  method: point-in-polygon;

With given point locations (through data statement), gstat will search which polygons the points are in. The output is a list of locations with a file number in the prediction column and a polygon number in the variance column of the output file. For points inside a polygon, this will be numbers counting from 0. Locations which are not inside of any polygon will have NA in both columns.