Question

我坚持认为不知道如何找到我的状态平面（阿拉斯加州立大学第4区NAD83脚）的正确投影输入以与函数sp::CRS一起使用的问题。

我正在遵循从州立飞机转换的指示： Convert latitude/longitude to state plane coordinates

我已经调查了？CRS从参考文献中把我带到了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使用OGR2OGR（在GDAL中，我看到你正在导入）。您应该能够利用这些工具从一个投影到另一个投影。

它应该是一个简单的命令行操作，例如 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/

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

1 个答案: