Advertisement
Guest User

Untitled

a guest
Feb 8th, 2016
47
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 0.88 KB | None | 0 0
  1. inputFolder = "/home/usr/Desktop/BufferFolder/"
  2. outputFolder = "/home/usr/Desktop//BufferFolder/Addition"
  3. template = gdal.Open("/home/usr/Desktop//BufferFolder/Template.tif")
  4. rows = template.RasterYSize
  5. cols = template.RasterXSize
  6. if template is None:
  7. print 'Unable to open template file'
  8. sys.exit(1)
  9. templateArray = np.array(template.GetRasterBand(1).ReadAsArray())
  10. mainArray = 0 * templateArray
  11.  
  12. os.chdir(inputFolder)
  13.  
  14. def findRasters(path, filter):
  15. for root, dirs, files in os.walk(path, filter):
  16. for file in fnmatch.filter(files, filter):
  17. yield os.path.join(root, file)
  18.  
  19. count = 0
  20. for raster in findRasters(inputFolder, '*.tif'):
  21. (infilepath, infilename) = os.path.split(raster)
  22. openraster = gdal.Open(infilename)
  23. rasterToArray = np.array(openraster.GetRasterBand(1).ReadAsArray())
  24. mainArray = mainArray + rasterToArray
  25. count = count + 1
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement