GDAL は Geospatial Data Abstraction Library の略であり、GIS データ機能の真の「スイス・アーミーナイフ」とも言える存在です。GDAL のサブセットには、標準フォーマットでのベクトル地理データの読み書きに特化した OGR Simple Features ライブラリがあります。
\nGeoDjango は、ベクトル空間データの読み込みと座標変換、ラスタ (画像) データに関する GDAL の機能の最小限のサポートを含む、OGR の機能の一部に対する高レベルの Python インタフェースを提供します。
\n\nここで説明する GDAL/OGR ツールは、地理空間データの読み込みを支援するように設計されています。もしあなたが使い始めたばかりで、まだ自分で使えるデータを持っていないのなら、GeoDjango のテストにはテストに使えるデータセットがたくさん含まれています。以下からダウンロードできます:
\n$ wget https://raw.githubusercontent.com/django/django/main/tests/gis_tests/data/cities/cities.{shp,prj,shx,dbf}\n$ wget https://raw.githubusercontent.com/django/django/main/tests/gis_tests/data/rasters/raster.tif\nDataSourceLink to this heading#DataSource は OGR データソースオブジェクトのラッパーであり、OGR がサポートする様々な地理空間ファイルフォーマットやデータソースから一貫したインターフェースを使用してデータを読み込むことをサポートします。各データソースは DataSource オブジェクトで表現され、1つ以上のデータのレイヤーを含みます。各レイヤーは Layer オブジェクトで表現され、いくつかの地理情報フィーチャ (Feature) 、そのレイヤーに含まれるフィーチャ (feture) の種類 (ポイント、ポリゴンなど) に関する情報、およびそのレイヤー内の各フィーチャに関連付けられた追加フィールド (Field) の名前とデータ型を含みます。
データソース DataSource のコンストラクタに必要なパラメータは1つだけです。しかし、OGRはデータベースを含むより複雑なデータソースもサポートしており、パスの代わりに特別な名前文字列を渡すことでアクセスできます。詳細は OGR Vector Formats のドキュメントを参照してください。 DataSource インスタンスの name プロパティは、使用するデータソースのOGR名を指定します。
オプションの encoding パラメータで、ソース中の文字列の非標準エンコードを指定できます。これは、フィールドの値を読み込む際に DjangoUnicodeDecodeError 例外が発生する場合に便利です。
データソース DataSource を作成したら、 layer_count プロパティにアクセスするか、len() 関数を使用することで、データのレイヤー数を知ることができます。データのレイヤー自体へのアクセスについては、次のセクションを参照してください:
>>> from django.contrib.gis.gdal import DataSource\n>>> ds = DataSource("/path/to/your/cities.shp")\n>>> ds.name\n'/path/to/your/cities.shp'\n>>> ds.layer_count # This file only contains one layer\n1\nデータソース内のレイヤーの数を返します。
\nデータソースの名前を返します。
\nLayerLink to this heading#Layer は DataSource オブジェクト内のデータ層をラップするためのものです。Layer オブジェクトを直接作成することはありません。代わりに、それらを DataSource オブジェクトから取得します。これは、本質的には Layer オブジェクトの標準的なPythonコンテナです。例えば、インデックスによって特定のレイヤーにアクセスできます (例: 最初のレイヤーにアクセスするには ds[0]) 、または for ループでコンテナ内のすべてのレイヤーをイテレートできます。Layer 自体が幾何学的な特徴のためのコンテナとして機能します。
通常、特定のレイヤー内のすべてのフィーチャは同じジオメトリタイプを持ちます。レイヤーの geom_type プロパティはフィーチャのタイプを識別する OGRGeomType であり、DataSource 内の各レイヤーに関する基本情報を出力するために使用できます。
>>> for layer in ds:\n... print('Layer "%s": %i %ss' % (layer.name, len(layer), layer.geom_type.name))\n...\nLayer "cities": 3 Points\n上記で読み込まれた都市データソースからの例の出力は、明らかに1つのレイヤーである "cities" を含み、その中に3つのポイント・フィーチャが含まれています。以下の例は簡略化のため、そのレイヤーを変数 layer に保存していると仮定しています:
>>> layer = ds[0]\nデータソース内のこのレイヤーの名前を返します。
\n>>> layer.name\n'cities'\nレイヤー内の要素数を返す。 len(layer) と同じです:
>>> layer.num_feat\n3\nレイヤーのジオメトリタイプを OGRGeomType オブジェクトとして返します:
>>> layer.geom_type.name\n'Point'\nレイヤーのフィールド数、つまりレイヤー内の各フィーチャに関連付けられたデータフィールドの数を返します:
\n>>> layer.num_fields\n4\nこのレイヤーの各フィールド名のリストを返します:
\n>>> layer.fields\n['Name', 'Population', 'Density', 'Created']\nこのレイヤーの各フィールドのデータ型のリストを返します。これらは、以下で説明される Field のサブクラスです:
>>> [ft.__name__ for ft in layer.field_types]\n['OFTString', 'OFTReal', 'OFTReal', 'OFTDate']\nこのレイヤーの各フィールドの最大フィールド幅のリストを返します:
\n>>> layer.field_widths\n[80, 11, 24, 10]\nこのレイヤーの各フィールドの数値精度のリストを返します。これは、数値フィールド以外では意味がありません (ゼロに設定されます):
\n>>> layer.field_precisions\n[0, 0, 15, 0]\nこのレイヤーの空間範囲を Envelope オブジェクトとして返します。
>>> layer.extent.tuple\n(-104.609252, 29.763374, -95.23506, 38.971823)\nこのレイヤーに関連付けられた SpatialReference を返すプロパティ:
>>> print(layer.srs)\nGEOGCS["GCS_WGS_1984",\n DATUM["WGS_1984",\n SPHEROID["WGS_1984",6378137,298.257223563]],\n PRIMEM["Greenwich",0],\n UNIT["Degree",0.017453292519943295]]\nLayer に空間参照情報が関連付けられていない場合、 None が返されます。
このレイヤーの空間フィルタを取得または設定するために使用できるプロパティ。空間フィルタは OGRGeometry インスタンス、4タプル、または None でのみ設定できます。None 以外で設定された場合、レイヤーをイテレートするときにはフィルタと交差するフィーチャのみが返されます。
>>> print(layer.spatial_filter)\nNone\n>>> print(len(layer))\n3\n>>> [feat.get("Name") for feat in layer]\n['Pueblo', 'Lawrence', 'Houston']\n>>> ks_extent = (-102.051, 36.99, -94.59, 40.00) # Extent for state of Kansas\n>>> layer.spatial_filter = ks_extent\n>>> len(layer)\n1\n>>> [feat.get("Name") for feat in layer]\n['Lawrence']\n>>> layer.spatial_filter = None\n>>> len(layer)\n3\nレイヤーの各フィーチャについて、指定されたフィールドの値のリストを返すメソッド:
\n>>> layer.get_fields("Name")\n['Pueblo', 'Lawrence', 'Houston']\nレイヤ内の各フィーチャのジオメトリを含むリストを返すメソッドです。オプションの引数 geos が True に設定されている場合、ジオメトリは GEOSGeometry オブジェクトに変換されます。そうでない場合は、 OGRGeometry オブジェクトとして返されます:
>>> [pt.tuple for pt in layer.get_geoms()]\n[(-104.609252, 38.255001), (-95.23506, 38.971823), (-95.363151, 29.763374)]\n指定された機能 (文字列) をサポートしているかどうかを示す真偽値を返します。 有効な機能文字列の例には、 'RandomRead', 'SequentialWrite', 'RandomWrite', 'FastSpatialFilter', 'FastFeatureCount', 'FastGetExtent', 'CreateField', 'Transactions', 'DeleteFeature', 'FastSetNextByIndex' などがあります。
FeatureLink to this heading#Feature は OGR の feture をラップします。Feature オブジェクトを直接作成することはありません。代わりに、 Layer オブジェクトからそれらを取得します。各フィーチャにはジオメトリと追加のプロパティを含むフィールドのセットがあります。フィールドのジオメトリは、その geom プロパティを介してアクセスできます。これは OGRGeometry オブジェクトを返します。Feature はフィールドの標準的な Python コンテナのように動作し、そのフィールドを Field オブジェクトとして返します。フィールドはインデックスまたは名前を指定して直接アクセスすることもでき、あるいは for ループ内でフィーチャのフィールドをイテレートすることもできます。
このフィーチャのジオメトリを OGRGeometry オブジェクトとして返します:
>>> city.geom.tuple\n(-104.609252, 38.255001)\nこのフィーチャーのフィールド (名前で指定) の値を返すメソッドで、 Field ラッパーオブジェクト ではありません :
>>> city.get("Population")\n102121\nこのフィーチャのジオメトリの種類を、 OGRGeomType オブジェクトとして返します。これは同じレイヤー内のすべてのフィーチャに対して同じであり、フィーチャの元となる Layer オブジェクトの Layer.geom_type プロパティと同等です。
フィーチャに関連付けられたデータフィールドの数を返します。これは、特定のレイヤー内のすべてのフィーチャに対して同じであり、フィーチャが属する Layer オブジェクトの Layer.num_fields プロパティと同等です。
フィーチャに関連付けられたデータのフィールド名のリストを返します。これは与えられたレイヤーの全てのフィーチャで同じであり、フィーチャの元となった Layer オブジェクトの Layer.fields プロパティと等価です。
レイヤー内のフィーチャ識別子を返します:
\n>>> city.fid\n0\nそのフィーチャーの Layer の名前を返します。これは指定されたレイヤーの全てのフィーチャに対して同じです:
>>> city.layer_name\n'cities'\n与えられたフィールド名のインデックスを返すメソッド。これは、指定されたレイヤーのすべてのフィーチャで同じになります:
\n>>> city.index("Population")\n1\nFieldLink to this heading#このフィールドの名前を返します:
\n>>> city["Name"].name\n'Name'\nこのフィールドの OGR 型を整数で返します。 FIELD_CLASSES 辞書は、これらの値を Field のサブクラスにマッピングします:
>>> city["Density"].type\n2\nこのフィールドのデータ型の名前を文字列で返します:
\n>>> city["Name"].type_name\n'String'\nこのフィールドの値を返します。Field クラス自体は値を文字列として返しますが、各サブクラスは最適な形式で値を返します:
>>> city["Population"].value\n102121\nこのフィールドの幅を返します:
\n>>> city["Name"].width\n80\nこのフィールドの数値精度を返します。非数値フィールドでは意味がありません (ゼロに設定されます) 。
\n>>> city["Density"].precision\n15\nフィールドの値を double (float) で返します:
\n>>> city["Density"].as_double()\n874.7\nフィールドの値を整数で返します:
\n>>> city["Population"].as_int()\n102121\nフィールドの値を文字列で返します:
\n>>> city["Name"].as_string()\n'Pueblo'\nフィールドの値を、日付と時刻のタプルで返します:
\n>>> city["Created"].as_datetime()\n(c_long(1999), c_long(5), c_long(23), c_long(0), c_long(0), c_long(0), c_long(0))\nDriverLink to this heading#Driver クラスは、OGR DataSource ドライバをラップするために内部で使用されます。
現在登録されているOGRベクタードライバーの数を返します。
\nOGRGeometryLink to this heading#OGRGeometry オブジェクトは GEOSGeometry オブジェクトと同様の機能を持ち、OGR の内部ジオメトリ表現の薄いラッパーです。そのため、 DataSource を使用する場合、より効率的にデータにアクセスできます。GEOSとは異なり、OGRGeometry は空間参照系 (spatial reference system) と座標変換をサポートしています。
>>> from django.contrib.gis.gdal import OGRGeometry\n>>> polygon = OGRGeometry("POLYGON((0 0, 5 0, 5 5, 0 5))")\nこのオブジェクトは 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、その他も同様です。
\n>>> polygon.dimension\n2\nこのジオメトリの座標の次元を返します。例えば、二次元のジオメトリであれば、その値は 2 になります。
\n\nこのジオメトリが Z 座標を持っているかどうかを示す真偽値。
\nZ 座標次元を追加または削除するメソッド。
\n>>> p = OGRGeometry("POINT (1 2 3)")\n>>> p.is_3d\nTrue\n>>> p.set_3d(False)\n>>> p.wkt\n"POINT (1 2)"\nこのジオメトリが M 座標を持っているかどうかを示す真偽値。
\nM 座標次元を追加または削除するメソッド。
\n>>> p = OGRGeometry("POINT (1 2)")\n>>> p.is_measured\nFalse\n>>> p.set_measured(True)\n>>> p.wkt\n"POINT M (1 2 0)"\nこのジオメトリの要素数を返します:
\n>>> polygon.geom_count\n1\nこのジオメトリがカーブジオメトリであるか、またはカーブジオメトリを含んでいるかどうかを示すブール値です。
\nジオメトリの線形バージョンを返します。変換ができない場合は、元のジオメトリがそのまま返されます。
\nジオメトリの曲線バージョンを返します。変換できない場合は、元のジオメトリがそのまま返されます。
\nこのジオメトリを記述するのに使用されたポイントの数を返します:
\n>>> polygon.point_count\n4\npoint_count のエイリアス。
point_count のエイリアス。
このジオメトリのタイプを OGRGeomType オブジェクトとして返します。
このジオメトリの型名を返します:
\n>>> polygon.geom_name\n'POLYGON'\nこのジオメトリの面積を返すか、面積を含まないジオメトリの場合は0を返します。
\n>>> polygon.area\n25.0\nこのジオメトリの外接矩形を Envelope オブジェクトとして返します。
このジオメトリの外接矩形を Envelope オブジェクトではなく、4値タプルとして返します。
>>> point.extent\n(0.0, 0.0, 5.0, 5.0)\nこのプロパティはこのジオメトリの空間参照系を制御します。空間参照系が割り当てられていない場合は None となります。割り当てられている場合、このプロパティにアクセスすると SpatialReference オブジェクトが返されます。このプロパティは別の SpatialReference オブジェクト、または SpatialReference が受け付ける任意の入力で設定できます。以下に例を示します。
>>> city.geom.srs.name\n'GCS_WGS_1984'\nこのジオメトリの SpatialReference に対応する空間参照系識別子を返すか設定します。このジオメトリに関連付けられている空間参照系情報がない場合、または SRID を決定できない場合は None を返します。
このジオメトリに対応する GEOSGeometry オブジェクトを返します。
このジオメトリの GML 形式の文字列表現を返します:
\n>>> OGRGeometry("POINT(1 2)").gml\n'<gml:Point><gml:coordinates>1,2</gml:coordinates></gml:Point>'\nこのジオメトリのHEX WKB形式での文字列表現を返します:
\n>>> OGRGeometry("POINT(1 2)").hex\n'0101000000000000000000F03F0000000000000040'\nこのジオメトリの文字列表現を JSON 形式で返します:
\n>>> OGRGeometry("POINT(1 2)").json\n'{ "type": "Point", "coordinates": [ 1.000000, 2.000000 ] }'\nこのジオメトリを KML 形式の文字列で返します。
\nこのジオメトリの WKB 表現を保持するために必要な WKB バッファのサイズを返します:
\n>>> OGRGeometry("POINT(1 2)").wkb_size\n21\nこのジオメトリのWKB表現を含む buffer を返します。
このジオメトリの文字列表現を WKT 形式で返します。
\nこのジオメトリの EWKT 表現を返します。
\nこのジオメトリオブジェクトの新しい OGRGeometry クローンを返します。
このジオメトリ内に閉じていないリングがある場合、このルーチンは始点を終点に追加することで閉じます:
\n>>> triangle = OGRGeometry("LINEARRING (0 0,0 1,1 0)")\n>>> triangle.close_rings()\n>>> triangle.wkt\n'LINEARRING (0 0,0 1,1 0,0 0)'\nこのジオメトリを異なる空間参照系に変換します。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 を返します。
点ジオメトリの座標をタプルとして、線ジオメトリの座標をタプルのタプルとして返す:
\n>>> OGRGeometry("POINT (1 2)").tuple\n(1.0, 2.0)\n>>> OGRGeometry("LINESTRING (1 2,3 4)").tuple\n((1.0, 2.0), (3.0, 4.0))\ntuple のエイリアス。
この点のX座標を返します:
\n>>> OGRGeometry("POINT (1 2)").x\n1.0\nこの点のY座標を返します:
\n>>> OGRGeometry("POINT (1 2)").y\n2.0\nこの点の Z 座標を返し、点が Z 座標を持たない場合は None を返します:
>>> OGRGeometry("POINT (1 2 3)").z\n3.0\nこの点の M 座標を返します。もし点に M 座標がない場合は None を返します。
>>> OGRGeometry("POINT ZM (1 2 3 4)").m\n4.0\nこの行のX座標のリストを返します:
\n>>> OGRGeometry("LINESTRING (1 2,3 4)").x\n[1.0, 3.0]\nこの行のY座標のリストを返します:
\n>>> OGRGeometry("LINESTRING (1 2,3 4)").y\n[2.0, 4.0]\nこのラインの Z 座標のリストを返すか、Z 座標を持たない場合は None を返します:
>>> OGRGeometry("LINESTRING (1 2 3,4 5 6)").z\n[3.0, 6.0]\nこの線分の M 座標のリストを返します。もし線分に M 座標がない場合は None を返します。
>>> OGRGeometry("LINESTRING(0 1 2 10, 1 2 3 11, 2 3 4 12)").m\n[10.0, 11.0, 12.0]\nこのポリゴンのシェルまたは外側のリングを LinearRing ジオメトリとして返します。
shell のエイリアス。
ジオメトリをこのジオメトリ コレクションに追加します。他のジオメトリタイプには適用できません。
\nOGRGeomTypeLink to this heading#このクラスは、OGR ジオメトリタイプをいくつかの方法で表現できます:
\n>>> from django.contrib.gis.gdal import OGRGeomType\n>>> gt1 = OGRGeomType(3) # Using an integer for the type\n>>> gt2 = OGRGeomType("Polygon") # Using a string\n>>> gt3 = OGRGeomType("POLYGON") # It's case-insensitive\n>>> print(gt1 == 3, gt1 == "Polygon") # Equivalence works w/non-OGRGeomType objects\nTrue True\nOGR Geometry の型の省略文字列形式を返します:
\n>>> gt1.name\n'Polygon'\nOGR ジオメトリタイプに対応する数値を返します:
\n>>> gt1.num\n3\nこの OGR 型を格納するために使用する Django フィールド型 (GeometryField のサブクラス) を返し、適切な Django 型がない場合は None を返します:
>>> gt1.django\n'PolygonField'\nEnvelopeLink to this heading#矩形のバウンディングボックスの最小、最大 X、Y 座標を含む OGR Envelope 構造体を表す。変数の名前は OGR Envelope C 構造体と互換性があります。
\nX座標の最小値。
\nX座標の最大値。
\nY座標の最小値。
\nY座標の最大値。
\n右上座標を表すタプル。
\n左下座標を表すタプル。
\nエンベロープを表すタプル。
\nWKT 形式でこのエンベロープを表すポリゴンを示す文字列。
\nSpatialReferenceLink to this heading#空間参照系 (spatial reference) オブジェクトは与えられた srs_input に基づいて初期化されます。この srs_input には以下のいずれかが指定されることがあります。
OGC Well Known Text (WKT) (文字列)
EPSG コード (整数または文字列)
PROJ 文字列
よく知られている標準の略記文字列 ('WGS84', 'WGS72', 'NAD27', 'NAD83')
例:
\n>>> wgs84 = SpatialReference("WGS84") # shorthand string\n>>> wgs84 = SpatialReference(4326) # EPSG code\n>>> wgs84 = SpatialReference("EPSG:4326") # EPSG string\n>>> proj = "+proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs "\n>>> wgs84 = SpatialReference(proj) # PROJ string\n>>> wgs84 = SpatialReference("""GEOGCS["WGS 84",\n... DATUM["WGS_1984",\n... SPHEROID["WGS 84",6378137,298.257223563,\n... AUTHORITY["EPSG","7030"]],\n... AUTHORITY["EPSG","6326"]],\n... PRIMEM["Greenwich",0,\n... AUTHORITY["EPSG","8901"]],\n... UNIT["degree",0.01745329251994328,\n... AUTHORITY["EPSG","9122"]],\n... AUTHORITY["EPSG","4326"]]""") # OGC WKT\n与えられた文字列属性ノードの値を返します。ノードが存在しない場合は None を返します。(target, child) のようなタプルもパラメータとして受け取ることができます。ここで、child は WKT 内の属性のインデックスです。例えば:
>>> wkt = 'GEOGCS["WGS 84", DATUM["WGS_1984, ... AUTHORITY["EPSG","4326"]]'\n>>> srs = SpatialReference(wkt) # could also use 'WGS84', or 4326\n>>> print(srs["GEOGCS"])\nWGS 84\n>>> print(srs["DATUM"])\nWGS_1984\n>>> print(srs["AUTHORITY"])\nEPSG\n>>> print(srs["AUTHORITY", 1]) # The authority value\n4326\n>>> print(srs["TOWGS84", 4]) # the fourth value in this wkt\n0\n>>> print(srs["UNIT|AUTHORITY"]) # For the units authority, have to use the pipe symbol.\nEPSG\n>>> print(srs["UNIT|AUTHORITY", 1]) # The authority value for the units\n9122\n指定されたターゲットノードの属性値 (例: 'PROJCS')。index キーワードは返す子ノードのインデックスを指定します。
指定された文字列の対象ノードに対するオーソリティ名を返します。
\n与えられた文字列の対象ノードのオーソリティコードを返します。
\nこの空間参照系 (SpatialReference) オブジェクトのクローンを返します。
\nこのメソッドは SpatialReference の WKT を検査し、EPSG 識別子が該当する場合は EPSG オーソリティノードを追加します。
この SpatialReference を ESRI のフォーマットから EPSG に変換します。
\nこの SpatialReference を ESRI のフォーマットに変換します。
\n与えられた空間参照系が有効かどうかをチェックし、有効でない場合は例外が発生します。
\nEPSGコードから空間参照系をインポートします。
\nPROJ文字列から空間参照系をインポートします。
\nWKT から空間参照系をインポートします。
\nXML から空間参照系をインポートします。
\nこの空間参照系の名前を返します。
\nトップレベルのオーソリティの SRID を返すか、未定義の場合は None を返します。
線形単位の名前を返します。
\n線形単位の値を返します。
\n角度の単位の名前を返します。
\n角度の単位の値を返します。
\n線形単位を返すか角度単位を返すかを自動的に決定し、単位値と単位名の2タプルを返します。
\nこの空間参照系の楕円体パラメータのタプル (半長軸、半短軸、逆扁平率) を返します。
\nこの空間参照系に対する楕円体の長軸を返します。
\nこの空間参照系における楕円体の短半径を返します。
\nこの空間参照系の楕円体の逆扁平率を返します。
\nこの空間参照系が地理座標系 (親ノードが GEOGCS) である場合に True を返します。
この空間参照系がローカルである場合 (ルートノードが LOCAL_CS である場合) には True を返します。
この空間参照系が投影座標系である場合 (ルートノードが PROJCS である場合) に True を返します。
この空間参照系の WKT 表現を返します。
\nWKT のきれい ("pretty") な表現を返します。
\nこの空間参照系の PROJ 表現を返します。
\nSpatialReference.proj のエイリアス。
この空間参照系の XML 表現を返します。
\nCoordTransformLink to this heading#座標系の変換を表します。ソースとターゲットの座標系をそれぞれ表す2つの SpatialReference で初期化されます。これらのオブジェクトは、異なるジオメトリに繰り返し同じ座標変換を行う場合に使用するべきです。
>>> ct = CoordTransform(SpatialReference("WGS84"), SpatialReference("NAD83"))\n>>> for feat in layer:\n... geom = feat.geom # getting clone of feature geometry\n... geom.transform(ct) # transforming\n...\nGDALRasterLink to this heading#GDALRaster は、GDALがサポートする様々な地理空間ファイル形式やデータソースからデータを読み込むための一貫したインターフェースを使用するGDALラスターソースオブジェクトのラッパーです。各データソースは、1つ以上のデータ層である「bands（バンド）」という名前のデータを含む GDALRaster オブジェクトによって表されます。各バンドは、ジオリファレンスされた画像データを含む GDALBand オブジェクトによって表されます。例えば、RGB画像は、赤、緑、青の各色に対応する3つのバンドで表されます。
GDALRaster のコンストラクターは 2 つのパラメーターを受け取ります。最初のパラメーターはラスターのソースを定義し、2 番目のパラメーターはラスターを書き込みモードで開くかどうかを定義します。新しく作成されたラスターに対しては、2 番目のパラメーターは無視され、常に書き込みモードで新しいラスターが作成されます。
最初のパラメータは、３つの形式を取ることができます。ファイルパスを表す文字列あるいは Path オブジェクト (ファイルシステムまたはGDAL仮想ファイルシステム) 、新しいラスタを定義する値が入った辞書、ラスタファイルを表すバイトオブジェクトです。
入力がファイルパスの場合は、そのパスからラスタがオープンされます。入力が辞書の生データの場合、パラメータ width、height、srid が必要です。入力がバイトオブジェクトの場合、GDALの仮想ファイルシステムを使用してオープンされます。
辞書入力を使用してラスタを作成する方法の詳細については、 データからラスターを作成する を参照してください。仮想ファイルシステム内でラスタを作成する方法の詳細については、 GDALの仮想ファイルシステムを使う を参照してください。
\n以下の例では、異なる入力ソースからどのようにラスタを作成できるかを示します (GeoDjango のテストのサンプルデータを使います。 サンプルデータ セクションも参照してください) 。
\n>>> from django.contrib.gis.gdal import GDALRaster\n>>> rst = GDALRaster("/path/to/your/raster.tif", write=False)\n>>> rst.name\n'/path/to/your/raster.tif'\n>>> rst.width, rst.height # This file has 163 x 174 pixels\n(163, 174)\n>>> rst = GDALRaster(\n... { # Creates an in-memory raster\n... "srid": 4326,\n... "width": 4,\n... "height": 4,\n... "datatype": 1,\n... "bands": [\n... {\n... "data": (2, 3),\n... "offset": (1, 1),\n... "size": (2, 2),\n... "shape": (2, 1),\n... "nodata_value": 5,\n... }\n... ],\n... }\n... )\n>>> rst.srs.srid\n4326\n>>> rst.width, rst.height\n(4, 4)\n>>> rst.bands[0].data()\narray([[5, 5, 5, 5],\n [5, 2, 3, 5],\n [5, 2, 3, 5],\n [5, 5, 5, 5]], dtype=uint8)\n>>> rst_file = open("/path/to/your/raster.tif", "rb")\n>>> rst_bytes = rst_file.read()\n>>> rst = GDALRaster(rst_bytes)\n>>> rst.is_vsi_based\nTrue\n>>> rst.name # Stored in a random path in the vsimem filesystem.\n'/vsimem/da300bdb-129d-49a8-b336-e410a9428dad'\n入力ファイルのパスまたはインスタンス化の際に指定された名前に相当するソースの名前。
\n>>> GDALRaster({"width": 10, "height": 10, "name": "myraster", "srid": 4326}).name\n'myraster'\n入力ファイルを処理するために使用されるGDALドライバの名前です。ファイルから作成された GDALRaster では、ドライバタイプが自動的に検出されます。ゼロからラスターを作成する場合、デフォルトではメモリ内のラスター ('MEM') ですが、必要に応じて変更できます。例えば、 GeoTiff ファイルには GTiff を使用します。ファイルタイプの一覧については、 GDAL Raster Formats リストも参照してください。
インメモリのラスタは次の例のように作成できます:
\n>>> GDALRaster({"width": 10, "height": 10, "srid": 4326}).driver.name\n'MEM'\nファイルベースの GeoTiff ラスターは、次の例のように作成できます:
\n>>> import tempfile\n>>> rstfile = tempfile.NamedTemporaryFile(suffix=".tif")\n>>> rst = GDALRaster(\n... {\n... "driver": "GTiff",\n... "name": rstfile.name,\n... "srid": 4326,\n... "width": 255,\n... "height": 255,\n... "nr_of_bands": 1,\n... }\n... )\n>>> rst.name\n'/tmp/tmp7x9H4J.tif' # The exact filename will be different on your computer\n>>> rst.driver.name\n'GTiff'\nピクセル単位のソースの幅(X軸)。
\n>>> GDALRaster({"width": 10, "height": 20, "srid": 4326}).width\n10\nピクセル単位のソースの高さ(Y軸)。
\n>>> GDALRaster({"width": 10, "height": 20, "srid": 4326}).height\n20\nラスターの空間参照系を SpatialReference インスタンスとして返します。SRSは他の SpatialReference に設定したり、 SpatialReference コンストラクタが受け付ける任意の入力を与えることで変更できます。
>>> rst = GDALRaster({"width": 10, "height": 20, "srid": 4326})\n>>> rst.srs.srid\n4326\n>>> rst.srs = 3086\n>>> rst.srs.srid\n3086\nラスタの空間参照システム識別子(SRID)。このプロパティは srs 属性を通してSRIDを取得または設定するためのショートカットです。
>>> rst = GDALRaster({"width": 10, "height": 20, "srid": 4326})\n>>> rst.srid\n4326\n>>> rst.srid = 3086\n>>> rst.srid\n3086\n>>> rst.srs.srid # This is equivalent\n3086\n以下のリレーションシップを使用して、ピクセル/ライン座標をジオリファレンス空間にマッピングする6つの係数のタプルとして表現されるアフィン変換行列:
\nXgeo = GT(0) + Xpixel * GT(1) + Yline * GT(2)\nYgeo = GT(3) + Xpixel * GT(4) + Yline * GT(5)\n同じ値は、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})\n>>> rst.geotransform\n[0.0, 1.0, 0.0, 0.0, 0.0, -1.0]\nソースの空間参照システム内で、ラスターの左上原点の座標を、x と y のメンバを持つ点オブジェクトとして表します。
>>> rst = GDALRaster({"width": 10, "height": 20, "srid": 4326})\n>>> rst.origin\n[0.0, 0.0]\n>>> rst.origin.x = 1\n>>> rst.origin\n[1.0, 0.0]\nラスタのジオリファレンスに使用されるピクセルの幅と高さ。 x と y のメンバーを持つポイントオブジェクトとして表されます。詳細については、geotransform を参照してください。
>>> rst = GDALRaster({"width": 10, "height": 20, "srid": 4326})\n>>> rst.scale\n[1.0, -1.0]\n>>> rst.scale.x = 2\n>>> rst.scale\n[2.0, -1.0]\nラスターを x と y のメンバを持つ点オブジェクトとしてジオリファレンスするために使用されるスキュー係数。ノースアップ画像の場合、これらの係数は両方とも 0 です。
>>> rst = GDALRaster({"width": 10, "height": 20, "srid": 4326})\n>>> rst.skew\n[0.0, 0.0]\n>>> rst.skew.x = 3\n>>> rst.skew\n[3.0, 0.0]\nラスターソースの範囲 (境界値)。ソースの空間参照システムにおける4タプル (xmin, ymin, xmax, ymax) で表されます。
>>> rst = GDALRaster({"width": 10, "height": 20, "srid": 4326})\n>>> rst.extent\n(0.0, -20.0, 10.0, 0.0)\n>>> rst.origin.x = 100\n>>> rst.extent\n(100.0, -20.0, 110.0, 0.0)\nソースの全てのバンドを、GDALBand のインスタンスとしてリストアップします。
>>> rst = GDALRaster(\n... {\n... "width": 1,\n... "height": 2,\n... "srid": 4326,\n... "bands": [{"data": [0, 1]}, {"data": [2, 3]}],\n... }\n... )\n>>> len(rst.bands)\n2\n>>> rst.bands[1].data()\narray([[ 2., 3.]], dtype=float32)\nこのラスターの変形されたバージョンを返します。
\n変化系パラメータは ds_input 引数で指定することができる。 ds_input の使用方法はクラスコンストラクタの対応する引数と同様です。 ds_input はターゲットとなるラスターの特性を辞書に格納したものです。指定できる辞書のキーは width, height, SRID, origin, scale, skew, datatype, driver, name (filename) です。
デフォルトでは、warp 関数はほとんどのパラメータを元のソースラスタの値と等しく保ちます。したがって、変更する必要のあるパラメータのみを指定すればよいです。ただし、これにはドライバも含まれるので、ファイルベースのラスタの場合、warp 関数は新しいラスタをディスク上に作成します。
\nソースラスタと異なるパラメータは名前のみです。ラスタ名のデフォルト値は、ソースラスタの名前に '_copy' + source_driver_name を追加したものです。ファイルベースのラスタの場合は、ターゲットとなるラスタのファイルパスを指定することを推奨します。
変形に使用するリサンプリングアルゴリズムは resampling 引数で指定します。デフォルトは NearestNeighbor で、他に指定できる値は Bilinear、Cubic、CubicSpline、Lanczos、Average、Mode です。
max_error 引数は、変換の近似において許容される入力ピクセルでの最大誤差を指定するために使用できます。デフォルト値は正確な計算のために 0.0 です。
GDAL に詳しいユーザーにとって、この関数は gdalwarp コマンドラインユーティリティと同様の機能を持っています。
たとえば、warp 関数は、ラスターを元のピクセルスケールの2倍に集約するために使用できます:
>>> rst = GDALRaster(\n... {\n... "width": 6,\n... "height": 6,\n... "srid": 3086,\n... "origin": [500000, 400000],\n... "scale": [100, -100],\n... "bands": [{"data": range(36), "nodata_value": 99}],\n... }\n... )\n>>> target = rst.warp({"scale": [200, -200], "width": 3, "height": 3})\n>>> target.bands[0].data()\narray([[ 7., 9., 11.],\n [ 19., 21., 23.],\n [ 31., 33., 35.]], dtype=float32)\nこのラスタを異なる空間参照システム (srs) に変換します。これは SpatialReference オブジェクトであり、または SpatialReference で受け入れられる他の入力 (空間参照 WKT や PROJ 文字列、整数の SRID を含む) でも構いません。
新しい空間参照系で現在のラスタの境界とスケールを計算し、 warp 関数を使用してラスタを変形します。
デフォルトでは、コピー元のラスタのドライバが使用され、ラスタ名はコピー元の名前に '_copy' + source_driver_name を追加したものになります。異なるドライバや名前を指定するには、driver および name 引数を使用します。
デフォルトのリサンプリングアルゴリズムは NearestNeighbour ですが、 resampling 引数を使用して変更できます。デフォルトのリサンプリングに許容される最大誤差は0.0であり、 max_error 引数を使用して変更できます。これらの引数の詳細については、 warp のドキュメントを参照してください。
>>> rst = GDALRaster(\n... {\n... "width": 6,\n... "height": 6,\n... "srid": 3086,\n... "origin": [500000, 400000],\n... "scale": [100, -100],\n... "bands": [{"data": range(36), "nodata_value": 99}],\n... }\n... )\n>>> target_srs = SpatialReference(4326)\n>>> target = rst.transform(target_srs)\n>>> target.origin\n[-82.98492744885776, 27.601924753080144]\nラスターの概要を含む文字列を返します。これは gdalinfo コマンドラインユーティリティに相当します。
\nこのラスタのメタデータ。ネストした辞書として表現されます。第1レベルのキーはメタデータ・ドメインです。第2レベルには、各ドメインのメタデータ項目名と値が含まれます。
\nメタデータ項目を設定または更新するには、上記の入れ子構造を使用して、対応するメタデータ項目をメソッドに渡します。指定された辞書にあるキーのみが更新され、メタデータの残りの部分は変更されません。
\nメタデータ項目を削除するには、メタデータの値として None を使用します。
>>> rst = GDALRaster({"width": 10, "height": 20, "srid": 4326})\n>>> rst.metadata\n{}\n>>> rst.metadata = {"DEFAULT": {"OWNER": "Django", "VERSION": "1.0"}}\n>>> rst.metadata\n{'DEFAULT': {'OWNER': 'Django', 'VERSION': '1.0'}}\n>>> rst.metadata = {"DEFAULT": {"OWNER": None, "VERSION": "2.0"}}\n>>> rst.metadata\n{'DEFAULT': {'VERSION': '2.0'}}\nA bytes representation of this raster. Returns None for rasters\nthat are not stored in GDAL's in-memory virtual filesystem.
このラスターがGDALの仮想ファイルシステムに保存されているかどうかを示す真偽値。
\nGDALBandLink to this heading#GDALBand のインスタンスは明示的に作成されるのではなく、代わりに GDALRaster オブジェクトから、その bands 属性を介して取得されます。 GDALBand にはラスターの実際のピクセル値が含まれています。
バンドの名前または説明 (あれば) 。
\nピクセル単位でのバンドの幅 (X軸)。
\nピクセル単位のバンドの高さ (Y軸)。
\nこのバンドに含まれるピクセルの総数。これは width * height に等しいです。
このバンドのピクセル値の統計量を計算します。戻り値は (minimum, maximum, mean, standard deviation) という構造を持つタプルです。
もし approximate 引数が True に設定されている場合、統計はオーバービューまたは画像タイルのサブセットに基づいて計算されます。
refresh 引数を True に設定すると、統計情報はデータから直接計算され、その結果でキャッシュが更新されます。
永続的なキャッシュ値が見つかった場合、その値が返されます。永続的補助メタデータ (PAM) サービスを使用するラスターフォーマットの場合、統計情報は補助ファイルにキャッシュされる可能性があります。場合によっては、このメタデータがピクセル値と同期していなかったり、 approximate 引数の値を反映しない以前の呼び出しの値が返されたりすることがあります。そのような場合は、 refresh 引数を使用して更新された値を取得し、キャッシュに保存します。
空のバンド (すべてのピクセル値が "no data "である) の場合、すべての統計は None として返されます。
バンドの最小ピクセル値 ("no data" 値を除く) 。
\nバンドの最大ピクセル値 ("no data" 値を除く) 。
\nバンドのすべてのピクセル値の平均 ("no data"値を除く) 。
\nバンドのすべてのピクセル値の標準偏差 ("no data" の値を除く)。
\nバンドの "no data" 値は通常、有効なデータではないピクセルを示す特別なマーカー値です。これらのピクセルは一般的に表示されるべきではなく、分析操作にも寄与すべきではありません。
\n既存の "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の利用を強くお勧めします。
\ndata パラメータが指定された場合、データは GDALBand に書き込まれます。入力には、パックされた文字列、バッファ、リスト、配列、NumPyの配列のいずれかを指定できます。入力に含まれるアイテムの数は、通常、バンド内の全ピクセル数に対応するか、 offset と size パラメータが指定されている場合は、ピクセル値の特定のブロックのピクセル数に対応します。
入力のアイテム数がターゲットピクセルブロックと異なる場合は、 shape パラメータを指定する必要があります。shapeは入力データの幅と高さをピクセル単位で指定するタプルです。その後、選択されたブロックのピクセル値を更新するためにデータが複製されます。これは、たとえば1つの値でバンド全体を塗りつぶすのに便利です。
例:
\n>>> rst = GDALRaster(\n... {"width": 4, "height": 4, "srid": 4326, "datatype": 1, "nr_of_bands": 1}\n... )\n>>> bnd = rst.bands[0]\n>>> bnd.data(range(16))\n>>> bnd.data()\narray([[ 0, 1, 2, 3],\n [ 4, 5, 6, 7],\n [ 8, 9, 10, 11],\n [12, 13, 14, 15]], dtype=int8)\n>>> bnd.data(offset=(1, 1), size=(2, 2))\narray([[ 5, 6],\n [ 9, 10]], dtype=int8)\n>>> bnd.data(data=[-1, -2, -3, -4], offset=(1, 1), size=(2, 2))\n>>> bnd.data()\narray([[ 0, 1, 2, 3],\n [ 4, -1, -2, 7],\n [ 8, -3, -4, 11],\n [12, 13, 14, 15]], dtype=int8)\n>>> bnd.data(data="\\x9d\\xa8\\xb3\\xbe", offset=(1, 1), size=(2, 2))\n>>> bnd.data()\narray([[ 0, 1, 2, 3],\n [ 4, -99, -88, 7],\n [ 8, -77, -66, 11],\n [ 12, 13, 14, 15]], dtype=int8)\n>>> bnd.data([1], shape=(1, 1))\n>>> bnd.data()\narray([[1, 1, 1, 1],\n [1, 1, 1, 1],\n [1, 1, 1, 1],\n [1, 1, 1, 1]], dtype=uint8)\n>>> bnd.data(range(4), shape=(1, 4))\narray([[0, 0, 0, 0],\n [1, 1, 1, 1],\n [2, 2, 2, 2],\n [3, 3, 3, 3]], dtype=uint8)\nこのバンドのメタデータ。機能は GDALRaster.metadata と同じです。
このセクションでは、ds_input パラメータを使用してゼロからラスタを作成する方法について説明します。
GDALRaster コンストラクタに dict を渡すことで新しいラスタが作成されます。この辞書には、新しいラスタの定義パラメータ (原点、サイズ、空間参照システムなど) が含まれています。さらに、辞書にはピクセルデータや新しいラスタのフォーマットに関する情報も含めることができます。したがって、作成されるラスタは指定したドライバに応じてファイルベースまたはメモリベースとなります。
ラスターデータを辞書や JSON で記述する標準はありません。そのため、 GDALRaster クラスへの辞書入力の定義は Django 固有のものです。これは geojson フォーマットに触発されたものですが、 geojson 標準は今のところベクタフォーマットに限定されています。
GDALRaster と GDALBand クラスの対応する属性とメソッドのドキュメントに、ラスタを作成する際に異なるキーを使用する例があります。
ds_input 辞書Link to this heading#ds_input 辞書内でラスターを作成するために必要なキーはわずかで、 width、height、および srid です。その他のパラメータにはデフォルト値が設定されています(以下の表を参照)。 ds_input 辞書に渡せるキーのリストは、GDALRaster プロパティに密接に関連していますが、完全には一致していません。多くのパラメータはこれらのプロパティに直接マッピングされており、その他のパラメータについては以下で説明します。
以下の表は、ds_input 辞書に設定できるすべてのキーを説明しています。
キー | \nデフォルト値 | \n使い方 | \n
|---|---|---|
| \n必須 | \n
| \n
| \n必須 | \n
| \n
| \n必須 | \n
| \n
| \n
| \n
| \n
| \n
| \n下記参照 | \n
| \n
| \n
| \n
| \n
| \n
| \n
| \n
| \n
| \n
| \n
| \n下記参照 | \n
| \n
| \n下記参照 | \n
| \n
| \n下記参照 | \n
| \n
| \n下記参照 | \n
ラスタの名前を表す文字列。ファイルベースのラスタを作成する場合、このパラメータは新しいラスタのファイルパスでなければなりません。名前が /vsimem/ で始まる場合、ラスタはGDALの仮想ファイルシステムに作成されます。
全バンドのデータ型を表す整数です。デフォルトは 6 (Float32) です。新しいラスターの全てのバンドは同じデータ型である必要があります。値のマッピングは以下の通りです:
値 | \nGDAL ピクセルタイプ | \n説明 | \n
|---|---|---|
1 | \nGDT_Byte | \n8 ビット符号なし整数 | \n
2 | \nGDT_UInt16 | \n16 ビット符号なし整数 | \n
3 | \nGDT_Int16 | \n16ビット符号付き整数 | \n
4 | \nGDT_UInt32 | \n32ビット符号なし整数 | \n
5 | \nGDT_Int32 | \n32ビット符号付き整数 | \n
6 | \nGDT_Float32 | \n32ビット浮動小数点 | \n
7 | \nGDT_Float64 | \n64ビット浮動小数点 | \n
12 | \nGDT_UInt64 | \n64ビット符号なし整数 (GDAL 3.5+) | \n
13 | \nGDT_Int64 | \n64ビット符号付き整数 (GDAL 3.5+) | \n
14 | \nGDT_Int8 | \n8ビット符号付き整数 (GDAL 3.7+) | \n
ラスターのバンド数を表す整数。作成時にバンドデータを渡さずにラスタを作成することもできます。バンド数が指定されていない場合は、入力 bands の長さから自動的に計算されます。作成後にバンド数を変更することはできません。
バンド入力データを持つ band_input 辞書のリスト。結果のバンドインデックスは、提供されたリストと同じです。バンド入力辞書の定義は以下に示します。バンドデータが提供されない場合、ラスターバンドの値はゼロの配列としてインスタンス化され、"no data" の値は None に設定されます。
ラスターを作成するためのオプションを持つ辞書。入力辞書のキーと値のペアは、ラスターの作成時にドライバーに渡されます。
\n利用可能なオプションは、ドライバごとのドキュメントに記載されています。
\n辞書内の値は大文字と小文字を区別しないため、作成時に適切な文字列形式に自動変換されます。
\n以下の例では、GTiff driver で利用可能なオプションのいくつかを使用しています。結果は、内部タイルスキームを持つ圧縮されたラスターです。内部タイルのブロックサイズは23×23です:
\n>>> GDALRaster(\n... {\n... "driver": "GTiff",\n... "name": "/path/to/new/file.tif",\n... "srid": 4326,\n... "width": 255,\n... "height": 255,\n... "nr_of_bands": 1,\n... "papsz_options": {\n... "compress": "packbits",\n... "tiled": "yes",\n... "blockxsize": 23,\n... "blockysize": 23,\n... },\n... }\n... )\nband_input 辞書Link to this heading#ds_input 辞書内の bands キーは band_input 辞書のリストです。各 band_input 辞書には、新しいラスタのバンドに設定されるピクセル値と "no data" 値が含まれている可能性があります。データ配列は新しいラスタの全サイズを持つことも、より小さくなることもあります。新しいラスタより小さい配列の場合、size、shape、および offset キーがピクセル値を制御します。これらの対応するキーは、data() メソッドに渡されます。これらの機能は、そのメソッドでバンドデータを設定するのと同じです。以下の表は使用可能なキーを説明しています。
キー | \nデフォルト値 | \n使い方 | \n
|---|---|---|
| \n
| \n
| \n
| \n
| \n
| \n
| \nラスターの | \n
| \n
| \nsize と同じ | \n
| \n
| \n
| \n
| \n
GDALはファイルシステムに保存されたファイルにアクセスするだけでなく、圧縮、暗号化、またはリモートファイルなど他の種類のファイルにアクセスを抽象化するための仮想ファイルシステムもサポートしています。
\nGDALには、メモリをファイルとして扱う内部メモリベースのファイルシステムがあります。これを使用して、 GDALRaster オブジェクトをバイナリファイルバッファーに読み込んだり書き込んだりすることができます。
これは、ラスターがリモートストレージからバッファとして取得されたり、ディスクに書き込まれることなくビューから返されたりするようなWebコンテキストで便利です。
\nbytes オブジェクトが入力として提供されるか、ファイルパスが /vsimem/ で始まる場合、仮想ファイルシステム内に GDALRaster オブジェクトが作成されます。
入力は、ファイルの bytes としての完全なバイナリ表現である必要があります。例えば、:
# Read a raster as a file object from a remote source.\n>>> from urllib.request import urlopen\n>>> dat = urlopen("https://example.com/raster.tif").read()\n# Instantiate a raster from the bytes object.\n>>> rst = GDALRaster(dat)\n# The name starts with /vsimem/, indicating that the raster lives in the\n# virtual filesystem.\n>>> rst.name\n'/vsimem/da300bdb-129d-49a8-b336-e410a9428dad'\n新しい仮想ファイルベースのラスタをゼロから作成するには、 ds_input 辞書表現を使用し、 /vsimem/ で始まる name 引数を指定します (辞書表現の詳細については、 データからラスターを作成する を参照してください) 。仮想ファイルベースのラスタの場合、 vsi_buffer 属性はラスタの bytes 表現を返します。
ここでは、ラスタを作成して HttpResponse にファイルとして返す方法を紹介します:
>>> from django.http import HttpResponse\n>>> rst = GDALRaster(\n... {\n... "name": "/vsimem/temporarymemfile",\n... "driver": "tif",\n... "width": 6,\n... "height": 6,\n... "srid": 3086,\n... "origin": [500000, 400000],\n... "scale": [100, -100],\n... "bands": [{"data": range(36), "nodata_value": 99}],\n... }\n... )\n>>> HttpResponse(rast.vsi_buffer, "image/tiff")\nGDAL のローカルビルドによっては、他の仮想ファイルシステムもサポートされているかもしれません。提供されたパスの先頭に適切な /vsi*/ 接頭辞を付けることで使用できます。詳細は GDAL Virtual Filesystems documentation を参照してください。
ファイルを解凍して結果のラスタをインスタンス化する代わりに、GDAL は /vsizip/, /vsigzip/, /vsitar/ 仮想ファイルシステムを使用して圧縮ファイルに直接アクセスできます:
>>> from django.contrib.gis.gdal import GDALRaster\n>>> rst = GDALRaster("/vsizip/path/to/your/file.zip/path/to/raster.tif")\n>>> rst = GDALRaster("/vsigzip/path/to/your/file.gz")\n>>> rst = GDALRaster("/vsitar/path/to/your/file.tar/path/to/raster.tif")\nGDALは、それらの機能がビルトインされている限り、オンラインリソースとストレージプロバイダを透過的にサポートできます。
\n認証なしで公開されたラスターファイルにアクセスするには、/vsicurl/ を使用できます。
>>> from django.contrib.gis.gdal import GDALRaster\n>>> rst = GDALRaster("/vsicurl/https://example.com/raster.tif")\n>>> rst.name\n'/vsicurl/https://example.com/raster.tif'\n商用のストレージプロバイダ (/vsis3/ など) の場合、認証やその他の設定を事前に行っておく必要があります (利用可能なオプションについては GDAL Virtual Filesystems documentation を参照してください) 。
Since GDALRaster always opens new rasters in write mode, it is\nessential to prevent instantiating one from untrusted input. Otherwise, an\nattacker might gain the ability to write a file or make a network request.
To mitigate this, spatial lookups prevent\nstr, pathlib.Path, and dict values from reaching\nGDALRaster altogether. To use these types with lookups, wrap them\nexplicitly with GDALRaster, indicating that the value is trusted.\nBytes are accepted without being wrapped in GDALRaster because they\nare opened through GDAL's memory-based virtual filesystem.
This protection applies only to spatial lookups. Assigning a dict value to\na RasterField will still open a new\nraster, and assigning a str or Path will still fetch and open the\nreferenced raster.
When validating geometry inputs, the\nGeometryField form field will reject raster\nvalues. When validating raster inputs, you should write custom validation.
For defense-in-depth strategies for limiting the available raster drivers, see\nGDAL security considerations.
\n\nGDAL_LIBRARY_PATHLink to this heading#GDALライブラリの場所を指定する文字列。通常、この設定はGDALライブラリが標準以外の場所にある場合にのみ使用されます (例えば /home/john/lib/libgdal.so など) 。
GDAL に関連するエラーを示す基本的な GDAL 例外。
\n空間参照システムオブジェクトの構築時または使用時にエラーが発生した場合に発生する例外。
\n