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

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

Crop image to path created in TikZ? Announcing the arrival of Valued Associate #679: Cesar Manara Planned maintenance scheduled April 17/18, 2019 at 00:00UTC (8:00pm US/Eastern)Crop an inserted image?TikZ pictures does not appear in posterImage behind and beyond crop marks?Tikz picture as large as possible on A4 PageTransparency vs image compression dilemmaHow to crop background from image automatically?Image does not cropTikzexternal capturing crop marks when externalizing pgfplots?How to include image path that contains a dollar signCrop image with left size given

Romeo and Juliet ContentsCharactersSynopsisSourcesDate and textThemes and motifsCriticism and interpretationLegacyScene by sceneSee alsoNotes and referencesSourcesExternal linksNavigation menu"Consumer Price Index (estimate) 1800–"10.2307/28710160037-3222287101610.1093/res/II.5.31910.2307/45967845967810.2307/2869925286992510.1525/jams.1982.35.3.03a00050"Dada Masilo: South African dancer who breaks the rules"10.1093/res/os-XV.57.1610.2307/28680942868094"Sweet Sorrow: Mann-Korman's Romeo and Juliet Closes Sept. 5 at MN's Ordway"the original10.2307/45957745957710.1017/CCOL0521570476.009"Ram Leela box office collections hit massive Rs 100 crore, pulverises prediction"Archived"Broadway Revival of Romeo and Juliet, Starring Orlando Bloom and Condola Rashad, Will Close Dec. 8"Archived10.1075/jhp.7.1.04hon"Wherefore art thou, Romeo? To make us laugh at Navy Pier"the original10.1093/gmo/9781561592630.article.O006772"Ram-leela Review Roundup: Critics Hail Film as Best Adaptation of Romeo and Juliet"Archived10.2307/31946310047-77293194631"Romeo and Juliet get Twitter treatment""Juliet's Nurse by Lois Leveen""Romeo and Juliet: Orlando Bloom's Broadway Debut Released in Theaters for Valentine's Day"Archived"Romeo and Juliet Has No Balcony"10.1093/gmo/9781561592630.article.O00778110.2307/2867423286742310.1076/enst.82.2.115.959510.1080/00138380601042675"A plague o' both your houses: error in GCSE exam paper forces apology""Juliet of the Five O'Clock Shadow, and Other Wonders"10.2307/33912430027-4321339124310.2307/28487440038-7134284874410.2307/29123140149-661129123144728341M"Weekender Guide: Shakespeare on The Drive""balcony"UK public library membership"romeo"UK public library membership10.1017/CCOL9780521844291"Post-Zionist Critique on Israel and the Palestinians Part III: Popular Culture"10.2307/25379071533-86140377-919X2537907"Capulets and Montagues: UK exam board admit mixing names up in Romeo and Juliet paper"Istoria Novellamente Ritrovata di Due Nobili Amanti2027/mdp.390150822329610820-750X"GCSE exam error: Board accidentally rewrites Shakespeare"10.2307/29176390149-66112917639"Exam board apologises after error in English GCSE paper which confused characters in Shakespeare's Romeo and Juliet""From Mariotto and Ganozza to Romeo and Guilietta: Metamorphoses of a Renaissance Tale"10.2307/37323537323510.2307/2867455286745510.2307/28678912867891"10 Questions for Taylor Swift"10.2307/28680922868092"Haymarket Theatre""The Zeffirelli Way: Revealing Talk by Florentine Director""Michael Smuin: 1938-2007 / Prolific dance director had showy career"The Life and Art of Edwin BoothRomeo and JulietRomeo and JulietRomeo and JulietRomeo and JulietEasy Read Romeo and JulietRomeo and Julieteeecb12003684p(data)4099369-3n8211610759dbe00d-a9e2-41a3-b2c1-977dd692899302814385X313670221313670221