R 将州平面坐标（阿拉斯加4）转换为纬度和经度

我一直在思考一个问题，就是不知道如何为我的州平面（阿拉斯加州平面4区NAD83英尺）找到正确的投影输入，以便与函数

sp:：CRS

一起使用

我将按照以下说明从状态平面进行转换：

我已经调查了？CRS，它把我从参考资料中带到了这里，但我甚至不知道这个网站是否会有我需要的东西

我正在使用Hadley Wickham的教程在以下位置绘制形状文件：

我可以绘制我喜欢的任何形状文件：

我的代码：

require(rgdal)
require(maptools)
require(ggplot2)
require(plyr)

  my_dsn <- "directory where the shapefile is"
  Assembly = readOGR(dsn=my_dsn, layer="assembly")
  Assembly@data$id = rownames(Assembly@data)
  Assembly.points = fortify(Assembly, region="id")
  Assembly.df = join(Assembly.points, Assembly@data, by="id")
  Assembly@data$id = rownames(Assembly@data)
  Assembly.points = fortify(Assembly, region="id")
  Assembly.df = join(Assembly.points, Assembly@data, by="id")
  #Assembly.df$DISTRICT <- factor(Assembly.df$DISTRICT)

  ggplot(Assembly.df) + 
  aes(long,lat,group=group) + 
  geom_path(color="black") +
  coord_equal() 

require（rgdal）
需要（地图工具）
需要（ggplot2）
需要（plyr）
我的_dsn我建议通过FWTools（在GDAL中，我看到您正在导入）使用OGR2OGR。您应该能够利用这些工具从一个投影投影到另一个投影

它应该是一个简单的命令行操作，例如
ogr2ogr-f“ESRI形状文件”original.shp wgs84.shp-s_srs EPSG:27700-t_srs EPSG:4326

proj4string（汇编）在readOGR行之后直接报告什么？如果你不同意的话，所有其他的东西都是没有帮助的结果输出是“+proj=tmerc+lat_0=54+lon_0=-150+k=0.9999+x_0=500000.0000000001+y_0=0+datum=NAD83+units=us ft+no_defs+ellps=GRS80+towgs84=0,0”我猜谷歌地图使用墨卡托投影。如何找到该投影的CRS（）输入？如果使用R，则可以避免使用ogr2ogr，这不一定更可取，但这里要求这样做。此外，FWTools已经完全过时了，最好使用OSGeo4W或其他最新的Windows版本。