Cleanup, logger output setup

Former-commit-id: 1bdc478329

modules/odm_dem2mesh/CMakeLists.txt.user

@@ -1,6 +1,6 @@
 <?xml version="1.0" encoding="UTF-8"?>
 <!DOCTYPE QtCreatorProject>
-<!-- Written by QtCreator 4.7.0, 2018-10-16T15:52:36. -->
+<!-- Written by QtCreator 4.7.0, 2018-10-17T12:29:47. -->
 <qtcreator>
  <data>
   <variable>EnvironmentId</variable>
@@ -452,7 +452,7 @@
      <value type="int">13</value>
      <value type="int">14</value>
     </valuelist>
-    <value type="QString" key="CMakeProjectManager.CMakeRunConfiguration.Arguments">-inputFile /data/drone/cmparks/mesh_dsm_12.tif -outputFile /data/drone/cmparks/mesh_dsm.ply -verbose</value>
+    <value type="QString" key="CMakeProjectManager.CMakeRunConfiguration.Arguments">-inputFile /data/drone/cmparks/mesh_dsm.tif -outputFile /data/drone/cmparks/mesh_dsm.ply -verbose</value>
     <value type="QString" key="CMakeProjectManager.CMakeRunConfiguration.UserWorkingDirectory"></value>
     <value type="QString" key="CMakeProjectManager.CMakeRunConfiguration.UserWorkingDirectory.default">/data/OpenDroneMap/modules/build-odm_dem2mesh-Desktop-Default</value>
     <value type="int" key="PE.EnvironmentAspect.Base">2</value>

modules/odm_dem2mesh/src/Logger.h

@@ -26,7 +26,7 @@ struct Logger{
             {
                     va_list args;
                     va_start( args , format );
-                    vfprintf( stderr, format , args );
+                    vfprintf( stdout, format , args );
                     va_end( args );
             }
     }

modules/odm_dem2mesh/src/main.cpp

@@ -273,8 +273,8 @@ void transform(const BoundingBox &extent){
     double ext_height = extent.max.y - extent.min.y;
 
     for(Simplify::Vertex &v : Simplify::vertices){
-        v.p.x = extent.min.x + (static_cast<float>(v.p.x) / static_cast<float>(arr_width)) * ext_width;
-        v.p.y = extent.max.y - (static_cast<float>(v.p.y) / static_cast<float>(arr_height)) * ext_height;;
+        v.p.x = extent.min.x + (static_cast<double>(v.p.x) / static_cast<double>(arr_width)) * ext_width;
+        v.p.y = extent.max.y - (static_cast<double>(v.p.y) / static_cast<double>(arr_height)) * ext_height;;
     }
 }
 
@@ -284,6 +284,7 @@ int main(int argc, char **argv) {
     if ( !MaxVertexCount.set ) MaxVertexCount.value = 100000;
 
     logWriter.verbose = Verbose.set;
+    logWriter.outputFile = "odm_dem2mesh.txt";
     logArgs(params, logWriter);
 
     GDALDataset  *dataset;
@@ -307,7 +308,7 @@ int main(int argc, char **argv) {
 
         GDALRasterBand *band = dataset->GetRasterBand(1);
 
-        int qtreeLevels = 2;
+        int qtreeLevels = 1;
         subdivisions = (int)pow(2, qtreeLevels);
         int numBlocks = subdivisions * subdivisions;
         blockSizeX = arr_width / subdivisions;
@@ -427,6 +428,10 @@ int main(int argc, char **argv) {
                     ss << OutputFile.value << "." << blockX << "-" << blockY << ".bin";
                     logWriter("Reading %s\n", ss.str().c_str());
                     readBin(ss.str(), blockX, blockY);
+
+                    if (std::remove(ss.str().c_str()) != 0){
+                        logWriter("Error while deleting intermediate file: %s\n", ss.str().c_str());
+                    }
                 }
             }
 

opendm/mesh.py

@@ -38,8 +38,8 @@ def create_25dmesh(inPointCloud, outMesh, dsm_radius=0.07, dsm_resolution=0.05,
     mesh = dem_to_mesh(os.path.join(tmp_directory, 'mesh_dsm.tif'), outMesh, maxVertexCount, verbose)
 
     # Cleanup tmp
-    # if os.path.exists(tmp_directory):
-    #     shutil.rmtree(tmp_directory)
+    if os.path.exists(tmp_directory):
+        shutil.rmtree(tmp_directory)
 
     return mesh
 

scripts/odm_meshing.py

@@ -88,8 +88,8 @@ class ODMeshingCell(ecto.Cell):
               if args.fast_orthophoto:
                   dsm_radius *= 2
               
-              # A good DSM size is 1/8 of the target orthophoto resolution
-              dsm_resolution = ortho_resolution * 8
+              # A good DSM size is 1/4 of the target orthophoto resolution
+              dsm_resolution = ortho_resolution * 4
 
               mesh.create_25dmesh(infile, tree.odm_25dmesh,
                     dsm_radius=dsm_radius,