أكثر

كيفية تغيير إسقاط ودقة ملف باستخدام gdal؟

كيفية تغيير إسقاط ودقة ملف باستخدام gdal؟


لدي ملف بهذه المعلومات:

lat_min -60.0 lat_max 80.0 lon_min -180.0 lon_max 179.9 lon_sampling 0.1 درجة أخذ عينات lat_max 0.1 درجة = 3600 سطر = 1401

أريد إعادة عرض الملف إلى

حاولت

gdalwarp -of "ENVI" -t_srs 'EPSG: 3410' -ts 1383586 ولكن خطأ gdalwarp -of "ENVI" -t_srs "+ proj = cea + lon_0 = 0 + lat_ts = 30 + x_0 = 0 + y_0 = 0 + a = 6371228 + b = 6371228 + الوحدات = m + no_defs "-ts 1383586 - ملف Float32. envi file_1383586.envi خطأ 1: لا توجد ترجمة PROJ.4 لمصدر SRS ، فشلت تهيئة تحويل الإحداثي

نظرًا لأن البيانات الوصفية للملفات غير مكتملة ، يجب عليك إضافتها يدويًا:

gdal_translate -a_srs epsg: 4326 -a_ullr -180 80179.9 -60 file.envi test.tif

والنتيجة تتناسب بشكل جيد مع ساحل الأرض الطبيعية ، وجاهزة لمزيد من الاعوجاج:

gdalwarp -of "ENVI" -t_srs "+ proj = cea + lon_0 = 0 + lat_ts = 30 + x_0 = 0 + y_0 = 0 + a = 6371228 + b = 6371228 + الوحدات = m + no_defs" -ts 1383586 -ot Float32 test.tif file_1383586.envi


يمكن أن يساعدك هذا الأمر في تحويل إسقاط ملف الشكل عندما يكون لديك مجموعتان مختلفتان من EPSG.

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

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

فيما يتعلق بـ GDAL OGR ، لا توجد مشكلة في الأداة نفسها. يحتوي ملف prj الذي تم إنشاؤه على إسقاط EPSG الذي قمت بتعيينه لملف الشكل. المشكلة في مثل هذه الحالات هي أن الخريطة وطبقة ملف الأشكال ليسا على نفس نظام الإحداثيات الذي يسبب التحول. لحل هذه المشكلة ، تأكد من أن ملف PRJ الذي تم إنشاؤه يحتوي على نظام الإحداثيات المستهدفة والأهم من ذلك هو التأكد من أن إحداثيات الخريطة متطابقة تمامًا.


استخراج المعلومات النقطية باستخدام GDAL و Python

مكتبة تجريد البيانات الجغرافية المكانية (GDAL) هي مكتبة لمعالجة البيانات النقطية. تتم صيانته من قبل مؤسسة الجغرافيا المكانية مفتوحة المصدر (OGF) وعادة ما تأتي مجمعة مع مكتبتها الشقيقة OGR. يسمح OGR بمعالجة بيانات المتجه. إذا كنت تستخدم بالفعل برنامجًا جغرافيًا مكانيًا ، فمن المحتمل أن يكون لديك GDAL مثبتًا بالفعل. يأتي GDAL مع مجموعة من الأدوات المساعدة لسطر الأوامر المفيدة وطريقة واحدة سريعة للتحقق مما إذا كان مثبتًا هو محاولة تشغيل أحدها. ومن الأمثلة الجيدة على ذلك هو gdalinfo الذي يعطي معلومات عن تنسيق الملف وإسقاطه ومدى انتشاره.

متطلبات

إذا لم يكن gdal & # 8217t مثبتًا ، فيمكن تنزيله من http://www.gdal.org/. في الأنظمة المستندة إلى نظام يونكس ، يمكن عادةً تثبيته عبر مدير الحزم الخاص بك. حاول البحث عن gdal أو gdal-bin.

للتفاعل مع GDAL من Python ، تحتاج إلى تثبيت حزمة Python GDAL. يمكن الحصول عليها من http://pypi.python.org/pypi/GDAL/. يتم توفير تعليمات التثبيت والوثائق الأخرى. يمكنك اختبار ما إذا كان مثبتًا عن طريق فتح موجه أوامر python والكتابة ،

استيراد osgeo

إذا لم & # 8217t تحصل على رسالة خطأ ، فكل شيء على ما يرام.

يمكنك أيضًا تثبيت Python و GDAL ومجموعة من الأدوات الجغرافية المكانية الأخرى باستخدام حزمة FWTools ، http://fwtools.maptools.org/. يقال إن هذه هي أسهل طريقة (خاصة إذا كنت تستخدم Windows) ، على الرغم من أنني لم أجربها بنفسي.

إذا لم يكن لديك & # 8217t صورة نقطية مناسبة لتلعب بها ، يمكنك تنزيل البيانات من & # 8220Global Land Cover Facility: Earth Science Data Interface & # 8221 على http://glcfapp.glcf.umd.edu:8080/esdi /index.jsp.

ما الذي يمكن لـ GDAL فعله؟

يتيح لك GDAL القياسي قراءة 81 تنسيقًا للملفات النقطية والكتابة إلى 41 تنسيقًا. يتيح لك الوصول إلى البيانات الموجودة في نطاقات الصور. يمكنك قراءة هذه البيانات في مصفوفة وتنفيذ معالجة الصور. يتيح لك GDAL أيضًا الوصول إلى المعلومات التالية:

  • تنسيق البيانات
  • عدد الصفوف في الصورة
  • عدد الأعمدة في الصورة
  • عدد النطاقات في الصورة
  • معلومات الإسقاط
  • إحداثيات X و Y في الزاوية اليسرى العليا من الصورة
  • دقة البكسل بين الشمال والجنوب
  • دقة البكسل بين الغرب والشرق
  • أي استدارة تم تطبيقها على الصورة

يمكن العثور على مزيد من التفاصيل وبعض أمثلة بايثون على http://www.gdal.org/gdal_tutorial.html.

تتيح لك مكتبة python GDAL أيضًا معالجة جوانب أخرى من الصورة مثل جدول الألوان. ستمنحك كتابة help (gdal) في موجه أوامر python مزيدًا من المعلومات حول الأساليب والسمات التي يدعمها كائن GDAL.

GDAL و Python

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

بادئ ذي بدء ، تحتاج إلى تحميل مكتبة gdal

تهيئة gdal وفتح ملف


"" حزمة osgeo تحتوي على مكتبات GDAL و OGR و OSR "" "

من osgeo import gdal، osr

فأنت بحاجة إلى فتح ملف مناسب ، في هذه الحالة & # 8220landsatETM.tif & # 8221. في بيثون ، يكتشف GDAL نوع الملف الذي تفتحه. يقوم بتحميل وتسجيل برنامج التشغيل الصحيح لك.

"" أنشئ كائن gdal وافتح ملف "LandsatETM.tif"

datafile = gdal.Open ("landatETM.tif")

استخراج المعلومات

الآن يمكننا استخراج بعض المعلومات من الملف. بادئ ذي بدء ، سنحصل على أبعاد الصورة وعدد النطاقات

الحصول على تفاصيل حجم الصورة


cols = ملف البيانات
الصفوف = ملف البيانات
العصابات = ملف البيانات

"" "اطبع المعلومات على الشاشة. تحويل الأرقام المعادة إلى السلاسل باستخدام str ()" ""

طباعة "عدد الأعمدة" + str (cols)
طباعة "عدد الصفوف" + str (صفوف)
طباعة "عدد النطاقات:" + str (نطاقات)

الحصول على موضع الزاوية XY

يمكننا أيضًا الحصول على الموضعين X و Y في الزاوية اليسرى العليا. أولاً ، نسمي طريقة GetGeoTransfrom التي تُرجع قائمة. يمكننا الوصول إلى عناصر القائمة للحصول على التفاصيل التي نريدها. يمكن العثور على مزيد من المعلومات حول قيم القائمة الأخرى على http://www.gdal.org/gdal_tutorial.html.


"" "نسمي أولاً طريقة GetGeoTransform لكائن ملف البيانات" ""

geoinformation = ملف البيانات. GetGeoTransform ()

"" "إحداثيات X و Y العلوية اليسرى في قائمة المواضع 0 و 3 على التوالي" ""

topLeftX = معلومات جغرافية [0]
topLeftY = معلومات جغرافية [3]

"" "طباعة هذه المعلومات إلى الشاشة" ""

طباعة "أعلى اليسار X:" + str (topLeftX)
طباعة "أعلى اليسار Y:" + str (topLeftY)

الحصول على نظام الإسناد المكاني

أخيرًا يمكننا الحصول على معلومات حول نظام الإسناد المكاني. يتم إرجاع هذا بتنسيق نص معروف (WKT). هذا تنسيق شائع يستخدم لتبادل المعلومات المكانية. يتم تعريفه من قبل Open Geospace Consortium. تنسيق WKT مشابه للإخراج الذي تحصل عليه من تشغيل gdalinfo في سطر الأوامر. يعد تنسيق WKT مفيدًا جدًا ولكن بالنسبة لبعض التطبيقات ، فأنت تحتاج إلى تنسيق proj4 الذي تستخدمه مكتبة عرض الخرائط proj4 (تأكد من ظهوره في المنشورات اللاحقة). يمكننا تحويل تنسيق WKT إلى تنسيق proj4 بسهولة بالغة.

"" نصل أولاً إلى معلومات الإسقاط داخل ملف البيانات الخاص بنا باستخدام طريقة GetProject (). هذا يعيد سلسلة بتنسيق WKT "" "

"" "ثم نستخدم وحدة osr التي تأتي مع GDAL لإنشاء كائن إسناد مكاني" ""

"" "نستورد سلسلة WKT الخاصة بنا إلى spatialRef" ""

"" "نستخدم طريقة ExportToProj4 () لإرجاع سلسلة إسناد مكاني بنمط proj4." ""

"" "يمكننا بعد ذلك طباعتها" ""

طباعة "WKT format:" + str (spatialRef)
طباعة "تنسيق Proj4:" + str (spatialRefProj)

التعليقات النهائية

لذلك ، فتحنا مجموعة بيانات نقطية واستخرجنا معلومات حول الملف. هذا مفيد في حد ذاته حيث يمكنك إنشاء بعض أدوات سطر الأوامر المخصصة البسيطة جدًا المشابهة لـ gdalinfo. يمكن العثور على معلومات حول الوصول إلى وسيطات سطر الأوامر على http://diveintopython.org/scripts_and_streams/command_line_arguments.html

تم إبقاء المثال بسيطًا عن عمد ، لذلك تم استبعاد أشياء مثل التحقق من الأخطاء

في المرة القادمة سأنظر في تحميل بيانات الصورة في مصفوفة ، وعرض صورة على الشاشة وربما بعض التلاعب الأساسي بالصور باستخدام مكتبات Python و Python Imaging العلمية.


GDAL و OGR: أدوات تحويل البيانات الجغرافية وإعادة الإسقاط

تكون العمليات على البيانات الجغرافية أكثر فاعلية عندما يكون لملفات الإدخال معلمات مكانية متطابقة: أي نظام الإحداثيات والمرجع / الإليسكوبويد. يمكن استخدام أدوات ومكتبات GDAL / OGR لـ "إعادة طرح" البيانات الجغرافية التي تم تنزيلها من نظام إحداثي إلى آخر ، لكل من بيانات المتجه والبيانات النقطية. يمكن العثور على مزيد من المعلومات حول هذا البرنامج هنا: http://www.remotesensing.org/gdal/index.html توضح الأمثلة التالية بعض الطرق التي يمكن من خلالها استخدام أدوات GDAL / OGR لمعالجة المعلمات المكانية لبعض الأنواع الشائعة من البيانات الجغرافية: SSURGO (تنسيق ESRI Shapefile) والصور الجوية (تنسيق GeoTiff).

كيف تقارن مع ArcToolbox؟

تم استخدام مضلع حدود SSURGO في اختبار سريع لمعرفة مدى تشابه نتائج إعادة إسقاط متجه بواسطة ArcMap و OGR. تم إسقاط متجه الإدخال (GCS NAD83) إلى منطقة UTM 11 NAD83 بواسطة كل من OGR و ArcMap 9.0 ، ثم تم استيراده إلى GRASS. ثم تم تصدير كلا المتجهين كـ ascii ، للوصول البسيط إلى إحداثيات قمة الخط. تم استخدام برنامج نصي AWK بسيط لمقارنة الاختلافات عدديًا في إحداثيات X و Y:

كانت قيم متوسط ​​dx و dy الناتجة هي: DX: 2.98838e-06 م و دى: -1.23948e-07 م مع 5035 العينات: تبدو جيدة بالنسبة لي. في المرة القادمة ، قد يكون تغيير قاعدة البيانات أكثر وضوحًا لأي اختلافات دقيقة يجب أن تكون على دراية بها. تحديث: NAD83 - & gt NAD27 تجربة إزاحة Datum تبدو كما لو أن تحولات الإسناد قابلة للمقارنة أيضًا.

أمثلة:
توجد أسماء الأوامر بتنسيق جريئة اكتب ، وملفات الإدخال باللون الأحمر ، وملفات الإخراج باللون الأخضر. يتم توفير الأرقام لسياق إضافي ، ويتم الرجوع إليها من خلال اسم ملف الإدخال / الإخراج في رمز المثال.

عمليات ناقلات

الحصول على معلومات حول ملف الشكل:

إعادة طرح بيانات وحدة خرائط SSURGO من الإحداثيات الجغرافية (خطوط الطول / العرض) إلى منطقة UTM 10:

استخرج جميع المضلعات من ملف أشكال SSURGO حيث يكون رمز mapunit هو "ScA":

تحويل ملف شكل إلى KML متوافق مع Google:

العمليات النقطية [برنامج تعليمي نقطي]

احصل على معلومات حول مجموعة البيانات النقطية

أعد إسقاط صورة جوية في منطقة المستوى 4 بولاية كاليفورنيا (إسقاط لامبرت المخروطي المطابق ، الوحدات = قدم) إلى UTM Zone 10 (الوحدات = متر) ، وأعد القياس إلى دقة إخراج تبلغ 1 متر ، واستخدم إعادة تشكيل شريحة رقيقة من الشرائح (-tps) :

لاحظ أنه من الجيد عادة "تحسين" الصورة الناتجة باستخدام gdal_translate.
تحويل ملف GeoTiff متعدد النطاقات إلى JPEG:

قراءة ملفات بتنسيق MrSid تفاصيل

أضف مسار DSDK للتكوين: - مع mrsid = / Geo_DSDK-6.0.7.1407 /
يوجد خطأ في Geo_DSDK-6.0.7.1407 // include / base / lti_sceneBuffer.h line 356:

علق على هذا السطر وأعد تجميعه. شكراً لـ Mateusz Loskot على قناة GDAL IRC للإشارة إلى ذلك.


إعادة طرح البيانات واختزالها¶

أوضح القسم السابق كيفية حفظ البيانات النقطية. ومع ذلك ، في كثير من الحالات ، هناك & # 8217s حاجة إلى إعادة طرح هذه البيانات وإعادة تشكيلها. الحل العملي يستخدم gdalwarp ويفعل ذلك على الغلاف. من ناحية ، يعد هذا مناسبًا ، ولكن في بعض الأحيان ، تحتاج إلى تنفيذ هذه المهمة كخطوة وسيطة ، ويكون إنشاء الملفات وحذفها مملاً وعرضة للخطأ. من الناحية المثالية ، سيكون لديك وظيفة بيثون من شأنها أن تؤدي الإسقاط نيابة عنك. يسمح GDAL بهذا من خلال تعريف الملفات النقطية في الذاكرة. هذه مجموعات بيانات GDAL عادية ، ولكن لا توجد على نظام الملفات ، فقط في ذاكرة الكمبيوتر & # 8217s. إنها لوحة مفاتيح مريحة & # 8220scratchpad & # 8221 لإجراء حسابات وسيطة سريعة. يوفر GDAL أيضًا وظيفة ، gdal.ReprojectImage الذي يفضح معظم قدرات gdalwarp. سنقوم بدمج هاتين الحيلتين لتنفيذ إعادة الإسقاط. على سبيل المثال ، سننظر في الحالة التي تحتاج فيها بيانات NDVI الخاصة بالجزر البريطانية المذكورة في القسم السابق إلى إعادة إسقاطها إلى الشبكة الوطنية لمسح الذخائر ، وهو توقع مناسب للمملكة المتحدة.

المضاعفات الرئيسية تأتي من الحاجة gdal.ReprojectImage للعمل على مجموعات بيانات GDAL. في القسم السابق ، قمنا بحفظ بيانات NDVI في ملف GeoTIFF ، وهذا يعطينا مجموعة بيانات البداية. ما زلنا بحاجة إلى إنشاء مجموعة بيانات الإخراج. هذا يعني أننا بحاجة إلى تحديد التحويل الجغرافي وحجم مجموعة بيانات الإخراج قبل إجراء الإسقاط. يستلزم ذلك جمع المعلومات حول مدى مجموعة البيانات الأصلية ، وإسقاطها على إسقاط الوجهة ، وحساب عدد وحدات البكسل ومعلمات التحويل الجغرافي من هناك. هذه وظيفة (تم التعليق عليها بشدة) تؤدي هذه المهمة فقط:

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

نقوم بتوسيع الجزء الرئيسي من البرنامج إلى (1) حفظ نتيجة إعادة الإسقاط كملف GeoTIFF ، (2) قراءة ملف البيانات الناتج و (3) رسمه:


قراءة متعمقة

  • ناسا DSCOVR: صفحة EPIC مع مناظر في الوقت الحقيقي للأرض المضاءة بنور الشمس.
  • مختبر مراقبة الأرض التابع لـ JSC للحصول على رؤية أقرب.
  • مجموعة خرائط David Rumsey لمنظور تاريخي.
  • Natural Earth للحصول على البيانات النقطية والأمبيرية الشاملة.
  • دليل الخرائط المواضيعية لخرائط المحور. . . .
  • حتى ورشة عمل Rouault's GDAL.
  • XKCD على خريطة الإسقاطات.

بفضل Frank Warmerdam و Kelsey Jordahl و Joe Kington للدعم الفني والتحرير.


24.2.5.2. استخراج الإسقاط¶

يستخرج إسقاط ملف نقطي ويكتبه في ملف العالمية ملف بملحق .wld.

هذه الخوارزمية مشتقة من الأداة المساعدة GDAL srsinfo.

القائمة الافتراضية: النقطية ► الإسقاطات

24.2.5.2.1. حدود¶

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

قم أيضًا بإنشاء ملف .prj

إذا تم تنشيط هذا ، فسيتم أيضًا إنشاء ملف .prj يحتوي على معلومات الإسقاط.

24.2.5.2.2. النواتج¶

ملف نصي بامتداد .wld يحتوي على معلمات التحويل لملف البيانات النقطية.

ESRI بحلقات ملف prj

ملف نصي بامتداد .prj يصف ملف CRS. سيكون بلا إذا كان إنشاء ملف .prj أيضًا خطأ.

24.2.5.2.3. كود بايثون¶

معرف الخوارزمية: gdal: extractproject

ال معرف الخوارزمية يتم عرضه عندما تحوم فوق الخوارزمية في صندوق أدوات المعالجة. ال قاموس المعلمات يوفر المعلمة NAMEs والقيم. راجع استخدام خوارزميات المعالجة من وحدة التحكم للحصول على تفاصيل حول كيفية تشغيل خوارزميات المعالجة من وحدة تحكم Python.


تحويل ملفات الشبكة

قدم PROJ 7 التنزيل عند الطلب (أجزاء) من شبكات التحويل من شبكة توصيل المحتوى إلى دليل قابل للكتابة بواسطة المستخدم على مسار بحث PROJ (عادةً ما يكون مكون المسار الأول). يتم تخزين حالة الشبكات التي تم تنزيلها في قاعدة بيانات SQLite أخرى ، cache.db. يتم تعيين دليل CDN القابل للكتابة بواسطة المستخدم PR ( phi ) J بمجرد الاستعلام عن مسار البحث الداخلي ، وبالنسبة لمعظم الاستخدامات ، ستسمح القيمة الافتراضية لجميع البرامج التي تستخدم PR ( phi ) J مثل R الحزم ، QGIS ، GRASS ، إلخ ، للوصول إلى أي شبكات تم تنزيلها. يتم فحص الشبكات للتأكد من خلوها من التلف على فترات منتظمة. يمكن تعيين هذا الدليل على قيمة غير افتراضية باستخدام متغير البيئة PROJ_USER_WRITABLE_DIRECTORY قبل رجدال (وأي حزمة أخرى تستخدم PR ( phi ) J) يتم تحميلها وإرفاقها ، من PR ( phi ) J & gt = 7.1.0. يتم هنا تكرار الرمز المستخدم في بداية هذه المقالة القصيرة للتسهيل:

دعونا نتحقق من ذلك رجدال قيد التشغيل مع تعطيل التنزيل عند الطلب:

يمكننا تمكين التنزيل عند الطلب باستخدام وظيفة تُبلغ عن قيمة دليل المستخدم القابل للكتابة في PR ( phi ) J CDN ، ونرى أنه باستخدام هذا الإعداد ، يتم تمكين تنزيل الشبكة:

سيتذكر القارئ أن مسار البحث تم تخزينه في shpr أعلاه ، والعنصر الأول هو الدليل القابل للكتابة بواسطة المستخدم في هذه المرحلة ، ولم يتم إنشاء قاعدة بيانات SQLite:

عندما نبحث بعد ذلك عن عمليات تنسيق المرشح ، نرى أن العملية التي تستخدم شبكة غائبة ترى الآن أن التنزيل ممكّن ، وتقترح دقة 1 مليون مرشح ، لأنه يمكن تنزيل الشبكة المطلوبة (مرة أخرى باستخدام منطقة الاهتمام):

عند إجراء التحويل ، قد نرى أن عملية الإحداثيات تستغرق وقتًا أطول من المتوقع ، لأنه عند المرور لأول مرة يتم تنزيل الشبكة من الشبكة:

تحدد عملية الإحداثيات المستخدمة الآن الشبكة في خط الأنابيب:

تختلف قيم الإحداثيات قليلاً عن خط أنابيب Helmert بدقة 2 متر:

كما نرى ، فإن نقطة الدقة 1 م هي 1.7 م من نقطة الدقة 2 م ، فقط إلى الغرب:

الآن تم إنشاء قاعدة بيانات ذاكرة التخزين المؤقت لشبكة SQLite وازداد حجمها:

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

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

يتم عرض مواضع النتيجة هنا عند التكبير / التصغير 18 ، يمكننا أن نرى أن النقطة الخضراء بدقة 1 متر تتطابق مع موقع Open Street Map الخاص بالمضخة جيدًا:


أمثلة¶

على سبيل المثال ، مشهد موضعي من ثماني بتات مخزّن في GeoTIFF مع نقاط تحكم ترسم الزوايا إلى خط العرض / الطول يمكن تشويهها إلى إسقاط UTM بأمر مثل هذا:

على سبيل المثال ، القناة الثانية لصورة ASTER المخزنة في HDF مع نقاط التحكم التي تعين الزوايا إلى خط العرض / الطول يمكن أن تتشوه إلى إسقاط UTM بأمر مثل هذا:

لتطبيق خط قطع على صورة غير محددة جغرافيًا ومقطع من بكسل (220،60) إلى بكسل (1160،690):

حيث يكون محتوى cutline.csv مثل:

لتحويل DEM من ارتفاعات الجيود (باستخدام EGM96) إلى ارتفاعات WGS84 الإهليلجية:


تلاحظ

ال r.in.gdal الأمر يدعم الميزات التالية ، طالما أن برنامج تشغيل التنسيق الأساسي يدعمها:

ستعمل فرق جدول الألوان مع العناصر الملونة المرتبطة على نقل جداول الألوان. لاحظ أنه في حالة عدم احتواء المصدر على خريطة ألوان ، فإن r.in.gdal في GRASS 5.0 لن يصدر أي مخطط ألوان. استخدم خريطة r.colors =. color = رمادي لتعيين مخطط ألوان بتدرج الرمادي. في إصدار مستقبلي من GRASS r.in.gdal ، من المحتمل أن تتم ترقيته ليصدر تلقائيًا خرائط ألوان بتدرج الرمادي.
أنواع البيانات يتم دعم معظم أنواع بيانات GDAL. يتم ترجمة نطاقي النوع Float32 و Float64 كخلايا فاصلة عائمة GRASS (ولكن ليس الدقة المزدوجة. يمكن إضافة هذا إذا لزم الأمر) ، ويتم ترجمة معظم الأنواع الأخرى كخلايا عدد صحيح GRASS. يتضمن ذلك مصادر بيانات ذات عدد صحيح 16 بت. تتم ترجمة نطاقات البيانات المعقدة (بعض تنسيقات بيانات إشارة SAR) إلى طبقتين من الخلايا النقطية العائمة (* .real و * .imaginary).
الإسناد الجغرافي إذا كانت مجموعة البيانات تحتوي على معلومات الإسناد الجغرافي ، فسيتم استخدام ذلك لتعيين الحواف الشمالية والجنوبية والشرقية والغربية. سيتم تجاهل معاملات الدوران ، مما يؤدي إلى تحديد موضع غير صحيح لمجموعات البيانات التي تم تدويرها.
الإسقاط سيتم استخدام إسقاط مجموعات البيانات للمقارنة بالموقع الحالي أو لتحديد موقع جديد. داخليًا ، يمثل GDAL الإسقاطات بتنسيق OpenGIS Well Known Text. يتم دعم مجموعة فرعية كبيرة من المجموعة الإجمالية لإسقاطات GRASS.
Null Values ​​النطاقات النقطية التي يتعرف GDAL على قيمة خالية لها ستتحول إلى وحدات بكسل فارغة إلى نمط GRASS فارغ أثناء الاستيراد. لا تملك العديد من التنسيقات العامة (والتنسيقات المدعومة بشكل سيئ بواسطة GDAL) طريقة للتعرف على وحدات البكسل الفارغة وفي هذه الحالة يجب استخدام r.null بعد الاستيراد.
مجموعات بيانات GCPs التي تحتوي على نقاط تحكم أرضية سيتم استيرادها كملف POINTS مرتبط بمجموعة الصور. مجموعات البيانات التي تحتوي على نطاق واحد فقط والتي كان من الممكن ترجمتها كخريطة نقطية بسيطة سيكون لها أيضًا مجموعة صور مرتبطة إذا كانت هناك نقاط تحكم أرضية. تم الإبلاغ عن نظام إحداثيات نقاط التحكم الأرضية بواسطة r.in.gdal ولكن لم يتم حفظه. يعود الأمر للمستخدم للتأكد من أن الموقع الذي تم إنشاؤه باستخدام i.target به نظام إحداثيات متوافق قبل استخدام النقاط مع i.rectify.
جداول البيانات النقطية r.in.gdal يمكن كتابة جداول البيانات النقطية كملفات CSV. علاوة على ذلك ، يتم استيراد المعلومات الموجودة في جداول البيانات النقطية تلقائيًا طالما أن تعريفات الحقول تحتوي على معلومات حول كيفية استخدام حقل ، على سبيل المثال للحصول على معلومات اللون أو التسميات.

التحسينات المخططة ل r.in.gdal في المستقبل ، قم بتضمين دعم للإبلاغ عن كل شيء معروف عن مجموعة البيانات إذا كان انتاج لم يتم تعيين المعلمة.

رسائل خاطئة

"خطأ: لا يبدو أن إسقاط مجموعة البيانات يطابق الموقع الحالي."
تحتاج إلى إنشاء موقع يتطابق عرضه مع البيانات التي ترغب في استيرادها. جرب استخدام موقعك المعلمة لإنشاء موقع جديد بناءً على معلومات الإسقاط في الملف. إذا رغبت في ذلك ، يمكنك إعادة عرضه على موقع آخر باستخدام r.proj. بدلاً من ذلك ، يمكنك تجاوز هذا الخطأ باستخدام ملف -o علم.

"تحذير: G_set_window (): خط عرض غير قانوني للشمال"
لا يمكن أن تحتوي مواقع خطوط الطول / العرض في GRASS على مناطق تتجاوز 90 درجة شمالًا أو جنوبًا. ستحتوي الصور غير المحددة جغرافيًا على إحداثيات بناءً على عدد وحدات البكسل للصور: 0،0 في الأعمدة اليسرى السفلية ، والصفوف في أعلى اليمين. عادةً ما يكون طول الصور أكثر من 90 بكسل ولذلك يرفض نظام المعلومات الجغرافية استيرادها. إذا كنت متأكدًا من أن البيانات مناسبة لموقع Lat / Lon الخاص بك وتريد إعادة تعيين حدود الخريطة باستخدام ص المنطقة الوحدة النمطية مباشرة بعد الاستيراد ، يمكنك استخدام علامة لتقييد إحداثيات الخريطة بالقيم القانونية. بينما من المحتمل أن تكون الحدود والدقة الناتجة خاطئة لخريطتك ، فإن بيانات الخريطة ستكون آمنة وغير قابلة للتغيير. بعد إعادة الضبط إلى الحدود المعروفة بـ ص المنطقة يجب عليك التحقق منها مرة أخرى r.info، مع إيلاء اهتمام خاص لدقة الخريطة. في معظم الحالات ، سترغب في الاستيراد إلى الإسقاط الأصلي لملف البيانات ، أو إلى موقع XY بسيط واستخدام أدوات التحليل الجغرافي (ط. تصحيح et al.) لإسقاط الموقع بشكل صحيح في الموقع المستهدف. ال يجب العلم فقط تُستخدم إذا كنت تعرف أن الإسقاط صحيح ولكن الإسناد الجغرافي الداخلي قد ضاع ، وأنت تعرف ما يجب أن تكون عليه حدود الخريطة ودقتها مسبقًا.


أمثلة

مدى الطباعة ودقة البيانات النقطية ثنائية وثلاثية الأبعاد

ز المنطقة-mp3 سيؤدي ذلك إلى طباعة المنطقة الحالية والمنطقة ثلاثية الأبعاد (المستخدمة مع وحدات البكسل) بالتنسيق:

ز المنطقة-ز ال -g الخيار يطبع المنطقة بتنسيق نمط البرنامج النصي التالي (مفتاح = القيمة):

ز المنطقة -bg ال -bg الخيار يطبع المنطقة بتنسيق نمط البرنامج النصي التالي (مفتاح = القيمة) بالإضافة إلى مربع الحدود في خط الطول والعرض / WGS84:

ز المنطقة -l ال الخيار يطبع المنطقة بالتنسيق التالي:

ز المنطقة - مساء سيؤدي هذا إلى طباعة المنطقة الحالية بالتنسيق (موقع خطوط الطول والعرض):

تغيير المدى ودقة البيانات النقطية باستخدام القيم

ز المنطقة n = 51: 36: 05 شمالاً e = 10: 10: 05 شرقًا = 51: 29: 55 شمالًا w = 9: 59: 55 درجة شرقًا = 0: 00: 01 سيعيد ضبط الاتجاه الشمالي ، والشرق ، والجنوب ، والغرب ، والدقة للمنطقة الحالية ، هنا في نمط خطوط الطول والعرض DMS (يمكن أيضًا استخدام الدرجات العشرية والدرجات مع الدقائق العشرية).

ز المنطقة -dp s = 698000 سيقوم بتعيين المنطقة الحالية من المنطقة الافتراضية لموقع قاعدة بيانات GRASS ، وإعادة تعيين الحافة الجنوبية إلى 698000 ، ثم طباعة النتيجة.

ز المنطقة n = n + 1000 w = w-500 ن =القيمة يمكن أيضًا تحديدها كدالة لقيمتها الحالية: n = n +القيمة يزيد من اتجاه الشمال الحالي ، بينما n = n-القيمة يقلل منه. هذا صحيح أيضًا بالنسبة لـ s =القيمة، البريد =القيمةو ث =القيمة. في هذا المثال ، يتم تمديد الحد الشمالي للمنطقة الحالية بمقدار 1000 وحدة ويتم تقليل الحد الغربي للمنطقة الحالية بمقدار 500 وحدة.

ز المنطقة n = s + 1000 e = w + 1000 يسمح هذا النموذج للمستخدم بتعيين قيم حدود المنطقة بالنسبة لبعضها البعض. هنا ، يتم تعيين إحداثيات الحدود الشمالية مساوية لـ 1000 وحدة أكبر من قيمة إحداثيات الحد الجنوبي ، ويتم تعيين قيمة إحداثيات الحد الشرقي مساوية لـ 1000 وحدة أكبر من قيمة إحداثيات الحد الغربي. الأشكال المقابلة s = n-القيمة و

ث = ه-القيمة يمكن استخدامها لتعيين قيم الحدود الجنوبية والغربية للمنطقة ، بالنسبة لقيم الحدود الشمالية والشرقية.

تغيير المدى ودقة البيانات النقطية باستخدام الخرائط

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

ز - خطوط المسح الإقليمية = التربة ال -u الخيار يمنع إعادة ضبط تعريف المنطقة الحالية. يمكن أن يكون هذا مفيدًا عندما يكون مطلوبًا لاستخراج معلومات المنطقة فقط. في هذه الحالة ، تتم طباعة ملف رأس الخلية لطبقة خريطة التربة دون تغيير إعدادات المنطقة الحالية.

ز- التكبير الإقليمي = حفظ التربة = التربة سيؤدي هذا إلى تكبير المنطقة الأصغر التي تشمل جميع قيم بيانات التربة غير الفارغة ، وحفظ إعدادات المنطقة الجديدة في ملف ليتم استدعاؤه التربة وتخزينها تحت دليل windows في مجموعة الخرائط الحالية للمستخدم. لم يتم تغيير إعدادات المنطقة الحالية.

تغيير المدى ودقة البيانات النقطية في الأبعاد الثلاثية

استخدام g.region في قذيفة بالاشتراك مع OGR

استخدام g.region في غلاف بالاشتراك مع GDAL


شاهد الفيديو: ArcMap - تغير أو تعريف نظام إحداثيات شيب فايل الى نظام آخر Define Projection u0026 Project