From 51e217bf052c364f1f96edf3ee4330e3b4100db7 Mon Sep 17 00:00:00 2001
From: Emmanuel Christophe <emmanuel.christophe@orfeo-toolbox.org>
Date: Fri, 5 Dec 2008 14:08:31 +0800
Subject: [PATCH] ENH: Adding lidar to image example
---
Examples/IO/LidarToImageExample.cxx | 25 +++++++++++++++++--------
1 file changed, 17 insertions(+), 8 deletions(-)
diff --git a/Examples/IO/LidarToImageExample.cxx b/Examples/IO/LidarToImageExample.cxx
index 441fbc2eb8..0ddc5ca492 100644
--- a/Examples/IO/LidarToImageExample.cxx
+++ b/Examples/IO/LidarToImageExample.cxx
@@ -31,13 +31,12 @@ int main( int argc, char* argv[] )
{
- if(argc!=3)
+ if(argc!=6)
{
std::cout << argv[0] <<" <input_lidar_filename> <output_image_filename>"
- << " <output_resolution> <upper_left_corner_latitude>"
- << " <size_x> <sizee_y> <number_of_stream_divisions>"
+ << " <output_resolution> <spline_order>"
+ << " <number_of_level>"
<< std::endl;
-
return EXIT_FAILURE;
}
@@ -65,21 +64,31 @@ int main( int argc, char* argv[] )
start[0] = 0;
start[1] = 0;
+// std::cout << std::setprecision(15);
+// std::cout << reader->GetMinX() << std::endl;
+// std::cout << reader->GetMaxX() << std::endl;
+// std::cout << reader->GetMinY() << std::endl;
+// std::cout << reader->GetMaxY() << std::endl;
+
ImageType::SizeType size;
size[0] = static_cast<long int >(
- ceil((reader->GetMaxX()-reader->GetMinX()+1)/resolution));
+ ceil((reader->GetMaxX()-reader->GetMinX()+1)/resolution))+1;
size[1] = static_cast<long int >(
- ceil((reader->GetMaxY()-reader->GetMinY()+1)/resolution));
+ ceil((reader->GetMaxY()-reader->GetMinY()+1)/resolution))+1;
ImageType::PointType origin;
- origin[0] = 0;//minX;
- origin[1] = 0;//maxY;
+ origin[0] = reader->GetMinX();//minX;
+ origin[1] = reader->GetMaxY();//maxY;
ImageType::SpacingType spacing;
spacing[0] = resolution;
spacing[1] = -resolution;
+// std::cout << "Size: " << size << std::endl;
+// std::cout << "Origin: " << origin << std::endl;
+// std::cout << "Spacing: " << spacing << std::endl;
+
typedef itk::BSplineScatteredDataPointSetToImageFilter
<PointSetType, VectorImageType> FilterType;
--
GitLab