Revisions to How do I change an existing srs in QGis?

added 1093 characters in body

Source Link

edited Jan 6, 2014 at 16:14

77.1k
5
88
163

Are you sure that your Postgis and QGIS share the same +towgs84 parameters?

My QGIS 2.0.1 Dufour comes with this definition:

+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

while my PostgreSQL has this stored in table 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

If your QGIS has another proj definition for EPSG:28992, your installation might be broken. Have a look at the setup log if postinstall was broken somewhere.

One other point: Postgis uses this string for the WKT definition:

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]]

which references to Amersfoort datum, but no +towgs84 are stored for that. I am not sure if QGIS therefore sets +towgs84 to zero, instead of using its own EPSG values, or those from the Postgis proj4 string.

QGIS writes this .qpj output:

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"]]

Are you sure that your Postgis and QGIS share the same +towgs84 parameters?

My QGIS 2.0.1 Dufour comes with this definition:

+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

while my PostgreSQL has this stored in table 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

If your QGIS has another proj definition for EPSG:28992, your installation might be broken. Have a look at the setup log if postinstall was broken somewhere.

One other point: Postgis uses this string for the WKT definition:

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]]

which references to Amersfoort datum, but no +towgs84 are stored for that. I am not sure if QGIS therefore sets +towgs84 to zero, instead of using its own EPSG values, or those from the Postgis proj4 string.

Are you sure that your Postgis and QGIS share the same +towgs84 parameters?

My QGIS 2.0.1 Dufour comes with this definition:

+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

while my PostgreSQL has this stored in table 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

If your QGIS has another proj definition for EPSG:28992, your installation might be broken. Have a look at the setup log if postinstall was broken somewhere.

One other point: Postgis uses this string for the WKT definition:

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]]

which references to Amersfoort datum, but no +towgs84 are stored for that. I am not sure if QGIS therefore sets +towgs84 to zero, instead of using its own EPSG values, or those from the Postgis proj4 string.

QGIS writes this .qpj output:

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"]]

added 924 characters in body

Source Link

edited Jan 6, 2014 at 15:59

AndreJ

77.1k
5
88
163

Are you sure that your Postgis and QGIS share the same +towgs84 parameters?

My QGIS 2.0.1 Dufour comes with this definition:

+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

while my PostgreSQL has this stored in table 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

If your QGIS has another proj definition for EPSG:28992, your installation might be broken. Have a look at the setup log if postinstall was broken somewhere.

One other point: Postgis uses this string for the WKT definition:

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]]

which references to Amersfoort datum, but no +towgs84 are stored for that. I am not sure if QGIS therefore sets +towgs84 to zero, instead of using its own EPSG values, or those from the Postgis proj4 string.

Are you sure that your Postgis and QGIS share the same +towgs84 parameters?

My QGIS 2.0.1 Dufour comes with this definition:

+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

while my PostgreSQL has this stored in table 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

If your QGIS has another proj definition for EPSG:28992, your installation might be broken. Have a look at the setup log if postinstall was broken somewhere.

Are you sure that your Postgis and QGIS share the same +towgs84 parameters?

My QGIS 2.0.1 Dufour comes with this definition:

+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

while my PostgreSQL has this stored in table 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

If your QGIS has another proj definition for EPSG:28992, your installation might be broken. Have a look at the setup log if postinstall was broken somewhere.

One other point: Postgis uses this string for the WKT definition:

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]]

which references to Amersfoort datum, but no +towgs84 are stored for that. I am not sure if QGIS therefore sets +towgs84 to zero, instead of using its own EPSG values, or those from the Postgis proj4 string.

Source Link

created Jan 6, 2014 at 15:43

AndreJ

77.1k
5
88
163

Are you sure that your Postgis and QGIS share the same +towgs84 parameters?

My QGIS 2.0.1 Dufour comes with this definition:

+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

while my PostgreSQL has this stored in table 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

If your QGIS has another proj definition for EPSG:28992, your installation might be broken. Have a look at the setup log if postinstall was broken somewhere.

Stack Exchange Network

Return to Answer