GeoDjango 教程¶
简介¶
GeoDjango 是 Django 的一个内置 contrib 模块，它将 Django 变成一个世界级的地理 Web 框架。GeoDjango 力求使创建地理 Web 应用程序（如基于位置的服务）尽可能简单。其功能包括
用于OGC 几何和栅格数据的 Django 模型字段。
Django ORM 的扩展，用于查询和操作空间数据。
松散耦合、高级 Python 接口，用于 GIS 几何和栅格操作以及不同格式的数据处理。
从管理界面编辑几何字段。
本教程假设您熟悉 Django；因此，如果您是 Django 新手，请先阅读常规教程，以先熟悉 Django。
注意
GeoDjango 除了 Django 的要求之外还有其他要求——请查阅安装文档了解更多详情。
本教程将指导您完成创建用于查看世界边界的地理 Web 应用程序。 [1] 本教程中使用的一些代码取自并/或受到GeoDjango 基本应用程序项目的启发。 [2]
注意
按顺序完成教程各个部分，以便逐步获得说明。
设置¶
创建空间数据库¶
通常不需要特殊设置，因此您可以像创建任何其他项目一样创建数据库。我们为选定的数据库提供了一些提示
创建一个新项目¶
使用标准的django-admin脚本创建一个名为geodjango的项目
$ django-admin startproject geodjango
...\> django-admin startproject geodjango
这将初始化一个新项目。现在，在geodjango项目中创建一个world Django 应用程序
$ cd geodjango
$ python manage.py startapp world
...\> cd geodjango
...\> py manage.py startapp world
配置settings.py¶
geodjango项目的设置存储在geodjango/settings.py文件中。编辑数据库连接设置以匹配您的设置
DATABASES = {
"default": {
"ENGINE": "django.contrib.gis.db.backends.postgis",
"NAME": "geodjango",
"USER": "geo",
},
}
此外，修改INSTALLED_APPS设置以包含django.contrib.admin、django.contrib.gis和world（您新创建的应用程序）
INSTALLED_APPS = [
"django.contrib.admin",
"django.contrib.auth",
"django.contrib.contenttypes",
"django.contrib.sessions",
"django.contrib.messages",
"django.contrib.staticfiles",
"django.contrib.gis",
"world",
]
地理数据¶
世界边界¶
世界边界数据可在此 zip 文件中获得。在world应用程序中创建一个data目录，下载世界边界数据并解压缩。在 GNU/Linux 平台上，使用以下命令
$ mkdir world/data
$ cd world/data
$ wget https://web.archive.org/web/20231220150759/https://thematicmapping.org/downloads/TM_WORLD_BORDERS-0.3.zip
$ unzip TM_WORLD_BORDERS-0.3.zip
$ cd ../..
...\> mkdir world\data
...\> cd world\data
...\> wget https://web.archive.org/web/20231220150759/https://thematicmapping.org/downloads/TM_WORLD_BORDERS-0.3.zip
...\> unzip TM_WORLD_BORDERS-0.3.zip
...\> cd ..\..
世界边界 ZIP 文件包含一组数据文件，统称为ESRI Shapefile，这是最流行的地理空间数据格式之一。解压缩后，世界边界数据集包含以下扩展名的文件：
.shp：保存世界边界几何的矢量数据。.shx：.shp中存储的几何的空间索引文件。.dbf：用于保存非几何属性数据（例如，整数和字符字段）的数据库文件。.prj：包含 shapefile 中存储的地理数据的空间参考信息。
使用ogrinfo检查空间数据¶
GDAL ogrinfo实用程序允许检查 shapefile 或其他矢量数据源的元数据
$ ogrinfo world/data/TM_WORLD_BORDERS-0.3.shp
INFO: Open of `world/data/TM_WORLD_BORDERS-0.3.shp'
using driver `ESRI Shapefile' successful.
1: TM_WORLD_BORDERS-0.3 (Polygon)
...\> ogrinfo world\data\TM_WORLD_BORDERS-0.3.shp
INFO: Open of `world/data/TM_WORLD_BORDERS-0.3.shp'
using driver `ESRI Shapefile' successful.
1: TM_WORLD_BORDERS-0.3 (Polygon)
ogrinfo告诉我们 shapefile 有一个图层，并且此图层包含多边形数据。要了解更多信息，我们将指定图层名称并使用-so选项仅获取重要的摘要信息
$ ogrinfo -so world/data/TM_WORLD_BORDERS-0.3.shp TM_WORLD_BORDERS-0.3
INFO: Open of `world/data/TM_WORLD_BORDERS-0.3.shp'
using driver `ESRI Shapefile' successful.
Layer name: TM_WORLD_BORDERS-0.3
Geometry: Polygon
Feature Count: 246
Extent: (-180.000000, -90.000000) - (180.000000, 83.623596)
Layer SRS WKT:
GEOGCS["GCS_WGS_1984",
DATUM["WGS_1984",
SPHEROID["WGS_1984",6378137.0,298.257223563]],
PRIMEM["Greenwich",0.0],
UNIT["Degree",0.0174532925199433]]
FIPS: String (2.0)
ISO2: String (2.0)
ISO3: String (3.0)
UN: Integer (3.0)
NAME: String (50.0)
AREA: Integer (7.0)
POP2005: Integer (10.0)
REGION: Integer (3.0)
SUBREGION: Integer (3.0)
LON: Real (8.3)
LAT: Real (7.3)
...\> ogrinfo -so world\data\TM_WORLD_BORDERS-0.3.shp TM_WORLD_BORDERS-0.3
INFO: Open of `world/data/TM_WORLD_BORDERS-0.3.shp'
using driver `ESRI Shapefile' successful.
Layer name: TM_WORLD_BORDERS-0.3
Geometry: Polygon
Feature Count: 246
Extent: (-180.000000, -90.000000) - (180.000000, 83.623596)
Layer SRS WKT:
GEOGCS["GCS_WGS_1984",
DATUM["WGS_1984",
SPHEROID["WGS_1984",6378137.0,298.257223563]],
PRIMEM["Greenwich",0.0],
UNIT["Degree",0.0174532925199433]]
FIPS: String (2.0)
ISO2: String (2.0)
ISO3: String (3.0)
UN: Integer (3.0)
NAME: String (50.0)
AREA: Integer (7.0)
POP2005: Integer (10.0)
REGION: Integer (3.0)
SUBREGION: Integer (3.0)
LON: Real (8.3)
LAT: Real (7.3)
此详细的摘要信息告诉我们图层中的要素数量 (246)、数据的地理范围、空间参考系统（“SRS WKT”）以及每个属性字段的类型信息。例如，FIPS: String (2.0)表示FIPS字符字段的最大长度为 2。同样，LON: Real (8.3)是一个浮点字段，最多可保存 8 位数字，最多三位小数。
地理模型¶
定义地理模型¶
现在您已经使用ogrinfo检查了您的数据集，请创建一个 GeoDjango 模型来表示此数据
from django.contrib.gis.db import models
class WorldBorder(models.Model):
# Regular Django fields corresponding to the attributes in the
# world borders shapefile.
name = models.CharField(max_length=50)
area = models.IntegerField()
pop2005 = models.IntegerField("Population 2005")
fips = models.CharField("FIPS Code", max_length=2, null=True)
iso2 = models.CharField("2 Digit ISO", max_length=2)
iso3 = models.CharField("3 Digit ISO", max_length=3)
un = models.IntegerField("United Nations Code")
region = models.IntegerField("Region Code")
subregion = models.IntegerField("Sub-Region Code")
lon = models.FloatField()
lat = models.FloatField()
# GeoDjango-specific: a geometry field (MultiPolygonField)
mpoly = models.MultiPolygonField()
# Returns the string representation of the model.
def __str__(self):
return self.name
请注意，models模块是从django.contrib.gis.db导入的。
几何字段的默认空间参考系统是 WGS84（这意味着SRID 为 4326）——换句话说，字段坐标以经度、纬度对表示，单位为度。要使用不同的坐标系，请使用srid参数设置几何字段的 SRID。使用表示坐标系 EPSG 代码的整数。
运行migrate¶
定义模型后，您需要将其与数据库同步。首先，创建一个数据库迁移
$ python manage.py makemigrations
Migrations for 'world':
world/migrations/0001_initial.py:
+ Create model WorldBorder
...\> py manage.py makemigrations
Migrations for 'world':
world/migrations/0001_initial.py:
+ Create model WorldBorder
让我们看看将为WorldBorder模型生成表的 SQL。
$ python manage.py sqlmigrate world 0001
...\> py manage.py sqlmigrate world 0001
此命令应产生以下输出
BEGIN;
--
-- Create model WorldBorder
--
CREATE TABLE "world_worldborder" (
"id" bigint NOT NULL PRIMARY KEY GENERATED BY DEFAULT AS IDENTITY,
"name" varchar(50) NOT NULL,
"area" integer NOT NULL,
"pop2005" integer NOT NULL,
"fips" varchar(2) NOT NULL,
"iso2" varchar(2) NOT NULL,
"iso3" varchar(3) NOT NULL,
"un" integer NOT NULL,
"region" integer NOT NULL,
"subregion" integer NOT NULL,
"lon" double precision NOT NULL,
"lat" double precision NOT NULL
"mpoly" geometry(MULTIPOLYGON,4326) NOT NULL
)
;
CREATE INDEX "world_worldborder_mpoly_id" ON "world_worldborder" USING GIST ("mpoly");
COMMIT;
如果看起来正确，请运行migrate以在数据库中创建此表
$ python manage.py migrate
Operations to perform:
Apply all migrations: admin, auth, contenttypes, sessions, world
Running migrations:
...
Applying world.0001_initial... OK
...\> py manage.py migrate
Operations to perform:
Apply all migrations: admin, auth, contenttypes, sessions, world
Running migrations:
...
Applying world.0001_initial... OK
导入空间数据¶
本节将向您展示如何使用LayerMapping 数据导入实用程序通过 GeoDjango 模型将世界边界 shapefile 导入数据库。
将数据导入空间数据库的方法有很多——除了 GeoDjango 中包含的工具外，您还可以使用以下工具：
ogr2ogr：GDAL 附带的一个命令行实用程序，可以将许多矢量数据格式导入 PostGIS、MySQL 和 Oracle 数据库。
shp2pgsql：此实用程序包含在 PostGIS 中，用于将 ESRI shapefile 导入 PostGIS。
GDAL 接口¶
前面，您使用ogrinfo检查了世界边界 shapefile 的内容。GeoDjango 还包含一个 Pythonic 接口，用于 GDAL 的强大 OGR 库，该库可以处理 OGR 支持的所有矢量数据源。
首先，调用 Django shell
$ python manage.py shell
...\> py manage.py shell
如果您在教程前面下载了世界边界数据，那么您可以使用 Python 的pathlib.Path确定其路径
>>> from pathlib import Path
>>> import world
>>> world_shp = Path(world.__file__).resolve().parent / "data" / "TM_WORLD_BORDERS-0.3.shp"
现在，使用 GeoDjango 的DataSource接口打开世界边界 shapefile
>>> from django.contrib.gis.gdal import DataSource
>>> ds = DataSource(world_shp)
>>> print(ds)
/ ... /geodjango/world/data/TM_WORLD_BORDERS-0.3.shp (ESRI Shapefile)
数据源对象可以具有不同的地理要素图层；但是，shapefile 只允许有一个图层
>>> print(len(ds))
1
>>> lyr = ds[0]
>>> print(lyr)
TM_WORLD_BORDERS-0.3
您可以看到图层的几何类型以及它包含的要素数量
>>> print(lyr.geom_type)
Polygon
>>> print(len(lyr))
246
注意
不幸的是，shapefile 数据格式不允许对几何类型进行更具体的指定。这个 shapefile，像许多其他 shapefile 一样，实际上包含的是 MultiPolygon 几何图形，而不是 Polygon。在模型中使用更通用的字段类型非常重要：GeoDjango 的 MultiPolygonField 可以接受 Polygon 几何图形，但 PolygonField 却不能接受 MultiPolygon 类型几何图形。这就是上面定义的 WorldBorder 模型使用 MultiPolygonField 的原因。
Layer 也可能具有与其关联的空间参考系统。如果它有，则 srs 属性将返回一个 SpatialReference 对象。
>>> srs = lyr.srs
>>> print(srs)
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.0174532925199433,
AUTHORITY["EPSG","9122"]],
AXIS["Latitude",NORTH],
AXIS["Longitude",EAST],
AUTHORITY["EPSG","4326"]]
>>> srs.proj # PROJ representation
'+proj=longlat +datum=WGS84 +no_defs'
此 shapefile 使用流行的 WGS84 空间参考系统——换句话说，数据使用经度、纬度对，单位为度。
此外，shapefile 还支持可能包含附加数据的属性字段。以下是世界边界图层上的字段。
>>> print(lyr.fields)
['FIPS', 'ISO2', 'ISO3', 'UN', 'NAME', 'AREA', 'POP2005', 'REGION', 'SUBREGION', 'LON', 'LAT']
以下代码将允许您检查与每个字段关联的 OGR 类型（例如整数或字符串）。
>>> [fld.__name__ for fld in lyr.field_types]
['OFTString', 'OFTString', 'OFTString', 'OFTInteger', 'OFTString', 'OFTInteger', 'OFTInteger64', 'OFTInteger', 'OFTInteger', 'OFTReal', 'OFTReal']
您可以迭代图层中的每个要素，并从要素的几何图形（通过 geom 属性访问）以及要素的属性字段（其**值**通过 get() 方法访问）中提取信息。
>>> for feat in lyr:
... print(feat.get("NAME"), feat.geom.num_points)
...
Guernsey 18
Jersey 26
South Georgia South Sandwich Islands 338
Taiwan 363
Layer 对象可以被切片。
>>> lyr[0:2]
[<django.contrib.gis.gdal.feature.Feature object at 0x2f47690>, <django.contrib.gis.gdal.feature.Feature object at 0x2f47650>]
并且可以通过其要素 ID 检索单个要素。
>>> feat = lyr[234]
>>> print(feat.get("NAME"))
San Marino
边界几何图形可以导出为 WKT 和 GeoJSON。
>>> geom = feat.geom
>>> print(geom.wkt)
POLYGON ((12.415798 43.957954,12.450554 ...
>>> print(geom.json)
{ "type": "Polygon", "coordinates": [ [ [ 12.415798, 43.957954 ], [ 12.450554, 43.979721 ], ...
LayerMapping¶
要导入数据，请在 Python 脚本中使用 LayerMapping。在 world 应用程序中创建一个名为 load.py 的文件，其中包含以下代码。
from pathlib import Path
from django.contrib.gis.utils import LayerMapping
from .models import WorldBorder
world_mapping = {
"fips": "FIPS",
"iso2": "ISO2",
"iso3": "ISO3",
"un": "UN",
"name": "NAME",
"area": "AREA",
"pop2005": "POP2005",
"region": "REGION",
"subregion": "SUBREGION",
"lon": "LON",
"lat": "LAT",
"mpoly": "MULTIPOLYGON",
}
world_shp = Path(__file__).resolve().parent / "data" / "TM_WORLD_BORDERS-0.3.shp"
def run(verbose=True):
lm = LayerMapping(WorldBorder, world_shp, world_mapping, transform=False)
lm.save(strict=True, verbose=verbose)
关于正在发生的事情的一些说明。
world_mapping字典中的每个键对应于WorldBorder模型中的一个字段。其值是从中加载数据的 shapefile 字段的名称。几何字段的键
mpoly为MULTIPOLYGON，这是 GeoDjango 将导入该字段的几何类型。即使 shapefile 中的简单多边形也会在插入数据库之前自动转换为集合。shapefile 的路径不是绝对路径——换句话说，如果您将
world应用程序（包含data子目录）移动到其他位置，脚本仍然可以工作。transform关键字设置为False，因为 shapefile 中的数据不需要转换——它已经位于 WGS84 (SRID=4326) 中。
之后，从 geodjango 项目目录调用 Django shell。
$ python manage.py shell
...\> py manage.py shell
接下来，导入 load 模块，调用 run 例程，并观察 LayerMapping 完成工作。
>>> from world import load
>>> load.run()
尝试使用 ogrinspect¶
既然您已经了解了如何使用 LayerMapping 数据导入实用程序 定义地理模型和导入数据，那么可以使用 ogrinspect 管理命令进一步自动化此过程。ogrinspect 命令会检查受 GDAL 支持的矢量数据源（例如 shapefile），并自动生成模型定义和 LayerMapping 字典。
该命令的一般用法如下所示。
$ python manage.py ogrinspect [options] <data_source> <model_name> [options]
...\> py manage.py ogrinspect [options] <data_source> <model_name> [options]
data_source 是受 GDAL 支持的数据源的路径，model_name 是要用于模型的名称。命令行选项可用于进一步定义模型的生成方式。
例如，以下命令几乎可以自动重现上面创建的 WorldBorder 模型和映射字典。
$ python manage.py ogrinspect world/data/TM_WORLD_BORDERS-0.3.shp WorldBorder \
--srid=4326 --mapping --multi
...\> py manage.py ogrinspect world\data\TM_WORLD_BORDERS-0.3.shp WorldBorder \
--srid=4326 --mapping --multi
关于上面给出的命令行选项的一些说明。
--srid=4326选项设置地理字段的 SRID。--mapping选项告诉ogrinspect也生成一个与LayerMapping一起使用的映射字典。--multi选项指定地理字段为MultiPolygonField而不是PolygonField。
该命令会生成以下输出，可以直接复制到 GeoDjango 应用程序的 models.py 中。
# This is an auto-generated Django model module created by ogrinspect.
from django.contrib.gis.db import models
class WorldBorder(models.Model):
fips = models.CharField(max_length=2)
iso2 = models.CharField(max_length=2)
iso3 = models.CharField(max_length=3)
un = models.IntegerField()
name = models.CharField(max_length=50)
area = models.IntegerField()
pop2005 = models.IntegerField()
region = models.IntegerField()
subregion = models.IntegerField()
lon = models.FloatField()
lat = models.FloatField()
geom = models.MultiPolygonField(srid=4326)
# Auto-generated `LayerMapping` dictionary for WorldBorder model
worldborders_mapping = {
"fips": "FIPS",
"iso2": "ISO2",
"iso3": "ISO3",
"un": "UN",
"name": "NAME",
"area": "AREA",
"pop2005": "POP2005",
"region": "REGION",
"subregion": "SUBREGION",
"lon": "LON",
"lat": "LAT",
"geom": "MULTIPOLYGON",
}
空间查询¶
空间查找¶
GeoDjango 将空间查找添加到 Django ORM。例如，您可以找到 WorldBorder 表中包含特定点的国家/地区。首先，启动管理 shell。
$ python manage.py shell
...\> py manage.py shell
现在，定义一个兴趣点 [3]
>>> pnt_wkt = "POINT(-95.3385 29.7245)"
pnt_wkt 字符串表示经度 -95.3385 度，纬度 29.7245 度的点。几何图形采用称为“众所周知文本 (WKT)”的格式，这是开放地理空间联盟 (OGC) 发布的一种标准。[4] 导入 WorldBorder 模型，并使用 pnt_wkt 作为参数执行 contains 查找。
>>> from world.models import WorldBorder
>>> WorldBorder.objects.filter(mpoly__contains=pnt_wkt)
<QuerySet [<WorldBorder: United States>]>
在这里，您检索到一个仅包含一个模型的 QuerySet：美国的边界（正是您所期望的）。
同样，您也可以使用 GEOS 几何对象。在这里，您可以将 intersects 空间查找与 get 方法结合使用，以仅检索圣马力诺的 WorldBorder 实例，而不是查询集。
>>> from django.contrib.gis.geos import Point
>>> pnt = Point(12.4604, 43.9420)
>>> WorldBorder.objects.get(mpoly__intersects=pnt)
<WorldBorder: San Marino>
contains 和 intersects 查找只是可用查询的一个子集——GeoDjango 数据库 API 文档中有更多内容。
自动空间转换¶
在执行空间查询时，如果几何图形位于不同的坐标系中，GeoDjango 会自动转换几何图形。在以下示例中，坐标将以 EPSG SRID 32140 表示，这是一个**仅限**德克萨斯州南部的坐标系，单位为**米**，而不是度。
>>> from django.contrib.gis.geos import GEOSGeometry, Point
>>> pnt = Point(954158.1, 4215137.1, srid=32140)
请注意，pnt 也可以使用 EWKT 构造，EWKT 是包含 SRID 的 WKT 的“扩展”形式。
>>> pnt = GEOSGeometry("SRID=32140;POINT(954158.1 4215137.1)")
GeoDjango 的 ORM 会自动将几何值包装在转换 SQL 中，从而允许开发人员在更高的抽象级别上工作。
>>> qs = WorldBorder.objects.filter(mpoly__intersects=pnt)
>>> print(qs.query) # Generating the SQL
SELECT "world_worldborder"."id", "world_worldborder"."name", "world_worldborder"."area",
"world_worldborder"."pop2005", "world_worldborder"."fips", "world_worldborder"."iso2",
"world_worldborder"."iso3", "world_worldborder"."un", "world_worldborder"."region",
"world_worldborder"."subregion", "world_worldborder"."lon", "world_worldborder"."lat",
"world_worldborder"."mpoly" FROM "world_worldborder"
WHERE ST_Intersects("world_worldborder"."mpoly", ST_Transform(%s, 4326))
>>> qs # printing evaluates the queryset
<QuerySet [<WorldBorder: United States>]>
原始查询
使用 原始查询 时，必须包装几何字段，以便 GEOS 可以识别字段值。
>>> from django.db import connection
>>> # or if you're querying a non-default database:
>>> from django.db import connections
>>> connection = connections["your_gis_db_alias"]
>>> City.objects.raw(
... "SELECT id, name, %s as point from myapp_city" % (connection.ops.select % "point")
... )
只有在确切知道自己在做什么的情况下，才应使用原始查询。
延迟几何¶
GeoDjango 以标准化的文本表示形式加载几何图形。第一次访问几何字段时，GeoDjango 会创建一个 GEOSGeometry 对象，从而提供强大的功能，例如流行地理空间格式的序列化属性。
>>> sm = WorldBorder.objects.get(name="San Marino")
>>> sm.mpoly
<MultiPolygon object at 0x24c6798>
>>> sm.mpoly.wkt # WKT
MULTIPOLYGON (((12.4157980000000006 43.9579540000000009, 12.4505540000000003 43.9797209999999978, ...
>>> sm.mpoly.wkb # WKB (as Python binary buffer)
<read-only buffer for 0x1fe2c70, size -1, offset 0 at 0x2564c40>
>>> sm.mpoly.geojson # GeoJSON
'{ "type": "MultiPolygon", "coordinates": [ [ [ [ 12.415798, 43.957954 ], [ 12.450554, 43.979721 ], ...
这包括对 GEOS 库提供的全部高级几何运算的访问。
>>> pnt = Point(12.4604, 43.9420)
>>> sm.mpoly.contains(pnt)
True
>>> pnt.contains(sm.mpoly)
False
地理标注¶
GeoDjango 还提供了一组地理标注，用于计算距离和其他多种操作（交集、差集等）。请参阅地理数据库函数文档。
在地图上显示您的数据¶
地理管理¶
Django 的管理应用程序支持编辑几何字段。
基础知识¶
Django 管理界面允许用户在 JavaScript 可平移地图上创建和修改几何图形（由OpenLayers提供支持）。
让我们开始吧。在world应用程序内创建一个名为admin.py的文件，其中包含以下代码
from django.contrib.gis import admin
from .models import WorldBorder
admin.site.register(WorldBorder, admin.ModelAdmin)
接下来，修改geodjango应用程序文件夹中的urls.py文件，如下所示
from django.contrib import admin
from django.urls import include, path
urlpatterns = [
path("admin/", admin.site.urls),
]
创建管理员用户
$ python manage.py createsuperuser
...\> py manage.py createsuperuser
接下来，启动 Django 开发服务器
$ python manage.py runserver
...\> py manage.py runserver
最后，浏览到http://localhost:8000/admin/，并使用您刚刚创建的用户登录。浏览任何WorldBorder条目——可以通过点击多边形并拖动顶点到所需位置来编辑边界。
GISModelAdmin¶
使用GISModelAdmin，GeoDjango 在管理界面中使用OpenStreetMap图层。这提供了比ModelAdmin（使用矢量地图级别 0 WMS 数据集，托管在OSGeo）更多的上下文（包括街道和道路详细信息）。
必须安装 PROJ 基准面转换文件（有关详细信息，请参阅PROJ 安装说明）。
如果您满足此要求，请在您的admin.py文件中使用GISModelAdmin选项类
admin.site.register(WorldBorder, admin.GISModelAdmin)
脚注
