Analisando dias atípicos e changesets suspeitos
Yuzen-it-id naoliv af 25 April 2015 s Brazilian Portuguese (Português do Brasil) Last updated on 10 June 2015.Para quem não conhece, temos um site com bastante tipo de estatísticas e visualizações de dados em http://resultmaps.neis-one.org/
Uma delas é a parte que mostra as estatísticas diárias de um determinado país, como a do Brasil
Dando uma breve olhada nos gráficos a gente acaba reparado em alguns picos bem atípicos para o padrão de edição em nosso país:
O pico maior de adição de nós ocorreu em 27 de fevereiro de 2014 (mais de 1 ano atrás), com 89404 nós adicionados (seguido de várias remoções).
Os dois próximos picos em verde ocorrem em 19 de janeiro desse ano, com 50133 nós adicionados, e em 9 de março, com 50862 nós.
Mas como encontrar o que foi alterado nesse dia sem precisar baixar o mundo inteiro ou ter uma base histórica do Brasil inteiro?
Através da API do OSM (devo uma cerveja pro pnorman pela sugestão)
Para obter os changesets de uma área precisamos do bbox e do período que desejamos.
O bbox conseguimos de forma bem fácil através do JOSM:
- baixar a relação do Brasil
- selecioná-la
Ctrl+i
- pegar o bbox calculado:
-73.9830625,-33.8689056,0.0,5.2842873
A data também temos (dia 27 de fevereiro de 2014), só precisando alterá-la para um formato reconhecido pela API.
Com um simples GET "http://api.openstreetmap.org/api/0.6/changesets/?bbox=-73.9830625,-33.8689056,0.0,5.2842873&time=20140226T000000,20140227T000000&closed=true"
temos então todos os changesets no Brasil de 00:00 do dia 26 de fevereiro de 2014 até 00:00 do dia 27 de fevereiro de 2014. Ou melhor, deveríamos ter…
A API retorna apenas 100 resultados (portanto se o período tiver mais de 100 edições, não retornará todas).
Quebrando as requisições em períodos menores, de 1 hora (ou, em alguns casos, de meia hora), temos então todos os changesets de um dia inteiro.
Com cada id de changeset fica fácil obter os dados deles e analisar o que aconteceu: quem editou no Brasil, que horas, qual editor utilizou, quantos objetos adicionou, removeu, modificou, etc.
O bom disso é que abre muitas possibilidades: podemos analisar qualquer dia anterior onde há picos no gráfico (seja por adição ou remoção de objetos) ou então, futuramente, analisar os dados diariamente na busca, também, de edições suspeitas.
Não é nada complicado e dá para ser todo automatizado o processo.
Por essa pequena análise pudemos ver que um usuário alemão é muito ativo no Brasil (ou era naquela época, pelo menos), que houve uma modificação automatizada de nomes de ruas e também, infelizmente, uma importação…
A importação precisará ser analisada, junto com os outros changesets do usuário.
Outros dias atípicos também serão verificados na medida do possível.
Aceito correções de ortografia, Português e pedidos para explicar algo melhor, desde que saiba explicar :-)
Discussion
Awennit n pizzaiolo di 25 April 2015 ɣef 02:58
Bacana! Bom saber.
Awennit n Alexandre Magno di 2 May 2015 ɣef 20:14
Como você obteve esse gráfico abrangendo vários meses?
Eu só consigo isto:
http://osmstats.neis-one.org/?item=countries&country=Brazil&date=2-4-2014
Awennit n naoliv di 2 May 2015 ɣef 20:20
Passa para a direita ou esquerda o segundo gráfico que chega nessa visualização.
Awennit n Alexandre Magno di 2 May 2015 ɣef 20:23
Na sua opinião, qual seria um critério quantificador adequado para detectar os picos? Olhando visualmente o gráfico, eu compreendo o que são os picos. Matematicamente, podemos ensaiar uma algoritmo na issue #3 do
osmdetective
.Awennit n naoliv di 2 May 2015 ɣef 20:25
Vão ser analizados todos os changesets, independente de ser pico ou não.
Awennit n Alexandre Magno di 2 May 2015 ɣef 20:39
Será que Pascal Neis não retornaria os hashmaps
data
do #graph10 que ele plota por Javascript?Pois isto aqui tem tudo que a mais simples análise automatizada precisaria:
var plot = $.plot($("#graph10"),[ { data: [[1320019201000,8041], ...
Awennit n Alexandre Magno di 2 May 2015 ɣef 20:40
Que tipo de coisas você quer perguntar nessa análise?
Awennit n Alexandre Magno di 2 May 2015 ɣef 22:07
Em Python,
datetime.datetime.fromtimestamp(1320019201000/1000).strftime('%Y-%m-%d %H:%M:%S')
retorna'2011-10-30 21:00:01'
.