Re: [Jts-topo-suite-user] Cut a polygon

Mon, 24 Mar 2014 04:59:34 -0700

In this case:
GEOMETRYCOLLECTION (POLYGON ((-189.34023087182823 -23.999596895664688, -189.7083100115792 -23.37233753317051, -189.81454603686237 -23.14082519789132, -199.50078909208844 -0.9798720034454546, -199.52178748107298 -0.9281806607285326, -211.22891653195984 27.400444524136535, -211.2997062270637 27.637290504825746, -213.71464580896264 35.76662932955155, -213.74479704612884 35.89522898299576, -215.2494697722787 42.22377265712712, -215.2767714671864 42.38967389312838, -216.47624293437477 49.60604976697388, -216.49214196462404 49.76963078682333, -216.83996462398034 53.458958645022264, -216.834999998697 53.757038455376595, -216.74128891653132 57.2953416235066, -216.70139296440243 57.589841196217215, -216.48379144706522 59.03487827929562, -216.4245391551226 59.280628758692266, -216.0524061392402 60.72445715131305, -215.96478648413984 60.97201636936241, -215.3515768533328 62.61308233211142, -215.23195700985696 62.861477357315756, -214.40998238046546 64.47385732834617, -214.2309728139015 64.7490111940272, -213.26661078576225 66.15523681395582, -213.08978567776563 66.36985557287719, -212.04434264211534 67.58687680621344, -211.9214774303339 67.71516372314322, -210.84919951552567 68.76249440345262, -210.5156749961542 69.01578286587261, -209.3072938831667 69.841950911113, -208.886443626351 70.0521306273389, -207.16958699858787 70.83218051705164, -206.95131583948574 70.91426057585146, -204.96590836642133 71.5879484311678, -204.64447828496196 71.66311854399945, -199.55858667889774 72.67321434593174, -199.3136447490366 72.7037412065061, -197.7445113520609 72.84644463661724, -197.51892653930102 72.85289089884859, -195.96952680096717 72.84671965984765, -195.7730466723533 72.83575663015081, -187.368412775841 72.10793599697728, -171.6401959525321 70.7664727153026, -171.42045119278444 70.7309914219543, -150.7728075365562 66.67364152959651, -150.7486290196615 66.66908398645607, -130.04471024865344 62.33072284949449, -130.03379894011644 62.32856518364399, -119.56955718999943 60.076174596068356, -119.55323279602226 60.07275249532549, -88.18913269054197 53.1256560251281, -55.61338730070978 46.956122970908645, -38.53201859550566 44.41037263272656, -23.339358952082787 43.29854039674695, -13.776360057861238 43.16130176334489, -4.888544850475359 43.505088861059164, 1.6297275014365256 44.179089494758, 7.63151865324651 45.13110079150486, 21.840670598718134 48.15470259714328, 27.778596817051202 49.508941348324186, 40.3754318263261 52.666291998520634, 40.91499362221456 52.68175112520325, 41.37046774472211 52.64853793356068, 42.38639249952542 52.38778273790826, 43.154260957485896 52.00213350029228, 43.8323364299301 51.47362592938373, 44.39426837850144 50.82250181425273, 44.818173173621325 50.07383707444774, 45.08752220581882 49.25656589366455, 45.12976093971702 48.91060507333589, 54.65814890704125 8.505516474622436, 63.28811294756317 -21.263547202396825, 63.30736615021132 -21.97188039539365, 63.26563339827509 -22.590883235934744, 62.92619878049044 -23.804589763013006, 62.429530264282604 -24.706225236754516, 61.75540065164076 -25.483987820655052, 60.93377932097545 -26.103483044649387, 60.001175040681986 -26.537430245971848, 58.9989620415941 -26.76687127437206, 58.66518257707229 -26.802854683839918, 26.272829786474368 -29.188383892125945, 26.24545582590859 -29.19069752576575, -6.3782356797584825 -31.24121064748633, -6.431708198875868 -31.24424880615893, -40.00673608081818 -32.47166165859356, -40.056036175676425 -32.472857457600405, -63.80038500832181 -32.678633587902326, -63.856770399612905 -32.67773329923304, -83.67546883819149 -32.1344330942304, -83.72015005994969 -32.13215520105822, -107.64182828773635 -30.868230868126098, -107.67150734173136 -30.866070312447455, -151.25771802234655 -28.121987705732177, -151.28052358250247 -28.12000020578541, -187.0037511260975 -25.56041964619766, -187.79748319078288 -25.2891786124156, -188.39282503436678 -24.940667998448973, -188.91422549801334 -24.488446922002574, -189.3437475829376 -23.947852499668663, -189.34023087182823 -23.999596895664688), (-183.74366557756184 -18.739878772142347, -150.7485654092785 -22.495365419428254, -107.36276669239462 -26.958861406537697, -83.52896053826457 -28.495125211394782, -63.804977597336766 -28.74845753342111, -40.19986198282159 -27.7507338132556, -6.832584608732528 -24.975533388689712, 25.520973085964023 -21.224887782716237, 54.189017693307484 -18.987834732311722, 45.93443926055564 6.024052613301043, 45.906274025073145 6.134953588074435, 38.339395706022046 44.442287803762405, 29.557961538294375 41.678366149750715, 29.498067331344803 41.66615370114021, 23.327138964485602 40.641671253655026, 23.296526860947093 40.636096342986825, 8.870405415466209 38.454212209143954, 8.771150552715861 38.439570958135796, 2.474393284635687 37.81357179841203, 2.375029250441676 37.80425642152594, -4.415658528724371 37.49409572147816, -4.527067285868672 37.49024817463117, -13.699157983132855 37.64796728242792, -13.78401092957036 37.64997950141432, -23.530456608543005 38.317507214350044, -23.614551351414224 38.32411124844026, -38.93234890919469 40.21406117902395, -39.02149122306413 40.22674153304924, -56.199038460920875 43.47161280620974, -56.2254772636307 43.47658810047885, -88.81134449565745 49.90372956652703, -88.82131521165785 49.90566039731954, -120.42457063543031 55.7430359088598, -130.95048832592494 57.48550654535399, -151.81966063939345 60.738653716759636, -172.5018396416465 63.63817152116468, -187.8335876763318 64.10452925620902, -196.1733387082745 64.36386586838687, -197.20029962490466 64.30370712571107, -198.1727281022146 64.1541190723478, -202.54424003533177 63.05515916697248, -203.73952351461017 62.56813188272883, -204.77965786322142 62.04511040510531, -204.98588859091916 61.869631160859676, -205.3398834965112 61.47116141963989, -206.03570047576105 60.626742496430495, -206.56238112446619 59.82511514821481, -206.9924953819437 58.9396506320999, -207.3501928159242 57.95968618165113, -207.55239805990203 57.17050741081968, -207.6761081958201 56.39658251512355, -207.78720889399727 53.63027251397561, -207.5558703731397 50.63975008069613, -206.58544242199187 43.82033469308253, -205.32045679401867 37.864413120444254, -203.2169543685503 30.059368823119506, -192.33143835472026 1.9827235332980355, -183.74366557756184 -18.739878772142347)), MULTILINESTRING ((-189.34023087182823 -23.999596895664688, -183.74366557756184 -18.739878772142347), (-191.2327182736261 1.9435105187995658, -200.23845107110347 -1.7886069931132687), (-200.09695434156373 24.778393721904354, -209.69004278199 21.05908478159426), (-205.9487409451637 46.44411032624282, -216.7095967435767 44.912800873526606), (-204.2720587167773 61.16931201162656, -212.145583032318 68.89381339475287), (-188.07167685810737 63.11144960580831, -187.75969930504144 73.14402887903553), (-164.37700832231837 61.48933939794763, -163.16611922251747 70.12657731244234), (-140.08236556499114 57.896933267656394, -138.94155179843858 65.21533536179686), (-115.89501749009426 53.889475540187135, -114.75995484472261 60.03464722440481), (-91.61852424623068 49.405411943966456, -90.64194761357712 54.69254838081611), (-67.42890548012129 44.66705397305445, -66.36969068153103 50.037315850654366), (-43.14317265533725 39.98764940408276, -42.00747077457889 45.99977471278514), (-17.99936703910044 36.93634384286121, -17.500966409299444 44.213397020710175), (7.793017300066136 37.33739911630528, 6.92879091797696 46.030419974492915), (33.78991162455544 41.96199264882447, 30.84844797872221 51.3075069564234), (38.411614128667026 38.916359904696236, 48.063745756099 40.82295317305089), (43.195431867293024 14.698274557783561, 53.73361646504761 16.77989057694918), (50.244230118862134 -10.225693887691547, 60.32637521706078 -6.898321733264217), (53.89897245750319 -18.03082179551973, 56.893603213707344 -27.911695606997462), (31.265778658173964 -19.773563511976143, 32.045514706073575 -29.765954080942624), (6.396971156154647 -22.435175745040425, 7.433185986004305 -31.37369988292409), (-17.911515323277328 -24.893529880813613, -17.26723487875295 -32.63995572769953), (-42.276369046538974 -26.837609229065336, -41.99522714673515 -33.48913701659018), (-66.67535223385984 -27.711508328442026, -66.75096609544008 -33.59866281381102), (-91.11920114646207 -27.003803252106476, -91.48784982868658 -32.723069177226556), (-115.54919334790411 -25.111370298597222, -116.18906568775478 -31.331017996598366), (-140.1027435788571 -22.585320850013503, -140.8428109081717 -29.778877712014648), (-164.56723087376625 -19.916074873881627, -165.49905434912606 -28.10292568611258)))

the method does not produce the correct output (see the missing cut in the lower right part of the shape).

In other cases, this for example:
GEOMETRYCOLLECTION (POLYGON ((-189.34023087182823 -23.999596895664688, -189.7083100115792 -23.37233753317051, -189.81454603686237 -23.14082519789132, -199.50078909208844 -0.9798720034454546, -199.52178748107298 -0.9281806607285326, -211.22891653195984 27.400444524136535, -211.2997062270637 27.637290504825746, -213.71464580896264 35.76662932955155, -213.74479704612884 35.89522898299576, -215.2494697722787 42.22377265712712, -215.2767714671864 42.38967389312838, -216.47624293437477 49.60604976697388, -216.49214196462404 49.76963078682333, -216.83996462398034 53.458958645022264, -216.834999998697 53.757038455376595, -216.74128891653132 57.2953416235066, -216.70139296440243 57.589841196217215, -216.48379144706522 59.03487827929562, -216.4245391551226 59.280628758692266, -216.0524061392402 60.72445715131305, -215.96478648413984 60.97201636936241, -215.3515768533328 62.61308233211142, -215.23195700985696 62.861477357315756, -214.40998238046546 64.47385732834617, -214.2309728139015 64.7490111940272, -213.26661078576225 66.15523681395582, -213.08978567776563 66.36985557287719, -212.04434264211534 67.58687680621344, -211.9214774303339 67.71516372314322, -210.84919951552567 68.76249440345262, -210.5156749961542 69.01578286587261, -209.3072938831667 69.841950911113, -208.886443626351 70.0521306273389, -207.16958699858787 70.83218051705164, -206.95131583948574 70.91426057585146, -204.96590836642133 71.5879484311678, -204.64447828496196 71.66311854399945, -199.55858667889774 72.67321434593174, -199.3136447490366 72.7037412065061, -197.7445113520609 72.84644463661724, -197.51892653930102 72.85289089884859, -195.96952680096717 72.84671965984765, -195.7730466723533 72.83575663015081, -187.368412775841 72.10793599697728, -171.6401959525321 70.7664727153026, -171.42045119278444 70.7309914219543, -150.7728075365562 66.67364152959651, -150.7486290196615 66.66908398645607, -130.04471024865344 62.33072284949449, -130.03379894011644 62.32856518364399, -119.56955718999943 60.076174596068356, -119.55323279602226 60.07275249532549, -88.18913269054197 53.1256560251281, -55.61338730070978 46.956122970908645, -38.53201859550566 44.41037263272656, -23.339358952082787 43.29854039674695, -13.776360057861238 43.16130176334489, -4.888544850475359 43.505088861059164, 1.6297275014365256 44.179089494758, 7.63151865324651 45.13110079150486, 21.840670598718134 48.15470259714328, 27.778596817051202 49.508941348324186, 40.3754318263261 52.666291998520634, 40.91499362221456 52.68175112520325, 41.37046774472211 52.64853793356068, 42.38639249952542 52.38778273790826, 43.154260957485896 52.00213350029228, 43.8323364299301 51.47362592938373, 44.39426837850144 50.82250181425273, 44.818173173621325 50.07383707444774, 45.08752220581882 49.25656589366455, 45.12976093971702 48.91060507333589, 54.65814890704125 8.505516474622436, 63.28811294756317 -21.263547202396825, 63.30736615021132 -21.97188039539365, 63.26563339827509 -22.590883235934744, 62.92619878049044 -23.804589763013006, 62.429530264282604 -24.706225236754516, 61.75540065164076 -25.483987820655052, 60.93377932097545 -26.103483044649387, 60.001175040681986 -26.537430245971848, 58.9989620415941 -26.76687127437206, 58.66518257707229 -26.802854683839918, 26.272829786474368 -29.188383892125945, 26.24545582590859 -29.19069752576575, -6.3782356797584825 -31.24121064748633, -6.431708198875868 -31.24424880615893, -40.00673608081818 -32.47166165859356, -40.056036175676425 -32.472857457600405, -63.80038500832181 -32.678633587902326, -63.856770399612905 -32.67773329923304, -83.67546883819149 -32.1344330942304, -83.72015005994969 -32.13215520105822, -107.64182828773635 -30.868230868126098, -107.67150734173136 -30.866070312447455, -151.25771802234655 -28.121987705732177, -151.28052358250247 -28.12000020578541, -187.0037511260975 -25.56041964619766, -187.79748319078288 -25.2891786124156, -188.39282503436678 -24.940667998448973, -188.91422549801334 -24.488446922002574, -189.3437475829376 -23.947852499668663, -189.34023087182823 -23.999596895664688), (-183.74366557756184 -18.739878772142347, -150.7485654092785 -22.495365419428254, -107.36276669239462 -26.958861406537697, -83.52896053826457 -28.495125211394782, -63.804977597336766 -28.74845753342111, -40.19986198282159 -27.7507338132556, -6.832584608732528 -24.975533388689712, 25.520973085964023 -21.224887782716237, 54.189017693307484 -18.987834732311722, 45.93443926055564 6.024052613301043, 45.906274025073145 6.134953588074435, 38.339395706022046 44.442287803762405, 29.557961538294375 41.678366149750715, 29.498067331344803 41.66615370114021, 23.327138964485602 40.641671253655026, 23.296526860947093 40.636096342986825, 8.870405415466209 38.454212209143954, 8.771150552715861 38.439570958135796, 2.474393284635687 37.81357179841203, 2.375029250441676 37.80425642152594, -4.415658528724371 37.49409572147816, -4.527067285868672 37.49024817463117, -13.699157983132855 37.64796728242792, -13.78401092957036 37.64997950141432, -23.530456608543005 38.317507214350044, -23.614551351414224 38.32411124844026, -38.93234890919469 40.21406117902395, -39.02149122306413 40.22674153304924, -56.199038460920875 43.47161280620974, -56.2254772636307 43.47658810047885, -88.81134449565745 49.90372956652703, -88.82131521165785 49.90566039731954, -120.42457063543031 55.7430359088598, -130.95048832592494 57.48550654535399, -151.81966063939345 60.738653716759636, -172.5018396416465 63.63817152116468, -187.8335876763318 64.10452925620902, -196.1733387082745 64.36386586838687, -197.20029962490466 64.30370712571107, -198.1727281022146 64.1541190723478, -202.54424003533177 63.05515916697248, -203.73952351461017 62.56813188272883, -204.77965786322142 62.04511040510531, -204.98588859091916 61.869631160859676, -205.3398834965112 61.47116141963989, -206.03570047576105 60.626742496430495, -206.56238112446619 59.82511514821481, -206.9924953819437 58.9396506320999, -207.3501928159242 57.95968618165113, -207.55239805990203 57.17050741081968, -207.6761081958201 56.39658251512355, -207.78720889399727 53.63027251397561, -207.5558703731397 50.63975008069613, -206.58544242199187 43.82033469308253, -205.32045679401867 37.864413120444254, -203.2169543685503 30.059368823119506, -192.33143835472026 1.9827235332980355, -183.74366557756184 -18.739878772142347)), MULTILINESTRING ((-189.34023087182823 -23.999596895664688, -183.74366557756184 -18.739878772142347), (-188.5657948589483 -4.491870893297033, -197.44800929687818 -8.172800451594462), (-195.06090859316748 11.789090324491482, -204.36790345256492 8.180701909963458), (-201.47042224208468 28.320933277299698, -211.14153625450135 24.57137326516025), (-205.594799278337 43.96941378071894, -216.1586361118554 41.725747270393256), (-206.5646956545791 56.99643675249693, -217.33470840088648 59.755954792645184), (-202.2949986505003 62.08671777376002, -205.0023189409107 72.60616275152768), (-187.4348102577381 63.09193681881397, -187.13069769885806 73.08979206289612), (-170.64101809956622 62.367516106219945, -169.38200443032684 71.34802522188485), (-153.10392146155024 59.908920039820956, -151.9796402168946 67.92838655784573), (-135.7081424329279 57.21506666578907, -134.60276034400096 64.30617200796495), (-118.38567054455416 54.34951919559204, -117.23434798603525 60.582721226212044), (-100.95639334383156 51.13019148548487, -99.91885561495641 56.74736821030223), (-83.58390792639308 47.853419564223806, -82.50951709550704 53.30062453948416), (-66.18396507853927 44.42150555720924, -65.112051989888 49.85614827910643), (-48.7400227789363 41.04490501425667, -47.62633572527875 46.9404892645863), (-30.899477092888002 38.21536155728028, -30.058123932329153 45.03441913096908), (-12.583209162751828 36.62863010064297, -12.452814399787528 44.21168461486262), (5.942115099995627 37.15338960591037, 5.088613885060564 45.738528890808304), (23.983165500167214 39.73689593179683, 22.39470738592889 49.3049087960316), (38.162934580730735 43.45798019368265, 39.964513039833164 53.507260931991595), (38.78624563662488 37.01978726562838, 48.50776606156939 38.94008695856811), (42.22078144897174 19.632443939639796, 52.57844246799403 21.678400992136115), (46.490533674430594 1.1482403916450807, 57.22541604227137 4.691032954528168), (52.041272758495914 -15.670845136051366, 61.90381212962846 -12.415948627835778), (53.7693764629456 -18.080144676301213, 58.250101466465054 -27.772017722009622), (39.43846859440403 -19.135824133435808, 40.22098715232908 -29.163872829291744), (21.440218799341935 -20.691259629078317, 22.568319562280866 -30.4224037982441), (3.8474373731058846 -22.73073513207625, 4.868079349529567 -31.534925807212588), (-13.53415080085013 -24.529458705784123, -12.872950985111197 -32.47931250955336), (-31.18924001831985 -25.99785573868433, -30.596280580270744 -33.12722984292165), (-48.6376604321453 -27.106483629426304, -48.365549658596535 -33.54434427135658), (-66.07863874862126 -27.719172440363664, -66.15436423427695 -33.61501776842519), (-83.57115859722508 -27.49033001543272, -83.93417234636664 -33.122174820322115), (-101.26841218685675 -26.349612080608992, -101.64463768448529 -32.18642575308495), (-118.73343666790491 -24.783777913504796, -119.38630286090094 -31.129727662353368), (-136.36162657730705 -22.97020395041905, -137.08642763453364 -30.015370534959935), (-153.9141805247203 -21.128600010482476, -154.79528638661813 -28.86985404827525), (-171.52018005435815 -19.12469349683492, -172.48510550308785 -27.602372909907114)))
the method (also with testbuilder) produce a topology exception on UnaryUnion. But if I add an MCIndexNoding operation before MergeLines/UnaryUnion/Polygonize operations, the exception disappear. Unfortunately i was not able to translate the MCIndexNoding operation in java code. Can you give me some advice? The sample code in TestBuilder was not helpful.

Ciao

Sandro


Il 20/03/14 19:00, Martin Davis ha scritto:
Currently that approach (extract lines/union/polygonize/keep the useful bits) is the best practice in JTS for cutting polygons with lines.


On Thu, Mar 20, 2014 at 7:04 AM, Sandro Zacchino <[email protected]> wrote:
Currently the best method I've found is the one explained here http://stackoverflow.com/questions/4129241/split-a-polygon-with-a-linestring-in-jts

To that algorithm I added a for-loop to remove the annulus "hole" geometry appearing together with slices in the result.
Here an excerpt of the code:

ArrayList<LineString> slices = new ArrayList<LineString>();

/* ...
    evaluating slices
    ... */


Geometry collection = Utils.createGeometryCollection(new Geometry[]{annulus, Utils.createMultiLineString(slices)});
LineMerger lineMerger = new LineMerger();
lineMerger.add(collection);

UnaryUnionOp op = new UnaryUnionOp(lineMerger.getMergedLineStrings());
Geometry geom1 = op.union();

Polygonizer polygonizer = new Polygonizer();
polygonizer.add(geom1);
Collection<Geometry> resultCollection = polygonizer.getPolygons();

Iterator<Geometry> iterator = resultCollection.iterator();
while(iterator.hasNext())
{
     Geometry g = (Geometry)iterator.next();
     if(g.contains(centroid))
         iterator.remove();
}
Geometry[] result = (Geometry[])resultCollection.toArray(new Geometry[0]);

Bye

Sandro

Nella citazione in data Thu Mar 20 09:31:36 2014, Sandro Zacchino ha scritto:
>
> Hi,
>
> I'm looking for an efficient way to cut an 2D annulus (that is a polygon with an hole, a donut) in slices. In particular I have a collection of segments that represents all the cuts I have to do on the annulus and, as a result, I need a collection of closed polygons, one for each portion of the annulus I obtain from the cutting operation.
>
> Thank you in advance
>
> Sandro
>
> ------------------------------------------------------------------------------
> Learn Graph Databases - Download FREE O'Reilly Book
> "Graph Databases" is the definitive new guide to graph databases and their
> applications. Written by three acclaimed leaders in the field,
> this first edition is now available. Download your free book today!
> http://p.sf.net/sfu/13534_NeoTech
> _______________________________________________
> Jts-topo-suite-user mailing list
> [email protected]
> https://lists.sourceforge.net/lists/listinfo/jts-topo-suite-user

------------------------------------------------------------------------------
Learn Graph Databases - Download FREE O'Reilly Book
"Graph Databases" is the definitive new guide to graph databases and their
applications. Written by three acclaimed leaders in the field,
this first edition is now available. Download your free book today!
http://p.sf.net/sfu/13534_NeoTech
_______________________________________________
Jts-topo-suite-user mailing list
[email protected]
https://lists.sourceforge.net/lists/listinfo/jts-topo-suite-user


------------------------------------------------------------------------------
Learn Graph Databases - Download FREE O'Reilly Book
"Graph Databases" is the definitive new guide to graph databases and their
applications. Written by three acclaimed leaders in the field,
this first edition is now available. Download your free book today!
http://p.sf.net/sfu/13534_NeoTech
_______________________________________________
Jts-topo-suite-user mailing list
[email protected]
https://lists.sourceforge.net/lists/listinfo/jts-topo-suite-user

Reply via email to