ArcGIS Yol Mesafesi aracıyla maliyet yüzeylerini kullanarak hareketi mi analiz ediyorsunuz?


9

Bir manzaradaki varsayımsal hareketi (yaya olarak) enerji harcamalarına dayalı olarak analiz etmek istiyorum, ancak bana yardımcı olabileceğinizi umduğum bir sorunla karşılaştım. Bunu, oluşturduğum maliyet yüzeylerini kullanarak Mekansal Analist'teki ArcGIS'in Yol Mesafesi aracını kullanarak yapmaya çalıştım, ancak sonuç beklediğim gibi değil.

Yükseklik yüzeyim böyle görünüyor (ASTER GDEM'den indirildi): Beyaz / purpe alanları en yüksek ve yeşil en düşük olmak üzere yükseklik verileri.

Yükseklik verilerine dayanarak, harita birimi (m) başına enerji harcaması (Watt cinsinden metabolik hız) içermesi gereken bir maliyet yüzeyi oluşturdum. Bunun için bu formülü kullandım: M = 1.5W + 2.0 (W + L) (L / W)2 + N (W + L) (1.5V2 + 0.35V * abs(G + 6))

Veya Raster Hesap Makinesi terimlerini koyun: (1.5 * 60) + (2.0 * (60 + 3) * Square((3 / 60))) + (1.2 * (60 + 3) * (Square((1.5 * "movementspeed")) + (0.35 * "movementspeed") * Abs(("slopeinpercent" + 6))))

M, Watt cinsinden metabolik hız olduğunda, W, modellenen bireyin ağırlığı, L, bireyin taşınan yüküdür, N, arazideki hareket kolaylığını (1.2'ye ayarlanan test amaçları için) tanımlayan bir faktördür, V, bireyin hareket hızı ve G yüzde olarak eğimdir. Bu, 90 ve 25000 arasında değerlere sahip bir yüzey yarattı, değerlerin çoğu 90 ile 1000 arasındadır (bu doğru görünüyor, saçma yüksek değerler büyük olasılıkla kolayca sabitlenebilen kusurlu eğim değerlerinin sonucudur).

Hareket hızı şu formül kullanılarak hesaplanmıştır: V = 6e^(-3.5 * |s + 0.05|burada s, derece cinsinden eğimdir.

Ya da Raster Hesaplayıcı terimlerini koyun: 6 * Exp( - 3.5 * Abs(Tan("slopeindegrees") + 0.05)) Bu, doğru görünen ve beklediğim ile tutarlı olan 0 ile 5.9 km / s arasında bir değere sahip bir yüzey yarattı.

Şimdi, bu yüzeyler Yol Mesafesi aracında girdi olarak kullanılmıştır; girdi yüzey raster olarak DEM (yani yüzey_raster), maliyet raster olarak enerji harcaması olan yüzey ve modellenmiş bireyin bir eğimde yukarı veya aşağı hareket edip etmediğini hesaplamak için aletin izin vermesi için dikey raster olarak DEM. Test amacıyla, DEM'in kuzey-batı ve güney-doğu köşelerinde iki nokta kaynak veri olarak kullanılmıştır (yani, in_source_data). Çıktı şuydu (kırmızı, istemeden en düşük değerler ve en yüksek mavi): Kırmızının en düşük değerler olduğu enerji harcaması

Çıktıya ilişkin yorumum, yükseklik farklarını hemen hemen görmezden gelmesi ve değerdeki farklılıkların sadece mesafe farklarıyla ilişkili olmasıdır. Yüzeyin bölgenin batı kesimindeki düz alanları takip etmesini beklerdim ve açıkça yapmadığı dağlık doğu kısımlarından kaçınırdım. Ancak, hala bu tür analizlerde yeniyim ve başkalarının yorumlarını takdir ediyorum. Öyleyse, metodolojim / formüllerimde garip çıktıya neden olabilecek kusurları gösterebilen var mı? Ya da çıktı bekleniyor ve ben sadece bir yol mesafe analizinden beklediğim yanlış anlıyorum?


"Yol mesafesi aracı" alanında "Giriş raster veya özellik kaynak verileri" alanı, bir yolu hesapladığınız BİR veya daha fazla noktayı temsil eder. ESRI şöyle diyor: "Bu, her çıkış hücresi konumu için en az birikmiş maliyet mesafesinin hesaplandığı hücreleri veya konumları tanımlayan bir tarama veya özellik veri kümesidir. Bunun için bir DEM kullanmak, IMO anlamına gelmez. "En Kısa Yol" hesaplaması (Araçlar Arc10'da karıştırılır ve farklı adlandırılır), böylece daha önce kaynak konumları tanımlayan bir konuma NEREDEN bir yol hesaplarsınız
G-sihirbazı

DEM'i kaynak katman olarak kullanmadım, sadece ARCGIS'in "in_surface_raster" olarak adlandırdığı şeyde girdi olarak kullanıldı. Kuzeybatı ve güneydoğuda kaynak katmanı olarak iki nokta kullandım. Karışıklık için özür dilerim, ikisi arasında düzgün bir şekilde ayrım yapmak için yazımı düzenleyeceğim.
Oulah

Yanıtlar:


4

Bu, bir dem, dikey raster ve dikey faktör spesifikasyonu içeren yol mesafesi aracından çıktımızla çok benzer (temel olarak direnç katmanınızla yapmaya çalıştığınız şeydir, ancak yokuş yukarı ve yokuş aşağı hareket arasında ayrım yapar). Yükseklik aralığınız ve direnç ağırlıklarınız göz önüne alındığında beklenen bu olabilir. Ancak, DEM'inize ve çıktınıza hızlı bir şekilde bakıldığında, sonuçlarınızın istediğiniz gibi görünmemesine neden olabilecek ve emin olmak için ikinci bir bakış atmak isteyebileceğiniz birkaç şey var gibi görünüyor.

1) Bölgenizin güneybatı kısmında nodata olarak kodlanmış gibi görünen büyük bir yığın var (DEM veya direnç katmanında). Bu funistonda, GIS nodata piksellerini esasen sonsuz bir dirence sahip olarak ele alır. (Bu yüzden bu ada şeyinin çok yüksek bir mesafe değeri vardır)

2) Yol mesafesi kullanıyorsanız ve dikey bir raster belirtiyor ancak dikey faktörler belirtmiyorsanız (veya tersi) veya bu iki parçadan herhangi biri yanlış belirtilmiş veya biçimlendirilmişse, işlev aletin bu bölümünü yürütemez ve geri kalanı çıktı üretmek için, ancak analizlerin dikey veya yatay yön kısımlarının düzgün yürütülmediğine dair herhangi bir uyarı veya gösterge vermeyecektir. Ayrıca, bazen program, bazı durumlarda ASCII dikey veya yatay faktör dosyasını kullanır ancak biçimlendirmeden bağımsız olarak diğerlerini (GUI kullanıyorsa çalışır, ancak python kullanmaz gibi) kullanmaz. Bu, bu aracın sorun gidermesini zorlaştırabilir. Genellikle gireriz ve farklı olup olmadıklarını görmek için bir koşudaki mesafe değerlerini dikey faktörlü ve dikey faktörlerle karşılaştırırız.

3) Test noktalarınızda birer birer çalıştırırsanız, aracın ne yaptığı hakkında daha fazla ayrıntı görebilirsiniz (şu anda işlev her bir pikselde yalnızca iki mesafeden daha kısa olanını görüntüleyebilirsiniz, çünkü işlev yalnızca her pikselden girişteki iki noktadan birine kadar olan mesafeler)

4) Bir çalışma alanı boyunca yüksek irtifa farklılıkları ve / veya VRMA faktörleri için geniş bir ağırlık aralığı olmadan, bir tepede yukarı ve aşağı hareket etmenin maliyetini de içeren bir analizden elde edilen çıktı genellikle pek farklı görünmemektedir. mesafenin öklid analizi. Bununla birlikte, aldığınız sayılar biraz farklı olacaktır ve bazı durumlarda en düşük maliyet yollarını eşlerseniz biraz farklı yollar alırlar.

5) Teknik olarak, dikey raster için giriş olarak DEM yerine z skorlu bir raster kullanmanız gerektiğini düşünüyorum, ancak her ikisi de forumlarda sıkça kullanılıyor ve en azından verilerimiz için çıktıdaki farklar minimum.

ESRI'nin bu konudaki belgeleri biraz dağınık, ancak dikey faktörlerin bu açıklaması oldukça iyi: http://webhelp.esri.com/arcgisdesktop/9.3/index.cfm?TopicName=Path%20Distance:%20adding%20more%20cost % 20complexity

Sitemizi kullandığınızda şunları okuyup anladığınızı kabul etmiş olursunuz: Çerez Politikası ve Gizlilik Politikası.
Licensed under cc by-sa 3.0 with attribution required.