Hi Freddy, to add borders to an object of class geodata just add them as an element to the list (class geodata basically is a list): obj$borders <- border.polygon border.polygon here must match coordinates as described in ?polygon, if I remember correctly. regards, Tom