أكثر

إعادة الإسقاط الجماعي للملفات النقطية في حلقة [Linux]


لدي مجلد به العديد.ascملفات نقطية يرافق كل منها ملف.prjملف. اسم الملف الأساسي هو نفسه لكلا الملفين. ومع ذلك ، يمكن أن يختلف الإسقاط بينهم جميعًا (هم حاليًا في مختلف إسقاط UTM مختلف).
ما أريد تحقيقه هو توحيد إسقاط كل تلك الملفات النقطية الفردية. توجهي الحالي هو هذاسحق-النصي

#! / bin / bash echo بدء التحويل إلى WGS84 لـ f في * .asc؛ القيام بتحويل الصدى "$ f" ؛ صدى إلى "$ {f٪ .asc} _wgs84.tif"؛ gdalwarp "$ f" "$ {f٪ .asc} _wgs84.tif" -t_srs "+ proj = longlat + ellps = WGS84" ؛ انتهى تحويل الصدى

ومع ذلك ، يفشل البرنامج النصي لكل ملف يعود:

"طلب غير صالح لـ pj_param ، فادح"

لذلك أود أن أرى مجموعة من ملفات .tif المعاد إسقاطها (أو أي تنسيق آخر) التي تحتوي جميعها على نفس الإسقاط.

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


يتصل:gdalinfo grid2000.ascيبدو أنه فشل حتى في النهاية؟

برنامج التشغيل: AAIGrid / Arc / Info ASCII Grid Files: grid2000.asc grid2000.prj الحجم 16 ، 16 نظام الإحداثيات هو: PROJCS ["WGS_1984_UTM_Zone_30N" ، GEOGCS ["GCS_WGS_1984" ، DATUM ["WGS_1984" ، SPHEROID ["WGS_1984" ، 6378137،298.257223563]]، PRIMEM ["Greenwich"، 0]، UNIT ["Degree"، 0.017453292519943295]]، PROJECTION ["Transverse_Mercator"]، PARAMETER ["latitude_of_origin"، 0]، PARAMETER ["central_meridian" ، PARAMETER ["scale_factor"، 0.9996]، PARAMETER ["false_easting"، 500000]، PARAMETER ["false_northing"، 0]، UNIT ["Meter"، 1]] الأصل = (565146.201264380943030،680158.85424000048001081) حجم البكسل = ( -30.000000000000000) طلب غير صالح لـ pj_param ، قاتل

call: gdalsrsinfo xyz.prj # (لاحظ أن xyz.prj قد يكون ملفًا مختلفًا)

PROJ.4: '+ proj = utm + zone = 17 + south + datum = WGS84 + وحدة = m + no_defs' OGC WKT: PROJCS ["WGS_1984_UTM_Zone_17S" ، GEOGCS ["GCS_WGS_1984" ، DATUM ["WGS_1984" " WGS_84 "، 6378137،298.257223563]]، PRIMEM [" Greenwich "، 0]، UNIT [" Degree "، 0.017453292519943295]]، PROJECTION [" Transverse_Mercator "]، PARAMETER [" latitude_of_origin "، 0]،" central_of_origin "، 0]" -81]، PARAMETER ["scale_factor"، 0.9996]، PARAMETER ["false_easting"، 500000]، PARAMETER ["false_northing"، 10000000]، UNIT ["Meter"، 1]]

قبل تشغيل البرنامج النصي في نظام دبيان الخاص بي:

تم تشغيله بشكل مثالي دون أي خطأ (وبدون تعديلات في أي مكان):

تم إسقاط كل صورة بشكل مناسب إلى WGS84.

ما هي نتيجة الإسقاط الفردي باستخدام الأمر في وحدة تحكم bash (بدون برنامج نصي)؟ هل هو ناجح؟


شاهد الفيديو: المحاضرة الخامسة GIS. تحديد وتغيير نظام الاسقاط Define Projection, Project, Project Raster (شهر اكتوبر 2021).