v.krigegrass - عبر الإنترنت في السحابة

برنامج:

اسم

ضد كريج - ينفذ عمليات عادية أو جماعية للخرائط المتجهة.

الكلمات الرئيسية

ناقلات ، الاستيفاء ، النقطية ، kriging

موجز

ضد كريج
ضد كريج --مساعدة
ضد كريج إدخال=الاسم عمود=الاسم [الناتج=الاسم] [صفقة=سلسلة]
[نموذج=سلسلة[,سلسلة، ...]] [منع=عدد صحيح] [نطاق=عدد صحيح] [كتلة صلبة=عدد صحيح]
[عتبة=عدد صحيح] [input_var=الاسم] [-اعادة الكتابة] [-مساعدة] [-مطنب] [-هدوء]
[-ui]

الأعلام:
--الكتابة فوق
السماح لملفات الإخراج بالكتابة فوق الملفات الموجودة

--مساعدة
طباعة ملخص الاستخدام

- الإسراف
إخراج وحدة مطول

--هادئ
إخراج وحدة هادئة

--ui
فرض إطلاق مربع حوار واجهة المستخدم الرسومية

المعلمات:
إدخال=الاسم [مطلوب]
اسم خريطة متجه الإدخال
اسم خريطة متجهية تحتوي على بيانات عينة

عمود=الاسم [مطلوب]
اسم عمود السمة الذي يحتوي على قيمة عددية ليتم تحريفها

الناتج=الاسم
اسم لخريطة نقطية الإخراج
إذا حذفت، سيكون _kriging

صفقة=سلسلة
حزمة R للاستخدام
خيارات: gstat
الافتراضي: gstat

نموذج=سلسلة [، سلسلة، ...]
نموذج (نماذج) فاريوجرام
اتركه فارغًا لاختبار جميع الطرز (يتطلب الأمر تلقائيًا)
خيارات: نوج، الخبرة، Sph ، غاو، باستثناء، حصيرة، ستي، سير، لين، بس، قلم جاف، لكل، هول ، سجل، أسير,
سبل، رجل، خطأ ، الباحث

منع=عدد صحيح
حجم الكتلة (كتلة مربعة)
مقاس الكتله. تستخدم من قبل بلوك كريجينج.

نطاق=عدد صحيح
قيمة النطاق
يتم إصلاحه تلقائيًا إذا لم يتم ضبطه

كتلة صلبة=عدد صحيح
قيمة الكتلة
يتم إصلاحه تلقائيًا إذا لم يتم ضبطه

عتبة=عدد صحيح
قيمة عتبة
يتم إصلاحه تلقائيًا إذا لم يتم ضبطه

input_var=الاسم
اسم الخريطة النقطية لتباين الإخراج
إذا حذفت، سيكون _kriging.var

الوصف

ضد كريج يسمح بإجراء عمليات Kriging في بيئة GRASS GIS ، باستخدام برنامج R.
وظائف في الخلفية.

الملاحظات

ضد كريج هي مجرد واجهة أمامية لـ R. الخيارات والمعلمات هي نفسها التي يقدمها
بكيج com.automap و gstat.

تعتمد Kriging ، مثل طرق الاستيفاء الأخرى ، بشكل كامل على ميزات بيانات الإدخال.
يتم تشجيع التحليل الاستكشافي للبيانات لمعرفة القيم المتطرفة، والاتجاهات، وتباين الخواص،
توزيعات غير متساوية وبالتالي اختر خوارزمية kriging التي ستعطي أكثر
نتيجة مقبولة. تعد المعرفة الجيدة بمجموعة البيانات أكثر قيمة من مئات منها
المعلمات أو الأجهزة القوية. انظر كتاب إسحاق وسريفاستافا ، الشامل والواضح
حتى لو عفا عليها الزمن بعض الشيء.

تبعيات
R نظام البرمجيات >= 2.x

rpy2
بايثون ملزم بـ R. Note! rpy الإصدار 1 غير مدعوم.

R بكيج أوتوماب غستات, rgrass7 و رجيو.
يعد automap اختياريًا (يوفر ملاءمة فوتوغرافي تلقائية). قم بتثبيت الحزم عبر R.
سطر الأوامر (أو واجهة المستخدم الرسومية المفضلة لديك):
install.packages ("rgeos" ، dep = T)
install.packages("gstat"، dep=T)
install.packages ("rgrass7"، dep = T)
install.packages("automap"، dep=T)

ملاحظة لـ ديبيان جنو / لينكس
تثبيت التبعيات. انتباه! بيثون-rpy IS لا مناسب.:
aptitude تثبيت R python-rpy2
لتثبيت حزم R ، استخدم وظائف R المذكورة أعلاه (كجذر أو كمستخدم) ، إما
حزم دبيان [5] ، أضفها إلى قائمة المستودعات 32 بت أو 64 بت (التقط ملف
الخط المناسب):
ديب http://debian.cran.r-project.org/cran2deb/debian-i386 اختبارات/
ديب http://debian.cran.r-project.org/cran2deb/debian-amd64 اختبارات/
واحصل على الحزم عبر aptitude:
aptitude install r-cran-gstat r-cran-rgrass7

ملاحظة لـ ويندوز
قم بتجميع GRASS GIS باتباع هذا الدليل. يمكنك أيضًا استخدام Linux في جهاز افتراضي. أو
قم بتثبيت Linux في قسم منفصل من HD. هذا ليس مؤلمًا كما يبدو ،
هناك الكثير من الأدلة عبر الإنترنت لمساعدتك.

حساب الوقت مسائل
يرجى ملاحظة أنه على الرغم من وجود عدد كبير من نقاط بيانات الإدخال و/أو دقة المنطقة العالية
المساهمة في إنتاج أفضل، وكلاهما سيؤدي أيضًا إلى إبطاء حساب kriging.

أمثلة

مثال Kriging على أساس خريطة الارتفاع (مجموعة بيانات Spearfish).

جزء 1: عشوائية أخذ العينات 2000 نقطة متجه من خريطة الارتفاع المعروفة. كل نقطة سوف
احصل على قيمة الارتفاع من بيانات الارتفاع النقطية، كما لو أنها جاءت من مسح نقطة.
# تقليل الدقة لهذا المثال
ز - خطوط المسح للمنطقة = الارتفاع - درجة الدقة = 100
v. الإخراج العشوائي = Rand2k_elev npoints = 2000
v.db.addtable Map=rand2k_elev columns="الدقة المزدوجة للارتفاع"
v.what.rast map = rand2k_elev raster = عمود الارتفاع = الارتفاع
جزء 2: إزالة نقطة تفتقر ارتفاع سمات. النقاط التي تم أخذ عينات منها عند حدود
لم تحصل خريطة الارتفاع على أي قيمة. v.krige ليس لديه أي إجراء مفضل للتعامل مع لا
قيم البيانات، لذلك يجب على المستخدم التحقق منها وتحديد ما يجب فعله (إزالة النقاط، وملء
بقيمة أقرب نقطة ، املأ بالمتوسط ​​العالمي / المحلي ...). في التالي
سطر من التعليمات البرمجية ، تتم إزالة النقاط التي لا تحتوي على بيانات من الخريطة.
v.extract المدخلات = إخراج rand2k_elev = rand2k_elev_filt حيث = "الارتفاع ليس فارغًا"
تحقق من نتيجة السطر السابق (يجب أن يكون "عدد السمات الفارغة" 0):
v.univar map = نوع rand2k_elev_filt = عمود النقطة = الارتفاع
جزء 3: إعادة DEM من خلال كريجينج. أبسط طريقة للتشغيل ضد كريج من CLI يستخدم
ملاءمة variogram تلقائية (ملاحظة: تتطلب حزمة R التلقائية). اسم خريطة الإخراج اختياري،
تقوم الوحدات بإنشائها تلقائيًا بإلحاق "_kriging" باسم خريطة الإدخال وأيضًا
الشيكات للكتابة. إذا تم تحديد output_var ، فسيتم أيضًا إنشاء خريطة التباين.
يتم توفير الملاءمة التلقائية للفاريوجرام بواسطة R package automap. تم اختبار نماذج varogram بواسطة
وظائف التركيب هي: أسي ، كروي ، غاوسي ، ماتيرن ، إم ستاين
المعلمات. تتوفر مجموعة واسعة من النماذج من حزمة gstat ويمكن اختبارها
على واجهة المستخدم الرسومية عبر الرسم المتغير. إذا تم تحديد نموذج في CLI ، فسيكون أيضًا عتبة ،
يجب توفير قيم الكتلة والنطاق ، وإلا فسيظهر خطأ (انظر الثانية
مثال على ضد كريج أمر).
# تناسب variogram التلقائي
v.krige input=rand2k_elev_filt column=elevation
الإخراج = Rand2k_elev_krigingoutput_var=rand2k_elev_kriging_var
# تحديد نموذج variogram، وإنشاء خريطة التباين أيضًا
v.krige input=rand2k_elev_filt column=elevation
الإخراج = Rand2k_elev_filt_krigingput_var=rand2k_elev_filt_kriging_var
الموديل = Mat sill = 2500 كتلة صلبة = 0 نطاق = 1000
أو قم بتشغيل wxGUI، لملاءمة المتغير بشكل تفاعلي واستكشاف الخيارات:
ضد كريج
احسب الشحن تنبؤ خطأ:
r.mapcalc "rand2k_elev_kriging_pe = sqrt (rand2k_elev_kriging_var)"
r.univar الخريطة = الارتفاع
r.univar Map=rand2k_elev_kriging
r.univar Map=rand2k_elev_kriging_pe
تظهر النتائج أخطاء عالية ، مثل تقنيات kriging (العادية و block kriging)
غير قادر على التعامل مع مجموعة بيانات ذات اتجاه، مثل ذلك المستخدم في هذا المثال: الارتفاع هو
أعلى في الزاوية الجنوبية الغربية وأسفل في الزاوية الشمالية الشرقية. يمكن أن يعطي kriging العالمي
نتائج أفضل بكثير في هذه الحالات لأنه يمكن التعامل مع هذا الاتجاه. وهي متوفرة في حزمة R
gstat وسيكون جزءًا من إصدار v.krige المستقبلي.

استخدم v.krigegrass عبر الإنترنت باستخدام خدمات onworks.net