أكثر

مقطع النقطية بواسطة shapefile في أجزاء؟

مقطع النقطية بواسطة shapefile في أجزاء؟


أحاول كتابة أمر نصي .sh لمقطع ملف الأشكال النقطية بواسطة مضلع. يحتوي ملف الشكل على حقل اسم بالاسم الفريد لعنصر المضلع. على الإخراج يجب أن تكون البيانات النقطية مع اسم المضلعات في شكل. أنا أستخدم Linux ، GDAL. الآن أنا استخدم الأمر

gdalwarp -dstnodata -9999 -cutline data / shapefile.shp data / input.tiff data / output.tiff

لكن في الإخراج أتلقى خط نقطي واحد.

مثال على الصورة:


أعد تنظيم ملف الشكل الخاص بك بحيث يحتوي ملف شكل واحد على ميزة واحدة (A ، B ، C في حالتك) فقط

ثم استخدم حلقة مثل

لأني في أ ب ج ؛ هل gdalwarp -cutline $ i.shp… $ i.tif تم

لإنشاء كل الناتج النقطي.

مثال على البرنامج النصي:

#! / bin / sh # "shp" - مجلد لملفات الأشكال # "outputraster" - مجلد لإخراج النقطيات cd / home / user / cliprasters / # حذف الملفات السابقة rm ./outputraster/*.tiff # تشغيل لـ fname1 في ./ shp / *. shp do fname2 = "$ (basename $ fname1 .shp)" gdalwarp -dstnodata -9999 -cutline ./shp/$fname2.shp -crop_to_cutline -of GTiff ./input.tiff ./outputraster/$fname2. شجار القيام به

حل أبسط ، على الرغم من أنه أقل أناقة هو الاستخدام ببساطة راستريو, جدال و حبيبي. لنفترض أنه في ملف الأشكال لدينا عمود يسمى "عمود" يحتوي على أسماء المعالم (مضلعاتنا) ، فلنفترض أن المضلعات لدينا هي "أ" و "ب" و "ج".

استيراد gdal، os، rasterio import numpy كـ np im = 'image.tif' shp = 'shapefile.shp' # افتح الصورة باستخدام rasterio وأعد تشكيلها إلى np.array a = rasterio.open (im) .read () a = a.reshape ((a.shape [1]، a.shape [2])) # gdal غريب عند العمل مع المجلدات في بعض الأحيان ، # لذلك ننتقل فقط إلى المجلد حيث الصورة os.chdir (مجلد) #folder with image #now نحن نشغل gdalwarp بالمعلمات المناسبة في حلقة # فوق مضلعاتنا. هذا هو الجزء الذي يمكن تحسينه # باستخدام fiona أو geopandas ، شيء لاستخراج الميزات # أسماء وإفراغها في قائمة. # الحيلة مع أسماء الميزات هي استخدام علامات الاقتباس المفردة والمزدوجة ، مثل "'this'" ، لأن gdal يريد الحصول على شيء مثل "this" كمعامل. polnames = ["'A'"، "'B'"، "'C'"] os.chdir (مجلد) لـ p في الأسماء المتعددة: os.system ('gdalwarp -cutline٪ s -cwhere "٪ s =٪ s "-crop_to_cutline٪ s٪ s_٪ s '٪ (shp، c، p، im، p، im))

الآن لديك الصورة مقطوعة إلى أجزاء وفقًا لميزاتك في ملف الأشكال. ثم نفتح الصور ونحصل على البيانات.

# افتح صورة واحدة باستخدام rasterio b = rasterio.open ('Aim.tif'). اقرأ (). reshape ((b.shape [1]، b.shape [2])) # إذا أردت ، يمكنك تسطيح و لديك صفيف 1d b = b.flatten ()

فعله! قلت إنه ليس أنيقًا ، هناك مجال كبير للتحسين ، لكنه يعمل.


نظم المعلومات الجغرافية

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

أرسل بريدًا إلكترونيًا إلى أمين المكتبة لمناقشة جدولة وقت خارج الأوقات المدرجة في المجدول حيث قد تتوفر أوقات أخرى.


كيفية استخراج المعلومات من البيانات النقطية إلى ملف الشكل

أنا طالب دكتوراه في الاقتصاد ومبتدئ في ArcGIS ، لذا سامحني لعدم كفاءتي.

لدي المشكلة التالية. لدي مجموعة بيانات نقطية تحتوي على معلومات حول عدد البرق على شبكة خلوية بحجم 50 كم × 50 كم فوق أوروبا.
تأتي المعلومات في الواقع من ملف ASCII ولكن ليس لدي مشكلة في التحويل إلى مجموعة بيانات نقطية.
لدي أيضًا ملف أشكال بحدود جميع مقاطعات أوروبا ، ويمكنني تركيب مجموعتي البيانات بشكل صحيح وتحديد إسقاط مشترك لـ 2.

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

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

لست متأكدًا مما يجب أن أفعله وما هي الأدوات أو أنواع الملفات التي يجب أن أعمل بها.
حاولت قص البيانات النقطية بملف الشكل ولكن لم يتمكن من استخدام البيانات النقطية كمدخل.
حاولت أيضًا حساب إحصائيات المنطقة كجدول (ضمن مربع أدوات التحليل المكاني) ، لكن النتيجة هي جدول من صف واحد فقط.
أخيرًا ، حاولت حساب إحصائيات المنطقة ولكن تعطل ArcGIS. أقوم بتشغيل ArcGIS بترخيص ArcInfo للطلاب ، راجع للشغل.

آمل أن أوضح نفسي: أخبرني إذا كان بإمكاني توضيح شيء ما.


القيود الهندسية

  • يوجد حد للحجم يبلغ 2 جيجا بايت لأي ملف مكون من ملفات الشكل ، والذي يترجم بحد أقصى 70 مليون نقطة كحد أقصى. يعتمد العدد الفعلي لميزات الخط أو المضلع التي يمكنك تخزينها في ملف شكل على عدد الرؤوس في كل سطر أو مضلع (الرأس يعادل نقطة).
  • لا تحتوي ملفات الأشكال على تفاوت x و y مثل فئات معالم قاعدة البيانات الجغرافية. تفاوت x و y هو الحد الأدنى للمسافة بين الإحداثيات قبل اعتبارهما متساويين. يتم استخدام تفاوت x و y عند تقييم العلاقات بين المعالم ضمن نفس فئة الميزة أو بين العديد من فئات الميزات المختلفة. يتم استخدامه أيضًا على نطاق واسع عند تحرير الميزات. إذا كنت تقوم بتنفيذ أي نوع من العمليات التي تتضمن مقارنة بين الميزات ، مثل استخدام أدوات التراكب أو أداة Clip أو أداة تحديد الطبقة حسب الموقع أو أي أداة تأخذ فئتين أو أكثر من فئات المعالم كمدخلات ، فيجب عليك استخدام ميزة قاعدة البيانات الجغرافية الفئات (التي لها تفاوت x ، y) بدلاً من ملفات الأشكال.
  • قد يستغرق ملف الشكل من ثلاثة إلى خمسة أضعاف مساحة قاعدة البيانات الجغرافية للملف أو SDE بسبب أساليب ضغط الشكل.
  • تدعم ملفات الأشكال التطابقات المتعددة ولكنها تفتقر إلى الدعم لإمكانيات الضربات المتعددة المتقدمة التالية:
    • إحداثيات الملمس
    • القوام ولون الجزء
    • معايير الإضاءة

    سؤال حول وقت تشغيل مقطع نقطي في ArcMap 10.3

    & # x27ve تم تشغيل مقطع نقطي على ملف 1 جيجا بايت بارتفاع الشجرة العالمي. أقوم & # x27m بقصها إلى ملف أشكال من مناطق اليابسة (بدون حدود داخلية) لأن البيانات النقطية الأصلية تشمل البحر أيضًا وهذا يجعل تطبيق نظام ألوان جيد أمرًا صعبًا.

    & # x27m تشغيل ArcMAP 10.3 و & # x27s على جهاز كمبيوتر جديد تمامًا مع تثبيتات جديدة لكل شيء. تم تشغيله ليلا ونهارا الآن لأكثر من 4 أيام بقليل. يُظهر Resource Monitor & amp Task Manager أن Runtimelocalserver.exe يبتعد باستخدام ما بين 12 و 18 ٪ من وحدة المعالجة المركزية ومقادير الذاكرة تتغير باستمرار. تعمل & # x27s كعملية في الخلفية ولا يزال بإمكاني معالجة الخريطة وحفظها وإضافة البيانات وما إلى ذلك.

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

    كم من الوقت أترك هذا يعمل؟ هل هو في الواقع عمل أي شيء ، أم أنه مجرد مظهر لفعل شيء ما؟


    حول QGIS

    كما هو موضح على موقع الويب الخاص بهم ، QGIS هو "نظام معلومات جغرافية مجاني ومفتوح المصدر" ، وأحدث إصدار (v2.6) من البرنامج متاح للتنزيل على موقع الويب الرسمي الخاص بهم. مثل منصات نظم المعلومات الجغرافية الأخرى ، يمكن استخدامه لتحرير مجموعة متنوعة من البيانات الجغرافية المكانية وإنشائها وتمثيلها بصريًا وتحليلها وتصديرها وهي مُرخصة بموجب رخصة جنو العمومية العامة. QGIS هو مشروع رسمي لمؤسسة Open Source Geospatial Foundation (OSGeo). يعمل على Linux و Unix و Mac OSX و Windows و Android ويدعم العديد من تنسيقات ووظائف المتجهات والنقطية وقواعد البيانات.


    تقديم نظم المعلومات الجغرافية باستخدام ArcGIS: منهج مصنف لتعلم نظم المعلومات الجغرافية ، الإصدار الثالث

    تقديم نظم المعلومات الجغرافية باستخدام ArcGISيدمج مقدمة واسعة إلى GIS مع كتاب مصنف خاص بالبرنامج لـ ArcGIS الخاص بـ Esri. حيث تقوم معظم الدورات باستخدام نصين منفصلين ، أحدهما يغطي نظم المعلومات الجغرافية والآخر البرنامج ، يمكّن هذا الكتاب الطلاب والمعلمين من استخدام نص واحد بنهج متكامل يغطي كلاهما في مجلد واحد مع مفردات مشتركة وأسلوب تعليمي.

    تركز هذه النسخة المنقحة على آخر تحديثات البرامج و mdashArcGIS 10.0 و 10.1. بالإضافة إلى تغطيته الناجحة بالفعل ، يسمح الكتاب للطلاب بتجربة نشر الخرائط على الإنترنت من خلال تمارين جديدة ، ويقدم فكرة البرمجة باللغة التي اختارتها Esri للتطبيقات (مثل Python). يتم حزم قرص DVD مع الكتاب ، كما في الطبعات السابقة ، يحتوي على بيانات للعمل على جميع التمارين.

    هذا الكتاب الدراسي الكامل سهل الاستخدام:

    • يتم تحديثه لأحدث إصدارات ArcGIS و mdashArcGIS 10.0 و 10.1
    • يقدم المفاهيم الأساسية لنظم المعلومات الجغرافية والموضوعات اللازمة لفهم تحليل المعلومات المكانية
    • يوفر قدرة كبيرة على تشغيل الأدوات المهمة في ArcGIS
    • يوضح الإمكانات الجديدة لـ ArcGIS 10.0 و 10.1
    • يوفر أساسًا للدراسة المتقدمة لنظم المعلومات الجغرافية ودراسة مجال GISc الناشئ حديثًا

    تقديم نظم المعلومات الجغرافية باستخدام ArcGIS، الطبعة الثالثة هو الدليل المثالي للطلاب الجامعيين الذين يأخذون دورات مثل مقدمة في نظم المعلومات الجغرافية ، وأساسيات نظم المعلومات الجغرافية ، ومقدمة إلى ArcGIS Desktop. وهو أيضًا دليل مهم للمحترفين الذين يتطلعون إلى تحديث مهاراتهم في ArcGIS 10.0 و 10.1.


    تعليمات

    1. افتح QGIS لإنشاء مشروع QQIS جديد يمكنك العمل فيه من أجل الإسناد الجغرافي لصور الخرائط الممسوحة ضوئيًا (راجع دليل تثبيت QGIS إذا لم يكن البرنامج مثبتًا لديك بالفعل). يرجى أن تضع في اعتبارك أن هذا الدليل قد تم تطويره لـ QGIS 3.x ، لذا إذا كان لديك QGIS 2.x مثبتًا على جهاز الكمبيوتر الخاص بك ، فستحتاج إلى ترقية التثبيت قبل المتابعة للتأكد من أن هذه الإرشادات تتوافق مع إصدار البرنامج الخاص بك. إذا كان QGIS مثبتًا لديك بالفعل ولكنك غير متأكد من الإصدار الذي لديك على جهاز الكمبيوتر الخاص بك ، فحاول فتح QGIS وانظر عن كثب إلى شاشة البداية التي تظهر أثناء تحميل البرنامج - يجب الإشارة إلى رقم الإصدار بشكل بارز. إذا كان لديك بالفعل QGIS مفتوحًا ، فيمكنك أيضًا التحقق من رقم الإصدار الخاص بك عن طريق الوصول إلى يساعد القائمة في الجزء العلوي من الشاشة واختيار تحقق من إصدار QGIS.
    2. بمجرد فتح QGIS بنجاح ، ستلاحظ أن مشروعك يبدو فارغًا - وهذا بسبب عدم وجود أي بيانات تم تحميلها عليه حاليًا. هذا جيد لأننا سنختار البيانات التي سيتم استخدامها في عملية الإسناد الجغرافي في الخطوة التالية. قبل بدء عملية الإسناد الجغرافي ، من الأفضل تحديد نظام إحداثيات مشروع QGIS الخاص بك لأن هذا سيجعل بعض الخطوات القادمة في هذه العملية أسهل قليلاً. للقيام بذلك ، انقر فوق مشروع في الجزء العلوي من نافذة QGIS (Windows) أو أعلى الشاشة (MacOS) وحدد ملكيات من القائمة التي تظهر. في نافذة Project Properties (خصائص المشروع) ، انقر فوق CRS في العمود الأيسر ثم حدد CRS الذي ترغب في أن يستخدمه مشروعك (والطبقات التي تم تحديدها جغرافيًا) في لوحة النافذة الرئيسية. WGS 84 (EPSG: 4326) هو النظام المرجعي القياسي للإحداثيات (CRS) المستخدم في خرائط الإسناد الجغرافي التي تعد جزءًا من مجموعة خرائط PCL ويوصى أيضًا بتحديد CRS هذا ما لم يكن لديك سبب محدد لاختيار خيار آخر. أسهل طريقة للعثور على ملف WGS 84 (EPSG: 4326) نظام الإحداثيات هو كتابة WGS 84 في شريط التصفية بالقرب من الجزء العلوي من نافذة محدد نظام الإحداثيات المرجعي ثم تصفحها في القائمة الضيقة الآن ضمن أنظمة الإحداثيات المرجعية في العالم (يجب أن تظهر بالقرب من أعلى هذه القائمة).
    3. لبدء الإسناد الجغرافي انقر فوق النقطية في الجزء العلوي من نافذة QGIS (نظام التشغيل Windows) أو أعلى الشاشة (نظام التشغيل MacOS). في القائمة المنسدلة التي تظهر ، انقر فوق Georeferencer و hellip لإحضار نافذة Georeference. إذا كنت لا ترى ملف Georeferencer و hellip الخيار في هذه القائمة ، يشير إلى أن ملف Georeferencer و hellip لم يتم تنشيط البرنامج المساعد بعد. لتنشيط البرنامج المساعد ، انقر فوق الإضافات في الجزء العلوي من نافذة QGIS الخاصة بك ثم حدد إدارة وتثبيت الإضافات. في القائمة المنسدلة لسحب نافذة مدير البرنامج المساعد التي تسرد جميع مكونات QGIS الإضافية المتاحة. قم بالتمرير لأسفل في قائمة المكونات الإضافية حتى تجدها Georeferencer GDAL. تأكد من تحديد المربع الموجود بجوار هذا المكون الإضافي في القائمة لتنشيط المكون الإضافي ثم إغلاق نافذة مدير البرنامج المساعد.

    4. الآن ، انقر فوق افتح النقطية في الزاوية اليسرى العلوية من نافذة Georeferencer واستعرض موقع نظام الملفات لملف صورة الخريطة الذي تريد تحديده جغرافيًا ، ثم انقر فوق فتح. إذا لم يكن لديك ملف صورة خريطة خاص بك للتدرب عليه ، أو ترغب فقط في جعل اتباع الإرشادات المتبقية في هذا الدليل سهلاً قدر الإمكان ، يمكنك استخدام نفس الخريطة الموضحة في لقطات الشاشة القادمة لعملية الإسناد الجغرافي. يمكنك تنزيل هذه الخريطة من UT Libraries & # 39 PCL Maps Collection على http://legacy.lib.utexas.edu/maps/ams/africa_1m/txu-pclmaps-oclc-6587819-na-32.jpg ، احفظها في جهاز الكمبيوتر الخاص بك ، ثم استعرض موقع ملف صورة الخريطة من نافذة Georeferencer كما هو موضح أعلاه. طالما أنك قد اخترت بالفعل بنجاح نظام مرجعي إحداثي لمشروعك ، يجب أن تكون جاهزًا للانتقال إلى الخطوة التالية. إذا كانت هناك أية مشكلات في تحديد نظام الإحداثي المرجعي ، فستلاحظ ظهور نافذة محدد نظام الإحداثي المرجعي بعد تحديد صورة الخريطة التي تريد الإشارة إليها جغرافيًا. إذا ظهرت لك هذه النافذة ، فيجب عليك تحديد ملف WGS 84 (EPSG: 4326) نظام الإحداثيات. ال WGS 84 (EPSG: 4326) يعد نظام الإحداثيات خيارًا جيدًا في معظم الحالات (وإذا كنت تقوم بإسناد الخرائط بشكل طوعي من مجموعة خرائط PCL ، فمن الأهمية بمكان أن يتم تحديد هذا النظام الإحداثي لضمان التوحيد القياسي). بمجرد تحديد النظام الإحداثي الخاص بك ، انقر فوق "موافق" وسترى خريطتك تظهر في نافذة Georeference.
    5. ستحتوي معظم الخرائط الموحدة مثل تلك التي تعد جزءًا من سلسلة AMS أو USGS أو JOG على قيم درجة نظام التنسيق الجغرافي (GCS) معروضة بشكل بارز في كل ركن من أركان الخريطة الأربعة على ورقة الخريطة. هذه المعلومات مفيدة للغاية عند الإسناد الجغرافي لأنها ستسمح لنا بتحديد مكان كل ركن من أركان الخريطة على الأرض بحيث يمكن رسمها بشكل صحيح فوق الطبقات الأخرى ذات المرجعية الجغرافية في برنامج نظم المعلومات الجغرافية. للبدء في إنشاء نقاط التحكم الأرضية (GCPs) التي سنستخدمها للإشارة إلى صورة الخريطة جغرافيًا ، قم بالتكبير إلى الزاوية اليسرى العلوية من خريطتك حيث نأمل أن ترى إحداثيات GCS. ثم انقر فوق إضافة زر GCP في شريط أدوات نافذة Georeferencer وانقر مرة أخرى على الزاوية عندما ينتهي نطاق الخريطة. هذا سوف يطرح ال أدخل إحداثيات الخريطة نافذة حيث يمكنك الآن إدخال قيم الإحداثي التي تراها معروضة في الزاوية التي يتم تكبيرها حاليًا. اقرأ إرشادات تنسيق الإحداثيات في هذه النافذة بعناية وتأكد من كتابة القيم الصحيحة في الحقلين X و Y لضمان إدخال الإحداثيات بشكل صحيح. كرر هذه العملية الآن حتى تقوم بإنشاء GCPs لجميع الزوايا الأربع للخريطة.
    6. تحقق مرة أخرى من إدخال إحداثيات نقاط GCP المدرجة في جدول GCP الخاص بك بشكل صحيح وإذا كان كل شيء يبدو جيدًا ، فابدأ وانقر على حفظ نقاط GCP باسم زر في شريط أدوات Georeferencer الخاص بك. سيظهر هذا نافذة منبثقة ستطالبك بتحديد موقع الحفظ واسم الملف لتصدير نقاط GCP. يجب أن تعطي هذا الملف نفس الاسم لملف صورة الخريطة الأصلي غير المرجعي (تأكد من إزالة .jpg أو أي امتداد صورة آخر في اسم الملف الافتراضي الذي تقترحه الأداة حتى لا يتم حفظه بامتداد مزدوج مثل file.jpg .نقاط). يجب عليك أيضًا حفظ هذا الملف في نفس الموقع الذي تخطط لحفظ صورة الخريطة المحددة جغرافيًا.
    7. بعد حفظ نقاط GCP الخاصة بك ، انقر فوق زر السهم الأخضر لملف ابدأ بالإشارة إلى البيانات الجغرافية الأداة الموجودة في شريط أدوات Georeferencer الخاص بك. في المرة الأولى التي تقوم فيها بالإشارة الجغرافية للخريطة ، سيظهر هذا نافذة إعدادات التحويل المنبثقة حيث يتعين عليك تحديد بعض المعلمات المهمة التي ستحدد كيفية تنفيذ عملية الإسناد الجغرافي. في الجزء العلوي من هذه النافذة ، اضبط نوع التحويل على متعدد الحدود 1، اضبط طريقة إعادة التشكيل على مكعب، والهدف SRS إلى WGS 84 (EPSG: 4326). ضمن إعدادات الإخراج ، انقر فوق الزر & hellip في صف Output Raster لتحديد موقع لحفظ البيانات النقطية التي تم الإشارة إليها جغرافيًا على جهاز الكمبيوتر الخاص بك. عند تحديد موقع الحفظ هذا ، تأكد من الاحتفاظ بالاسم الأصلي للملف ثم إلحاق المعلومات التالية باسم الملف _transformationtype_resamplingmethod_compressionmethod_coordinatesystem باستخدام الشرطة السفلية لفصل كل سمة مهمة. وبالتالي ، يمكن حفظ ملف صورة الخريطة المسمى txu-pclmaps-oclc-6587819-na-32 كـ txu-pclmaps-oclc-6587819-na-32_polynomial1_cubic_lzw_wgs84.tif. تأكد أيضًا من تحديد ملف LZW نوع الضغط. بمجرد تحديد المعلمات بشكل صحيح ، انقر فوق الزر "موافق" لتطبيق إعدادات التحويل وإغلاق نافذة إعدادات التحويل. في هذه المرحلة ، ستحتاج إلى النقر مرة أخرى على زر السهم الأخضر لبدء عملية الإسناد الجغرافي باستخدام إعدادات التحويل هذه. يجب أن تكتمل عملية الإسناد الجغرافي في غضون ثوانٍ قليلة وسترى بعد ذلك صورة الخريطة التي تم تحديدها جغرافيًا حديثًا تظهر في لوحة مشروع QGIS.
    8. للتحقق من أن خريطتك قد تم تحديدها جغرافيًا بشكل صحيح ، يمكنك إضافة خريطة أساس OpenStreetMap إلى مشروع QGIS الخاص بك والذي يمكنك استخدامه للتحقق من وضع صورة الخريطة عن طريق مقارنة مواقع المعالم البارزة في كلتا الخريطتين. للقيام بذلك ، انقر فوق XYZ البلاط في لوحة المتصفح ثم اسحب ملف OpenStreetMap طبقة في لوحة الطبقات الخاصة بك ، مع الحرص على وضعها أسفل طبقة صورة الخريطة التي تم تصديرها مؤخرًا.
    9. الآن وقد تم تحديد موقع خريطتك جغرافيًا ، قد ترغب أيضًا في قص صورة الخريطة الممسوحة ضوئيًا لإزالة الطوق (المساحة الهامشية البيضاء حول الخريطة التي تحتوي على معلومات مثل عنوان الخريطة و rsquos والمقياس والمفتاح وما إلى ذلك) لتسهيل استخدام صورة نقطية في برنامج GIS. للقيام بذلك انقر فوق النقطية في الجزء العلوي من نافذة QGIS (نظام التشغيل Windows) أو أعلى الشاشة (نظام التشغيل MacOS). في القائمة المنسدلة التي تظهر ، انقر فوق قص النقطية حسب المدى لإحضار نافذة الأداة.
    10. في نافذة Clip Raster by Extent المنبثقة ، انقر فوق الزر & hellip بجوار حقل مدى القطع (xmin ، xmax ، ymin ، ymax) وحدديختارمدى على قماش من القائمة المنسدلة التي تظهر. سيعيدك هذا للحظات إلى لوحة الخريطة حيث يجب عليك الآن النقر مع الاستمرار على المؤشر في الزاوية العلوية اليسرى من نطاق الخريطة ثم الاستمرار في الضغط على زر الماوس لأسفل أثناء سحب المؤشر إلى الزاوية اليمنى السفلية من الخريطة مدى. يجب أن تلاحظ مستطيلًا أحمر كبير يتم رسمه فوق الخريطة أثناء سحب المؤشر. بمجرد أن يكون المؤشر فوق الزاوية اليمنى السفلية من الخريطة ، توقف عن الضغط باستمرار على زر الماوس الأيسر لإنهاء تحديد النطاق الخاص بك والعودة إلى نافذة Clip Raster by Extent.
    11. مرة أخرى في نافذة Clip Raster by Extent ، اضبط معلمة ملف التعريف (الموجودة ضمن Advanced Parameters) على ضغط JPEG واضبط قيمة JPEG_QUALITY على 99. بعد ذلك ، انقر فوق الزر & hellip الموجود بجوار الحقل المحول لتحديد Save to File & hellip والذي سيسمح لك لاختيار مكان حفظ الصورة المضغوطة. بعد ذلك ، انقر فوق الزر & hellip الموجود بجوار الحقل & ldquoClipped (مدى) & rdquo لفتح نافذة حفظ الملف حيث يجب عليك استعراض موقع صورة الخريطة ذات المرجعية الجغرافية التي تحاول ضغطها والنقر عليها لتحديدها. سيؤدي هذا تلقائيًا إلى ملء اسم هذه الصورة في حقل اسم الملف. يوصى بتعديل اسم الملف هذا عن طريق إضافة _c إلى نهاية اسم الملف للإشارة إلى أن هذا الإصدار من الصورة مضغوط. أثناء وجودك في نافذة Save File هذه ، ستحتاج أيضًا إلى تغيير قيمة الحقل Save as Type إلى ملفات JPG (* .jpg) بحيث يتم ضغط الملف بشكل صحيح. بمجرد الانتهاء من إدخال هذه المعلمات ، امض قدمًا وانقر فوق تشغيل.
    12. تهانينا ، لقد نجحت الآن في إحداث إشارة جغرافية لصورة خريطة ممسوحة ضوئيًا. إذا كنت ترغب في إنشاء مرجع جغرافي لصورة خريطة أخرى ، فارجع إلى نافذة Georeferencer حيث لا يزال من المفترض أن ترى آخر صورة قمت بمعالجتها. لاستبدال هذا بصورة خريطة جديدة تحتاج إلى مرجع جغرافي ، انقر فوق افتح النقطية في أقصى يسار شريط الأدوات ، استعرض وصولاً إلى موقع الصورة الجديدة التي تريد معالجتها ، ثم انقر فوق "موافق". بعد بضع ثوان من المفترض أن ترى صورة الخريطة الجديدة تحل محل الصورة السابقة التي أكملتها للتو. يمكنك الآن متابعة تحديد الموقع الجغرافي لهذه الخريطة الجديدة عن طريق بدء هذه التعليمات مرة أخرى من الخطوة 4.

    أهداف

    بنهاية دورة GL O MACS التدريبية ، سيكون المشاركون قادرين على:

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

    افتح طبقة المتجهات

    بعد ذلك ، افتح طبقة متجه تحتوي على مدى الاقتصاص الذي تريد استخدامه لاقتصاص بياناتك. لفتح ملف شكل ، يمكنك استخدام وظيفة gpd.read_file () من geopandas.

    بعد ذلك ، استكشف النظام المرجعي للإحداثيات (CRS) لكل من مجموعتي البيانات. تذكر أنه من أجل إجراء أي تحليل باستخدام مجموعتي البيانات هاتين معًا ، يجب أن يكونا في نفس CRS.

    مخطط لطبقة مدى القطع التي ستستخدمها لاقتصاص البيانات النقطية. يمثل الامتداد المكاني لملف الشكل "الحافة" الجغرافية أو الموقع الذي يقع في أقصى الشمال والجنوب الشرقي والغرب. وهكذا يمثل التغطية الجغرافية الشاملة للجسم المكاني. مصدر الصورة: Colin Williams، NEON.

    الآن بعد أن قمت باستيراد ملف الشكل. ارسم الطبقتين معًا لضمان التداخل مع بعضهما البعض. إذا كان ملف الشكل لا يتداخل مع الخطوط النقطية ، فلا يمكنك استخدامه للقص!

    تم تراكب مدى القطع أعلى البيانات النقطية. عندما تقوم بقص البيانات النقطية ، ستتم إزالة جميع البيانات الموجودة خارج مدى القطع.

    لاقتصاص البيانات ، استخدم وظيفة crop_image في earthpy.spatial.

    هنا يمكنك رؤية نتائج وظيفة المحاصيل الخاصة بك.

    شاهد الفيديو: Clipping Raster DEM by a Polygon in QGIS