I'm using Small Area shapefiles in Ireland from the Central Statistics Office (CSO) website. The link to these is here http://census.cso.ie/censusasp/saps/boundaries/ED_SA%20Disclaimer1.htm (Look for the small area ones)
I can't seem to convert the coordinates to latitude longitude, can anyone see what I'm doing wrong?
I've loaded these in using R maptools package using the following code:
But when I check the
cte@proj4string it reurns NA. Then I found this document which suggests the projection used is EPSG:3035
I then found a site which gives the proj4 string for EPSG:3035. I then tried to set the coorinates as this string.
cte@proj4string <- CRS("+proj=laea +lat_0=52 +lon_0=10 +x_0=4321000 +y_0=3210000 +ellps=GRS80 +units=m +no_defs" )
Finally I tried to transform the object to latlong coordinates using the following:
cte <- spTransform(cte, CRS("+proj=longlat +datum=WGS84"))
I also tried:
coordinates(cte) <- project(cte, "+proj=longlat +datum=WGS84")
The first transformed it to new coordinates, but they don't seem to be correct. The second one just put everything as infinite.
Can anybody let me know what I'm doing wrong?