def Show(TheDataset, width=800, height=800): # IF a file path is specified, load the dataset if (type(TheDataset) is str): filename, file_extension = os.path.splitext(TheDataset) file_extension = file_extension.lower() if (file_extension == ".shp"): TheDataset = SpaVectors.Load( TheDataset) # load the contents of the layer else: TheDataset = SpaRasters.Load(TheDataset) # Create the view TheView = SpaView(width, height) # Set the bounds of the view to match the dataset TheBounds = TheDataset.GetBounds() TheView.SetBounds(TheBounds) # Render the dataset into the view if (isinstance(TheDataset, SpaRasters.SpaDatasetRaster)): TheView.RenderRaster(TheDataset) elif (isinstance(TheDataset, SpaVectors.SpaDatasetVector)): TheView.RenderVectors(TheDataset) TheView.Show()
############################################################################ # Globals ############################################################################ OutputFolderPath = "../Temp/" RasterFilePath = "../Data/MtStHelens/Mt St Helens Post Eruption DEM Int16.tif" #RasterFilePath2="../Data/MtStHelens/Mt St Helens PreEruption DEM Float32.tif" Path1 = "../Data/MtStHelens/Mt St Helens PreEruption DEM Float32.tif" #Path2="../Data/MtStHelens/Mt St Helens Post Eruption DEM.tif" ############################################################################ # SpaView Tests ############################################################################ TheDataset = SpaRasters.Load(Path1) SpaTopo.Contour(Path1, OutputFilePath=OutputFolderPath + "Countours.shp") #SpaView.Show(OutputFolderPath+"Countours.shp") SpaTopo.TRI(Path1, OutputFilePath=OutputFolderPath + "TRI.tif") SpaView.Show(OutputFolderPath + "TRI.tif") # deal with no data values #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")
TheBounds = TheDataset.GetBounds() TheView.SetBounds(TheBounds) TheView.RenderRaster(TheDataset) TheView.Save(OutputFolderPath + "Raster.png") TheView.Show() ####################################################################### # find missing area of mt st helens if (True): import math SpaView.Show(Path1) # deal with no data values TheDataset = SpaRasters.Load(Path1) SpaTopo.Contour(Path1, OutputFilePath=OutputFolderPath + "Countours.shp") TRIDataset = SpaTopo.TRI(Path1) SpaView.Show(TRIDataset) # deal with no data values TRIDataset.Save(OutputFolderPath + "TRI.tif") SlopeDataset = SpaTopo.Slope(TheDataset) SpaView.Show(SlopeDataset) AspectDataset = SpaTopo.Aspect(TheDataset) SpaView.Show(AspectDataset) SteepDataset = SlopeDataset > 20 SpaView.Show(SteepDataset)