img2grdgmt - عبر الإنترنت في السحابة

برنامج:

اسم

img2grd - استخراج مجموعة فرعية من ملف img بتنسيق Mercator أو Geographic

موجز

img2grd ملف grdfile منطقة نوع [] [[مينلات / ماكسلات]] [] [ دقائق ] [] [ التنقل
] [[مقياس]] [[مستوى]] [ ماكسلون ] [ -n]

ملحوظة: لا يسمح بمسافة بين علامة الخيار والوسيطات المرتبطة بها.

الوصف

img2grd يقرأ ملف تنسيق img ويستخرج مجموعة فرعية ويكتبها في ملف شبكة. ال -M
يحدد الخيار ما إذا كان إسقاط Spherical Mercator لملف img هو الإسقاط أم لا
تم الاحتفاظ بها أو إذا كان يجب كتابة شبكة جغرافية عن طريق التراجع عن إسقاط مركاتور. لو
يتم تحديد الشبكة الجغرافية ، يمكنك أيضًا طلب إعادة تشكيلها على وجه الدقة -R معطى.

مطلوبة الحجج

ملف
ملف تنسيق img مركاتور مثل حقول الجاذبية البحرية أو تضاريس قاع البحر
المقدرة من بيانات مقياس الارتفاع عبر الأقمار الصناعية بواسطة Sandwell و Smith. إذا قام المستخدم بتعيين
متغير البيئة GMT_DATADIR دولار، ثم img2grd سيحاول العثور عليه ملف in
GMT_DATADIR دولار؛ وإلا فإنه سيحاول فتح ملف مباشرة.

-Ggrdfile
grdfile هو اسم ملف الشبكة الناتج.

-R [وحدة]غرب/الشرق/جنوب/شمالي[/zmin/com.zmax] [ص]
غرب, الشرق, جنوبو شمالي تحديد المنطقة محل الاهتمام ، ويمكنك تحديد
عليها بالدرجات العشرية أو بتنسيق [+ -] dd: mm [: ss.xxx] [W | E | S | N]. ألحق r إذا كانت أقل
يتم إعطاء إحداثيات الخريطة اليمنى واليسرى العلوية بدلاً من w / e / s / n. الاثنان
الاختصارات -Rg و -بحث وتطوير تقف للمجال العالمي (0/360 و -180 / + 180 في خط الطول
على التوالي ، مع -90 / + 90 في خط العرض). بدلا من ذلك لإنشاء الشبكة ، أعط
Rكودلون/اللات/nx/ny، حيث الكود عبارة عن تركيبة مكونة من حرفين من L و C و R (لليسار ،
الوسط أو اليمين) و T أو M أو B للأعلى أو الوسط أو الأسفل. على سبيل المثال ، BL لأسفل اليسار.
يشير هذا إلى أي نقطة على منطقة مستطيلة الشكل LON/اللات تنسيق يشير
و أبعاد الشبكة nx و ny مع تباعد الشبكة عبر -I يستخدم لخلق
المنطقة المقابلة. بدلاً من ذلك ، حدد اسم ملف الشبكة الموجود
و -R يتم نسخ الإعدادات (وتباعد الشبكة ، إن أمكن) من الشبكة.
باستخدام -Rوحدة تتوقع الإحداثيات المتوقعة (الديكارتية) المتوافقة مع المختار -J
ونقوم عكسيًا بتحديد المنطقة الجغرافية المستطيلة الفعلية. ل
وجهة نظر (-p) ، قم بإلحاق اختياريًا /zmin/com.zmax. في حالة عرض المنظور
(-p) ، نطاق z (zmin, com.zmax) للإشارة إلى البعد الثالث. هذه
يجب القيام به فقط عند استخدام ملف -جيز الخيار ، وليس عند استخدام فقط -p الخيار.
في الحالة الأخيرة ، يتم رسم عرض منظور للمستوى ، بدون ثلث
البعد.

اختياري الحجج

-C عيّن إحداثيات x و y Mercator المتعلقة بمركز الإسقاط [الافتراضي هو
نسبة إلى الزاوية اليسرى السفلية للشبكة]. يتطلب -M.

-د[مينلات / ماكسلات]
استخدم نطاق خط العرض الممتد -80.738 / + 80.738. بدلاً من ذلك ، قم بإلحاق
مينلات / ماكسلات كمدى خط العرض لملف img المدخلات. [الافتراضي هو
-72.006 / 72.006]. ليس مطلوبًا عادةً نظرًا لأنه يمكننا تحديد المدى من
فحص حجم الملف.

-E يمكن استخدامها عندما -M لم يتم تعيينه لإجبار الشبكة النهائية على أن يكون لها نفس الشيء تمامًا
المنطقة كما هو مطلوب مع -R. بشكل افتراضي ، المنطقة النهائية هي إسقاط مباشر لـ
منطقة مركاتور الأصلية وستمتد عادةً إلى ما بعد نطاق
نطاق خط العرض المطلوب ، وعلاوة على ذلك ، فإن الزيادة في الشبكة في خط العرض لا تفعل ذلك
تطابق زيادة خط الطول. ومع ذلك ، فإن إعادة التشكيل الإضافية تقدم صغيرة
أخطاء الاستيفاء ويجب استخدامها فقط إذا كان يجب أن تتطابق شبكة الإخراج مع
المنطقة المطلوبة ولها x_inc = y_inc. في هذه الحالة المنطقة التي حددتها -R يجب أن تكون
تعطى بمضاعفات الزيادة (.eg ، -R0/45/45/72).

-I تشير دقائق كعرض إدخال img بكسل في دقائق من خط الطول.
[الافتراضي هو 2.0]. ليس مطلوبًا عادةً نظرًا لأنه يمكننا تحديد حجم البكسل من
فحص الحجم.

-M إخراج شبكة مركاتور كروية [الافتراضي هو شبكة جغرافية / خطوط الطول]. ال
يتم الاحتفاظ بإسقاط Spherical Mercator لملف img ، بحيث يتم الاحتفاظ بالمنطقة -R
تم تعديله من قبل المستخدم بشكل طفيف ؛ المنطقة المعدلة تتوافق مع الحواف
عدد البكسل [أو مجموعات من التنقل بكسل]. يتم تعيين رأس ملف الشبكة بحيث يتم تعيين x و
تمثل أطوال المحور y المسافة من الحواف الغربية والجنوبية للصورة ،
تقاس بالوحدات الافتراضية للمستخدم ، مع -جم1 والمعدلة -R. عن طريق تحديد
الافتراضي مشروع_ بيضاوي = المجال ، يمكن للمستخدم عمل تراكبات مع التعديل -R
حتى يتطابقوا. يرى أمثلة أقل. المعدل -R هو مكتوب أيضا في الشبكة
ملاحظة العنوان ، لذلك يمكن العثور عليها لاحقًا. يرى -C لتعيين الإحداثيات بالنسبة إلى
مركز الإسقاط.

-Nالتنقل متوسط ​​القيم في إدخال img بكسل إلى التنقل by التنقل المربعات وخلق
بكسل إخراج واحد لكل مربع من هذا القبيل. إذا استخدمت مع -T3 سيبلغ عن المتوسط
القيد بين 0 و 1. إذا تم استخدامه مع -T2 سيكون الناتج متوسط ​​قيمة البيانات
أو NaN وفقًا لما إذا كان متوسط ​​القيد> 0.5. التنقل يجب أن تقسم بالتساوي
في أبعاد ملف img بالبكسل. [تقصير 1 لا متوسط].

-س[مقياس]
اضرب قيم ملف img في مقياس قبل التخزين في ملف الشبكة. [الافتراضي هو
1.0]. لملفات img الحديثة: يتم تخزين ملفات img topo في (مصححة) متر [-S1] ؛
ملفات الجاذبية الجوية الحرة في mGal * 10 [-S0.1 للحصول على mGal] ؛ ملفات الانحراف العمودي بتنسيق
ميكرو راديان * 10 [-S0.1 للحصول على ميكرو راديان] ، ملفات التدرج الجاذبية العمودية بتنسيق
يوتفوس * 50 [-S0.02 للحصول على Eotvos أو -S0.002 للحصول على mGal / km]). إذا كان الجواب لا مقياس معطى
نحاول تحديد المقياس من خلال فحص اسم الملف بحثًا عن أدلة.

-Tنوع نوع يعالج ترميز معلومات القيد. نوع = 0 يشير إلى أن لا
يتم ترميز هذه المعلومات في ملف img (المستخدم لإصدارات ما قبل 1995 من
بيانات الجاذبية) ويحصل على جميع البيانات. نوع > 0 يشير إلى أن معلومات القيد هي
المشفرة (1995 والإصدارات الأحدث (الحالية) من ملفات img) بحيث يمكن للمرء أن
إنتاج ملف شبكي على النحو التالي: -T1 يحصل على قيم البيانات في جميع النقاط ، -T2 يحصل على البيانات
القيم عند نقاط مقيدة و NaN عند نقاط محرفة ؛ -T3 يحصل على 1 في
نقاط مقيدة و 0 عند نقاط محرفة [الافتراضي هو 1].

-الخامس[مستوى] (أكثر ...)
حدد مستوى الإسهاب [ج]. موصى به بشكل خاص هنا ، لأنه من المفيد رؤيته
كيف يتم تعديل الإحداثيات.

-Wماكسلون
تشير ماكسلون كحد أقصى لخط الطول لملف img الإدخال. إصدارات
منذ عام 1995 ماكسلون = 360.0 ، بينما كانت بعض الملفات السابقة ماكسلون = 390.0.
[الافتراضي هو 360.0].

-^ or م -
اطبع رسالة قصيرة حول صيغة الأمر ، ثم قم بالخروج (ملاحظة: في Windows
استخدم فقط -).

-+ or م +
اطبع رسالة (تعليمات) استخدام شاملة ، بما في ذلك شرح أي منها
خيار خاص بالوحدة النمطية (ولكن ليس خيارات GMT الشائعة) ، ثم يخرج.

-? or لا الحجج
اطبع رسالة استخدام (مساعدة) كاملة ، بما في ذلك شرح الخيارات ، بعد ذلك
المخارج.

--الإصدار
اطبع نسخة GMT واخرج.

--إظهار- datadir
طباعة المسار الكامل إلى دليل مشاركة GMT والخروج.

جغرافيا أمثلة

• -M يجب استبعاد الخيار إذا كنت تريد أن تكون شبكة الإخراج في منطقة جغرافية
إحداثيات. لاستخراج البيانات في المنطقة -R-40 / 40 / -70 / -30 من World_grav.img.7.2 و
إعادة طرح لإحداث إحداثيات جغرافية ، يمكنك المحاولة

img2grd world_grav.img.16.1 -Gmerc_grav.nc -R-40/40 / -70 / -30 -V

نظرًا لأن التباعد بين خطوط العرض في ملف img متساوي الأبعاد في وحدات Mercator ، فإن ملف
الشبكة الناتجة لن تتطابق مع المحدد -R بالضبط ، ولن يتم التباعد بين خطوط العرض
تساوي المسافات بين خطوط الطول. إذا كنت بحاجة إلى تطابق تام مع ملف -R ونفس التباعد
في خطوط الطول والعرض ، استخدم -E الخيار:

img2grd world_grav.img.16.1 -Gmerc_grav.nc -R-40/40 / -70 / -30 -E -V

ميركاتور أمثلة

نظرًا لأن ملفات img موجودة في إسقاط Mercator ، فلا يجب عليك استخراج شبكة جغرافية
إذا كانت خطتك هي إنشاء خريطة مركاتور. إذا فعلت ذلك فإنك ستنتهي من الإسقاط و
إعادة إسقاط الشبكة ، وفقدان تفاصيل الطول الموجي القصير. أفضل للاستخدام -M ورسم الشبكة
باستخدام إسقاط خطي بنفس مقياس إسقاط Mercator المطلوب (انظر GMT
مثال 29). لاستخراج البيانات في المنطقة -R-40 / 40 / -70 / -30 من World_grav.img.7.2، يركض

gmt img2grd -M world_grav.img.7.2 -Gmerc_grav.nc -R-40/40 / -70 / -30 -V

نلاحظ أن -V يخبرنا الخيار أنه تم تعديل النطاق إلى
-R-40 / 40 / -70.0004681551 / -29.9945810754. يمكننا أيضًا استخدام ملفات grdinfo للعثور على ملف الشبكة
يظهر header منطقته المراد أن تكون -R0/80/0 / 67.9666667 هذا هو نطاق x ، y الذي سنحصل عليه
من إسقاط مركاتور كروي باستخدام -R-40 / 40 / -70.0004681551 / -29.9945810754 و
-جم1. لذلك ، لأخذ ship.lonlatgrav واستخدامه لأخذ عينات من merc_grav.nc ، يمكننا القيام بذلك:

تعيين gmt المجال PROJ_ELLIPSOID

gmt mapproject -R-40/40 / -70.0004681551 / -29.9945810754 -Jm1i ship.lonlatgrav | \
gmt grdtrack -Gmerc_grav.nc | مخطط gmt \
-R-40/40 / -70.0004681551 / -29.9945810754 -Jm1i -I> ship.lonlatgravsat

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

لإنشاء شبكة lon ، lat من الشبكة أعلاه يمكننا استخدامها

gmt grdproject merc_grav.nc -R-40/40 / -70.0004681551 / -29.9945810754 -Jm1i -I -D2m -Ggrav.nc

في بعض الحالات ، لن يكون هذا سهلاً مثل -R في نظامي الإحداثيات قد لا تتم محاذاة
حسنًا. عندما يحدث هذا ، يمكننا أيضًا استخدام (في الواقع ، قد يكون من الأفضل دائمًا استخدام)

gmt grd2xyz merc_grav.nc | مخطط gmt \
-R-40/40 / -70.0004681551 / -29.994581075 -Jm1i -I | \
السطح GT-R-40/40 / -70 / 70 -I2m -Ggrav.nc

لعمل خريطة مركاتور للمنطقة أعلاه ، افترض قيمة gmt.conf الخاصة بنا لـ
PROJ_LENGTH_UNIT بوصة. ثم منذ أن تم عرض ملف merc_grav.nc أعلاه بامتداد -جم1i
يبلغ عرضها 80 بوصة. يمكننا إنشاء خريطة بعرض 8 بوصات باستخدام -جكس0.1i على أي برامج خرائط
المطبقة على هذه الشبكة (على سبيل المثال ، com.grdcontour, grdimage, com.grdview) ، ثم للتراكبات التي
العمل في lon ، lat (على سبيل المثال ، com.psxy, بسكوست) يمكننا استخدام التعديل أعلاه -R و -جم0.1 للحصول عليه
النظامين لمطابقة.

ومع ذلك ، يمكننا أن نكون أكثر ذكاءً من هذا. إدراك أن ملف img الإدخال يحتوي على وحدات بكسل 2.0
دقيقة (أو التحقق من nx و ny مع grdinfo merc_grav.nc) ندرك ذلك
استخدم merc_grav.nc الدقة الكاملة لملف img ولديه 2400 × 2039 بكسل ، و
بعرض 8 بوصات ، هذا 300 بكسل لكل بوصة. قررنا أننا لسنا بحاجة إلى هذا العدد الكبير ونحن
سيكون راضيًا عن 100 بكسل لكل بوصة ، لذلك نريد أن نقسم البيانات إلى 3 × 3
مربعات. (إذا أردنا رسمًا كفافيًا ، فسنختار على الأرجح متوسط ​​البيانات أكثر من ذلك بكثير
(على سبيل المثال ، 6 في 6) للحصول على ملامح سلسة.) نظرًا لأن عام 2039 غير قابل للقسمة على 3 ، فسنحصل على a
تعديل مختلف -R هذا الوقت:

gmt img2grd -M world_grav.img.7.2 -Gmerc_grav_2.nc -R-40/40 / -70 / -30 -N3 -V

هذه المرة نجد أن المنطقة المعدلة هي -R-40 / 40 / -70.023256525 / -29.9368261101 و
الناتج 800 × 601 بكسل ، حجم أفضل بالنسبة لنا. الآن يمكننا إنشاء مصطنعة
ملف الإضاءة لهذا الاستخدام grdgradient:

gmt grdgradient merc_grav_2.nc -Gillum.nc -A0 / 270 -Ne0.6

وإذا كان لدينا أيضًا ملف CPT يسمى "grav.cpt" ، فيمكننا إنشاء خريطة إغاثة مظللة بالألوان
مثله:

gmt grdimage merc_grav_2.nc -Iillum.nc -Cgrav.cpt -Jx0.1i -K> map.ps
gmt psbasemap -R-40/40 / -70.023256525 / -29.9368261101 -Jm0.1i -Ba10 -O >> map.ps

افترض أنك تريد الحصول على قيم البيانات المقيدة فقط من ملف img ، في خطوط الطول / العرض
إحداثيات. ثم اركض img2grd مع الالجائزة -T2 خيار ، استخدام grd2xyz لتفريغ القيم ، الأنابيب
من خلال grep -v NaN للقضاء على NaNs ، والأنابيب من خلال مخطط مع العكس
الإسقاط على النحو الوارد أعلاه.

استخدم img2grdgmt عبر الإنترنت باستخدام خدمات onworks.net