أكثر

توليد شكل الشبكة في R؟

توليد شكل الشبكة في R؟


لدي ملف الأشكال التالي لشبكة خطوط الطول / العرض 10x10 درجة التي أنشأتها في QGIS. يمكنني قراءته إلى R باستخدام الحزمة rgdal.

الشبكة <-readOGR ("."، "GridShapeFile")

لديها السمات والهيكل التالي.

ملخص (شبكة) كائن من فئة SpatialPolygonsDataFrame إحداثيات: min max x -180 190 y -100 90 مُسقط: FALSE proj4string: [+ proj = longlat + datum = WGS84 + no_defs + ellps = WGS84 + towgs84 = 0،0،0] سمات البيانات: ID XMIN XMAX YMIN YMAX Min. : 0.0 دقيقة. : -180 دقيقة. : -170 دقيقة. : -100 دقيقة. : -90 1st Qu.:175.5 1st Qu .: -90 1st Qu: -80 1st Qu .: -60 1st Qu.:-50 Median: 351.0 Median: 0 Median: 10 Median: -10 Median: 0 المتوسط: 351.0 المتوسط: 0 المتوسط: 10 المتوسط: - 10 المتوسط: 0 الفصل الثالث: 526.5 الفصل الثالث: 90 الفصل الثالث: 100 الفصل الثالث: 40 الدور الثالث: 50 الحد الأقصى. : 702.0 كحد أقصى. : 180 كحد أقصى. : 190 كحد أقصى. : 80 كحد أقصى. : 90 # مثال صف من البيانات الشبكة [50 ،] الفئة: ميزات SpatialPolygonsDataFrame: 1 مدى: -60 ، -50 ، 70 ، 80 (xmin ، xmax ، ymin ، ymax) تنسيق. المرجع. : + proj = longlat + datum = WGS84 + no_defs + ellps = WGS84 + towgs84 = 0،0،0 متغيرات: 5 أسماء: ID، XMIN، XMAX، YMIN، YMAX min القيم: 49، -60، -50، 70، 80 قيمة قصوى: 49 ، -60 ، -50 ، 70 ، 80

أود أن أتعلم كيفية إنشاء نفس الشبكة / ملف الأشكال بالكامل داخل R ، بدلاً من QGIS ، بالإضافة إلى شبكات بأحجام أخرى (1x1،5x5 ، إلخ). الهدف النهائي هو استخدام over () في حزمة sp لتراكب هذه الشبكة على ملف آخر من المضلعات وحساب عدد الشبكات التي تتقاطع مع كل مضلع. أعرف بالفعل كيفية القيام بهذه الخطوة التالية (على ما أعتقد) ، أريد فقط معرفة كيفية إنشاء طبقة الشبكة في R.


ألق نظرة على وظيفة البيانات النقطية في حزمة البيانات النقطية. سيتيح لك إنشاء نقطية بمدى محدد وعدد الصفوف / الأعمدة ودقة الوضوح.

سأستخدم هنا خصائص ملخص البيانات الخاصة بك لإنشاء خطوط نقطية 100 × 100 ضمن النطاق المحدد. أقوم بتمرير كائن مدى لتحديد حدود x و y. يمكنك أيضًا استخدام الوسيطات المحددة (xmn ، xmx ، ymn ، ymx) داخل دالة البيانات النقطية.

مكتبة (نقطية) مكتبة (sp) r <- خطوط نقطية (المدى (مصفوفة (c (-180 ، -100 ، 190 ، 90) ، nrow = 2)) ، nrow = 100 ، ncol = 100 ، crs = "+ proj = longlat + datum = WGS84 + no_defs + ellps = WGS84 + towgs84 = 0،0،0 ") r [] <- 1: ncell (r) Summary (r) print (r) plot (r)

من السهل إجبار الكائنات النقطية على كائن sp شبكي باستخدام ؛

sp.r <- as (r، "SpatialPixelsDataFrame") class (sp.r) spplot (sp.r، "layer")

شاهد الفيديو: تعلم بايثون - دورة كاملة للمبتدئين