Obawy zgłoszone przez @Adam Cara są bardzo ważne. Jednak gdy spojrzałem na samo skrzyżowanie, okazało się, że większość z nich to albo chude bity na granicy między 2 masywnymi wielokątami, albo odłączone bity / wyspy dużego wielokąta siedzącego wewnątrz innego. W tym przypadku kasowanie nakładek i przypisywanie rozpuszczonych nakładek może działać atrybut najbliższego sąsiada (w zależności od tego, co nastąpi wcześniej).
Proces:
arcpy.Intersect_analysis("Landuse2b #","D:/Scratch/mpart.shp","ALL","#","INPUT")
arcpy.MultipartToSinglepart_management("mpart","D:/Scratch/spart.shp")
arcpy.Dissolve_management("spart","D:/Scratch/dissolved_1.shp","#","#","SINGLE_PART","DISSOLVE_LINES")
arcpy.Erase_analysis("Landuse2b","dissolved_1","D:/Scratch/erased.shp","#")
arcpy.DeleteField_management("spart","FID_Landus;OBJECTID;Lnduse2013;Shape_Leng;Shape_Area")
arcpy.SpatialJoin_analysis("spart","erased","D:/Scratch/sj_1.shp","JOIN_ONE_TO_ONE","KEEP_ALL"…..)
arcpy.Merge_management("sj;erased","D:/Scratch/merge_1.shp"….)
arcpy.Dissolve_management("merge_1","D:/Scratch/dissolved.shp","Lnduse2013","#","SINGLE_PART","DISSOLVE_LINES")
Test:
arcpy.Intersect_analysis("dissolved #","D:/Scratch/test.shp","ALL","#","INPUT")
Wiadomość:
- Wykonywanie: Przecięcie „rozwiązany #” D: \ Scratch \ test.shp ALL # INPUT
- Czas rozpoczęcia: Pt 22 kwietnia 10:16:07 2016
- Funkcje czytania ...
- Funkcje pękania ...
- Składanie funkcji ...
- OSTRZEŻENIE 000117: Ostrzeżenie wygenerowane puste wyjście.
- Udało się w piątek 22 kwietnia 10:16:09 2016 (Czas, który upłynął: 2,34 sekundy)
Zauważ, że ten komunikat oznacza, że nie ma przecinających się danych wyjściowych z końcowego rozpadu
Wynik:
https://drive.google.com/file/d/0B6B33FrGNitPT0FvbVdHQm0xMDg/view?usp=sharing