Touching Polygon Intersection The Next CEO of Stack OverflowCalculating intersection between points and polygon in shapefileDetermine if point is within an irregular polygon using PythonAdding a Property (attribute) to a geometry in Shapely/FionaIntersection postgis TopologyExceptionLine vs. Polygon Intersection CoordinatesMulti-part Geometrical Intersection errorPolygon intersection with LineString in PostGisEditing touching polygons simultaneously?Geopandas Line Polygon IntersectionHow is intersection implemented in Shapely?
How to make a software documentation "officially" citable?
Are there languages with no euphemisms?
Is it okay to store user locations?
Inappropriate reference requests from Journal reviewers
How can I get through very long and very dry, but also very useful technical documents when learning a new tool?
How long to clear the 'suck zone' of a turbofan after start is initiated?
How to be diplomatic in refusing to write code that breaches the privacy of our users
How to safely derail a train during transit?
What happens if you roll doubles 3 times then land on "Go to jail?"
Too much space between section and text in a twocolumn document
Increase performance creating Mandelbrot set in python
Visit to the USA with ESTA approved before trip to Iran
What is the purpose of the Evocation wizard's Potent Cantrip feature?
Why didn't Theresa May consult with Parliament before negotiating a deal with the EU?
If the heap is initialized for security, then why is the stack uninitialized?
What makes a siege story/plot interesting?
Apart from "berlinern", do any other German dialects have a corresponding verb?
Customer Requests (Sometimes) Drive Me Bonkers!
How do I construct this japanese bowl?
The King's new dress
Solution of this Diophantine Equation
When airplanes disconnect from a tanker during air to air refueling, why do they bank so sharply to the right?
Is the concept of a "numerable" fiber bundle really useful or an empty generalization?
Why do professional authors make "consistency" mistakes? And how to avoid them?
Touching Polygon Intersection
The Next CEO of Stack OverflowCalculating intersection between points and polygon in shapefileDetermine if point is within an irregular polygon using PythonAdding a Property (attribute) to a geometry in Shapely/FionaIntersection postgis TopologyExceptionLine vs. Polygon Intersection CoordinatesMulti-part Geometrical Intersection errorPolygon intersection with LineString in PostGisEditing touching polygons simultaneously?Geopandas Line Polygon IntersectionHow is intersection implemented in Shapely?
I have two touching Polygons (two neighbouring counties) and there intersection is a MultiLineString, where:
len(list(counties["King's Co."].intersection(counties["Queen's Co."]).geoms))
list(counties["King's Co."].intersection(counties["Queen's Co."]).geoms )[0]
list(counties["King's Co."].intersection(counties["Queen's Co."]).geoms )[1]
gives:
924
LINESTRING (-7.1060859 53.1680546, -7.1061766 53.1678101)
LINESTRING (-7.1061766 53.1678101, -7.1072205 53.1666679)
I would have expected instead:
len(list(counties["King's Co."].intersection(counties["Queen's Co."]).geoms))
len(list(counties["King's Co."].intersection(counties["Queen's Co."]).geoms)[0])
to give:
1
924
why is intersection giving the result like this? I wish to simplify the borders between a set of bordering counties by reducing the number of points, while at the same time preserving topology.
topology shapely matplotlib
bumped to the homepage by Community♦ yesterday
This question has answers that may be good or bad; the system has marked it active so that they can be reviewed.
add a comment |
I have two touching Polygons (two neighbouring counties) and there intersection is a MultiLineString, where:
len(list(counties["King's Co."].intersection(counties["Queen's Co."]).geoms))
list(counties["King's Co."].intersection(counties["Queen's Co."]).geoms )[0]
list(counties["King's Co."].intersection(counties["Queen's Co."]).geoms )[1]
gives:
924
LINESTRING (-7.1060859 53.1680546, -7.1061766 53.1678101)
LINESTRING (-7.1061766 53.1678101, -7.1072205 53.1666679)
I would have expected instead:
len(list(counties["King's Co."].intersection(counties["Queen's Co."]).geoms))
len(list(counties["King's Co."].intersection(counties["Queen's Co."]).geoms)[0])
to give:
1
924
why is intersection giving the result like this? I wish to simplify the borders between a set of bordering counties by reducing the number of points, while at the same time preserving topology.
topology shapely matplotlib
bumped to the homepage by Community♦ yesterday
This question has answers that may be good or bad; the system has marked it active so that they can be reviewed.
Without knowing your data, this is guesswork.
– inc42
Aug 4 '18 at 9:11
add a comment |
I have two touching Polygons (two neighbouring counties) and there intersection is a MultiLineString, where:
len(list(counties["King's Co."].intersection(counties["Queen's Co."]).geoms))
list(counties["King's Co."].intersection(counties["Queen's Co."]).geoms )[0]
list(counties["King's Co."].intersection(counties["Queen's Co."]).geoms )[1]
gives:
924
LINESTRING (-7.1060859 53.1680546, -7.1061766 53.1678101)
LINESTRING (-7.1061766 53.1678101, -7.1072205 53.1666679)
I would have expected instead:
len(list(counties["King's Co."].intersection(counties["Queen's Co."]).geoms))
len(list(counties["King's Co."].intersection(counties["Queen's Co."]).geoms)[0])
to give:
1
924
why is intersection giving the result like this? I wish to simplify the borders between a set of bordering counties by reducing the number of points, while at the same time preserving topology.
topology shapely matplotlib
I have two touching Polygons (two neighbouring counties) and there intersection is a MultiLineString, where:
len(list(counties["King's Co."].intersection(counties["Queen's Co."]).geoms))
list(counties["King's Co."].intersection(counties["Queen's Co."]).geoms )[0]
list(counties["King's Co."].intersection(counties["Queen's Co."]).geoms )[1]
gives:
924
LINESTRING (-7.1060859 53.1680546, -7.1061766 53.1678101)
LINESTRING (-7.1061766 53.1678101, -7.1072205 53.1666679)
I would have expected instead:
len(list(counties["King's Co."].intersection(counties["Queen's Co."]).geoms))
len(list(counties["King's Co."].intersection(counties["Queen's Co."]).geoms)[0])
to give:
1
924
why is intersection giving the result like this? I wish to simplify the borders between a set of bordering counties by reducing the number of points, while at the same time preserving topology.
topology shapely matplotlib
topology shapely matplotlib
asked Mar 17 '18 at 20:51
BazBaz
1022
1022
bumped to the homepage by Community♦ yesterday
This question has answers that may be good or bad; the system has marked it active so that they can be reviewed.
bumped to the homepage by Community♦ yesterday
This question has answers that may be good or bad; the system has marked it active so that they can be reviewed.
Without knowing your data, this is guesswork.
– inc42
Aug 4 '18 at 9:11
add a comment |
Without knowing your data, this is guesswork.
– inc42
Aug 4 '18 at 9:11
Without knowing your data, this is guesswork.
– inc42
Aug 4 '18 at 9:11
Without knowing your data, this is guesswork.
– inc42
Aug 4 '18 at 9:11
add a comment |
1 Answer
1
active
oldest
votes
In your case;
mls = a.intersection(b)
returnsMultiLineString
.len(mls)
gives you lines' count inmls
.gs = mls.geoms
returnsGeometrySequence
andlen(gs)
equalslen(mls)
.
So, len(mls) = len(mls.geoms) = len(list(mls)) = len(list(mls.geoms))
.
More specifically, it's all about __len__()
special class method definition.
Note that len(list(mls))
equals len(list(mls.geoms))
, however, list(mls)==list(mls.geoms)
returns False
(but, of course, lines in the lists have similar geometry)
But why do I get these 2 point LINESTRINGs? I would like a LINESTRING containing all the boarding points, presuming the counties have a single border with each other.
– Baz
Mar 18 '18 at 8:46
There is probably a small gap between them or the points that you think are overlapping in the two counties' polygon are not actually overlapping.
– Kadir Şahbaz
Mar 18 '18 at 18:51
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%2f276186%2ftouching-polygon-intersection%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
In your case;
mls = a.intersection(b)
returnsMultiLineString
.len(mls)
gives you lines' count inmls
.gs = mls.geoms
returnsGeometrySequence
andlen(gs)
equalslen(mls)
.
So, len(mls) = len(mls.geoms) = len(list(mls)) = len(list(mls.geoms))
.
More specifically, it's all about __len__()
special class method definition.
Note that len(list(mls))
equals len(list(mls.geoms))
, however, list(mls)==list(mls.geoms)
returns False
(but, of course, lines in the lists have similar geometry)
But why do I get these 2 point LINESTRINGs? I would like a LINESTRING containing all the boarding points, presuming the counties have a single border with each other.
– Baz
Mar 18 '18 at 8:46
There is probably a small gap between them or the points that you think are overlapping in the two counties' polygon are not actually overlapping.
– Kadir Şahbaz
Mar 18 '18 at 18:51
add a comment |
In your case;
mls = a.intersection(b)
returnsMultiLineString
.len(mls)
gives you lines' count inmls
.gs = mls.geoms
returnsGeometrySequence
andlen(gs)
equalslen(mls)
.
So, len(mls) = len(mls.geoms) = len(list(mls)) = len(list(mls.geoms))
.
More specifically, it's all about __len__()
special class method definition.
Note that len(list(mls))
equals len(list(mls.geoms))
, however, list(mls)==list(mls.geoms)
returns False
(but, of course, lines in the lists have similar geometry)
But why do I get these 2 point LINESTRINGs? I would like a LINESTRING containing all the boarding points, presuming the counties have a single border with each other.
– Baz
Mar 18 '18 at 8:46
There is probably a small gap between them or the points that you think are overlapping in the two counties' polygon are not actually overlapping.
– Kadir Şahbaz
Mar 18 '18 at 18:51
add a comment |
In your case;
mls = a.intersection(b)
returnsMultiLineString
.len(mls)
gives you lines' count inmls
.gs = mls.geoms
returnsGeometrySequence
andlen(gs)
equalslen(mls)
.
So, len(mls) = len(mls.geoms) = len(list(mls)) = len(list(mls.geoms))
.
More specifically, it's all about __len__()
special class method definition.
Note that len(list(mls))
equals len(list(mls.geoms))
, however, list(mls)==list(mls.geoms)
returns False
(but, of course, lines in the lists have similar geometry)
In your case;
mls = a.intersection(b)
returnsMultiLineString
.len(mls)
gives you lines' count inmls
.gs = mls.geoms
returnsGeometrySequence
andlen(gs)
equalslen(mls)
.
So, len(mls) = len(mls.geoms) = len(list(mls)) = len(list(mls.geoms))
.
More specifically, it's all about __len__()
special class method definition.
Note that len(list(mls))
equals len(list(mls.geoms))
, however, list(mls)==list(mls.geoms)
returns False
(but, of course, lines in the lists have similar geometry)
answered Mar 18 '18 at 0:33
Kadir ŞahbazKadir Şahbaz
4,48721231
4,48721231
But why do I get these 2 point LINESTRINGs? I would like a LINESTRING containing all the boarding points, presuming the counties have a single border with each other.
– Baz
Mar 18 '18 at 8:46
There is probably a small gap between them or the points that you think are overlapping in the two counties' polygon are not actually overlapping.
– Kadir Şahbaz
Mar 18 '18 at 18:51
add a comment |
But why do I get these 2 point LINESTRINGs? I would like a LINESTRING containing all the boarding points, presuming the counties have a single border with each other.
– Baz
Mar 18 '18 at 8:46
There is probably a small gap between them or the points that you think are overlapping in the two counties' polygon are not actually overlapping.
– Kadir Şahbaz
Mar 18 '18 at 18:51
But why do I get these 2 point LINESTRINGs? I would like a LINESTRING containing all the boarding points, presuming the counties have a single border with each other.
– Baz
Mar 18 '18 at 8:46
But why do I get these 2 point LINESTRINGs? I would like a LINESTRING containing all the boarding points, presuming the counties have a single border with each other.
– Baz
Mar 18 '18 at 8:46
There is probably a small gap between them or the points that you think are overlapping in the two counties' polygon are not actually overlapping.
– Kadir Şahbaz
Mar 18 '18 at 18:51
There is probably a small gap between them or the points that you think are overlapping in the two counties' polygon are not actually overlapping.
– Kadir Şahbaz
Mar 18 '18 at 18:51
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%2f276186%2ftouching-polygon-intersection%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
Without knowing your data, this is guesswork.
– inc42
Aug 4 '18 at 9:11