Як я можу програмно створювати та додавати функції до шару пам'яті в 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/isissue/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 API:

# 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)