יותר

נתוני נקודת XYZ לא סדירים לרסטר (TIF) - חלופה של LASTools?

נתוני נקודת XYZ לא סדירים לרסטר (TIF) - חלופה של LASTools?


אני מנסה להמיר את נתוני הנקודות הבלתי סדירות שלי (xyz) ל- DEM ולאחר מכן לרסטר (geotiff). בגלל המרווח הלא סדיר אני צריך טריאנגולציה קודם כל, נכון (אז, אין רווחים ברסטר הסופי)?

אז הדרך היא: XYZ -> Triangulate -> Interpolate To Raster. יש לי הרבה נקודות (מינימום 150.000), גודל הרסטר שלי צריך להיות 0.25 על 0.25 מ 'או קטן יותר. אני צריך דרך לעשות את זה עם קוד פתוח / כלים בחינם!

מה שניסיתי עד כה:

  • LAST כלים (las2dem) - עובד מצוין, מהיר במיוחד < - זה בדיוק מה שאני צריך, אבל זה לא בחינם.

  • היתוך (TINSurfaceCreateלאחר מכןDTM2ASCII) - עובד טוב, אבל מאט במיוחד מכיוון שגודל הרסטר מתכווץ (1 x 1 מ 'ו- 0.5 x 0.5 מ' זה בסדר, הכל למטה ממש איטי).

  • SAGA כלים /gdal_grid(באמצעות QGIS) -> גם איטי במיוחד.

יש עוד הצעות שכדאי לנסות?


אני חושב שהעיקר עבור לידר הוא האם אתה צריך לסנן את ההחזרים הלא קרקעיים - אני מניח שכבר עשית זאת.

אתה לא צריך טריאנגולציה בהכרח, יש גם kriging, IDW, צלחת דקה ופלטות אחרות. עיין בחבילות gstat ו- raster ב- R לקבלת מידע נוסף.

ובכל זאת, חבילת RTriangle ב- R מציעה טריאנגולציה מהירה במיוחד (עם אילוצים) ולכן הייתי מתחיל שם. חבילות אחרות מציעות טריאנגולציה של Delaunay והן יהיו טובות לא פחות במקרה זה.

הכנתי כאן מסמך סימון גס מאוד: https://github.com/mdsumner/xyztrisurf/blob/master/xyztrisurf.md

ניתן יהיה להגדיל זאת לנתונים שלך, אך ייתכן שיהיה עליך טיפול זהיר למדי - לא עשיתי דבר כדי לראות כיצד הדבר יתבצע כאן.

FWIW, זה "אינטרפולט לרשת ממשולשים" הוא בעצם מה שהגרידטה של ​​מטלאב עושה כברירת מחדל (אפייק), ומה שעושה מניפולד GIS בשיטת "הטריאנגולציה" ליצירת משטח מנקודות.


PDAL יכול לעשות זאת באמצעות points2grid. עיין ב http://www.pdal.io/stages/writers.p2g.html לפרטים על אופן פעולתו, ועיין ב https://github.com/Applied-GeoSolutions/lidar2dems לקבלת פתרון צינורות עיבוד אם עליך לעשות זאת הרבה מזה.

שים לב שעדיין אין קורא ASCII ל- PDAL, והוא מניח שאתה משתמש בפורמט LiDAR כלשהו.


אז הדרך לעקיפת הבעיה הייתה כדלקמן: אני משתמשהיתוךלמשולש ולחישוב רסטר. זה שם א*.dtmקובץ החוצה (בינארי) -> ההמרה עםהיתוךכליםDTM2ASCIIהוא צוואר הבקבוק.

אז שמתי את ידי על הקורא הבינארי של C#, קרא את כל*.dtmקובץ וכתבתי פלט tif משלי עםlibtiff.net(http://bitmiracle.com/libtiff/).

זה לא מהיר כמו LAStools, אבל לצרכים שלי זה בסדר ועובד די טוב :).


צפו בסרטון: QGIS and LiDAR LAStools