### Home ranges of Australian petcats
### Kays R, Dunn RR, Parsons AW, Mcdonald B, Perkins T, Powers S, Shell L, McDonald JL, Cole H, Kikillus H, Woods L, Tindle H, Roetman P (2020) The small home ranges and large local ecological impacts of pet cats. Animal Conservation. https://doi.org/10.1111/acv.12563
###actual data from https://www.movebank.org/cms/webapp?gwt_fragment=page=studies,path=study64283289


library(rworldmap)
data(countryExData) 
sPDF<-joinCountryData2Map(countryExData,joinCode = "ISO3",nameJoinColumn="ISO3V10")


ozzy_cat<-read.csv("Pet-Cats-Australia.csv")
attach(ozzy_cat)
long_range<-range(location.long)
lat_range<-range(location.lat)
kitty<-individual.local.identifier


mapCountryData(sPDF,mapTitle="Australian Moggies", xlim=long_range, ylim=lat_range, colourPalette="white2Black",borderCol="black") 

for(i in 1:length(levels(individual.local.identifier))){
rows<-which(ozzy_cat$individual.local.identifier==levels(individual.local.identifier)[i])
lines(ozzy_cat$location.long[rows], ozzy_cat$location.lat[rows],col=as.numeric(ozzy_cat$individual.local.identifier[rows]),lwd=2)
}

#### Scale too large, limit it to part around median region


medlong<-median(location.long)
medlat<-median(location.lat)
newxlim<-c(medlong-0.4,medlong+0.4)
newylim<-c(medlat-0.4,medlat+0.4)

mapCountryData(sPDF,mapTitle="Australian Moggies", xlim=newxlim, ylim=newylim, colourPalette="white2Black",borderCol="black") 
for(i in 1:length(levels(individual.local.identifier))){
rows<-which(ozzy_cat$individual.local.identifier==levels(individual.local.identifier)[i])
lines(ozzy_cat$location.long[rows], ozzy_cat$location.lat[rows],col=as.numeric(ozzy_cat$individual.local.identifier[rows]),lwd=2)
}



detach(ozzy_cat)