#TRIDataset.Save(OutputFolderPath+"TRI.tif") SpaTopo.Slope(TheDataset, OutputFilePath=OutputFolderPath + "Slope.tif") SpaView.Show(OutputFolderPath + "Slope.tif") SpaTopo.Aspect(TheDataset, OutputFilePath=OutputFolderPath + "Aspect.tif") SpaView.Show(OutputFolderPath + "Aspect.tif") SlopeDataset = SpaRasters.Load(OutputFolderPath + "Slope.tif") AspectDataset = SpaRasters.Load(OutputFolderPath + "Aspect.tif") SteepDataset = SlopeDataset > 20 SpaView.Show(SteepDataset) GreaterThan90Degrees = AspectDataset > 90 SpaView.Show(GreaterThan90Degrees) LessThan270Degres = AspectDataset < 270 SpaView.Show(LessThan270Degres) SouthDataset = SpaRasters.And(GreaterThan90Degrees, LessThan270Degres) SpaView.Show(SouthDataset) NewDataset3 = SpaRasters.And(SouthDataset, SteepDataset) SpaView.Show(NewDataset3) # works NewDataset = SpaTopo.Hillshade(Path1) NewDataset.Save(OutputFolderPath + "Hillshade.tif") SpaView.Show(NewDataset)
NewDataset=SpaRasters.Multiply(RasterFilePath,RasterFile2) NewDataset.Save(TempFolderPath+"Multiple.tif") NewDataset=SpaRasters.Multiply(RasterFilePath,3) NewDataset.Save(TempFolderPath+"Multiple1.tif") NewDataset=SpaRasters.Multiply(1.0,RasterFile2) NewDataset.Save(TempFolderPath+"Multiple2.tif") #Test Logical Operators if (True): #create some boolean rasters to work with RasterFile2 =SpaRasters.LessThan(RasterFilePath,5) BoolRaster1=SpaRasters.GreaterThan(RasterFilePath,RasterFile2) BoolRaster2=SpaRasters.LessThan(RasterFilePath,RasterFile2) NewDataset=SpaRasters.And(BoolRaster1,BoolRaster2) NewDataset.Save(TempFolderPath + "And.tif") NewDataset=SpaRasters.Or(BoolRaster1,BoolRaster2) NewDataset.Save(TempFolderPath + "Or.tif") NewDataset=SpaRasters.Not(BoolRaster1) NewDataset.Save(TempFolderPath + "Not.tif") #Test Rounding Functions: if (True): #create a raster with decimal numbers DecimalRaster=SpaRasters.Divide(RasterFilePath,RasterFile2) DecimalRaster.Save(TempFolderPath + "Divide.tif") NewDataset=SpaRasters.Round(DecimalRaster, 1)