Adds linear interpolation to grid_tile:ing of dtms. Hopefully results look a bit better as a result

dvlp/grid_tile.sh

@@ -140,7 +140,8 @@ function grid_shape() {
     xmax=$5
     ymax=$6
     ncells=$7
-    method="invdist"
+    method=$8
+    
 
     if [ ! -f $tiffile ]; then
 	gdal_grid \
@@ -319,7 +320,7 @@ fi
 
 # interpolate raster from dtm points
 if [ ! -f roofs.tif ]; then
-    grid_shape  $UNCLASSIFIED_POINTS_IN_BUILDINGS_SHAPE roofs.tif $x1 $y1 $x2 $y2 $ncells
+    grid_shape  $UNCLASSIFIED_POINTS_IN_BUILDINGS_SHAPE roofs.tif $x1 $y1 $x2 $y2 $ncells "invdist"
 fi
 
 # fix for bug in gdal_grid causing rasters with strange orientation, that cannot be handled by gdal_calc.py
@@ -327,7 +328,7 @@ gdalwarp roofs.tif roofs_corr.tif
 
 # interpolate raster from dem points
 if [ ! -f $dem ]; then
-    grid_shape  $DEM_POINTS_OUTSIDE_BUILDINGS_SHAPE dem.tif $x1 $y1 $x2 $y2 $ncells
+    grid_shape  $DEM_POINTS_OUTSIDE_BUILDINGS_SHAPE dem.tif $x1 $y1 $x2 $y2 $ncells "linear"
 fi
 
 # fix for bug in gdal_grid causing rasters with strange orientation, that cannot be handled by gdal_calc.py