LasDetails.h

Go to the documentation of this file.

// ----------------------------------------------------------------------------
// -                        CloudViewer: www.cloudViewer.org                  -
// ----------------------------------------------------------------------------
// Copyright (c) 2018-2024 www.cloudViewer.org
// SPDX-License-Identifier: MIT
// ----------------------------------------------------------------------------
 
#pragma once
 
// ##########################################################################
// #                                                                        #
// #                CloudViewer PLUGIN: LAS-IO Plugin                      #
// #                                                                        #
// #  This program is free software; you can redistribute it and/or modify  #
// #  it under the terms of the GNU General Public License as published by  #
// #  the Free Software Foundation; version 2 of the License.               #
// #                                                                        #
// #  This program is distributed in the hope that it will be useful,       #
// #  but WITHOUT ANY WARRANTY; without even the implied warranty of        #
// #  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the         #
// #  GNU General Public License for more details.                          #
// #                                                                        #
// #                   COPYRIGHT: Thomas Montaigu                           #
// #                                                                        #
// ##########################################################################
 
// CC
#include <CVTypes.h>
 
// Qt
#include <QtGlobal>
 
// System
#include <array>
#include <cmath>
#include <limits>
#include <string>
#include <vector>
 
class ccPointCloud;
class ccScalarField;
 
class QDataStream;
 
struct laszip_header;
struct laszip_vlr;
typedef laszip_vlr laszip_vlr_struct;
 
constexpr size_t LAS_VLR_HEADER_SIZE = 54;
constexpr double SCAN_ANGLE_SCALE    = 0.06;
 
namespace LasNames
{
    constexpr const char* Intensity         = "Intensity";
    constexpr const char* ReturnNumber      = "Return Number";
    constexpr const char* NumberOfReturns   = "Number Of Returns";
    constexpr const char* ScanDirectionFlag = "Scan Direction Flag";
    constexpr const char* EdgeOfFlightLine  = "EdgeOfFlightLine";
    constexpr const char* Classification    = "Classification";
    constexpr const char* SyntheticFlag     = "Synthetic Flag";
    constexpr const char* KeypointFlag      = "Keypoint Flag";
    constexpr const char* WithheldFlag      = "Withheld Flag";
    constexpr const char* ScanAngleRank     = "Scan Angle Rank";
    constexpr const char* UserData          = "User Data";
    constexpr const char* PointSourceId     = "Point Source ID";
    constexpr const char* GpsTime           = "Gps Time";
 
    // 1.4 point format 6 stuff
    constexpr const char* ScanAngle      = "Scan Angle";
    constexpr const char* ScannerChannel = "Scanner Channel";
    constexpr const char* OverlapFlag    = "Overlap Flag";
    constexpr const char* NearInfrared   = "Near Infrared";
} // namespace LasNames
 
namespace LasDetails
{
    // The position of the overlap flag in the classification flags
    // (valid for fmt >= 6)
    constexpr unsigned OVERLAP_FLAG_BIT_POS  = 3;
    constexpr unsigned OVERLAP_FLAG_BIT_MASK = 1 << OVERLAP_FLAG_BIT_POS;
 
    struct EvlrHeader
    {
        static constexpr size_t SIZE             = 60;
        static constexpr size_t USER_ID_SIZE     = 16;
        static constexpr size_t DESCRIPTION_SIZE = 32;
 
        char     userID[USER_ID_SIZE];
        uint16_t recordID{0};
        uint64_t recordLength{0};
        char     description[DESCRIPTION_SIZE];
 
        EvlrHeader() = default;
 
        static EvlrHeader Waveform();
 
        bool isWaveFormDataPackets() const;
 
        friend QDataStream& operator>>(QDataStream& stream, EvlrHeader& hdr);
        friend QDataStream& operator<<(QDataStream& stream, const EvlrHeader& hdr);
    };
 
    uint16_t PointFormatSize(unsigned pointFormat);
 
    uint16_t HeaderSize(unsigned versionMinor);
 
    inline bool HasGpsTime(unsigned pointFormatId)
    {
        return pointFormatId == 1
               || pointFormatId >= 3;
    }
 
    inline bool HasRGB(unsigned pointFormatId)
    {
        return pointFormatId == 2
               || pointFormatId == 3
               || pointFormatId == 5
               || pointFormatId == 7
               || pointFormatId == 8
               || pointFormatId == 10;
    }
 
    inline bool HasWaveform(unsigned pointFormatId)
    {
        return pointFormatId == 4
               || pointFormatId == 5
               || pointFormatId >= 9;
    }
 
    inline bool HasNearInfrared(unsigned pointFormatId)
    {
        return pointFormatId == 8
               || pointFormatId == 10;
    }
 
    unsigned SizeOfVlrs(const laszip_vlr_struct* vlrs, unsigned numVlrs);
 
    bool IsLaszipVlr(const laszip_vlr_struct&);
 
    bool IsExtraBytesVlr(const laszip_vlr_struct&);
 
    const std::vector<unsigned>* PointFormatsAvailableForVersion(QString version);
 
    const std::array<const char*, 3>& AvailableVersions();
 
    struct LasVersion
    {
        int pointFormat  = 3;
        int minorVersion = 2;
    };
 
    LasVersion SelectBestVersion(const ccPointCloud& cloud);
 
    void CloneVlrInto(const laszip_vlr_struct& src, laszip_vlr_struct& dst);
 
} // namespace LasDetails