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: