def createRaster():
rows = 40
cols = 60
xorg = 0.
yorg = 0.
xp = 5
yp = 5
nodata = -999.999
cellsize = 1.
levels = 4
datahi = 100.
datalow = 0
randpercent = 0.2
resampleFactorA = 1
raster = createRanRasterSlope(rows, cols, cellsize, xorg, yorg, nodata,
levels, datahi, datalow, xp, yp, randpercent)
data = raster.getData()
rainraster = createRanRasterSlope(rows // resampleFactorA,
cols // resampleFactorA,
cellsize * resampleFactorA, xorg, yorg,
nodata, levels, 4000, 1, 36, 4, .1)
return raster, data, rainraster
@author: nrjh
"""
from RasterHandler import createRanRasterSlope
import matplotlib.pyplot as mp
from Flow import flowRaster
rows=10
cols=15
xorg=0.
yorg=0.
nodata=-999.999
cellsize=1.
#raster=createRanRaster()
raster=createRanRasterSlope(rows,cols,cellsize,xorg,yorg,nodata,3,100.,0.,1,1,0.2)
data=raster.getData()
mp.matshow(data)
mp.colorbar()
fr=flowRaster(raster.getData(),raster.getOrgs()[0],raster.getOrgs()[0],raster.getCellsize(),raster.getNoData())
fr.setDownCells()
pointlist=fr.getPointList()
for p in pointlist:
mp.scatter(p.get_x(),p.get_y())
for p in pointlist:
############# step 1 to 4 #######################################
# Create Random Raster
rows = 40
cols = 40
xorg = 0.
yorg = 0.
xp = 5
yp = 5
nodata = -999.999
cellsize = 1.
levels = 4
datahi = 100.
datalow = 0
randpercent = 0.1
resampleFactorA = 1
elevationRasterA = createRanRasterSlope(rows, cols, cellsize, xorg, yorg,
nodata, levels, datahi, datalow, xp,
yp, randpercent)
rainrasterA = createRanRasterSlope(rows // resampleFactorA,
cols // resampleFactorA,
cellsize * resampleFactorA, xorg, yorg,
nodata, levels, 4000, 1, 36, 4, .1)
##random raster
calculateFlowsAndPlot(elevationRasterA, rainrasterA, resampleFactorA)
############# step 5 #######################################
#calculateFlowsAndPlot(readRaster('ascifiles/dem_hack.txt'), readRaster('ascifiles/rain_small_hack.txt'), 10)
cols=30
xorg=0.
yorg=0.
xp=5
yp=5
nodata=-999.999
cellsize=1.
levels=4
datahi=100.
datalow=0
randpercent=0.2
colours=["black","red","blue","yellow","green","cyan","white","orange","grey","brown"]
#raster=createRanRaster()
raster=createRanRasterSlope(rows,cols,cellsize,xorg,yorg,nodata,levels,datahi,datalow,xp,yp,randpercent)
data=raster.getData()
mp.matshow(data)
mp.colorbar()
fr=flowRaster(raster.getData(),raster.getOrgs()[0],raster.getOrgs()[0],raster.getCellsize(),raster.getNoData())
fr.setDownCells()
pointlist=fr.getPointList()
for p in pointlist:
mp.scatter(p.get_x(),p.get_y())
#for p in pointlist: