Viktigt

Översättning är en gemenskapsinsats du kan gå med i. Den här sidan är för närvarande översatt till 100.00%.

8. Stöd för projektioner

Råd

Kodsnuttarna på den här sidan behöver följande import om du befinner dig utanför pyqgis-konsolen:

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

8.1. Referenssystem för koordinater

Koordinatreferenssystem (CRS) är inkapslade av QgsCoordinateReferenceSystem-klassen. Instanser av denna klass kan skapas på flera olika sätt:

• ange CRS med dess ID

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

  True
  

  QGIS stöder olika CRS-identifierare med följande format:

  • EPSG:<code> — ID tilldelat av EPSG-organisationen - hanteras med createFromOgcWms()

  • POSTGIS:<srid>— ID som används i PostGIS-databaser - hanteras med createFromSrid()

  • INTERNAL:<srsid> — ID som används i den interna QGIS-databasen - hanteras med createFromSrsId()

  • PROJ:<proj> - hanterad med createFromProj()

  • WKT:<wkt> - hanterad med createFromWkt()

Om inget prefix anges antas WKT-definitionen.

• ange CRS genom dess välkända 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
  

• skapa en ogiltig CRS och sedan använda en av funktionerna create* för att initialisera den. I följande exempel använder vi en Proj-sträng för att initiera projektionen.

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

  True
  

Det är klokt att kontrollera om skapandet (d.v.s. uppslagningen i databasen) av CRS har lyckats: isValid() måste returnera True.

Observera att för initialisering av spatiala referenssystem måste QGIS leta upp lämpliga värden i sin interna databas srs.db. Om du skapar en oberoende applikation måste du därför ställa in sökvägar korrekt med QgsApplication.setPrefixPath(), annars kommer den inte att hitta databasen. Om du kör kommandona från QGIS Python-konsolen eller utvecklar ett tillägg behöver du inte bry dig: allt är redan förberett åt dig.

Tillgång till information om spatialt referenssystem:

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())

Utdata:

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. CRS-omvandling

Du kan göra transformationer mellan olika spatiala referenssystem genom att använda klassen QgsCoordinateTransform. Det enklaste sättet att använda den är att skapa ett CRS för källa och destination och konstruera en QgsCoordinateTransform-instans med dem och det aktuella projektet. Sedan anropar du bara upprepade gånger transform()-funktionen för att göra transformationen. Som standard gör den framåtriktad transformation, men den kan också göra invers transformation.

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)

Utdata:

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