[Qgis-developer] Datum Transformation - parameters for mainland Portugal

Pedro Venâncio pedrongvenancio at yahoo.com
Wed Jun 18 06:13:39 PDT 2014


Hi all,

I have an update to the parameters that are in Datum Transformation for mainland Portugal.

After an analysis of the parameters that are currently inserted in srs.db, I think we should keep only the latest parameters that are provided by the Portuguese Surveying Authority (Directorate-General of the Territory), to avoid confusion among users.

So, we should keep only the transformation parameters of Molodensky, Bursa-Wolf, and two grids NTv2 that are in use here in Portugal, one from the Directorate-General of the Territory, and another from the University of Porto.

I already have the SQL queries with the changes that I suggest:


delete from tbl_datum_transform
where source_crs_code = 4207 or source_crs_code = 4274;


insert into tbl_datum_transform (epsg_nr,coord_op_code,source_crs_code, target_crs_code,coord_op_method_code,p1,p2,p3,p4,p5,p6,p7,remarks,preferred,deprecated,area_of_use_code)
values (30000,100014,4803,4326,9606,-283.088,-70.693,117.445,-1.157,0.059,-0.652,-4.058,'Parâmetros de Transformação de Bursa-Wolf do Datum Lisboa para PT-TM06-ETRS89 (Direção-Geral do Território, 2011)',0,0,1294);

insert into tbl_datum_transform (epsg_nr,coord_op_code,source_crs_code, target_crs_code,coord_op_method_code,p1,p2,p3,remarks,preferred,deprecated,area_of_use_code)
values (30001,100015,4803,4326,9603,-303.861,-60.693,103.607,'Parâmetros de Transformação de Molodensky do Datum Lisboa para PT-TM06-ETRS89 (Direção-Geral do Território, 2011)',0,0,1294);

insert into tbl_datum_transform (epsg_nr,coord_op_code,source_crs_code, target_crs_code,coord_op_method_code,p1,remarks,preferred,deprecated,area_of_use_code)
values (30002,100016,4803,4258,9615,'ptLX_e89.gsb','Grelhas no formato NTv2 para a Transformação de Coordenadas do Sistema HGDLx para o Sistema PT-TM06/ETRS89 (José Alberto Gonçalves - FCUP)',1,0,1294);

insert into tbl_datum_transform (epsg_nr,coord_op_code,source_crs_code, target_crs_code,coord_op_method_code,p1,remarks,preferred,deprecated,area_of_use_code)
values (30003,100017,4803,4258,9615,'DLX_ETRS89_geo.gsb','Grelhas no formato NTv2 para a Transformação de Coordenadas do Sistema HGDLx para o Sistema PT-TM06/ETRS89 (Direção-Geral do Território)',0,0,1294);


insert into tbl_datum_transform (epsg_nr,coord_op_code,source_crs_code, target_crs_code,coord_op_method_code,p1,p2,p3,p4,p5,p6,p7,remarks,preferred,deprecated,area_of_use_code)
values (30004,100018,4274,4326,9606,-230.994,102.591,25.199,0.633,-0.239,0.900,1.950,'Parâmetros de Transformação de Bursa-Wolf do Datum 73 para PT-TM06-ETRS89 (Direção-Geral do Território, 2011)',0,0,1294);

insert into tbl_datum_transform (epsg_nr,coord_op_code,source_crs_code, target_crs_code,coord_op_method_code,p1,p2,p3,remarks,preferred,deprecated,area_of_use_code)
values (30005,100019,4274,4326,9603,-223.150,110.132,36.711,'Parâmetros de Transformação de Molodensky do Datum 73 para PT-TM06-ETRS89 (Direção-Geral do Território, 2011)',0,0,1294);

insert into tbl_datum_transform (epsg_nr,coord_op_code,source_crs_code, target_crs_code,coord_op_method_code,p1,remarks,preferred,deprecated,area_of_use_code)
values (30006,100020,4274,4258,9615,'pt73_e89.gsb','Grelhas no formato NTv2 para a Transformação de Coordenadas do Sistema HGD73 para o Sistema PT-TM06/ETRS89 (José Alberto Gonçalves - FCUP)',1,0,1294);

insert into tbl_datum_transform (epsg_nr,coord_op_code,source_crs_code, target_crs_code,coord_op_method_code,p1,remarks,preferred,deprecated,area_of_use_code)
values (30007,100021,4274,4258,9615,'D73_ETRS89_geo.gsb','Grelhas no formato NTv2 para a Transformação de Coordenadas do Sistema HGD73 para o Sistema PT-TM06/ETRS89 (Direção-Geral do Território)',0,0,1294);


Regarding the queries above, there are some codes (epsg_nr, coord_op_code) that I assigned with a sequential number, because I did not know what I should put. In coord_op_method_code, I put what seemed right to me, but please confirm.



Regarding NTv2 grids (http://goo.gl/jDhxJB) I have another doubt here. On Linux, these files are placed in /usr/share/proj. On Windows, using the OSGeo4W installer, they are placed in ..\OSGeo4W64\share\proj. With QGIS standalone, they are placed at ..\Program Files\QGIS Valmiera\share\proj.
There is no possibility of copy the grid files, automatically during the installation? This will make it easier for less experienced users.



However I came across something quite odd. In the Datum Transformation window, automatically arise some strange transformations (red in the image http://goo.gl/iIWOuw). These transformations give huge errors, and I do not understand what they represent.

May it be a bug?

The transformation operation in question (epsg_nr 1571) uses as source_crs_code the 4258 (Geodetic CRS ETRS89), which is used in some target_crs_code transformations for Portugal. May it be an incorrect match between source and target codes? It is because this transformation, and these parameters, makes no sense to me.

Thank you very much!

Best regards,
Pedro
-------------- next part --------------
An embedded and charset-unspecified text was scrubbed...
Name: qgis_datum_transformation_portugal.sql
URL: <http://lists.osgeo.org/pipermail/qgis-developer/attachments/20140618/22c5a7b7/attachment.ksh>


More information about the Qgis-developer mailing list