QGISで大量データを扱う方法

layer = iface.activeLayer()
provider = layer.dataProvider()

# 空間インデックスを作成
if provider.capabilities() & QgsVectorDataProvider.CreateSpatialIndex:
    provider.createSpatialIndex()
    print("空間インデックスを作成しました")

import processing

result = processing.run("native:simplifygeometries", {
    'INPUT': layer,
    'METHOD': 0,  # Douglas-Peucker
    'TOLERANCE': 10,  # メートル
    'OUTPUT': 'memory:'
})

# 頂点数の確認
for feature in layer.getFeatures():
    geom = feature.geometry()
    if geom:
        vertex_count = sum(1 for _ in geom.vertices())
        print(f"ID {feature.id()}: {vertex_count} 頂点")

# 必要なフィールドのみ保持
fields_to_keep = ['id', 'name', 'area']

processing.run("native:retainfields", {
    'INPUT': layer,
    'FIELDS': fields_to_keep,
    'OUTPUT': 'memory:'
})

# 使っていないレイヤを削除
for layer_id, layer in QgsProject.instance().mapLayers().items():
    if not layer.name() in ['必要なレイヤ1', '必要なレイヤ2']:
        QgsProject.instance().removeMapLayer(layer_id)

# 特定範囲のみ読み込み
uri = f"{filepath}|layername={layername}|subset=ST_Within(geometry, ST_GeomFromText('POLYGON(...)'))"
layer = QgsVectorLayer(uri, "layer", "ogr")

# 処理対象の範囲を限定
extent = QgsRectangle(139.7, 35.6, 139.8, 35.7)
request = QgsFeatureRequest().setFilterRect(extent)

for feature in layer.getFeatures(request):
    # 処理
    pass

# 空間インデックスを作成
spatial_index = QgsSpatialIndex(layer.getFeatures())

# 特定範囲内のフィーチャIDを取得（高速）
ids = spatial_index.intersects(QgsRectangle(139.7, 35.6, 139.8, 35.7))

# IDからフィーチャを取得
for fid in ids:
    feature = layer.getFeature(fid)
    # 処理

# 悪い例：毎回コミット
for feature in layer.getFeatures():
    layer.startEditing()
    layer.changeAttributeValue(feature.id(), 0, "value")
    layer.commitChanges()  # 毎回コミット → 遅い

# 良い例：まとめてコミット
layer.startEditing()
for feature in layer.getFeatures():
    layer.changeAttributeValue(feature.id(), 0, "value")
layer.commitChanges()  # 最後に1回 → 速い

from qgis.core import QgsTask

class LargeDataTask(QgsTask):
    def run(self):
        total = self.layer.featureCount()
        for i, feature in enumerate(self.layer.getFeatures()):
            if self.isCanceled():
                return False

            # 処理
            self.setProgress(i / total * 100)

        return True