أكثر

النقطية المحاصيل من مشكلة المتجهات

النقطية المحاصيل من مشكلة المتجهات


لذا لدي ملف شكل لناقل لخطوط مدينة فيلادلفيا وصورة نقطية لمقاطعة فيلادلفيا حصلت عليها من أشخاص رائعين من وزارة الزراعة الأمريكية: NRCS. لقد قمت بتحويل البيانات النقطية من ملف sid إلى TIFF (باستخدام gdal_translate) ثم قمت بإعادة إسقاط صيغة TIFF النقطية المحولة حديثًا إلى 3857 (باستخدام gdalwarp).

يحتوي هذا النقطي الجديد الآن على نفس الإسقاط الأصلي مثل المتجه لحدود المدينة: EPSG: 3857. لذلك أقوم بتحميل كل من ملف الشكل لحدود المدينة وهذا tif النقطي الجديد الذي أعددته في QGIS وأطلق أداة المقص المفيدة (Raster-> Extraction-> Clipper) وأقوم بعملي المتمثل في تعيين البيانات النقطية كمصدر واختيار المتجه كطبقة القناع واضغط على "موافق". #thatwassimple. يبدأ الكمبيوتر في العمل بشكل متواصل أثناء القيام بأشياءه ويتم تحفيزي بشأن وصول البيانات النقطية الجديدة.

دينغ! انتهى الكمبيوتر ... فقط جنوب فيلي وجزء صغير من وسط فيلي - يبدو أن فيلي الغربية والشمالية والشرقية قد تهربوا من الحفلة وكلها سوداء.

حسنًا ... هذا غريب ، ربما يمكنني استخدام سطر أوامر GDAL لفرض ذلك ، لأنك تعلم ، إذا تم ذلك في سطر الأوامر ، فأنا متسلل حقيقي وليس مثل QGIS يفعل ذلك بالفعل في الخلفية ، لكنني استطراداً.

لذلك أنا كمد تلك mofos (باستخدام gdalwarp مرة أخرى!) 0… 10… 100! دينغ! حسنًا ، حسنًا ، حسنًا ، لأنني ماكونهي بنفسي انتقلت إلى الدليل الذي أخرجت ملف TIFF المقطوع. افتح هذا المصاص ، وأنا مع نفس الناتج بالضبط من ساوثسايد وسنترال بقيت فقط لتمثيلها.

لقد أرفقت صورتين ، الأولى تظهر مضلعًا وردي اللون أريد قصه من الخطوط النقطية ، والثانية تُظهر المخرجات في المرتين.

الصورة 1:

الصورة 2:

ملفات المصدر

الشكل المستخدم (في الرمز البريدي)

تم استخدام TIFF (في ملف zip)

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

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


يبدو أن هناك فرقًا بسيطًا ولكنه مهم بين الإسقاطات التي تجعل gdalwarp تتعثر.

جريgdalsrsinfo philly_3857.tif> city_limits.prjوأدى إزالة الأسطر الأولى منه إلى الحيلة بالنسبة لي.

للمقارنة:

philly_3857.tif:

PROJCS ["WGS 84 / Pseudo-Mercator"، GEOGCS ["WGS 84"، DATUM ["WGS_1984"، SPHEROID ["WGS 84"، 6378137،298.257223563، AUTHORITY ["EPSG"، "7030"]]، AUTHORITY [" EPSG "،" 6326 "]]، PRIMEM [" Greenwich "، 0]، UNIT [" degree "، 0.0174532925199433]، AUTHORITY [" EPSG "،" 4326 "]]، PROJECTION [" Mercator_1SP "]، PARAMETER [" central_meridian "، 0]، PARAMETER [" scale_factor "، 1]، PARAMETER [" false_easting "، 0]، PARAMETER [" false_northing "، 0]، UNIT [" meter "، 1، AUTHORITY [" EPSG "،" 9001 "] ]، EXTENSION ["PROJ4"، "+ proj = merc + a = 6378137 + b = 6378137 + lat_ts = 0.0 + lon_0 = 0.0 + x_0 = 0.0 + y_0 = 0 + k = 1.0 + الوحدات = m + nadgrids = @ فارغ + wktext + no_defs "]، AUTHORITY [" EPSG "،" 3857 "]]

city_limits.prj (الأصل):

PROJCS ["WGS_84_Pseudo_Mercator"، GEOGCS ["GCS_WGS_1984"، DATUM ["D_WGS_1984"، SPHEROID ["WGS_1984"، 6378137،298.257223563]]، PRIMEM ["Greenwich" ["0]"، 0.02599 ["Mercator"]، PARAMETER ["central_meridian"، 0]، PARAMETER ["false_easting"، 0]، PARAMETER ["false_northing"، 0]، UNIT ["Meter"، 1]، PARAMETER ["standard_parallel_1"، 0.0] ]

بالنسبة للمستخدم العادي ، يبدو أنهما متماثلان. ولكن إذا حذفت ملفات .q ** من ملف zip ، فسيتم تحميل ملف الشكل في QGIS بتنسيقEPSG: 54004، وسوف ليس تتماشى مع خطوط المسح. لا يبحث gdalwarp في ملفات .q ** ، وبالتالي ينتج عنه ملف مهمل.


شاهد الفيديو: Vector Space شرح كامل بالتفصيل مع امثله كثير وطرق مختلفه للحل