Advertisement
Guest User

Untitled

a guest
Jan 23rd, 2017
101
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 1.26 KB | None | 0 0
  1. ##Cell Statistics=name
  2. ##Select_directory=Folder
  3. ##statistic_name=string average
  4. ##Output_raster=Output raster
  5. import glob, os
  6. from PyQt4.QtCore import QFileInfo
  7. from qgis.core import QgsRasterLayer, QgsRectangle
  8.  
  9. os.chdir(Select_directory)
  10. raster_list = []
  11. extent = QgsRectangle()
  12. extent.setMinimal()
  13.  
  14. statistic_dic = {'average':0, 'count':1, 'median':2,'mode':3, 'minimum':4, 'min_raster':5,'maximum':6, 'max_raster':7, 'std':8, 'range':9}
  15.  
  16.  
  17. for raster in glob.glob("*.tif"):
  18. fileInfo = QFileInfo(raster)
  19. baseName = fileInfo.baseName()
  20. rlayer = QgsRasterLayer(raster, baseName)
  21. # Combine raster layers to list
  22. raster_list .append(rlayer)
  23. # Combine raster extents
  24. extent.combineExtentWith(rlayer.extent())
  25.  
  26. # Get extent
  27. xmin = extent.xMinimum()
  28. xmax = extent.xMaximum()
  29. ymin = extent.yMinimum()
  30. ymax = extent.yMaximum()
  31. # Run algorithm and set relevant parameters
  32. processing.runalg("grass7:r.series",
  33. {"input":raster_list ,
  34. "-n":False,
  35. "method":statistic_dic[statistic_name],
  36. "range":'-10000000000,10000000000',
  37. "GRASS_REGION_PARAMETER":"%f,%f,%f,%f"% (xmin, xmax, ymin, ymax),
  38. "GRASS_REGION_CELLSIZE_PARAMETER":0,
  39. "output":Output_raster})
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement