/usr/lib/python3/dist-packages/pyshp-1.2.1.egg-info/PKG-INFO is in python3-pyshp 1.2.1-1.
This file is owned by root:root, with mode 0o644.
The actual contents of the file can be viewed below.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 153 154 155 156 157 158 159 160 161 162 163 164 165 166 167 168 169 170 171 172 173 174 175 176 177 178 179 180 181 182 183 184 185 186 187 188 189 190 191 192 193 194 195 196 197 198 199 200 201 202 203 204 205 206 207 208 209 210 211 212 213 214 215 216 217 218 219 220 221 222 223 224 225 226 227 228 229 230 231 232 233 234 235 236 237 238 239 240 241 242 243 244 245 246 247 248 249 250 251 252 253 254 255 256 257 258 259 260 261 262 263 264 265 266 267 268 269 270 271 272 273 274 275 276 277 278 279 280 281 282 283 284 285 286 287 288 289 290 291 292 293 294 295 296 297 298 299 300 301 302 303 304 305 306 307 308 309 310 311 312 313 314 315 316 317 318 319 320 321 322 323 324 325 326 327 328 329 330 331 332 333 334 335 336 337 338 339 340 341 342 343 344 345 346 347 348 349 350 351 352 353 354 355 356 357 358 359 360 361 362 363 364 365 366 367 368 369 370 371 372 373 374 375 376 377 378 379 380 381 382 383 384 385 386 387 388 389 390 391 392 393 394 395 396 397 398 399 400 401 402 403 404 405 406 407 408 409 410 411 412 413 414 415 416 417 418 419 420 421 422 423 424 425 426 427 428 429 430 431 432 433 434 435 436 437 438 439 440 441 442 443 444 445 446 447 448 449 450 451 452 453 454 455 456 457 458 459 460 461 462 463 464 465 466 467 468 469 470 471 472 473 474 475 476 477 478 479 480 481 482 483 484 485 486 487 488 489 490 491 492 493 494 495 496 497 498 499 500 501 502 503 504 505 506 507 508 509 510 511 512 513 514 515 516 517 518 519 520 521 522 523 524 525 526 527 528 529 530 531 532 533 534 535 536 537 538 539 540 541 542 543 544 545 546 547 548 549 550 551 552 553 554 555 556 557 558 559 560 561 562 563 564 565 566 567 568 569 570 571 572 573 574 575 576 577 578 579 580 581 582 583 584 585 586 587 588 589 590 591 592 593 594 595 596 597 598 599 600 601 602 603 604 605 606 607 608 609 610 611 612 613 614 615 616 617 618 619 620 621 622 623 624 625 626 627 628 629 630 631 632 633 634 635 636 637 638 639 640 641 642 643 644 645 646 647 648 649 650 651 652 653 654 655 656 657 658 659 660 661 662 663 664 665 666 667 668 669 670 | Metadata-Version: 1.1
Name: pyshp
Version: 1.2.1
Summary: Pure Python read/write support for ESRI Shapefile format
Home-page: http://code.google.com/p/pyshp
Author: Joel Lawhead
Author-email: jlawhead@geospatialpython.com
License: MIT
Description: PyShp
========================
:Author: Joel Lawhead - jlawhead@geospatialpython.com
:Version 1.2.1
:Revised: May 11, 2014
.. contents::
Overview
--------
PyShp provides read and write support for the Esri
Shapefile format. The Shapefile format is a popular Geographic Information
System vector data format created by Esri. For more information about this format
please read the well-written "ESRI Shapefile Technical Description - July 1998"
located at http://www.esri.com/library/whitepapers/pdfs/shapefile.pdf.
The Esri document describes the shp and shx file formats. However a third file
format called dbf is also required. This format is documented on the web as the
"XBase File Format Description" and is a simple file-based database format created
in the 1960's. For more on this specification see:
http://www.clicketyclick.dk/databases/xbase/format/index.html
Both the Esri and XBase file-formats are very simple in design and
memory efficient which is part of the reason the shapefile format remains popular
despite the numerous ways to store and exchange GIS data available today.
Pyshp is compatible with Python 2.4-3.x.
This document provides examples for using PyShp to read and write shapefiles.
Currently the sample census blockgroup shapefile referenced in the examples is
only available on the google code project site at http://code.google.com/p/pyshp.
These examples are straight-forward and you can also easily run them against your
own shapefiles manually with minimal modification. Other examples for specific
topics are continually added to the pyshp wiki on google code and the blog
http://GeospatialPython.com.
Important: For information about map projections, shapefiles,
and Python please visit: http://code.google.com/p/pyshp/wiki/MapProjections
I sincerely hope this library eliminates the mundane distraction of simply
reading and writing data, and allows you to focus on the challenging and FUN
part of your geospatial project.
Examples
--------
Before doing anything you must import the library.
>>> import shapefile
The examples below will use a shapefile created from the U.S. Census Bureau
Blockgroups data set near San Francisco, CA and available in the subversion
repository of the pyshp google code site.
Reading Shapefiles
++++++++++++++++++
To read a shapefile create a new "Reader" object and pass it the name of an
existing shapefile. The shapefile format is acutally a collection of three
files. You specify the base filename of the shapefile or the complete filename
of any of the shapefile component files.
>>> sf = shapefile.Reader("shapefiles/blockgroups")
OR
>>> sf = shapefile.Reader("shapefiles/blockgroups.shp")
OR
>>> sf = shapefile.Reader("shapefiles/blockgroups.dbf")
OR any of the other 5+ formats which are potentially part of a shapefile.
The library does not care about extensions.
Reading Shapefiles from File-Like Objects
.........................................
You can also load shapefiles from any Python file-like object using keyword
arguments to specify any of the three files. This feature is very powerful
and allows you to load shapefiles from a url, from a zip file, serialized
object, or in some cases a database.
>>> myshp = open("shapefiles/blockgroups.shp", "rb")
>>> mydbf = open("shapefiles/blockgroups.dbf", "rb")
>>> r = shapefile.Reader(shp=myshp, dbf=mydbf)
Notice in the examples above the shx file is never used. The shx file is a
very simple fixed-record index for the variable length records in the shp file.
This file is optional for reading. If it's available pyshp will use the shx file
to access shape records a little faster but will do just fine without it.
Reading Geometry
................
A shapefile's geometry is the collection of points or shapes made from verticies
and implied arcs representing physical locations. All types of shapefiles
just store points. The metadata about the points determine how they are handled by
software.
You can get the a list of the shapefile's geometry by calling the shapes()
method.
>>> shapes = sf.shapes()
The shapes method returns a list of Shape objects describing the
geometry of each shape record.
>>> len(shapes)
663
You can iterate through the shapefile's geometry using the iterShapes() method.
>>> len(list(sf.iterShapes()))
663
Each shape record contains the following attributes:
>>> for name in dir(shapes[3]):
... if not name.startswith('__'):
... name
'bbox'
'parts'
'points'
'shapeType'
- shapeType: an integer representing the type of shape as defined by the
shapefile specification.
>>> shapes[3].shapeType
5
- bbox: If the shape type contains multiple points this tuple describes the
lower left (x,y) coordinate and upper right corner coordinate creating a
complete box around the points. If the shapeType is a Null
(shapeType == 0) then an AttributeError is raised.
>>> # Get the bounding box of the 4th shape.
>>> # Round coordinates to 3 decimal places
>>> bbox = shapes[3].bbox
>>> ['%.3f' % coord for coord in bbox]
['-122.486', '37.787', '-122.446', '37.811']
- parts: Parts simply group collections of points into shapes. If the shape record
has multiple parts this attribute contains the index of the first point of each part.
If there is only one part then a list containing 0 is returned.
>>> shapes[3].parts
[0]
- points: The points attribute contains a list of tuples containing an (x,y)
coordinate for each point in the shape.
>>> len(shapes[3].points)
173
>>> # Get the 8th point of the fourth shape
>>> # Truncate coordinates to 3 decimal places
>>> shape = shapes[3].points[7]
>>> ['%.3f' % coord for coord in shape]
['-122.471', '37.787']
To read a single shape by calling its index use the shape() method. The index
is the shape's count from 0. So to read the 8th shape record you would
use its index which is 7.
>>> s = sf.shape(7)
>>> # Read the bbox of the 8th shape to verify
>>> # Round coordinates to 3 decimal places
>>> ['%.3f' % coord for coord in s.bbox]
['-122.450', '37.801', '-122.442', '37.808']
Reading Records
................
A record in a shapefile contains the attributes for each shape in the
collection of geometry. Records are stored in the dbf file. The link
between geometry and attributes is the foundation of Geographic Information
Systems. This critical link is implied by the order of shapes and
corresponding records in the shp geometry file and the dbf attribute file.
The field names of a shapefile are available as soon as you read a shapefile.
You can call the "fields" attribute of the shapefile as a Python list. Each
field is a Python list with the following information:
- Field name: the name describing the data at this column index.
- Field type: the type of data at this column index. Types can be: Character, Numbers, Longs, Dates, or Memo.
The "Memo" type has no meaning within a GIS and is part of the xbase spec instead.
- Field length: the length of the data found at this column index. Older GIS software may truncate this
length to 8 or 11 characters for "Character" fields.
- Decimal length: the number of decimal places found in "Number" fields.
To see the fields for the Reader object above (sf) call the "fields" attribute:
>>> fields = sf.fields
>>> assert 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], ["WHITE", "N", 9, 0],
... ["BLACK", "N", 8, 0], ["AMERI_ES", "N", 7, 0], ["ASIAN_PI", "N", 8, 0],
... ["OTHER", "N", 8, 0], ["HISPANIC", "N", 8, 0], ["AGE_UNDER5", "N", 8, 0],
... ["AGE_5_17", "N", 8, 0], ["AGE_18_29", "N", 8, 0], ["AGE_30_49", "N", 8, 0],
... ["AGE_50_64", "N", 8, 0], ["AGE_65_UP", "N", 8, 0],
... ["NEVERMARRY", "N", 8, 0], ["MARRIED", "N", 9, 0], ["SEPARATED", "N", 7, 0],
... ["WIDOWED", "N", 8, 0], ["DIVORCED", "N", 8, 0], ["HSEHLD_1_M", "N", 8, 0],
... ["HSEHLD_1_F", "N", 8, 0], ["MARHH_CHD", "N", 8, 0],
... ["MARHH_NO_C", "N", 8, 0], ["MHH_CHILD", "N", 7, 0],
... ["FHH_CHILD", "N", 7, 0], ["HSE_UNITS", "N", 9, 0], ["VACANT", "N", 7, 0],
... ["OWNER_OCC", "N", 8, 0], ["RENTER_OCC", "N", 8, 0],
... ["MEDIAN_VAL", "N", 7, 0], ["MEDIANRENT", "N", 4, 0],
... ["UNITS_1DET", "N", 8, 0], ["UNITS_1ATT", "N", 7, 0], ["UNITS2", "N", 7, 0],
... ["UNITS3_9", "N", 8, 0], ["UNITS10_49", "N", 8, 0],
... ["UNITS50_UP", "N", 8, 0], ["MOBILEHOME", "N", 7, 0]]
You can get a list of the shapefile's records by calling the records() method:
>>> records = sf.records()
>>> len(records)
663
Similar to the geometry methods, you can iterate through dbf records using the
recordsIter() method.
>>> len(list(sf.iterRecords()))
663
Each record is a list containing an attribute corresponding to each field in the
field list.
For example in the 4th record of the blockgroups shapefile the 2nd and 3rd
fields are the blockgroup id and the 1990 population count of
that San Francisco blockgroup:
>>> records[3][1:3]
['060750601001', 4715]
To read a single record call the record() method with the record's index:
>>> rec = sf.record(3)
>>> rec[1:3]
['060750601001', 4715]
Reading Geometry and Records Simultaneously
...........................................
You way want to examine both the geometry and the attributes for a record at the
same time. The shapeRecord() and shapeRecords() method let you do just that.
Calling the shapeRecords() method will return the geometry and attributes for
all shapes as a list of ShapeRecord objects. Each ShapeRecord instance has a
"shape" and "record" attribute. The shape attribute is a ShapeRecord object as
dicussed in the first section "Reading Geometry". The record attribute is a
list of field values as demonstrated in the "Reading Records" section.
>>> shapeRecs = sf.shapeRecords()
Let's read the blockgroup key and the population for the 4th blockgroup:
>>> shapeRecs[3].record[1:3]
['060750601001', 4715]
Now let's read the first two points for that same record:
>>> points = shapeRecs[3].shape.points[0:2]
>>> len(points)
2
The shapeRec() method reads a single shape/record pair at the specified index.
To get the 4th shape record from the blockgroups shapfile use the third index:
>>> shapeRec = sf.shapeRecord(3)
The blockgroup key and population count:
>>> shapeRec.record[1:3]
['060750601001', 4715]
>>> points = shapeRec.shape.points[0:2]
>>> len(points)
2
Writing Shapefiles
++++++++++++++++++
PyShp tries to be as flexible as possible when writing shapefiles while
maintaining some degree of automatic validation to make sure you don't
accidentally write an invalid file.
PyShp can write just one of the component files such as the shp or dbf file
without writing the others. So in addition to being a complete
shapefile library, it can also be used as a basic dbf (xbase) library. Dbf files are
a common database format which are often useful as a standalone simple
database format. And even shp files occasionaly have uses as a standalone
format. Some web-based GIS systems use an user-uploaded shp file to specify
an area of interest. Many precision agriculture chemical field sprayers also
use the shp format as a control file for the sprayer system (usually in
combination with custom database file formats).
To create a shapefile you add geometry and/or attributes using methods in the
Writer class until you are ready to save the file.
Create an instance of the Writer class to begin creating a shapefile:
>>> w = shapefile.Writer()
Setting the Shape Type
......................
The shape type defines the type of geometry contained in the shapefile. All of
the shapes must match the shape type setting.
Shape types are represented by numbers between 0 and 31 as defined by the
shapefile specification. It is important to note that numbering system has
several reserved numbers which have not been used yet therefore the numbers of
the existing shape types are not sequential.
You can reference shape types by the numbers or by constants defined by PyShp:
shapefile.NULL = 0
shapefile.POINT = 1
shapefile.POLYLINE = 3
shapefile.POLYGON = 5
shapefile.MULTIPOINT = 8
shapefile.POINTZ = 11
shapefile.POLYLINEZ = 13
shapefile.POLYGONZ = 15
shapefile.MULTIPOINTZ = 18
shapefile.POINTM = 21
shapefile.POLYLINEM = 23
shapefile.POLYGONM = 25
shapefile.MULTIPOINTM = 28
shapefile.MULTIPATCH = 31
There are three ways to set the shape type:
- Set it when creating the class instance.
- Set it by assigning a value to an existing class instance.
- Set it automatically to the type of the first shape by saving the shapefile.
To manually set the shape type for a Writer object when creating the Writer:
>>> w = shapefile.Writer(shapeType=1)
or we can use the constants as explained above:
>>> w = shapefile.Writer(shapefile.POINT)
As you can see, specifying the shapeType argument explicitly isn't necessary.
>>> w.shapeType
1
OR you can set it after the Writer is created by changing the property:
>>> w.shapeType = 3
>>> w.shapeType
3
Geometry and Record Balancing
.............................
Because every shape must have a corresponding record it is critical that the
number of records equals the number of shapes to create a valid shapefile. To
help prevent accidental misalignment PyShp has an "auto balance" feature to
make sure when you add either a shape or a record the two sides of the
equation line up. This feature is NOT turned on by default. To activate it
set the attribute autoBalance to 1 (True):
>>> w.autoBalance = 1
You also have the option of manually calling the balance() method each time you
add a shape or a record to ensure the other side is up to date. When balancing
is used null shapes are created on the geometry side or a record with a value of
"NULL" for each field is created on the attribute side.
The balancing option gives you flexibility in how you build the shapefile.
Without auto balancing you can add geometry or records at anytime. You can
create all of the shapes and then create all of the records or vice versa. You
can use the balance method after creating a shape or record each time and make
updates later. If you do not use the balance method and forget to manually
balance the geometry and attributes the shapefile will be viewed as corrupt by
most shapefile software.
With auto balanacing you can add either shapes or geometry and update blank
entries on either side as needed. Even if you forget to update an entry the
shapefile will still be valid and handled correctly by most shapefile software.
Adding Geometry
...............
Geometry is added using one of three methods: "null", "point", or "poly". The "null"
method is used for null shapes, "point" is used for point shapes, and "poly" is
used for everything else.
**Adding a Null shape**
Because Null shape types (shape type 0) have no geometry the "null" method is
called without any arguments.
>>> w = shapefile.Writer()
>>> w.null()
The writer object's shapes list will now have one null shape:
>>> assert w.shapes()[0].shapeType == shapefile.NULL
**Adding a Point shape**
Point shapes are added using the "point" method. A point is specified by an
x, y, and optional z (elevation) and m (measure) value.
>>> w = shapefile.Writer(shapefile.POINT)
>>> w.point(122, 37) # No elevation or measure values
>>> w.shapes()[0].points
[[122, 37, 0, 0]]
>>> w.point(118, 36, 4, 8)
>>> w.shapes()[1].points
[[118, 36, 4, 8]]
**Adding a Poly Shape**
"Poly" shapes can be either polygons or lines. Shapefile polygons must have at
least 4 points and the last point must be the same as the first. PyShp automatically
enforces closed polygons.
A line must have at least two points.
Because of the similarities between these two shape types they are created using
a single method called "poly".
>>> w = shapefile.Writer(shapefile.POLYGON)
>>> w.poly(shapeType=3, parts=[[[122,37,4,9], [117,36,3,4]], [[115,32,8,8],
... [118,20,6,4], [113,24]]])
**Adding a Polygon with Rings**
Polygons consist of rings which mean they are closed. The first point and last point
of a ring must be the same. PyShp enforces ring closure if the ring is incomplete when
you add the shape. Polygons can have inner rings which create holes. Holes are defined
by the order of the points. Normally points in a ring run clockwise. If the points
run counter-clockwise then they form a hole. If you don't order the points correctly
you'll just have overlapping polygons.
>>> w = shapefile.Writer(shapefile.POLYGON)
>>> outer_ring = [[10,10],[50,50],[100,10],[50,-50],[10,10]]
>>> inner_ring = [[40,10],[50,30],[70,10],[50,-30],[40,10]]
>>> inner_ring.reverse()
You can use the "shapefile.signed_area()" method to determine if a ring is clockwise
or counter-clockwise. A value >= 0 means the ring is counter-clockwise and < 0 means
the ring is clockwise. The value returned is also the area of the polygon.
>>> # Clockwise ring
... shapefile.signed_area(outer_ring)
-4500.0
>>> # Counter-clockwise ring
... shapefile.signed_area(inner_ring)
900.0
**Creating 3D Polygons**
Elevation values, known as "Z" values allow you to create 3-dimensional shapefiles. The
z value is an extra value specified as part of a point.
>>> w = shapefile.Writer(shapeType=shapefile.POLYGONZ)
>>> w.poly([[[-89.0, 33, 12], [-90, 31, 11], [-91, 30, 12]]], shapeType=15)
>>> w.field("NAME")
>>> w.record("PolyZTest")
>>> w.save("shapefiles/test/MyPolyZ")
The z values are stored in a seperate shape attribute.
>>> r = shapefile.Reader("shapefiles/test/MyPolyZ")
>>> s = r.shape(0)
>>> s.points
[[-89.0, 33.0], [-90.0, 31.0], [-91.0, 30.0], [-89.0, 33.0]]
>>> s.z
[12.0, 11.0, 12.0, 12.0]
Creating Attributes
...................
Creating attributes involves two steps. Step 1 is to create fields to contain
attribute values and step 2 is to populate the fields with values for each
shape record.
The following attempts to create a complete shapefile:
>>> w = shapefile.Writer(shapefile.POINT)
>>> w.point(1,1)
>>> w.point(3,1)
>>> w.point(4,3)
>>> w.point(2,2)
>>> 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')
>>> 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')
>>> w = shapefile.Writer(shapefile.POLYLINE)
>>> w.line(parts=[[[1,5],[5,5],[5,1],[3,3],[1,1]]])
>>> w.poly(parts=[[[1,3],[5,3]]], shapeType=shapefile.POLYLINE)
>>> w.field('FIRST_FLD','C','40')
>>> w.field('SECOND_FLD','C','40')
>>> w.record('First','Line')
>>> w.record('Second','Line')
>>> w.save('shapefiles/test/line')
You can also add attributes using keyword arguments where the keys are field names.
>>> w = shapefile.Writer(shapefile.POLYLINE)
>>> w.line(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_FLD='First', SECOND_FLD='Line')
>>> w.save('shapefiles/test/line')
File Names
..........
File extensions are optional when reading or writing shapfiles. If you specify them Pyshp
ignores them anyway. When you save files you can specify a base file name that is used for
all three file types. Or you can specify a nmae for one or more file types. In that case,
any file types not assigned will not save and only file types with file names will be saved.
If you do not specify any file names (i.e. save()), then a unique file name is generated with
the prefix "shapefile\_" followed by random characters which is used for all three files. The
unique file name is returned as a string.
>>> targetName = w.save()
>>> assert("shapefile_" in targetName)
Saving to File-Like Objects
...........................
Just as you can read shapefiles from python file-like objects you can also write them.
>>> try:
... from StringIO import StringIO
... except ImportError:
... from io import BytesIO as StringIO
>>> shp = StringIO()
>>> shx = StringIO()
>>> dbf = StringIO()
>>> w.saveShp(shp)
>>> w.saveShx(shx)
>>> w.saveDbf(dbf)
>>> # Normally you would call the "StringIO.getvalue()" method on these objects.
>>> shp = shx = dbf = None
Editing Shapefiles
++++++++++++++++++
The Editor class attempts to make changing existing shapefiles easier by handling the reading and writing details behind the scenes.
Let's add shapes to existing shapefiles:
Add a point to a point shapefile
>>> e = shapefile.Editor(shapefile="shapefiles/test/point.shp")
>>> e.point(0,0,10,2)
>>> e.record("Appended","Point")
>>> # We added z and m values so
>>> # change the shapetype
>>> e.shapeType = shapefile.POINTZ
>>> e.save('shapefiles/test/point')
Edit the appended point to change the "y" and "z" value
>>> e = shapefile.Editor(shapefile="shapefiles/test/point.shp")
>>> # Find the point by the attribute
>>> for s in enumerate(e.records):
... i, record = s
... if record[0] == "Appended":
... geom = e._shapes[i]
... # Change the y value to 5
... geom.points[0][1] = 5
... # Change the z value to 9
... if hasattr(geom, "z"):
... geom.z = (9,)
... else:
... geom.points[0][2] = 9
>>> e.save('shapefiles/test/point')
Add a new line to a line shapefile:
>>> e = shapefile.Editor(shapefile="shapefiles/test/line.shp")
>>> e.line(parts=[[[10,5],[15,5],[15,1],[13,3],[11,1]]])
>>> e.record('Appended','Line')
>>> e.save('shapefiles/test/line')
>>> e = None
Add a new polygon to a polygon shapefile:
>>> e = shapefile.Editor(shapefile="shapefiles/test/polygon.shp")
>>> e.poly(parts=[[[5.1,5],[9.9,5],[9.9,1],[7.5,3],[5.1,1]]])
>>> e.record("Appended","Polygon")
>>> e.save('shapefiles/test/polygon')
>>> e = None
Remove the first point in each shapefile - for a point shapefile that is
the first shape and record"
>>> e = shapefile.Editor(shapefile="shapefiles/test/point.shp")
>>> e.delete(0)
>>> e.save('shapefiles/test/point')
>>> e = None
Remove the last shape in the polygon shapefile.
>>> e = shapefile.Editor(shapefile="shapefiles/test/polygon.shp")
>>> e.delete(-1)
>>> e.save('shapefiles/test/polygon')
>>> e = None
Python __geo_interface__
++++++++++++++++++++++++
The Python __geo_interface__ convention provides a data interchange interface
among geospatial Python libraries. The interface returns data as GeoJSON.
More information on the __geo_interface__ protocol can be found at:
https://gist.github.com/sgillies/2217756.
More information on GeoJSON is available at http://geojson.org http://geojson.org.
>>> s = sf.shape(0)
>>> s.__geo_interface__["type"]
'MultiPolygon'
Keywords: gis geospatial geographic shapefile shapefiles
Platform: UNKNOWN
Classifier: Programming Language :: Python
Classifier: Topic :: Scientific/Engineering :: GIS
Classifier: Topic :: Software Development :: Libraries
Classifier: Topic :: Software Development :: Libraries :: Python Modules
|