3

我从如何转换 Shapefile 的坐标开始?.

那里的回应使我开始[我认为]正确的轨道,但我仍然无法解决我的问题。

一个问题是我还没有找到正确的投影:https ://gis.stackexchange.com/questions/13330/how-can-i-correctly-transform-unproject-from-lcc

编辑: gis 网站上的那个问题已经得到解答,我能够使用 PROJ 命令行工具 cs2cs 重现正确的转换。它看起来像这样:

larry$  cs2cs -f "%.8f" +proj=lcc +lat_1=37.06666666666667 +lat_2=38.43333333333333 +lat_0=36.5 +lon_0=-120.5 +x_0=2000000 +y_0=500000.0000000002 +ellps=GRS80 +datum=NAD83 +to_meter=0.3048006096012192 +no_defs +to +proj=lonlat +datum=WGS84 +ellps=WGS84
6011287.4999795845 2100857.2499904726
-122.40375492   37.74919006 0.00000000

现在,我进行了正确的转换,我可以使用 RGeo 以简单的形式尝试同样的事情:

ruby-1.9.2-p180 :001 >     projection_str = ' +proj=lcc +lat_1=37.06666666666667 +lat_2=38.43333333333333 +lat_0=36.5 +lon_0=-120.5 +x_0=2000000 +y_0=500000.0000000002 +ellps=GRS80 +datum=NAD83 +to_meter=0.3048006096012192 +no_defs'
 => " +proj=lcc +lat_1=37.06666666666667 +lat_2=38.43333333333333 +lat_0=36.5 +lon_0=-120.5 +x_0=2000000 +y_0=500000.0000000002 +ellps=GRS80 +datum=NAD83 +to_meter=0.3048006096012192 +no_defs" 
ruby-1.9.2-p180 :002 >     projection = RGeo::CoordSys::Proj4.new(projection_str)
 => #<RGeo::CoordSys::Proj4:0x805cba18 " +proj=lcc +lat_1=37.06666666666667 +lat_2=38.43333333333333 +lat_0=36.5 +lon_0=-120.5 +x_0=2000000 +y_0=500000.0000000002 +ellps=GRS80 +datum=NAD83 +to_meter=0.3048006096012192 +no_defs +towgs84=0,0,0"> 
ruby-1.9.2-p180 :003 >     desired_str = '+proj=lonlat +datum=WGS84 +ellps=WGS84'
 => "+proj=lonlat +datum=WGS84 +ellps=WGS84" 
ruby-1.9.2-p180 :004 >     desired = RGeo::CoordSys::Proj4.new(desired_str)
 => #<RGeo::CoordSys::Proj4:0x805271ac " +proj=lonlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0"> 
ruby-1.9.2-p180 :005 >     RGeo::CoordSys::Proj4::transform_coords(projection, desired, 6011287.4999795845, 2100857.2499904726 )
 => [-140.92282523143973, 30.16981659183029] 
  1. 为什么 RGeo 和 cs2cs 的结果不同?
  2. 一旦我可以让 RGeo 执行正确的转换,有没有办法可以创建合适的工厂来转换完整的几何而不是点?
  3. 是否有一个命令行工具可以用作解决方法来转换我的 shapefile 中的所有点,以便我可以继续我的生活?

一般来说:有人可以指导我如何正确使用这个库吗?

非常感谢您的关注。

4

1 回答 1

2

作为黑暗中的狂野刺,因为我不知道 RGeo 甚至 Ruby,请尝试将您的坐标以英尺替换为等效的米:1832244.0944819663048746863094224、640342.57048223700783128534419392(尽管您可能不需要那个小数位数......)另一种可能性是交换坐标——也许 RGeo 做了一些非常规的假设。

如果您能够从 Ruby 调用可执行文件,您可以简单地使用ogr2ogr来转换您的 shapefile。

于 2011-08-10T17:06:19.320 回答