Determine if a polygon is not enclosed by other polygons

.everyoneloves__top-leaderboard:empty,.everyoneloves__mid-leaderboard:empty margin-bottom:0;

up vote
2
down vote

favorite

enter image description here Does anyone know of a routine/function or transformer that would identify if a polygon is not fully enclosed/surrounded by other polygons?

The polygons have no gaps or slivers between them so if a certain section of polygon is not attached to another polygon then these are the ones I want to keep.

I am happy to use any of PostGIS/FME/QGIS I just cannot seem to find a function or transformer that would do it.

EDIT

I have added a picture to help explain:

I am after the purple polygons and ignore the yellow ones as they are not fully enclosed by other squares

edited Aug 6 at 17:04

Spacedman

19.6k13047

asked Aug 6 at 8:46

tjmgis

2,0951325

Are the polygons on the same layer? And if you worked with USA data then Texas would be selected but Colorado not, or?
â€“Â user30184
Aug 6 at 9:59

Yes the same layer. Thats exactly what I am after.
â€“Â tjmgis
Aug 6 at 10:44

I undestood that if a certain section of polygon is not attached to another polygon then these are the ones I want to keep means the purple polygons.
â€“Â user30184
Aug 6 at 11:05

sorry - yes your are right, i just edited. Mondays!
â€“Â tjmgis
Aug 6 at 11:11

1

So your question is "How can I tell if a feature in a set of features has any external boundaries?". Test the linear intersection of A with the dissolved features.
â€“Â Spacedman
Aug 6 at 13:18

Â |Â
show 3 more comments

up vote
2
down vote

favorite

enter image description here Does anyone know of a routine/function or transformer that would identify if a polygon is not fully enclosed/surrounded by other polygons?

The polygons have no gaps or slivers between them so if a certain section of polygon is not attached to another polygon then these are the ones I want to keep.

I am happy to use any of PostGIS/FME/QGIS I just cannot seem to find a function or transformer that would do it.

EDIT

I have added a picture to help explain:

I am after the purple polygons and ignore the yellow ones as they are not fully enclosed by other squares

edited Aug 6 at 17:04

Spacedman

19.6k13047

asked Aug 6 at 8:46

tjmgis

2,0951325

Are the polygons on the same layer? And if you worked with USA data then Texas would be selected but Colorado not, or?
â€“Â user30184
Aug 6 at 9:59

Yes the same layer. Thats exactly what I am after.
â€“Â tjmgis
Aug 6 at 10:44

I undestood that if a certain section of polygon is not attached to another polygon then these are the ones I want to keep means the purple polygons.
â€“Â user30184
Aug 6 at 11:05

sorry - yes your are right, i just edited. Mondays!
â€“Â tjmgis
Aug 6 at 11:11

1

So your question is "How can I tell if a feature in a set of features has any external boundaries?". Test the linear intersection of A with the dissolved features.
â€“Â Spacedman
Aug 6 at 13:18

Â |Â
show 3 more comments

up vote
2
down vote

favorite

enter image description here Does anyone know of a routine/function or transformer that would identify if a polygon is not fully enclosed/surrounded by other polygons?

The polygons have no gaps or slivers between them so if a certain section of polygon is not attached to another polygon then these are the ones I want to keep.

I am happy to use any of PostGIS/FME/QGIS I just cannot seem to find a function or transformer that would do it.

EDIT

I have added a picture to help explain:

I am after the purple polygons and ignore the yellow ones as they are not fully enclosed by other squares

edited Aug 6 at 17:04

Spacedman

19.6k13047

asked Aug 6 at 8:46

tjmgis

2,0951325

enter image description here Does anyone know of a routine/function or transformer that would identify if a polygon is not fully enclosed/surrounded by other polygons?

The polygons have no gaps or slivers between them so if a certain section of polygon is not attached to another polygon then these are the ones I want to keep.

I am happy to use any of PostGIS/FME/QGIS I just cannot seem to find a function or transformer that would do it.

EDIT

I have added a picture to help explain:

I am after the purple polygons and ignore the yellow ones as they are not fully enclosed by other squares

edited Aug 6 at 17:04

Spacedman

19.6k13047

asked Aug 6 at 8:46

tjmgis

2,0951325

edited Aug 6 at 17:04

Spacedman

19.6k13047

edited Aug 6 at 17:04

Spacedman

19.6k13047

edited Aug 6 at 17:04

Spacedman

19.6k13047

asked Aug 6 at 8:46

tjmgis

2,0951325

asked Aug 6 at 8:46

tjmgis

2,0951325

asked Aug 6 at 8:46

tjmgis

2,0951325

Are the polygons on the same layer? And if you worked with USA data then Texas would be selected but Colorado not, or?
â€“Â user30184
Aug 6 at 9:59

Yes the same layer. Thats exactly what I am after.
â€“Â tjmgis
Aug 6 at 10:44

I undestood that if a certain section of polygon is not attached to another polygon then these are the ones I want to keep means the purple polygons.
â€“Â user30184
Aug 6 at 11:05

sorry - yes your are right, i just edited. Mondays!
â€“Â tjmgis
Aug 6 at 11:11

1

So your question is "How can I tell if a feature in a set of features has any external boundaries?". Test the linear intersection of A with the dissolved features.
â€“Â Spacedman
Aug 6 at 13:18

Â |Â
show 3 more comments

Are the polygons on the same layer? And if you worked with USA data then Texas would be selected but Colorado not, or?
â€“Â user30184
Aug 6 at 9:59

Yes the same layer. Thats exactly what I am after.
â€“Â tjmgis
Aug 6 at 10:44

I undestood that if a certain section of polygon is not attached to another polygon then these are the ones I want to keep means the purple polygons.
â€“Â user30184
Aug 6 at 11:05

sorry - yes your are right, i just edited. Mondays!
â€“Â tjmgis
Aug 6 at 11:11

1

So your question is "How can I tell if a feature in a set of features has any external boundaries?". Test the linear intersection of A with the dissolved features.
â€“Â Spacedman
Aug 6 at 13:18

Are the polygons on the same layer? And if you worked with USA data then Texas would be selected but Colorado not, or?
â€“Â user30184
Aug 6 at 9:59

Yes the same layer. Thats exactly what I am after.
â€“Â tjmgis
Aug 6 at 10:44

I undestood that if a certain section of polygon is not attached to another polygon then these are the ones I want to keep means the purple polygons.
â€“Â user30184
Aug 6 at 11:05

sorry - yes your are right, i just edited. Mondays!
â€“Â tjmgis
Aug 6 at 11:11

So your question is "How can I tell if a feature in a set of features has any external boundaries?". Test the linear intersection of A with the dissolved features.
â€“Â Spacedman
Aug 6 at 13:18

Â |Â
show 3 more comments

4 Answers
4

active

oldest

votes

up vote
4
down vote

Assuming your polygons are topologically correct (no gaps, or overlaps), then a polygon A is enclosed by other polygons if its perimeter is equal to the sum of intersection length between its boundary and the boundary of other polygons. This method lends itself to use of a spatial index and avoids performing a relatively costly union of all other polygons.

SELECT a.id
FROM my_data a 
INNER JOIN my_data b ON (ST_Intersects(a.geom, b.geom) AND a.id != b.id) 
GROUP BY a.id
HAVING 1e-6 > 
 abs(ST_Length(ST_ExteriorRing(a.geom)) - 
 sum(ST_Length(ST_Intersection(ST_Exteriorring(a.geom), ST_ExteriorRing(b.geom)))));

edited Aug 6 at 20:33

answered Aug 6 at 19:52

dbaston

6,88722045

add a commentÂ |Â

up vote
2
down vote

I have an answer that should be logically correct but the SQL part is not complete. The idea is based on intersections of adjacent polygons. The intesections are either lines if polygons share segment(s), or points if polygons meet at one vertex.

The first step is to compute the intersection of one polygon and the union of all the polygons on the layer, except the selected polygon.

enter image description here

The intersection of polygon 1 and union of all the rest is a multilinestring with two members. The intersection of polygon 6 and anything else but polygon 6 consists of four linestrings members which are connected.

enter image description here

The second step is to try if it is possible to build a new polygon from the intersections. For the connected lines from polygon 6 it is possible but not for the two linestrings from polygon 1.

This is the SQL that I used for testing.

select ST_MakePolygon(ST_LineMerge(c.intersection)) from 
(select 
ST_Intersection(sub_a.a,sub_b.b) as intersection from
(select ST_Union("GEOMETRY") as a from outline_test where gid!=6) as sub_a,
(select "GEOMETRY" as b from outline_test where gid=6) as sub_b) as c;

With gid=1 the result is an error:

ERROR: lwpoly_from_lwlines: shell must have at least 4 points

and with gid=2:

ERROR: lwpoly_from_lwlines: shell must be closed

answered Aug 6 at 12:51

user30184

26.2k22749

add a commentÂ |Â

up vote
0
down vote

Another approach: The polygons you want to find have 1-dimensional intersection (lines) with the outer ring of the whole area as an union. If polygon is only touching the outer boundary the intersection is a point with dimension=0, and polygons which are totally inside to not intersect as all (so they are disjoint).

So you can select the features with ST_Dimension=1 from SQL query like

select outline_test.gid,ST_Dimension
(ST_Intersection(outline_test."GEOMETRY",sub_a.geometry))
from 
(select ST_ExteriorRing(ST_Union("GEOMETRY")) as geometry from outline_test)
as sub_a, outline_test;

The query does not make difference between "touches" and "disjoint" as documented in https://postgis.net/docs/ST_Dimension.html If the dimension is unknown (empty GEOMETRYCOLLECTION) 0 is returned but that does not matter in your use case. Rereffing to the data and image in my other answer, other polygons than numbers 5, 6, 8, and 10 should be found and that happens.

1;1
2;1
3;1
4;1
5;0
6;0
7;1
8;0
9;1
10;0
11;1
12;1
13;1
14;1

I am sure it is possible to write more elegant queries based on the outer boundary and DE-9IM relations http://postgis.net/docs/using_postgis_dbmanagement.html#DE-9IM.

answered Aug 6 at 13:30

user30184

26.2k22749

I think you've expressed the inverse problem with "The polygons you want to find have 1-dimensional intersection (lines) with the outer ring of the whole area as an union" - the required polygons are the yellow ones in the diagram and they don't have 1-d intersection with the external.
â€“Â Spacedman
Aug 6 at 15:45

OP wants the purple polygons.
â€“Â user30184
Aug 6 at 15:57

Oh then the title is misleading: "Determine if a polygon is enclosed by other polygons"... I'll edit..
â€“Â Spacedman
Aug 6 at 17:04

add a commentÂ |Â

up vote
0
down vote

So, in FME you would probably use the NeighborFinder and count the number of neighbors. If you included diagonal polygons, then you just need to check if there are 8 neighbors, but since you don't seem to need that, we need to check for 4 neighbors, but only horizontally/vertically.

So, here's a workspace to do that:

enter image description here

Basically:

Use a CentrePointReplacer to convert polygons to centre points

Use a NeighborFinder:
- Set Candidate Only Mode
- Neighbors to Find = 4
- Maximum distance = width/height of your polygons
- Closest Candidate List Name = list

Use a ListElementCounter to count the number of neighbors found

Use a Tester to test if _element_count < 4

Use a SpatialFilter to filter out the original data
- Tester:Passed = Filter, Original Data = Candidates
- Predicate to Test = "Filter is Within Candidate"

The SpatialFilter:Passed output is the data you need.

I put my demo workspace on Dropbox if you want to check it out.

edited Aug 7 at 19:26

answered Aug 7 at 14:50

Mark Ireland

8,99022259

add a commentÂ |Â

Your Answer

StackExchange.ready(function()
var channelOptions =
tags: "".split(" "),
id: "79"
;
initTagRenderer("".split(" "), "".split(" "), channelOptions);

StackExchange.using("externalEditor", function()
// Have to fire editor after snippets, if snippets enabled
if (StackExchange.settings.snippets.snippetsEnabled)
StackExchange.using("snippets", function()
createEditor();
);

else
createEditor();

);

function createEditor()
StackExchange.prepareEditor(
heartbeatType: 'answer',
convertImagesToLinks: false,
noModals: false,
showLowRepImageUploadWarning: true,
reputationToPostImages: null,
bindNavPrevention: true,
postfix: "",
onDemand: true,
discardSelector: ".discard-answer"
,immediatelyShowMarkdownHelp:true
);

);

draft saved

draft discarded

StackExchange.ready(
function ()
StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fgis.stackexchange.com%2fquestions%2f291824%2fdetermine-if-a-polygon-is-not-enclosed-by-other-polygons%23new-answer', 'question_page');

);

Post as a guest

Name

4 Answers
4

active

oldest

votes

4 Answers
4

active

oldest

votes

up vote
4
down vote

SELECT a.id
FROM my_data a 
INNER JOIN my_data b ON (ST_Intersects(a.geom, b.geom) AND a.id != b.id) 
GROUP BY a.id
HAVING 1e-6 > 
 abs(ST_Length(ST_ExteriorRing(a.geom)) - 
 sum(ST_Length(ST_Intersection(ST_Exteriorring(a.geom), ST_ExteriorRing(b.geom)))));

edited Aug 6 at 20:33

answered Aug 6 at 19:52

dbaston

6,88722045

add a commentÂ |Â

up vote
4
down vote

SELECT a.id
FROM my_data a 
INNER JOIN my_data b ON (ST_Intersects(a.geom, b.geom) AND a.id != b.id) 
GROUP BY a.id
HAVING 1e-6 > 
 abs(ST_Length(ST_ExteriorRing(a.geom)) - 
 sum(ST_Length(ST_Intersection(ST_Exteriorring(a.geom), ST_ExteriorRing(b.geom)))));

edited Aug 6 at 20:33

answered Aug 6 at 19:52

dbaston

6,88722045

add a commentÂ |Â

up vote
4
down vote

SELECT a.id
FROM my_data a 
INNER JOIN my_data b ON (ST_Intersects(a.geom, b.geom) AND a.id != b.id) 
GROUP BY a.id
HAVING 1e-6 > 
 abs(ST_Length(ST_ExteriorRing(a.geom)) - 
 sum(ST_Length(ST_Intersection(ST_Exteriorring(a.geom), ST_ExteriorRing(b.geom)))));

edited Aug 6 at 20:33

answered Aug 6 at 19:52

dbaston

6,88722045

SELECT a.id
FROM my_data a 
INNER JOIN my_data b ON (ST_Intersects(a.geom, b.geom) AND a.id != b.id) 
GROUP BY a.id
HAVING 1e-6 > 
 abs(ST_Length(ST_ExteriorRing(a.geom)) - 
 sum(ST_Length(ST_Intersection(ST_Exteriorring(a.geom), ST_ExteriorRing(b.geom)))));

edited Aug 6 at 20:33

answered Aug 6 at 19:52

dbaston

6,88722045

edited Aug 6 at 20:33

answered Aug 6 at 19:52

dbaston

6,88722045

answered Aug 6 at 19:52

dbaston

6,88722045

answered Aug 6 at 19:52

dbaston

6,88722045

add a commentÂ |Â

up vote
2
down vote

The first step is to compute the intersection of one polygon and the union of all the polygons on the layer, except the selected polygon.

enter image description here

The second step is to try if it is possible to build a new polygon from the intersections. For the connected lines from polygon 6 it is possible but not for the two linestrings from polygon 1.

This is the SQL that I used for testing.

select ST_MakePolygon(ST_LineMerge(c.intersection)) from 
(select 
ST_Intersection(sub_a.a,sub_b.b) as intersection from
(select ST_Union("GEOMETRY") as a from outline_test where gid!=6) as sub_a,
(select "GEOMETRY" as b from outline_test where gid=6) as sub_b) as c;

With gid=1 the result is an error:

ERROR: lwpoly_from_lwlines: shell must have at least 4 points

and with gid=2:

ERROR: lwpoly_from_lwlines: shell must be closed

answered Aug 6 at 12:51

user30184

26.2k22749

add a commentÂ |Â

up vote
2
down vote

The first step is to compute the intersection of one polygon and the union of all the polygons on the layer, except the selected polygon.

enter image description here

The second step is to try if it is possible to build a new polygon from the intersections. For the connected lines from polygon 6 it is possible but not for the two linestrings from polygon 1.

This is the SQL that I used for testing.

select ST_MakePolygon(ST_LineMerge(c.intersection)) from 
(select 
ST_Intersection(sub_a.a,sub_b.b) as intersection from
(select ST_Union("GEOMETRY") as a from outline_test where gid!=6) as sub_a,
(select "GEOMETRY" as b from outline_test where gid=6) as sub_b) as c;

With gid=1 the result is an error:

ERROR: lwpoly_from_lwlines: shell must have at least 4 points

and with gid=2:

ERROR: lwpoly_from_lwlines: shell must be closed

answered Aug 6 at 12:51

user30184

26.2k22749

add a commentÂ |Â

up vote
2
down vote

The first step is to compute the intersection of one polygon and the union of all the polygons on the layer, except the selected polygon.

enter image description here

The second step is to try if it is possible to build a new polygon from the intersections. For the connected lines from polygon 6 it is possible but not for the two linestrings from polygon 1.

This is the SQL that I used for testing.

select ST_MakePolygon(ST_LineMerge(c.intersection)) from 
(select 
ST_Intersection(sub_a.a,sub_b.b) as intersection from
(select ST_Union("GEOMETRY") as a from outline_test where gid!=6) as sub_a,
(select "GEOMETRY" as b from outline_test where gid=6) as sub_b) as c;

With gid=1 the result is an error:

ERROR: lwpoly_from_lwlines: shell must have at least 4 points

and with gid=2:

ERROR: lwpoly_from_lwlines: shell must be closed

answered Aug 6 at 12:51

user30184

26.2k22749

The first step is to compute the intersection of one polygon and the union of all the polygons on the layer, except the selected polygon.

enter image description here

The second step is to try if it is possible to build a new polygon from the intersections. For the connected lines from polygon 6 it is possible but not for the two linestrings from polygon 1.

This is the SQL that I used for testing.

select ST_MakePolygon(ST_LineMerge(c.intersection)) from 
(select 
ST_Intersection(sub_a.a,sub_b.b) as intersection from
(select ST_Union("GEOMETRY") as a from outline_test where gid!=6) as sub_a,
(select "GEOMETRY" as b from outline_test where gid=6) as sub_b) as c;

With gid=1 the result is an error:

ERROR: lwpoly_from_lwlines: shell must have at least 4 points

and with gid=2:

ERROR: lwpoly_from_lwlines: shell must be closed

answered Aug 6 at 12:51

user30184

26.2k22749

answered Aug 6 at 12:51

user30184

26.2k22749

answered Aug 6 at 12:51

user30184

26.2k22749

answered Aug 6 at 12:51

user30184

26.2k22749

add a commentÂ |Â

up vote
0
down vote

So you can select the features with ST_Dimension=1 from SQL query like

select outline_test.gid,ST_Dimension
(ST_Intersection(outline_test."GEOMETRY",sub_a.geometry))
from 
(select ST_ExteriorRing(ST_Union("GEOMETRY")) as geometry from outline_test)
as sub_a, outline_test;

1;1
2;1
3;1
4;1
5;0
6;0
7;1
8;0
9;1
10;0
11;1
12;1
13;1
14;1

I am sure it is possible to write more elegant queries based on the outer boundary and DE-9IM relations http://postgis.net/docs/using_postgis_dbmanagement.html#DE-9IM.

answered Aug 6 at 13:30

user30184

26.2k22749

I think you've expressed the inverse problem with "The polygons you want to find have 1-dimensional intersection (lines) with the outer ring of the whole area as an union" - the required polygons are the yellow ones in the diagram and they don't have 1-d intersection with the external.
â€“Â Spacedman
Aug 6 at 15:45

OP wants the purple polygons.
â€“Â user30184
Aug 6 at 15:57

Oh then the title is misleading: "Determine if a polygon is enclosed by other polygons"... I'll edit..
â€“Â Spacedman
Aug 6 at 17:04

add a commentÂ |Â

up vote
0
down vote

So you can select the features with ST_Dimension=1 from SQL query like

select outline_test.gid,ST_Dimension
(ST_Intersection(outline_test."GEOMETRY",sub_a.geometry))
from 
(select ST_ExteriorRing(ST_Union("GEOMETRY")) as geometry from outline_test)
as sub_a, outline_test;

1;1
2;1
3;1
4;1
5;0
6;0
7;1
8;0
9;1
10;0
11;1
12;1
13;1
14;1

I am sure it is possible to write more elegant queries based on the outer boundary and DE-9IM relations http://postgis.net/docs/using_postgis_dbmanagement.html#DE-9IM.

answered Aug 6 at 13:30

user30184

26.2k22749

I think you've expressed the inverse problem with "The polygons you want to find have 1-dimensional intersection (lines) with the outer ring of the whole area as an union" - the required polygons are the yellow ones in the diagram and they don't have 1-d intersection with the external.
â€“Â Spacedman
Aug 6 at 15:45

OP wants the purple polygons.
â€“Â user30184
Aug 6 at 15:57

Oh then the title is misleading: "Determine if a polygon is enclosed by other polygons"... I'll edit..
â€“Â Spacedman
Aug 6 at 17:04

add a commentÂ |Â

up vote
0
down vote

So you can select the features with ST_Dimension=1 from SQL query like

select outline_test.gid,ST_Dimension
(ST_Intersection(outline_test."GEOMETRY",sub_a.geometry))
from 
(select ST_ExteriorRing(ST_Union("GEOMETRY")) as geometry from outline_test)
as sub_a, outline_test;

1;1
2;1
3;1
4;1
5;0
6;0
7;1
8;0
9;1
10;0
11;1
12;1
13;1
14;1

I am sure it is possible to write more elegant queries based on the outer boundary and DE-9IM relations http://postgis.net/docs/using_postgis_dbmanagement.html#DE-9IM.

answered Aug 6 at 13:30

user30184

26.2k22749

So you can select the features with ST_Dimension=1 from SQL query like

select outline_test.gid,ST_Dimension
(ST_Intersection(outline_test."GEOMETRY",sub_a.geometry))
from 
(select ST_ExteriorRing(ST_Union("GEOMETRY")) as geometry from outline_test)
as sub_a, outline_test;

1;1
2;1
3;1
4;1
5;0
6;0
7;1
8;0
9;1
10;0
11;1
12;1
13;1
14;1

I am sure it is possible to write more elegant queries based on the outer boundary and DE-9IM relations http://postgis.net/docs/using_postgis_dbmanagement.html#DE-9IM.

answered Aug 6 at 13:30

user30184

26.2k22749

answered Aug 6 at 13:30

user30184

26.2k22749

answered Aug 6 at 13:30

user30184

26.2k22749

answered Aug 6 at 13:30

user30184

26.2k22749

I think you've expressed the inverse problem with "The polygons you want to find have 1-dimensional intersection (lines) with the outer ring of the whole area as an union" - the required polygons are the yellow ones in the diagram and they don't have 1-d intersection with the external.
â€“Â Spacedman
Aug 6 at 15:45

OP wants the purple polygons.
â€“Â user30184
Aug 6 at 15:57

Oh then the title is misleading: "Determine if a polygon is enclosed by other polygons"... I'll edit..
â€“Â Spacedman
Aug 6 at 17:04

add a commentÂ |Â

I think you've expressed the inverse problem with "The polygons you want to find have 1-dimensional intersection (lines) with the outer ring of the whole area as an union" - the required polygons are the yellow ones in the diagram and they don't have 1-d intersection with the external.
â€“Â Spacedman
Aug 6 at 15:45

OP wants the purple polygons.
â€“Â user30184
Aug 6 at 15:57

Oh then the title is misleading: "Determine if a polygon is enclosed by other polygons"... I'll edit..
â€“Â Spacedman
Aug 6 at 17:04

I think you've expressed the inverse problem with "The polygons you want to find have 1-dimensional intersection (lines) with the outer ring of the whole area as an union" - the required polygons are the yellow ones in the diagram and they don't have 1-d intersection with the external.
â€“Â Spacedman
Aug 6 at 15:45

OP wants the purple polygons.
â€“Â user30184
Aug 6 at 15:57

Oh then the title is misleading: "Determine if a polygon is enclosed by other polygons"... I'll edit..
â€“Â Spacedman
Aug 6 at 17:04

add a commentÂ |Â

up vote
0
down vote

So, here's a workspace to do that:

enter image description here

Basically:

Use a CentrePointReplacer to convert polygons to centre points

Use a NeighborFinder:
- Set Candidate Only Mode
- Neighbors to Find = 4
- Maximum distance = width/height of your polygons
- Closest Candidate List Name = list

Use a ListElementCounter to count the number of neighbors found

Use a Tester to test if _element_count < 4

Use a SpatialFilter to filter out the original data
- Tester:Passed = Filter, Original Data = Candidates
- Predicate to Test = "Filter is Within Candidate"

The SpatialFilter:Passed output is the data you need.

I put my demo workspace on Dropbox if you want to check it out.

edited Aug 7 at 19:26

answered Aug 7 at 14:50

Mark Ireland

8,99022259

add a commentÂ |Â

up vote
0
down vote

So, here's a workspace to do that:

enter image description here

Basically:

Use a CentrePointReplacer to convert polygons to centre points

Use a NeighborFinder:
- Set Candidate Only Mode
- Neighbors to Find = 4
- Maximum distance = width/height of your polygons
- Closest Candidate List Name = list

Use a ListElementCounter to count the number of neighbors found

Use a Tester to test if _element_count < 4

Use a SpatialFilter to filter out the original data
- Tester:Passed = Filter, Original Data = Candidates
- Predicate to Test = "Filter is Within Candidate"

The SpatialFilter:Passed output is the data you need.

I put my demo workspace on Dropbox if you want to check it out.

edited Aug 7 at 19:26

answered Aug 7 at 14:50

Mark Ireland

8,99022259

add a commentÂ |Â

up vote
0
down vote

So, here's a workspace to do that:

enter image description here

Basically:

Use a CentrePointReplacer to convert polygons to centre points

Use a NeighborFinder:
- Set Candidate Only Mode
- Neighbors to Find = 4
- Maximum distance = width/height of your polygons
- Closest Candidate List Name = list

Use a ListElementCounter to count the number of neighbors found

Use a Tester to test if _element_count < 4

Use a SpatialFilter to filter out the original data
- Tester:Passed = Filter, Original Data = Candidates
- Predicate to Test = "Filter is Within Candidate"

The SpatialFilter:Passed output is the data you need.

I put my demo workspace on Dropbox if you want to check it out.

edited Aug 7 at 19:26

answered Aug 7 at 14:50

Mark Ireland

8,99022259

So, here's a workspace to do that:

enter image description here

Basically:

Use a CentrePointReplacer to convert polygons to centre points

Use a NeighborFinder:
- Set Candidate Only Mode
- Neighbors to Find = 4
- Maximum distance = width/height of your polygons
- Closest Candidate List Name = list

Use a ListElementCounter to count the number of neighbors found

Use a Tester to test if _element_count < 4

Use a SpatialFilter to filter out the original data
- Tester:Passed = Filter, Original Data = Candidates
- Predicate to Test = "Filter is Within Candidate"

The SpatialFilter:Passed output is the data you need.

I put my demo workspace on Dropbox if you want to check it out.

edited Aug 7 at 19:26

answered Aug 7 at 14:50

Mark Ireland

8,99022259

edited Aug 7 at 19:26

answered Aug 7 at 14:50

Mark Ireland

8,99022259

answered Aug 7 at 14:50

Mark Ireland

8,99022259

answered Aug 7 at 14:50

Mark Ireland

8,99022259

add a commentÂ |Â

draft saved

draft discarded

draft saved

draft discarded

Post as a guest

Name

Search This Blog

Iyfjky