Gunakan Fiona dari Sean Gillies, pembungkus yang sangat sederhana dari perpustakaan OGR ( The Fiona User Manual )
Semua elemen dari shapefile (skema, catatan) diproses menggunakan kamus Python:
skema salah satu shapefile saya sebagai contoh:
{'geometry': 'LineString', 'properties': {u'faille': 'str:20', u'type': 'str:20', u'id': 'int'}}
satu catatan di shapefile:
{'geometry': {'type': 'LineString', 'coordinates': [(269884.20917418826, 151805.1917153612), (270409.89083992655, 153146.21637285672), (272298.05355768028, 154047.38494269375), (272941.74539327814, 155484.96337552898), (272169.31519056071, 156117.92701386689)]}, 'id': '1', 'properties': {'faille': u'de Salinas', 'type': u'normale'}}
jadi untuk menduplikasi sebuah shapefile:
from shapely.geometry import mapping, shape
import fiona
# Read the original Shapefile
with fiona.collection('original.shp', 'r') as input:
# The output has the same schema
schema = input.schema.copy()
# write a new shapefile
with fiona.collection(''file1.shp', 'w', 'ESRI Shapefile', schema) as output:
for elem in input:
output.write({'properties': elem['properties'],'geometry': mapping(shape(elem['geometry']))})
Jika Anda ingin membagi shapefile besar menjadi banyak yang lebih kecil, semuanya terjadi di loop for tetapi semua skema dari shapefile asli disimpan dalam kamus dengan schema = input.schema.copy()
dan{'properties': elem['properties']
lihat Bagaimana cara menemukan bantalan garis vektor di QGIS atau GRASS? untuk contoh
- membagi shapefile
- pertahankan atribut dari shapefile asli dalam shapefile yang terpecah
- dan tambahkan bidang baru di shapefile yang terpecah
Untuk pengguna Mac OS X atau Linux, mudah untuk menginstal. Untuk pengguna Windows, gunakan versi Christoph Gohlke Binari Windows Tidak Resmi untuk Paket Ekstensi Python
Bagi mereka yang kesulitan untuk menginstal perpustakaan di atas:
sumber