Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- # script to apply fuzzy membership to climate grids and map malaria in Ethiopia
- import arcpy
- arcpy.CheckOutExtension ("Spatial")
- arcpy.env.overwriteOutput = True
- arcpy.env.workspace = r"H:\MRes\BatchprocessingGuide\Ethiopia\Climate"
- # processing average temperature
- for inRaster in arcpy.ListRasters("tavg*.tif"):
- fuzzyOutput = r"H:\MRes\BatchprocessingGuide\Malaria\MinAvg" + inRaster [5:]
- fuzzyAlgorithm = "LINEAR 18 22"
- arcpy.gp.FuzzyMembership_sa(inRaster, fuzzyOutput, fuzzyAlgorithm)
- print ("fuzzyOutput")
- # processing minimum temperature
- for inRaster in arcpy.ListRasters("tmin*.tif"):
- fuzzyOutput = r"H:\MRes\BatchprocessingGuide\Malaria\MinT" + inRaster [5:]
- fuzzyAlgorithm = "LINEAR 4 6"
- arcpy.gp.FuzzyMembership_sa(inRaster, fuzzyOutput, fuzzyAlgorithm)
- print ("fuzzyOutput")
- # processing rainfall
- for inRaster in arcpy.ListRasters("prec*.tif"):
- fuzzyOutput = r"H:\MRes\BatchprocessingGuide\Malaria\rain" + inRaster [5:]
- fuzzyAlgorithm = "LINEAR 0 80"
- arcpy.gp.FuzzyMembership_sa(inRaster, fuzzyOutput, fuzzyAlgorithm)
- print ("fuzzyOutput")
- #monthly suitability
- 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")
- 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")
- 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")
- 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")
- 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")
- 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")
- 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")
- 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")
- 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")
- 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")
- 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")
- 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")
- print ("suitable months done")
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement