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

Question

我对我认为不知道如何为我的状态平面（阿拉斯加州平面4区NAD83 Feet）找到正确的投影输入以与功能sp::CRS使用感到困扰。

我在以下位置遵循从状态平面转换的说明： 将纬度/经度转换为状态平面坐标

并且我研究了？CRS，它使我从References转到http://trac.osgeo.org/proj/ ，但我什至无法确定该站点是否将满足我的需求。

我正在使用Hadley Wickham的教程来绘制shapefile， 网址为： https : //github.com/hadley/ggplot2/wiki/plotting-polygon-shapefiles

而且我可以从以下位置绘制我喜欢的任何shapefile： http ： //munimaps.muni.org/moagis/download.htm

我的代码：

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() 

请帮助我找到转换投影所需的输入。 我是使用地图投影的新手。

Answer 1

我建议您通过FWTools（在GDAL中，我看到您正在导入）中使用OGR2OGR。 您应该能够利用这些工具从一个投影投射到另一个投影。

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

http://www.mercatorgeosystems.com/blog-articles/2008/05/30/using-ogr2ogr-to-re-project-a-shape-file/