Orijinal set:
Pseudo-copy (TOC içinde CNTRL-drag) oluşturun ve klonla mekansal bire birleştirme yapın. Bu durumda 500m mesafe kullandım. Çıktı tablosu:
PAR_ID = PAR_ID_1 - kolay olan bu tablodan kayıtları kaldırın.
Tabloda yineleme yapın ve üzerindeki kayıtları (PAR_ID, PAR_ID_1) = (PAR_ID_1, PAR_ID) burada kaldırın. O kadar kolay değil, acrpy kullanın.
Havza sentroidlerini hesaplayın (UniqID = PAR_ID). Bunlar düğüm veya ağdır. Uzamsal birleştirme tablosu kullanarak bunları çizgilerle bağlayın. Bu, kesinlikle bu forumda bir yerde ele alınan ayrı bir konudur.
Aşağıdaki komut dosyası düğüm tablosunun şöyle göründüğünü varsayar:
MUID'nin parsellerden geldiği yerde, P2013 özetlenecek bir alandır. Bu durumda sadece sayım için = 1. [rcvnode] - tanımlanan grup / kümedeki ilk düğümün grup kimliğine eşit NODEREC'i depolayan komut dosyası çıktısı.
Önemli tabloların vurgulandığı bağlantılar tablosu yapısı
Times, bağlantı / kenar ağırlığını, yani düğümden düğüme seyahat maliyetini depolar. Bu durumda 1'e eşittir, böylece tüm komşulara seyahat maliyeti aynıdır. [fi] ve [ti] bağlı düğümlerin sıralı sayısıdır. Bu tabloyu doldurmak için, bağlantı kurmak için düğümlerden ve düğümlere nasıl atayacağınız konusunda bu forumda arama yapın.
Kendi tezgah mxd için özelleştirilmiş komut dosyası. Alanlar ve kaynaklar isimlendirilmeniz için değiştirilmeli, kodlanmalıdır:
import arcpy, traceback, os, sys,time
import itertools as itt
scriptsPath=os.path.dirname(os.path.realpath(__file__))
os.chdir(scriptsPath)
import COMMON
sys.path.append(r'C:\Users\felix_pertziger\AppData\Roaming\Python\Python27\site-packages')
import networkx as nx
RATIO = int(arcpy.GetParameterAsText(0))
try:
def showPyMessage():
arcpy.AddMessage(str(time.ctime()) + " - " + message)
mxd = arcpy.mapping.MapDocument("CURRENT")
theT=COMMON.getTable(mxd)
DÜĞÜM KATMANI BUL
theNodesLayer = COMMON.getInfoFromTable(theT,1)
theNodesLayer = COMMON.isLayerExist(mxd,theNodesLayer)
LAYKS LAYER ALIN
theLinksLayer = COMMON.getInfoFromTable(theT,9)
theLinksLayer = COMMON.isLayerExist(mxd,theLinksLayer)
arcpy.SelectLayerByAttribute_management(theLinksLayer, "CLEAR_SELECTION")
linksFromI=COMMON.getInfoFromTable(theT,14)
linksToI=COMMON.getInfoFromTable(theT,13)
G=nx.Graph()
arcpy.AddMessage("Adding links to graph")
with arcpy.da.SearchCursor(theLinksLayer, (linksFromI,linksToI,"Times")) as cursor:
for row in cursor:
(f,t,c)=row
G.add_edge(f,t,weight=c)
del row, cursor
pops=[]
pops=arcpy.da.TableToNumPyArray(theNodesLayer,("P2013"))
length0=nx.all_pairs_shortest_path_length(G)
nNodes=len(pops)
aBmNodes=[]
aBig=xrange(nNodes)
host=[-1]*nNodes
while True:
RATIO+=-1
if RATIO==0:
break
aBig = filter(lambda x: x not in aBmNodes, aBig)
p=itt.combinations(aBig, 2)
pMin=1000000
small=[]
for a in p:
S0,S1=0,0
for i in aBig:
p=pops[i][0]
p0=length0[a[0]][i]
p1=length0[a[1]][i]
if p0<p1:
S0+=p
else:
S1+=p
if S0!=0 and S1!=0:
sMin=min(S0,S1)
sMax=max(S0,S1)
df=abs(float(sMax)/sMin-RATIO)
if df<pMin:
pMin=df
aBest=a[:]
arcpy.AddMessage('%s %i %i' %(aBest,sMax,sMin))
if df<0.005:
break
lSmall,lBig,S0,S1=[],[],0,0
arcpy.AddMessage ('Ratio %i' %RATIO)
for i in aBig:
p0=length0[aBest[0]][i]
p1=length0[aBest[1]][i]
if p0<p1:
lSmall.append(i)
S0+=p0
else:
lBig.append(i)
S1+=p1
if S0<S1:
aBmNodes=lSmall[:]
for i in aBmNodes:
host[i]=aBest[0]
for i in lBig:
host[i]=aBest[1]
else:
aBmNodes=lBig[:]
for i in aBmNodes:
host[i]=aBest[1]
for i in lSmall:
host[i]=aBest[0]
with arcpy.da.UpdateCursor(theNodesLayer, "rcvnode") as cursor:
i=0
for row in cursor:
row[0]=host[i]
cursor.updateRow(row)
i+=1
del row, cursor
except:
message = "\n*** PYTHON ERRORS *** "; showPyMessage()
message = "Python Traceback Info: " + traceback.format_tb(sys.exc_info()[2])[0]; showPyMessage()
message = "Python Error Info: " + str(sys.exc_type)+ ": " + str(sys.exc_value) + "\n"; showPyMessage()
6 grup için çıktı örneği:
NETWORKX site paketine ihtiyacınız olacak
http://networkx.github.io/documentation/development/install.html
Komut dosyası gerekli sayıda kümeyi parametre olarak alır (yukarıdaki örnekte 6). Eşit ağırlık / seyahat kenarları mesafesine sahip bir grafik yapmak için düğümler ve bağlantı tabloları kullanır (Times = 1). Tüm düğümlerin kombinasyonunu 2 ile ele alır ve iki komşu grupta [P2013] toplamını hesaplar. Gerekli oran elde edildiğinde, örneğin ilk yinelemede (6-1) / 1, oran 1, yani 4 vb. İle 1'e kadar devam eder. Başlangıç noktaları çok önemlidir, bu yüzden 'uç' düğümlerinizin üstte oturduğundan emin olun (sıralama?) Örnek çıktıdaki ilk 3 gruba bakın. Bir sonraki yinelemede 'dal kesmeyi' önlemeye yardımcı olur.
MXD'den çalışmak için komut dosyası özelleştirme:
- ithalat COMMON gerekmez. Bu benim kendi şeyim, burada kendi ortam tabloumu okuyan, buradaNodesLayer, theLinksLayer, linksFromI, linksToI belirtildi. İlgili satırları kendi düğüm ve bağlantı katmanları adınızla değiştirin.
- P2013 alanının, kiracı sayısı veya parsel alanı gibi her şeyi saklayabileceğini unutmayın. Öyleyse çokgenleri yaklaşık eşit sayıda kişi vb. Alacak şekilde kümeleyebilirsiniz.
giriş şeklinin kesiştiği balık ağı 