diff --git a/src/Bridge.cpp b/src/Bridge.cpp
index 84339dae..593da5f5 100644
--- a/src/Bridge.cpp
+++ b/src/Bridge.cpp
@@ -27,7 +27,6 @@
 */
 
 #include "Bridge.h"
-#include "io.h"
 
 float Bridge::_heightref;
 bool Bridge::_flatten;
diff --git a/src/Bridge.h b/src/Bridge.h
index 83bc348d..914b62a2 100644
--- a/src/Bridge.h
+++ b/src/Bridge.h
@@ -30,6 +30,7 @@
 #define Bridge_h
 
 #include "TopoFeature.h"
+#include "io.h"
 
 class Bridge: public Boundary3D {
 public:
diff --git a/src/Building.cpp b/src/Building.cpp
index c21c861d..1b95a477 100644
--- a/src/Building.cpp
+++ b/src/Building.cpp
@@ -27,8 +27,6 @@
 */
 
 #include "Building.h"
-#include "io.h"
-#include <algorithm>
 
 //-- static variable
 float Building::_heightref_top;
diff --git a/src/Building.h b/src/Building.h
index ac90ba56..a309006d 100644
--- a/src/Building.h
+++ b/src/Building.h
@@ -30,6 +30,8 @@
 #define Building_h
 
 #include "TopoFeature.h"
+#include "io.h"
+#include <algorithm>
 
 class Building: public Flat {
 public:
diff --git a/src/Forest.cpp b/src/Forest.cpp
index 3ccdd48f..33107f68 100644
--- a/src/Forest.cpp
+++ b/src/Forest.cpp
@@ -27,7 +27,6 @@
 */
 
 #include "Forest.h"
-#include "io.h"
 
 Forest::Forest(char *wkt, std::string layername, AttributeMap attributes, std::string pid, int simplification, double simplification_tinsimp, float innerbuffer)
   : TIN(wkt, layername, attributes, pid, simplification, simplification_tinsimp, innerbuffer) {}
diff --git a/src/Forest.h b/src/Forest.h
index bbd9ec9d..4fbe3681 100644
--- a/src/Forest.h
+++ b/src/Forest.h
@@ -30,6 +30,7 @@
 #define Forest_h
 
 #include "TopoFeature.h"
+#include "io.h"
 
 class Forest: public TIN {
 public:
diff --git a/src/Map3d.cpp b/src/Map3d.cpp
index 0db65ce5..489a2f0d 100644
--- a/src/Map3d.cpp
+++ b/src/Map3d.cpp
@@ -27,12 +27,6 @@
 */
 
 #include "Map3d.h"
-#include "io.h"
-#include "boost/locale.hpp"
-#include "boost/tokenizer.hpp"
-#include "boost/algorithm/string.hpp"
-#include "nlohmann-json/json.hpp"
-#include "codecvt"
 
 Map3d::Map3d() {
   OGRRegisterAll();
diff --git a/src/Map3d.h b/src/Map3d.h
index 63953c4a..add8beee 100644
--- a/src/Map3d.h
+++ b/src/Map3d.h
@@ -31,6 +31,7 @@
 
 #include "definitions.h"
 #include "geomtools.h"
+#include "io.h"
 #include "TopoFeature.h"
 #include "Building.h"
 #include "Terrain.h"
@@ -39,9 +40,7 @@
 #include "Road.h"
 #include "Separation.h"
 #include "Bridge.h"
-
-#include <set>
-#include <array>
+#include "boost/locale.hpp"
 
 typedef std::pair<Box2, TopoFeature*> PairIndexed;
 
diff --git a/src/Road.cpp b/src/Road.cpp
index 97adfc83..697c7fb2 100644
--- a/src/Road.cpp
+++ b/src/Road.cpp
@@ -27,7 +27,6 @@
 */
 
 #include "Road.h"
-#include "io.h"
 
 float Road::_heightref;
 bool  Road::_filter_outliers;
diff --git a/src/Road.h b/src/Road.h
index 509f1749..de661f3a 100644
--- a/src/Road.h
+++ b/src/Road.h
@@ -30,6 +30,7 @@
 #define Road_h
 
 #include "TopoFeature.h"
+#include "io.h"
 
 class Road: public Boundary3D {
 public:
diff --git a/src/Separation.cpp b/src/Separation.cpp
index d6e610b5..cd544f43 100644
--- a/src/Separation.cpp
+++ b/src/Separation.cpp
@@ -27,7 +27,6 @@
 */
 
 #include "Separation.h"
-#include "io.h"
 
 float Separation::_heightref;
 
diff --git a/src/Separation.h b/src/Separation.h
index 8207afee..67df13a8 100644
--- a/src/Separation.h
+++ b/src/Separation.h
@@ -30,6 +30,7 @@
 #define Separation_h
 
 #include "TopoFeature.h"
+#include "io.h"
 
 class Separation: public Boundary3D {
 public:
diff --git a/src/Terrain.cpp b/src/Terrain.cpp
index 522ee634..03a6ec1a 100644
--- a/src/Terrain.cpp
+++ b/src/Terrain.cpp
@@ -27,8 +27,6 @@
 */
 
 #include "Terrain.h"
-#include "io.h"
-#include <algorithm>
 
 Terrain::Terrain(char *wkt, std::string layername, AttributeMap attributes, std::string pid, int simplification, double simplification_tinsimp, float innerbuffer)
   : TIN(wkt, layername, attributes, pid, simplification, simplification_tinsimp, innerbuffer) {}
diff --git a/src/Terrain.h b/src/Terrain.h
index e995f739..4b69df27 100644
--- a/src/Terrain.h
+++ b/src/Terrain.h
@@ -30,6 +30,7 @@
 #define Terrain_h
 
 #include "TopoFeature.h"
+#include "io.h"
 
 class Terrain: public TIN {
 public:
diff --git a/src/TopoFeature.cpp b/src/TopoFeature.cpp
index ae12b623..5bad49c7 100644
--- a/src/TopoFeature.cpp
+++ b/src/TopoFeature.cpp
@@ -27,8 +27,6 @@
 */
 
 #include "TopoFeature.h"
-#include "io.h"
-#include "polyfit.hpp"
 
 TopoFeature::TopoFeature(char *wkt, std::string layername, AttributeMap attributes, std::string pid) {
   _id = pid;
diff --git a/src/TopoFeature.h b/src/TopoFeature.h
index c7237aa2..c75c005c 100644
--- a/src/TopoFeature.h
+++ b/src/TopoFeature.h
@@ -31,8 +31,9 @@
 
 #include "definitions.h"
 #include "geomtools.h"
+#include "io.h"
+#include "polyfit.hpp"
 #include "nlohmann-json/json.hpp"
-#include <random>
 
 class TopoFeature {
 public:
diff --git a/src/Water.cpp b/src/Water.cpp
index 65616a9c..3f8d91d5 100644
--- a/src/Water.cpp
+++ b/src/Water.cpp
@@ -27,7 +27,6 @@
 */
 
 #include "Water.h"
-#include "io.h"
 
 float Water::_heightref;
 
diff --git a/src/Water.h b/src/Water.h
index 0ef21099..2993ce89 100644
--- a/src/Water.h
+++ b/src/Water.h
@@ -30,6 +30,7 @@
 #define Water_h
 
 #include "TopoFeature.h"
+#include "io.h"
 
 class Water: public Flat {
 public:
diff --git a/src/definitions.h b/src/definitions.h
index a918c3b1..c3a56d57 100644
--- a/src/definitions.h
+++ b/src/definitions.h
@@ -2,11 +2,8 @@
 #define __3DFIER__Definitions__
 
 #include <iostream>
-#include <iomanip>
 #include <cmath>
 #include <vector>
-#include <set>
-#include <map>
 #include <unordered_map>
 
 #include <liblas/liblas.hpp>
@@ -15,17 +12,13 @@
 #include <boost/geometry.hpp>
 #include <boost/geometry/geometries/point_xy.hpp>
 #include <boost/geometry/geometries/box.hpp>
-#include <boost/geometry/geometries/linestring.hpp>
 #include <boost/geometry/geometries/polygon.hpp>
 #include <boost/geometry/index/rtree.hpp>
-#include "boost/filesystem.hpp"
-#include <boost/filesystem/operations.hpp>
 
 namespace bg = boost::geometry;
 namespace bgi = boost::geometry::index;
 typedef bg::model::d2::point_xy<double> Point2;
 typedef bg::model::segment<Point2> Segment2;
-typedef bg::model::linestring<Point2> Linestring2;
 typedef bg::model::polygon<Point2, true, false> Polygon2; //-- cw, first!=last
 typedef bg::model::ring<Point2, true, false> Ring2; //-- cw, first!=last
 typedef bg::model::box<Point2> Box2;
diff --git a/src/geomtools.cpp b/src/geomtools.cpp
index 2fd42306..c13cad6b 100644
--- a/src/geomtools.cpp
+++ b/src/geomtools.cpp
@@ -26,8 +26,8 @@
   Julianalaan 134, Delft 2628BL, the Netherlands
 */
 
-#include "io.h"
 #include "geomtools.h"
+#include "io.h"
 #include <CGAL/Exact_predicates_inexact_constructions_kernel.h>
 #include <CGAL/Constrained_Delaunay_triangulation_2.h>
 #include <CGAL/Projection_traits_xy_3.h>
@@ -37,8 +37,6 @@
 
 #include <vector>
 #include <unordered_set>
-#include <iterator>
-#include <memory>
 #include <boost/heap/fibonacci_heap.hpp>
 
 // fibonacci heap for greedy insertion code
diff --git a/src/main.cpp b/src/main.cpp
index 258aef98..70d8c3ee 100644
--- a/src/main.cpp
+++ b/src/main.cpp
@@ -35,7 +35,6 @@
 #include "boost/chrono.hpp"
 #include <boost/program_options.hpp>
 #include <boost/filesystem.hpp>
-#include <map>
 
 std::string VERSION = "1.0";