python - Create new shapely polygon by subtracting the intersection with another polygon -


i have 2 shapely multipolygon instances (made of lon,lat points) intersect @ various parts. i'm trying loop through, determine if there's intersection between 2 polygons, , create new polygon excludes intersection. attached image, don't want red circle overlap yellow contour, want edge yellow contour starts.

i've tried following instructions here doesn't change output @ all, plus don't want merge them 1 cascading union. i'm not getting error messages, when add these multipolygons kml file (just raw text manipulation in python, no fancy program) they're still showing circles without modifications. 

# multipol1 , multipol2 shapely multipolygons shapely.ops import cascaded_union itertools import combinations shapely.geometry import polygon,multipolygon  outmulti = [] pol in multipoly1:     pol2 in multipoly2:         if pol.intersects(pol2)==true:             # if intersect, create new polygon             # pol minus intersection             intersection = pol.intersection(pol2)             nonoverlap = pol.difference(intersection)             outmulti.append(nonoverlap)          else:             # otherwise, keep initial polygon is.             outmulti.append(pol)  finalpol = multipolygon(outmulti)

polygon overlap

i guess can use symmetric_difference between theses 2 polygons, combined difference second polygon achieve want (the symmetric difference brings non-overlapping parts 2 polygons, on removed parts of polygon 2 difference). haven't tested might :

# multipol1 , multipol2 shapely multipolygons shapely.ops import cascaded_union itertools import combinations shapely.geometry import polygon,multipolygon  outmulti = [] pol in multipoly1:     pol2 in multipoly2:         if pol.intersects(pol2)==true:             # if intersect, create new polygon             # pol minus intersection             nonoverlap = (pol.symmetric_difference(pol2)).difference(pol2)             outmulti.append(nonoverlap)          else:             # otherwise, keep initial polygon is.             outmulti.append(pol)  finalpol = multipolygon(outmulti)

-

Comments

Post a Comment

Popular posts from this blog

javascript - Thinglink image not visible until browser resize -

Image
i trying repurpose pre-existing theme, includes drop down menu displays image. trying embed thinglink drop down menu's content section, not visible until resize browser window. here menu looks like: then when clicked, can see empty space image should show: then when browser resized, image appears: this code content of menu: <div class="ajax_accordion_content" style="display: none;"> <div class="report-detail"> <table border="0" cellpadding="0" cellspacing="0" width="100%" class="table_basic"> <tbody> <!-- thinglink image embedded --> <img style="max-width: 100%;" src="http://cdn.thinglink.me/api/image/832705713204625409/1024/10/scaletowidth#tl-832705713204625409;1043138249'" class="alwaysthinglink"/> <script async charset="utf-8" src="http://cdn.thinglink.me/jse/embed.j
Read more

firebird - Error "invalid transaction handle (expecting explicit transaction start)" executing script from Delphi -

i'm developing update app delphi 10. i'm running windows 7 64x, firebird 2.5.1.26351 32x. execution order: checks current version(select statement firebird database) downloads update(via ftp) apply it(.exe files, , execute .sql) error code -901 pops out when try execute .sql files. know there incompatibility issue w7 64x , fdb 32x, but, can everything: backup , restore(gbak), db commands(all). tried fdb 64x instead, , not working. code: //components fdwaitcursor: tfdguixwaitcursor; fddriverlink: tfdphysfbdriverlink; fdconnection:tfdconnection; fdscript: tfdscript; function tthr_script.createcomponents:boolean; begin try fddriverlink := tfdphysfbdriverlink.create(fddriverlink); fdwaitcursor := tfdguixwaitcursor.create(fdwaitcursor); fdconnection := tfdconnection.create(fdconnection); fdconnection begin drivername := 'fb'; params.username := thrbanco.bduser; params.password := thrbanco.bdpass; params.database :=
Read more

mongodb - How to keep track of users making Stripe Payments -

i'm trying figure out way mean stack application keep track of users have paid grant them access portion of webpage. i've considered several options: stripe customer id, mongodb record, , html attribute can update. my mean stack keeps track of users jwt, , appears stripe assigns them own customer id isn't ideal. can done jwt opposed forced cutomer id? mongodb record. i'm thinking might best option. when new user has been created, i'll give attribute of haspaid = no. update record of customer when payment submitted. guess run script set unpaid each day? html element/attribute. don't know if possible; cool create key carried during html session after payment received. if person closers browser session closed? i'm looking guidance on 3 options determine if they're best solution. also, if has suggestions alternatives, i'm ears! thanks in advance. speaking generally, common approach second one: use attribute in data model indicates
Read more