Intersection function returns multiple coordinates Planned maintenance scheduled April 17/18, 2019 at 00:00UTC (8:00pm US/Eastern) Announcing the arrival of Valued Associate #679: Cesar Manara Unicorn Meta Zoo #1: Why another podcast?Removing intersect of one polygon from another using R?How to do a “fuzzy” intersection?Oracle Spatial Returns Different Intersection Results with SRID DefinedLine vs. Polygon Intersection CoordinatesHow to get the Area of two intersecting polygons on PostGIS?Getting intersection of multiple polygons efficiently in PythonRecover M-values from PostGIS `ST_Intersection`?Aggregate version of ST_IntersectionHow to write multiple .las files from polygon clipping with PDAL?Dealing with a hole (or inner ring) at the intersection of two MultiPolygon polygons in Python Shapely
How come Sam didn't become Lord of Horn Hill?
What do you call the main part of a joke?
Compare a given version number in the form major.minor.build.patch and see if one is less than the other
Most bit efficient text communication method?
What is the longest distance a player character can jump in one leap?
Denied boarding although I have proper visa and documentation. To whom should I make a complaint?
List of Python versions
Is there any way for the UK Prime Minister to make a motion directly dependent on Government confidence?
How to find 'n' nodes where all distances between them are greater than 'k'?
What does the "x" in "x86" represent?
How to answer "Have you ever been terminated?"
Muzzle or not? HOw would you name this thing?
Short story by Isaac Asimov about a mental chef
The logistics of corpse disposal
Fundamental Solution of the Pell Equation
Sci-Fi book where patients in a coma ward all live in a subconscious world linked together
Why did the US and UK choose different solutions to the problem of an undemocratic upper house?
How to override model in magento2?
Why wasn't DOSKEY integrated with command.com?
How to tell that you are a giant?
Delete nth line from bottom
For a new assistant professor in CS, how to build/manage a publication pipeline
Can an alien society believe that their star system is the universe?
Using audio cues to encourage good posture
Intersection function returns multiple coordinates
Planned maintenance scheduled April 17/18, 2019 at 00:00UTC (8:00pm US/Eastern)
Announcing the arrival of Valued Associate #679: Cesar Manara
Unicorn Meta Zoo #1: Why another podcast?Removing intersect of one polygon from another using R?How to do a “fuzzy” intersection?Oracle Spatial Returns Different Intersection Results with SRID DefinedLine vs. Polygon Intersection CoordinatesHow to get the Area of two intersecting polygons on PostGIS?Getting intersection of multiple polygons efficiently in PythonRecover M-values from PostGIS `ST_Intersection`?Aggregate version of ST_IntersectionHow to write multiple .las files from polygon clipping with PDAL?Dealing with a hole (or inner ring) at the intersection of two MultiPolygon polygons in Python Shapely
.everyoneloves__top-leaderboard:empty,.everyoneloves__mid-leaderboard:empty,.everyoneloves__bot-mid-leaderboard:empty margin-bottom:0;
I am trying to find intersection vertices of two polygons defined below
PG = Polygon([(0.0,0.0), (0.0,4.0), (2.0, 0.0)])
QG = shapely.affinity.rotate(PG, 10.0, origin="centroid")
print(PG.intersection(QG))
It returns a polygon which represents the intersection vertices but I get multiple vertices belonging to polygon PG and QG and not the exact intersection vertices something like this
POLYGON ((0 1.275007557649384, 0 3.233068398870882, 1.176669689726421 1.646660620547158, 1.901423900453083 0.1971521990938349, 0.7833182180345657 0, 0.2248182330207899 0, 0 1.275007557649384))
The graph looks something like this
How to get the exact coordinates of intersection vertices A, B and C as shown above?
python coordinates intersection shapely
asked Apr 10 at 6:21
Ashutosh MishraAshutosh Mishra
103
|
show 2 more comments
I am trying to find intersection vertices of two polygons defined below
PG = Polygon([(0.0,0.0), (0.0,4.0), (2.0, 0.0)])
QG = shapely.affinity.rotate(PG, 10.0, origin="centroid")
print(PG.intersection(QG))
It returns a polygon which represents the intersection vertices but I get multiple vertices belonging to polygon PG and QG and not the exact intersection vertices something like this
POLYGON ((0 1.275007557649384, 0 3.233068398870882, 1.176669689726421 1.646660620547158, 1.901423900453083 0.1971521990938349, 0.7833182180345657 0, 0.2248182330207899 0, 0 1.275007557649384))
The graph looks something like this
How to get the exact coordinates of intersection vertices A, B and C as shown above?
python coordinates intersection shapely
asked Apr 10 at 6:21
Ashutosh MishraAshutosh Mishra
103
1
Intersection of two overlapping polygons is a polygon. Are you perhaps more interested in finding the intersection of the boundaries of the polygons (linear rings) that would give a multipoint in this case?
– user30184
Apr 10 at 7:29
I am interested in getting just the intersection coordinates of polygons PG and QG.
– Ashutosh Mishra
Apr 10 at 7:38
the vertices of the polygon are the points you want, aren't they? other wise can you mark on your image which points you do want
– Ian Turton♦
Apr 10 at 7:42
There are 3 intersection coordinates if you see the above graph where both the polygons intersect. I want to find the coordinates of those 3 points.
– Ashutosh Mishra
Apr 10 at 7:48
I can see 6 points where the boundaries of the two polygons intersect. Could you mark your favorites? Polygons intersect topologically also in all interior points which are common for both and obviously you are not interested in that.
– user30184
Apr 10 at 8:11
|
show 2 more comments
I am trying to find intersection vertices of two polygons defined below
PG = Polygon([(0.0,0.0), (0.0,4.0), (2.0, 0.0)])
QG = shapely.affinity.rotate(PG, 10.0, origin="centroid")
print(PG.intersection(QG))
It returns a polygon which represents the intersection vertices but I get multiple vertices belonging to polygon PG and QG and not the exact intersection vertices something like this
POLYGON ((0 1.275007557649384, 0 3.233068398870882, 1.176669689726421 1.646660620547158, 1.901423900453083 0.1971521990938349, 0.7833182180345657 0, 0.2248182330207899 0, 0 1.275007557649384))
The graph looks something like this
How to get the exact coordinates of intersection vertices A, B and C as shown above?
python coordinates intersection shapely
asked Apr 10 at 6:21
Ashutosh MishraAshutosh Mishra
103
I am trying to find intersection vertices of two polygons defined below
PG = Polygon([(0.0,0.0), (0.0,4.0), (2.0, 0.0)])
QG = shapely.affinity.rotate(PG, 10.0, origin="centroid")
print(PG.intersection(QG))
It returns a polygon which represents the intersection vertices but I get multiple vertices belonging to polygon PG and QG and not the exact intersection vertices something like this
POLYGON ((0 1.275007557649384, 0 3.233068398870882, 1.176669689726421 1.646660620547158, 1.901423900453083 0.1971521990938349, 0.7833182180345657 0, 0.2248182330207899 0, 0 1.275007557649384))
The graph looks something like this
How to get the exact coordinates of intersection vertices A, B and C as shown above?
python coordinates intersection shapely
python coordinates intersection shapely
asked Apr 10 at 6:21
Ashutosh MishraAshutosh Mishra
103
asked Apr 10 at 6:21
Ashutosh MishraAshutosh Mishra
103
edited Apr 10 at 8:08
Ashutosh Mishra
asked Apr 10 at 6:21
Ashutosh MishraAshutosh Mishra
103
asked Apr 10 at 6:21
Ashutosh MishraAshutosh Mishra
103
asked Apr 10 at 6:21
Ashutosh MishraAshutosh Mishra
103
103
1
Intersection of two overlapping polygons is a polygon. Are you perhaps more interested in finding the intersection of the boundaries of the polygons (linear rings) that would give a multipoint in this case?
– user30184
Apr 10 at 7:29
I am interested in getting just the intersection coordinates of polygons PG and QG.
– Ashutosh Mishra
Apr 10 at 7:38
the vertices of the polygon are the points you want, aren't they? other wise can you mark on your image which points you do want
– Ian Turton♦
Apr 10 at 7:42
There are 3 intersection coordinates if you see the above graph where both the polygons intersect. I want to find the coordinates of those 3 points.
– Ashutosh Mishra
Apr 10 at 7:48
I can see 6 points where the boundaries of the two polygons intersect. Could you mark your favorites? Polygons intersect topologically also in all interior points which are common for both and obviously you are not interested in that.
– user30184
Apr 10 at 8:11
|
show 2 more comments
1
Intersection of two overlapping polygons is a polygon. Are you perhaps more interested in finding the intersection of the boundaries of the polygons (linear rings) that would give a multipoint in this case?
– user30184
Apr 10 at 7:29
I am interested in getting just the intersection coordinates of polygons PG and QG.
– Ashutosh Mishra
Apr 10 at 7:38
the vertices of the polygon are the points you want, aren't they? other wise can you mark on your image which points you do want
– Ian Turton♦
Apr 10 at 7:42
There are 3 intersection coordinates if you see the above graph where both the polygons intersect. I want to find the coordinates of those 3 points.
– Ashutosh Mishra
Apr 10 at 7:48
I can see 6 points where the boundaries of the two polygons intersect. Could you mark your favorites? Polygons intersect topologically also in all interior points which are common for both and obviously you are not interested in that.
– user30184
Apr 10 at 8:11
1
1
Intersection of two overlapping polygons is a polygon. Are you perhaps more interested in finding the intersection of the boundaries of the polygons (linear rings) that would give a multipoint in this case?
– user30184
Apr 10 at 7:29
Intersection of two overlapping polygons is a polygon. Are you perhaps more interested in finding the intersection of the boundaries of the polygons (linear rings) that would give a multipoint in this case?
– user30184
Apr 10 at 7:29
I am interested in getting just the intersection coordinates of polygons PG and QG.
– Ashutosh Mishra
Apr 10 at 7:38
I am interested in getting just the intersection coordinates of polygons PG and QG.
– Ashutosh Mishra
Apr 10 at 7:38
the vertices of the polygon are the points you want, aren't they? other wise can you mark on your image which points you do want
– Ian Turton♦
Apr 10 at 7:42
the vertices of the polygon are the points you want, aren't they? other wise can you mark on your image which points you do want
– Ian Turton♦
Apr 10 at 7:42
There are 3 intersection coordinates if you see the above graph where both the polygons intersect. I want to find the coordinates of those 3 points.
– Ashutosh Mishra
Apr 10 at 7:48
There are 3 intersection coordinates if you see the above graph where both the polygons intersect. I want to find the coordinates of those 3 points.
– Ashutosh Mishra
Apr 10 at 7:48
I can see 6 points where the boundaries of the two polygons intersect. Could you mark your favorites? Polygons intersect topologically also in all interior points which are common for both and obviously you are not interested in that.
– user30184
Apr 10 at 8:11
I can see 6 points where the boundaries of the two polygons intersect. Could you mark your favorites? Polygons intersect topologically also in all interior points which are common for both and obviously you are not interested in that.
– user30184
Apr 10 at 8:11
|
show 2 more comments
1 Answer
1
active
oldest
votes
As user30184 says, the intersection of two overlapping polygons is a polygon and with your result (light blue polygon)
The intersection of the boundaries of the polygons (linearRings) gives a multipoint in you case
one = LineString(list(PG.exterior.coords))
two = LineString(list(QG.exterior.coords))
print(one.intersection(two).wkt)
MULTIPOINT (0 1.275007557649384, 0 3.233068398870882, 0.2248182330207899 0, 0.7833182180345657 0, 1.176669689726421 1.646660620547158, 1.901423900453083 0.1971521990938349)
These points are also the vertices of the resulting polygon
points = [Point(pt) for pt in list(PG.intersection(QG).exterior.coords)[:-1]]
for i in points:
print(i.wkt)
POINT (0 1.275007557649384)
POINT (0 3.233068398870882)
POINT (1.176669689726421 1.646660620547158)
POINT (1.901423900453083 0.1971521990938349)
POINT (0.7833182180345657 0)
POINT (0.2248182330207899 0)
But shapely alone does not allow to select the A,B,C points from the others.
answered Apr 10 at 17:05
genegene
37.6k157122
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',
autoActivateHeartbeat: false,
convertImagesToLinks: false,
noModals: true,
showLowRepImageUploadWarning: true,
reputationToPostImages: null,
bindNavPrevention: true,
postfix: "",
imageUploader:
brandingHtml: "Powered by u003ca class="icon-imgur-white" href="https://imgur.com/"u003eu003c/au003e",
contentPolicyHtml: "User contributions licensed under u003ca href="https://creativecommons.org/licenses/by-sa/3.0/"u003ecc by-sa 3.0 with attribution requiredu003c/au003e u003ca href="https://stackoverflow.com/legal/content-policy"u003e(content policy)u003c/au003e",
allowUrls: true
,
onDemand: true,
discardSelector: ".discard-answer"
,immediatelyShowMarkdownHelp:true
);
);
Sign up or log in
StackExchange.ready(function ()
StackExchange.helpers.onClickDraftSave('#login-link');
);
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Post as a guest
Required, but never shown
StackExchange.ready(
function ()
StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fgis.stackexchange.com%2fquestions%2f318325%2fintersection-function-returns-multiple-coordinates%23new-answer', 'question_page');
);
Post as a guest
Required, but never shown
1 Answer
1
active
oldest
votes
1 Answer
1
active
oldest
votes
active
oldest
votes
active
oldest
votes
As user30184 says, the intersection of two overlapping polygons is a polygon and with your result (light blue polygon)
The intersection of the boundaries of the polygons (linearRings) gives a multipoint in you case
one = LineString(list(PG.exterior.coords))
two = LineString(list(QG.exterior.coords))
print(one.intersection(two).wkt)
MULTIPOINT (0 1.275007557649384, 0 3.233068398870882, 0.2248182330207899 0, 0.7833182180345657 0, 1.176669689726421 1.646660620547158, 1.901423900453083 0.1971521990938349)
These points are also the vertices of the resulting polygon
points = [Point(pt) for pt in list(PG.intersection(QG).exterior.coords)[:-1]]
for i in points:
print(i.wkt)
POINT (0 1.275007557649384)
POINT (0 3.233068398870882)
POINT (1.176669689726421 1.646660620547158)
POINT (1.901423900453083 0.1971521990938349)
POINT (0.7833182180345657 0)
POINT (0.2248182330207899 0)
But shapely alone does not allow to select the A,B,C points from the others.
answered Apr 10 at 17:05
genegene
37.6k157122
add a comment |
As user30184 says, the intersection of two overlapping polygons is a polygon and with your result (light blue polygon)
The intersection of the boundaries of the polygons (linearRings) gives a multipoint in you case
one = LineString(list(PG.exterior.coords))
two = LineString(list(QG.exterior.coords))
print(one.intersection(two).wkt)
MULTIPOINT (0 1.275007557649384, 0 3.233068398870882, 0.2248182330207899 0, 0.7833182180345657 0, 1.176669689726421 1.646660620547158, 1.901423900453083 0.1971521990938349)
These points are also the vertices of the resulting polygon
points = [Point(pt) for pt in list(PG.intersection(QG).exterior.coords)[:-1]]
for i in points:
print(i.wkt)
POINT (0 1.275007557649384)
POINT (0 3.233068398870882)
POINT (1.176669689726421 1.646660620547158)
POINT (1.901423900453083 0.1971521990938349)
POINT (0.7833182180345657 0)
POINT (0.2248182330207899 0)
But shapely alone does not allow to select the A,B,C points from the others.
answered Apr 10 at 17:05
genegene
37.6k157122
add a comment |
As user30184 says, the intersection of two overlapping polygons is a polygon and with your result (light blue polygon)
The intersection of the boundaries of the polygons (linearRings) gives a multipoint in you case
one = LineString(list(PG.exterior.coords))
two = LineString(list(QG.exterior.coords))
print(one.intersection(two).wkt)
MULTIPOINT (0 1.275007557649384, 0 3.233068398870882, 0.2248182330207899 0, 0.7833182180345657 0, 1.176669689726421 1.646660620547158, 1.901423900453083 0.1971521990938349)
These points are also the vertices of the resulting polygon
points = [Point(pt) for pt in list(PG.intersection(QG).exterior.coords)[:-1]]
for i in points:
print(i.wkt)
POINT (0 1.275007557649384)
POINT (0 3.233068398870882)
POINT (1.176669689726421 1.646660620547158)
POINT (1.901423900453083 0.1971521990938349)
POINT (0.7833182180345657 0)
POINT (0.2248182330207899 0)
But shapely alone does not allow to select the A,B,C points from the others.
answered Apr 10 at 17:05
genegene
37.6k157122
As user30184 says, the intersection of two overlapping polygons is a polygon and with your result (light blue polygon)
The intersection of the boundaries of the polygons (linearRings) gives a multipoint in you case
one = LineString(list(PG.exterior.coords))
two = LineString(list(QG.exterior.coords))
print(one.intersection(two).wkt)
MULTIPOINT (0 1.275007557649384, 0 3.233068398870882, 0.2248182330207899 0, 0.7833182180345657 0, 1.176669689726421 1.646660620547158, 1.901423900453083 0.1971521990938349)
These points are also the vertices of the resulting polygon
points = [Point(pt) for pt in list(PG.intersection(QG).exterior.coords)[:-1]]
for i in points:
print(i.wkt)
POINT (0 1.275007557649384)
POINT (0 3.233068398870882)
POINT (1.176669689726421 1.646660620547158)
POINT (1.901423900453083 0.1971521990938349)
POINT (0.7833182180345657 0)
POINT (0.2248182330207899 0)
But shapely alone does not allow to select the A,B,C points from the others.
answered Apr 10 at 17:05
genegene
37.6k157122
edited Apr 10 at 17:12
answered Apr 10 at 17:05
genegene
37.6k157122
answered Apr 10 at 17:05
genegene
37.6k157122
answered Apr 10 at 17:05
genegene
37.6k157122
37.6k157122
add a comment |
add a comment |
Thanks for contributing an answer to Geographic Information Systems Stack Exchange!
- Please be sure to answer the question. Provide details and share your research!
But avoid …
- Asking for help, clarification, or responding to other answers.
- Making statements based on opinion; back them up with references or personal experience.
To learn more, see our tips on writing great answers.
Sign up or log in
StackExchange.ready(function ()
StackExchange.helpers.onClickDraftSave('#login-link');
);
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Post as a guest
Required, but never shown
StackExchange.ready(
function ()
StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fgis.stackexchange.com%2fquestions%2f318325%2fintersection-function-returns-multiple-coordinates%23new-answer', 'question_page');
);
Post as a guest
Required, but never shown
Sign up or log in
StackExchange.ready(function ()
StackExchange.helpers.onClickDraftSave('#login-link');
);
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Post as a guest
Required, but never shown
Sign up or log in
StackExchange.ready(function ()
StackExchange.helpers.onClickDraftSave('#login-link');
);
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Post as a guest
Required, but never shown
Sign up or log in
StackExchange.ready(function ()
StackExchange.helpers.onClickDraftSave('#login-link');
);
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Post as a guest
Required, but never shown
Required, but never shown
Required, but never shown
Required, but never shown
Required, but never shown
Required, but never shown
Required, but never shown
Required, but never shown
Required, but never shown
1
Intersection of two overlapping polygons is a polygon. Are you perhaps more interested in finding the intersection of the boundaries of the polygons (linear rings) that would give a multipoint in this case?
– user30184
Apr 10 at 7:29
I am interested in getting just the intersection coordinates of polygons PG and QG.
– Ashutosh Mishra
Apr 10 at 7:38
the vertices of the polygon are the points you want, aren't they? other wise can you mark on your image which points you do want
– Ian Turton♦
Apr 10 at 7:42
There are 3 intersection coordinates if you see the above graph where both the polygons intersect. I want to find the coordinates of those 3 points.
– Ashutosh Mishra
Apr 10 at 7:48
I can see 6 points where the boundaries of the two polygons intersect. Could you mark your favorites? Polygons intersect topologically also in all interior points which are common for both and obviously you are not interested in that.
– user30184
Apr 10 at 8:11