¿Cómo cambiar un srs existente en QGis?

Question

Yo estoy mostrando a postgis datos con srid 28992 (Amersfoort / RD Nueva). Pero cuando me la trama de datos geoespaciales en un mapa (con openlayers plugin), se muestran los datos de un hunderd metros de frente porque el dato de turno no es tomado en cuenta.

Cuando hago una costumbre CRS con el datum correcto desplazamiento y puse mi capa para que CRS es mostrado correctamente.

Quiero actualizar QGis' interpretación estándar de 28992 por lo que toma en cuenta este dato cambio (+towgs84 parámetro) correctamente y cuando abro de nuevo los datos se muestran correctamente inmediatamente.

He encontrado que QGis tiene una sqlite db llamado srs.db en /usr/share/qgis/recursos/. Pero incluso cuando yo actualice el registro para 28992 todavía no muestra los datos correctamente en el mapa. Cuando me inspeccionar la capa en QGis veo que todavía no tiene +towgs84 parámetro.

¿Significa esto que QGis no uso de este db sqlite? Debido a que los datos que muestro es de una base de datos postgresql también he editado spatial_ref_sys.sql para asegurarse de que también se conoce acerca del cambio de datum.

También he editado el epsg es esri archivos en /usr/share/proj/ sólo para estar seguro, pero ningún efecto.

En este momento no se me ocurre otro lugar donde QGis obtiene su información de proyección.

Answer 1

Está seguro de que su Postgis y QGIS comparten el mismo +towgs84 parámetros?

Mi QGIS 2.0.1 Dufour viene con esta definición:

+proj=sterea +lat_0=52.15616055555555 +lon_0=5.38763888888889 +k=0.9999079 +x_0=155000 +y_0=463000 +ellps=bessel +towgs84=565.417,50.3319,465.552,-0.398957,0.343988,-1.8774,4.0725 +units=m +no_defs

mientras mi PostgreSQL tiene esta almacenado en la tabla spatial_ref_sys

+proj=sterea +lat_0=52.15616055555555 +lon_0=5.38763888888889 +k=0.9999079 +x_0=155000 +y_0=463000 +ellps=bessel +towgs84=565.237,50.0087,465.658,-0.406857,0.350733,-1.87035,4.0812 +units=m +no_defs

Si su QGIS tiene otro proj definición para EPSG:28992, su instalación puede ser roto. Tener una mirada en el registro de instalación si postinstall fue roto en alguna parte.

Otro punto: Postgis utiliza esta cadena para el WKT definición:

PROJCS["Amersfoort / RD New",
  GEOGCS["Amersfoort",
    DATUM["Amersfoort",
      SPHEROID["Bessel 1841",6377397.155,299.1528128,
      AUTHORITY["EPSG","7004"]],
    AUTHORITY["EPSG","6289"]],
    PRIMEM["Greenwich",0,AUTHORITY["EPSG","8901"]],
    UNIT["degree",0.01745329251994328,AUTHORITY["EPSG","9122"]],
    AUTHORITY["EPSG","4289"]],
    UNIT["metre",1,AUTHORITY["EPSG","9001"]],
  PROJECTION["Oblique_Stereographic"],
    PARAMETER["latitude_of_origin",52.15616055555555],
    PARAMETER["central_meridian",5.38763888888889],
    PARAMETER["scale_factor",0.9999079],
    PARAMETER["false_easting",155000],
    PARAMETER["false_northing",463000],
    AUTHORITY["EPSG","28992"],
    AXIS["X",EAST],
    AXIS["Y",NORTH]]

que las referencias a Amersfoort dato, pero no +towgs84 se almacenan para que. No estoy seguro de si QGIS, por tanto, los conjuntos de +towgs84 a cero, en lugar de utilizar su propio EPSG valores, o a aquellas desde las Postgis proj4 cadena.

QGIS escribe esto .qpj de salida:

PROJCS["Amersfoort / RD New",
  GEOGCS["Amersfoort",
    DATUM["Amersfoort",
      SPHEROID["Bessel 1841",6377397.155,299.1528128,
      AUTHORITY["EPSG","7004"]],
    TOWGS84[565.417,50.3319,465.552,-0.398957,0.343988,-1.8774,4.0725],
    AUTHORITY["EPSG","6289"]],
    PRIMEM["Greenwich",0,AUTHORITY["EPSG","8901"]],
    UNIT["degree",0.0174532925199433,AUTHORITY["EPSG","9122"]],
    AUTHORITY["EPSG","4289"]],
PROJECTION["Oblique_Stereographic"],
  PARAMETER["latitude_of_origin",52.15616055555555],
  PARAMETER["central_meridian",5.38763888888889],
  PARAMETER["scale_factor",0.9999079],
  PARAMETER["false_easting",155000],
  PARAMETER["false_northing",463000],
  UNIT["metre",1,AUTHORITY["EPSG","9001"]],
  AXIS["X",EAST],
  AXIS["Y",NORTH],
  AUTHORITY["EPSG","28992"]]