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

Export of reprojected layer from GEE fails 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?Iterate over features and years to export images for each feature-yearEarth engine reproject: why does reprojecting a pixel at 30m scale gives me pixels with an area of ~550?Exporting slope raster for entire country in GEEGEE Landsat SR year composite // mask cloud / shadow out w/ other method than quality pixel"export region contains no valid (un-masked) pixels - Google Earth EgineGoogle Earth Engine, how to distinguish between rivers/streams and ponds/lakes in a water maskPython script tool fails when processing rainfall data from Google Earth EngineGoogle Earth Engine Error: Number of pixels requested from Image.load exceeds the maximum allowedCloudfree images in small area from Sentinel-2Using Image exportToDrive in Google Earth Engine?

Creating closest line along the point''s azimuth using PostgreSQL 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?Drawing line between points at specific distance in PostGIS?How to efficiently find the closest point over the dateline?How to find the nearest point by using PostGIS function?PostGIS nearest point with LATERAL JOIN in PostgreSQL 9.3+Creating a table and inserting selected streets using plpgsql functionsCreating a table that stores Distances and other columnSaving select query results (year wise) from PostgreSQL/PostGIS to text filesWhat is the information behind this geometry?How to give start and end vertex ids dynamically in pgr_dijkstra?Point to Polygon nearest distance DS_distance is not using geography index & knn <-> or <#> does not give result in orderLine to point conversion with start point and end point detection?

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