0
У меня есть два полигональных шейп-файла, и я хотел бы закрепить их друг за другом. Я ищу в google, но я мог найти только обрезание ограничивающей рамкой или обрезками точек полигонами, и это не то, что мне нужно. Я также нахожу что-то в другом программировании, кроме R (http://rosettacode.org/wiki/Sutherland-Hodgman_polygon_clipping#Python). Не могли бы вы мне помочь?Как скопировать файл формы многоугольника с помощью другого файла формы полигона в R?
Благодаря Tiago
Попробуйте 'rgeos :: gIntersection()' как demo'd в (например) принятом ответе на [этот вопрос] (http://stackoverflow.com/questions/13982773/crop-for-spatialpolygonsdataframe/13986029# 13986029). Альтернативно, 'растр :: intersect()' также должен делать трюк. –
@ JoshO'Brien Спасибо за ваш ответ. Это помогает мне получить почти то, что я хочу. У меня еще одна небольшая проблема. Когда я делаю 'rgeos :: gIntersection()', я потерял все строки, хранящиеся в таблице атрибутов. Есть ли способ, которым я могу сохранить его в соответствии с тем, что было пересечено? – Tiago
Возможно, попробуйте 'raster :: intersect()', который был специально разработан для хранения связанных атрибутов data.frame. (И единственная причина, по которой я говорю «возможно», заключается в том, что иногда это не работает так же, как «rgeos :: gIntersection()», как, например, с данными в ответе, который я связал выше. Но он обычно делает работа очень приятная.) Надеюсь, что это поможет! –