ArcGIS olmadan Python üzerinden mevcut Shapefile'a nitelik alanı nasıl eklenir?

24

Eğer yoksa bir Shapefile'a nitelik alanı ekleyen bir Python betiğim var. Bu ArcGIS (grafiksel ya da Python üzerinden) ile yapmak kolaydır, ancak ArcGIS'e bağlı olmayan bir şey arıyorum.

Shapefile özellik içeriyor çünkü bu başarısız OGR ile denedim .

Pyshp'e baktım , ancak benzer şekilde, oluşturulduktan sonra şemayı değiştirmenin bir yolu yok. Shapefile (Python için) ile hiç denemedim , ancak bu özelliğin reklamını görmedim. Ayrıca bunun DBF dosyası ile dbfpy aracılığıyla düzeltilerek nasıl yapılacağını da göremiyorum .

Bir fikri olan var mı?

— Mike T
kaynak

Mevcut şekil dosyası yapısını klonlamak, yeni bir sütun eklemek ve ardından orijinal şekil dosyasını temel alarak doldurmak kabul edilebilir mi?

— DavidF

Özel Özellik özniteliklerini Python ile ESRI Shapefile Ekleme işleminin

— whuber

Bu soru, gis.stackexchange.com/q/3623/664'ün bir kopyası olarak kapatılmalıdır .

— whuber

Evet, aslında aynı. Baktım ama görmedim.

— Mike T

8

Zaten cevaplandığından beri bu sorulara bir göz atmalısınız: Python kullanarak Shapefile'a nasıl özel özellik eklenir?

/programming/4215658/adding-custom-feature-attributes-to-esri-shapefile-with-python

Sonuç olarak, yalnızca bir şekil dosyası, giriş dosyalarınızı komut dosyanızın sonunda silin.

— simo
kaynak

Aynı zamanda gis.stackexchange.com/questions/3623/… 'e benzer. Bunu kazdığın için teşekkürler

— Mike T

5

İşte bir çözüm. Geç olsun güç olmasın:

http://geospatialpython.com/2013/04/add-field-to-existing-shapefile.html?m=1

— GeospatialPython.com
kaynak

4

DBF olarak adlandırılan oldukça beyinsiz bir format sayesinde, mevcut özellik verileri ile şekil dosyalarına alan eklemek, yeniden yazma veya DBF'ye dolgu eklemeden mümkün değildir. Hazır bir çözümü bilmiyorum, fakat yapacağım şey var olanı temel alan yeni bir şekil dosyası oluşturmak için bir komut dosyası yazmak ve fazladan alanları yeni şekil dosyasına eklemek. Ardından geometri / özellik verilerini eskiden yeni şekil dosyasına kopyalayın. Ve son adım olarak, eski şekil dosyasını kaldırın ve yenisini yeniden adlandırın. Tüm bunlar, OGR piton bağları kullanılarak oldukça kolay bir şekilde gerçekleştirilir.

Alternatif olarak, yukarıda yalnızca DBF dosyasıyla yapmak için dbfpy komutunu kullanabilirsiniz. Adımların sırası aynı kalır:

Özgünle aynı yapıya sahip yeni bir DBF oluşturun
Yeni DBF'de yeni özellik alanları oluşturun
Verileri orijinal DBF'den yeni DBF'ye kopyala
Eski DBF'yi kaldırın, yeni DBF'yi eski DBF olarak yeniden adlandırın

DBF'de yer alan özellik bilgilerine referans vermediklerinden, shapefile (.shp) dosyasının kendisinde veya diğer dosyalarda herhangi bir değişiklik yapmanız gerekmez. Bununla birlikte, eski ve yeni DBF'deki kayıtların sırasını tamamen aynı tutmanız gerekir.

— Sasa Ivetic
kaynak

3

DBFpy bunun için çalışmalıdır. Daha sonra bu sayfada örnek gördünüz mü:

http://dbfpy.sourceforge.net/

Kilitleme sırasında sorunlara yol açabileceğinden, şekil dosyasının ArcGIS dahil başka bir uygulama tarafından düzenlenmediğinden emin olun.

— Rob Clark
kaynak

DBF dosyasında zaten verileriniz varsa, bunun işe yaradığını sanmıyorum. Daha önce bakmıştım, ancak bir hata alıyorum: "En az bir kayıt eklendi, yapı değiştirilemez". Aklında belirli bir örnek var mı?

— Mike T

Ah, şimdi evet, Sasa'nın yeni bir DBF'nin oluşturulması gerekeceğini söylediğini hatırlıyorum. Şemayı (alanlardaki vb.) Kopyalayın, sonra eklemeler yapın, sonra kayıtları kopyalayın. "Harika" DBF ... :(

— Rob Clark

@Mike Tek yapmak istediğiniz alan bir alan eklemek olduğunda, bir kayıt nasıl eklendi ? Kayıt eklemek bir hatadır çünkü özellikler ve şekiller arasındaki bağlantıyı mahveder. Bir alan eklemek hiçbir şekilde zarar vermez. DBF dosyalarını düzenleyebilen herhangi bir kütüphane işi doğru şekilde yapar.

— whuber

@ whuber: Bu onların hata mesajı. Verileri olan mevcut bir dbf'yi açın ve şunu görün:from dbfpy import dbf; db = dbf.Dbf('my.dbf'); db.addField(("FOO", "C", 15))

— Mike T

@Mike Durumu açıklığa kavuşturduğun için teşekkürler. Dbfpy :-( 'de gereksiz bir sınırlamanın sonucu gibi geliyor. Neden olduğunu tahmin edebilirim: boş olmayan bir veritabanına bir alan eklemek tüm kayıtların fiziksel olarak okunmasını, genişletilmesini ve tekrar yazılmasını gerektirir. farklı bir dBase kütüphanesi bulmak veya başka bir yazılımı kullanmak için ;-).

— whuber

1

OGR kullanarak bir çözüm buldum ve önceki bir sorudan dolayı teşekkürler . İşte tam bir örnek:

from osgeo import ogr

# Open a Shapefile, and get field names
source = ogr.Open('my.shp', update=True)
layer = source.GetLayer()
layer_defn = layer.GetLayerDefn()
field_names = [layer_defn.GetFieldDefn(i).GetName() for i in range(layer_defn.GetFieldCount())]
print len(field_names), 'MYFLD' in field_names

# Add a new field
new_field = ogr.FieldDefn('MYFLD', ogr.OFTInteger)
layer.CreateField(new_field)

# Close the Shapefile
source = None

Benim sorunum, onun layer_defn.AddFieldDefn(new_field)yerine kullanmamdı layer.CreateField(new_field). Yardım için çok teşekkürler, ve benzer diğer soruları iyice kontrol etmediğim için üzgünüm.

— Mike T
kaynak