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

Sound is not coming out while implementing Text-to-speech in Android activity -

i using texttospeak in 1 of app. while adding speech string texttospeech object instance, doesn't produce sound. though not getting error in log still sound not coming out. want device speak value of textview gets set resultmsg.settext() method. public class resultactivity extends activity { @override protected void oncreate(bundle savedinstancestate){ ............ resultmsg.settext(resultmsg_str); texttospeech = new texttospeech(getapplicationcontext(),newtexttospeech.oninitlistener(){ @override public void oninit(int status) { system.out.println("txt speech status = "+status); system.out.println("txt speech error status = "+texttospeech.error); system.out.println("txt speech success status = "+texttospeech.success); if(status == texttospeech.success){ texttospeech.setlanguage(locale.english); ...
Read more