أكثر

إنشاء نوافذ متعددة الأضلاع للتحليل المكاني في r

إنشاء نوافذ متعددة الأضلاع للتحليل المكاني في r


أحاول تطبيع نمط النقطة المكانية باستخدام كثافة سكانية نقطية في r. مشكلتي هي أن الوظيفة التي أحاول تشغيلها هي التعرف على "المسافة البيضاء" حول البيانات النقطية التي أستخدمها كـ "NA" ، مما يعطيني رسالة الخطأ: قيم المتغير المشترك "Y" كانت NA أو غير محددة عند 32٪ (19493 من 60241) من النقاط التربيعية. حدث أثناء التنفيذ: ppm.ppp (X، ~ offset (log (Y))، المتغيرات المشتركة = list (Y)).

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

مكتبة (spatstat) مكتبة (maptools) مكتبة (نقطية) مكتبة (geostatsp) (rgdal) X <- readShapePoints ("~ / Documents / Data / Cyber ​​/ US / CONTIG_US / CONUS02_FREQ") Y <- نقطية ("~ / المستندات /Data/Cyber/US/US_POP1.tif ") ## تعيين المدى e <- مدى (-125، -67،25، 50) ## إنشاء نافذة لـ spatstat WIN <- مثل (e ،" SpatialPoints ") X.ppp <- as (X، "ppp") ## تحويل البيانات النقطية إلى كائن im لـ spatstat Y <- asImRaster (Y) Y <- Eval.im ((Y)) Y <- as.im (Y) ## Unmark X للنموذج: X <- unmark (X.ppp) Z <- ppm.ppp (X، ~ offset (log (Y))، المتغيرات المشتركة = list (Y))

تبدو الحبكة النقطية كما يلي:


توجد قيم NA في مخطط التربيع الخاص بك بسبب الطريقة التي تحدد بها نافذة الملاحظة (اوين). ال اوين لأحداثك (X) تم تعريفه على أنه النطاق الكامل للخطوط النقطية. تذكر أن أي خطوط نقطية ، بشكل افتراضي ، هي بنية أو مصفوفة من NxN الصفوف والأعمدة. في المثال الخاص بك ، الخلايا ذات الكثافة السكانية لها قيمة تتراوح من 0 إلى 41 ، ولكن هناك أيضًا خلايا ذات قيم NA ، (المساحات الفارغة حيث توجد المحيطات أو كندا أو المكسيك).

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

US <- readShapePoly ("USBORDER.shp") USw <- مثل (US، "owin") Xpoints <- readShapePoints ("~ / Documents / Data / Cyber ​​/ US / CONTIG_US / CONUS02_FREQ") X <- مثل (Xpoints، "ppp") [USw]

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

بالإضافة إلى ذلك ، حتى بعد استخدام مضلع لتعريف owin ، قد تحصل على قيم NA ، وقد يكون هذا بسبب دقة البيانات النقطية و / أو مقياس أصل المضلع. إذا حدث هذا ، يمكنك:

1) قم بتشغيل النموذج

Z <- جزء في المليون (** X ** ، ~ PD ، المتغيرات المشتركة = قائمة (PD = Y))

سيعطيك نفس التحذير ، "قيم المتغير" PD "كانت NA أو غير محددة عند n٪"

2) استخراج مخطط التربيع من Z

** Qz ** <- quad.ppm (Z، drop = TRUE) # ستسقط هذه الخطوة نقاطًا داخل مخطط التربيع الذي يحتوي على قيم NA

3) قم بتشغيل النموذج مرة أخرى ، ولكن هذه المرة ، استخدم مخطط التربيع المصحح

Z1 <- جزء في المليون (** Qz ** ، ~ PD ، المتغيرات المشتركة = القائمة (PD = Y))


شاهد الفيديو: طرق الاستكمال المكاني