project logo

Εγχειρίδιο Γρήγορης Εκκίνησης libLAS

Εκτέλεση

Το εργαλείο τερματικού libLAS προμηθεύει μια πληθώρα από εύχρηστα λειτουργικά λογισμικά για το libLAS, αν και η θεμελιώδης libLAS βιβλιοθήκη δίνει πολλές δυνατότητες για αυτά. Παρακάτω παρουσιάζεται μια λίστα με τις βασικές λειτουργίες, που μπορεί να χρειαστείτε για εφαρμογή σε LAS δεδομένα, τις χρήσεις και τις προσεγγίσεις για την ολοκλήρωση των συγκεκριμένων θεμάτων.

Πραγματοποιώντας ορθή προβολή ενός αρχείου LAS

Όλα τα LAS δεδομένα είναι σε ένα σύστημα συντεταγμένων, ακόμα και αν το σύστημα αυτό δεν αναφέρεται στο αρχείο LAS. Για τα LAS δεδομένα που αφορούν την γή, οι περιγραφές του συστήματος συνταταγμένων συχνά αποτυπώνονται σε ένα σύστημα που περιγράφεται από το EPSG database. Άλλη πηγή πληροφοριών σχετικά με τα συστήματα συντεταγμένων μπορεί να βρεθεί στο http://spatialreference.org.

lasinfo --no-check srs.las

Note

Η επιλογή –no-check ορίζει στο lasinfο να εκτυπώσει μόνο τις επικεφαλίδες των πληροφοριών του αρχείου και να μην σαρώσει όλα τα σημεία. Για ένα αρχείο που περιέχει μόνο 10 σημεία δεν έχει μεγάλη σημασία, αλλά για ένα αρχείο με 50 και 500 χιλιάδες σημεία δεν αξίζει η αναμονή για την πλήρη σάρωση όλων των σημείων εάν αρκούν οι πληροφορίες από την επικεφαλίδα του αρχείου.

To ‘lasinfo’ ορίζει ότι το αρχείο srs.las είναι στο σύστημα συντεταγμένων UTM North Zone 17:

PROJCS["WGS 84 / UTM zone 17N",
    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["Transverse_Mercator"],
    PARAMETER["latitude_of_origin",0],
    PARAMETER["central_meridian",-81],
    PARAMETER["scale_factor",0.9996],
    PARAMETER["false_easting",500000],
    PARAMETER["false_northing",0],
    UNIT["metre",1,
        AUTHORITY["EPSG","9001"]],
    AUTHORITY["EPSG","32617"]]

Τώρα που είναι γνωστό το σύστημα συντεταγμένων που εισάγεται, μπορεί να παρθεί η απόφαση σε τι θα προβληθούν τα δεδομένα. Στο πρώτο παράδειγμα, θα χρησιμοποιηθεί το σύστημα Plate Carrée μη-συντεταγμένων EPSG:4326.

las2las srs.las --t_srs EPSG:4326

Η όλη διαδικασία ολοκληρώθηκε, αλλά μετά από μια γρήγορη επιθεώρηση των δεδομένων με το lasinfo output.las παρατηρήθηκε το πρόβλημα:

...
Scale Factor X Y Z:          0.01 0.01 0.01
Offset X Y Z:                -0.00 -0.00 -0.00
...
Min X, Y, Z:                -83.43, 39.01, 170.58,
Max X, Y, Z:                -83.43, 39.01, 170.76,

Το αρχείο srs.las είχε κλίμακα 0.01, ή δύο δεκαδικά ψηφία ακρίβεια για τις συντεταγμένες X, Y και Z. Για τα UTM δεδομένα, αυτό αρκεί, καθώς υποδηλώνει μια υπονοούμενη ακρίβεια του 1 cm. Για τα δεκαδικά δεδομένα του Plate Carrée συστήματος συντεταγμένων, το γεγονός αυτό μας κάνει να χάσουμε ακρίβεια. Πρέπει λοιπόν να οριστεί η τιμή της κλίμακας σε κάτι που έχει μεγαλύτερη ακρίβεια, στην περίπτωσή μας:

las2las --t_srs EPSG:4326 srs.las --scale 0.000001 0.000001 0.01

Μια άλλη γρήγορη επισκόπηση στο ‘lasinfo’ μας δίνει κάτι με το οποίο είμαστε πιο εξοικειωμένοι:

...
Scale Factor X Y Z:          0.000001 0.000001 0.01
Offset X Y Z:                -0.000000 -0.000000 -0.00
...
Min X, Y, Z:                -83.427598, 39.012599, 170.58
Max X, Y, Z:                -83.427548, 39.012618, 170.76

Αποτέλεσμα του LAS αρχείο στο κείμενο

las2txt input.las --parse xyzti