I had an earlier correspondence with someone on this list that I haven't
recently been able to contact. I am still having issues gluing a projection
on to the shapefiles I am writing from R regarding KUDs. Everytime I try to
load one onto ArcGIS it gives me a message saying the file is missing
spatial reference data so it won't overlay on to my map. I attached the
script, a head of my data, and the shapefile that I have exported from R.
If anyone sees anything I am doing wrong I would really appreciate it.
data <- read.table(file.choose(), header=T, na.strings=c("", "NA"),
sep="\t")
lapply(c('adehabitatLT','adehabitatHR','sp','rgdal'), library,
character.only=T)
#==convert data to spatial points dataframe
xy=data[9:10] # coords
df=data[-9:-10] # other data variables
spdf.x=SpatialPointsDataFrame(coords=xy,data=df)
#==Estimate utilization distribution
kud <-
kernelUD(spdf.x[,4],h="href",hlim=c(0,1000),same4all=T,kern='bivnorm')
ud.1=estUDm2spixdf(kud)
writeOGR(ud.1,"KUD",driver="ESRI Shapefile")
#==home range from UD using LSCV
home=getverticeshr(kud,percent=c(95))
proj4string(home)=CRS("+init=epsg:27573 +proj=utm +zone=11 +ellps=WGS84
+units=m +no_defs")
writePolyShape(home,"HOME")
Date TimeGMT DTHI DTConvertHI TConvertHI Frequency Signal
Gain (dB) Latitude Longitude Diel
6/14/2013 59:59.7 06142013 23:59:59.719 06142013 13:59:59.719
14:00:00.000 63 71 6 21.4326 -157.79103 AM
6/14/2013 59:57.7 06142013 23:59:57.678 06142013 13:59:57.678
13:59:58.000 63 70 6 21.4326 -157.79103 AM
6/14/2013 59:56.7 06142013 23:59:56.658 06142013 13:59:56.658
13:59:57.000 63 73 6 21.4326 -157.79103 AM
6/14/2013 59:54.6 06142013 23:59:54.618 06142013 13:59:54.618
13:59:55.000 63 75 6 21.4326 -157.79103 AM
6/14/2013 59:53.6 06142013 23:59:53.598 06142013 13:59:53.598
13:59:54.000 63 73 6 21.4326 -157.79103 AM
6/14/2013 59:52.6 06142013 23:59:52.578 06142013 13:59:52.578
13:59:53.000 63 74 6 21.4326 -157.79103 AM
6/14/2013 59:51.6 06142013 23:59:51.558 06142013 13:59:51.558
13:59:52.000 63 78 6 21.4326 -157.79103 AM
6/14/2013 59:21.0 06142013 23:59:20.951 06142013 13:59:20.951
13:59:21.000 63 72 12 21.43263 -157.79091 AM
6/14/2013 59:19.9 06142013 23:59:19.932 06142013 13:59:19.932
13:59:20.000 63 70 12 21.43263 -157.79091 AM
6/14/2013 59:00.6 06142013 23:59:00.562 06142013 13:59:00.562
13:59:01.000 63 74 12 21.43262 -157.7908 AM
6/14/2013 58:58.5 06142013 23:58:58.522 06142013 13:58:58.522
13:58:59.000 63 80 12 21.43262 -157.7908 AM
6/14/2013 58:57.5 06142013 23:58:57.502 06142013 13:58:57.502
13:58:58.000 63 79 12 21.43262 -157.7908 AM
6/14/2013 58:56.5 06142013 23:58:56.482 06142013 13:58:56.482
13:58:56.000 63 82 12 21.43262 -157.7908 AM
6/14/2013 58:55.5 06142013 23:58:55.462 06142013 13:58:55.462
13:58:55.000 63 70 12 21.43262 -157.7908 AM
6/14/2013 58:48.3 06142013 23:58:48.324 06142013 13:58:48.324
13:58:48.000 63 75 12 21.43262 -157.7908 AM
6/14/2013 58:47.3 06142013 23:58:47.305 06142013 13:58:47.305
13:58:47.000 63 80 12 21.43262 -157.7908 AM
_______________________________________________
R-sig-ecology mailing list
[email protected]
https://stat.ethz.ch/mailman/listinfo/r-sig-ecology
