Add raster shadedrelief command

Author: jericks

README.rst

@@ -255,6 +255,8 @@ raster divide raster --name1 topo --name2 reduce --output-format topo_new --outp
 
 raster stylize --name raster --output-format raster_stylized --output-name raster_stylized
 
+raster shadedrelief --name raster --output-format shaded --output-name shaded --scale 1.0 --altitude 25.0 --azimuth 35.0
+
 style
 -----
 style create --params "stroke=navy stroke=width=5" --file earth_outline.sld

examples/raster_shadedrelief.txt

@@ -0,0 +1,17 @@
+format open --name raster --input raster.tif
+
+raster open --format raster --raster raster --name raster
+
+format open --name shaded --input shaded.tif
+
+raster shadedrelief --name raster --output-format shaded --output-name shaded --scale 1.0 --altitude 25.0 --azimuth 35.0
+
+map open --name map
+
+map add raster --name map --raster shaded
+
+map draw --name map
+
+map close --name map
+
+open --file image.png

src/main/groovy/org/geoshell/raster/RasterCommands.groovy

@@ -603,4 +603,36 @@ class RasterCommands implements CommandMarker {
         }
     }
 
+    @CliCommand(value = "raster shadedrelief", help = "Create a shaded relief raster")
+    String shadedRelief(
+            @CliOption(key = "name", mandatory = true, help = "The Raster name") RasterName name,
+            @CliOption(key = "output-format", mandatory = true, help = "The output Format Workspace") FormatName formatName,
+            @CliOption(key = "output-name", mandatory = false, help = "The output Raster name") String outputRasterName,
+            @CliOption(key = "scale", mandatory = true, help = "The scale") double scale,
+            @CliOption(key = "altitude", mandatory = true, help = "The altitude") double altitude,
+            @CliOption(key = "azimuth", mandatory = true, help = "The azimuth") double azimuth,
+            @CliOption(key = "resx", mandatory = false, specifiedDefaultValue = "0.5", unspecifiedDefaultValue = "0.5", help = "The x resolution") double resX,
+            @CliOption(key = "resy", mandatory = false, specifiedDefaultValue = "0.5", unspecifiedDefaultValue = "0.5", help = "The y resolution") double resY,
+            @CliOption(key = "zetafactory", mandatory = false, specifiedDefaultValue = "1.0", unspecifiedDefaultValue = "1.0", help = "The zeta factory") double zetaFactory,
+            @CliOption(key = "algorithm", mandatory = false, specifiedDefaultValue = "DEFAULT", unspecifiedDefaultValue = "DEFAULT", help = "The x resolution") String algorithm
+    ) throws Exception {
+        Raster raster = catalog.rasters[name]
+        if (raster) {
+            Format format = catalog.formats[formatName]
+            if (format) {
+                Raster shadedReliefRaster = raster.createShadedRelief(scale, altitude, azimuth, resX: resX, resY: resY, zetaFactory: zetaFactory, algorithm: algorithm)
+                format.write(shadedReliefRaster)
+                if (!outputRasterName) {
+                    outputRasterName = formatName.name
+                }
+                catalog.rasters[new RasterName(outputRasterName)] = format.read(outputRasterName)
+                "Create shaded relief ${outputRasterName} from ${name}!"
+            } else {
+                "Unable to find Raster Format ${formatName}"
+            }
+        } else {
+            "Unable to find Raster ${name}"
+        }
+    }
+
 }

src/test/groovy/org/geoshell/raster/RasterCommandsTest.groovy

@@ -485,4 +485,28 @@ class RasterCommandsTest {
         assertEquals(214.0, outRaster.getValue(10,10,0), 0.1)
         assertEquals(212.0, outRaster.getValue(30,20,0), 0.1)
     }
+
+    @Test void shadedRelief() {
+        Catalog catalog = new Catalog()
+        File file = new File(getClass().getClassLoader().getResource("raster.tif").toURI())
+        Format format = Format.getFormat(file)
+        catalog.formats[new FormatName("raster")] = format
+        catalog.rasters[new RasterName("raster")] = catalog.formats[new FormatName("raster")].read()
+
+        File outFile = new File(temporaryFolder.newFolder("shaded"), "shaded.tif")
+        Format outFormat = Format.getFormat(outFile)
+        catalog.formats[new FormatName("shaded")] = outFormat
+
+        RasterCommands cmds = new RasterCommands(catalog: catalog)
+        cmds.open(new FormatName("raster"), new RasterName("raster"), "raster")
+
+        String result = cmds.shadedRelief(new RasterName("raster"), new FormatName("shaded"), "shaded", 1.0, 15.0, 15.0, 0.5, 0.5, 1.0, "DEFAULT")
+        assertEquals("Create shaded relief shaded from raster!", result)
+        Raster outRaster = catalog.rasters[new RasterName("shaded")]
+        assertNotNull outRaster
+        assertEquals(67.0, outRaster.getValue(0,0,0),  0.1)
+        assertEquals(67.0, outRaster.getValue(10,10,0), 0.1)
+        assertEquals(1.0, outRaster.getValue(30,20,0), 0.1)
+    }
+
 }