Spatial Indexing #3
There are no files selected for viewing
This file contains hidden or bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
This file contains hidden or bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
This file contains hidden or bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
This file contains hidden or bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
This file contains hidden or bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| Original file line number | Diff line number | Diff line change |
|---|---|---|
| @@ -0,0 +1,91 @@ | ||
| /* | ||
| * SPDX-FileCopyrightText: (c) 2025-2026 Trailblaze Software, all rights reserved | ||
| * SPDX-License-Identifier: MIT | ||
| */ | ||
|
|
||
| #include <cmath> | ||
| #include <fstream> | ||
| #include <iomanip> | ||
| #include <iostream> | ||
| #include <string> | ||
|
|
||
| #include "las_reader.hpp" | ||
| #include "spatial_index.hpp" | ||
| #include "utilities/assert.hpp" | ||
| #include "utilities/printing.hpp" | ||
| #include "vlr.hpp" | ||
|
|
||
| using namespace laspp; | ||
|
|
||
| int main(int argc, char* argv[]) { | ||
|
Check warning on line 20 in apps/inspect_vlrs.cpp
|
||
| if (argc != 2) { | ||
| std::cerr << "Usage: " << argv[0] << " <laz_file>" << std::endl; | ||
| return 1; | ||
| } | ||
|
|
||
| std::string filename = argv[1]; | ||
| std::filesystem::path file_path(filename); | ||
| LASReader reader(file_path); | ||
|
|
||
| std::cout << "=== VLRs ===" << std::endl; | ||
| const auto& vlrs = reader.vlr_headers(); | ||
| for (size_t i = 0; i < vlrs.size(); ++i) { | ||
| const auto& vlr = vlrs[i]; | ||
| std::cout << "VLR " << i << ":" << std::endl; | ||
| std::cout << indented(vlr, " "); | ||
| std::cout << " Is LAZ VLR: " << (vlr.is_laz_vlr() ? "yes" : "no") << std::endl; | ||
| std::cout << " Is LAStools spatial index VLR: " | ||
| << (vlr.is_lastools_spatial_index_vlr() ? "yes" : "no") << std::endl; | ||
| std::cout << std::endl; | ||
| } | ||
|
|
||
| std::cout << "=== EVLRs ===" << std::endl; | ||
| const auto& evlrs = reader.evlr_headers(); | ||
| for (size_t i = 0; i < evlrs.size(); ++i) { | ||
| const auto& evlr = evlrs[i]; | ||
| std::cout << "EVLR " << i << ":" << std::endl; | ||
| std::cout << indented(evlr, " "); | ||
| std::cout << " Is LAZ VLR: " << (evlr.is_laz_vlr() ? "yes" : "no") << std::endl; | ||
| std::cout << " Is LAStools spatial index EVLR: " | ||
| << (evlr.is_lastools_spatial_index_evlr() ? "yes" : "no") << std::endl; | ||
|
|
||
| // If it's a LAStools spatial index, show more details | ||
| if (evlr.is_lastools_spatial_index_evlr()) { | ||
| LASPP_ASSERT(reader.has_lastools_spatial_index(), | ||
| "Failed to load LAStools spatial index from EVLR."); | ||
| std::cout << "Spatial index: " << std::endl; | ||
| const auto& index = reader.lastools_spatial_index(); | ||
| const auto& quadtree = index.quadtree_header(); | ||
|
|
||
| std::cout << " Spatial index header:" << std::endl; | ||
| std::cout << indented(quadtree, " "); | ||
| std::cout << " Number of cells: " << index.num_cells() << std::endl; | ||
|
|
||
| // Show cell details | ||
| size_t total_intervals = 0; | ||
| for (const auto& [cell_index, cell] : index.cells()) { | ||
| total_intervals += cell.intervals.size(); | ||
| } | ||
| std::cout << " Total intervals: " << total_intervals << std::endl; | ||
| std::cout << " Average intervals per cell: " << std::fixed << std::setprecision(1) | ||
| << (index.num_cells() > 0 ? static_cast<double>(total_intervals) / | ||
| static_cast<double>(index.num_cells()) | ||
| : 0.0) | ||
| << std::endl; | ||
|
|
||
| // Show first few cells as examples | ||
| std::cout << indented(limited_map(index.cells(), 3), " ") << std::endl; | ||
| } | ||
| std::cout << std::endl; | ||
| } | ||
|
|
||
| if (reader.has_lastools_spatial_index()) { | ||
| const auto& spatial_index = reader.lastools_spatial_index(); | ||
| std::cout << "Loaded LAStools spatial index with " << spatial_index.num_cells() << " cells." | ||
| << std::endl; | ||
| std::cout << "Spatial index details:" << std::endl; | ||
| std::cout << indented(spatial_index, " ") << std::endl; | ||
| } | ||
|
|
||
| return 0; | ||
| } | ||
Oops, something went wrong.
Oops, something went wrong.
Add this suggestion to a batch that can be applied as a single commit.
This suggestion is invalid because no changes were made to the code.
Suggestions cannot be applied while the pull request is closed.
Suggestions cannot be applied while viewing a subset of changes.
Only one suggestion per line can be applied in a batch.
Add this suggestion to a batch that can be applied as a single commit.
Applying suggestions on deleted lines is not supported.
You must change the existing code in this line in order to create a valid suggestion.
Outdated suggestions cannot be applied.
This suggestion has been applied or marked resolved.
Suggestions cannot be applied from pending reviews.
Suggestions cannot be applied on multi-line comments.
Suggestions cannot be applied while the pull request is queued to merge.
Suggestion cannot be applied right now. Please check back later.
Uh oh!
There was an error while loading. Please reload this page.