Sto cercando di creare una sovrapposizione spaziale tra due poligoni ("grd" e "range") usando la funzione over (rgeos).Come posso gestire un "foro orfano" in uno shapefile poligonale in R o ArcGIS?
Questo codice funziona per tutti i file tranne uno:
range.grd.overlay <- grd[!is.na(over(grd, range)), ]@data$grdID # get gridcells that overlap with range
e ottenere questo messaggio di errore:
Error in createPolygonsComment(p) :
rgeos_PolyCreateComment: orphaned hole, cannot find containing polygon for hole at index 9
Capisco che c'è un buco, senza un poligono che circonda da qualche parte nel shapefile, ma Non riesco a capire come aggirarlo. Posso semplicemente cancellare il foro nell'indice 9? Come ? Se si dispone di una soluzione per correggere lo shapefile in ArcGIS, anche per me andrebbe bene.
Grazie mille in anticipo per il vostro aiuto !!
saluti, Carsten
nessuno può aiutare con questo, dal momento che non abbiamo alcuna idea sulla struttura di GRD o gamma. Molti oggetti possono avere uno slot @ data $ grdID. Questo ha bisogno di lavoro per fare un esempio riproducibile. Utilizzare QGIS (gratuito) con Geometry Tools o Collettore (economico) con Topology Factory per un controllo geometrico della geometria. – mdsumner
Sembra che cancellare il buco all'indice 9 dovrebbe funzionare. Cosa succede quando ci provi? 'subset' nel pacchetto' taRifx' dovrebbe farlo se non si dispone di un modo migliore (non sono aggiornato sui miei metodi spaziali in questi giorni). –
Ok, vedo, scusa per non essere molto specifico sul mio problema: – user1012744