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)