將狀態平面坐標（阿拉斯加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/