open source - إحصائيات المناطق قيس




open-source statistics (3)

هل هناك أي بديل مفتوح المصدر لأداة الإحصاءات الإقليمية (أرتجيس سبيشيال أناليست)؟ ما هي أفضل أداة (التي يمكنني استخدامها في البرنامج النصي) دور عد إحصاءات الملفات النقطية؟


تصحيح لي إذا أنا مخطئ، روبرث، ولكن أعتقد زونال () يتطلب أن المناطق هي بالفعل "النقطية" في بعض المعنى، في حين عدة مرات واحد سوف تريد إحصاءات الخلايا النقطية التي تقع ضمن المضلعات. طرق تراكب مختلفة في R داخل حزمة سب (انظر:؟ "تراكب-طرق") ضرورية لهذا، على الرغم من أنني إذا كنت مخطئا سأكون سعيدا لسماع ذلك. أنا أفضل تماما حزمة النقطية أكثر من استخدام سباتيالغريدزداتافرامس، ولكن أعتقد أن المرء يجب أن تعتمد على الطبقات سب لخلط المضلعات والبيانات الشبكية. ما هو موافق، إلا يصبح مشكلة لأنه يفتقر إلى إدارة الذاكرة كبيرة من حزمة النقطية، الأمر الذي يجعل نقطة في المضلعات على غرار عمليات من الصعب حقا القيام به في R على راسترز كبيرة.

كما أنني أؤمن، ولكن لم أحاول، أن هذا يمكن أن يتم داخل غراس، و / أو من خلال قيس، مع الإصدار التالي من قيس (1.7) أن يكون نوعا من المدمج في ميزة المناطق.


حزمة راسترستاتس هي أداة مفتوحة المصدر لطيفة التي عملت بشكل جيد بالنسبة لي: http://blog.perrygeo.net/2013/09/24/python-raster-stats/

بدأت استخدامه كعمل حول لأن طريقة أركبي زونالستاتيستيكش كان ينتج النقطية إشكالية التي تؤدي إلى خطأ غريب عند محاولة تحويل النقطية إلى مصفوفة ( https://gis.stackexchange.com/questions/110274/save- فشل على النقطية، وجوه خلق من-نومبياريتراوستر ). عملت راسترستاتس بشكل جيد وقدمت حل فعال لمشكلتي.


يمكنك القيام بذلك مع غراس باستخدام أساليب مختلفة. أي واحد هو الأكثر ملاءمة سوف تعتمد على البيانات الخاصة بك والناتج المطلوب. لاحظ أنه يمكنك استخدام غراس أيضا من داخل قيس باستخدام مربع الأدوات غراس أو سيكستانت.

لنفترض أن لديك:

  • خريطة متجهة، مثل vector_zones مع المناطق المحددة في myzones العمود في جدول الخصائص.
  • "قيم" طبقة النقطية التي تريد حساب إحصاءات منطقتك

r.statistics

لاستخدام r.statistics، تحتاج أولا إلى تحويل خريطة المتجه إلى طبقة النقطية، والتي يمكنك القيام به مع v.to.rast . بعد ذلك، استخدم r.statistics لحساب إحصاءات المناطق.

v.to.rast input=vector_zones output=zones column=myzones
r.statistics base=zones cover=values out=outputmap method=average

هذا سيعطيك طبقة جديدة مع إحصائية المنطقة المحددة، والتي يمكن أن تكون متوسطة، واسطة، وسيط، والتباين، وما إلى ذلك (انظر الرابط صفحة الرجل أعلاه).

r.univar

وظيفة r.univar أيضا يعمل على طبقات النقطية.

v.to.rast input=vector_zones output=zones column=myzones    
r.univar map=values zones=zones output=output.file fs=;

الإخراج هو جدول مع الإحصائيات المحلية.

v.rast.stats

هذا لا يتطلب منك تحويل طبقة المتجه إلى طبقة النقطية (يتم ذلك داخليا). تقوم الدالة بحساب الإحصائيات الأساسية المتغيرة أحادية المتغير لكل فئة متجه (كات) من الخريطة النقطية.

v.rast.stats vector=vector_zones layer=1 raster=values column_prefix=val

يتم تحميل النتائج إلى جدول سمة خريطة المتجهات.