Advertisement
Guest User

Untitled

a guest
Aug 20th, 2019
90
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 0.51 KB | None | 0 0
  1. import arcpy
  2. import os
  3. ws = arcpy.env.workspace = r'D:BRB Snow CoverSampleTerra'
  4. Terra_Rasters = arcpy.ListRasters("MOD*")
  5. matches = 0
  6. for item1 in Terra_Rasters:
  7. for item2 in Terra_Rasters:
  8. print item1
  9. print item2
  10. if item1.endswith(item2[32:37]):
  11. print "Match"
  12. output = "D:\BRB Snow Cover\Sample\output""\" + str(item1) + '.tif'
  13. arcpy.gp.Con_sa(item1,item2,output,item1,""Value" =8")
  14. else:
  15. print "No"
  16. matches = matches+1
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement