PyQGISでフィールドを追加して式を計算しますか？

PyQGISを使用して新しいフィールドを追加し、各フィーチャの値を計算したいと思います。Field Calculatorオプションに似ています。

'Field Calculator'式は次のとおりです。 y(start_point($geometry))

from PyQt4.QtCore import QVariant
from qgis.core import QgsField, QgsExpression, QgsFeature
vl = iface.activeLayer()

vl.startEditing()

#step 1
myField = QgsField( 'myNewColumn', QVariant.Float )
vl.addAttribute( myField )
idx = vl.fieldNameIndex( 'myNewColumn' )

#step 2
e = QgsExpression( 'y(start_point($geometry))' )
e.prepare( vl.pendingFields() )

for f in vl.getFeatures():
    f[idx] = e.evaluate( f )
    vl.updateFeature( f )

vl.commitChanges()

これは私が得るエラーです：

Traceback (most recent call last):
  File "<input>", line 1, in <module>
  File "/var/folders/0r/q6nxg6w54fv4l7c8gksb2t500000gn/T/tmp9dosIe.py", line 30, in <module>
    f[idx] = e.evaluate( f )
KeyError: '-1'

表示されるエラーでは、フィールドインデックスが-1であるため、新しいフィールドが属性テーブルに見つかりませんでした。

これは、次の理由で発生している可能性があります。

1. のQVariant.Double代わりに使用する必要がありますQVariant.Float。
2. 新しい列インデックスを要求する前に、新しいフィールドをレイヤープロバイダーにコミットしていません。
3. のインデックスを要求してmyNewColumnいますが、プロバイダーはフィールド名として10文字しか保存できないため、保存されていますmyNewColum（finalがありませんn）。（答えを試してみているときに私に起こっただけです）

代わりにこれを試してください：

#step 1
myField = QgsField( 'newColumn', QVariant.Double )
vl.dataProvider().addAttributes([myField])
vl.updateFields()
idx = vl.fieldNameIndex( 'newColumn' )

これidxで、とは異なるはず-1です。確認できます。

if idx != -1:
    print "Field found!"

ちなみに、#step 1編集ブロックからコードを実行できます。