Como criar programaticamente e adicionar recursos a uma camada de memória no QGIS 1.9?

Eu tinha um plug-in funcional no QGIS 1.8 que lia dados de um banco de dados MSAccess e o adicionava a uma série de camadas de memória. Há algum processamento envolvido nesse ínterim, então não acho que simplesmente usar o QGIS para ler diretamente do banco de dados seja uma opção.

Gostaria de passar do QGIS 1.8 para 1.9 (principalmente devido à qualidade aprimorada do compositor de impressão). O plug-in não funciona com a nova API.

Eu tentei uma variedade de métodos que surgiram nas pesquisas do Google. Um, modificando o código abaixo - de http://www.qgis.org/pyqgis-cookbook/vector.html#memory-provider , isto é, adicionando geometria e atributos ao dataprovider e atualizando a camada - para se adequar à nova API trabalhada um pouco, mas os atributos não estavam visíveis até que eu entrei no modo de edição manualmente (semelhante a http://hub.qgis.org/issues/3713 ). Uma abordagem alternativa, detalhada na resposta nº 1 do link acima, adicionou a camada e os atributos corretamente, mas não foi possível adicionar recursos à camada.

Dado que deve ser uma tarefa bastante simples, espero que alguém aqui possa oferecer um exemplo prático de como isso deve ser feito. (PS: Eu não sou um programador profissional e a maior parte da minha codificação é bastante grosseira - agradeço qualquer orientação, mas peço que você desculpe alguma ignorância da minha parte)

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

Graças a gsherman acima, o exemplo do plugin PinPoint é perfeito.

Pelo que entendi, o processo é:

1. Crie a camada com os atributos na cadeia de construção
2. Adicione a referida camada ao registro do mapa
3. Comece a editar nessa camada
4. Adicionar recursos e confirmar as alterações

Aqui está um extrato do meu código que agora funciona.

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

Com base na resposta de Adam Bioletti, testes adicionais ao processo descrito mostram que o único requisito essencial é começar a editar a camada de memória antes de fazer alterações, como criar atributos e recursos, e depois confirmar as alterações. Isso pode ser feito antes de adicionar a camada ao mapeamento do registro.

Aqui está uma atualização do código do Cookbook que funciona com a 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)