using System;

using System.Collections.Generic;

using NetTopologySuite.IO;

using NetTopologySuite.Geometries;

using NetTopologySuite.Operation.OverlayNG;

​

class Program

{

    static void Main(string[] args)

    {

        // Řetězec WKT polygonu

        string wktParcela = "POLYGON ((16.81280733 49.07273129, 16.81292522 49.07275114, 16.81299744 49.07277767, 16.81305914 49.07281305, 16.81309458 49.07284182, 16.81311965 49.07286496, 16.81313119 49.07292478, 16.81315019 49.07302325, 16.81315893 49.07306905, 16.81316668 49.07312144, 16.81317303 49.07321012, 16.81317639 49.07329721, 16.81317763 49.07346075, 16.81320286 49.07371568, 16.81321103 49.07376085, 16.81321129 49.07376153, 16.81321134 49.07376186, 16.8132674 49.07390993, 16.81335772 49.07405905, 16.8134846 49.07419145, 16.81348811 49.07419511, 16.81348814 49.07419514, 16.81348861 49.07419563, 16.81373393 49.07440992, 16.81386098 49.0745379, 16.81398696 49.07467835, 16.8140437 49.07475395, 16.81407921 49.07480593, 16.81410941 49.07485487, 16.81415209 49.07494396, 16.81418903 49.07503138, 16.8142261 49.07512036, 16.81425679 49.07519422, 16.81431425 49.07540981, 16.81431749 49.07542209, 16.81383771 49.0754737, 16.81380656 49.07534883, 16.81379286 49.0753068, 16.81375044 49.07520471, 16.81371365 49.0751164, 16.8136792 49.07503487, 16.81364498 49.07496344, 16.81362561 49.07493205, 16.81359746 49.07489085, 16.81355115 49.07482913, 16.81344625 49.07470808, 16.81337851 49.07463947, 16.81333668 49.07459711, 16.813097 49.07439139, 16.81307155 49.07436738, 16.81307109 49.07436689, 16.81307079 49.07436659, 16.81307058 49.07436636, 16.81306728 49.07436293, 16.8129404 49.07423053, 16.81289955 49.07417773, 16.81280922 49.07402861, 16.81278888 49.07398705, 16.81273283 49.07383898, 16.8127217 49.07379905, 16.81271353 49.07375388, 16.81271113 49.07373669, 16.8126859 49.07348176, 16.81268487 49.07346236, 16.81268366 49.07330211, 16.81268056 49.07322181, 16.81267504 49.07314472, 16.81265014 49.0729764, 16.81264162 49.07292904, 16.81260543 49.07273872, 16.81280206 49.07273136, 16.81280733 49.07273129))";

        string wktBlok = "POLYGON ((16.80459502 49.07918068, 16.80459144 49.07913664, 16.80447417 49.07769977, 16.8043569 49.07626301, 16.80435528 49.07623939, 16.80446511 49.07623378, 16.80496536 49.07621282, 16.80528282 49.07621077, 16.80548629 49.07621609, 16.80562906 49.07622073, 16.80578848 49.07622723, 16.80613022 49.07624096, 16.80646558 49.07626238, 16.80665541 49.07626405, 16.80680863 49.07625725, 16.80712148 49.07625424, 16.80731558 49.07626192, 16.80757042 49.07626565, 16.80787648 49.07627472, 16.80839379 49.07627619, 16.80904509 49.07624896, 16.80980512 49.07620373, 16.81008632 49.0761764, 16.81040025 49.07614181, 16.81047816 49.07613066, 16.81056709 49.07611607, 16.81061032 49.07610292, 16.81064308 49.07608894, 16.81066267 49.0760773, 16.81067691 49.07606311, 16.81068814 49.07603681, 16.8106977 49.07599619, 16.81070615 49.07594613, 16.81070967 49.07589746, 16.81070381 49.07584368, 16.81068447 49.0757282, 16.81067513 49.07567449, 16.81065593 49.0756135, 16.81062323 49.07549472, 16.81050417 49.07509942, 16.81050244 49.07505327, 16.81050193 49.07503178, 16.8104914 49.0750018, 16.81047343 49.07497196, 16.81043549 49.07491863, 16.81039663 49.07485162, 16.8101722 49.07437433, 16.80988806 49.07376375, 16.80987233 49.07373363, 16.80983696 49.07367353, 16.80979663 49.07360846, 16.80975071 49.0735433, 16.80970362 49.07347723, 16.80965931 49.07341486, 16.80949465 49.07322294, 16.80948011 49.07319748, 16.80947491 49.07318309, 16.80947263 49.07316852, 16.80947477 49.07314849, 16.80948254 49.07313088, 16.80949592 49.07312032, 16.80951913 49.07310579, 16.80954652 49.07309442, 16.80958167 49.07308304, 16.80962319 49.07307508, 16.80966719 49.07307011, 16.80970908 49.07306673, 16.80994937 49.07306242, 16.81023466 49.07305659, 16.81045834 49.07305219, 16.81072761 49.0730338, 16.8111058 49.07296956, 16.81139429 49.07289754, 16.81156154 49.07286016, 16.81175393 49.07282451, 16.81226697 49.0727514, 16.81280206 49.07273136, 16.81309313 49.07272749, 16.81313119 49.07292478, 16.81315019 49.07302325, 16.81315968 49.07307296, 16.813167 49.07312246, 16.81317335 49.07321113, 16.81317671 49.07329822, 16.81317795 49.07346176, 16.81320318 49.07371669, 16.81321134 49.07376186, 16.8132674 49.07390993, 16.81335772 49.07405905, 16.8134846 49.07419145, 16.81348811 49.07419511, 16.81348814 49.07419514, 16.81348861 49.07419563, 16.81373393 49.07440992, 16.81386098 49.0745379, 16.81398696 49.07467835, 16.8140437 49.07475395, 16.81407921 49.07480593, 16.81410941 49.07485487, 16.81415209 49.07494396, 16.81418903 49.07503138, 16.8142261 49.07512036, 16.81425679 49.07519422, 16.81431425 49.07540981, 16.81440961 49.07577101, 16.81442727 49.07585524, 16.81444167 49.07594098, 16.81447402 49.07613225, 16.81448645 49.07622321, 16.81449562 49.07629207, 16.81449816 49.07636811, 16.8144948 49.07641, 16.81449031 49.07645347, 16.81446436 49.0765487, 16.81445916 49.07657095, 16.81444338 49.07663161, 16.81441908 49.07669868, 16.81436093 49.07683065, 16.81434392 49.076888, 16.81434048 49.076899, 16.8143267 49.07697003, 16.81431053 49.07702466, 16.81429021 49.07709475, 16.81428348 49.07712968, 16.81427689 49.07716885, 16.81427583 49.0771956, 16.81427446 49.07723014, 16.81427729 49.07725868, 16.81428386 49.07731155, 16.81430342 49.07740268, 16.81431793 49.07751314, 16.81432812 49.07754587, 16.81435439 49.07760894, 16.81436583 49.07765566, 16.81438674 49.0777519, 16.81439433 49.07781334, 16.81439497 49.07788151, 16.81438945 49.07793468, 16.81437694 49.07803316, 16.81436726 49.07812488, 16.81436878 49.07819405, 16.81436593 49.07825128, 16.81439342 49.07842424, 16.81446273 49.07866127, 16.81447767 49.07870585, 16.81451684 49.07888481, 16.81450292 49.07922454, 16.81446487 49.07921329, 16.81439016 49.07918388, 16.81432168 49.07914803, 16.81424543 49.07910012, 16.8140179 49.07893153, 16.81401456 49.0789292, 16.81388638 49.07884487, 16.81388115 49.07884171, 16.81378769 49.07879004, 16.81378279 49.07878753, 16.81371845 49.07875714, 16.81371311 49.07875483, 16.81363665 49.07872435, 16.81363511 49.07872375, 16.81363018 49.07872202, 16.81355666 49.07869837, 16.81354923 49.07869626, 16.81345156 49.07867234, 16.81345096 49.0786722, 16.81344361 49.07867069, 16.81329455 49.07864513, 16.81328784 49.07864418, 16.81299367 49.07861122, 16.8127637 49.07857865, 16.81261452 49.0785512, 16.81250319 49.07852616, 16.8122999 49.07846953, 16.81204195 49.07839057, 16.81170704 49.07827698, 16.81143352 49.07817122, 16.81143292 49.078171, 16.81142609 49.07816868, 16.8113236 49.0781378, 16.81131745 49.07813615, 16.81131618 49.07813585, 16.81123741 49.07811797, 16.81123338 49.07811714, 16.81108583 49.07808913, 16.81108294 49.07808862, 16.81074332 49.07803289, 16.81073842 49.07803219, 16.81056118 49.07801055, 16.81056043 49.07801046, 16.8105573 49.07801014, 16.81032716 49.07798947, 16.81018517 49.0779705, 16.81002026 49.07793466, 16.80977876 49.07786185, 16.80977344 49.07786039, 16.80959197 49.07781552, 16.80959112 49.07781531, 16.80958294 49.07781367, 16.80943919 49.07779008, 16.80943355 49.07778929, 16.80943071 49.077789, 16.80933438 49.07778011, 16.80932746 49.07777967, 16.80925594 49.07777709, 16.80925144 49.07777701, 16.80924937 49.07777703, 16.80918145 49.07777815, 16.80917458 49.07777845, 16.8091104 49.07778307, 16.80910508 49.07778357, 16.80910173 49.077784, 16.8090051 49.07779793, 16.80900002 49.07779878, 16.80899503 49.07779984, 16.80885881 49.07783186, 16.80885284 49.07783345, 16.80864701 49.07789461, 16.80864092 49.07789663, 16.80855413 49.07792836, 16.808553 49.07792877, 16.80854935 49.07793024, 16.80849212 49.07795465, 16.80848715 49.07795694, 16.80842748 49.07798673, 16.80842421 49.07798845, 16.80832162 49.07804558, 16.80826156 49.07807492, 16.80810272 49.07813571, 16.8080969 49.07813815, 16.80803131 49.07816837, 16.80802666 49.07817068, 16.80791564 49.07822997, 16.80782837 49.07826535, 16.80776744 49.07828148, 16.80766628 49.07830024, 16.80745625 49.07832997, 16.8074521 49.07833064, 16.80737634 49.0783441, 16.80737119 49.07834514, 16.80725993 49.07837016, 16.80725845 49.07837051, 16.80725213 49.07837222, 16.80715261 49.07840238, 16.80714591 49.07840466, 16.80707689 49.0784309, 16.8070719 49.07843296, 16.80700188 49.07846418, 16.80699652 49.07846679, 16.80688626 49.07852508, 16.80688496 49.07852579, 16.80688143 49.07852785, 16.80673111 49.0786209, 16.80672877 49.0786224, 16.80651167 49.07876795, 16.80640551 49.07883512, 16.80630253 49.0788933, 16.80619492 49.0789442, 16.80610778 49.07897718, 16.80601365 49.07900523, 16.8056903 49.07908004, 16.80568949 49.07908023, 16.80568516 49.07908136, 16.80544146 49.07915014, 16.8053119 49.0791829, 16.80523314 49.07919744, 16.80515855 49.07920577, 16.80509242 49.0792089, 16.8050121 49.0792085, 16.80492479 49.07920495, 16.80470157 49.07918959, 16.80459502 49.07918068), (16.81294122 49.07786349, 16.81291337 49.07785888, 16.81288071 49.07786458, 16.81285586 49.07787917, 16.81284081 49.07789262, 16.81288309 49.07797206, 16.81294474 49.07804509, 16.81299885 49.07803847, 16.81304256 49.07801823, 16.81300193 49.07794009, 16.81294122 49.07786349), (16.81004379 49.07727585, 16.81005086 49.07727578, 16.81006057 49.07727591, 16.81007026 49.07727628, 16.81007993 49.07727691, 16.81008955 49.07727779, 16.81080804 49.07735283, 16.81104116 49.07736528, 16.81096116 49.07732733, 16.81054129 49.07720495, 16.81022466 49.07714154, 16.80991102 49.07710705, 16.80951912 49.0771026, 16.80911318 49.07713892, 16.80869425 49.07721092, 16.80738866 49.07751889, 16.80741572 49.07752126, 16.80773681 49.07749862, 16.80880569 49.07735304, 16.80881292 49.07735213, 16.80882008 49.07735137, 16.80948581 49.07728747, 16.80948664 49.0772874, 16.80948753 49.07728731, 16.80949108 49.07728703, 16.80949468 49.07728673, 16.80949557 49.07728667, 16.80949644 49.0772866, 16.80949999 49.07728641, 16.80950358 49.0772862, 16.80950446 49.07728616, 16.80950539 49.07728611, 16.80950899 49.077286, 16.8095125 49.07728588, 16.80951345 49.07728586, 16.80951436 49.07728584, 16.81004379 49.07727585))";

​

        // Vytvoření instance třídy Geometry z WKT

        var reader = new WKTReader();

        Geometry geometryParcela = reader.Read(wktParcela);

        Geometry geometryBlok = reader.Read(wktBlok);

        var factory = new GeometryFactory();

        var ubionBLokParcela = OverlayNG.Overlay(geometryParcela, geometryBlok, OverlayNG.UNION);

        var intersectionBLokParcela = OverlayNG.Overlay(geometryParcela, ubionBLokParcela, OverlayNG.SYMDIFFERENCE);

        // Určení vzdálenosti pro zarovnání

        double distance = 0.001;

​

        // Vytvoření bufferu kolem většího polygonu

        var buffer = geometryParcela.Buffer(distance);

​

        // Najít průnik bufferu a menšího polygonu

        var intersec = buffer.Intersection(geometryBlok);

​

        // Zarovnání hranic menšího polygonu na hranici většího polygonu

        var alignedPolygon = OverlayNG.Overlay(geometryParcela, intersec, OverlayNG.UNION);

        var symDiff = OverlayNG.Overlay(geometryParcela, intersec, OverlayNG.SYMDIFFERENCE);

        // Najít části hranice menšího polygonu, které jsou v určité vzdálenosti od většího polygonu

        var segmentsToBuffer = new List<LineString>();

        foreach (var segment in GetLineSegments(geometryParcela))

        {

            if (segment.Distance(geometryBlok) <= distance)

            {

                segmentsToBuffer.Add(segment);

            }

        }

​

        // Vytvoření bufferu kolem vybraných částí hranice menšího polygonu

        var buffer2 = factory.CreateMultiLineString(segmentsToBuffer.ToArray()).Buffer(distance);

​

        // Najít průnik bufferu a většího polygonu

        var intersection2 = buffer2.Intersection(geometryBlok);

​

        // Zarovnání hranic menšího polygonu na hranici většího polygonu

        var alignedPolygon2 = OverlayNG.Overlay(geometryParcela, intersection2, OverlayNG.UNION);

​

​

        Console.WriteLine("buffer: " + buffer);

        Console.WriteLine("intersection: " + intersec);

        Console.WriteLine("alignedPolygon: " + alignedPolygon);

        Console.WriteLine("symDiff: " + symDiff);

        Console.WriteLine("ubionBLokParcela: " + ubionBLokParcela);

        Console.WriteLine("intersectionBLokParcela: " + intersectionBLokParcela);

        // Vytvoření dvou sousedících polygonů

        var polygon1 = factory.CreatePolygon(new Coordinate[]

        {

            new Coordinate(0, 0),

            new Coordinate(4, 0),

            new Coordinate(4, 4),

            new Coordinate(0, 4),

            new Coordinate(0, 0)

        });

​

        var polygon2 = factory.CreatePolygon(new Coordinate[]

        {

            new Coordinate(4, 0),

            new Coordinate(8, 0),

            new Coordinate(8, 4),

            new Coordinate(4, 4),

            new Coordinate(4, 0)

        });

​

        // Zarovnání geometrie polygonů

        var union = OverlayNGRobust.Union(new[] { polygon1, polygon2 });

        var intersection = OverlayNG.Overlay(polygon1, polygon2, OverlayNG.INTERSECTION);

​

        // Vytvoření nového polygonu, který převezme společnou část

        var alignedPolygon1 = OverlayNG.Overlay(polygon1, intersection, OverlayNG.DIFFERENCE);

        //var alignedPolygon2 = OverlayNG.Overlay(alignedPolygon1, intersection, OverlayNG.UNION);

​

        // Výstup výsledných polygonů

        Console.WriteLine("Polygon 1: " + alignedPolygon1);

        //Console.WriteLine("Polygon 2: " + alignedPolygon2);

        static List<LineString> GetLineSegments(Geometry polygon)

        {

            var segments = new List<LineString>();

            var coordinates = polygon.Coordinates;

            for (int i = 0; i < coordinates.Length - 1; i++)

            {

                segments.Add(new LineString(new[] { coordinates[i], coordinates[i + 1] }));

            }

            return segments;

        }

    }

}

​