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?










0















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.










share|improve this question














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















0















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.










share|improve this question














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













0












0








0








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.










share|improve this question














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






share|improve this question













share|improve this question











share|improve this question




share|improve this question










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

















  • 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










1 Answer
1






active

oldest

votes


















0














In your case;




  • mls = a.intersection(b) returns MultiLineString.


  • len(mls)gives you lines' count in mls.


  • gs = mls.geomsreturns GeometrySequenceand len(gs) equals len(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)






share|improve this answer























  • 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












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
);



);













draft saved

draft discarded


















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









0














In your case;




  • mls = a.intersection(b) returns MultiLineString.


  • len(mls)gives you lines' count in mls.


  • gs = mls.geomsreturns GeometrySequenceand len(gs) equals len(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)






share|improve this answer























  • 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
















0














In your case;




  • mls = a.intersection(b) returns MultiLineString.


  • len(mls)gives you lines' count in mls.


  • gs = mls.geomsreturns GeometrySequenceand len(gs) equals len(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)






share|improve this answer























  • 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














0












0








0







In your case;




  • mls = a.intersection(b) returns MultiLineString.


  • len(mls)gives you lines' count in mls.


  • gs = mls.geomsreturns GeometrySequenceand len(gs) equals len(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)






share|improve this answer













In your case;




  • mls = a.intersection(b) returns MultiLineString.


  • len(mls)gives you lines' count in mls.


  • gs = mls.geomsreturns GeometrySequenceand len(gs) equals len(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)







share|improve this answer












share|improve this answer



share|improve this answer










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


















  • 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


















draft saved

draft discarded
















































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.




draft saved


draft discarded














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





















































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







Popular posts from this blog

រឿង រ៉ូមេអូ និង ហ្ស៊ុយលីយេ សង្ខេបរឿង តួអង្គ បញ្ជីណែនាំ

QGIS export composer to PDF scale the map [closed] Planned maintenance scheduled April 23, 2019 at 23:30 UTC (7:30pm US/Eastern) Announcing the arrival of Valued Associate #679: Cesar Manara Unicorn Meta Zoo #1: Why another podcast?Print Composer QGIS 2.6, how to export image?QGIS 2.8.1 print composer won't export all OpenCycleMap base layer tilesSave Print/Map QGIS composer view as PNG/PDF using Python (without changing anything in visible layout)?Export QGIS Print Composer PDF with searchable text labelsQGIS Print Composer does not change from landscape to portrait orientation?How can I avoid map size and scale changes in print composer?Fuzzy PDF export in QGIS running on macSierra OSExport the legend into its 100% size using Print ComposerScale-dependent rendering in QGIS PDF output

PDF-ში გადმოწერა სანავიგაციო მენიუproject page