Merges compatible points into multipoints, closes systemed#719

src/tile_worker.cpp

@@ -275,15 +275,38 @@ void ProcessObjects(
 		OutputObjectID oo = *jt;
 		if (zoom < oo.oo.minZoom) { continue; }
 
-		if (oo.oo.geomType == POINT_) {
+                if (oo.oo.geomType == POINT_) {
+			// The very first point; below we check if there are more compatible points
+			// so that we can write a multipoint instead of many point features
+
+			std::vector<std::pair<int, int>> multipoint;
+
+			LatpLon pos = source->buildNodeGeometry(jt->oo.objectID, bbox);
+			pair<int,int> xy = bbox.scaleLatpLon(pos.latp/10000000.0, pos.lon/10000000.0);
+			multipoint.push_back(xy);
+
+			while (jt<(ooSameLayerEnd-1) && oo.oo.compatible((jt+1)->oo)) {
+				jt++;
+				LatpLon pos = source->buildNodeGeometry(jt->oo.objectID, bbox);
+				pair<int,int> xy = bbox.scaleLatpLon(pos.latp/10000000.0, pos.lon/10000000.0);
+				multipoint.push_back(xy);
+			}
+
 			vtzero::point_feature_builder fbuilder{vtLayer};
 			if (sharedData.config.includeID && oo.id) fbuilder.set_id(oo.id);
 
-			LatpLon pos = source->buildNodeGeometry(oo.oo.objectID, bbox);
-			pair<int,int> xy = bbox.scaleLatpLon(pos.latp/10000000.0, pos.lon/10000000.0);
-			fbuilder.add_point(xy.first, xy.second);
+			fbuilder.add_points(multipoint.size());
+
+			if (verbose && multipoint.size() > 1)
+				std::cout << "Merging " << multipoint.size() << " points into a multipoint" << std::endl;
+
+			for (const auto point : multipoint)
+				fbuilder.set_point(point.first, point.second);
+
 			oo.oo.writeAttributes(attributeStore, fbuilder, zoom);
 			fbuilder.commit();
+
+			oo = *jt;
 		} else {
 			Geometry g;
 			try {