Guest User

Untitled

a guest
Feb 23rd, 2018
87
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 0.58 KB | None | 0 0
  1. temp = '/Users/lisa/Desktop/GIS/FISH/DATA/01.03.TIFF'
  2. chl = '/Users/lisa/Desktop/GIS/FISH/DATA/01.03CA.TIFF'
  3.  
  4. lyrt = processing.getObjectFromUri(temp)
  5. lyrc = processing.getObjectFromUri(chl)
  6.  
  7. entries = []
  8.  
  9. temp = QgsRasterCalculatorEntry()
  10. chl = QgsRasterCalculatorEntry()
  11.  
  12. temp.ref = 'lyr@t'
  13. chl.ref = 'lyr@c'
  14.  
  15. temp.raster = lyrt
  16. chl.raster = lyrc
  17.  
  18. entries.append( temp )
  19. entries.append( chl )
  20.  
  21. calc = QgsRasterCalculator( 'lyr@t > 26 AND lyr@t < 30 AND lyr@c > 0.3 AND lyr@c < 2.8', output_path, 'GTiff', lyrt.extent(), lyrt.width(), lyrt.height(), entries )
  22. calc.processCalculation()
Add Comment
Please, Sign In to add comment