[QGIS-it-user] Integrare i grigliati NTv2 in QGIS 3.10 (futura LTR)

Andrea Giudiceandrea andreaerdna a libero.it
Lun 6 Gen 2020 04:55:08 PST


Ciao Salvatore,
ho finalmente potuto installare QGIS 3.10.1 su Windows 7 64 bit tramite lo
standalone installer QGIS-OSGeo4W-3.10.1-1-Setup-x86_64.exe con GDAL 3.0.2
(gdal-3.0.2-3) e PROJ 6.2.1 (proj-6.2.1-2).

Dopo aver copiato il file ItalyRome40ToWGS84_NTV2_GN.gsb nella cartella
C:\Program Files\QGIS 3.10\share\proj ho eseguito la sequente query sql sul
db proj.db nella stessa cartella per fare in modo che venga utilizzato nelle
transformazioni che coinvolgono EPSG:4265 Monte Mario ed EPSG:4258 ETRS89:

INSERT INTO `grid_transformation`
VALUES ('PROJ','EPSG_4265_TO_EPSG_4258','Monte Mario (ROMA40) TO
ETRS89/ETRF89 (GN)',
NULL,NULL,
'EPSG','9615','NTv2',
'EPSG','4265','EPSG','4258',
'EPSG','1127',
NULL,
'EPSG','8656','Latitude and longitude difference file',
'ItalyRome40ToWGS84_NTV2_GN.gsb',
NULL,NULL,NULL,NULL,NULL,NULL,NULL,0);


Ho anche eseguito la seguente query, ma ho poi verificato che ciò non è
necessario (l'uso del grigliato funziona anche senza modificare la tabella
grid_alternatives, in quanto queste informazioni dovrebbero solo servire a
far capire quale pacchetto proj-datumgrid contiene il file del grigliato e
quindi quale pacchetto deve essere scaricato nel caso serva quel grigliato):

INSERT INTO `grid_alternatives`(
`original_grid_name`,`proj_grid_name`,
`proj_grid_format`,`proj_method`,`inverse_direction`,`package_name`,
`url`,`direct_download`,`open_license`,`directory`) 
VALUES (
'ItalyRome40ToWGS84_NTV2_GN.gsb','ItalyRome40ToWGS84_NTV2_GN.gsb',
'NTv2','hgridshift',0,'proj-datumgrid-europe',
NULL,NULL,NULL,NULL);


Comunque sia, il grigliato ItalyRome40ToWGS84_NTV2_GN.gsb viene utilizzato
automaticamente in QGIS 3.10.1 senza ulteriori necessità di impostare in
QGIS stesso alcuna Default Datum Trasformation (la presenza di una
trasformazione tramite grigliati NTv2 sembra annullare tutte le altre
trasformazioni considerate intrinsecamente meno precise).

Tuttavia puoi verificare che effettivamente la trasformazione sia presente
in QGIS, utilizzando la finestra di impostazioni "Select Datum
Transformations" (sempre in Settings->Options->CRS o
Project->Properties->CRS) per l'aggiunta di una trasformazione di datum
preferita:

impostando per esempio come Source CRS l'EPSG:3004 Monte Mario / Italy zone
2 e come Destination CRS l'EPSG:25833 ETRS89 / UTM zone 33N, visualizzo
l'unica trasformazione
"Inverse of Italy zone 2 + Monte Mario (ROMA40) TO ETRS89/ETRF89 (GN) + UTM
zone 33N"
che, nella nuova stringa proj contiene, come accennavo in un precedente
messaggio, i nuovi parametri "+proj=hgridshift
+grids=ItalyRome40ToWGS84_NTV2_GN.gsb" al posto del vecchio "+nadgrids=".


Il non funzionamento della tua query per l'inserimento nella tabella
grid_transformation penso sia molto probabilemente dovuto al fatto che hai
indicato il codice EPSG dell'area di uso della trasformazione come 3004 che
però corrisponde alla zona antartica (vedi tabella "area"), mentre dovrebbe
funzionare utilizzando il codice EPSG 1127 che corrisponde all'intera
Italia.

A presto.

Andrea



--
Sent from: http://osgeo-org.1560.x6.nabble.com/QGIS-Italian-User-f5250612.html


Maggiori informazioni sulla lista QGIS-it-user