GDAL stands for Geospatial Data Abstraction Library, and is a veritable "Swiss army knife" of GIS data functionality. A subset of GDAL is the OGR Simple Features Library, which specializes in reading and writing vector geographic data in a variety of standard formats.
GeoDjango は、ベクトル空間データの読み込みと座標変換、ラスタ (画像) データに関する GDAL の機能の最小限のサポートを含む、OGR の機能の一部に対する高レベルの Python インタフェースを提供します。
注釈
モジュールの名前は gdal ですが、GeoDjangoは現時点ではOGRの一部機能とGDALのラスタ機能のみをサポートしています。
The GDAL/OGR tools described here are designed to help you read in your geospatial data, in order for most of them to be useful you have to have some data to work with. If you're starting out and don't yet have any data of your own to use, GeoDjango tests contain a number of data sets that you can use for testing. You can download them here:
$ 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
DataSource¶DataSource is a wrapper for the OGR data source object that
supports reading data from a variety of OGR-supported geospatial file
formats and data sources using a consistent interface. Each
data source is represented by a DataSource object which contains
one or more layers of data. Each layer, represented by a Layer
object, contains some number of geographic features (Feature),
information about the type of features contained in that layer (e.g.
points, polygons, etc.), as well as the names and types of any
additional fields (Field) of data that may be associated with
each feature in that layer.
データソース DataSource のコンストラクタに必要なパラメータは1つだけです。しかし、OGRはデータベースを含むより複雑なデータソースもサポートしており、パスの代わりに特別な名前文字列を渡すことでアクセスできます。詳細は OGR Vector Formats のドキュメントを参照してください。 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¶Layer は DataSource オブジェクト内のデータ層をラップするためのものです。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
上記で読み込まれた都市データソースからの例の出力は、明らかに1つのレイヤーである "cities" を含み、その中に3つのポイント・フィーチャが含まれています。以下の例は簡略化のため、そのレイヤーを変数 layer に保存していると仮定しています:
>>> layer = ds[0]
データソース内のこのレイヤーの名前を返します。
>>> layer.name
'cities'
レイヤー内の要素数を返す。 len(layer) と同じです:
>>> layer.num_feat
3
レイヤーのジオメトリタイプを OGRGeomType オブジェクトとして返します:
>>> layer.geom_type.name
'Point'
レイヤーのフィールド数、つまりレイヤー内の各フィーチャに関連付けられたデータフィールドの数を返します:
>>> layer.num_fields
4
このレイヤーの各フィールド名のリストを返します:
>>> layer.fields
['Name', 'Population', 'Density', 'Created']
このレイヤーの各フィールドのデータ型のリストを返します。これらは、以下で説明される Field のサブクラスです:
>>> [ft.__name__ for ft in layer.field_types]
['OFTString', 'OFTReal', 'OFTReal', 'OFTDate']
このレイヤーの各フィールドの最大フィールド幅のリストを返します:
>>> layer.field_widths
[80, 11, 24, 10]
このレイヤーの各フィールドの数値精度のリストを返します。これは、数値フィールド以外では意味がありません (ゼロに設定されます):
>>> layer.field_precisions
[0, 0, 15, 0]
このレイヤーの空間範囲を Envelope オブジェクトとして返します。
>>> layer.extent.tuple
(-104.609252, 29.763374, -95.23506, 38.971823)
このレイヤーに関連付けられた 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 が返されます。
このレイヤーの空間フィルタを取得または設定するために使用できるプロパティ。空間フィルタは 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
レイヤーの各フィーチャについて、指定されたフィールドの値のリストを返すメソッド:
>>> layer.get_fields("Name")
['Pueblo', 'Lawrence', 'Houston']
レイヤ内の各フィーチャのジオメトリを含むリストを返すメソッドです。オプションの引数 geos が True に設定されている場合、ジオメトリは GEOSGeometry オブジェクトに変換されます。そうでない場合は、 OGRGeometry オブジェクトとして返されます:
>>> [pt.tuple for pt in layer.get_geoms()]
[(-104.609252, 38.255001), (-95.23506, 38.971823), (-95.363151, 29.763374)]
Returns a boolean indicating whether this layer supports the given
capability (a string). Examples of valid capability strings include:
'RandomRead', 'SequentialWrite', 'RandomWrite',
'FastSpatialFilter', 'FastFeatureCount', 'FastGetExtent',
'CreateField', 'Transactions', 'DeleteFeature', and
'FastSetNextByIndex'.
Feature¶Feature は OGR の feture をラップします。Feature オブジェクトを直接作成することはありません。代わりに、 Layer オブジェクトからそれらを取得します。各フィーチャにはジオメトリと追加のプロパティを含むフィールドのセットがあります。フィールドのジオメトリは、その geom プロパティを介してアクセスできます。これは OGRGeometry オブジェクトを返します。Feature はフィールドの標準的な Python コンテナのように動作し、そのフィールドを Field オブジェクトとして返します。フィールドはインデックスまたは名前を指定して直接アクセスすることもでき、あるいは for ループ内でフィーチャのフィールドをイテレートすることもできます。
このフィーチャのジオメトリを OGRGeometry オブジェクトとして返します:
>>> city.geom.tuple
(-104.609252, 38.255001)
このフィーチャーのフィールド (名前で指定) の値を返すメソッドで、 Field ラッパーオブジェクト ではありません :
>>> city.get("Population")
102121
このフィーチャのジオメトリの種類を、 OGRGeomType オブジェクトとして返します。これは同じレイヤー内のすべてのフィーチャに対して同じであり、フィーチャの元となる Layer オブジェクトの Layer.geom_type プロパティと同等です。
フィーチャに関連付けられたデータフィールドの数を返します。これは、特定のレイヤー内のすべてのフィーチャに対して同じであり、フィーチャが属する Layer オブジェクトの Layer.num_fields プロパティと同等です。
フィーチャに関連付けられたデータのフィールド名のリストを返します。これは与えられたレイヤーの全てのフィーチャで同じであり、フィーチャの元となった Layer オブジェクトの Layer.fields プロパティと等価です。
レイヤー内のフィーチャ識別子を返します:
>>> city.fid
0
そのフィーチャーの Layer の名前を返します。これは指定されたレイヤーの全てのフィーチャに対して同じです:
>>> city.layer_name
'cities'
与えられたフィールド名のインデックスを返すメソッド。これは、指定されたレイヤーのすべてのフィーチャで同じになります:
>>> city.index("Population")
1
Field¶このフィールドの名前を返します:
>>> city["Name"].name
'Name'
このフィールドの OGR 型を整数で返します。 FIELD_CLASSES 辞書は、これらの値を Field のサブクラスにマッピングします:
>>> city["Density"].type
2
このフィールドのデータ型の名前を文字列で返します:
>>> city["Name"].type_name
'String'
このフィールドの値を返します。Field クラス自体は値を文字列として返しますが、各サブクラスは最適な形式で値を返します:
>>> city["Population"].value
102121
このフィールドの幅を返します:
>>> city["Name"].width
80
このフィールドの数値精度を返します。非数値フィールドでは意味がありません (ゼロに設定されます) 。
>>> city["Density"].precision
15
フィールドの値を double (float) で返します:
>>> city["Density"].as_double()
874.7
フィールドの値を整数で返します:
>>> city["Population"].as_int()
102121
フィールドの値を文字列で返します:
>>> city["Name"].as_string()
'Pueblo'
フィールドの値を、日付と時刻のタプルで返します:
>>> 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¶Driver クラスは、OGR DataSource ドライバをラップするために内部で使用されます。
現在登録されているOGRベクタードライバーの数を返します。
OGRGeometry¶OGRGeometry オブジェクトは GEOSGeometry オブジェクトと同様の機能を持ち、OGR の内部ジオメトリ表現の薄いラッパーです。そのため、 DataSource を使用する場合、より効率的にデータにアクセスできます。GEOSとは異なり、OGRGeometry は空間参照系 (spatial reference system) と座標変換をサポートしています。
>>> from django.contrib.gis.gdal import OGRGeometry
>>> polygon = OGRGeometry("POLYGON((0 0, 5 0, 5 5, 0 5))")
このオブジェクトは OGR Geometry クラスのラッパーです。このオブジェクトは、与えられた geom_input パラメータから直接インスタンス化されます。このパラメータには、WKT、HEX、GeoJSON、WKB データを含む buffer、または OGRGeomType オブジェクトを含む文字列を指定できます。これらのオブジェクトは Layer (DataSource の一部) からベクトルデータを読み込む際に、 Feature.geom 属性からも返されます。
与えられたGML文字列から OGRGeometry を構築する。
与えられたバウンディングボックス (4要素タプル) から Polygon を構築します。
LineString の点の数、Polygon のリングの数、または GeometryCollection のジオメトリの数を返します。他のジオメトリタイプには適用されません。
LineString の点、Polygon のリング、または GeometryCollection のジオメトリをイテレートします。他のジオメトリタイプには適用されません。
LineString オブジェクトにおける指定されたインデックスのポイント、Polygon オブジェクトにおける指定されたインデックスの内部リング、または GeometryCollection における指定されたインデックスのジオメトリを返します。他のジオメトリタイプには適用されません。
ジオメトリの座標次元数を返します。例えば、点なら0、線なら1、その他も同様です。
>>> polygon.dimension
2
このジオメトリが Z 座標を持っているかどうかを示す真偽値。
Z 座標次元を追加または削除するメソッド。
>>> p = OGRGeometry("POINT (1 2 3)")
>>> p.is_3d
True
>>> p.set_3d(False)
>>> p.wkt
"POINT (1 2)"
このジオメトリが M 座標を持っているかどうかを示す真偽値。
M 座標次元を追加または削除するメソッド。
>>> p = OGRGeometry("POINT (1 2)")
>>> p.is_measured
False
>>> p.set_measured(True)
>>> p.wkt
"POINT M (1 2 0)"
このジオメトリの要素数を返します:
>>> polygon.geom_count
1
このジオメトリがカーブジオメトリであるか、またはカーブジオメトリを含んでいるかどうかを示すブール値です。
ジオメトリの線形バージョンを返します。変換ができない場合は、元のジオメトリがそのまま返されます。
ジオメトリの曲線バージョンを返します。変換できない場合は、元のジオメトリがそのまま返されます。
このジオメトリを記述するのに使用されたポイントの数を返します:
>>> polygon.point_count
4
point_count のエイリアス。
point_count のエイリアス。
このジオメトリのタイプを OGRGeomType オブジェクトとして返します。
このジオメトリの型名を返します:
>>> polygon.geom_name
'POLYGON'
このジオメトリの面積を返すか、面積を含まないジオメトリの場合は0を返します。
>>> polygon.area
25.0
このジオメトリの外接矩形を Envelope オブジェクトとして返します。
このジオメトリの外接矩形を Envelope オブジェクトではなく、4値タプルとして返します。
>>> point.extent
(0.0, 0.0, 5.0, 5.0)
This property controls the spatial reference for this geometry, or
None if no spatial reference system has been assigned to it.
If assigned, accessing this property returns a SpatialReference
object. It may be set with another SpatialReference object,
or any input that SpatialReference accepts. Example:
>>> city.geom.srs.name
'GCS_WGS_1984'
Returns or sets the spatial reference identifier corresponding to
SpatialReference of this geometry. Returns None if
there is no spatial reference information associated with this
geometry, or if an SRID cannot be determined.
このジオメトリに対応する GEOSGeometry オブジェクトを返します。
このジオメトリの GML 形式の文字列表現を返します:
>>> OGRGeometry("POINT(1 2)").gml
'<gml:Point><gml:coordinates>1,2</gml:coordinates></gml:Point>'
このジオメトリのHEX WKB形式での文字列表現を返します:
>>> OGRGeometry("POINT(1 2)").hex
'0101000000000000000000F03F0000000000000040'
このジオメトリの文字列表現を JSON 形式で返します:
>>> OGRGeometry("POINT(1 2)").json
'{ "type": "Point", "coordinates": [ 1.000000, 2.000000 ] }'
このジオメトリを KML 形式の文字列で返します。
このジオメトリの WKB 表現を保持するために必要な WKB バッファのサイズを返します:
>>> OGRGeometry("POINT(1 2)").wkb_size
21
このジオメトリのWKB表現を含む buffer を返します。
このジオメトリの文字列表現を WKT 形式で返します。
このジオメトリの EWKT 表現を返します。
このジオメトリオブジェクトの新しい OGRGeometry クローンを返します。
このジオメトリ内に閉じていないリングがある場合、このルーチンは始点を終点に追加することで閉じます:
>>> triangle = OGRGeometry("LINEARRING (0 0,0 1,1 0)")
>>> triangle.close_rings()
>>> triangle.wkt
'LINEARRING (0 0,0 1,1 0,0 0)'
このジオメトリを異なる空間参照系に変換します。CoordTransform オブジェクト、SpatialReference オブジェクト、または SpatialReference が受け付けるその他の入力 (空間参照系の WKT や PROJ 文字列、整数の SRID など) を受け付けます。
デフォルトでは何も返されず、ジオメトリはインプレースで変換されます。しかし、 clone キーワードが True に設定されている場合は、このジオメトリの変換されたクローンが返されます。
このジオメトリがもう一方のジオメトリと交差している場合は True を返し、そうでない場合は False を返します。
このジオメトリがもう一方のジオメトリと等価であれば True を返し、そうでなければ False を返します。
このジオメトリがもう一方のジオメトリと空間的に不連続であれば (つまり交差していなければ) True を返し、そうでなければ False を返します。
このジオメトリがもう一方のジオメトリに接している場合は True を返し、接していない場合は False を返します。
このジオメトリがもう一方のジオメトリと交差している場合は True を返し、交差していない場合は False を返します。
このジオメトリがもう一方のジオメトリに含まれる場合は True を返し、そうでない場合は False を返します。
このジオメトリがもう一方のジオメトリを含んでいる場合は True を返し、そうでない場合は False を返します。
このジオメトリがもう一方のジオメトリと重なっている場合は True を返し、重なっていない場合は False を返します。
このジオメトリの境界を、新しい OGRGeometry オブジェクトとして返します。
このジオメトリを含む最小の凸ポリゴンを、新しい OGRGeometry オブジェクトとして返します。
このジオメトリともう一方のジオメトリの差分からなる領域を、新しい OGRGeometry オブジェクトとして返します。
このジオメトリともう一方のジオメトリの交点からなる領域を、新しい OGRGeometry オブジェクトとして返します。
このジオメトリともう一方のジオメトリの対称差からなる領域を、新しい OGRGeometry オブジェクトとして返します。
このジオメトリともう一方のジオメトリの和からなる領域を、新しい OGRGeometry オブジェクトとして返します。
このジオメトリの重心を表す Point を返します。
点ジオメトリの座標をタプルとして、線ジオメトリの座標をタプルのタプルとして返す:
>>> OGRGeometry("POINT (1 2)").tuple
(1.0, 2.0)
>>> OGRGeometry("LINESTRING (1 2,3 4)").tuple
((1.0, 2.0), (3.0, 4.0))
tuple のエイリアス。
この点のX座標を返します:
>>> OGRGeometry("POINT (1 2)").x
1.0
この点のY座標を返します:
>>> OGRGeometry("POINT (1 2)").y
2.0
この点の Z 座標を返し、点が Z 座標を持たない場合は None を返します:
>>> OGRGeometry("POINT (1 2 3)").z
3.0
この点の M 座標を返します。もし点に M 座標がない場合は None を返します。
>>> OGRGeometry("POINT ZM (1 2 3 4)").m
4.0
この行のX座標のリストを返します:
>>> OGRGeometry("LINESTRING (1 2,3 4)").x
[1.0, 3.0]
この行のY座標のリストを返します:
>>> OGRGeometry("LINESTRING (1 2,3 4)").y
[2.0, 4.0]
このラインの Z 座標のリストを返すか、Z 座標を持たない場合は None を返します:
>>> OGRGeometry("LINESTRING (1 2 3,4 5 6)").z
[3.0, 6.0]
この線分の M 座標のリストを返します。もし線分に M 座標がない場合は None を返します。
>>> OGRGeometry("LINESTRING(0 1 2 10, 1 2 3 11, 2 3 4 12)").m
[10.0, 11.0, 12.0]
OGRGeomType¶このクラスは、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
OGR Geometry の型の省略文字列形式を返します:
>>> gt1.name
'Polygon'
OGR ジオメトリタイプに対応する数値を返します:
>>> gt1.num
3
この OGR 型を格納するために使用する Django フィールド型 (GeometryField のサブクラス) を返し、適切な Django 型がない場合は None を返します:
>>> gt1.django
'PolygonField'
Envelope¶SpatialReference¶空間参照系 (spatial reference) オブジェクトは与えられた srs_input に基づいて初期化されます。この srs_input には以下のいずれかが指定されることがあります。
OGC Well Known Text (WKT) (文字列)
EPSG コード (整数または文字列)
PROJ 文字列
よく知られている標準の略記文字列 ('WGS84', 'WGS72', 'NAD27', 'NAD83')
例:
>>> 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
与えられた文字列属性ノードの値を返します。ノードが存在しない場合は 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
指定されたターゲットノードの属性値 (例: 'PROJCS')。index キーワードは返す子ノードのインデックスを指定します。
指定された文字列の対象ノードに対するオーソリティ名を返します。
与えられた文字列の対象ノードのオーソリティコードを返します。
この空間参照系 (SpatialReference) オブジェクトのクローンを返します。
このメソッドは SpatialReference の WKT を検査し、EPSG 識別子が該当する場合は EPSG オーソリティノードを追加します。
この SpatialReference を ESRI のフォーマットから EPSG に変換します。
この SpatialReference を ESRI のフォーマットに変換します。
与えられた空間参照系が有効かどうかをチェックし、有効でない場合は例外が発生します。
EPSGコードから空間参照系をインポートします。
PROJ文字列から空間参照系をインポートします。
WKT から空間参照系をインポートします。
XML から空間参照系をインポートします。
この空間参照系の名前を返します。
トップレベルのオーソリティの SRID を返すか、未定義の場合は None を返します。
線形単位の名前を返します。
線形単位の値を返します。
角度の単位の名前を返します。
角度の単位の値を返します。
線形単位を返すか角度単位を返すかを自動的に決定し、単位値と単位名の2タプルを返します。
この空間参照系の楕円体パラメータのタプル (半長軸、半短軸、逆扁平率) を返します。
この空間参照系に対する楕円体の長軸を返します。
この空間参照系における楕円体の短半径を返します。
この空間参照系の楕円体の逆扁平率を返します。
この空間参照系が地理座標系 (親ノードが GEOGCS) である場合に True を返します。
この空間参照系がローカルである場合 (ルートノードが LOCAL_CS である場合) には True を返します。
この空間参照系が投影座標系である場合 (ルートノードが PROJCS である場合) に True を返します。
この空間参照系の WKT 表現を返します。
WKT のきれい ("pretty") な表現を返します。
この空間参照系の PROJ 表現を返します。
SpatialReference.proj のエイリアス。
この空間参照系の XML 表現を返します。
CoordTransform¶座標系の変換を表します。ソースとターゲットの座標系をそれぞれ表す2つの 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 is a wrapper for the GDAL raster source object that
supports reading data from a variety of GDAL-supported geospatial file formats
and data sources using a consistent interface. Each data source is represented
by a GDALRaster object which contains one or more layers of data named
bands. Each band, represented by a GDALBand object, contains
georeferenced image data. For example, an RGB image is represented as three
bands: one for red, one for green, and one for blue.
注釈
ラスターデータにおいては、ラスターインスタンスとそのデータソースとの間に差異はありません。ジオメトリオブジェクトと異なり、GDALRaster オブジェクトは常にデータソースです。一時的なラスターは、対応するドライバーを使用してメモリ内にインスタンス化できますが、ファイルベースのラスターソースと同じクラスになります。
GDALRaster のコンストラクターは 2 つのパラメーターを受け取ります。最初のパラメーターはラスターのソースを定義し、2 番目のパラメーターはラスターを書き込みモードで開くかどうかを定義します。新しく作成されたラスターに対しては、2 番目のパラメーターは無視され、常に書き込みモードで新しいラスターが作成されます。
最初のパラメータは、３つの形式を取ることができます。ファイルパスを表す文字列あるいは Path オブジェクト (ファイルシステムまたはGDAL仮想ファイルシステム) 、新しいラスタを定義する値が入った辞書、ラスタファイルを表すバイトオブジェクトです。
入力がファイルパスの場合は、そのパスからラスタがオープンされます。入力が辞書の生データの場合、パラメータ width、height、srid が必要です。入力がバイトオブジェクトの場合、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'
入力ファイルのパスまたはインスタンス化の際に指定された名前に相当するソースの名前。
>>> GDALRaster({"width": 10, "height": 10, "name": "myraster", "srid": 4326}).name
'myraster'
入力ファイルを処理するために使用されるGDALドライバの名前です。ファイルから作成された GDALRaster では、ドライバタイプが自動的に検出されます。ゼロからラスターを作成する場合、デフォルトではメモリ内のラスター ('MEM') ですが、必要に応じて変更できます。例えば、 GeoTiff ファイルには GTiff を使用します。ファイルタイプの一覧については、 GDAL Raster Formats リストも参照してください。
インメモリのラスタは次の例のように作成できます:
>>> 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'
ラスターの空間参照系を SpatialReference インスタンスとして返します。SRSは他の SpatialReference に設定したり、 SpatialReference コンストラクタが受け付ける任意の入力を与えることで変更できます。
>>> rst = GDALRaster({"width": 10, "height": 20, "srid": 4326})
>>> rst.srs.srid
4326
>>> rst.srs = 3086
>>> rst.srs.srid
3086
ラスタの空間参照システム識別子(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
以下のリレーションシップを使用して、ピクセル/ライン座標をジオリファレンス空間にマッピングする6つの係数のタプルとして表現されるアフィン変換行列:
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]
ソースの空間参照システム内で、ラスターの左上原点の座標を、x と y のメンバを持つ点オブジェクトとして表します。
>>> rst = GDALRaster({"width": 10, "height": 20, "srid": 4326})
>>> rst.origin
[0.0, 0.0]
>>> rst.origin.x = 1
>>> rst.origin
[1.0, 0.0]
ラスタのジオリファレンスに使用されるピクセルの幅と高さ。 x と y のメンバーを持つポイントオブジェクトとして表されます。詳細については、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]
ラスターを x と y のメンバを持つ点オブジェクトとしてジオリファレンスするために使用されるスキュー係数。ノースアップ画像の場合、これらの係数は両方とも 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]
ラスターソースの範囲 (境界値)。ソースの空間参照システムにおける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)
ソースの全てのバンドを、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)
このラスターの変形されたバージョンを返します。
変化系パラメータは ds_input 引数で指定することができる。 ds_input の使用方法はクラスコンストラクタの対応する引数と同様です。 ds_input はターゲットとなるラスターの特性を辞書に格納したものです。指定できる辞書のキーは width, height, SRID, origin, scale, skew, datatype, driver, name (filename) です。
デフォルトでは、warp 関数はほとんどのパラメータを元のソースラスタの値と等しく保ちます。したがって、変更する必要のあるパラメータのみを指定すればよいです。ただし、これにはドライバも含まれるので、ファイルベースのラスタの場合、warp 関数は新しいラスタをディスク上に作成します。
ソースラスタと異なるパラメータは名前のみです。ラスタ名のデフォルト値は、ソースラスタの名前に '_copy' + source_driver_name を追加したものです。ファイルベースのラスタの場合は、ターゲットとなるラスタのファイルパスを指定することを推奨します。
変形に使用するリサンプリングアルゴリズムは resampling 引数で指定します。デフォルトは NearestNeighbor で、他に指定できる値は Bilinear、Cubic、CubicSpline、Lanczos、Average、Mode です。
max_error 引数は、変換の近似において許容される入力ピクセルでの最大誤差を指定するために使用できます。デフォルト値は正確な計算のために 0.0 です。
GDAL に詳しいユーザーにとって、この関数は gdalwarp コマンドラインユーティリティと同様の機能を持っています。
たとえば、warp 関数は、ラスターを元のピクセルスケールの2倍に集約するために使用できます:
>>> 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)
このラスタを異なる空間参照システム (srs) に変換します。これは SpatialReference オブジェクトであり、または SpatialReference で受け入れられる他の入力 (空間参照 WKT や PROJ 文字列、整数の SRID を含む) でも構いません。
新しい空間参照系で現在のラスタの境界とスケールを計算し、 warp 関数を使用してラスタを変形します。
デフォルトでは、コピー元のラスタのドライバが使用され、ラスタ名はコピー元の名前に '_copy' + source_driver_name を追加したものになります。異なるドライバや名前を指定するには、driver および name 引数を使用します。
デフォルトのリサンプリングアルゴリズムは 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_srs = SpatialReference(4326)
>>> target = rst.transform(target_srs)
>>> target.origin
[-82.98492744885776, 27.601924753080144]
このラスタのメタデータ。ネストした辞書として表現されます。第1レベルのキーはメタデータ・ドメインです。第2レベルには、各ドメインのメタデータ項目名と値が含まれます。
メタデータ項目を設定または更新するには、上記の入れ子構造を使用して、対応するメタデータ項目をメソッドに渡します。指定された辞書にあるキーのみが更新され、メタデータの残りの部分は変更されません。
メタデータ項目を削除するには、メタデータの値として 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'}}
GDALBand¶GDALBand のインスタンスは明示的に作成されるのではなく、代わりに GDALRaster オブジェクトから、その bands 属性を介して取得されます。 GDALBand にはラスターの実際のピクセル値が含まれています。
バンドの名前または説明 (あれば) 。
ピクセル単位でのバンドの幅 (X軸)。
ピクセル単位のバンドの高さ (Y軸)。
このバンドに含まれるピクセルの総数。これは width * height に等しいです。
このバンドのピクセル値の統計量を計算します。戻り値は (minimum, maximum, mean, standard deviation) という構造を持つタプルです。
もし approximate 引数が True に設定されている場合、統計はオーバービューまたは画像タイルのサブセットに基づいて計算されます。
refresh 引数を True に設定すると、統計情報はデータから直接計算され、その結果でキャッシュが更新されます。
永続的なキャッシュ値が見つかった場合、その値が返されます。永続的補助メタデータ (PAM) サービスを使用するラスターフォーマットの場合、統計情報は補助ファイルにキャッシュされる可能性があります。場合によっては、このメタデータがピクセル値と同期していなかったり、 approximate 引数の値を反映しない以前の呼び出しの値が返されたりすることがあります。そのような場合は、 refresh 引数を使用して更新された値を取得し、キャッシュに保存します。
空のバンド (すべてのピクセル値が "no data "である) の場合、すべての統計は None として返されます。
バンドの最小ピクセル値 ("no data" 値を除く) 。
バンドの最大ピクセル値 ("no data" 値を除く) 。
バンドのすべてのピクセル値の平均 ("no data"値を除く) 。
バンドのすべてのピクセル値の標準偏差 ("no data" の値を除く)。
バンドの "no data" 値は通常、有効なデータではないピクセルを示す特別なマーカー値です。これらのピクセルは一般的に表示されるべきではなく、分析操作にも寄与すべきではありません。
既存の "no data" 値を削除するには、このプロパティを None に設定します。
バンドに含まれるデータ型。0 (Unknown) から 14 までの整数定数。 as_string が True の場合、データ型は文字列として返されます。指定可能な値については、 datatype の値テーブル の "GDAL ピクセルタイプ" カラムを参照してください。
0 から 16 までの整数値。 as_string が True の場合、データ型は文字列として返されます: 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, GCI_YCbCr_CrBand 。 GCI_YCbCr_CrBand は GCI_Max も表します。どちらも 16 の整数に対応しますが、文字列として返されるのは GCI_YCbCr_CrBand だけです。
GDALBand のピクセル値へのアクセサ。パラメータを指定しない場合は完全なデータ配列を返す。オフセットとブロックサイズをタプルで指定することで、ピクセル配列のサブセットをリクエストできます。
NumPyが利用可能な場合、データはNumPy配列として返されます。パフォーマンス上の理由から、NumPyの利用を強くお勧めします。
data パラメータが指定された場合、データは GDALBand に書き込まれます。入力には、パックされた文字列、バッファ、リスト、配列、NumPyの配列のいずれかを指定できます。入力に含まれるアイテムの数は、通常、バンド内の全ピクセル数に対応するか、 offset と size パラメータが指定されている場合は、ピクセル値の特定のブロックのピクセル数に対応します。
入力のアイテム数がターゲットピクセルブロックと異なる場合は、 shape パラメータを指定する必要があります。shapeは入力データの幅と高さをピクセル単位で指定するタプルです。その後、選択されたブロックのピクセル値を更新するためにデータが複製されます。これは、たとえば1つの値でバンド全体を塗りつぶすのに便利です。
例:
>>> 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)
このバンドのメタデータ。機能は GDALRaster.metadata と同じです。
このセクションでは、ds_input パラメータを使用してゼロからラスタを作成する方法について説明します。
GDALRaster コンストラクタに dict を渡すことで新しいラスタが作成されます。この辞書には、新しいラスタの定義パラメータ (原点、サイズ、空間参照システムなど) が含まれています。さらに、辞書にはピクセルデータや新しいラスタのフォーマットに関する情報も含めることができます。したがって、作成されるラスタは指定したドライバに応じてファイルベースまたはメモリベースとなります。
ラスターデータを辞書や JSON で記述する標準はありません。そのため、 GDALRaster クラスへの辞書入力の定義は Django 固有のものです。これは geojson フォーマットに触発されたものですが、 geojson 標準は今のところベクタフォーマットに限定されています。
GDALRaster と GDALBand クラスの対応する属性とメソッドのドキュメントに、ラスタを作成する際に異なるキーを使用する例があります。
ds_input 辞書¶ds_input 辞書内でラスターを作成するために必要なキーはわずかで、 width、height、および srid です。その他のパラメータにはデフォルト値が設定されています(以下の表を参照)。 ds_input 辞書に渡せるキーのリストは、GDALRaster プロパティに密接に関連していますが、完全には一致していません。多くのパラメータはこれらのプロパティに直接マッピングされており、その他のパラメータについては以下で説明します。
以下の表は、ds_input 辞書に設定できるすべてのキーを説明しています。
キー |
デフォルト値 |
使い方 |
|---|---|---|
|
必須 |
|
|
必須 |
|
|
必須 |
|
|
|
|
|
|
下記参照 |
|
|
|
|
|
|
|
|
|
|
|
下記参照 |
|
|
下記参照 |
|
|
下記参照 |
|
|
下記参照 |
ラスタの名前を表す文字列。ファイルベースのラスタを作成する場合、このパラメータは新しいラスタのファイルパスでなければなりません。名前が /vsimem/ で始まる場合、ラスタはGDALの仮想ファイルシステムに作成されます。
全バンドのデータ型を表す整数です。デフォルトは 6 (Float32) です。新しいラスターの全てのバンドは同じデータ型である必要があります。値のマッピングは以下の通りです:
値 |
GDAL ピクセルタイプ |
説明 |
|---|---|---|
1 |
GDT_Byte |
8 ビット符号なし整数 |
2 |
GDT_UInt16 |
16 ビット符号なし整数 |
3 |
GDT_Int16 |
16ビット符号付き整数 |
4 |
GDT_UInt32 |
32ビット符号なし整数 |
5 |
GDT_Int32 |
32ビット符号付き整数 |
6 |
GDT_Float32 |
32ビット浮動小数点 |
7 |
GDT_Float64 |
64ビット浮動小数点 |
12 |
GDT_UInt64 |
64ビット符号なし整数 (GDAL 3.5+) |
13 |
GDT_Int64 |
64ビット符号付き整数 (GDAL 3.5+) |
14 |
GDT_Int8 |
8ビット符号付き整数 (GDAL 3.7+) |
ラスターのバンド数を表す整数。作成時にバンドデータを渡さずにラスタを作成することもできます。バンド数が指定されていない場合は、入力 bands の長さから自動的に計算されます。作成後にバンド数を変更することはできません。
バンド入力データを持つ band_input 辞書のリスト。結果のバンドインデックスは、提供されたリストと同じです。バンド入力辞書の定義は以下に示します。バンドデータが提供されない場合、ラスターバンドの値はゼロの配列としてインスタンス化され、"no data" の値は None に設定されます。
ラスターを作成するためのオプションを持つ辞書。入力辞書のキーと値のペアは、ラスターの作成時にドライバーに渡されます。
利用可能なオプションは、ドライバごとのドキュメントに記載されています。
辞書内の値は大文字と小文字を区別しないため、作成時に適切な文字列形式に自動変換されます。
以下の例では、GTiff driver で利用可能なオプションのいくつかを使用しています。結果は、内部タイルスキームを持つ圧縮されたラスターです。内部タイルのブロックサイズは23×23です:
>>> GDALRaster(
... {
... "driver": "GTiff",
... "name": "/path/to/new/file.tif",
... "srid": 4326,
... "width": 255,
... "height": 255,
... "nr_of_bands": 1,
... "papsz_options": {
... "compress": "packbits",
... "tiled": "yes",
... "blockxsize": 23,
... "blockysize": 23,
... },
... }
... )
band_input 辞書¶ds_input 辞書内の bands キーは band_input 辞書のリストです。各 band_input 辞書には、新しいラスタのバンドに設定されるピクセル値と "no data" 値が含まれている可能性があります。データ配列は新しいラスタの全サイズを持つことも、より小さくなることもあります。新しいラスタより小さい配列の場合、size、shape、および offset キーがピクセル値を制御します。これらの対応するキーは、data() メソッドに渡されます。これらの機能は、そのメソッドでバンドデータを設定するのと同じです。以下の表は使用可能なキーを説明しています。
キー |
デフォルト値 |
使い方 |
|---|---|---|
|
|
|
|
|
|
|
ラスターの |
|
|
size と同じ |
|
|
|
|
GDALはファイルシステムに保存されたファイルにアクセスするだけでなく、圧縮、暗号化、またはリモートファイルなど他の種類のファイルにアクセスを抽象化するための仮想ファイルシステムもサポートしています。
GDALには、メモリをファイルとして扱う内部メモリベースのファイルシステムがあります。これを使用して、 GDALRaster オブジェクトをバイナリファイルバッファーに読み込んだり書き込んだりすることができます。
これは、ラスターがリモートストレージからバッファとして取得されたり、ディスクに書き込まれることなくビューから返されたりするようなWebコンテキストで便利です。
bytes オブジェクトが入力として提供されるか、ファイルパスが /vsimem/ で始まる場合、仮想ファイルシステム内に GDALRaster オブジェクトが作成されます。
入力は、ファイルの bytes としての完全なバイナリ表現である必要があります。例えば、:
# 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'
新しい仮想ファイルベースのラスタをゼロから作成するには、 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 のローカルビルドによっては、他の仮想ファイルシステムもサポートされているかもしれません。提供されたパスの先頭に適切な /vsi*/ 接頭辞を付けることで使用できます。詳細は GDAL Virtual Filesystems documentation を参照してください。
ファイルを解凍して結果のラスタをインスタンス化する代わりに、GDAL は /vsizip/, /vsigzip/, /vsitar/ 仮想ファイルシステムを使用して圧縮ファイルに直接アクセスできます:
>>> 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")
GDALは、それらの機能がビルトインされている限り、オンラインリソースとストレージプロバイダを透過的にサポートできます。
認証なしで公開されたラスターファイルにアクセスするには、/vsicurl/ を使用できます。
>>> from django.contrib.gis.gdal import GDALRaster
>>> rst = GDALRaster("/vsicurl/https://example.com/raster.tif")
>>> rst.name
'/vsicurl/https://example.com/raster.tif'
商用のストレージプロバイダ (/vsis3/ など) の場合、認証やその他の設定を事前に行っておく必要があります (利用可能なオプションについては GDAL Virtual Filesystems documentation を参照してください) 。
GDAL_LIBRARY_PATH¶A string specifying the location of the GDAL library. Typically,
this setting is only used if the GDAL library is in a non-standard
location (e.g., /home/john/lib/libgdal.so).
12月 03, 2025