0
J'ai deux shapefiles de polygones et j'aimerais les couper les uns par les autres. Je fais une recherche sur google mais je ne peux trouver que l'écrêtage par une boîte de délimitation ou des points de découpage par polygones, et ce n'est pas ce dont j'ai besoin. Je trouve aussi quelque chose dans d'autres langages de programmation, sauf dans R (http://rosettacode.org/wiki/Sutherland-Hodgman_polygon_clipping#Python). Pourriez-vous m'aider?Comment couper un fichier de formes polygone par un autre fichier de formes polygone dans R?
Merci Tiago
Essayez 'rgeos :: gIntersection()' comme demo'd dans (par exemple) la réponse acceptée à [cette question] (http://stackoverflow.com/questions/13982773/crop-for-spatialpolygonsdataframe/13986029# 13986029). Alternativement, 'raster :: intersect()' devrait aussi faire l'affaire. –
@ JoshO'Brien Merci pour votre réponse. Cela m'aide à obtenir presque exactement ce que je veux. Je n'ai toujours qu'un petit problème. Quand je fais le 'rgeos :: gIntersection()' j'ai perdu toutes les lignes stockées dans la table attributaire. Y a-t-il un moyen de le garder correspondant à ce qui a été recoupé? – Tiago
Essayez peut-être 'raster :: intersect()', qui a été expressément conçu pour conserver les attributs data.frame associés. (Et la seule raison pour laquelle je dis "peut-être" est que parfois, cela ne marche pas aussi bien que "rgeos :: gIntersection()", comme par exemple, avec les données de la réponse à laquelle je suis lié ci-dessus. travailler assez bien.) J'espère que ça aide! –