Advertisement
Guest User

Untitled

a guest
Dec 10th, 2019
109
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 4.12 KB | None | 0 0
  1. # script to apply fuzzy membership to climate grids and map malaria in Ethiopia
  2.  
  3. import arcpy
  4. arcpy.CheckOutExtension ("Spatial")
  5. arcpy.env.overwriteOutput = True
  6. arcpy.env.workspace = r"H:\MRes\BatchprocessingGuide\Ethiopia\Climate"
  7.  
  8. # processing average temperature
  9.  
  10. for inRaster in arcpy.ListRasters("tavg*.tif"):
  11. fuzzyOutput = r"H:\MRes\BatchprocessingGuide\Malaria\MinAvg" + inRaster [5:]
  12. fuzzyAlgorithm = "LINEAR 18 22"
  13. arcpy.gp.FuzzyMembership_sa(inRaster, fuzzyOutput, fuzzyAlgorithm)
  14. print ("fuzzyOutput")
  15.  
  16. # processing minimum temperature
  17.  
  18. for inRaster in arcpy.ListRasters("tmin*.tif"):
  19. fuzzyOutput = r"H:\MRes\BatchprocessingGuide\Malaria\MinT" + inRaster [5:]
  20. fuzzyAlgorithm = "LINEAR 4 6"
  21. arcpy.gp.FuzzyMembership_sa(inRaster, fuzzyOutput, fuzzyAlgorithm)
  22. print ("fuzzyOutput")
  23.  
  24. # processing rainfall
  25.  
  26. for inRaster in arcpy.ListRasters("prec*.tif"):
  27. fuzzyOutput = r"H:\MRes\BatchprocessingGuide\Malaria\rain" + inRaster [5:]
  28. fuzzyAlgorithm = "LINEAR 0 80"
  29. arcpy.gp.FuzzyMembership_sa(inRaster, fuzzyOutput, fuzzyAlgorithm)
  30. print ("fuzzyOutput")
  31.  
  32. #monthly suitability
  33.  
  34. arcpy.gp.RasterCalculator_sa(r' "H:\MRes\BatchprocessingGuide\Malaria\MinAvg01.tif" * "H:\MRes\BatchprocessingGuide\Malaria\MinT01.tif" * "H:\MRes\BatchprocessingGuide\Malaria\rain01.tif" ', "H:\MRes\BatchprocessingGuide\Malaria\Monthsuitable01.tif")
  35. arcpy.gp.RasterCalculator_sa(r' "H:\MRes\BatchprocessingGuide\Malaria\MinAvg02.tif" * "H:\MRes\BatchprocessingGuide\Malaria\MinT02.tif" * "H:\MRes\BatchprocessingGuide\Malaria\rain02.tif" ', "H:\MRes\BatchprocessingGuide\Malaria\Monthsuitable02.tif")
  36. arcpy.gp.RasterCalculator_sa(r' "H:\MRes\BatchprocessingGuide\Malaria\MinAvg03.tif" * "H:\MRes\BatchprocessingGuide\Malaria\MinT03.tif" * "H:\MRes\BatchprocessingGuide\Malaria\rain03.tif" ', "H:\MRes\BatchprocessingGuide\Malaria\Monthsuitable03.tif")
  37. arcpy.gp.RasterCalculator_sa(r' "H:\MRes\BatchprocessingGuide\Malaria\MinAvg04.tif" * "H:\MRes\BatchprocessingGuide\Malaria\MinT04.tif" * "H:\MRes\BatchprocessingGuide\Malaria\rain04.tif" ', "H:\MRes\BatchprocessingGuide\Malaria\Monthsuitable04.tif")
  38. arcpy.gp.RasterCalculator_sa(r' "H:\MRes\BatchprocessingGuide\Malaria\MinAvg05.tif" * "H:\MRes\BatchprocessingGuide\Malaria\MinT05.tif" * "H:\MRes\BatchprocessingGuide\Malaria\rain05.tif" ', "H:\MRes\BatchprocessingGuide\Malaria\Monthsuitable05.tif")
  39. arcpy.gp.RasterCalculator_sa(r' "H:\MRes\BatchprocessingGuide\Malaria\MinAvg06.tif" * "H:\MRes\BatchprocessingGuide\Malaria\MinT06.tif" * "H:\MRes\BatchprocessingGuide\Malaria\rain06.tif" ', "H:\MRes\BatchprocessingGuide\Malaria\Monthsuitable06.tif")
  40. arcpy.gp.RasterCalculator_sa(r' "H:\MRes\BatchprocessingGuide\Malaria\MinAvg07.tif" * "H:\MRes\BatchprocessingGuide\Malaria\MinT07.tif" * "H:\MRes\BatchprocessingGuide\Malaria\rain07.tif" ', "H:\MRes\BatchprocessingGuide\Malaria\Monthsuitable07.tif")
  41. arcpy.gp.RasterCalculator_sa(r' "H:\MRes\BatchprocessingGuide\Malaria\MinAvg08.tif" * "H:\MRes\BatchprocessingGuide\Malaria\MinT08.tif" * "H:\MRes\BatchprocessingGuide\Malaria\rain08.tif" ', "H:\MRes\BatchprocessingGuide\Malaria\Monthsuitable08.tif")
  42. arcpy.gp.RasterCalculator_sa(r' "H:\MRes\BatchprocessingGuide\Malaria\MinAvg09.tif" * "H:\MRes\BatchprocessingGuide\Malaria\MinT09.tif" * "H:\MRes\BatchprocessingGuide\Malaria\rain09.tif" ', "H:\MRes\BatchprocessingGuide\Malaria\Monthsuitable09.tif")
  43. arcpy.gp.RasterCalculator_sa(r' "H:\MRes\BatchprocessingGuide\Malaria\MinAvg10.tif" * "H:\MRes\BatchprocessingGuide\Malaria\MinT10.tif" * "H:\MRes\BatchprocessingGuide\Malaria\rain10.tif" ', "H:\MRes\BatchprocessingGuide\Malaria\Monthsuitable10.tif")
  44. arcpy.gp.RasterCalculator_sa(r' "H:\MRes\BatchprocessingGuide\Malaria\MinAvg11.tif" * "H:\MRes\BatchprocessingGuide\Malaria\MinT11.tif" * "H:\MRes\BatchprocessingGuide\Malaria\rain11.tif" ', "H:\MRes\BatchprocessingGuide\Malaria\Monthsuitable11.tif")
  45. arcpy.gp.RasterCalculator_sa(r' "H:\MRes\BatchprocessingGuide\Malaria\MinAvg12.tif" * "H:\MRes\BatchprocessingGuide\Malaria\MinT12.tif" * "H:\MRes\BatchprocessingGuide\Malaria\rain12.tif" ', "H:\MRes\BatchprocessingGuide\Malaria\Monthsuitable12.tif")
  46. print ("suitable months done")
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement