This library reads and writes ESRI Shapefiles in pure Python. You can read and write shp, shx, and dbf files with all types of geometry. Everything in the public ESRI shapefile specification is implemented. This library is compatible with Python versions 2.4 to 3.x.
Overview
This library reads and writes ESRI Shapefiles in pure Python. You can read and write shp, shx, and dbf files with all types of geometry. Everything in the public ESRI shapefile specification is implemented. This library is compatible with Python versions 2.4 to 3.x.
Get Started Instantly
- Download shapefile.py
- Start Python
- import shapefile
- Try one of the examples below
OR
Just run: easy_install pyshp
OR
pip install pyshp
If you are looking for information on .sbn and .sbx file formats some documentation is available here.
License
This library is released under the MIT license allowing it to be used for both commercial and non-commercial use. If your organization requires a different license for legal or policy reasons please contact the author through GeospatialPython.com.
Usage
The "Source" tab contains the SVN trunk with the latest release. Tagged releases representing relatively stable versions are also in the repository. The python module "shapefile.py" contains the complete library.
Detailed documentation and examples can be found in the Wiki as well as the "Download" tab. These files are different formats for the same information.
Below are minimal examples to give you a quick idea of what using the library "feels" like. More examples can also be found onGeospatialPython.com.
Important: For information about map projections, shapefiles, and Python please visit: http://code.google.com/p/pyshp/wiki/MapProjections
Reading Shapefiles
Reading Points in Shapes
>>> import shapefile>>> sf = shapefile.Reader("shapefiles/blockgroups") >>> shapes = sf.shapes() >>> # Read the bounding box from the 4th shape >>> shapes[3].bbox[-122.485792, 37.786931000000003, -122.446285, 37.811019000000002] >>># Read the 8th point in the 4th shape >>> shapes[3].points[7] [-122.471063, 37.787402999999998]
Reading Database Attributes
>>> # Read the field descriptors for the database file >>> sf.fields[("DeletionFlag", "C", 1, 0), ["AREA", "N", 18, 5], ... ["BKG_KEY", "C", 12, 0], ["POP1990", "N", 9, 0], ["POP90_SQMI", "N", 10, 1], ... ["HOUSEHOLDS", "N", 9, 0], ... ["MALES", "N", 9, 0], ["FEMALES", "N", 9, 0]] >>> # Read the 2nd and 3rd field values of the 4th database record >>> sf.records[3][1:3] ['060750601001', 4715]
Writing Shapefiles
>>> import shapefile>>> # Make a point shapefile >>> w = shapefile.Writer(shapefile.POINT) >>> w.point(90.3, 30) >>> w.point(92, 40) >>> w.point(-122.4, 30) >>> w.point(-90, 35.1) >>> w.field('FIRST_FLD') >>> w.field('SECOND_FLD','C','40') >>> w.record('First','Point') >>> w.record('Second','Point') >>> w.record('Third','Point') >>> w.record('Fourth','Point') >>> w.save('shapefiles/test/point') >>> # Create a polygon shapefile >>> w = shapefile.Writer(shapefile.POLYGON) >>> w.poly(parts=[[[1,5],[5,5],[5,1],[3,3],[1,1]]]) >>> w.field('FIRST_FLD','C','40') >>> w.field('SECOND_FLD','C','40') >>> w.record('First','Polygon') >>> w.save('shapefiles/test/polygon')
Tests
The above examples are very simple. There are many other shortcuts/features listed in the longer usage examples.
The file "README.txt" is a doctest module containing basic tests which are called if you run "shapefile.py" directly instead of importing it.
The sample shapefile used in the test and the output directory structure are in the "shapefiles" directory in the trunk.
Bug reports, fixes, improvements are welcome
Commercial support is available from NVision Solutions
