[Github]ArcPy and ArcGIS Second Edition by Packt

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

本帖隐藏的内容

ArcPy-and-ArcGIS-Second-Edition-master.zip (366.88 KB)

2017-7-7 04:25:20 上传

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：Edition Second GitHub editio dition

1. # -*- coding: utf-8 -*-
   
2. # ---------------------------------------------------------------------------
   
3. # Chapter2Model1.py
   
4. # Created on: 2017-01-26 04:26:31.00000
   
5. #   (generated by ArcGIS/ModelBuilder)
   
6. # Description:
   
7. # ---------------------------------------------------------------------------
   
8. 
   
9. #Import libraries
   
10. import arcpy
    
11. import csv
    
12. 
    
13. # if uncommented, the next line allows for overwriting existing files with output. Be careful!
    
14. arcpy.env.overwriteOutput = True   
    
15. 
    
16. # Local variables:
    
17. Bus_Stops = r"C:\Projects\SanFrancisco.gdb\SanFrancisco\Bus_Stops"
    
18. CensusBlocks2010 = r"C:\Projects\SanFrancisco.gdb\SanFrancisco\CensusBlocks2010"
    
19. Inbound71 = r"C:\Projects\SanFrancisco.gdb\Chapter2Results\Inbound71"
    
20. Inbound71_400ft_buffer = r"C:\Projects\SanFrancisco.gdb\Chapter2Results\Inbound71_400ft_buffer"
    
21. Intersect71Census = r"C:\Projects\SanFrancisco.gdb\Chapter2Results\Intersect71Census"
    
22. 
    
23. # Process: Select
    
24. arcpy.Select_analysis(Bus_Stops,
    
25.                       Inbound71,
    
26.                       "NAME = '71 IB' AND BUS_SIGNAG = 'Ferry Plaza'")
    
27. 
    
28. # Process: Buffer
    
29. arcpy.Buffer_analysis(Inbound71,
    
30.                       Inbound71_400ft_buffer,
    
31.                       "400 Feet", "FULL", "ROUND", "NONE", "")
    
32. 
    
33. # Process: Intersect
    
34. arcpy.Intersect_analysis([Inbound71_400ft_buffer,CensusBlocks2010],
    
35.                          Intersect71Census, "ALL", "", "INPUT")
    
36. 
    
37. # Organize the census block population using a dictionary
    
38. # The stop ID is the key and the census block populations are added to a list as the value
    
39. dataDictionary = {}
    
40. with arcpy.da.SearchCursor(Intersect71Census, ["STOPID","POP10"]) as cursor:
    
41.     for row in cursor:
    
42.         busStopID = row[0]
    
43.         pop10 = row[1]
    
44.         if busStopID not in dataDictionary.keys():
    
45.             dataDictionary[busStopID] = [pop10]
    
46.         else:
    
47.             dataDictionary[busStopID].append(pop10)
    
48. 
    
49. #Create an output CSV spreadsheet. Check the folder path to be sure it's legal
    
50. with open(r'C:\Projects\Averages.csv', 'wb') as csvfile:
    
51.     csvwriter = csv.writer(csvfile, delimiter=',')
    
52.     for busStopID in dataDictionary.keys():
    
53.         popList = dataDictionary[busStopID]
    
54.         averagePop = sum(popList)/len(popList)
    
55.         data = [busStopID, averagePop]
    
56.         csvwriter.writerow(data)
    
57. 
    
58. print "Data Analysis Complete"

复制代码

ReneeBK 发表于 2017-7-7 04:26
**** 本内容被作者隐藏 ****

谢谢分享

thanks a lot