Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- import arcpy
- from arcpy.sa import *
- from arcpy import env
- arcpy.CheckOutExtension("Spatial")
- env.workspace = r"G:Place project 8 changex2tempShp"
- shplist = arcpy.ListFeatureClasses()
- for each_shapefile in shplist:
- name_shp = arcpy.Describe(each_shapefile).basename
- try: #PointDensity
- outPdens=PointDensity(each_shapefile,"NONE",10,"","")
- # Reciprocal
- outRaster=1+Raster(outPdens)
- outRaster.save(r"G:Place project"+"\"+str(name_shp))
- except:
- print "error in PointDensity"
Add Comment
Please, Sign In to add comment