¿Qué código EPGS es el CRS real en este GeoPDF?

Question

Estoy teniendo problemas reales, tratando de adaptarse a un GeoPDF a un conjunto de archivos de forma a partir gdaltransform y gdalwarp. El problema me parece que no puede solucionar es que cada vez que ejecuto los comandos, ya sea en QGIS o de la línea de comandos, el nuevo GeoTiff es siempre por un par de cientos de metros. He leído acerca del uso de la +towgs84 opciones a la hora de establecer el CRS, pero la documentación en línea para esta opción es muy vaga acerca de cómo configurar el 3-7 parámetros y que el lugar para ponerlos.

En completo y absoluto desaliento, empecé a pensar que tal vez no soy la correcta utilización de la CRS en el primer lugar. Cuando se ejecuta gdalinfo, veo a los múltiples códigos EPSG. Yo había estado usando EPSG:4267, pensando que esta era la real CRS de la GeoPDF. Pero si EPSG:4267 es la correcta CRS, entonces ¿por qué el archivo de seguir mostrando un desplazamiento de ~188 metros, cuando tratando de igualar a mi shapefiles en EPSG:4326?

Aquí está el gdalinfo salida de mi GeoPDF archivo. Por favor alguien puede identificar cuál es la correcta CRS del archivo, y cómo puedo usar el CRS para eliminar el offset?

Driver: PDF/Geospatial PDF Files: ~/Desktop/myfile.pdf Size is 4470, 3375 Coordinate System is: PROJCS["UTM Zone 18, Northern Hemisphere", GEOGCS["NAD27", DATUM["North_American_Datum_1927", SPHEROID["Clarke 1866",6378206.4,294.9786982138982, AUTHORITY["EPSG","7008"]], AUTHORITY["EPSG","6267"]], PRIMEM["Greenwich",0, AUTHORITY["EPSG","8901"]], UNIT["degree",0.0174532925199433, AUTHORITY["EPSG","9122"]], AUTHORITY["EPSG","4267"]], PROJECTION["Transverse_Mercator"], PARAMETER["latitude_of_origin",0], PARAMETER["central_meridian",-75], PARAMETER["scale_factor",0.9996], PARAMETER["false_easting",500000], PARAMETER["false_northing",0], UNIT["Meter",1]] GeoTransform = 786521.192808, 8.469778080000001, 0.15067296 2179199.881801, 0.14403264, -8.469838560000001 GCP Projection = PROJCS["UTM Zone 18, Northern Hemisphere", GEOGCS["NAD27", DATUM["North_American_Datum_1927", SPHEROID["Clarke 1866",6378206.4,294.9786982138982, AUTHORITY["EPSG","7008"]], AUTHORITY["EPSG","6267"]], PRIMEM["Greenwich",0, AUTHORITY["EPSG","8901"]], UNIT["degree",0.0174532925199433, AUTHORITY["EPSG","9122"]], AUTHORITY["EPSG","4267"]], PROJECTION["Transverse_Mercator"], PARAMETER["latitude_of_origin",0], PARAMETER["central_meridian",-75], PARAMETER["scale_factor",0.9996], PARAMETER["false_easting",500000], PARAMETER["false_northing",0], UNIT["Meter",1]] GCP[ 0]: Id=1, Info= (207.8125,2466.92708333333) -> (788652.79799,2158337.734581,0) GCP[ 1]: Id=2, Info= (211.197916666667,287.760416666667) -> (788355.62363,2176795.854646,0) GCP[ 2]: Id=3, Info= (3308.59375,288.020833333333) -> (814587.658583,2177239.223702,0) GCP[ 3]: Id=4, Info= (3307.8125,2467.70833333333) -> (814911.952376,2158777.947799,0) GCP[ 4]: Id=5, Info= (2275.78125,1377.34375) -> (806001.09572,2167857.005616,0) GCP[ 5]: Id=6, Info= (1242.70833333333,1377.08333333333) -> (797252.559747,2167709.741736,0) Metadata: AUTHOR=NGA CREATION_DATE=D:20071008152059-04'00' KEYWORDS=003 NEATLINE=POLYGON ((788648.68490375 2158330.85305763,788592.08154525 2161760.771554,788380.2253165 2174662.13561613,788346.57604175 2176801.72714438,802820.27008 2177043.44654713,812895.73768325 2177216.99106225,814591.898970625 2177245.83509875,814733.508010375 2169285.51026463,814856.815543875 2162105.92365575,814916.143021875 2158770.92472275,802630.671176375 2158555.38531638,793697.623632125 2158407.88726563,789672.19489725 2158343.84562738,788648.68490375 2158330.85305763,788648.68490375 2158330.85305763,788648.68490375 2158330.85305763,788648.68490375 2158330.85305763)) TITLE=Map Edition Corner Coordinates: Upper Left ( 786521.193, 2179199.882) ( 72d16' 1.60"W, 19d41'19.09"N) Lower Left ( 787029.714, 2150614.177) ( 72d15'59.83"W, 19d25'49.84"N) Upper Right ( 824381.101, 2179843.708) ( 71d54'22.52"W, 19d41'18.91"N) Lower Right ( 824889.622, 2151258.003) ( 71d54'22.82"W, 19d25'49.93"N) Center ( 805705.407, 2165228.942) ( 72d 5'11.67"W, 19d33'34.77"N) Band 1 Block=4470x1 Type=Byte, ColorInterp=Red Band 2 Block=4470x1 Type=Byte, ColorInterp=Green Band 3 Block=4470x1 Type=Byte, ColorInterp=Blue

Answer 1

Sus datos parece estar en Haití, el uso de la North American Datum de 1927.

QGIS y GDAL el uso de los datos de la cuadrícula de cambio de archivos para NAD27 para Estados Unidos y Canadá, pero por desgracia sus datos está justo fuera de la conus de la cuadrícula de archivo. por lo tanto, todos sus reprojections obtener un valor null cambio, que está definitivamente mal.

Clifford J. Mugnier publicado un Cuadrículas Y las Referencias del artículo en Haití en septiembre de 2011, dando a cambio de estos valores:

from NAD1927 to WGS84: DeltaX =  –3 m ± 3 m, DeltaY = +142 m ± 9 m, 
DeltaZ = +183 m ± 12 m. Curiously, the shift parameters that work out 
for  the single Haïtian point, "FORT NATIONAL"  from NAD1927  to 
WGS84 are: DeltaX =  –13 m, DeltaY =  –95 m, DeltaZ =  –197 m. 

Sus datos están en NAD27/UTM 18N, pero EPSG:26718 falla por las razones dadas anteriormente, y usted tiene que construir -s_srs y -t_srs por cada GDAL-basado en la operación de forma manual:

+proj=utm +zone=18 +ellps=clrk66 +towgs84=-3,142,183,0,0,0 +units=m +no_defs
+proj=utm +zone=18 +ellps=clrk66 +towgs84=-13,-95,-197,0,0,0 +units=m +no_defs

Tomando la Trou du Nord mapa de http://www.lib.utexas.edu/maps/topo/haiti/ y el primer dato cambio, se adapta muy bien a un OpenStreetMap fondo, con el borde interior en rojo:

Answer 2

Pegar el WKT en http://prj2epsg.org da el sonido plausible EPSG: 2032

 PROJCS["NAD27(CGQ77) / UTM zone 18N", 
  GEOGCS["NAD27(CGQ77)", 
    DATUM["North American Datum 1927 (CGQ77)", 
      SPHEROID["Clarke 1866", 6378206.4, 294.9786982138982, AUTHORITY["EPSG","7008"]], 
      AUTHORITY["EPSG","6609"]], 
    PRIMEM["Greenwich", 0.0, AUTHORITY["EPSG","8901"]], 
    UNIT["degree", 0.017453292519943295], 
    AXIS["Geodetic longitude", EAST], 
    AXIS["Geodetic latitude", NORTH], 
    AUTHORITY["EPSG","4609"]], 
  PROJECTION["Transverse Mercator", AUTHORITY["EPSG","9807"]], 
  PARAMETER["central_meridian", -75.0], 
  PARAMETER["latitude_of_origin", 0.0], 
  PARAMETER["scale_factor", 0.9996], 
  PARAMETER["false_easting", 500000.0], 
  PARAMETER["false_northing", 0.0], 
  UNIT["m", 1.0], 
  AXIS["Easting", EAST], 
  AXIS["Northing", NORTH], 
  AUTHORITY["EPSG","2032"]]