Add raster mosaic command

Author: jericks

README.rst

@@ -257,6 +257,8 @@ raster stylize --name raster --output-format raster_stylized --output-name raste
 
 raster shadedrelief --name raster --output-format shaded --output-name shaded --scale 1.0 --altitude 25.0 --azimuth 35.0
 
+raster mosaic --name1 raster1 --name2 raster2 --output-format mosaic --output-name mosaic
+
 style
 -----
 style create --params "stroke=navy stroke=width=5" --file earth_outline.sld

examples/alki2.tif

@@ -0,0 +1,21 @@
+format open --input alki2.tif --name alki2
+
+raster open --format alki2 --raster alki2 --name alki2
+
+format open --input alki3.tif --name alki3
+
+raster open --format alki3 --raster alki3 --name alki3
+
+format open --input mosaic.tif --name mosaic
+
+raster mosaic --name1 alki2 --name2 alki3 --output-format mosaic --output-name mosaic
+
+map open --name map
+
+map add raster --name map --raster mosaic
+
+map draw --name map
+
+map close --name map
+
+open --file image.png

examples/alki3.tif

@@ -635,4 +635,28 @@ class RasterCommands implements CommandMarker {
         }
     }
 
+    @CliCommand(value = "raster mosaic", help = "Mosaic two Rasters together")
+    String mosaic(
+            @CliOption(key = "name1", mandatory = true, help = "The Raster name") RasterName name1,
+            @CliOption(key = "name2", mandatory = true, help = "The Raster name") RasterName name2,
+            @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
+    ) throws Exception {
+        Raster raster1 = catalog.rasters[name1]
+        Raster raster2 = catalog.rasters[name2]
+        Format format = catalog.formats[formatName]
+        if (format) {
+            Raster mosaicRaster = Raster.mosaic([raster1, raster2])
+            format.write(mosaicRaster)
+            if (!outputRasterName) {
+                outputRasterName = formatName.name
+            }
+            catalog.rasters[new RasterName(outputRasterName)] = format.read(outputRasterName)
+            "Mosaic ${name1} and ${name2} to create ${outputRasterName}!"
+        } else {
+            "Unable to find Raster Format ${formatName}"
+        }
+
+    }
+
 }

examples/raster_mosaic.txt

@@ -1,5 +1,6 @@
 package org.geoshell.raster
 
+import geoscript.geom.Bounds
 import geoscript.layer.Format
 import geoscript.layer.Layer
 import geoscript.layer.Raster
@@ -509,4 +510,31 @@ class RasterCommandsTest {
         assertEquals(1.0, outRaster.getValue(30,20,0), 0.1)
     }
 
+    @Test void mosaic() {
+
+        Catalog catalog = new Catalog()
+        File file1 = new File(getClass().getClassLoader().getResource("mosaic/alki2.tif").toURI())
+        Format format1 = Format.getFormat(file1)
+        catalog.formats[new FormatName("raster1")] = format1
+        catalog.rasters[new RasterName("raster1")] = catalog.formats[new FormatName("raster1")].read()
+
+        File file2 = new File(getClass().getClassLoader().getResource("mosaic/alki3.tif").toURI())
+        Format format2 = Format.getFormat(file2)
+        catalog.formats[new FormatName("raster2")] = format2
+        catalog.rasters[new RasterName("raster2")] = catalog.formats[new FormatName("raster2")].read()
+
+        File outFile = new File(temporaryFolder.newFolder("mosaic"), "mosaic.tif")
+        Format outFormat = Format.getFormat(outFile)
+        catalog.formats[new FormatName("mosaic")] = outFormat
+
+        RasterCommands cmds = new RasterCommands(catalog: catalog)
+        cmds.open(new FormatName("raster1"), new RasterName("raster1"), "raster1")
+        cmds.open(new FormatName("raster2"), new RasterName("raster2"), "raster2")
+
+        String result = cmds.mosaic(new RasterName("raster1"), new RasterName("raster2"), new FormatName("mosaic"), "mosaic")
+        assertEquals("Mosaic raster1 and raster2 to create mosaic!", result)
+        Raster outRaster = catalog.rasters[new RasterName("mosaic")]
+        assertNotNull outRaster
+    }
+
 }