GDAL APILink to this heading#

ExempeldataLink to this heading#

GDAL/OGR-verktygen som beskrivs här är utformade för att hjälpa dig att läsa in dina geospatiala data, för att de flesta av dem ska vara användbara måste du ha några data att arbeta med. Om du börjar och ännu inte har några egna data att använda, innehåller GeoDjango-tester ett antal datauppsättningar som du kan använda för testning. Du kan ladda ner dem här:

$ wget https://raw.githubusercontent.com/django/django/main/tests/gis_tests/data/cities/cities.{shp,prj,shx,dbf}
$ wget https://raw.githubusercontent.com/django/django/main/tests/gis_tests/data/rasters/raster.tif

DataSourceLink to this heading#

DataSource är ett omslag för OGR:s datakällsobjekt som stöder läsning av data från en mängd olika OGR-stödda geospatiala filformat och datakällor med hjälp av ett konsekvent gränssnitt. Varje datakälla representeras av ett DataSource-objekt som innehåller ett eller flera datalager. Varje skikt, som representeras av ett Layer-objekt, innehåller ett antal geografiska egenskaper (Feature), information om vilken typ av egenskaper som finns i skiktet (t.ex. punkter, polygoner etc.) samt namn och typ för eventuella ytterligare datafält (Field) som kan associeras med varje egenskap i skiktet.

class DataSource(ds_input, encoding='utf-8')Link to this definition#

Konstruktören för DataSource kräver bara en parameter: sökvägen till den fil du vill läsa. OGR stöder dock också en mängd mer komplexa datakällor, inklusive databaser, som kan nås genom att skicka en speciell namnsträng istället för en sökväg. För mer information, se dokumentationen `OGR Vector Formats`__. Egenskapen name för en DataSource-instans ger OGR-namnet på den underliggande datakällan som den använder.

Med den valfria parametern encoding kan du ange en icke-standardiserad kodning av strängarna i källan. Detta är vanligtvis användbart när du får DjangoUnicodeDecodeError-undantag när du läser fältvärden.

När du har skapat din DataSource kan du ta reda på hur många datalager den innehåller genom att använda egenskapen layer_count eller (på motsvarande sätt) genom att använda funktionen len(). Information om hur du kommer åt själva datalagren finns i nästa avsnitt:

Python console

>>> from django.contrib.gis.gdal import DataSource
>>> ds = DataSource("/path/to/your/cities.shp")
>>> ds.name
'/path/to/your/cities.shp'
>>> ds.layer_count  # This file only contains one layer
1

layer_countLink to this definition#

Returnerar antalet lager i datakällan.

nameLink to this definition#

Returnerar namnet på datakällan.

LagerLink to this heading#

class LayerLink to this definition#

Layer är ett omslag för ett lager av data i ett DataSource-objekt. Du skapar aldrig ett Layer-objekt direkt. Istället hämtar du dem från ett DataSource-objekt, som i huvudsak är en standard Python-behållare för Layer-objekt. Du kan till exempel komma åt ett specifikt lager genom dess index (t.ex. ds[0] för att komma åt det första lagret), eller så kan du iterera över alla lager i behållaren i en for-loop. Själva Layer fungerar som en behållare för geometriska egenskaper.

Vanligtvis har alla objekt i ett visst skikt samma geometrityp. Egenskapen geom_type för ett skikt är en OGRGeomType som identifierar funktionstypen. Vi kan använda den för att skriva ut grundläggande information om varje lager i en DataSource:

Python console

>>> for layer in ds:
...     print('Layer "%s": %i %ss' % (layer.name, len(layer), layer.geom_type.name))
...
Layer "cities": 3 Points

Exemplet är från datakällan cities, som laddades ovan och som uppenbarligen innehåller ett skikt, kallat "cities", som innehåller tre punktfunktioner. För enkelhetens skull antar exemplen nedan att du har lagrat det skiktet i variabeln layer:

Python console

>>> layer = ds[0]

nameLink to this definition#

Returnerar namnet på detta lager i datakällan.

Python console

>>> layer.name
'cities'

num_featLink to this definition#

Returnerar antalet funktioner i lagret. Samma som len(layer):

Python console

>>> layer.num_feat
3

geom_typeLink to this definition#

Returnerar geometritypen för skiktet, som ett OGRGeomType-objekt:

Python console

>>> layer.geom_type.name
'Point'

num_fieldsLink to this definition#

Returnerar antalet fält i skiktet, dvs. antalet fält med data som är associerade med varje funktion i skiktet:

Python console

>>> layer.num_fields
4

fieldsLink to this definition#

Returnerar en lista med namnen på alla fält i detta skikt:

Python console

>>> layer.fields
['Name', 'Population', 'Density', 'Created']

Returnerar en lista med datatyperna för varje fält i detta skikt. Dessa är subklasser till Field, som beskrivs nedan:

Python console

>>> [ft.__name__ for ft in layer.field_types]
['OFTString', 'OFTReal', 'OFTReal', 'OFTDate']

field_widthsLink to this definition#

Returnerar en lista med de maximala fältbredderna för varje fält i detta skikt:

Python console

>>> layer.field_widths
[80, 11, 24, 10]

field_precisionsLink to this definition#

Returnerar en lista med de numeriska precisionsvärdena för varje fält i detta skikt. Detta är meningslöst (och sätts till noll) för icke-numeriska fält:

Python console

>>> layer.field_precisions
[0, 0, 15, 0]

extentLink to this definition#

Returnerar den rumsliga utbredningen av detta lager, som ett Envelope-objekt:

Python console

>>> layer.extent.tuple
(-104.609252, 29.763374, -95.23506, 38.971823)

srsLink to this definition#

Egenskap som returnerar den SpatialReference som är associerad med detta lager:

Python console

>>> print(layer.srs)
GEOGCS["GCS_WGS_1984",
    DATUM["WGS_1984",
        SPHEROID["WGS_1984",6378137,298.257223563]],
    PRIMEM["Greenwich",0],
    UNIT["Degree",0.017453292519943295]]

Om Layer inte har någon spatial referensinformation associerad med sig, returneras None.

spatial_filterLink to this definition#

Egenskap som kan användas för att hämta eller ställa in ett rumsligt filter för detta skikt. Ett rumsligt filter kan endast ställas in med en OGRGeometry-instans, en 4-tupel extent eller None. När det anges med något annat än None kommer endast objekt som skär filtret att returneras vid iterering över skiktet:

Python console

>>> print(layer.spatial_filter)
None
>>> print(len(layer))
3
>>> [feat.get("Name") for feat in layer]
['Pueblo', 'Lawrence', 'Houston']
>>> ks_extent = (-102.051, 36.99, -94.59, 40.00)  # Extent for state of Kansas
>>> layer.spatial_filter = ks_extent
>>> len(layer)
1
>>> [feat.get("Name") for feat in layer]
['Lawrence']
>>> layer.spatial_filter = None
>>> len(layer)
3

get_fields()Link to this definition#

En metod som returnerar en lista med värdena för ett visst fält för varje objekt i skiktet:

Python console

>>> layer.get_fields("Name")
['Pueblo', 'Lawrence', 'Houston']

get_geoms(geos=False)Link to this definition#

En metod som returnerar en lista som innehåller geometrin för varje objekt i skiktet. Om det valfria argumentet geos är satt till True konverteras geometrierna till GEOSGeometry-objekt. Annars returneras de som OGRGeometry-objekt:

Python console

>>> [pt.tuple for pt in layer.get_geoms()]
[(-104.609252, 38.255001), (-95.23506, 38.971823), (-95.363151, 29.763374)]

test_capability(capability)Link to this definition#

Returnerar ett booleanskt värde som anger om lagret stöder den angivna egenskapen (en sträng). Exempel på giltiga kapacitetssträngar är: 'RandomRead', 'SequentialWrite', 'RandomWrite', 'FastSpatialFilter', 'FastFeatureCount', 'FastGetExtent', 'CreateField', 'Transactions', 'DeleteFeature' och 'FastSetNextByIndex'.

FunktionLink to this heading#

class FeatureLink to this definition#

Feature omsluter en OGR-funktion. Du skapar aldrig ett Feature-objekt direkt. Istället hämtar du dem från ett Layer-objekt. Varje feature består av en geometri och en uppsättning fält som innehåller ytterligare egenskaper. Geometrin för ett fält är tillgänglig via dess egenskap geom, som returnerar ett OGRGeometry-objekt. En Feature beter sig som en standard Python-behållare för sina fält, som den returnerar som Field-objekt: du kan komma åt ett fält direkt via dess index eller namn, eller så kan du iterera över en features fält, t.ex. i en for-loop.

geomLink to this definition#

Returnerar geometrin för den här funktionen som ett OGRGeometry-objekt:

Python console

>>> city.geom.tuple
(-104.609252, 38.255001)

getLink to this definition#

En metod som returnerar värdet för det angivna fältet (specificerat med namn) för denna funktion, inte ett Field omslagsobjekt:

Python console

>>> city.get("Population")
102121

geom_typeLink to this definition#

Returnerar geometritypen för den här funktionen, som ett OGRGeomType-objekt. Detta kommer att vara detsamma för alla objekt i ett visst skikt och motsvarar egenskapen Layer.geom_type i objektet Layer som objektet kommer från.

num_fieldsLink to this definition#

Returnerar antalet datafält som är associerade med objektet. Detta kommer att vara detsamma för alla funktioner i ett visst skikt och motsvarar egenskapen Layer.num_fields i objektet Layer som funktionen kommer från.

fieldsLink to this definition#

Returnerar en lista med namnen på de datafält som är associerade med objektet. Detta kommer att vara detsamma för alla objekt i ett visst skikt och motsvarar egenskapen Layer.fields i objektet Layer som objektet kommer från.

fidLink to this definition#

Returnerar identifieraren för objektet i skiktet:

Python console

>>> city.fid
0

layer_nameLink to this definition#

Returnerar namnet på Layer som objektet kommer från. Detta kommer att vara detsamma för alla funktioner i ett visst lager:

Python console

>>> city.layer_name
'cities'

indexLink to this definition#

En metod som returnerar indexet för det angivna fältnamnet. Detta kommer att vara detsamma för alla objekt i ett visst lager:

Python console

>>> city.index("Population")
1

FältLink to this heading#

class FieldLink to this definition#

nameLink to this definition#

Returnerar namnet på detta fält:

Python console

>>> city["Name"].name
'Name'

typeLink to this definition#

Returnerar OGR-typen för detta fält, som ett heltal. Ordboken FIELD_CLASSES mappar dessa värden till subklasser av Field:

Python console

>>> city["Density"].type
2

type_nameLink to this definition#

Returnerar en sträng med namnet på datatypen för detta fält:

Python console

>>> city["Name"].type_name
'String'

valueLink to this definition#

Returnerar värdet för detta fält. Klassen Field returnerar själv värdet som en sträng, men varje underklass returnerar värdet i den mest lämpliga formen:

Python console

>>> city["Population"].value
102121

widthLink to this definition#

Returnerar bredden på detta fält:

Python console

>>> city["Name"].width
80

precisionLink to this definition#

Returnerar den numeriska precisionen för detta fält. Detta är meningslöst (och sätts till noll) för icke-numeriska fält:

Python console

>>> city["Density"].precision
15

as_double()Link to this definition#

Returnerar fältets värde som en dubbel (float):

Python console

>>> city["Density"].as_double()
874.7

as_int()Link to this definition#

Returnerar fältets värde som ett heltal:

Python console

>>> city["Population"].as_int()
102121

as_string()Link to this definition#

Returnerar fältets värde som en sträng:

Python console

>>> city["Name"].as_string()
'Pueblo'

as_datetime()Link to this definition#

Returnerar fältets värde som en tupel av datum- och tidskomponenter:

Python console

>>> city["Created"].as_datetime()
(c_long(1999), c_long(5), c_long(23), c_long(0), c_long(0), c_long(0), c_long(0))

DrivrutinLink to this heading#

class Driver(dr_input)Link to this definition#

Klassen Driver används internt för att omsluta en OGR DataSource-drivrutin.

driver_countLink to this definition#

Returnerar antalet OGR-vektorförare som för närvarande är registrerade.

OGRGeometriLink to this heading#

OGRGeometry-objekt delar liknande funktionalitet med GEOSGeometry-objekt och är tunna omslag runt OGR:s interna geometrirepresentation. De möjliggör således effektivare åtkomst till data när man använder DataSource. Till skillnad från sin GEOS-motsvarighet har OGRGeometry stöd för rumsliga referenssystem och koordinattransformation:

>>> from django.contrib.gis.gdal import OGRGeometry
>>> polygon = OGRGeometry("POLYGON((0 0, 5 0, 5 5, 0 5))")

class OGRGeometry(geom_input, srs=None)Link to this definition#

Detta objekt är en omslagsform för klassen `OGR Geometry`__. Dessa objekt instansieras direkt från den givna parametern geom_input, som kan vara en sträng som innehåller WKT, HEX, GeoJSON, en buffer som innehåller WKB-data, eller ett OGRGeomType-objekt. Dessa objekt returneras också från attributet Feature.geom när vektordata läses från Layer (som i sin tur är en del av en DataSource).

classmethod from_gml(gml_string)Link to this definition#

Konstruerar en OGRGeometry från den givna GML-strängen.

classmethod from_bbox(bbox)Link to this definition#

Konstruerar en Polygon från den givna bounding-boxen (en 4-tupel).

__len__()Link to this definition#

Returnerar antalet punkter i en LineString, antalet ringar i en Polygon eller antalet geometrier i en GeometryCollection. Gäller inte för andra geometrityper.

__iter__()Link to this definition#

Itererar över punkterna i en LineString, ringarna i en Polygon eller geometrierna i en GeometryCollection. Gäller inte för andra geometrityper.

__getitem__()Link to this definition#

Returnerar punkten vid det angivna indexet för en LineString, den inre ringen vid det angivna indexet för en Polygon, eller geometrin vid det angivna indexet i en GeometryCollection. Gäller inte för andra geometrityper.

dimensionLink to this definition#

Returnerar antalet koordinerade dimensioner i geometrin, dvs. 0 för punkter, 1 för linjer osv:

Python console

>>> polygon.dimension
2

coord_dimLink to this definition#

Returnerar koordinatdimensionen för denna geometri. Värdet är t.ex. 2 för tvådimensionella geometrier.

is_3dLink to this definition#

En boolean som anger om geometrin har Z-koordinater.

set_3d(value)Link to this definition#

En metod som lägger till eller tar bort Z-koordinatdimensionen.

Python console

>>> p = OGRGeometry("POINT (1 2 3)")
>>> p.is_3d
True
>>> p.set_3d(False)
>>> p.wkt
"POINT (1 2)"

is_measuredLink to this definition#

En boolean som anger om geometrin har M-koordinater.

set_measured(value)Link to this definition#

En metod för att lägga till eller ta bort M-koordinatdimensionen.

Python console

>>> p = OGRGeometry("POINT (1 2)")
>>> p.is_measured
False
>>> p.set_measured(True)
>>> p.wkt
"POINT M (1 2 0)"

geom_countLink to this definition#

Returnerar antalet element i denna geometri:

Python console

>>> polygon.geom_count
1

has_curveLink to this definition#

En boolean som anger om geometrin är eller innehåller en kurvgeometri.

get_linear_geometry()Link to this definition#

Returnerar en linjär version av geometrin. Om ingen konvertering kan göras returneras originalgeometrin.

get_curve_geometry()Link to this definition#

Returnerar en krökt version av geometrin. Om ingen konvertering kan göras returneras originalgeometrin.

point_countLink to this definition#

Returnerar antalet punkter som används för att beskriva denna geometri:

Python console

>>> polygon.point_count
4

num_pointsLink to this definition#

Alias för point_count.

num_coordsLink to this definition#

Alias för point_count.

geom_typeLink to this definition#

Returnerar typen av denna geometri, som ett OGRGeomType-objekt.

geom_nameLink to this definition#

Returnerar namnet på typen för denna geometri:

Python console

>>> polygon.geom_name
'POLYGON'

areaLink to this definition#

Returnerar arean för denna geometri, eller 0 för geometrier som inte innehåller någon area:

Python console

>>> polygon.area
25.0

envelopeLink to this definition#

Returnerar kuvertet för denna geometri, som ett Envelope-objekt.

extentLink to this definition#

Returnerar kuvertet för denna geometri som en 4-tupel, istället för som ett Envelope-objekt:

Python console

>>> point.extent
(0.0, 0.0, 5.0, 5.0)

srsLink to this definition#

Denna egenskap styr den rumsliga referensen för denna geometri, eller None om inget rumsligt referenssystem har tilldelats den. Om den är tilldelad returnerar åtkomst till denna egenskap ett SpatialReference-objekt. Den kan ställas in med ett annat SpatialReference-objekt, eller någon annan indata som SpatialReference accepterar. Exempel på detta:

Python console

>>> city.geom.srs.name
'GCS_WGS_1984'

sridLink to this definition#

Returnerar eller anger den identifierare för rumslig referens som motsvarar SpatialReference för denna geometri. Returnerar None om det inte finns någon rumslig referensinformation associerad med denna geometri, eller om en SRID inte kan fastställas.

geosLink to this definition#

Returnerar ett GEOSGeometry-objekt som motsvarar denna geometri.

gmlLink to this definition#

Returnerar en strängrepresentation av denna geometri i GML-format:

Python console

>>> OGRGeometry("POINT(1 2)").gml
'<gml:Point><gml:coordinates>1,2</gml:coordinates></gml:Point>'

hexLink to this definition#

Returnerar en strängrepresentation av denna geometri i HEX WKB-format:

Python console

>>> OGRGeometry("POINT(1 2)").hex
'0101000000000000000000F03F0000000000000040'

jsonLink to this definition#

Returnerar en strängrepresentation av denna geometri i JSON-format:

Python console

>>> OGRGeometry("POINT(1 2)").json
'{ "type": "Point", "coordinates": [ 1.000000, 2.000000 ] }'

kmlLink to this definition#

Returnerar en strängrepresentation av denna geometri i KML-format.

wkb_sizeLink to this definition#

Returnerar storleken på den WKB-buffert som behövs för att hålla en WKB-representation av denna geometri:

Python console

>>> OGRGeometry("POINT(1 2)").wkb_size
21

wkbLink to this definition#

Returnerar en buffer som innehåller en WKB-representation av denna geometri.

wktLink to this definition#

Returnerar en strängrepresentation av denna geometri i WKT-format.

ewktLink to this definition#

Returnerar EWKT-representationen av denna geometri.

clone()Link to this definition#

Returnerar en ny OGRGeometry-klon av detta geometriobjekt.

close_rings()Link to this definition#

Om det finns några ringar inom denna geometri som inte har stängts, kommer denna rutin att göra det genom att lägga till startpunkten till slutet:

Python console

>>> triangle = OGRGeometry("LINEARRING (0 0,0 1,1 0)")
>>> triangle.close_rings()
>>> triangle.wkt
'LINEARRING (0 0,0 1,1 0,0 0)'

transform(coord_trans, clone=False)Link to this definition#

Transformerar denna geometri till ett annat rumsligt referenssystem. Kan ta ett CoordTransform-objekt, ett SpatialReference-objekt, eller någon annan indata som accepteras av SpatialReference (inklusive rumsliga referens WKT- och PROJ-strängar, eller ett heltal SRID).

Som standard returneras ingenting och geometrin transformeras på plats. Men om nyckelordet clone är satt till True returneras istället en transformerad klon av denna geometri.

intersects(other)Link to this definition#

Returnerar True om denna geometri skär den andra, annars returneras False.

equals(other)Link to this definition#

Returnerar True om denna geometri är likvärdig med den andra, annars returneras False.

disjoint(other)Link to this definition#

Returnerar True om denna geometri är rumsligt disjunkt med (dvs. inte skär) den andra, annars returneras False.

touches(other)Link to this definition#

Returnerar True om denna geometri berör den andra, annars returneras False.

crosses(other)Link to this definition#

Returnerar True om denna geometri korsar den andra, annars returneras False.

within(other)Link to this definition#

Returnerar True om denna geometri är innesluten i den andra, annars returneras False.

contains(other)Link to this definition#

Returnerar True om denna geometri innehåller den andra, annars returneras False.

overlaps(other)Link to this definition#

Returnerar True om denna geometri överlappar den andra, annars returneras False.

boundary()Link to this definition#

Gränsen för denna geometri, som ett nytt OGRGeometry-objekt.

convex_hullLink to this definition#

Den minsta konvexa polygonen som innehåller denna geometri, som ett nytt OGRGeometry-objekt.

difference()Link to this definition#

Returnerar regionen som består av skillnaden mellan den här geometrin och den andra, som ett nytt OGRGeometry-objekt.

intersection()Link to this definition#

Returnerar regionen som består av skärningspunkten mellan denna geometri och den andra, som ett nytt OGRGeometry-objekt.

sym_difference()Link to this definition#

Returnerar den region som består av den symmetriska skillnaden mellan denna geometri och den andra, som ett nytt OGRGeometry-objekt.

union()Link to this definition#

Returnerar den region som består av föreningen av denna geometri och den andra, som ett nytt OGRGeometry-objekt.

centroidLink to this definition#

Returnerar en Point som representerar geometrins centroid.

tupleLink to this definition#

Returnerar koordinaterna för en punktgeometri som en tupel, koordinaterna för en linjegeometri som en tupel av tuplar osv:

Python console

>>> OGRGeometry("POINT (1 2)").tuple
(1.0, 2.0)
>>> OGRGeometry("LINESTRING (1 2,3 4)").tuple
((1.0, 2.0), (3.0, 4.0))

coordsLink to this definition#

Ett alias för tuple.

class PointLink to this definition#

xLink to this definition#

Returnerar X-koordinaten för denna punkt:

Python console

>>> OGRGeometry("POINT (1 2)").x
1.0

yLink to this definition#

Returnerar Y-koordinaten för denna punkt:

Python console

>>> OGRGeometry("POINT (1 2)").y
2.0

zLink to this definition#

Returnerar Z-koordinaten för denna punkt, eller None om punkten inte har någon Z-koordinat:

Python console

>>> OGRGeometry("POINT (1 2 3)").z
3.0

mLink to this definition#

Returnerar M-koordinaten för denna punkt, eller None om Point inte har någon M-koordinat:

Python console

>>> OGRGeometry("POINT ZM (1 2 3 4)").m
4.0

class LineStringLink to this definition#

xLink to this definition#

Returnerar en lista med X-koordinater i denna linje:

Python console

>>> OGRGeometry("LINESTRING (1 2,3 4)").x
[1.0, 3.0]

yLink to this definition#

Returnerar en lista med Y-koordinater i denna linje:

Python console

>>> OGRGeometry("LINESTRING (1 2,3 4)").y
[2.0, 4.0]

zLink to this definition#

Returnerar en lista med Z-koordinater i linjen, eller None om linjen inte har några Z-koordinater:

Python console

>>> OGRGeometry("LINESTRING (1 2 3,4 5 6)").z
[3.0, 6.0]

mLink to this definition#

Returnerar en lista med M-koordinater i den här linjen eller None om linjen inte har M-koordinater:

Python console

>>> OGRGeometry("LINESTRING(0 1 2 10, 1 2 3 11, 2 3 4 12)").m
[10.0, 11.0, 12.0]

class PolygonLink to this definition#

shellLink to this definition#

Returnerar skalet eller den yttre ringen för denna polygon, som en LinearRing geometri.

exterior_ringLink to this definition#

Ett alias för shell.

class GeometryCollectionLink to this definition#

add(geom)Link to this definition#

Lägger till en geometri till denna geometrisamling. Gäller inte för andra geometrityper.

OGRGeomTypeLink to this heading#

class OGRGeomType(type_input)Link to this definition#

Denna klass gör det möjligt att representera en OGR-geometrityp på flera olika sätt:

Python console

>>> from django.contrib.gis.gdal import OGRGeomType
>>> gt1 = OGRGeomType(3)  # Using an integer for the type
>>> gt2 = OGRGeomType("Polygon")  # Using a string
>>> gt3 = OGRGeomType("POLYGON")  # It's case-insensitive
>>> print(gt1 == 3, gt1 == "Polygon")  # Equivalence works w/non-OGRGeomType objects
True True

nameLink to this definition#

Returnerar en kortfattad strängform av OGR Geometry-typen:

Python console

>>> gt1.name
'Polygon'

numLink to this definition#

Returnerar det nummer som motsvarar OGR-geometritypen:

Python console

>>> gt1.num
3

djangoLink to this definition#

Returnerar Django-fälttypen (en underklass till GeometryField) som ska användas för att lagra denna OGR-typ, eller None om det inte finns någon lämplig Django-typ:

Python console

>>> gt1.django
'PolygonField'

KuvertLink to this heading#

class Envelope(*args)Link to this definition#

Representerar en OGR Envelope-struktur som innehåller de lägsta och högsta X- och Y-koordinaterna för en rektangels begränsningsbox. Namngivningen av variablerna är kompatibel med OGR Envelope C-strukturen.

min_xLink to this definition#

Värdet på den minsta X-koordinaten.

min_yLink to this definition#

Värdet på den maximala X-koordinaten.

max_xLink to this definition#

Värdet på den minsta Y-koordinaten.

max_yLink to this definition#

Värdet på den maximala Y-koordinaten.

urLink to this definition#

Den övre högra koordinaten, som en tupel.

llLink to this definition#

Den nedre vänstra koordinaten, som en tupel.

tupleLink to this definition#

En tupel som representerar kuvertet.

wktLink to this definition#

En sträng som representerar detta kuvert som en polygon i WKT-format.

expand_to_include(*args)Link to this definition#

SpatialReferenceLink to this heading#

class SpatialReference(srs_input)Link to this definition#

Rumsliga referensobjekt initialiseras på den givna srs_input, som kan vara en av följande:

• OGC Välkänd text (WKT) (en sträng)

• EPSG-kod (heltal eller sträng)

• PROJ sträng

• En kortfattad sträng för välkända standarder ('WGS84', 'WGS72', 'NAD27', 'NAD83')

Exempel:

Python console

>>> wgs84 = SpatialReference("WGS84")  # shorthand string
>>> wgs84 = SpatialReference(4326)  # EPSG code
>>> wgs84 = SpatialReference("EPSG:4326")  # EPSG string
>>> proj = "+proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs "
>>> wgs84 = SpatialReference(proj)  # PROJ string
>>> wgs84 = SpatialReference("""GEOGCS["WGS 84",
... DATUM["WGS_1984",
...      SPHEROID["WGS 84",6378137,298.257223563,
...          AUTHORITY["EPSG","7030"]],
...      AUTHORITY["EPSG","6326"]],
... PRIMEM["Greenwich",0,
...      AUTHORITY["EPSG","8901"]],
... UNIT["degree",0.01745329251994328,
...      AUTHORITY["EPSG","9122"]],
... AUTHORITY["EPSG","4326"]]""")  # OGC WKT

__getitem__(target)Link to this definition#

Returnerar värdet för den angivna strängattributnoden, None om noden inte finns. Kan också ta en tupel som parameter, (target, child), där child är attributets index i WKT. Till exempel: (target)

Python console

>>> wkt = 'GEOGCS["WGS 84", DATUM["WGS_1984, ... AUTHORITY["EPSG","4326"]]'
>>> srs = SpatialReference(wkt)  # could also use 'WGS84', or 4326
>>> print(srs["GEOGCS"])
WGS 84
>>> print(srs["DATUM"])
WGS_1984
>>> print(srs["AUTHORITY"])
EPSG
>>> print(srs["AUTHORITY", 1])  # The authority value
4326
>>> print(srs["TOWGS84", 4])  # the fourth value in this wkt
0
>>> print(srs["UNIT|AUTHORITY"])  # For the units authority, have to use the pipe symbol.
EPSG
>>> print(srs["UNIT|AUTHORITY", 1])  # The authority value for the units
9122

attr_value(target, index=0)Link to this definition#

Attributvärdet för den angivna målnoden (t.ex. 'PROJCS'). Nyckelordet index anger ett index för den underordnade nod som ska returneras.

auth_name(target)Link to this definition#

Returnerar auktoritetsnamnet för den givna strängens målnod.

auth_code(target)Link to this definition#

Returnerar auktoritetskoden för den givna strängens målnod.

clone()Link to this definition#

Returnerar en klon av detta rumsliga referensobjekt.

identify_epsg()Link to this definition#

Denna metod inspekterar WKT för denna SpatialReference och lägger till EPSG-auktoritetsnoder där en EPSG-identifierare är tillämplig.

from_esri()Link to this definition#

Omvandlar denna SpatialReference från ESRI:s format till EPSG

to_esri()Link to this definition#

Morphar denna SpatialReference till ESRI:s format.

validate()Link to this definition#

Kontrollerar om den angivna rumsliga referensen är giltig, om så inte är fallet kommer ett undantag att uppstå.

import_epsg(epsg)Link to this definition#

Importera rumslig referens från EPSG-kod.

import_proj(proj)Link to this definition#

Importera rumslig referens från PROJ-sträng.

import_user_input(user_input)Link to this definition#

import_wkt(wkt)Link to this definition#

Importera rumslig referens från WKT.

import_xml(xml)Link to this definition#

Importera rumslig referens från XML.

nameLink to this definition#

Returnerar namnet på denna rumsliga referens.

sridLink to this definition#

Returnerar SRID för högsta nivåns auktoritet, eller None om den inte är definierad.

linear_nameLink to this definition#

Returnerar namnet på de linjära enheterna.

linear_unitsLink to this definition#

Returnerar värdet för de linjära enheterna.

angular_nameLink to this definition#

Returnerar namnet på vinkelenheterna.”

angular_unitsLink to this definition#

Returnerar värdet för vinkelenheterna.

unitsLink to this definition#

Returnerar en 2-tupel av enhetsvärdet och enhetsnamnet och avgör automatiskt om det är linjära enheter eller vinkelenheter som ska returneras.

ellipsoidLink to this definition#

Returnerar en tupel av ellipsoidparametrarna för denna rumsliga referens: (semimajoraxel, semiminoraxel och invers avplaning).

semi_majorLink to this definition#

Returnerar ellipsoidens halva huvudaxel för denna rumsliga referens.

semi_minorLink to this definition#

Returnerar ellipsoidens halva mindre axel för denna rumsliga referens.

inverse_flatteningLink to this definition#

Returnerar ellipsoidens inversa tillplattning för denna rumsliga referens.

geographicLink to this definition#

Returnerar True om denna rumsliga referens är geografisk (rotnoden är GEOGCS).

localLink to this definition#

Returnerar True om denna rumsliga referens är lokal (rotnoden är LOCAL_CS).

projectedLink to this definition#

Returnerar True om denna rumsliga referens är ett projicerat koordinatsystem (rotnoden är PROJCS).

wktLink to this definition#

Returnerar WKT-representationen av denna rumsliga referens.

pretty_wktLink to this definition#

Returnerar den ”vackra” representationen av WKT.

projLink to this definition#

Returnerar PROJ-representationen för denna rumsliga referens.

proj4Link to this definition#

Alias för SpatialReference.proj.

xmlLink to this definition#

Returnerar XML-representationen av denna rumsliga referens.

CoordTransformLink to this heading#

class CoordTransform(source, target)Link to this definition#

Representerar en koordinatsystemtransformation. Den initialiseras med två SpatialReference, som representerar käll- respektive målkoordinatsystemen. Dessa objekt bör användas när samma koordinattransformation utförs upprepade gånger på olika geometrier:

>>> ct = CoordTransform(SpatialReference("WGS84"), SpatialReference("NAD83"))
>>> for feat in layer:
...     geom = feat.geom  # getting clone of feature geometry
...     geom.transform(ct)  # transforming
...

GDALRasterLink to this heading#

GDALRaster är ett omslag för GDAL:s rasterkällobjekt som stöder läsning av data från en mängd olika GDAL-stödda geospatiala filformat och datakällor med hjälp av ett konsekvent gränssnitt. Varje datakälla representeras av ett GDALRaster-objekt som innehåller ett eller flera datalager som kallas band. Varje band, som representeras av ett GDALBand-objekt, innehåller georefererade bilddata. En RGB-bild representeras t.ex. av tre band: ett för rött, ett för grönt och ett för blått.

class GDALRaster(ds_input, write=False)Link to this definition#

Konstruktören för GDALRaster accepterar två parametrar. Den första parametern definierar rasterkällan och den andra parametern definierar om ett raster ska öppnas i skrivläge. För nyskapade raster ignoreras den andra parametern och det nya rastret skapas alltid i skrivläge.

Den första parametern kan ha tre former: en sträng eller Path som representerar en filsökväg (filsystem eller GDAL:s virtuella filsystem), en ordbok med värden som definierar ett nytt raster eller ett bytesobjekt som representerar en rasterfil.

Om indata är en sökväg till en fil öppnas rastret därifrån. Om indata är rådata i en ordbok krävs parametrarna width, height och rid. Om indata är ett bytesobjekt kommer det att öppnas med hjälp av ett virtuellt GDAL-filsystem.

För en detaljerad beskrivning av hur man skapar raster med hjälp av ordboksinmatning, se Skapa raster från data. För en detaljerad beskrivning av hur man skapar raster i det virtuella filsystemet, se Använda GDAL:s virtuella filsystem.

Följande exempel visar hur raster kan skapas från olika ingångskällor (med hjälp av exempeldata från GeoDjango-testerna; se även avsnittet Exempeldata).

Python console

>>> from django.contrib.gis.gdal import GDALRaster
>>> rst = GDALRaster("/path/to/your/raster.tif", write=False)
>>> rst.name
'/path/to/your/raster.tif'
>>> rst.width, rst.height  # This file has 163 x 174 pixels
(163, 174)
>>> rst = GDALRaster(
...     {  # Creates an in-memory raster
...         "srid": 4326,
...         "width": 4,
...         "height": 4,
...         "datatype": 1,
...         "bands": [
...             {
...                 "data": (2, 3),
...                 "offset": (1, 1),
...                 "size": (2, 2),
...                 "shape": (2, 1),
...                 "nodata_value": 5,
...             }
...         ],
...     }
... )
>>> rst.srs.srid
4326
>>> rst.width, rst.height
(4, 4)
>>> rst.bands[0].data()
array([[5, 5, 5, 5],
       [5, 2, 3, 5],
       [5, 2, 3, 5],
       [5, 5, 5, 5]], dtype=uint8)
>>> rst_file = open("/path/to/your/raster.tif", "rb")
>>> rst_bytes = rst_file.read()
>>> rst = GDALRaster(rst_bytes)
>>> rst.is_vsi_based
True
>>> rst.name  # Stored in a random path in the vsimem filesystem.
'/vsimem/da300bdb-129d-49a8-b336-e410a9428dad'

nameLink to this definition#

Namnet på källan, vilket motsvarar sökvägen till indatafilen eller det namn som anges vid instansiering.

Python console

>>> GDALRaster({"width": 10, "height": 10, "name": "myraster", "srid": 4326}).name
'myraster'

driverLink to this definition#

Namnet på den GDAL-drivrutin som används för att hantera indatafilen. För GDALRaster som skapas från en fil identifieras drivrutinstypen automatiskt. När raster skapas från grunden är standardinställningen ett raster i minnet ('MEM), men det kan ändras efter behov. Använd t.ex. GTiff för en GeoTiff-fil. För en lista över filtyper, se även listan `GDAL Raster Formats`__.

Ett raster i minnet skapas med hjälp av följande exempel:

Python console

>>> GDALRaster({"width": 10, "height": 10, "srid": 4326}).driver.name
'MEM'

Ett filbaserat GeoTiff-raster skapas med hjälp av följande exempel:

Python console

>>> import tempfile
>>> rstfile = tempfile.NamedTemporaryFile(suffix=".tif")
>>> rst = GDALRaster(
...     {
...         "driver": "GTiff",
...         "name": rstfile.name,
...         "srid": 4326,
...         "width": 255,
...         "height": 255,
...         "nr_of_bands": 1,
...     }
... )
>>> rst.name
'/tmp/tmp7x9H4J.tif'           # The exact filename will be different on your computer
>>> rst.driver.name
'GTiff'

widthLink to this definition#

Källans bredd i pixlar (X-axeln).

Python console

>>> GDALRaster({"width": 10, "height": 20, "srid": 4326}).width
10

heightLink to this definition#

Källans höjd i pixlar (Y-axeln).

Python console

>>> GDALRaster({"width": 10, "height": 20, "srid": 4326}).height
20

srsLink to this definition#

Rastrets rumsliga referenssystem, som en SpatialReference-instans. SRS kan ändras genom att ställa in det till en annan SpatialReference eller genom att tillhandahålla någon indata som accepteras av SpatialReference konstruktören.

Python console

>>> rst = GDALRaster({"width": 10, "height": 20, "srid": 4326})
>>> rst.srs.srid
4326
>>> rst.srs = 3086
>>> rst.srs.srid
3086

sridLink to this definition#

Rastrets SRID (Spatial Reference System Identifier). Denna egenskap är en genväg till att hämta eller ställa in SRID genom attributet srs.

Python console

>>> rst = GDALRaster({"width": 10, "height": 20, "srid": 4326})
>>> rst.srid
4326
>>> rst.srid = 3086
>>> rst.srid
3086
>>> rst.srs.srid  # This is equivalent
3086

geotransformLink to this definition#

Den affina transformationsmatris som används för att georeferera källan, som en tupel av sex koefficienter som mappar pixel-/linjekoordinater till georefererat utrymme med hjälp av följande relation:

Code

Xgeo = GT(0) + Xpixel * GT(1) + Yline * GT(2)
Ygeo = GT(3) + Xpixel * GT(4) + Yline * GT(5)

Samma värden kan hämtas genom att använda egenskaperna origin (index 0 och 3), scale (index 1 och 5) och skew (index 2 och 4).

Standardvärdet är [0,0, 1,0, 0,0, 0,0, 0,0, -1,0].

Python console

>>> rst = GDALRaster({"width": 10, "height": 20, "srid": 4326})
>>> rst.geotransform
[0.0, 1.0, 0.0, 0.0, 0.0, -1.0]

originLink to this definition#

Koordinater för rastrets övre vänstra ursprung i källans rumsliga referenssystem, som ett punktobjekt med x och y medlemmar.

Python console

>>> rst = GDALRaster({"width": 10, "height": 20, "srid": 4326})
>>> rst.origin
[0.0, 0.0]
>>> rst.origin.x = 1
>>> rst.origin
[1.0, 0.0]

scaleLink to this definition#

Pixelbredd och -höjd som används för georeferering av rastret, som ett punktobjekt med x och y medlemmar. Se geotransform för mer information.

Python console

>>> rst = GDALRaster({"width": 10, "height": 20, "srid": 4326})
>>> rst.scale
[1.0, -1.0]
>>> rst.scale.x = 2
>>> rst.scale
[2.0, -1.0]

skewLink to this definition#

Skevhetskoefficienter som används för att georeferera rastret, som ett punktobjekt med x och y medlemmar. När det gäller bilder från norr är båda dessa koefficienter 0.

Python console

>>> rst = GDALRaster({"width": 10, "height": 20, "srid": 4326})
>>> rst.skew
[0.0, 0.0]
>>> rst.skew.x = 3
>>> rst.skew
[3.0, 0.0]

extentLink to this definition#

Rasterkällans utbredning (gränsvärden), som en 4-tupel (xmin, ymin, xmax, ymax) i källans rumsliga referenssystem.

Python console

>>> rst = GDALRaster({"width": 10, "height": 20, "srid": 4326})
>>> rst.extent
(0.0, -20.0, 10.0, 0.0)
>>> rst.origin.x = 100
>>> rst.extent
(100.0, -20.0, 110.0, 0.0)

bandsLink to this definition#

Lista över alla band i källan, som GDALBand-instanser.

Python console

>>> rst = GDALRaster(
...     {
...         "width": 1,
...         "height": 2,
...         "srid": 4326,
...         "bands": [{"data": [0, 1]}, {"data": [2, 3]}],
...     }
... )
>>> len(rst.bands)
2
>>> rst.bands[1].data()
array([[ 2.,  3.]], dtype=float32)

warp(ds_input, resampling='NearestNeighbour', max_error=0.0)Link to this definition#

Returnerar en förvrängd version av detta raster.

Warpingparametrarna kan specificeras genom argumentet ds_input. Användningen av ds_input är analog med motsvarande argument i klasskonstruktorn. Det är en ordbok med egenskaperna hos målrastret. Tillåtna nyckelvärden i ordboken är bredd, höjd, SRID, ursprung, skala, skevhet, datatyp, drivrutin och namn (filnamn).

Som standard behåller warp-funktionen de flesta parametrar till samma värden som i det ursprungliga källrastret, så endast parametrar som ska ändras behöver anges. Observera att detta inkluderar drivrutinen, så för filbaserade raster skapar warp-funktionen ett nytt raster på disken.

Den enda parameter som skiljer sig från källrastret är namnet. Standardvärdet för rasternamnet är namnet på källrastret med tillägget '_copy' + source_driver_name. För filbaserade raster rekommenderas att man anger filvägen till målrastret.

Den resamplingsalgoritm som används för warping kan specificeras med argumentet resampling. Standardvärdet är NearestNeighbor och de andra tillåtna värdena är Bilinear, Cubic, CubicSpline, Lanczos, Average och Mode.

Argumentet max_error kan användas för att ange det maximala felet mätt i ingående pixlar som tillåts vid approximering av transformationen. Standardvärdet är 0,0 för exakta beräkningar.

För användare som är bekanta med GDAL har denna funktion en liknande funktionalitet som kommandoradsverktyget gdalwarp.

Warp-funktionen kan t.ex. användas för att aggregera ett raster till det dubbla av den ursprungliga pixelskalan:

Python console

>>> rst = GDALRaster(
...     {
...         "width": 6,
...         "height": 6,
...         "srid": 3086,
...         "origin": [500000, 400000],
...         "scale": [100, -100],
...         "bands": [{"data": range(36), "nodata_value": 99}],
...     }
... )
>>> target = rst.warp({"scale": [200, -200], "width": 3, "height": 3})
>>> target.bands[0].data()
array([[  7.,   9.,  11.],
       [ 19.,  21.,  23.],
       [ 31.,  33.,  35.]], dtype=float32)

transform(srs, driver=None, name=None, resampling='NearestNeighbour', max_error=0.0)Link to this definition#

Transformerar detta raster till ett annat rumsligt referenssystem (srs), som kan vara ett SpatialReference-objekt, eller någon annan indata som accepteras av SpatialReference (inklusive rumsliga referens WKT- och PROJ-strängar, eller ett heltal SRID).

Den beräknar gränserna och skalan för det aktuella rastret i det nya rumsliga referenssystemet och förvränger rastret med hjälp av funktionen warp.

Som standard används drivrutinen för källrastret och namnet på rastret är det ursprungliga namnet med tillägget '_copy' + source_driver_name. En annan drivrutin eller ett annat namn kan anges med argumenten driver och name.

Standardalgoritmen för omsampling är NearestNeighbour men kan ändras med argumentet resampling. Det högsta tillåtna standardfelet för omsampling är 0,0 och kan ändras med argumentet max_error. Se dokumentationen warp för mer information om dessa argument.

Python console

>>> rst = GDALRaster(
...     {
...         "width": 6,
...         "height": 6,
...         "srid": 3086,
...         "origin": [500000, 400000],
...         "scale": [100, -100],
...         "bands": [{"data": range(36), "nodata_value": 99}],
...     }
... )
>>> target_srs = SpatialReference(4326)
>>> target = rst.transform(target_srs)
>>> target.origin
[-82.98492744885776, 27.601924753080144]

infoLink to this definition#

Returnerar en sträng med en sammanfattning av rastret. Detta motsvarar kommandoradsverktyget `gdalinfo`__.

metadataLink to this definition#

Metadata för detta raster, representerat som en nästlad ordbok. Nyckeln på första nivån är metadatadomänen. Den andra nivån innehåller metadataelementens namn och värden från varje domän.

För att ställa in eller uppdatera ett metadataobjekt skickar du motsvarande metadataobjekt till metoden med hjälp av den nästlade struktur som beskrivs ovan. Endast nycklar som finns i den angivna ordlistan uppdateras; resten av metadata förblir oförändrade.

Om du vill ta bort ett metadataobjekt använder du None som metadatavärde.

Python console

>>> rst = GDALRaster({"width": 10, "height": 20, "srid": 4326})
>>> rst.metadata
{}
>>> rst.metadata = {"DEFAULT": {"OWNER": "Django", "VERSION": "1.0"}}
>>> rst.metadata
{'DEFAULT': {'OWNER': 'Django', 'VERSION': '1.0'}}
>>> rst.metadata = {"DEFAULT": {"OWNER": None, "VERSION": "2.0"}}
>>> rst.metadata
{'DEFAULT': {'VERSION': '2.0'}}

vsi_bufferLink to this definition#

En bytes``representation av detta raster. Returnerar ``None för raster som inte lagras i GDAL:s virtuella filsystem.

is_vsi_basedLink to this definition#

En boolean som anger om detta raster lagras i GDAL:s virtuella filsystem.

GDALBandLink to this heading#

class GDALBandLink to this definition#

instanser av GDALBand skapas inte explicit utan erhålls från ett GDALRaster-objekt genom dess bands-attribut. GDALBands innehåller de faktiska pixelvärdena i rastret.

descriptionLink to this definition#

Namn eller beskrivning av bandet, om sådant finns.

widthLink to this definition#

Bandets bredd i pixlar (X-axeln).

heightLink to this definition#

Bandets höjd i pixlar (Y-axeln).

pixel_countLink to this definition#

Det totala antalet pixlar i detta band. Är lika med bredd * höjd.

statistics(refresh=False, approximate=False)Link to this definition#

Beräknar statistik över pixelvärdena i detta band. Returvärdet är en tupel med följande struktur: (minimum, maximum, mean, standard deviation).

Om argumentet approximate är inställt på True kan statistiken beräknas baserat på översikter eller en delmängd av bildplattor.

Om argumentet refresh är satt till True, kommer statistiken att beräknas direkt från data och cachen kommer att uppdateras med resultatet.

Om det finns ett värde i den permanenta cachen returneras det värdet. För rasterformat som använder PAM-tjänster (Persistent Auxiliary Metadata) kan statistiken cachelagras i en hjälpfil. I vissa fall kan dessa metadata vara osynkroniserade med pixelvärdena eller leda till att värden från ett tidigare anrop returneras som inte återspeglar värdet på argumentet approximate. I sådana fall ska du använda argumentet refresh för att hämta uppdaterade värden och lagra dem i cacheminnet.

För tomma band (där alla pixelvärden är ”inga data”) returneras all statistik som None.

Statistiken kan också hämtas direkt genom att använda egenskaperna min, max, mean och std.

minLink to this definition#

Bandets lägsta pixelvärde (exklusive värdet för ”inga data”).

maxLink to this definition#

Det maximala pixelvärdet för bandet (exklusive värdet ”inga data”).

meanLink to this definition#

Medelvärdet av alla pixelvärden i bandet (exklusive värdet ”inga data”).

stdLink to this definition#

Standardavvikelsen för alla pixelvärden i bandet (exklusive värdet ”inga data”).

nodata_valueLink to this definition#

Värdet ”no data” för ett band är i allmänhet ett speciellt markörvärde som används för att markera pixlar som inte är giltiga data. Sådana pixlar ska i allmänhet inte visas och inte heller bidra till analysoperationer.

För att ta bort ett befintligt ”no data”-värde, sätt denna egenskap till None.

datatype(as_string=False)Link to this definition#

Den datatyp som bandet innehåller, som en heltalskonstant mellan 0 (okänd) och 14. Om as_string är True returneras datatypen som en sträng. Kolla in kolumnen ”GDAL Pixel Type” i datatype value table för möjliga värden.

color_interp(as_string=False)Link to this definition#

Färgtolkningen för bandet, som ett heltal mellan 0 och 16. Om as_string är True returneras datatypen som en sträng med följande möjliga värden: GCI_Undefined, GCI_GrayIndex, GCI_PaletteIndex, GCI_RedBand, GCI_GreenBand, GCI_BlueBand, GCI_AlphaBand, GCI_HueBand, GCI_SaturationBand, GCI_LightnessBand, GCI_CyanBand, GCI_MagentaBand, GCI_YellowBand, GCI_BlackBand, GCI_YCbCr_YBand, GCI_YCbCr_CbBand och GCI_YCbCr_CrBand. GCI_YCbCr_CrBand representerar också GCI_Max eftersom båda motsvarar heltalet 16, men endast GCI_YCbCr_CrBand returneras som en sträng.

data(data=None, offset=None, size=None, shape=None)Link to this definition#

Åtkomst till pixelvärdena i GDALBand. Returnerar den fullständiga dataarrayen om inga parametrar anges. En delmängd av pixelmatrisen kan begäras genom att ange en offset och blockstorlek som tupler.

Om NumPy är tillgängligt returneras data som en NumPy-array. Av prestandaskäl rekommenderas det starkt att använda NumPy.

Data skrivs till GDALBand om parametern data anges. Indata kan vara av någon av följande typer - packad sträng, buffert, lista, array och NumPy-array. Antalet objekt i indata bör normalt motsvara det totala antalet pixlar i bandet, eller antalet pixlar för ett specifikt block av pixelvärden om parametrarna offset och size anges.

Om antalet objekt i indata skiljer sig från målpixelblocket måste parametern shape anges. Formen är en tupel som anger bredden och höjden på indata i pixlar. Data replikeras sedan för att uppdatera pixelvärdena i det valda blocket. Detta är användbart om man t.ex. vill fylla ett helt band med ett enda värde.

Till exempel:

Python console

>>> rst = GDALRaster(
...     {"width": 4, "height": 4, "srid": 4326, "datatype": 1, "nr_of_bands": 1}
... )
>>> bnd = rst.bands[0]
>>> bnd.data(range(16))
>>> bnd.data()
array([[ 0,  1,  2,  3],
       [ 4,  5,  6,  7],
       [ 8,  9, 10, 11],
       [12, 13, 14, 15]], dtype=int8)
>>> bnd.data(offset=(1, 1), size=(2, 2))
array([[ 5,  6],
       [ 9, 10]], dtype=int8)
>>> bnd.data(data=[-1, -2, -3, -4], offset=(1, 1), size=(2, 2))
>>> bnd.data()
array([[ 0,  1,  2,  3],
       [ 4, -1, -2,  7],
       [ 8, -3, -4, 11],
       [12, 13, 14, 15]], dtype=int8)
>>> bnd.data(data="\x9d\xa8\xb3\xbe", offset=(1, 1), size=(2, 2))
>>> bnd.data()
array([[  0,   1,   2,   3],
       [  4, -99, -88,   7],
       [  8, -77, -66,  11],
       [ 12,  13,  14,  15]], dtype=int8)
>>> bnd.data([1], shape=(1, 1))
>>> bnd.data()
array([[1, 1, 1, 1],
       [1, 1, 1, 1],
       [1, 1, 1, 1],
       [1, 1, 1, 1]], dtype=uint8)
>>> bnd.data(range(4), shape=(1, 4))
array([[0, 0, 0, 0],
       [1, 1, 1, 1],
       [2, 2, 2, 2],
       [3, 3, 3, 3]], dtype=uint8)

metadataLink to this definition#

Metadata för detta band. Funktionaliteten är identisk med GDALRaster.metadata.

Skapa raster från dataLink to this heading#

I detta avsnitt beskrivs hur du skapar raster från grunden med hjälp av parametern ds_input.

Ett nytt raster skapas när en dict skickas till GDALRaster-konstruktören. Dictionariet innehåller definierande parametrar för det nya rastret, t.ex. ursprung, storlek eller rumsligt referenssystem. Dictionariet kan också innehålla pixeldata och information om formatet för det nya rastret. Det resulterande rastret kan därför vara filbaserat eller minnesbaserat, beroende på vilken drivrutin som anges.

Det finns ingen standard för att beskriva rasterdata i en ordbok eller JSON-smak. Definitionen av ordboksinmatningen till GDALRaster-klassen är därför specifik för Django. Den är inspirerad av `geojson`__-formatet, men geojson-standarden är för närvarande begränsad till vektorformat.

Exempel på hur man använder de olika nycklarna när man skapar raster finns i dokumentationen för motsvarande attribut och metoder i klasserna GDALRaster och GDALBand.

Använda GDAL:s virtuella filsystemLink to this heading#

GDAL kan komma åt filer som lagras i filsystemet, men stöder också virtuella filsystem för att abstrahera åtkomst till andra typer av filer, t.ex. komprimerade, krypterade eller fjärrfiler.

# Read a raster as a file object from a remote source.
>>> from urllib.request import urlopen
>>> dat = urlopen("https://example.com/raster.tif").read()
# Instantiate a raster from the bytes object.
>>> rst = GDALRaster(dat)
# The name starts with /vsimem/, indicating that the raster lives in the
# virtual filesystem.
>>> rst.name
'/vsimem/da300bdb-129d-49a8-b336-e410a9428dad'

>>> from django.http import HttpResponse
>>> rst = GDALRaster(
...     {
...         "name": "/vsimem/temporarymemfile",
...         "driver": "tif",
...         "width": 6,
...         "height": 6,
...         "srid": 3086,
...         "origin": [500000, 400000],
...         "scale": [100, -100],
...         "bands": [{"data": range(36), "nodata_value": 99}],
...     }
... )
>>> HttpResponse(rast.vsi_buffer, "image/tiff")

>>> from django.contrib.gis.gdal import GDALRaster
>>> rst = GDALRaster("/vsizip/path/to/your/file.zip/path/to/raster.tif")
>>> rst = GDALRaster("/vsigzip/path/to/your/file.gz")
>>> rst = GDALRaster("/vsitar/path/to/your/file.tar/path/to/raster.tif")

>>> from django.contrib.gis.gdal import GDALRaster
>>> rst = GDALRaster("/vsicurl/https://example.com/raster.tif")
>>> rst.name
'/vsicurl/https://example.com/raster.tif'

GDAL_BIBLIOTEK_VÄGLink to this heading#

En sträng som anger platsen för GDAL-biblioteket. Vanligtvis används denna inställning endast om GDAL-biblioteket finns på en plats som inte är standard (t.ex. /home/john/lib/libgdal.so).