Hi list,<br><br>is it correct that multipolygons are not implemented for the Intersects spatial filter? I can not get it to work, also not with 5.0, but I assume GEOS is able to perform this operation or not?<br><br>in mapogcfilter.c
 I see:<br><br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; psGMLElement = CPLGetXMLNode(psXMLNode, &quot;Polygon&quot;);<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; if (psGMLElement)<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; bPolygon = 1;<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; else if ((psGMLElement= CPLGetXMLNode(psXMLNode, &quot;LineString&quot;)))
<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; {<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; if (psGMLElement)<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; bLine = 1;<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; }<br>&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; else<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; {<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; psGMLElement = CPLGetXMLNode(psXMLNode, &quot;Point&quot;);
<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; if (!psGMLElement)<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; psGMLElement =&nbsp; CPLGetXMLNode(psXMLNode, &quot;PointType&quot;);<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; if (psGMLElement)<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; bPoint =1;<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; } 
<br><br>which leads to me to believe multipolygons are not implemented.<br><br>This is the filter I am trying:<br><br>&lt;Filter xmlns=&quot;<a href="http://www.opengis.net/ogc">http://www.opengis.net/ogc</a>&quot; xmlns:xsi=&quot;
<a href="http://www.w3.org/2001/XMLSchema-instance">http://www.w3.org/2001/XMLSchema-instance</a>&quot; xsi:schemaLocation=&quot;<a href="http://www.opengis.net/ogc">http://www.opengis.net/ogc</a> <a href="http://schemas.opengis.net/filter/1.0.0/filter.xsd">
http://schemas.opengis.net/filter/1.0.0/filter.xsd</a>&quot;&gt;<br>&nbsp;&nbsp; &lt;Intersects&gt;<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &lt;PropertyName&gt;Geometry&lt;/PropertyName&gt;<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &lt;gml:MultiPolygon xmlns:gml=&quot;<a href="http://www.opengis.net/gml">
http://www.opengis.net/gml</a>&quot; srsName=&quot;EPSG:28992&quot;&gt;<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &lt;gml:polygonMember&gt;<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &lt;gml:Polygon&gt;<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &lt;gml:outerBoundaryIs&gt;<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &lt;gml:LinearRing&gt;
<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &lt;gml:coordinates decimal=&quot;.&quot; cs=&quot;,&quot; ts=&quot; &quot;&gt;142031.375,451226.6875 142456.53125,449448.6875 142596.390625,448863.8125 142778.140625,448103.78125 142940.78125,447422.5
 142940.921875,447421.9375 142937.578125,447422.90625 141994.703125,447696.65625 141994.46875,447695.6875 141895.9375,447619.375 141846.328125,447576.84375 141832.140625,447520.15625 141828.609375,447424.46875 141839.234375
,447318.1875 141835.6875,447261.46875 141800.25,447151.625 141750.640625,447087.84375 141701.03125,446999.25 141669.140625,446949.65625 141495.515625,446471.28125 141445.140625,446360.96875 141441.546875,446361.65625 141273.34375
,446394 141269.921875,446394.65625 141026.8125,446440.71875 141023.53125,446441.8125 139836.65625,446835.46875 139836.046875,446834.5 139606.984375,446467.96875 139604.34375,446469.40625 138646.375,446996.3125 137102.234375
,448524.375 136506.828125,449142.25 136557.09375,449410.28125 136448.859375,449561 136130.171875,450003.71875 135952.984375,450464.1875 136569.34375,450617.625 137113.34375,450733.3125 137186.03125,450748.78125 137828.765625
,450879.875 138255.171875,450979.21875 138843.046875,451103.34375 139248.375,451197.1875 139280.15625,451202.28125 142031.375,451226.6875&lt;/gml:coordinates&gt;<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &lt;/gml:LinearRing&gt;<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &lt;/gml:outerBoundaryIs&gt;
<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &lt;/gml:Polygon&gt;<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &lt;/gml:polygonMember&gt;<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &lt;/gml:MultiPolygon&gt;<br>&nbsp;&nbsp; &lt;/Intersects&gt;<br>&lt;/Filter&gt;<br><br>TIA.<br><br>Best regards,<br>Bart<br>