在GeoDjango中导入空间数据时出错 - 对于mpoly字段的KeyError

时间:2015-06-28 05:19:19

标签: python django python-2.7 geodjango

我正在关注https://docs.djangoproject.com/en/1.8/ref/contrib/gis/tutorial/#importing-spatial-data上的教程,以便在我的机器上设置GeoDjango。但似乎那里存在一些问题。通过运行load.run()使用LayerMapping导入数据时,我收到以下错误:

Traceback (most recent call last):
  File "<console>", line 1, in <module>
  File "/home/ubuntu/src/django/world/load.py", line 23, in run
    lm = LayerMapping(WorldBorder, world_shp, world_mapping, transform=False, encoding='iso-8859-1')
  File "/home/ubuntu/Envs/vir-env/local/lib/python2.7/site-packages/django/contrib/gis/utils/layermapping.py", line 105, in __init__
    self.check_layer()
  File "/home/ubuntu/Envs/vir-env/local/lib/python2.7/site-packages/django/contrib/gis/utils/layermapping.py", line 178, in check_layer
    ogr_field_types = self.layer.field_types
  File "/home/ubuntu/Envs/vir-env/local/lib/python2.7/site-packages/django/contrib/gis/gdal/layer.py", line 153, in field_types
    for i in range(self.num_fields)]
KeyError: 12

然后我发现,没有&#39; MULTIPOLYGON&#39; .shp文件中的字段:

>>> from django.contrib.gis.gdal import DataSource
>>> ds = DataSource('world/data/TM_WORLD_BORDERS-0.3.shp')
>>> layer = ds[0]
>>> layer.fields
[u'FIPS', u'ISO2', u'ISO3', u'UN', u'NAME', u'AREA', u'POP2005', u'REGION', u'SUBREGION', u'LON', u'LAT']

但它在WorldBorder模型中,类型为MultiPolygonField。 因此,绝对在world_mapping文件中,'mpoly': 'MULTIPOLYGON'映射的导入将失败。还有其他人遇到过这个问题吗?我希望如此,因为我一步一步地按照教程。但它对此类问题没有任何说明。如果我通过删除mpoly映射加载数据会产生什么影响?

这是我的load.py文件:

  1 import os
  2 from django.contrib.gis.utils import LayerMapping
  3 from models import WorldBorder
  4
  5 world_mapping = {
  6     'fips' : 'FIPS',
  7     'iso2' : 'ISO2',
  8     'iso3' : 'ISO3',
  9     'un' : 'UN',
 10     'name' : 'NAME',
 11     'area' : 'AREA',
 12     'pop2005' : 'POP2005',
 13     'region' : 'REGION',
 14     'subregion' : 'SUBREGION',
 15     'lon' : 'LON',
 16     'lat' : 'LAT',
 17     'mpoly' : 'MULTIPOLYGON',
 18 }
 19
 20 world_shp = os.path.abspath(os.path.join(os.path.dirname(__file__), 'data/TM_WORLD_BORDERS-0.3.shp'))
 21
 22 def run(verbose=True):
 23     lm = LayerMapping(WorldBorder, world_shp, world_mapping, transform=False, encoding='iso-8859-1')
 24
 25     lm.save(strict=True, verbose=verbose)

只是更新: 在浏览源代码后,通过堆栈跟踪,我发现我无法访问field_types模块的layer属性。所以,从python shell,当我访问该属性时,我得到同样的错误:

>>> from django.contrib.gis.gdal import DataSource
>>> ds = DataSource(wshp)
>>> layer = ds[0]
>>> layer.fields
[u'FIPS', u'ISO2', u'ISO3', u'UN', u'NAME', u'AREA', u'POP2005', u'REGION', u'SUBREGION', u'LON', u'LAT']
>>> layer.field_types
Traceback (most recent call last):
  File "<console>", line 1, in <module>
  File "/home/ubuntu/Envs/rj-venv/local/lib/python2.7/site-packages/django/contrib/gis/gdal/layer.py", line 153, in field_types
    for i in range(self.num_fields)]
KeyError: 12

现在,这很奇怪,因为我现在还从mpoly模型中删除了WorldBorder字段。

更新2:

在深入了解源代码后,我发现,在我的gdal版本中,OGDFieldTypes可能没有密钥12,如此处的源代码:https://github.com/django/django/blob/master/django/contrib/gis/gdal/field.py。 但它表示密钥12和13将可用于GDAL 2,这就是我已安装的密码。现在看来,图书馆之间似乎存在一些冲突。

我已经安装了以下库:

  • GEOS-3.4.2.tar.bz2
  • PROJ-datumgrid-1.5.tar.gz
  • PROJ-4.8.0.tar.gz
  • GDAL-2.0.0.tar.gz

PostGIS版本2.1.5安装在Amazon RDS实例中。

1 个答案:

答案 0 :(得分:7)

这里的问题是Django的version 1.8.2(撰写本文时的最新版本)不支持GDAL 2.0.0字段类型。对此的修复是在Django主代码分支中 - 但尚未在发行版中。

在编写本文时,您有两个选择 - 从Github上的最新源构建Django(即升级),或降级GDAL,以便本机代码called by this line不能将值12作为字段返回类型。

因为GDAL 2.0 can return 12 as a field type,使用Django 1.8.2代码看到的this dictionary查找失败了KeyErrorOn masterdict包含GDAL 2.0的正确字段 - 由this commit修正(在撰写本文时仅9天!)。

为了完整性 - 在SOPython聊天室中讨论了这个错误的诊断 - 指向书签对话的链接is here

<强> N.B。来自聊天 如果您安装了多个GDAL版本,则可能需要在django设置文件中添加GDAL_LIBRARY_PATH以指向libgdal.so的正确版本(或者我猜是相应的.dll,如果你在Windows上)