{"id":7776,"date":"2018-12-04T07:30:14","date_gmt":"2018-12-04T10:30:14","guid":{"rendered":"http:\/\/www.fernandoquadro.com.br\/html\/?p=7776"},"modified":"2018-12-05T08:06:55","modified_gmt":"2018-12-05T11:06:55","slug":"desenvolvimento-de-plugins-do-qgis-vizinho-mais-proximo-parte-3","status":"publish","type":"post","link":"https:\/\/www.fernandoquadro.com.br\/html\/2018\/12\/04\/desenvolvimento-de-plugins-do-qgis-vizinho-mais-proximo-parte-3\/","title":{"rendered":"Desenvolvimento de plugins do QGIS: Vizinho mais pr\u00f3ximo"},"content":{"rendered":"

No \u00faltimo post<\/a> foi descrito o b\u00e1sico da manipula\u00e7\u00e3o de camadas vetoriais. Como o nosso objetivo \u00e9 criar um plug-in personalizado totalmente funcional com base em uma dist\u00e2ncia, gostaria de discutir a indexa\u00e7\u00e3o espacial e a an\u00e1lise do vizinho mais pr\u00f3ximo.<\/p>\n

\"\"<\/center><\/p>\n

A figura acima ilustra a tarefa que pode ser resolvida apenas usando a API do QGIS. Imagine que voc\u00ea tenha uma camada de origem com um atributo preenchido com valores. Voc\u00ea recebe uma camada de destino tamb\u00e9m, infelizmente, os valores nesta camada est\u00e3o faltando (n\u00e3o t\u00e3o raros no mundo GIS, certo?). No entanto, voc\u00ea sabe que o valor ausente de cada recurso na camada de destino pode ser preenchido pelo valor de seu vizinho mais pr\u00f3ximo da camada de origem. Como voc\u00ea faz isso?<\/p>\n

1. GERANDO DADOS FICT\u00cdCIOS<\/strong><\/p>\n

Vamos criar dois conjuntos de dados na mem\u00f3ria com atributos id e value. Ambos ter\u00e3o dez recursos.<\/p>\n

\r\nfrom qgis.core import QgsMapLayerRegistry, QgsVectorLayer, QgsFeature, QgsGeometry, QgsPoint, QgsSpatialIndex\r\nfrom qgis.utils import iface\r\n\r\nsource_layer = QgsVectorLayer(\"point?crs=epsg:4326&field=id:integer&field=value:integer\", \"Source layer\", \"memory\")\r\ntarget_layer = QgsVectorLayer(\"point?crs=epsg:4326&field=id:integer&field=value:integer\", \"Target layer\", \"memory\")\r\n\r\ndef create_dummy_data():\r\n\r\n    source_layer.startEditing()\r\n    target_layer.startEditing()\r\n\r\n    feature = QgsFeature(source_layer.pendingFields())\r\n\r\n    for i in range(10):\r\n        feature.setGeometry(QgsGeometry.fromPoint(QgsPoint(i, i)))\r\n        feature.setAttribute(\"id\", i)\r\n        feature.setAttribute(\"value\", i)\r\n        source_layer.addFeature(feature)\r\n\r\n    feature = QgsFeature(source_layer.pendingFields())\r\n\r\n    for i in range(10):\r\n        feature.setGeometry(QgsGeometry.fromPoint(QgsPoint(i + i, i)))\r\n        feature.setAttribute(\"id\", i)\r\n        target_layer.addFeature(feature)\r\n\r\n    source_layer.commitChanges()\r\n    target_layer.commitChanges()\r\n\r\n    QgsMapLayerRegistry.instance().addMapLayer(source_layer)\r\n    QgsMapLayerRegistry.instance().addMapLayer(target_layer)\r\n\r\ncreate_dummy_data()\r\n<\/pre>\n
2. ESCREVENDO VALORES DO VIZINHO MAIS PR\u00d3XIMO<\/strong><\/p>\n
A an\u00e1lise real do vizinho mais pr\u00f3ximo pode ser feita em dez linhas de c\u00f3digo! Primeiro, o qgis.core.QgsSpatialIndex \u00e9 constru\u00eddo a partir de todos os recursos source_layer<\/strong>. Ent\u00e3o, voc\u00ea itera sobre os recursos target_layer<\/strong> e para cada um deles, obt\u00e9m apenas um ( nearestNeighbor(f.geometry().asPoint(), 1)[0]) vizinho mais pr\u00f3ximo. Por fim, voc\u00ea apenas grava o valor do atributo do vizinho mais pr\u00f3ximo na camada de destino e confirma as altera\u00e7\u00f5es. <\/p>\n
\r\ndef write_values_from_nn():\r\n    source_layer_index = QgsSpatialIndex(source_layer.getFeatures())\r\n    source_layer_features = {feature.id(): feature for (feature) in source_layer.getFeatures()}\r\n    target_layer_features = target_layer.getFeatures()\r\n\r\n    target_layer.startEditing()\r\n\r\n    for f in target_layer_features:\r\n        nearest = source_layer_index.nearestNeighbor(f.geometry().asPoint(), 1)[0]\r\n        value = source_layer_features[nearest].attribute(\"value\")\r\n        target_layer.changeAttributeValue(f.id(), 1, value)\r\n\r\n    target_layer.commitChanges()\r\n\r\nwrite_values_from_nn()\r\n<\/pre>\n
3. O QUE VEM POR AI<\/strong><\/p>\n
Estamos um pouco mais perto do nosso objetivo. O que esta faltando<\/a>?<\/p>\n