Tüm ülkeye yayılmış 130 milyondan fazla kaydı bulunan müşterilere ait konum verilerim var. Bu PostGIS'te Mekansal Bir Tablodur

Şimdi tüm ülkeyi 'hizmet alanlarına' bölmek istiyorum, böylece her hizmet alanı:

• En az 30 Müşteri
• En fazla 1000 km2'dir.

ArcGIS'in Toplu aracını araştırdım, ancak kriterlerimi karşılamıyor, çünkü:

• Minimum sayımı dikkate almaz
• Elimde sadece Qgis & Grass var.

Bunu başarmak için hangi aracı / süreci kullanabilirim?

Küçük bir deneme yaptım:

• 900 rastgele puan üretti
• onları minimum yayılan ağaç kullanarak bağladılar
• ortasında bir yerde lavabo aldı

• yönlü grafik oluşturmak için bağlantıları çevirdi, yani her bir düğüm sonunda aynı lavaboya akıyor

• 30 ortalama "boyutu" almaya çalışırken büyük olanlar içine küçük havza toplamak için kullandığım komut dosyası koştu.

SONUÇ:

İlginizi çekiyorsa bana bildirin, betiği kazabilirim (zaten bu sitede bir yerde yayınlanmış olabilir) veya iş akışını daha ayrıntılı olarak açıklayabilirim. ArcGIS kullandım.

-------------------------------------

** GÜNCELLEME Eylül 2017

Yukarıdakine çok benzer, ancak raster yaklaşımı, ArcGIS'ten Mesafe ve Hidroloji araçları kullanılır. Zor kısmı, minimum yayılan ağacın oluşturulması ve raster formatında 'yakalama' çıkışlarının bulunmasıdır.

İzlediğim adımlar şunlardır:

• Düğümleri üçgenleyin (kümeye işaret eder), üçgen ağ kenarlarını ayıklayın, düğümleri AĞIRLIK rasterine dönüştürün :

• Kenarlar tablosundaki hesaplama alanı: (ShapeLength ^ 3 / 1e6) kullanıyorum. Raster dönüştürmek, COST raster oluşturmak için çok yüksek değer boşlukları doldurun . Bu, birbirine yakın noktalar arasındaki akışı teşvik eder. Akış yollarının lavaboya doğru giden düz çizgilere yakın olmak yerine minimum yayılma ağacına (MST) benzeyeceği umulmaktadır.
• Herhangi bir düğümü (OUTLET / SINK) seçin ve kaynak olarak COST yüzeyini ve SINK kullanarak Cost Back Link raster'ı oluşturun. Geri bağlantı raster kullanarak iyi Akış Yön raster dönüştürün Int(Power(2,"backlink"-1)). Akış yönü ve ağırlık rasterini kullanarak akışı biriktirin. Maliyet atama hilesi gerçekten MST'ye benzer bir şey üretir:

Puanları 50'ye göre gruplamaya karar verdim. Ağın fraktal yapısını akılda tutarak, sınırı 'Yüksek Akış Akımları' başlangıcında, yani Con ("FlowAccum"> 45,1) başında biraz daha düşük INLET puanları belirledim. Girişler, hücreler olarak tanımlanır.

Yüksek Akış var ve odak istatistikleri = 2 ve hücre bir EVYE değil (Akış Yönünde veri değeri yok) :

Havzaları tanımlamak için girişleri akma noktaları ve akış yönü olarak kullanın. Resimde elde edilen 115 havza gösterilmektedir:

İstatistikleri: ortalama = 50.33, min = 46 ve maks = 74.

İkinci havza setini elde etmek için, önceden tanımlanmış havzalarda AĞIRLIK rasterini silmek, Akış Toplama'yı hesaplamak vb.

Yöntem, raster tabanlı olduğu için milyonlarca puan için işe yarayacak, bu sayıların üçgenlenmesi de bir sorun olmayacak.