Guest User

Untitled

a guest
Nov 23rd, 2017
80
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 0.50 KB | None | 0 0
  1. import arcpy
  2. from arcpy.sa import *
  3. from arcpy import env
  4. arcpy.CheckOutExtension("Spatial")
  5. env.workspace = r"G:Place project8 changex2tempShp"
  6.  
  7. shplist = arcpy.ListFeatureClasses()
  8. for each_shapefile in shplist:
  9. name_shp = arcpy.Describe(each_shapefile).basename
  10. try: #PointDensity
  11. outPdens=PointDensity(each_shapefile,"NONE",10,"","")
  12. # Reciprocal
  13. outRaster=1+Raster(outPdens)
  14. outRaster.save(r"G:Place project"+"\"+str(name_shp))
  15. except:
  16. print "error in PointDensity"
Add Comment
Please, Sign In to add comment