GDAL API — Django 文档

来自菜鸟教程
Django/docs/2.2.x/ref/contrib/gis/gdal
跳转至:导航、​搜索

GDAL API

GDAL代表地理空间数据抽象库,是名副其实的GIS数据功能“瑞士军刀”。 GDAL 的一个子集是 OGR 简单特征库,它专门用于读取和写入各种标准格式的矢量地理数据。

GeoDjango 为 OGR 的一些功能提供了一个高级 Python 接口,包括矢量空间数据的读取和坐标转换以及对 GDAL 与栅格(图像)数据相关的功能的最小支持。

笔记

虽然该模块被命名为gdal,但GeoDjango目前仅支持OGR和GDAL的光栅特性的部分功能。


概况

样本数据

此处描述的 GDAL/OGR 工具旨在帮助您读入您的地理空间数据,为了使它们中的大多数有用,您必须有一些数据可以使用。 如果您刚开始并且还没有任何自己的数据可供使用,GeoDjango 测试包含许多可用于测试的简单数据集。 你可以在这里下载它们:

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

矢量数据源对象

DataSource

DataSource 是 OGR 数据源对象的包装器,支持使用简单、一致的接口从各种 OGR 支持的地理空间文件格式和数据源读取数据。 每个数据源由一个 DataSource 对象表示,该对象包含一层或多层数据。 每个图层,由 Layer 对象表示,包含一定数量的地理特征 (Feature),关于该图层中包含的特征类型的信息(例如 点、多边形等),以及可能与该图层中的每个要素相关联的任何其他数据字段 (Field) 的名称和类型。

class DataSource(ds_input, encoding='utf-8')

DataSource 的构造函数只需要一个参数:要读取的文件的路径。 但是,OGR 还支持各种更复杂的数据源,包括数据库,可以通过传递特殊名称字符串而不是路径来访问这些数据源。 有关详细信息,请参阅 OGR 矢量格式 文档。 DataSource 实例的 name 属性给出了它正在使用的底层数据源的 OGR 名称。

可选的 encoding 参数允许您指定源中字符串的非标准编码。 当您在读取字段值时获得 DjangoUnicodeDecodeError 异常时,这通常很有用。

创建 DataSource 后,您可以通过访问 layer_count 属性或(等效地)使用 len() 函数找出它包含多少层数据. 有关访问数据层本身的信息,请参阅下一节:

>>> 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_count

返回数据源中的层数。

name

返回数据源的名称。


Layer

class Layer

LayerDataSource 对象中一层数据的包装器。 您永远不会直接创建 Layer 对象。 相反,您可以从 DataSource 对象中检索它们,该对象本质上是 Layer 对象的标准 Python 容器。 例如,您可以通过索引访问特定图层(例如 ds[0] 访问第一层),或者您可以在 for 循环中遍历容器中的所有层。 Layer 本身充当几何特征的容器。

通常,给定图层中的所有要素都具有相同的几何类型。 图层的 geom_type 属性是标识要素类型的 OGRGeomType。 我们可以用它打印出DataSource中每一层的一些基本信息:

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

示例输出来自上面加载的城市数据源,其中显然包含一层,称为 "cities",其中包含三个点要素。 为简单起见,以下示例假设您已将该图层存储在变量 layer 中:

>>> layer = ds[0]
name

返回数据源中该图层的名称。

>>> layer.name
'cities'
num_feat

返回图层中的要素数量。 与 len(layer) 相同:

>>> layer.num_feat
3
geom_type

返回层的几何类型,作为 OGRGeomType 对象:

>>> layer.geom_type.name
'Point'
num_fields

返回图层中的字段数,即与图层中每个要素关联的数据的字段数:

>>> layer.num_fields
4
fields

返回该层中每个字段的名称列表:

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

返回该层中每个字段的数据类型列表。 这些是 Field 的子类,讨论如下:

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

返回该层中每个字段的最大字段宽度列表:

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

返回该层中每个字段的数值精度列表。 这对于非数字字段毫无意义(并设置为零):

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

返回该层的空间范围,作为 Envelope 对象:

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

返回与此图层关联的 SpatialReference 的属性:

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

如果 Layer 没有与之关联的空间参考信息,则返回 None

spatial_filter

可用于检索或设置此图层的空间过滤器的属性。 空间过滤器只能使用 OGRGeometry 实例、4 元组范围或 None 设置。 当设置为 None 以外的其他内容时,在迭代图层时只会返回与过滤器相交的特征:

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

为图层中的每个要素返回给定字段的值列表的方法:

>>> layer.get_fields('Name')
['Pueblo', 'Lawrence', 'Houston']
get_geoms(geos=False)

一种返回包含图层中每个要素几何的列表的方法。 如果可选参数 geos 设置为 True,则几何图形将转换为 GEOSGeometry 对象。 否则,它们将作为 OGRGeometry 对象返回:

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

返回一个布尔值,指示该层是否支持给定的功能(字符串)。 有效功能字符串的示例包括:'RandomRead''SequentialWrite''RandomWrite''FastSpatialFilter''FastFeatureCount''FastGetExtent'、[ X112X]、'Transactions''DeleteFeature''FastSetNextByIndex'


Feature

class Feature

Feature 包装了一个 OGR 功能。 您永远不会直接创建 Feature 对象。 相反,您可以从 Layer 对象中检索它们。 每个要素都由一个几何图形和一组包含附加属性的字段组成。 字段的几何形状可通过其 geom 属性访问,该属性返回一个 OGRGeometry 对象。 Feature 的行为类似于其字段的标准 Python 容器,它作为 Field 对象返回:您可以通过其索引或名称直接访问字段,或者您可以迭代功能的字段,例如 在 for 循环中。

geom

返回此特征的几何图形,作为 OGRGeometry 对象:

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

返回此功能的给定字段(由名称指定)的值的方法,not 一个 Field 包装器对象:

>>> city.get('Population')
102121
geom_type

返回此特征的几何类型,作为 OGRGeomType 对象。 对于给定图层中的所有要素,这将相同,并且等效于该要素来自的 Layer 对象的 Layer.geom_type 属性。

num_fields

返回与特征关联的数据字段数。 对于给定图层中的所有要素,这将相同,并且等效于该要素来自的 Layer 对象的 Layer.num_fields 属性。

fields

返回与要素关联的数据字段的名称列表。 对于给定图层中的所有要素,这将相同,并且等效于该要素来自的 Layer 对象的 Layer.fields 属性。

fid

返回图层内的特征标识符:

>>> city.fid
0
layer_name

返回特征来自的 的名称。 对于给定层中的所有特征,这将是相同的:

>>> city.layer_name
'cities'
index

返回给定字段名称索引的方法。 对于给定层中的所有特征,这将是相同的:

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


Field

class Field
name

返回此字段的名称:

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

以整数形式返回此字段的 OGR 类型。 FIELD_CLASSES 字典将这些值映射到 Field 的子类上:

>>> city['Density'].type
2
type_name

返回具有此字段数据类型名称的字符串:

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

返回此字段的值。 Field 类本身以字符串形式返回值,但每个子类都以最合适的形式返回值:

>>> city['Population'].value
102121
width

返回此字段的宽度:

>>> city['Name'].width
80
precision

返回此字段的数字精度。 这对于非数字字段毫无意义(并设置为零):

>>> city['Density'].precision
15
as_double()

以双精度(浮点数)形式返回该字段的值:

>>> city['Density'].as_double()
874.7
as_int()

以整数形式返回字段的值:

>>> city['Population'].as_int()
102121
as_string()

以字符串形式返回字段的值:

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

以日期和时间组件的元组形式返回字段的值:

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


Driver

class Driver(dr_input)

Driver 类在内部用于包装 OGR DataSource 驱动程序。

driver_count

返回当前注册的 OGR 矢量驱动程序的数量。


OGR 几何图形

OGRGeometry

OGRGeometry 对象与 GEEOSGeometry 对象共享相似的功能,并且是 OGR 内部几何表示的薄包装。 因此,当使用 DataSource 时,它们允许更有效地访问数据。 与其对应的 GEOS 不同,OGRGeometry 支持空间参考系统和坐标变换:

>>> from django.contrib.gis.gdal import OGRGeometry
>>> polygon = OGRGeometry('POLYGON((0 0, 5 0, 5 5, 0 5))')
class OGRGeometry(geom_input, srs=None)

该对象是 OGR Geometry 类的包装器。 这些对象直接从给定的 geom_input 参数实例化,该参数可能是包含 WKT、HEX、GeoJSON 的字符串、包含 WKB 数据的 bufferOGRGeomType 对象。 当从 Layer(又是 DataSource 的一部分)读取矢量数据时,这些对象也从 Feature.geom 属性返回。

classmethod from_gml(gml_string)

从给定的 GML 字符串构造一个 OGRGeometry

classmethod from_bbox(bbox)

从给定的边界框(一个 4 元组)构造一个 Polygon

__len__()

返回 LineString 中的点数、Polygon 中的环数或 GeometryCollection 中的几何数。 不适用于其他几何类型。

__iter__()

迭代 LineString 中的点、Polygon 中的环或 GeometryCollection 中的几何。 不适用于其他几何类型。

__getitem__()

返回 LineString 指定索引处的点、Polygon 指定索引处的内环或 GeometryCollection 中指定索引处的几何. 不适用于其他几何类型。

dimension

返回几何的协调维数,即 0 代表点,1 代表线,依此类推:

>> polygon.dimension
2
coord_dim

返回或设置此几何图形的坐标尺寸。 例如,二维几何图形的值为 2。

geom_count

返回此几何中的元素数:

>>> polygon.geom_count
1
point_count

返回用于描述此几何图形的点数:

>>> polygon.point_count
4
num_points

point_count 的别名。

num_coords

point_count 的别名。

geom_type

返回此几何的类型,作为 OGRGeomType 对象。

geom_name

返回此几何类型的名称:

>>> polygon.geom_name
'POLYGON'
area

返回此几何图形的面积,对于不包含面积的几何图形,返回 0:

>>> polygon.area
25.0
envelope

返回此几何体的包络,作为 Envelope 对象。

extent

将此几何的包络作为 4 元组返回,而不是作为 Envelope 对象:

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

此属性控制此几何图形的空间参考,如果尚未为其指定空间参考系统,则控制 None。 如果已分配,访问此属性将返回一个 SpatialReference 对象。 可以使用另一个 SpatialReference 对象或 SpatialReference 接受的任何输入来设置它。 例子:

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

返回或设置与此几何的 SpatialReference 对应的空间参考标识符。 如果没有与此几何相关联的空间参考信息,或者无法确定 SRID,则返回 None

geos

返回与此几何体对应的 GEOSGeometry 对象。

gml

以 GML 格式返回此几何图形的字符串表示形式:

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

以 HEX WKB 格式返回此几何图形的字符串表示形式:

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

以 JSON 格式返回此几何图形的字符串表示形式:

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

以 KML 格式返回此几何图形的字符串表示形式。

wkb_size

返回保存此几何图形的 WKB 表示所需的 WKB 缓冲区的大小:

>>> OGRGeometry('POINT(1 2)').wkb_size
21
wkb

返回包含此几何图形的 WKB 表示的 buffer

wkt

以 WKT 格式返回此几何图形的字符串表示形式。

ewkt

返回此几何的 EWKT 表示。

clone()

返回此几何对象的新 OGRGeometry 克隆。

close_rings()

如果此几何图形中有任何尚未闭合的环,则此例程将通过将起点添加到终点来完成:

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

将此几何图形转换为不同的空间参考系统。 可以采用 CoordTransform 对象、SpatialReference 对象或 SpatialReference 接受的任何其他输入(包括空间参考 WKT 和 PROJ.4 字符串,或整数 SRID )。

默认情况下,不返回任何内容,并且几何图形就地转换。 但是,如果 clone 关键字设置为 True,则返回此几何体的转换克隆。

intersects(other)

如果此几何图形与另一个几何图形相交,则返回 True,否则返回 False

equals(other)

如果此几何图形与另一个几何图形等效,则返回 True,否则返回 False

disjoint(other)

如果此几何体在空间上与(即 不相交)另一个,否则返回 False

touches(other)

如果此几何体与另一个几何体接触,则返回 True,否则返回 False

crosses(other)

如果此几何图形与另一个几何图形交叉,则返回 True,否则返回 False

within(other)

如果此几何图形包含在另一个几何图形中,则返回 True,否则返回 False

contains(other)

如果此几何包含另一个,则返回 True,否则返回 False

overlaps(other)

如果此几何图形与另一个几何图形重叠,则返回 True,否则返回 False

boundary()

此几何体的边界,作为一个新的 OGRGeometry 对象。

convex_hull

包含此几何体的最小凸多边形,作为新的 OGRGeometry 对象。

difference()

以新的 OGRGeometry 对象的形式返回由该几何体与另一个几何体的差异组成的区域。

intersection()

返回由这个几何体和另一个几何体的交集组成的区域,作为一个新的 OGRGeometry 对象。

sym_difference()

返回由该几何图形和另一个几何图形的对称差异组成的区域,作为新的 OGRGeometry 对象。

union()

返回由这个几何体和另一个几何体的联合组成的区域,作为一个新的 OGRGGeometry 对象。

tuple

将点几何的坐标作为元组返回,线几何的坐标作为元组的元组返回,依此类推:

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

tuple 的别名。

class Point
x

返回该点的 X 坐标:

>>> OGRGeometry('POINT (1 2)').x
1.0
y

返回该点的 Y 坐标:

>>> OGRGeometry('POINT (1 2)').y
2.0
z

返回该点的 Z 坐标,如果该点没有 Z 坐标,则返回 None

>>> OGRGeometry('POINT (1 2 3)').z
3.0
class LineString
x

返回此行中的 X 坐标列表:

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

返回此行中的 Y 坐标列表:

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

返回此行中的 Z 坐标列表,如果该行没有 Z 坐标,则返回 None

>>> OGRGeometry('LINESTRING (1 2 3,4 5 6)').z
[3.0, 6.0]
class Polygon
shell

返回此多边形的壳或外环,作为 LinearRing 几何图形。

exterior_ring

shell 的别名。

centroid

返回表示此多边形质心的 Point

class GeometryCollection
add(geom)

向此几何集合添加几何。 不适用于其他几何类型。


OGRGeomType

class OGRGeomType(type_input)

此类允许以多种方式中的任何一种表示 OGR 几何类型:

>>> 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
name

返回 OGR Geometry 类型的简写字符串形式:

>>> gt1.name
'Polygon'
num

返回对应于 OGR 几何类型的数字:

>>> gt1.num
3
django

返回用于存储此 OGR 类型的 Django 字段类型(GeometryField 的子类),如果没有合适的 Django 类型,则返回 None

>>> gt1.django
'PolygonField'


Envelope

class Envelope(*args)

表示包含矩形边界框的最小和最大 X、Y 坐标的 OGR Envelope 结构。 变量的命名与 OGR Envelope C 结构兼容。

min_x

最小 X 坐标值。

min_y

最大 X 坐标值。

max_x

最小 Y 坐标值。

max_y

最大 Y 坐标值。

ur

右上角的坐标,作为一个元组。

ll

左下角坐标,作为元组。

tuple

代表信封的元组。

wkt

将此信封表示为 WKT 格式的多边形的字符串。

expand_to_include(*args)


坐标系对象

SpatialReference

class SpatialReference(srs_input)

空间参考对象在给定的 srs_input 上初始化,它可能是以下之一:

  • OGC 众所周知的文本 (WKT)(字符串)

  • EPSG 代码(整数或字符串)

  • PROJ.4 字符串

  • 知名标准的简写字符串('WGS84''WGS72''NAD27''NAD83'

例子:

>>> wgs84 = SpatialReference('WGS84') # shorthand string
>>> wgs84 = SpatialReference(4326) # EPSG code
>>> wgs84 = SpatialReference('EPSG:4326') # EPSG string
>>> proj4 = '+proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs '
>>> wgs84 = SpatialReference(proj4) # PROJ.4 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)

返回给定字符串属性节点的值,如果节点不存在,则返回 None。 也可以带一个元组作为参数,(target, child),其中child是属性在WKT中的索引。 例如:

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

给定目标节点的属性值(例如 'PROJCS')。 index 关键字指定要返回的子节点的索引。

auth_name(target)

返回给定字符串目标节点的权限名称。

auth_code(target)

返回给定字符串目标节点的权限代码。

clone()

返回此空间参考对象的副本。

identify_epsg()

此方法检查此 SpatialReference 的 WKT,并将在 EPSG 标识符适用的地方添加 EPSG 授权节点。

from_esri()

将此 SpatialReference 从 ESRI 的格式变形为 EPSG

to_esri()

将此 SpatialReference 变形为 ESRI 的格式。

validate()

检查给定的空间参考是否有效,否则将引发异常。

import_epsg(epsg)

从 EPSG 代码导入空间参考。

import_proj(proj)

从 PROJ.4 字符串导入空间参考。

import_user_input(user_input)
import_wkt(wkt)

从 WKT 导入空间参考。

import_xml(xml)

从 XML 导入空间参考。

name

返回此空间参考的名称。

srid

返回顶级权限的 SRID,如果未定义,则返回 None

linear_name

返回线性单位的名称。

linear_units

返回线性单位的值。

angular_name

返回角度单位的名称。”

angular_units

返回角度单位的值。

units

返回单位值和单位名称的 2 元组,并将自动确定是返回线性单位还是角度单位。

ellipsoid

返回此空间参考的椭球参数元组:(长半轴、短半轴和反向展平)。

semi_major

返回此空间参考的椭球的半长轴。

semi_minor

返回此空间参考的椭球的短半轴。

inverse_flattening

返回此空间参考的椭圆体的逆展平。

geographic

如果此空间参考是地理参考(根节点为 GEOGCS),则返回 True

local

如果此空间参考是本地的(根节点是 LOCAL_CS),则返回 True

projected

如果此空间参考是投影坐标系(根节点为 PROJCS),则返回 True

wkt

返回此空间参考的 WKT 表示。

pretty_wkt

返回 WKT 的“漂亮”表示。

proj

返回此空间参考的 PROJ.4 表示。

proj4

SpatialReference.proj 的别名。

xml

返回此空间参考的 XML 表示。


CoordTransform

class CoordTransform(source, target)

表示坐标系变换。 它用两个 SpatialReference 初始化,分别代表源坐标系和目标坐标系。 在不同几何图形上重复执行相同的坐标变换时,应使用这些对象:

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

栅格数据对象

GDALRaster

GDALRaster 是 GDAL 栅格源对象的包装器,支持使用简单、一致的接口从各种 GDAL 支持的地理空间文件格式和数据源读取数据。 每个数据源都由一个 GDALRaster 对象表示,该对象包含一个或多个名为 band 的数据层。 每个波段由 GDALBand 对象表示,包含地理参考图像数据。 例如,RGB 图像表示为三个波段:一个用于红色,一个用于绿色,一个用于蓝色。

笔记

对于栅格数据,栅格实例与其数据源之间没有区别。 与 Geometry 对象不同,GDALRaster 对象始终是数据源。 临时栅格可以使用相应的驱动程序在内存中实例化,但它们将与基于文件的栅格源属于同一类。


class GDALRaster(ds_input, write=False)

GDALRaster 的构造函数接受两个参数。 第一个参数定义栅格源,第二个参数定义是否应以写入模式打开栅格。 对于新创建的栅格,第二个参数将被忽略,并且新栅格始终以写入模式创建。

第一个参数可以采用三种形式:表示文件路径的字符串、具有定义新光栅的值的字典或表示光栅文件的字节对象。

如果输入是文件路径,则从那里打开栅格。 如果输入是字典中的原始数据,则需要参数widthheightsrid。 如果输入是字节对象,它将使用 GDAL 虚拟文件系统打开。

有关如何使用字典输入创建栅格的详细说明,请参阅 从数据创建栅格 。 有关如何在虚拟文件系统中创建栅格的详细说明,请参阅 使用 GDAL 的虚拟文件系统

以下示例显示了如何从不同的输入源创建栅格(使用来自 GeoDjango 测试的示例数据;另见 示例数据 部分)。

>>> 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'
name

源的名称,相当于输入文件路径或实例化时提供的名称。

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

用于处理输入文件的 GDAL 驱动程序的名称。 对于从文件创建的 GDALRaster,会自动检测驱动程序类型。 默认情况下,从头开始创建的栅格是内存中的栅格 ('MEM'),但可以根据需要进行更改。 例如,对 GeoTiff 文件使用 GTiff。 有关文件类型列表,另请参阅 GDAL 光栅格式 列表。

通过以下示例创建内存中栅格:

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

通过以下示例创建基于文件的 GeoTiff 栅格:

>>> 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'
width

源的宽度(以像素为单位)(X 轴)。

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

源的高度(以像素为单位)(Y 轴)。

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

栅格的空间参考系统,作为 SpatialReference 实例。 可以通过将 SRS 设置为其他 SpatialReference 或提供 SpatialReference 构造函数接受的任何输入来更改 SRS。

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

栅格的空间参考系统标识符 (SRID)。 此属性是通过 srs 属性获取或设置 SRID 的快捷方式。

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

用于对源进行地理配准的仿射变换矩阵,作为六个系数的元组,使用以下关系将像素/线坐标映射到地理配准空间:

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

可以通过访问 origin(索引 0 和 3)、scale(索引 1 和 5)和 skew(索引 2 和 4)来检索相同的值特性。

默认值为 [0.0, 1.0, 0.0, 0.0, 0.0, -1.0]

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

源空间参考系统中栅格左上角原点的坐标,作为具有 xy 成员的点对象。

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

用于对栅格进行地理配准的像素宽度和高度,作为具有 xy 成员的点对象。 有关详细信息,请参阅 geotransform

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

用于对栅格进行地理配准的倾斜系数,作为具有 xy 成员的点对象。 在北上图像的情况下,这些系数都是 0

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

栅格源的范围(边界值),作为源空间参考系统中的 4 元组 (xmin, ymin, xmax, ymax)

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

源的所有波段列表,作为 GDALBand 实例。

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

返回此栅格的扭曲版本。

可以通过 ds_input 参数指定变形参数。 ds_input 的使用类似于类构造函数的相应参数。 它是具有目标栅格特征的字典。 允许的字典键值为宽度、高度、SRID、原点、比例、倾斜、数据类型、驱动程序和名称(文件名)。

默认情况下,warp 函数使大多数参数与原始源栅格的值相同,因此只需要指定应该更改的参数。 请注意,这包括驱动程序,因此对于基于文件的栅格,warp 函数将在磁盘上创建一个新栅格。

与源栅格设置不同的唯一参数是名称。 栅格名称的默认值是源栅格名称后附加 '_copy' + source_driver_name。 对于基于文件的栅格,建议提供目标栅格的文件路径。

可以使用 resampling 参数指定用于变形的重采样算法。 默认为 NearestNeighbor,其他允许值为 BilinearCubicCubicSplineLanczosAverage , 和 Mode

max_error 参数可用于指定在近似变换中允许的输入像素中测量的最大误差。 精确计算的默认值为 0.0。

对于熟悉 GDAL 的用户,此功能具有与 gdalwarp 命令行实用程序类似的功能。

例如,warp 函数可用于将栅格聚合为其原始像素比例的两倍:

>>> 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(srid, driver=None, name=None, resampling='NearestNeighbour', max_error=0.0)

返回具有指定 SRID 的此栅格的转换版本。

此函数将当前栅格转换为可使用 srid 指定的新空间参考系统。 它计算新空间参考系统中当前栅格的边界和比例,并使用 warp 函数扭曲栅格。

默认情况下,使用源栅格的驱动程序,栅格名称为原始名称后加'_copy' + source_driver_name。 可以使用 drivername 参数指定不同的驱动程序或名称。

默认重采样算法为 NearestNeighbour,但可以使用 resampling 参数进行更改。 重采样的默认最大允许误差为 0.0,可以使用 max_error 参数进行更改。 有关这些参数的详细信息,请参阅 warp 文档。

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

返回带有栅格摘要的字符串。 这相当于 gdalinfo 命令行实用程序。

metadata

此栅格的元数据,表示为嵌套字典。 第一级键是元数据域。 第二级包含来自每个域的元数据项名称和值。

要设置或更新元数据项,请将相应的元数据项传递给使用上述嵌套结构的方法。 只更新指定字典中的键; 其余元数据保持不变。

要删除元数据项,请使用 None 作为元数据值。

>>> 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_buffer

此栅格的 bytes 表示。 对于未存储在 GDAL 的虚拟文件系统中的栅格,返回 None

is_vsi_based

一个布尔值,指示此栅格是否存储在 GDAL 的虚拟文件系统中。


GDALBand

class GDALBand

GDALBand 实例不是显式创建的,而是通过 bands 属性从 GDALRaster 对象中获取的。 GDALBands 包含光栅的实际像素值。

description

乐队的名称或描述(如果有)。

width

带的宽度(以像素为单位)(X 轴)。

height

带的高度(以像素为单位)(Y 轴)。

pixel_count

此波段中的像素总数。 等于 width * height

statistics(refresh=False, approximate=False)

计算此波段的像素值的统计信息。 返回值是一个具有以下结构的元组:(minimum, maximum, mean, standard deviation)

如果 approximate 参数设置为 True,则可以根据概览或图像块的子集计算统计数据。

如果 refresh 参数设置为 True,则将直接从数据计算统计信息,并将结果更新缓存。

如果找到持久缓存值,则返回该值。 对于使用持久辅助元数据 (PAM) 服务的栅格格式,统计信息可能会缓存在辅助文件中。 在某些情况下,此元数据可能与像素值不同步,或者导致上次调用的值返回,而这些值不反映 approximate 参数的值。 在这种情况下,使用 refresh 参数来获取更新的值并将它们存储在缓存中。

对于空波段(其中所有像素值都是“无数据”),所有统计数据都以 None 形式返回。

还可以通过访问 minmaxmeanstd 属性直接检索统计信息。

min

波段的最小像素值(不包括“无数据”值)。

max

波段的最大像素值(不包括“无数据”值)。

mean

波段所有像素值的平均值(不包括“无数据”值)。

std

波段所有像素值的标准偏差(不包括“无数据”值)。

nodata_value

波段的“无数据”值通常是一个特殊的标记值,用于标记无效数据的像素。 通常不应显示此类像素,也不应参与分析操作。

要删除现有的“无数据”值,请将此属性设置为 None(需要 GDAL ≥ 2.1)。

datatype(as_string=False)

带中包含的数据类型,作为 0(未知)和 11 之间的整数常量。 如果 as_stringTrue,则数据类型以字符串形式返回,可能的值如下:GDT_UnknownGDT_ByteGDT_UInt16GDT_Int16GDT_UInt32GDT_Int32GDT_Float32GDT_Float64GDT_CInt16、 X215X] 和 GDT_CFloat64

color_interp(as_string=False)

带的颜色解释,为 0 到 16 之间的整数。 如果 as_stringTrue,则数据类型以字符串形式返回,可能的值如下:GCI_UndefinedGCI_GrayIndexGCI_PaletteIndexGCI_RedBandGCI_GreenBandGCI_BlueBandGCI_AlphaBandGCI_HueBandGCI_SaturationBand、 X215X]、GCI_MagentaBandGCI_YellowBandGCI_BlackBandGCI_YCbCr_YBandGCI_YCbCr_CbBandGCI_YCbCr_CrBand 也代表 GCI_Max,因为两者都对应整数 16,但只有 GCI_YCbCr_CrBand 作为字符串返回。

data(data=None, offset=None, size=None, shape=None)

GDALBand 像素值的访问器。 如果没有提供参数,则返回完整的数据数组。 可以通过将偏移量和块大小指定为元组来请求像素阵列的子集。

如果 NumPy 可用,则数据作为 NumPy 数组返回。 出于性能原因,强烈建议使用 NumPy。

如果提供 data 参数,数据将写入 GDALBand。 输入可以是以下类型之一 - 压缩字符串、缓冲区、列表、数组和 NumPy 数组。 如果 offsetsize 参数是假如。

如果输入中的项数与目标像素块不同,则必须指定 shape 参数。 形状是一个元组,以像素为单位指定输入数据的宽度和高度。 然后复制数据以更新所选块的像素值。 例如,这对于用单个值填充整个波段很有用。

例如:

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

该乐队的元数据。 该功能与 GDALRaster.metadata 相同。


从数据创建栅格

本节介绍如何使用 ds_input 参数从头开始创建栅格。

dict 传递给 GDALRaster 构造函数时,会创建一个新栅格。 字典包含新栅格的定义参数,例如原点、大小或空间参考系统。 字典还可以包含像素数据和有关新栅格格式的信息。 因此,生成的栅格可以是基于文件的或基于内存的,具体取决于指定的驱动程序。

在字典或 JSON 风格中描述栅格数据没有标准。 GDALRaster 类的字典输入的定义因此特定于 Django。 它的灵感来自 geojson 格式,但 geojson 标准目前仅限于矢量格式。

创建栅格时使用不同键的示例可以在 GDALRasterGDALBand 类的相应属性和方法的文档中找到。

ds_input 字典

ds_input 字典中只需要几个键即可创建栅格:widthheightsrid。 所有其他参数都有默认值(见下表)。 可以在 ds_input 字典中传递的键列表与 GDALRaster 属性密切相关但不完全相同。 许多参数直接映射到这些属性; 其他的描述如下。

下表描述了 ds_input 字典中可以设置的所有键。

默认 如何使用
srid 必须的 映射到 srid 属性
width 必须的 映射到 width 属性
height 必须的 映射到 height 属性
driver MEM 映射到 driver 属性
name 见下文
origin 0 映射到 origin 属性
scale 0 映射到 scale 属性
skew 0 映射到 width 属性
bands [] 见下文
nr_of_bands 0 见下文
datatype 6 见下文
papsz_options {} 见下文
name
表示栅格名称的字符串。 创建基于文件的栅格时,此参数必须是新栅格的文件路径。 如果名称以 /vsimem/ 开头,则栅格将在 GDAL 的虚拟文件系统中创建。
datatype

表示所有波段的数据类型的整数。 默认为 6 (Float32)。 新栅格的所有波段都需要具有相同的数据类型。 值映射为:

GDAL 像素类型

描述

1

GDT_Byte

八位无符号整数

2

GDT_UInt16

十六位无符号整数

3

GDT_Int16

十六位有符号整数

4

GDT_UInt32

三十二位无符号整数

5

GDT_Int32

三十二位有符号整数

6

GDT_Float32

三十二位浮点数

7

GDT_Float64

六十四位浮点数

nr_of_bands
表示栅格波段数的整数。 可以在创建时不传递波段数据的情况下创建栅格。 如果未指定波段数,则会根据 bands 输入的长度自动计算。 创建后无法更改波段数。
bands
带有波段输入数据的 band_input 字典列表。 生成的波段索引与提供的列表中的相同。 带输入字典的定义如下。 如果未提供波段数据,则栅格波段值将实例化为零数组,并且“无数据”值设置为 None
papsz_options

带有栅格创建选项的字典。 输入字典的键值对在创建栅格时传递给驱动程序。

可用选项是特定于驱动程序的,并在每个驱动程序的文档中进行了描述。

字典中的值不区分大小写,并在创建时自动转换为正确的字符串格式。

以下示例使用了一些可用于 GTiff 驱动程序 的选项。 结果是带有内部切片方案的压缩有符号字节栅格。 内部磁贴的块大小为 23 x 23:

>>> GDALRaster({
...    'driver': 'GTiff',
...    'name': '/path/to/new/file.tif',
...    'srid': 4326,
...    'width': 255,
...    'height': 255,
...    'nr_of_bands': 1,
...    'papsz_options': {
...        'compress': 'packbits',
...        'pixeltype': 'signedbyte',
...        'tiled': 'yes',
...        'blockxsize': 23,
...        'blockysize': 23,
...    }
... })


band_input 字典

ds_input 字典中的 bands 键是 band_input 字典的列表。 每个 band_input 字典都可以包含像素值和要在新栅格波段上设置的“无数据”值。 数据数组可以具有新栅格的完整大小或更小。 对于小于完整光栅的阵列,sizeshapeoffset 键控制像素值。 相应的键被传递给 data() 方法。 它们的功能与使用该方法设置波段数据相同。 下表描述了可以使用的键。

默认 如何使用
nodata_value None 映射到 nodata_value 属性
data nodata_value0 相同 传递给 data() 方法
size (with, height) 光栅 传递给 data() 方法
shape 大小一样 传递给 data() 方法
offset (0, 0) 传递给 data() 方法


使用 GDAL 的虚拟文件系统

GDAL 有一个基于内部内存的文件系统,它允许将内存块视为文件。 它可用于从二进制文件缓冲区读取和写入 GDALRaster 对象。

这在 Web 上下文中很有用,其中栅格可以作为缓冲区从远程存储获取或从视图返回而不写入磁盘。

GDALRaster 对象是在提供 bytes 对象作为输入时,或文件路径以 /vsimem/ 开头时在虚拟文件系统中创建的。

作为 bytes 提供的输入必须是文件的完整二进制表示。 例如:

# Read a raster as a file object from a remote source.
>>> from urllib.request import urlopen
>>> dat = urlopen('http://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'

要从头开始创建新的基于虚拟文件的栅格,请使用 ds_input 字典表示并提供以 /vsimem/ 开头的 name 参数(有关字典表示的详细信息,请参阅从数据创建栅格)。 对于基于虚拟文件的栅格,vsi_buffer 属性返回栅格的 bytes 表示。

以下是创建栅格并将其作为 HttpResponse 中的文件返回的方法:

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

配置

GDAL_LIBRARY_PATH

指定 GDAL 库位置的字符串。 通常,仅当 GDAL 库位于非标准位置(例如,/home/john/lib/libgdal.so)时才使用此设置。


例外

exception GDALException
基本 GDAL 异常,指示与 GDAL 相关的错误。
exception SRSException
构造或使用空间参考系统对象时发生错误时引发的异常。