Как я могу программно создавать и добавлять объекты в слой памяти в QGIS 1.9?

У меня был рабочий плагин в QGIS 1.8, который считывал данные из базы данных MSAccess и добавлял их в ряд слоев памяти. В промежутке между тем происходит некоторая обработка, поэтому я не думаю, что просто использовать QGIS для чтения непосредственно из базы данных.

Я бы хотел перейти с QGIS 1.8 на 1.9 (в первую очередь из-за улучшенного качества композитора печати). Плагин не работает с новым API.

Я пробовал различные методы, которые были найдены в поиске Google. Во-первых, изменение кода ниже - с http://www.qgis.org/pyqgis-cookbook/vector.html#memory-provider , то есть добавление геометрии и атрибутов к провайдеру данных, а затем обновление слоя - в соответствии с новым API немного, но атрибуты не были видны, пока я вручную не вошел в режим редактирования (похоже на http://hub.qgis.org/issues/3713 ). Альтернативный подход, подробно описанный в ответе 1 вышеупомянутой ссылки, правильно добавил слой и атрибуты, но я не смог добавить объекты в слой.

Учитывая, что это должно быть довольно простой задачей, я надеюсь, что кто-то здесь может предложить рабочий пример того, как это должно быть сделано. (PS Я не профессиональный программист, и большая часть моего кодирования довольно грубая - я приветствую любые рекомендации, но прошу вас извинить за некоторое невежество с моей стороны)

# Receivers = a list of lists returned from a database query

# create layer
vl = QgsVectorLayer("Point", item, "memory")
pr = vl.dataProvider()

# add fields
pr.addAttributes( [ QgsField("Rec_No", QVariant.Int), QgsField("Include",  QVariant.String), QgsField("Label",  QVariant.String), QgsField("X", QVariant.Double),
                    QgsField("Y", QVariant.Double), QgsField("Z", QVariant.Double), QgsField("Height", QVariant.Double),
                    QgsField("Project_Re", QVariant.String), QgsField("NCA", QVariant.String),
                    QgsField("DayCrit", QVariant.Int), QgsField("EveCrit", QVariant.Int), QgsField("NightCrit", QVariant.Int) ] )

for i in range(len(Receivers)):          
  # add a feature
  fet = QgsFeature()
  X = Receivers[i][3]
  Y = Receivers[i][4]
  fet.setGeometry( QgsGeometry.fromPoint(QgsPoint(X,Y)) )

  # Details = a list of results returned from a database query specific to each result in 'Receivers'

  if Receivers[i][3] != 0:
    Include = 'Yes'
  else:
    Include = 'No'

  fet.setAttributeMap( { 0 : QVariant(Receivers[i][0]), 1 : QVariant(Include), 2 : QVariant(Receivers[i][2]),
                         3 : QVariant(Receivers[i][3]), 4 : QVariant(Receivers[i][4]), 5 : QVariant(Receivers[i][5]), 6 : QVariant(Receivers[i][6]),
                         7 : QVariant(Details[0]), 8 : QVariant(Details[1]), 9 : QVariant(Details[2]), 10 : QVariant(Details[3]), 11 : QVariant(Details[4]) } )
  pr.addFeatures( [ fet ] )

# add a style
vl.loadNamedStyle('C:/OSGeo4W/apps/qgis/python/plugins/Gopher2QGIS/styles/Receiver_Style.qml')

# update layer's extent when new features have been added
# because change of extent in provider is not propagated to the layer
vl.commitChanges()
vl.updateExtents()
vl.updateFieldMap()

QgsMapLayerRegistry.instance().addMapLayer(vl)

Благодаря gsherman выше пример плагина PinPoint идеален.

Насколько я понимаю, процесс таков:

1. Создайте слой с атрибутами в строительной строке
2. Добавить указанный слой в реестр карт
3. Начать редактирование на этом слое
4. Добавить функции и зафиксировать изменения

Вот выдержка из моего кода, который сейчас работает.

layer =  QgsVectorLayer(
          "Point?field=Rec_No:integer&field=Include:string(120)&field=Label:string(120)&field=X:double&field=Y:double&field=Z:double&field=Height:double&field=Project_Re:string(120)&field=NCA:string(120)&field=DayCrit:integer&field=EveCrit:integer&field=NightCrit:integer",
          item,
          "memory")
QgsMapLayerRegistry.instance().addMapLayer(layer)

# Receivers = as in the above example 'Receivers' is a list of results
for i in range(len(Receivers)):

  # add a feature
  feature = QgsFeature()

  X = Receivers[i][3]
  Y = Receivers[i][4]
  feature.setGeometry( QgsGeometry.fromPoint(QgsPoint(X,Y)) )

  # Details = as in the above example 'Details' is a list of results

  if Receivers[i][1] != 0:
    Include = 'Yes'
  else:
    Include = 'No'

  values = [ QVariant(Receivers[i][0]), QVariant(Include), QVariant(Receivers[i][2]),
                         QVariant(Receivers[i][3]), QVariant(Receivers[i][4]), QVariant(Receivers[i][5]), QVariant(Receivers[i][6]),
                         QVariant(Details[0]), QVariant(Details[1]), QVariant(Details[2]), QVariant(Details[3]), QVariant(Details[4]) ]

  feature.setAttributes(values)
  layer.startEditing()
  layer.addFeature(feature, True)
  layer.commitChanges()

Основываясь на ответе Адама Биолетти, дальнейшие тесты описанного процесса показывают, что единственное существенное требование - начать редактирование слоя памяти перед внесением каких-либо изменений, таких как создание атрибутов и функций, а затем зафиксировать изменения. Это можно сделать перед добавлением слоя на карту реестра.

Вот обновление кода поваренной книги, который работает с API 2.0:

# create layer
vl = QgsVectorLayer("Point", "temporary_points", "memory")
pr = vl.dataProvider()

# changes are only possible when editing the layer
vl.startEditing()
# add fields
pr.addAttributes([QgsField("name", QVariant.String),QgsField("age", QVariant.Int),QgsField("size", QVariant.Double)])

# add a feature
fet = QgsFeature()
fet.setGeometry(QgsGeometry.fromPoint(QgsPoint(10,10)))
fet.setAttributes(["Johny", 2, 0.3])
pr.addFeatures([fet])

# commit to stop editing the layer
vl.commitChanges()

# update layer's extent when new features have been added
# because change of extent in provider is not propagated to the layer
vl.updateExtents()

# add layer to the legend
QgsMapLayerRegistry.instance().addMapLayer(vl)