8. Supporto alle proiezioni

Suggerimento

I frammenti di codice di questa pagina necessitano delle seguenti importazioni se sei è al di fuori della console pyqgis:

from qgis.core import (
    QgsCoordinateReferenceSystem,
    QgsCoordinateTransform,
    QgsProject,
    QgsPointXY,
)

8.1. Sistemi di riferimento delle coordinate

I sistemi di riferimento delle coordinate (SR) sono incapsulati dalla classe QgsCoordinateReferenceSystem. Le istanze di questa classe possono essere create in diversi modi:

• Specificare SR con il suo ID

  # EPSG 4326 is allocated for WGS84
  crs = QgsCoordinateReferenceSystem("EPSG:4326")
  print(crs.isValid())
  

  True
  

  QGIS supporta diversi identificatori SR con i seguenti formati:

  • EPSG:<code> — ID assegnato dall’organizzazione EPSG - gestito con createFromOgcWms()

  • POSTGIS:<srid>— ID utilizzato nei database PostGIS - gestito con createFromSrid()

  • INTERNAL:<srsid> — ID utilizzato nel database interno di QGIS - gestito con createFromSrsId()

  • PROJ:<proj> - gestito con createFromProj()

  • WKT:<wkt> - gestito con createFromWkt()

Se non viene specificato alcun prefisso, viene assunta la definizione WKT.

• specifica il SR con il suo well-known text (WKT)

  wkt = 'GEOGCS["WGS84", DATUM["WGS84", SPHEROID["WGS84", 6378137.0, 298.257223563]],' \
        'PRIMEM["Greenwich", 0.0], UNIT["degree",0.017453292519943295],' \
        'AXIS["Longitude",EAST], AXIS["Latitude",NORTH]]'
  crs = QgsCoordinateReferenceSystem(wkt)
  print(crs.isValid())
  

  True
  

• crea un SR non valido e poi utilizzare una delle funzioni create* per inizializzarlo. Nell’esempio seguente si utilizza una stringa Proj per inizializzare la proiezione.

  crs = QgsCoordinateReferenceSystem()
  crs.createFromProj("+proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs")
  print(crs.isValid())
  

  True
  

È opportuno verificare se la creazione (cioè la ricerca nel database) del SR è avvenuta con successo: isValid() deve restituire True.

Si noti che per l’inizializzazione dei sistemi di riferimento spaziali QGIS deve cercare i valori appropriati nel suo database interno srs.db. Pertanto, se creai un’applicazione indipendente, devi impostare correttamente i percorsi con QgsApplication.setPrefixPath(), altrimenti non riuscirà a trovare il database. Se esegui i comandi dalla console QGIS Python o sviluppi un plugin, non è importante: tutto è già stato impostato.

Accesso alle informazioni del sistema di riferimento spaziale:

crs = QgsCoordinateReferenceSystem("EPSG:4326")

print("QGIS CRS ID:", crs.srsid())
print("PostGIS SRID:", crs.postgisSrid())
print("Description:", crs.description())
print("Projection Acronym:", crs.projectionAcronym())
print("Ellipsoid Acronym:", crs.ellipsoidAcronym())
print("Proj String:", crs.toProj())
# check whether it's geographic or projected coordinate system
print("Is geographic:", crs.isGeographic())
# check type of map units in this CRS (values defined in QGis::units enum)
print("Map units:", crs.mapUnits())

In uscita:

QGIS CRS ID: 3452
PostGIS SRID: 4326
Description: WGS 84
Projection Acronym: longlat
Ellipsoid Acronym: EPSG:7030
Proj String: +proj=longlat +datum=WGS84 +no_defs
Is geographic: True
Map units: 6

8.2. Trasformazione SR

Puoi effettuare trasformazioni tra sistemi di riferimento spaziali diversi utilizzando la classe QgsCoordinateTransform. Il modo più semplice per usarla è creare un SR di origine e uno di destinazione e costruire un’istanza QgsCoordinateTransform con essi e il progetto corrente. Quindi basta chiamare ripetutamente la funzione transform() per eseguire la trasformazione. Per impostazione predefinita esegue la trasformazione in avanti, ma è in grado di eseguire anche la trasformazione inversa.

crsSrc = QgsCoordinateReferenceSystem("EPSG:4326")    # WGS 84
crsDest = QgsCoordinateReferenceSystem("EPSG:32633")  # WGS 84 / UTM zone 33N
transformContext = QgsProject.instance().transformContext()
xform = QgsCoordinateTransform(crsSrc, crsDest, transformContext)

# forward transformation: src -> dest
pt1 = xform.transform(QgsPointXY(18,5))
print("Transformed point:", pt1)

# inverse transformation: dest -> src
pt2 = xform.transform(pt1, QgsCoordinateTransform.ReverseTransform)
print("Transformed back:", pt2)

In uscita:

Transformed point: <QgsPointXY: POINT(832713.79873844375833869 553423.98688333143945783)>
Transformed back: <QgsPointXY: POINT(18 4.99999999999999911)>