Comecei comComo posso transformar as coordenadas de um Shapefile? .

A resposta que me levou a seguir [o que eu acho que é] o caminho certo, mas ainda não consegui resolver meu problema.

Uma questão é que ainda não encontrei a projeção correta:https: //gis.stackexchange.com/questions/13330/how-can-i-correctly-transform-unproject-from-lc

EDIT: Essa pergunta no site gis foi respondida e eu pude reproduzir uma transformação correta usando a ferramenta de linha de comando PROJ cs2cs. Se parece com isso

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

gora, que eu tive a transformação correta, pude tentar a mesma coisa de uma forma simples usando o RGe

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] 

Por que os resultados são diferentes entre o RGeo e o cs2cUma vez que posso fazer o RGeo executar a tradução correta, existe uma maneira de criar a fábrica adequada para transformar uma Geometria completa em vez de um ponto? Existe uma ferramenta de linha de comando que eu possa usar como solução alternativa para transformar todos os pontos no meu shapefile para que eu possa seguir adiante com minha vid

Em geral: alguém poderia me instruir sobre como usar corretamente esta biblioteca?

Muito obrigado para olhar