<div dir="ltr">
<div>Dear all, </div><div>I'm working with QGIS and PostGIS. As input, I
 have 25 polygonal layers covering a large area (multicities area). One 
of these data is a very large dataset (1 million objects). The other 24 
are much smaller (a maximum of a hundred objects).</div>For information,
 I should point out that some of these polygonal datasets are in 
"multi-part features" mode and others in "single-part features" mode. I 
imagine this may ultimately have a slight impact on the method/result. 
These 25 polygonal .shp files have highly variable, 
non-homogeneous/non-harmonized data structures. Each layer has a 
"data_id" field that allows  to define/link/reference, for each feature,
 its membership in the layer. For example, all values in the "data_id"
 field for the first layer have a value of '1'. For the second layer, 
the field values are '2', etc.<br><br>My goal would be to be able to apply/adapt the existing QGIS geoprocessing tool called "Multiple Union":<br><a href="https://docs.qgis.org/3.40/en/docs/user_manual/processing_algs/qgis/vectoroverlay.html#union-multiple">https://docs.qgis.org/3.40/en/docs/user_manual/processing_algs/qgis/vectoroverlay.html#union-multiple</a><br><br>Below a screenshot from the QGIS documentation :<br><br><div><img src="cid:ii_mcswuu7s0" alt="image.png" width="385" height="185" style="margin-right: 0px;"></div><div><br></div>My goal would be to have an output file:<br><br><ul><li> Which
 would be the result of the union/overlay of the 25 input data. To use 
the terms of the QGIS documentation, the processing should check for 
overlaps between features within the 25 layers and create separate 
features for the overlapping and non-overlapping parts. This "multiple 
union" geoprocessing seems interesting for my goal where there is no 
overlap (a, NULL; b, NULL; c, NULL).</li></ul><ul><li>For areas where 
there is an overlap, the QGIS union geoprocessing creates as many 
identical overlapping features as there are features participating in 
this overlap. This doesn't bother me. But since, ultimately, I'd like a 
field in the result/output file to allow, for each feature, to retrieve 
the list of input layers that participate/contribute to this result 
feature (in order to retrieve the origin/source of the data). I was 
wondering/thinking it might be better if only one feature was created 
per overlapping area?</li></ul><ul><li> I'd like a field in the result 
file to allow, for each feature, to retrieve the list of input layers 
that participate/contribute to this result feature. In order to retrieve
 the origin/source of the data.</li></ul><ul><li>Ideally, a field that 
allows you to retrieve the number (COUNT) of layers that contribute to 
this feature (at least 1 layer, at most 25 layers).</li></ul><ul><li>Regarding
 the non-geometric attributes/fields, I would like to be able to specify
 the selection/define the list of fields I ultimately want to keep. I 
don't want to keep all of the fields, but rather just some of the fields
 for each of the 25 input layers.</li></ul><br>I imagine it's 
recommended to do this processing in PostGIS rather than QGIS? I can, if
 necessary, import my 25 SHP files into a PostGIS database. I also 
imagine it's important to keep in mind that the "multi-part features" / 
"single-part pieces/features" mode of the input layers can affect the 
result. If I'm using a PostGIS query, I was thinking it might be helpful
 to force all features to be in single-part mode (using the PostGIS 
'st_dump' function?).<br><br><div>In advance, Thanks so much for your help, guidance.</div><br>

<br></div>