Uso de proj para convertir EPSG:2263 a WGS84

Question

Tengo datos de entrada en EPSG:2263 con coordenadas X e Y como:

X,Y
987043,206407
982441,212359
990907,209020

(estos lugares están dentro de NYC).

Necesito convertirlos a WGS84. Con QGIS, si importo estas coordenadas y le digo que los datos de entrada están en "EPSG:2263, NAD83 / New York Long Island (ftUS)", importa correctamente las coordenadas y las representa (por ejemplo, si se superponen a un mapa WGS84).

Me da esto como las opciones de conversión:

+proj=lcc +lat_1=41.03333333333333 +lat_2=40.66666666666666 +lat_0=40.16666666666666  +lon_0=-74 +x_0=300000.0000000001 +y_0=0 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=us-ft +no_defs

Sin embargo, si intento convertir manualmente un punto en la línea de comandos, por ejemplo:

cs2cs +proj=lcc +lat_1=41.03333333333333 +lat_2=40.66666666666666 +lat_0=40.16666666666666 +lon_0=-74 +x_0=300000.0000000001 +y_0=0 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=us-ft +no_defs -f '%.8f'
987043,206407
-73.99000712    40.16666623 0.00000000,206407

El punto resultante es obviamente no correcta y (si se muestra en Google Earth) está frente a la costa de Nueva Jersey.

Admito que las proyecciones y los sistemas de coordenadas y la conversión entre ellos todavía me desconciertan bastante...

Sugerencias sobre cómo debería ajustarlos cs2cs opciones para conseguir una conversión correcta?

Answer 1

La sintaxis correcta es:

cs2cs +proj=lcc +lat_1=41.03333333333333 +lat_2=40.66666666666666 +lat_0=40.16666666666666 +lon_0=-74 +x_0=300000 +y_0=0 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=us-ft +no_defs +to +init=epsg:4326 -f %%.8f<nad.txt >out.txt

con

987043 206407

en el archivo de entrada nad.txt para obtener este resultado:

-73.98992216    40.73321572 0.00000000

Tenga en cuenta que los valores de entrada deben estar delimitados por espacios en blanco, no por comas.