Python Geospatial Analysis Cookbook

谢谢分享

kankan

谢谢！

要看内容就要回复，既然是这样的话，那我就回复一下...

好书啊！谢谢！

1. How to do it...
   
2. 1. Now, let's get into the actual importing of our Shapefile into a PostGIS database that
   
3. will automatically create a new table from our Shapefile:
   
4. #!/usr/bin/env python
   
5. # -*- coding: utf-8 -*-
   
6. import subprocess
   
7. # database options
   
8. db_schema = "SCHEMA=geodata"
   
9. overwrite_option = "OVERWRITE=YES"
   
10. geom_type = "MULTILINESTRING"
    
11. output_format = "PostgreSQL"
    
12. # database connection string
    
13. db_connection = """PG:host=localhost port=5432
    
14. user=pluto dbname=py_test password=stars"""
    
15. # input shapefile
    
16. input_shp = "../geodata/bikeways.shp"
    
17. # call ogr2ogr from python
    
18. subprocess.call(["ogr2ogr","-lco", db_schema, "-lco",
    
19. overwrite_option,
    
20. "-nlt", geom_type, "-f", output_format, db_connection,
    
21. input_shp])
    
22. 2. Next, we'll call our script from the command line:
    
23. $ python ch03-01_shp2pg.py

复制代码