Analisando dias atípicos e changesets suspeitos
naoliv님이 Brazilian Portuguese (Português do Brasil)로 2015년 4월 25일에 게시함. 최근 2015년 6월 10일에 업데이트됨.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 :-)
토론
2015년 4월 25일 02:58에 pizzaiolo님의 의견
Bacana! Bom saber.
2015년 5월 2일 20:14에 Alexandre Magno님의 의견
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
2015년 5월 2일 20:20에 naoliv님의 의견
Passa para a direita ou esquerda o segundo gráfico que chega nessa visualização.
2015년 5월 2일 20:23에 Alexandre Magno님의 의견
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
.2015년 5월 2일 20:25에 naoliv님의 의견
Vão ser analizados todos os changesets, independente de ser pico ou não.
2015년 5월 2일 20:39에 Alexandre Magno님의 의견
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], ...
2015년 5월 2일 20:40에 Alexandre Magno님의 의견
Que tipo de coisas você quer perguntar nessa análise?
2015년 5월 2일 22:07에 Alexandre Magno님의 의견
Em Python,
datetime.datetime.fromtimestamp(1320019201000/1000).strftime('%Y-%m-%d %H:%M:%S')
retorna'2011-10-30 21:00:01'
.