My goal is to link all the attributes of an input layer (that has a UTM CRS projection matrix and point coordinates but no geometries) to the geometries of another layer based on the spatial analysis. In particular, I want to join them if they intersect or if the geometry contains the dot equal to the coordinates in the input layer. Using the plugin MMQGIS, I am able to merge the two dataset. The problem is following: often multiple observations lies on the same geometry. MMQGIS only returns a random attribute of the input layer but reports correctly the number of points that intersect the geometry of the second one. I would like to have instead a new layer in which all the input attributes are reported. This means that the geometries' id will be duplicated in those cases with multiple inputs on the same geometry. Is there a way to easily perform this task? Both PyQGIS or icons solutions are useful.