BR112017020838B1 - Método para estimar um número de pessoas reunidas em uma área de interesse, e, sistema acoplado com uma rede de telecomunicação sem fio para estimar um número de pessoas reunidas em uma área de interesse - Google Patents

Método para estimar um número de pessoas reunidas em uma área de interesse, e, sistema acoplado com uma rede de telecomunicação sem fio para estimar um número de pessoas reunidas em uma área de interesse Download PDF

Info

Publication number
BR112017020838B1
BR112017020838B1 BR112017020838-5A BR112017020838A BR112017020838B1 BR 112017020838 B1 BR112017020838 B1 BR 112017020838B1 BR 112017020838 A BR112017020838 A BR 112017020838A BR 112017020838 B1 BR112017020838 B1 BR 112017020838B1
Authority
BR
Brazil
Prior art keywords
interest
area
user equipment
radius
surface elements
Prior art date
Application number
BR112017020838-5A
Other languages
English (en)
Other versions
BR112017020838A2 (pt
Inventor
Massimo Colonna
Marco MAMEI
Original Assignee
Telecom Italia S.P.A.
Filing date
Publication date
Application filed by Telecom Italia S.P.A. filed Critical Telecom Italia S.P.A.
Priority claimed from PCT/EP2015/056878 external-priority patent/WO2016155769A1/en
Publication of BR112017020838A2 publication Critical patent/BR112017020838A2/pt
Publication of BR112017020838B1 publication Critical patent/BR112017020838B1/pt

Links

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04WWIRELESS COMMUNICATION NETWORKS
    • H04W4/00Services specially adapted for wireless communication networks; Facilities therefor
    • H04W4/02Services making use of location information
    • H04W4/021Services related to particular areas, e.g. point of interest [POI] services, venue services or geofences
    • H04W4/022Services related to particular areas, e.g. point of interest [POI] services, venue services or geofences with dynamic range variability
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04WWIRELESS COMMUNICATION NETWORKS
    • H04W4/00Services specially adapted for wireless communication networks; Facilities therefor
    • H04W4/02Services making use of location information
    • H04W4/029Location-based management or tracking services

Abstract

MÉTODO E SISTEMA PARA ESTIMAR UM NÚMERO DE PESSOAS REUNIDAS EM UMA ÁREA DE INTERESSE. Um método para estimar um número de pessoas reunidas em uma Área de Interesse (107) durante um intervalo de tempo em um dia, em que a dita Área de Interesse (107) é definida por um centro de Área de Interesse e um raio de Área de Interesse e é coberto por uma rede de telecomunicação móvel (105) tendo uma pluralidade de estações de comunicação (105a) cada uma das quais adaptada para gerenciar comunicações de equipamento de usuário em uma ou mais áreas atendidas (105b) em uma região geográfica coberta (300; 300') sobre a qual a rede de telecomunicação móvel (105) se estende, o método compreendendo as etapas de: a) subdividir a região geográfica coberta (300; 300') em uma pluralidade de elementos de superfície (205); b) definir (604) uma pluralidade de valores de raio calculados do raio de Área de Interesse, e, para cada valor de raio calculado: c) identificar (606) um número de elementos de superfície relevantes (505a-d) da região geográfica coberta (300; 300') compreendida dentro da Área de Interesse; d) computar (610; 612) um primeiro número de Equipamento de Usuário atendido pela rede de comunicação móvel (105)durante o intervalo de tempo no dia dentro da Área de Interesse (...).

Description

Fundamentos da Invenção Campo da Invenção
[001] A presente invenção refere-se à contagem de multidão, isto é, às técnicas para contagem ou estimativa do número de pessoas em uma multidão. Na presente descrição e com os propósitos da presente invenção, por “multidão” entende-se uma reunião de um certo número de pessoas, reunidas em um certo local para, por exemplo, estar presente em eventos ou acontecimentos públicos, de naturezas mais disparates, como, por exemplo, (e não exaustivamente) acontecimentos públicos de televisão ao vivo, performances artísticas/de entretenimento, exibições culturais, peças teatrais, competições esportivas, concertos, filmes, demonstrações e/ou para visitar um lugar de interesse em particular, tais como, por exemplo, um museu, um monumento, um prédio e congêneres.
[002] Particularmente, a presente invenção refere-se a técnicas de contagem de multidão que exploram a informação provida por redes de telecomunicação sem fio ou móveis.
Visão Geral da Técnica Relacionada
[003] Nas tarefas de planejamento urbano, de gerenciamento de atividades (por exemplo, gerenciamento de sistemas de transporte e gerenciamento de emergências) e de marketing para turismo e local, é útil ter um conhecimento de quantidades de pessoas que se reuniram em certos locais ou Áreas de Interesse (AoI em abreviatura, por exemplo, um prédio, tais como, por exemplo, um estádio ou um teatro ou uma sala de projeções, as cercanias dos mesmos, uma praça ou uma(s) rua(s) de uma cidade ou município ou vila, um distrito, etc.), por exemplo, em virtude de eles terem estado presentes em acontecimentos públicos, como espetáculos (por exemplo, relacionados a cultura, entretenimento, política ou esportes) que ocorreram na área de interesse e/ou para visitar um lugar de interesse (também denotado como ponto de interesse) na área de interesse.
[004] No caso de uma ou mais reuniões de pessoas relacionadas a acontecimentos públicos, embora as seguintes considerações também se apliquem a reuniões de pessoas relacionadas a pontos de interesse, este conhecimento permite, por exemplo, um planejamento mais efetivo de subsequentes acontecimentos públicos do mesmo tipo. Particularmente, este conhecimento permite um planejamento e um gerenciamento mais efetivos de recursos e atividades (tais como infraestruturas, sistema de transporte e segurança) diretamente ou indiretamente relacionados a acontecimentos públicos similares que podem ocorrer no futuro (tais como, por exemplo, jogos esportivos que ocorrem regularmente em um estádio). Além do mais, a partir de um ponto de vista comercial, este conhecimento permite um melhor gerenciamento de atividades de marketing projetadas para promover eventos similares que podem ocorrer no futuro.
[005] Hoje em dia, os dispositivos de comunicação móvel (referidos como telefones celulares ou UE a seguir, incluindo telefones celulares, telefones inteligentes, tablets e congêneres) têm alcançado uma completa difusão entre a população de muitos países, e possuidores de telefone celular quase sempre carregam seus telefones celulares consigo. Já que telefones celulares comunicam com uma pluralidade de estações bases das redes de telefonia celular, e cada estação base cobre (isto é, atende) uma ou mais áreas de atendimento predeterminadas, ou células, que são conhecidas pelo provedor de serviços de comunicação móvel (por exemplo, possuidor da rede de telefonia celular ou provedor de serviços do telefone celular virtual), os telefones celulares resultam como candidatos ideais como dispositivos de rastreamento para coletar dados usados para identificar a quantidade de pessoas que estiveram presentes em um ou mais acontecimentos públicos.
[006] Na técnica, muitos sistemas e métodos foram propostos a fim de coletar informação sobre horário e locais em que, e nos quais, um Equipamento de Usuário (UE, por exemplo, um telefone celular, um telefone inteligente, um tablet, etc.) de um indivíduo conecta na rede de telefonia celular (por exemplo, para realizar uma chamada de voz ou enviar uma mensagem de texto), e usa tal informação coletada a fim de derivar informação relacionada a quantos presentes um certo acontecimento público teve.
[007] Por exemplo, F. Calabrese, F. C. Pereira, G. Di Lorenzo, L. Liu, C. Ratti, "The Geography of Taste: Analyzing Cell-Phone Mobility in Social Events," Pervasive Computing, LNCS 6030, Springer, 2010, pp. 22-37, descreve a análise de mobilidade de multidão durante eventos especiais. Quase 1 milhão de rastros de telefone celular foram analisados e associados com seus destinos com eventos sociais. Foi observado que as origens de pessoas que estão presentes em um evento são fortemente correlacionadas ao tipo de evento, com implicações em gerenciamento da cidade, já que o conhecimento de fluxos aditivos pode ser uma informação crítica sobre a qual tomam-se decisões sobre gerenciamento de eventos e mitigação de congestionamento.
[008] Traag, V. A.; Browet, A.; Calabrese, F.; Morlot, F., "Social Event Detection in Massive Mobile Phone Data Using Probabilistic Location Inference", 2011 IEEE Third International Conference on Privacy, Security, Risk and Trust (Passat), and 2011 IEEE Third International Conference on Social Computing (Socialcom), pp.625, 628, 9-11 Oct. 2011, focam em reuniões de pessoas incomumente grandes, isto é, eventos sociais incomuns. A metodologia de detecção de tais eventos sociais em massivos dados de telefone celular é introduzida, com base em uma estrutura de inferência de local Bayesiana. Mais especificamente, uma estrutura para decidir quem está presente em um evento também é desenvolvida. O método, em poucos exemplos, é demonstrado. Finalmente, algumas possíveis futuras abordagens para detecção de evento, e algumas possíveis análises dos eventos sociais detectados são discutidas.
[009] Francesco Calabrese, Carlo Ratti, "Real Time Rome", Networks and Communications Studies 20(3-4), pages 247-258, 2006, descrevem o projeto Real Time Rome, apresentado na 10a Exibição de Arquitetura Internacional em Veneza, Itália. O projeto Real Time Rome é o primeiro exemplo de um sistema de monitoramento urbano em tempo real que coleta e processa dados providos por redes de telecomunicações e sistemas de transporte a fim de entender padrões da vida diária em Roma. Observar a vida diária em tempo real em um município torna-se um meio de entendimento do presente e previsão do ambiente urbano futuro.
[0010] F. Manfredini, P. Pucci, P. Secchi, P. Tagliolato, S. Vantini, V. Vitelli, "Treelet decomposition of mobile phone data for deriving city usage and mobility pattern in the Milan urban region", MOX-Report No. 25/2012, MOX, Departamento de Matemática "F. Brioschi", Politecnico di Milano, disponível em http://mox.polimi.it, descrevem uma técnica de aprendizado geoestatístico não supervisionado visada a identificar informação útil em padrões ocultos de uso de telefone celular. Estes padrões ocultos consideram diferentes usos da cidade no tempo e no espaço que são relacionados a mobilidade individual, esquematizando o potencial desta tecnologia para a comunidade de planejamento urbano. A metodologia permite obter uma tendência de referência que reporta o efeito específico de algumas atividades nos dados Erlang gravados e um conjunto de mapas que mostram a contribuição de cada atividade para o sinal Erlang local. Os resultados que são significativos para explicar mobilidade e padrões de usos da cidade específicos (deslocamentos para o trabalho, atividades noturnas, distribuição de residências, mobilidade não sistemática) foram selecionados e sua significância e sua interpretação a partir de uma perspectiva de análise e planejamento urbanos na escala da região urbana de Milão foram testadas.
[0011] Ramon Caceres, James Rowland, Christopher Small, and Simon Urbanek, "Exploring the Use of Urban Greenspace through Cellular Network Activity", 2nd Workshop on Pervasive Urban Applications (PURBA), June 2012, descrevem o uso de gravações anônimas de atividade da rede celular para estudar os padrões espaço-temporais da densidade humana em uma área urbana. Este documento apresenta a visão e alguns resultados anteriores deste esforço. Primeiramente, um conjunto de dados de seis meses de atividade na área metropolitana de Nova Iorque é descrito. Em segundo lugar, uma técnica para estimar áreas de cobertura da rede é apresentada. Em terceiro lugar, a abordagem usada na análise de mudanças nos volumes de atividade nestas áreas é descrita. Finalmente, os resultados preliminares em relação a mudanças na densidade humana ao redor do Central Park são apresentados.
Sumário da Invenção
[0012] O Requerente observou que, no geral, o método e os sistemas conhecidos na técnica proveem resultados insatisfatórios, já que eles não são capazes de determinar (ou têm uma capacidade limitada de determinação) se um possuidor de UE esteve em uma área de interesse (AoI) onde um ou mais acontecimentos públicos foram organizados, por estar presentes nos mesmos ou por outros motivos (por exemplo, em virtude de o possuidor de UE residir ou ter uma empresa na proximidade da, ou na, área de interesse). Além do mais, os resultados providos pelas soluções conhecidas são fortemente influenciados pelo tamanho da área de interesse selecionada para a análise da quantidade de presentes nos um ou mais acontecimentos públicos. Em outras palavras, se a área de interesse tiver um grande tamanho, um certo número de possuidores de UE que não são realmente parte da multidão será considerado na avaliação do número de pessoas na multidão. Inversamente, se a área de interesse tiver pequeno tamanho, um certo número de possuidores de UE realmente parte da multidão será excluído da avaliação do número de pessoas na multidão.
[0013] Portanto, subsequentes planejamento e gerenciamento de recursos e atividades (do tipo supramencionado) com base nos resultados obtidos pelos métodos e sistemas conhecidos na técnica irão atingir uma eficiência limitada devido à precisão limitada dos mesmos.
[0014] Além do mais, conhecidos métodos e sistemas com base no uso de informação em relação às posições ocupadas por cada UE enquanto conectado na rede de telefonia celular (informação que é coletada por redes de telefonia celular que atendem o UE) podem ser intrusivos de uma privacidade dos possuidores do UE.
[0015] De fato, tal informação permite conhecer hábitos, rotinas e locais (por exemplo, locais de casa e trabalho) diariamente frequentados pelos possuidores de UE.
[0016] Desta maneira, o uso de tal informação é, assim, usualmente, altamente restrito (até o limite de ser proibido) por leis emitidas por muitas Autoridades Nacionais a fim de proteger a privacidade dos possuidores de UE.
[0017] Neste aspecto, as técnicas de "anonimização" conhecidas na técnica e usadas para anonimizar a informação sobre os possuidores de UE, a fim de burlar questões de privacidade, não concedem uma proteção satisfatória da privacidade dos mesmos.
[0018] No geral, as técnicas de anonimização compreendem mascarar todos os identificadores (tais como, por exemplo, uma Identidade de Equipamento Móvel Internacional - IMEI, uma Identidade de Assinante Móvel Internacional - IMSI, ou um Número ISDN de Assinante Móvel) associados com o UE e/ou o possuidor de UE com identificadores cifrados.
[0019] Contudo, uma análise da informação coletada em uma pluralidade de dias pode ser intrusiva de possuidores de UE, já que ela, de qualquer modo, permite identificar informação sensível em relação a hábitos, locais de casa e trabalho dos possuidores de UE e, possivelmente, os próprios possuidores de UE pela análise de tal informação sensível assim obtida.
[0020] O Requerente, assim, lidou com o problema de conceber um sistema e método adaptados para superar os problemas que afetam as soluções da técnica anterior.
[0021] O Requerente descobriu que é possível determinar o tamanho de uma área de interesse ideal com base nos dados relacionados ao UE durante o curso dos um ou mais acontecimentos públicos e em um certo número de dias que precedem os um ou mais acontecimentos públicos.
[0022] Além do mais, o Requerente descobriu que é possível proteger a privacidade dos possuidores de UE pela exploração de dados agregados em relação ao uso de pelo menos uma rede de telefonia celular.
[0023] Por exemplo, os dados agregados exploráveis pela presente invenção compreendem um número de UEs atendidos pela rede de telefonia celular em um ou mais intervalos de tempo (isto é, nenhuma informação sobre UE individual é provida que pode infringir na privacidade do possuidor de UE).
[0024] Preferivelmente, pelo uso de dados agregados em relação separadamente a uma ou mais áreas atendidas da rede de telefonia celular, é possível determinar o tamanho de uma área de interesse ideal e, então, um número de pessoas que se reuniu na mesma com uma alta precisão.
[0025] Deve-se notar que o conhecimento do número de UEs atendidos pela rede de telefonia celular em um ou mais intervalos de tempo na área ideal, no geral, não corresponde ao número de pessoas na multidão. De fato, o número de UEs atendidos pela rede de telefonia celular em um ou mais intervalos de tempo na área ideal compreende os UEs possuídos por pessoas na área ideal por motivos (por exemplo, trabalho, pessoas simplesmente cruzando a área ideal) diferentes da reunião na multidão.
[0026] O Requerente descobriu adicionalmente que é possível discernir o número de pessoas na multidão na área ideal das pessoas que estão na área ideal, mas não estão na multidão com base na análise de dados agregados referidos ao uso da rede de telefonia celular durante a reunião e durante um número de dias anteriores ao dia no qual a reunião de pessoas ocorreu.
[0027] Particularmente, um aspecto da presente invenção propõe um método para estimar um número An de pessoas reunidas em uma área de interesse durante um intervalo de tempo em um dia, em que a dita área de interesse é definida por um centro da área de interesse C e um raio da Área de Interesse Rs e é coberta por uma rede de telecomunicação móvel que tem uma pluralidade de estações de comunicação, cada uma das quais adaptada para gerenciar as comunicações do equipamento de usuário em uma ou mais áreas atendidas em uma região geográfica coberta sobre a qual a rede de telecomunicação móvel provê serviços, o método compreendendo as etapas de: a) subdividir a região geográfica coberta em uma pluralidade de elementos de superfície; b) definir uma pluralidade de valores de raio calculados Rk do raio da Área de Interesse Rs e, para cada valor de raio calculado: c) identificar um número de elementos de superfície relevantes da região geográfica coberta compreendidos dentro da área de interesse; d) computar um primeiro número Uk de Equipamento de Usuário atendido pela rede de comunicação móvel durante o intervalo de tempo [Ts, Te] no dia g na área de interesse com base nos dados agregados uq,t em relação a um uso da rede de comunicação móvel; e) computar um segundo número Upk de Equipamento de Usuário atendido pela rede de comunicação móvel durante o intervalo de tempo [Ts, Te] para cada dia gp de um número predeterminado de dias anteriores P precedendo o dia g na área de interesse com base nos dados agregados uq,t em relação ao uso da rede de comunicação móvel; f) combinar o primeiro número Uk de Equipamento de Usuário e os segundos números Upk de Equipamento de Usuário para obter uma quantidade estatística Zk; g) computar uma quantidade estatística normalizada Z'k pela normalização da quantidade estatística Zk em relação aos raios Ra das áreas atendidas relevantes; h) computar um valor de raio ideal Ro do raio da Área de Interesse como a média dos valores de raio calculados pesados pela quantidade estatística normalizada Z'k; i) estimar o número An de pessoas reunidas na área de interesse que tem o raio da Área de Interesse Rs igual ao valor de raio ideal Ro.
[0028] Os recursos preferidos da presente invenção são apresentados nas reivindicações dependentes.
[0029] Em uma modalidade da invenção, os dados agregados uq,t em relação a um uso da rede de comunicação móvel compreendem um número de carga de tráfego do Equipamento de Usuário atendido, número de chamadas de voz, número de SMS transmitido e/ou volume de dados binários trocados em, preferivelmente, cada uma das estações de comunicação da rede de comunicação móvel.
[0030] Em uma modalidade da invenção, um elemento de superfície é identificado como um elemento de superfície relevante se ele verifica a seguinte desigualdade: Dist(C, B) < IRs + Rkl, onde C é o centro da área de interesse, B é o centro do elemento de superfície, Dist(C, B) é a distância geográfica entre o centro da área de interesse C e o centro do elemento de superfície B, Rs é o raio do elemento de superfície, e Rk é o valor de raio calculado.
[0031] Em uma modalidade da invenção, o método compreende adicionalmente, para cada valor de raio calculado: j) receber uma pluralidade de dados agregados uq,t em relação a um uso da rede de comunicação móvel referida a cada um dos ditos elementos de superfície.
[0032] Em uma modalidade da invenção, a etapa i) de receber uma pluralidade de dados agregados uq,t em relação a um uso da rede de comunicação móvel para cada um dos ditos elementos de superfície compreende receber um conjunto {uq,t} de dados agregados, cada um dos dados agregados uq,t do conjunto {uq,t} dos dados agregados sendo referido como um respectivo intervalo de tempo de referência dt que é uma porção de um período de aquisição ΔT durante o qual os dados agregados uq,t são coletados.
[0033] Em uma modalidade da invenção, a etapa c) de computar um primeiro número Uk de Equipamento de Usuário atendido pela rede de comunicação móvel durante o intervalo de tempo de observação [Ts, Te] no dia g na área de interesse com base nos dados agregados uq,t em relação a um uso da rede de comunicação móvel compreende computar o primeiro número Uk de Equipamento de Usuário com base nos conjuntos {uq,t} de dados agregados referidos aos respectivos intervalos de tempo de referência dt compreendidos dentro do intervalo de tempo de observação [Ts, Te] no dia g, e em que a etapa d) de computar um segundo número Upk de Equipamento de Usuário que foi atendido pela rede de comunicação móvel durante o intervalo de tempo de observação [Ts, Te] para cada dia gpn de um número predeterminado P de dias anteriores precedentes ao dia g na área de interesse com base nos dados agregados uq,t em relação ao uso da rede de comunicação móvel compreende computar cada segundo número Upk de Equipamento de Usuário com base nos conjuntos {uq,t} de dados agregados referidos aos respectivos intervalos de tempo de referência dt compreendidos dentro do intervalo de tempo de observação [Ts, Te] do respectivo dia anterior gpn do número predeterminado P de dias anteriores precedentes ao dia g.
[0034] Em uma modalidade da invenção, o primeiro número Uk de Equipamento de Usuário e/ou cada segundo número Upk de Equipamento de Usuário podem ser computados como um número total, um número médio ou um número máximo de Equipamento de Usuário nos elementos de superfície relevantes compreendidos dentro da área de interesse.
[0035] Em uma modalidade da invenção, a etapa e) de combinar o primeiro número Uk de Equipamento de Usuário e os segundos números Upk de Equipamento de Usuário para obter uma quantidade estatística Zk compreende: combinar o segundo número Upk de Equipamento de Usuário de cada um dos dias anteriores gpn a fim de determinar um número médio de equipamentos de usuário k e um desvio padrão do número de equipamentos de usuário ok.
[0036] Em uma modalidade da invenção, a etapa e) de combinar o primeiro número Uk de Equipamento de Usuário e os segundos números Upk de Equipamento de Usuário para obter uma quantidade estatística compreende adicionalmente: computar a quantidade estatística como: Zk = (Uk - μk) / ok, onde Uk é o primeiro número de Equipamento de Usuário, μk é o número médio de equipamentos de usuário e ok é o desvio padrão do número de equipamentos de usuário.
[0037] Em uma modalidade da invenção, a pluralidade de valores de raio calculados Rk varia de um valor de raio mínimo Rmin até um valor de raio máximo Rmax, cada valor de raio calculado sendo separado de um próximo valor de raio por uma largura de iteração Δ.
[0038] Em uma modalidade da invenção, a etapa h) de estimar um número An de pessoas reunidas na área de interesse que tem o raio da Área de Interesse Rs igual ao valor de raio ideal Ro compreende: k) definir um número de elementos de superfície relevantes dentre os elementos de superfície compreendidos na região geográfica coberta, em que os ditos elementos de superfície relevantes são elementos de superfície pelo menos parcialmente sobrepostos na Área de Interesse que tem o raio da Área de Interesse igual ao valor de raio ideal.
[0039] Em uma modalidade da invenção, um elemento de superfície é identificado como um elemento de superfície relevante se ele verifica a seguinte desigualdade: Dist(C, B) < IRs + Rol, onde C é o centro da área de interesse, B é o centro do elemento de superfície, Dist(C, B) é a distância geográfica entre o centro da área de interesse C e o centro do elemento de superfície B, Rs é o raio do elemento de superfície, e Ro é o valor de raio ideal.
[0040] Em uma modalidade da invenção, a etapa h) de estimar um número An de pessoas reunidas na área de interesse que tem o raio da Área de Interesse Rs igual ao valor de raio ideal Ro compreende adicionalmente: computar um terceiro número UIAOI de Equipamento de Usuário como um número de equipamentos de usuário compreendidos dentro dos elementos de superfície relevantes compreendidos dentro da área de interesse que tem o raio da Área de Interesse Rs igual ao valor de raio ideal Ro durante o intervalo de tempo com base nos dados agregados uq,t em relação a um uso da rede de comunicação móvel.
[0041] Em uma modalidade da invenção, a computação de um terceiro número UIAOI de Equipamento de Usuário compreende computar o terceiro número UIAOI de Equipamento de Usuário com base nos conjuntos {uq,t} de dados agregados referidos aos respectivos intervalos de tempo de referência dt compreendidos dentro do intervalo de tempo de observação [Ts, Te] no dia g.
[0042] Em uma modalidade da invenção, a etapa h) de estimar um número An de pessoas reunidas na área de interesse que tem o raio da Área de Interesse Rs igual ao valor de raio ideal Ro compreende adicionalmente: computar um quarto número UpIAOI de Equipamento de Usuário como um número de equipamentos de usuário compreendidos dentro dos elementos de superfície relevantes compreendidos dentro da área de interesse que tem o raio da Área de Interesse Rs igual ao valor de raio ideal Ro para cada dia gpn do número predeterminado P de dias anteriores precedentes ao dia com base nos dados agregados uq,t em relação a um uso da rede de comunicação móvel.
[0043] Em uma modalidade da invenção, a computação de um quarto número Up|AOI de Equipamento de Usuário compreende computar cada quarto número Up|AOI de Equipamento de Usuário com base nos conjuntos {uq,t} de dados agregados referidos aos respectivos intervalos de tempo de referência dt compreendidos dentro do intervalo de tempo de observação [Ts, Te] do respectivo dia anterior gpn do número predeterminado P de dias anteriores precedentes ao dia g.
[0044] Em uma modalidade da invenção, o terceiro número U|AOI de Equipamento de Usuário e/ou cada quarto número Up|AOI de Equipamento de Usuário podem ser computados como um número total, um número médio ou um número máximo de Equipamento de Usuário nos elementos de superfície relevantes compreendidos dentro da área de interesse que tem o raio da Área de Interesse igual ao valor de raio ideal.
[0045] Em uma modalidade da invenção, a etapa h) de estimar um número An de pessoas reunidas na área de interesse que tem o raio da Área de Interesse Rs igual ao valor de raio ideal Ro compreende adicionalmente: combinar o quarto número Up|AOI de Equipamentos de Usuário de cada um dos dias anteriores gpn a fim de determinar um número de equipamentos de usuário médio adicional μlAOI, o número de equipamentos de usuário médio adicional μlAOI provendo uma indicação de um número médio de pessoas normalmente compreendidas dentro da área de interesse que tem o raio da Área de Interesse Rs igual ao valor de raio ideal Ro durante o intervalo de tempo de observação considerado [Ts, Te] em qualquer dia.
[0046] Em uma modalidade da invenção, a etapa h) de estimar um número An de pessoas reunidas na área de interesse que tem o raio da Área de Interesse Rs igual ao valor de raio ideal Ro compreende adicionalmente: combinar o terceiro número U|AOI de Equipamento de Usuário e o número de equipamentos de usuário médio adicional μlAOI a fim de obter o número An de pessoas reunidas na área de interesse que tem o raio da Área de Interesse Rs igual ao valor de raio ideal Ro.
[0047] Em uma modalidade da invenção, a combinação do terceiro número U|AOI de Equipamento de Usuário e do número de equipamentos de usuário médio adicional μ|AOI compreende subtrair o número de equipamentos de usuário médio adicional μ|AOI do terceiro número UPIAOI de Equipamento de Usuário.
[0048] Um outro aspecto da presente invenção propõe um sistema acoplado com uma rede de telecomunicação sem fio para estimar um número de pessoas reunidas em uma área de interesse, o sistema compreendendo: um motor de computação adaptado para processar dados recuperados a partir de uma rede de telefonia móvel; um repositório adaptado para armazenar dados relativos a interações entre o Equipamento de Usuário e a rede de telefonia móvel, resultados de computação gerados pelo motor de computação e, possivelmente, quaisquer dados de processamento gerados por e/ou providos para o sistema, e uma interface de administrador operável para modificar parâmetros e/ou algoritmos usados pelo motor de computação e/ou acessar dados armazenados no repositório. Além do mais, o sistema compreende adicionalmente um elemento de memória que armazena um produto de programa de software configurado para implementar o método exposto através do sistema.
[0049] Em uma modalidade da invenção, o sistema compreende adicionalmente pelo menos uma interface de usuário adaptada para receber entradas a partir de, e para prover saída para, um usuário do sistema, o usuário compreendendo um ou mais seres humanos e/ou um ou mais sistemas de computação externos assinantes dos serviços providos pelo sistema.
[0050] Uma das vantagens da solução de acordo com a presente invenção é que ela é computacionalmente simples, envolvendo somente operações de contagem e operações algébricas.
Breve Descrição dos Desenhos
[0051] Estes e ainda outros recursos e vantagens da solução de acordo com a presente invenção serão melhor entendidos pela leitura da seguinte descrição detalhada de uma modalidade da mesma, provida meramente a título de exemplos não limitantes, a ser lida em conjunto com os desenhos anexos, em que: a figura 1 é uma representação esquemática de um sistema de estimativa de multidão de acordo com uma modalidade da presente invenção; as figuras 2A - 2E são formas exemplares que os elementos de superfície podem tomar de acordo com uma modalidade da presente invenção; as figuras 3A - 3B são exemplos de regiões geográficas cobertas associadas com uma rede de comunicação móvel subdividida em correspondentes conjuntos de elementos de superfície de acordo com uma modalidade da presente invenção; as figuras 4A - 4E são formas exemplares que a AoI a ser determinada pode tomar de acordo com uma modalidade da presente invenção; as figuras 5A - 5D são elementos de superfície relevantes dentre elementos de superfície nos quais a região geográfica coberta AoI é subdividida de acordo com uma modalidade da invenção, e as figuras 6A - 6C são um fluxograma esquemático de um algoritmo de estimativa de multidão de acordo com uma modalidade da presente invenção.
Descrição Detalhada da Invenção
[0052] Em relação aos desenhos, a figura 1 é uma representação esquemática de um sistema de estimativa de multidão, simplesmente denotado como sistema 100 a seguir, de acordo com uma modalidade exemplar da presente invenção.
[0053] O sistema de estimativa de multidão e o método descritos a seguir permitem realizar uma estimativa de um número de pessoas em uma multidão reunida, por exemplo, a fim de estar presente em um ou mais acontecimentos públicos, da natureza mais disparate, como, por exemplo, (e não exaustivamente) acontecimentos públicos de televisão ao vivo, performances artísticas/de entretenimento, exibições culturais, peças teatrais, competições esportivas, concertos, filmes, demonstrações e congêneres.
[0054] Além do mais, entende-se que o sistema e o método de estimativa de multidão descritos a seguir também permitem realizar uma estimativa de um número de pessoas em uma multidão reunida para visitar um lugar de interesse em particular, tais como, por exemplo, um museu, um monumento, um prédio histórico e congêneres.
[0055] O sistema 100 é acoplado em uma rede de comunicação móvel 105, tal como uma rede de telefonia móvel (2G, 3G, 4G ou geração superior).
[0056] A rede de comunicação móvel 105 pode prover recursos de comunicação (por exemplo, uma porção de uma largura de banda de comunicação disponível) para o Equipamento de Usuário, UE a seguir (por exemplo, um telefone celular, um telefone inteligente, um tablet com conectividade 2G-3G-4G, etc.) que solicita os mesmos em uma região geográfica coberta (não detalhada na figura 1, mas descrita a seguir em relação às figuras 3A e 3B). Em outras palavras, o UE na região geográfica coberta pode ser atendido pela rede de comunicação móvel.
[0057] Em uma modalidade da presente invenção, a região geográfica coberta pode compreender a íntegra de um território coberto (isto é, atendido) pela rede de comunicação móvel 105, mesmo embora, em outras modalidades da presente invenção, uma região geográfica coberta que compreende apenas uma porção da íntegra do território coberto pela rede de comunicação móvel 105 possa ser considerada.
[0058] A rede de comunicação móvel 105 compreende uma pluralidade (isto é, duas ou mais) de estações de comunicação 105a (por exemplo, estações bases de rádio da rede de telefonia móvel) implementadas na região geográfica coberta.
[0059] Cada estação de comunicação 105a é adaptada para gerenciar as comunicações do UE (não mostrado, tais como, por exemplo, telefones celulares) em uma ou mais áreas ou células atendidas 105b (no exemplo em questão, três células são atendidas por cada estação de comunicação 105a).
[0060] Desta maneira, a região geográfica coberta compreende a área de uma pluralidade das células 105b da rede de comunicação móvel 105, por exemplo, em uma modalidade da invenção, a soma das áreas de todas as células 105b da rede de comunicação móvel 105 constitui a região geográfica coberta.
[0061] Diferentemente, uma Área de Interesse (geográfica), AoI em resumo, 107 esquematizada na figura 1 como a área na linha de traço e ponto 107 se estende sobre uma ou mais células 105b da rede de comunicação móvel 105. A AoI 107 é uma área na qual as pessoas reunidas em uma multidão, por exemplo, a fim de estar presentes em um ou mais acontecimentos públicos e cuja extensão é determinada pelo algoritmo de estimativa de multidão da presente invenção (da forma descrita a seguir).
[0062] O sistema 100 é configurado para receber a partir da rede de comunicação móvel 105 dados agregados em relação ao uso da rede de comunicação móvel 105 em um ou mais intervalos de tempo de referência (da forma descrita a seguir).
[0063] O termo 'dados agregados', da forma usada na presente descrição, indica dados relativos à operação da rede de comunicação móvel 105, tais como, por exemplo, um número de UEs atendidos, a carga de tráfego, o número de chamadas de voz, p número de SMS transmitido, o volume de dados binários trocados, etc. Os dados agregados são tipicamente usados por um operador que gerencia a rede de comunicação móvel 105 para analisar tendências ou valores gerais (por exemplo, mudanças de um número de UEs atendidos durante o tempo) na exploração de recursos (por exemplo, largura de banda e/ou capacidades computacionais) da rede de comunicação móvel 105 sem identificar cada UE individual (e, portanto, os possuidores do UE) que interagiu com a rede de comunicação móvel 105.
[0064] Particularmente, os dados agregados não compreendem nenhum identificador do UE atendido pela rede de comunicação móvel 105, portanto, nenhuma indicação sobre as identidades, hábitos ou locais frequentados dos possuidores de UE (tais como, por exemplo, locais de casa e trabalho) pode ser obtida a partir dos dados agregados providos pela rede de comunicação móvel 105, assim, a privacidade dos possuidores de UE é garantida.
[0065] Em uma modalidade da invenção, os dados agregados adquiridos pelo sistema 100 da rede de comunicação móvel 105 compreendem uma indicação em relação a um número de UE, isto é, indicativa, também, de um número de indivíduos (os possuidores de UE), localizados na região geográfica coberta.
[0066] A região geográfica coberta recentemente descrita compreende a área de interesse, AoI em resumo, esquematizada na figura 1 como a área na linha de traço e ponto 107.
[0067] A AoI 107 (adicionalmente descrita a seguir) pode, no geral, compreender um local central (por exemplo, um estádio, um teatro, uma praça da cidade e congêneres) onde um acontecimento público, que atrai uma respectiva multidão, está ocorrendo e, possivelmente, as cercanias (por exemplo, estacionamentos nas proximidades, ruas nas proximidades, estações de transporte nas proximidades e congêneres) do local central.
[0068] Deve-se notar que, já que a AoI 107 está compreendida na região geográfica coberta da rede de comunicação móvel 105, assim, o UE na AoI 107 pode ser atendido pela rede de comunicação móvel 105.
[0069] Preferivelmente, os dados agregados recebidos a partir da rede de comunicação móvel 105 pelo sistema 100 compreendem uma indicação de um número de UEs atendidos pela rede de comunicação móvel 105.
[0070] Mais preferivelmente, os dados agregados recebidos a partir da rede de comunicação móvel 105 pelo sistema 100 compreendem as indicações do número do UE atendido pela rede de comunicação 105 em uma pluralidade de subpartes, ou elementos de superfície (descritos a seguir), da região geográfica coberta. Por exemplo, cada elemento de superfície pode compreender uma das células 105b, um grupo de duas ou mais células 105b, e/ou partes das células 105b da rede de comunicação móvel 105 (por exemplo, em uma modalidade da invenção, os elementos de superfície podem ser modelados como quadrados que têm um lado de 150 m, portanto, cada uma das células 105b compreende mais do que um elemento de superfície especialmente em áreas extraurbanas onde as células 105b usualmente têm uma maior extensão em relação às células 105b nas áreas urbanas).
[0071] No geral, cada estação de comunicação 105a da rede de comunicação móvel 105 é adaptada para interagir com qualquer UE localizado em uma das células 105b atendidas por tal estação de comunicação 105a (por exemplo, interações no ligamento/desligamento, na atualização da área local, em chamadas de entrada/saída, no envio/recepção de SMS e/ou de MMS, no acesso à Internet, etc.). Tais interações entre o UE e a rede de comunicação móvel 105 serão, no geral, denotadas como eventos ei (i = 1,... , I; em que I é um número inteiro) a seguir.
[0072] Portanto, os dados agregados que compreendem a indicação de um número de UEs podem ser computados simplesmente pela contagem de um número de UEs que teve pelo menos uma interação com a rede de comunicação móvel 105. Em outras palavras, os dados agregados compreendem uma indicação do número de UEs que produziu pelo menos um evento ei com uma das estações de comunicação 105a que provê serviços sobre respectivas células 105b da rede de comunicação móvel 105.
[0073] Alternativamente, a indicação de um número de UEs pode ser com base na carga de tráfego experimentada pela rede de comunicação móvel 105. De fato, cada evento ei, a fim de ser realizado, exige uma porção de recursos de comunicação (por exemplo, partes de uma largura de banda de comunicação) gerenciados pela rede de comunicação móvel 105, isto é, cada evento ei produz uma certa quantidade de carga de tráfego. Desta maneira, os dados agregados, preferivelmente, proveem uma estimativa de um número de UEs com base em uma carga de tráfego dividida por uma carga de tráfego média do UE (isto é, uma carga de tráfego média gerada por cada UE associado com a rede de comunicação móvel 105).
[0074] Na presente descrição, considera-se que os dados agregados são providos pela rede de comunicação móvel 105 periodicamente, por exemplo, no lapso de intervalos de tempo de referência predeterminados (por exemplo, a cada certo número de minutos, horas, em uma base diária ou semanal) de acordo com uma capacidade da rede de comunicação móvel 105 de coletar, processar e prover tais dados agregados.
[0075] Por exemplo, em uma modalidade da presente invenção, os dados agregados são coletados pela rede de comunicação móvel 105 com uma periodicidade igual a, ou inferior a, quinze (15) minutos (isto é, os intervalos de tempo de referência têm uma duração de quinze minutos cada), que é uma periodicidade sustentável por capacidades computacionais das presentes redes de comunicação móveis.
[0076] Contudo, a rede de comunicação móvel 105 pode ser configurada para prover dados agregados de forma assíncrona ou, alternativamente, os dados agregados podem ser providos por um terceiro módulo de processamento (não mostrado) associado com a rede de comunicação móvel 105 para receber os dados relativos à operação da rede de comunicação móvel 105 e com o sistema 100 para prover os dados agregados sem fugir do escopo da presente invenção.
[0077] O sistema 100 compreende um motor de computação 110 configurado para processar dados agregados recuperados da rede de comunicação móvel 105, e um repositório 115 (tais como uma base de dados, um sistema de arquivo, etc.) configurado para armazenar dados agregados recebidos a partir da rede de comunicação móvel 105, resultados de computação gerados pelo motor de computação 110 e, possivelmente, quaisquer dados de processamento gerados por e/ou providos para o sistema 100 (no geral, em um formato binário). O sistema 100 é provido com uma interface de administrador 120 (por exemplo, um computador) configurada e operável para modificar parâmetros e/ou algoritmos usados pelo motor de computação 110 e/ou para acessar dados armazenados no repositório 115.
[0078] Preferivelmente, o sistema 100 compreende uma ou mais interfaces de usuário 125 (por exemplo, um terminal de usuário, um software em execução em um terminal remoto conectado no sistema 100) adaptadas para receber entradas a partir de, e para prover saídas para, um usuário do sistema 100. O termo "usuário do sistema" da forma usada na presente descrição pode se referir a um ou mais seres humanos e/ou a sistemas de computação externos (tal como uma rede de computador, não mostrada) de um terceiro que é assinante de serviços providos pelo sistema 100 e habilitado para acessar o sistema 100 - por exemplo, sob assinatura de um contrato com um provedor de serviço possuidor do sistema 100 e, tipicamente, com reduzido direito de acesso ao sistema 100, se comparado com o direito de acesso mantido por um administrador do sistema 100 que opera através da interface de administrador 120.
[0079] Deve ser percebido que o sistema 100 pode ser implementado de qualquer maneira conhecida; por exemplo, o sistema 100 pode compreender um único computador, ou uma rede de computadores distribuídos, tanto do tipo físico (por exemplo, com uma ou mais máquinas principais que implementam o motor de computação 110 e o repositório 115 conectado em outras máquinas que implementam o administrador e as interfaces de usuário 120 e 125) ou de tipo virtual (por exemplo, pela implementação de uma ou mais máquinas virtuais em uma rede de computador).
[0080] Preferivelmente, o motor de computação 110 processa dados agregados (da forma discutida com detalhes a seguir) de acordo com um algoritmo de estimativa de multidão (descrito a seguir) implementado por um produto de programa de software armazenado em um elemento de memória 110a do sistema 100, por exemplo, compreendido no motor de computação 110 no exemplo da figura 1, mesmo embora o produto de programa de software também possa ser armazenado no repositório 115 (ou em qualquer outro elemento de memória provido no sistema 100).
[0081] Em operação, os dados agregados podem ser continuamente recuperados pelo sistema 100 da rede de comunicação móvel 105. Por exemplo, os dados agregados podem ser transferidos da rede de comunicação móvel 105 para o sistema 100 à medida que eles são gerados, em uma espécie de modalidade tipo "empurrar".
[0082] Alternativamente, os dados agregados podem ser coletados pela rede de comunicação móvel 105 com uma primeira periodicidade (por exemplo, cada 15 minutos) e, então, empacotados e transferidos para o sistema 100 com uma segunda periodicidade inferior à primeira periodicidade (por exemplo, cada hora), ou apenas mediante solicitação pelo sistema 100.
[0083] Os dados agregados recuperados da rede de comunicação móvel 105 são preferivelmente armazenados no repositório 115, onde eles são feitos disponíveis para o motor de computação 110 para processamento.
[0084] Os dados agregados são processados de acordo com as instruções providas pelo administrador do sistema (através da interface de administrador 120), por exemplo, armazenadas no repositório 115 e, possivelmente, de acordo com as instruções providas por um usuário (através da interface de usuário 125), da forma descrita a seguir. Finalmente, o motor de computação 110 provê os resultados do processamento realizado nos dados agregados para o usuário através da interface de usuário 125 e, opcionalmente, armazena tais resultados de processamento no repositório 115.
[0085] Deve ser adicionalmente notado que o método descrito na presente descrição pode ser implementado pelo uso de qualquer fonte de dados (por exemplo, providos por um ou mais dentre redes WiFi, WiMax, Bluetooth ou combinações das mesmas com as redes de telefonia móvel) a partir do que, é possível obter dados agregados que podem ser relacionados a um número de pessoas em uma área predeterminada (por exemplo, os elementos de superfície, ou a AoI 107).
[0086] Voltando, agora, para as figuras 2A - 2E, elas são formas exemplares nas quais os elementos de superfície 205 da região geográfica coberta associada com a rede de comunicação móvel 105 (isto é, a região geográfica coberta na qual a rede de comunicação móvel 105 pode atender o UE) podem ser modelados de acordo com uma modalidade da presente invenção.
[0087] Com os propósitos da presente invenção, cada elemento de superfície 205 da região geográfica coberta pela rede de comunicação móvel 105 pode ser modelado como uma superfície (da forma mostrada na figura 2A) que tem um respectivo centro de superfície B e um respectivo raio de superfície Rs.
[0088] De acordo com uma modalidade da presente invenção, no geral, o centro de superfície B e o raio de superfície Rs do elemento de superfície 205 não são relacionados com uma posição geográfica das uma ou mais estações de comunicação 105a ou as posições das células 105b da rede de comunicação móvel 105.
[0089] Da forma previamente notada, os elementos de superfície 205 podem se estender sobre uma ou mais células 105b da rede de comunicação móvel 105 ou, inversamente, os elementos de superfície 205 podem ser menores do que uma célula 105b da rede de comunicação móvel 105.
[0090] Deve-se notar que os elementos de superfície 205 não são limitados a uma forma tipo disco, de fato, os elementos de superfície 205 podem ter a forma de um polígono, preferivelmente, embora não estritamente necessariamente, regular. Neste caso, o centro de superfície B corresponde a um centro de massa (ou centroide) do polígono, ao mesmo tempo em que o raio de superfície Rs corresponde a um segmento adjacente ao centro de massa do polígono, isto é, o centro de superfície B, com um vértice do polígono (da forma mostrada nas figuras 2B e 2D) ou com um ponto médio de um lado do polígono (da forma mostrada nas figuras 2C e 2E).
[0091] Alternativamente, a rede de comunicação móvel 105 pode ser modelada por meio de um diagrama de infusão de mosaicos de Voronoi, em que cada célula de Voronoi corresponde a um elemento de superfície 205 da região geográfica coberta associada com a rede de comunicação móvel 105 (já que diagramas de infusão de mosaicos de Voronoi são bem conhecidos na técnica, eles não são aqui adicionalmente discutidos).
[0092] Preferivelmente, a modelagem, a lista e o número de elementos de superfície 205 da região geográfica coberta associada com a rede de comunicação móvel 105 podem ser predeterminados pela rede de comunicação móvel 105, pelo sistema 100, ou são inseridos no sistema 100 pelo administrador através da interface de administrador 120.
[0093] Considerando as figuras 3A e 3B, elas são exemplos de regiões geográficas cobertas 300 e 300' associadas com a rede de comunicação móvel 105 subdivididas em correspondentes conjuntos de elementos de superfície 2051-9 e 205’1-9, respectivamente, de acordo com uma modalidade da presente invenção.
[0094] A região geográfica coberta 300 mostrada na figura 3A foi subdividida em nove elementos de superfície 2051-9 que têm a forma de um polígono regular, isto é, um quadrado. Inversamente, a região geográfica 300' mostrada na figura 3B foi subdividida em nove elementos de superfície 205’19 que têm a forma de polígonos irregulares.
[0095] No geral, os recursos geométricos dos elementos de superfície 2051-9 e 2051-9 podem ser com base em um número de parâmetros da região geográfica, tais como, por exemplo, recursos urbanos (presença e distribuição de ruas, seções administrativas, etc.) e/ou recursos naturais (presença e distribuições de rios, colinas, etc.). Além do mais, outros sistemas de referências e/ou de mapeamento (tais como, por exemplo, bem conhecidas ferramentas de software de planejamento de rede usadas pelo provedor da rede de comunicação móvel 105) podem ser considerados para definir as formas e os tamanhos dos elementos de superfície 2051-9 e 205’1-9, além dos, ou como uma alternativa aos, recursos urbanos e naturais.
[0096] Deve-se notar que nada impede a definição de formas e tamanhos dos elementos de superfície 2051-9 e 205’1-9 de acordo com uma distribuição das células 105b da rede de comunicação móvel 105.
[0097] Além do mais, no geral, não há relacionamentos entre número, formas e tamanhos dos elementos de superfícies 2051-9 e 205’1-9 e a AoI 107.
[0098] Em uma modalidade da presente invenção, os elementos de superfície quadrados são preferivelmente usados, tais como os elementos de superfície 2051-9 da região geográfica coberta 300. Ainda mais preferivelmente, os elementos de superfície 2051-9 correspondem a pixels determinados durante o planejamento de rede pelas ferramentas de software de planejamento de rede supramencionadas.
[0099] De fato, os elementos de superfície quadrados 2051-9 permitem a simples subdivisão da região geográfica coberta 300, por exemplo, tendo determinado um ponto de referência (tal como, por exemplo, o centro de superfície B) de um elemento de superfície genérico 205q (q = 1, Q; em que Q é um número inteiro positivo, Q = 9 no exemplo da figura 3A) e o tamanho dos lados dos elementos de superfície quadrados 2051-9, é simplesmente possível determinar os vértices e os centros de superfície B de todos os elementos de superfície 2051-9.
[00100] Por exemplo, os pixels usados como elementos de superfície podem ser modelados como quadrados que têm um lado que tem um tamanho compreendido entre 200 m e 50 m, tal como 150 m. Isto garante uma boa proporcionalidade entre o nível de detalhes da região geográfica coberta 300 e a complexidade computacional exigida para analisar os dados agregados referidos à região geográfica coberta 300.
[00101] De acordo com uma modalidade da presente invenção, os dados agregados providos pela rede de comunicação móvel 105 compreendem uma indicação em relação a um número de UEs (e, portanto, de possuidores de UE) para cada um dos elementos de superfícies 2051-9 e 205’1-9.
[00102] A seguir, referência é feita apenas à área geográfica de cobertura 300 e aos respectivos elementos de superfície 2051-9 da figura 3A a título de simplicidade e concisão; entretanto, deve-se notar que considerações similares podem ser aplicadas na área geográfica de cobertura 300' e nos respectivos elementos de superfície 205’1-9 da figura 3B também.
[00103] Em uma modalidade da invenção, para cada elemento de superfície genérico 205q, o sistema 100 recebe a partir da rede de comunicação móvel 105 um respectivo número de UEs agregados uq.
[00104] Preferivelmente, o sistema 100 recebe da rede de comunicação móvel 105 uma pluralidade de números de UEs agregados uq,t para cada elemento de superfície genérico 205q. Cada número de UEs uq,t dos números de UEs agregados uq,t é referido a um t-ésimo intervalo de tempo de referência de uma pluralidade de consecutivos intervalos de tempo de referência dt (t = 1, ..., T; em que T é um número inteiro positivo).
[00105] Em outras palavras, o sistema 100 recebe um conjunto {uq,t} de Q números de UEs agregados uq,t (um para cada um dos elementos de superfície 205q; isto é, nove números de UEs agregados uq,t no exemplo da figura 3A), cada conjunto sendo referido como um intervalo de tempo de referência dt de T intervalos de tempo de referência consecutivos dt, por exemplo, cada conjunto {uq,t} de Q números de UEs agregados uq,t é gerado em instantes de tempo tt correspondentes ao final de um respectivo intervalo de tempo dt.
[00106] Por exemplo, pela consideração de um período de aquisição ΔT de vinte e quatro horas de duração (ΔT = 24 hr), durante o qual T conjuntos {uq,t} de Q números de UEs agregados uq,t são recebidos pelo sistema 100 para armazenamento e/ou processamento (da forma descrita a seguir), e intervalo de tempo dt de quinze minutos de duração (dt = 15 min.), no final do período de aquisição ΔT, o sistema 100 tem disponível T = 96 conjuntos {uq,t} de Q = 9 números de UEs agregados uq,t, um para cada intervalo de tempo de referência dt. De fato, o período de aquisição ΔT é subdividido em 96 consecutivos intervalos de tempo de referência dt que têm a seguinte estrutura: d1 = [00:00, 00:15), d2 = [00:15, 00:30), d95 = [23:30, 23:45) e d96 = [23:45, 00:00).
[00107] De acordo com uma modalidade da presente invenção, os números de UEs agregados uq,t são computados com base na carga de tráfego de cada célula 105b da rede de comunicação móvel 105 durante o correspondente intervalo de tempo de referência dt.
[00108] Preferivelmente, o número de UEs agregados uq,t é computado pela combinação da carga de tráfego (por exemplo, em Erlang) medida em cada uma das células 105b compreendidas dentro da área geográfica de cobertura 300 durante o correspondente t-ésimo intervalo de tempo de referência dt com a carga de tráfego média do UE estimada para o UE nas células 105b (isto é, uma carga de tráfego média gerada por cada UE associado com a rede de comunicação móvel 105).
[00109] A carga de tráfego de cada célula 105b é dividida pela carga de tráfego média do UE, assim, obtendo uma estimativa do número de UEs atendidos por cada célula 105b durante o intervalo de tempo de referência dt. Subsequentemente, o número de UEs atendidos por células 105b é distribuído entre os elementos de superfície 205q da região geográfica coberta 300.
[00110] Por exemplo, a determinação (isto é, a distribuição) do número de UEs em cada um dos elementos de superfície 205q da região geográfica coberta 300 pode ser com base no método descrito no documento de Francesco Calabrese, Carlo Ratti: "Real Time Rome", Networks and Communications Studies 20(3-4), pages 247-258, 2006 supramencionado com base no número de UEs atendidos por células 105b.
[00111] É apontado que a presente invenção é independente da quantidade usada para determinar o número de UEs agregados uq,t. De fato, à parte da carga de tráfego medida em Erlangs, outras formas de contagem da carga de tráfego podem ser usadas, por exemplo, um número de chamadas por célula 105b, um número de conexões por célula 105b ou o número de UEs exclusivos conectados por célula 105b (isto é, cada qual referido ao correspondente t-ésimo intervalo de tempo de referência dt).
[00112] Também é apontado que a presente invenção é independente de uma função usada para distribuir os números de UEs atendidos por células 105b entre os elementos de superfície 205q da região geográfica coberta 300.
[00113] Em uma outra modalidade da presente invenção, a rede de comunicação móvel 105 tem o conhecimento do número de UEs conectados em cada uma das células 105b, portanto, não há necessidade de determinar o número de UEs através da carga de tráfego, e os números de UEs agregados uq,t podem ser diretamente determinados pela adição em conjunto do número de UEs atendidos pelas células 105b compreendidas no elemento de superfície 205q durante o respectivo intervalo de tempo de referência dt.
[00114] Voltando, agora, para as figuras 4A - 4E, elas são formas exemplares que a AoI 107 a ser determinada pode tomar de acordo com uma modalidade da presente invenção.
[00115] No geral, a AoI 107 do acontecimento público pode ser modelada como uma área que tem um centro da AoI C e um raio da AoI Ra. Por exemplo, a AoI 107 pode ser delimitada por uma circunferência centralizada no centro da AoI C e que tem o raio da AoI Ra como o raio da circunferência (da forma mostrada na figura 4A).
[00116] Deve-se notar que a AoI 107 pode ter formas diferentes da circunferência. Por exemplo, a AoI 107 pode ter a forma de um polígono, preferivelmente, regular. Neste caso, o centro da AoI C corresponde a um centro de massa (ou centroide) do polígono, ao mesmo tempo em que o raio da AoI Ra corresponde a um segmento adjacente ao centro de massa do polígono com um vértice do polígono (da forma mostrada nas figuras 4B e 4D) ou com um ponto médio de um lado do polígono (da forma mostrada nas figuras 3C e 3E) de uma maneira similar aos elementos de superfície 205, da forma discutida anteriormente.
[00117] O centro da AoI C pode ser definido (por exemplo, por um usuário através da interface de usuário 125 ou por um administrador do sistema através da interface de administrador 120) como um ponto central (geográfico) da AoI 107 (por exemplo, um ponto central geográfico do local central), como um endereço do local central do acontecimento público, como um ponto provido por um software de mapeamento, tais como serviços de mapeamento na Internet (por exemplo, Google maps™, OpenStreetMap™, etc.).
[00118] Como será descrito com mais detalhes a seguir, o raio da AoI Ra pode tomar valores zero ou negativos juntamente com valores positivos. No caso em que o raio da AoI Ra tomar valores zero ou negativos, a AoI 107 é limitada ao centro da AoI C (isto é, o local central dos um ou mais acontecimentos públicos). O significado de valores zero ou negativos para o raio da AoI Ra será adicionalmente esclarecido pela referência a tais valores zero ou negativos nas modalidades descritas a seguir.
[00119] O algoritmo descrito a seguir é configurado para determinar um valor de raio ideal Ro para o raio de AoI Ra da AoI 107. Em uma modalidade da invenção, o valor de raio ideal Ro é determinado por meio de etapas iterativas que começam a partir de um valor de raio mínimo Rmin até um valor de raio máximo Rmax (da forma descrita a seguir). Preferivelmente, o valor de raio mínimo Rmin e o valor de raio máximo Rmax são definidos pelo administrador do sistema 100 através da interface de administrador 120.
[00120] Em uma modalidade da presente invenção, com base na análise estatística de dados empíricos relativos a uma pluralidade de acontecimentos públicos passados realizados pelo Requerente, o valor de raio mínimo Rmin é definido igual a -1.500 m (Rmin = -1.500 m), ao mesmo tempo em que o valor de raio máximo Rmax é definido igual a 1.500 m (Rmax = 1.500 m).
[00121] Tendo definido a forma dos elementos de superfície 205 da região geográfica coberta 300 coberta pela rede de comunicação móvel 105 e a forma da AoI 107, o conceito de elementos de superfície relevantes, isto é, um elemento de superfície 205q da região geográfica coberta 300 que é considerado pertencendo pelo menos parcialmente à AoI 107 de acordo com uma modalidade da invenção será agora introduzido.
[00122] As figuras 5A - 5D são áreas atendidas ou elementos de superfície 505a-d relevantes entre os elementos de superfície 205q da região geográfica coberta 300 em relação à AoI 107 de acordo com uma modalidade da invenção.
[00123] Em uma modalidade da invenção, dada a AoI 107 que tem o centro da AoI C e os elementos de superfície genéricos 205q que têm o centro de superfície B e o raio de superfície Rs, a célula genérica 105b pode ser considerada um elemento de superfície relevante 505a-d para a AoI 107 se a seguinte desigualdade for verificada: onde Dist(C, B) é a distância geográfica entre o centro da AoI C e o centro de superfície B.
[00124] De acordo com o valor do raio de AoI Ra da AoI 107, a desigualdade (1) pode tomar três diferentes significados.
[00125] A saber, se o raio de AoI Ra da AoI 107 for maior do que zero (isto é, Ra > 0), a desigualdade (1) reduz para: e o elemento de superfície genérico 205 é considerado um elemento de superfície relevante (tal como o caso dos elementos de superfície relevantes 505a na figura 4A) para a AoI 107 que tem um raio da AoI Ra maior do que zero se a área da AoI 107 e o elemento de superfície genérico 205 forem pelo menos parcialmente sobrepostos (mesmo se o centro da AoI C cair fora do elemento de superfície genérico 205).
[00126] Se o raio de AoI Ra da AoI 107 for igual a zero (isto é, Ra = 0), a desigualdade (1) reduz para: e o elemento de superfície genérico 205q é considerado um elemento de superfície relevante (tal como o caso dos elementos de superfície relevantes 505b e 505c das figuras 4B e 4C) para a AoI 107 que tem um raio da AoI Ra igual a zero se o centro de AoI C da AoI 107 estiver compreendido no elemento de superfície genérico 205q.
[00127] Finalmente, se o raio de AoI Ra da AoI 107 for menor do que zero (isto é, Ra < 0), o elemento de superfície genérico 205q é considerado um elemento de superfície relevante (tal como o caso do elemento de superfície relevante 505d na figura 4D) para a AoI 107 que tem um raio da AoI Ra menor do que zero se o centro de AoI C da AoI 107 estiver compreendido no elemento de superfície genérico 205q em uma distância do centro de superfície B igual a ou menor do que Rs - |Ra|.
[00128] Um acontecimento público (genérico) S, além de ser mantido em um local específico (isto é, a AoI 107), tem um tempo de início Ts e um tempo de fim Te. Consequentemente, com os propósitos da presente invenção, o acontecimento público S tem uma duração relevante igual a um intervalo de tempo de observação [Ts, Te] (isto é, um intervalo de tempo que inicia em um tempo de início Ts e termina no tempo de fim Te, durando Te - Ts unidades de tempo, por exemplo, segundos, minutos ou horas).
[00129] Tanto o tempo de início Ts quanto o tempo de fim Te podem ser definidos para corresponder aos tempos de início e de fim oficiais (oficialmente anunciados) agendados para este acontecimento público S. Contudo, o Requerente observou que, pela antecipação do tempo de início Ts em relação ao tempo de início oficial do acontecimento público S, é possível levar em conta o fato de que pessoas (isto é, possuidores de UE que estão presentes no acontecimento público S) chegam na AoI 107 antes do tempo de início oficial do acontecimento público S, o que pode ser útil para coletar dados sobre uma tendência no tempo de um fluxo de presentes que chegam no acontecimento público S. Por exemplo, com base nos dados empíricos de acontecimentos públicos anteriores, o Requerente descobriu que o tempo de início Ts pode ser utilmente antecipado para 60 minutos antes do tempo de início oficial do acontecimento público S a fim de levar em conta a tendência de presentes que chegam no acontecimento público S.
[00130] Similarmente, o Requerente observou que o tempo de fim Te pode ser atrasado em relação ao tempo de fim oficial do acontecimento público S a fim de levar em conta o fato de que pessoas deixam a AoI 107 depois do tempo de fim oficial do acontecimento público, o que pode ser útil para coletar dados sobre uma tendência no tempo de um fluxo de presentes que deixam o acontecimento público S. Por exemplo, com base nos dados empíricos de acontecimentos públicos anteriores, o Requerente descobriu que o tempo de fim Ts pode ser utilmente atrasado em 30 minutos depois do tempo de fim oficial do acontecimento público S a fim de levar em conta a tendência de presentes que deixam o acontecimento público S.
[00131] Todo modo, o administrador através da interface de administrador 120, e/ou o usuário através da interface de usuário 125, pode definir quaisquer tempo de início Ts e tempo de fim Te customizados para o acontecimento público S. Por exemplo, o tempo de início Ts e o tempo de fim Te podem ser definidos a fim de definir o intervalo de tempo de observação [Ts, Te] mais curto do que a efetiva duração do acontecimento público S (mais curto do que a duração da íntegra do evento) a fim de analisar um número ou uma variação de presentes em relação ao acontecimento público S apenas durante uma subparte da íntegra da duração de tempo do acontecimento público S.
[00132] Além do mais, durante o curso do acontecimento público, o administrador pode mudar em tempo real o tempo de fim Te a fim de obter um número instantâneo de presentes no acontecimento público S. Por exemplo, o administrador pode definir o tempo de fim Te em um instante de tempo atual a fim de determinar o número de presentes no acontecimento público S até o instante de tempo atual. Além do mais, o administrador pode determinar o número de presentes no acontecimento público S até uma pluralidade de subsequentes tempos de fim Te a fim de identificar uma tendência no número de presentes no acontecimento público S durante o tempo.
[00133] Tendo descrito o sistema 100 e as características temporal (isto é, o tempo de início Ts e o tempo de fim Te) e espacial (isto é, o centro de AoI C e o raio de AoI Ra da AoI 107) do acontecimento público S, um algoritmo de estimativa de multidão (ou algoritmo de contagem de multidão) de acordo com uma modalidade da presente invenção será agora descrito, pela referência às figuras 6A -6C, que são um diagrama de blocos esquemático da mesma.
[00134] Uma primeira parte do algoritmo de estimativa de multidão é configurada para determinar o valor de raio ideal Ro para o raio de AoI Ra da AoI 107 com base nos dados relativos ao acontecimento público S considerado.
[00135] Inicialmente (etapa 602), o centro da AoI C, o dia de observação g e o tempo de início Ts e tempos de fim Te do acontecimento público S são inseridos no sistema 100, por exemplo, por um usuário através da interface de usuário 125 ou pelo administrador através da interface de administrador 120.
[00136] Posteriormente (etapa 604), uma variável de iteração k é inicializada em zero (isto é, k = 0) e um valor de raio calculado Rk é inicialmente definido no valor de raio mínimo Rmin (isto é, Rk = Rmin). A variável de iteração k leva em conta o número de iterações da primeira parte do algoritmo.
[00137] A seguir (etapa 606), os elementos de superfície relevantes 505a-d para a AoI 107 que têm um raio da AoI Ra igual ao valor de raio calculado Rk (Ra = Rk) são identificados por meio da desigualdade (1), como exposto.
[00138] Posteriormente (etapa 608), um raio total Rtotk é computado pela combinação dos raios Rs dos elementos de superfície relevantes 505a-d identificados na etapa anterior 606. Por exemplo, um raio total é definido simplesmente como: onde Rsr denota o raio do r-ésimo elemento de superfície relevante 505a-d (1 < r < R, R sendo um número inteiro positivo) identificado na etapa 606.
[00139] Todos os conjuntos {uq,t} de Q números de UEs agregados uq,t referidos ao dia de observação g durante um intervalo de tempo de observação [Ts, Te] (isto é, referidos em intervalos de tempo dt compreendidos dentro do intervalo de tempo de observação [Ts, Te]) e referidos aos elementos de superfície relevantes 505a-d determinados na etapa 606 são recuperados (etapa 610) a partir do repositório 115.
[00140] Subsequentemente (etapa 612), um primeiro número de UEs Uk é computado como o número de UEs nos elementos de superfície relevantes 505a-d que foram recuperados na etapa anterior 606 com base nos conjuntos {uq,t} de Q números de UEs agregados uq,t (o primeiro número de UEs Uk depende dos elementos de superfície relevantes e, portanto, do valor de raio calculado Rk).
[00141] O primeiro número de UEs Uk pode ser computado como um número total, um número médio ou um número máximo (pico) de UEs nos elementos de superfície relevantes 505a-d, por exemplo, de acordo com uma definição selecionada pelo administrador do sistema 100 através da interface de administrador 120 e/ou pelo usuário do sistema 100 através da interface de usuário 125.
[00142] Por exemplo, o primeiro número de UEs Uk como o número total de UEs nos elementos de superfície relevantes 505a-d pode ser computado da seguinte maneira. Primeiramente, a soma dos valores dos números de UEs agregados uq,t (determinada na etapa 610) em cada um dos intervalos de tempo de referência dt compreendidos dentro do intervalo de tempo de observação [Ts, Te] em todos os elementos de superfície relevantes 505a-d (determinados na etapa 606) é computada. Subsequentemente, uma soma dos valores recentemente computados para cada um dos intervalos de tempo de referência dt compreendidos dentro do intervalo de tempo de observação [Ts, Te] é realizada. Em outras palavras, o primeiro número de UEs Uk como o número total de UEs nos elementos de superfície relevantes 505a-d pode ser computado como:
[00143] Similarmente, o primeiro número de UEs Uk como o número médio de UEs nos elementos de superfície relevantes 505a-d pode ser computado da seguinte maneira. Primeiramente, a soma dos valores dos números de UEs agregados uq,t (determinada na etapa 610) em cada um dos intervalos de tempo de referência dt compreendidos dentro do intervalo de tempo de observação [Ts, Te] em todos os elementos de superfície relevantes 505a-d (determinados na etapa 606) é computada. Subsequentemente, uma soma dos valores recentemente computados para cada um dos intervalos de tempo de referência dt compreendidos dentro do intervalo de tempo de observação [Ts, Te] é realizada. Finalmente, o valor recentemente obtido é dividido pelo número T de intervalos de tempo de referência dt compreendidos dentro do intervalo de tempo de observação [Ts, Te]. Em outras palavras, o primeiro número de UEs Uk como o número médio de UEs nos elementos de superfície relevantes 505a-d pode ser computado como:
[00144] Inversamente, o primeiro número de UEs Uk como o número máximo (pico) de UEs nos elementos de superfície relevantes 505a-d pode ser computado da seguinte maneira. Primeiramente, a soma dos valores dos números de UEs agregados uq,t (determinada na etapa 610) em cada um dos intervalos de tempo de referência dt compreendidos dentro do intervalo de tempo de observação [Ts, Te] em todos os elementos de superfície relevantes 505a-d (determinados na etapa 606) é computada. Subsequentemente, o valor máximo entre os valores recentemente computados é selecionado como o primeiro número de UEs Uk. Em outras palavras, o primeiro número de UEs Uk como o máximo (pico) UE nos elementos de superfície relevantes 505a-d pode ser computado como:
[00145] Deve-se notar que o primeiro número de UEs Uk, independente se é computado como um número total, um número médio ou um número máximo (pico) de UEs nos elementos de superfície relevantes 505a-d, é sempre dependente do valor de raio calculado Rk usado para determinar os elementos de superfície relevantes 505a-d aos quais os números de UEs agregados uq,t são referidos.
[00146] Similarmente, todos os conjuntos {u’q,t} de Q números de UEs agregados u’q,t referidos aos dias anteriores gp precedendo o dia de observação g durante o intervalo de tempo de observação [Ts, Te] e tendo ocorrido nos elementos de superfície relevantes 505a-d determinados na etapa 606 são recuperados (etapa 614) a partir do repositório 115.
[00147] Em uma modalidade da invenção, para o acontecimento público S, um conjunto de P dias anteriores gp (em que 1 < p < P e P é um número inteiro) precedendo o dia de observação g é considerado. O número P de dias anteriores gp considerado é preferivelmente definido pelo administrador (através da interface de administrador 120). Em uma modalidade da presente invenção, o administrador define o número P de dias anteriores gp de acordo com as capacidades de armazenamento do repositório 115 (isto é, a fim de poder armazenar todos os dados relativos aos P dias anteriores gp) e/ou com base nas capacidades computacionais do motor de computação 110 (isto é, a fim de poder processar todos os dados relativos aos P dias anteriores gp). Preferivelmente, o administrador define o número P de dias anteriores gp também com base em uma análise estatística de acontecimentos públicos passados do mesmo tipo (isto é, acontecimentos culturais, de entretenimento, políticos ou esportivos públicos).
[00148] O Requerente descobriu que a definição do número P de dias anteriores gp igual a 6 (isto é, P = 6) provê bons resultados para a maior parte de tipos de acontecimentos públicos (embora isto não deva ser interpretado como limitante para a presente invenção).
[00149] Então (etapa 616), um segundo número de UEs Upk é computado, para cada um dos dias anteriores gp, como o número de UEs nos elementos de superfície relevantes 505a-d com base nos conjuntos {u’q,t} de Q números de UEs agregados u’q,t referidos aos elementos de superfície relevantes 505a-d que foram recuperados na etapa anterior 606 (o segundo número de UEs Upk depende dos elementos de superfície relevantes 505a-d e, portanto, do valor de raio calculado Rk).
[00150] Similarmente ao primeiro número de UEs Uk, cada um dos segundos números de UE Upk pode ser computado como um número total, um número médio ou um número máximo (pico) de UEs nos elementos de superfície relevantes 505a-d, por exemplo, de acordo com uma definição selecionada pelo administrador do sistema 100 através da interface de administrador 120 e/ou pelo usuário do sistema 100 através da interface de usuário 125.
[00151] Por exemplo, os segundos números de UE Upk como os números totais de UE nos elementos de superfície relevantes 505a-d podem ser computados da seguinte maneira. Primeiramente, para cada um dos dias anteriores gp, a soma dos valores dos números de UEs agregados u’q,t (determinada na etapa 614) em cada um dos intervalos de tempo de referência dt compreendidos dentro do intervalo de tempo de observação [Ts, Te] durante o dia anterior gp considerado em todos os elementos de superfície relevantes 505a-d (determinados na etapa 606) é computada. Subsequentemente, uma soma dos valores recentemente computados para cada um dos intervalos de tempo de referência dt compreendidos dentro do intervalo de tempo de observação [Ts, Te] durante o dia anterior gp considerado é realizada. Em outras palavras, o segundo número de UEs Upk para o dia anterior gp como o número total de UEs nos elementos de superfície relevantes 505a-d pode ser computado como:
[00152] Os segundos números de UE Upk como os números médios de UE nos elementos de superfície relevantes 505a-d podem ser computados da seguinte maneira. Primeiramente, para cada um dos dias anteriores gp, a soma dos valores dos números de UEs agregados u’q,t (determinada na etapa 614) em cada um dos intervalos de tempo de referência dt compreendidos dentro do intervalo de tempo de observação [Ts, Te] durante o dia anterior gp considerado em todos os elementos de superfície relevantes 505a-d (determinados na etapa 606) é computada. Subsequentemente, uma soma dos valores recentemente computados para cada um dos intervalos de tempo de referência dt compreendidos dentro do intervalo de tempo de observação [Ts, Te] durante o dia anterior gp considerado é realizada. Finalmente, o valor recentemente obtido é dividido pelo número T de intervalos de tempo de referência dt compreendidos dentro do intervalo de tempo de observação [Ts, Te] durante o dia anterior gp considerado. Em outras palavras, os segundos números de UE Upk para o dia anterior gp como os números médios de UE nos elementos de superfície relevantes 505a-d podem ser computados como:
[00153] Os segundos números de UE Upk como os números máximos (pico) de UE nos elementos de superfície relevantes 505a-d podem ser computados da seguinte maneira. Primeiramente, para cada um dos dias anteriores gp, a soma dos valores dos números de UEs agregados u’q,t (determinada na etapa 614) em cada um dos intervalos de tempo de referência dt compreendidos dentro do intervalo de tempo de observação [Ts, Te] computados é selecionado como o segundo número de UEs Upk para o dia anterior gp considerado. Em outras palavras, o segundo número de UEs Upk para o dia anterior gp considerado como o máximo (pico) UE nos elementos de superfície relevantes 505a-d pode ser computado como:
[00154] Também neste caso, os segundos números de UE Upk, independente se são computados como um número total, um número médio ou um número máximo (pico) de UEs nos elementos de superfície relevantes 505a-d, são sempre dependentes do valor de raio calculado Rk usado para determinar os elementos de superfície relevantes 505a-d aos quais os números de UEs agregados uq,t são referidos.
[00155] Os segundos números de UE Upk recentemente computados são combinados (etapa 618) a fim de determinar um número médio de UEs μk (com assim, o número médio de UEs μk é claramente diferente dos segundos números de UE Upk mesmo se eles forem computados como o número médio de UEs nos elementos de superfície relevantes 505a-d) e um desvio padrão do número de UEs ok (com dos números de UE nos elementos de superfície relevantes 505a-d durante os P dias anteriores gp considerados.
[00156] O número médio de UEs μk e o desvio padrão do número de UEs ok são combinados (etapa 620) com o primeiro número de UEs Uk a fim de obter uma quantidade (estatística) definida z-escore Zk (que depende do valor de raio calculado Rk):
[00157] O z-escore Zk recentemente computado é normalizado (etapa 622) em relação ao raio total Rtotk computado na etapa 608:
[00158] A seguir, a variável k é aumentada em unidade (etapa 624; isto é, k = k + 1) e o valor de raio calculado Rk é aumentado (etapa 626): onde Δ é uma etapa ou largura da iteração que podem ser definidas pelo administrador (por exemplo, Δ = 100 m), assim, cada valor de raio calculado Rk sendo separado do próximo valor de raio calculado por uma largura de iteração Δ. Deve-se notar que a largura de iteração Δ define um máximo valor de iteração kmax para a variável de iteração k - e, portanto, um número máximo de iteração para determinar o valor de raio ideal Ro - como:
[00159] Deve-se notar que a largura de iteração Δ pode ser usada pelo administrador do sistema para ajustar uma granularidade (isto é, fineza) com a qual o valor de raio ideal Ro é determinado, isto é, quanto menor for a largura de iteração Δ definida pelo administrador, mais alto é o número de iterações definido pelo máximo valor de iteração kmax e, assim, mais fina é uma granularidade do algoritmo.
[00160] Em uma modalidade da presente invenção, já que o valor de raio mínimo Rmin é definido em - 1.500 m, o valor de raio máximo Rmax é definido em 1.500 m e a largura de iteração Δ é definida em 100 m, o máximo valor de iteração kmax para a variável de iteração k resulta como igual a 30 e, portanto, o número máximo de iterações para determinar o valor de raio ideal Ro é limitado a 30.
[00161] Posteriormente, é verificado (etapa 628) se o valor de raio calculado Rk é inferior ao, ou igual ao, valor de raio máximo Rmax:
[00162] No caso afirmativo (ramificação de saída Y do bloco de decisão 628), isto é, o valor de raio calculado Rk é inferior ao, ou igual ao, valor de raio máximo Rmax (isto é, Rk < Rmax), a operação retorna para a etapa 606 para iniciar uma nova iteração da primeira parte do algoritmo com base no valor de raio calculado Rk recentemente aumentado (na etapa 626) por uma k-ésima largura de iteração Δ adicional.
[00163] No caso negativo (ramificação de saída N do bloco de decisão 628), isto é, o valor de raio calculado Rk é maior do que o valor de raio máximo Rmax (isto é, Rk > Rmax), o valor de raio ideal Ro é computado (etapa 630) como a média dos valores de raio calculados Rk (com 1 < k < kmax) pesada pelo z-escore normalizado Z'k computado em:
[00164] As etapas 606 até 628 da primeira parte do algoritmo são iteradas até que o valor de raio calculado Rk seja maior do que o valor de raio máximo Rmax (isto é, Rk > Rmax), e o valor de raio ideal Ro seja computado (na etapa 630).
[00165] Com a computação do valor de raio ideal Ro na etapa 630, a primeira parte do algoritmo termina e uma segunda parte do algoritmo inicia (na etapa 632, descrita a seguir). No final da primeira parte do algoritmo, a AoI 107 é apropriadamente definida pelo centro da AoI C e pelo raio da AoI Ra definidos iguais ao valor de raio ideal Ro (Ra = Ro).
[00166] A segunda parte do algoritmo de acordo com uma modalidade da presente invenção é configurada para determinar um número de pessoas que está presente no acontecimento público S.
[00167] Depois que o valor de raio ideal Ro tiver sido computado, na etapa 630, um conjunto de elementos de superfície realmente relevantes 505a- d é definido (etapa 632). Este conjunto inclui todos os elementos de superfície 205 da rede de comunicação móvel 105 para os quais a desigualdade (1) é verificada quando o raio da AoI Ra for definido igual ao valor de raio ideal Ro, ou:
[00168] Então (etapa 634), todos os conjuntos {uq,t} de Q números de UEs agregados uq,t referidos ao dia de observação g durante o intervalo de tempo de observação [Ts, Te] e tendo ocorrido nos elementos de superfície realmente relevantes 505a-d determinados na etapa 632 são recuperados a partir do repositório 115.
[00169] Subsequentemente (etapa 636), um terceiro número de UEs UIAOI é computado como um número de UEs compreendidos dentro dos elementos de superfície relevantes 505a-d compreendidos na AoI 107 que tem o raio da AoI Ra igual ao valor de raio ideal Ro durante o intervalo de tempo de observação [Ts, Te] com base nos conjuntos {uq,t} que foram recuperados na etapa 632.
[00170] Similarmente ao primeiro número de UEs Uk, o terceiro número de UEs UIAOI pode ser computado como um número total, um número médio ou um número máximo (pico) de UEs nos elementos de superfície relevantes 505a-d, por exemplo, de acordo com uma definição selecionada pelo administrador do sistema 100 através da interface de administrador 120 e/ou pelo usuário do sistema 100 através da interface de usuário 125. Deve-se notar que, neste caso, o terceiro número de UEs UIAOI é computado tanto como um número total, como um número médio quanto como um número máximo (pico) de UEs nos elementos de superfície relevantes 505a-d na AoI 107 que tem o raio da AoI Ra igual ao valor de raio ideal Ro, assim, o terceiro número de UEs UIAOI é dependente do valor de raio ideal Ro em vez do valor de raio calculado Rk (do qual o primeiro número de UEs Uk é dependente).
[00171] Uma vez que o terceiro número de UEs U|AOI foi computado, um número de pessoas A é inicializado em zero (isto é, A = 0) (etapa 638). O número de pessoas A leva em conta o número de pessoas na multidão, por exemplo, presentes no acontecimento público S (da forma descrita a seguir).
[00172] Todos os conjuntos {u’q,t} referidos a cada um dos dias anteriores gp no intervalo de tempo de observação [Ts, Te] e tendo ocorrido no elemento de superfície realmente relevante 505a-d determinado na etapa 632 são recuperados (etapa 640) a partir do repositório 115.
[00173] Então (etapa 642), um quarto número de UEs Up|AOI é computado para cada um dos P dias anteriores gp como um número de UEs compreendidos dentro dos elementos de superfície relevantes 505a-d compreendidos na AoI 107 que tem o raio da AoI Ra igual ao valor de raio ideal Ro durante o intervalo de tempo de observação [Ts, Te] com base nos conjuntos {u’q,t} que foram recuperados na etapa 640.
[00174] Também neste caso, similarmente aos segundos números de UE Upk, cada um dos quartos números de UE Up|AOI pode ser computado como um número total, um número médio ou um número máximo (pico) de UEs nos elementos de superfície relevantes 505a-d, por exemplo, de acordo com uma definição selecionada pelo administrador do sistema 100 através da interface de administrador 120 e/ou pelo usuário do sistema 100 através da interface de usuário 125. Deve-se notar que, neste caso, os quartos números de UE Up|AOI são computados tanto como um número total, como um número médio quanto como um número máximo (pico) de UEs nos elementos de superfície relevantes 505a-d na AoI 107 que tem o raio da AoI Ra igual ao valor de raio ideal Ro; assim, os quartos números de UE Up|AOI são dependentes do valor de raio ideal Ro em vez do valor de raio calculado Rk (do qual os segundos números de UE Upk são dependentes).
[00175] Os quartos números de UE Up|AOI recentemente computados são combinados (etapa 644) a fim de determinar um número médio adicional de UEs μlAoi do número de UEs nos elementos de superfície relevantes 505ad. Por exemplo, o número médio adicional de UEs μ|AOI pode ser computado como:
[00176] O número médio adicional de UEs μlAOI provê uma indicação de um número médio de pessoas normalmente compreendidas na AoI 107 que tem o raio da AoI Ra igual ao valor de raio ideal Ro durante o intervalo de tempo de observação considerado [Ts, Te] (isto é, pessoas que não se reúnem na multidão).
[00177] É apontado que, embora o número médio adicional de UEs μnlAOI computado da forma supradescrita possa prover um tipo de precisão limitada (já que duas ou mais atividades de um mesmo UE no intervalo de tempo de observação considerado [Ts, Te] podem ser consideradas como, cada qual, pertencendo a diferentes UEs), o número médio adicional de UEs μnlAOI provê uma estimativa do número médio de pessoas normalmente compreendidas na AoI 107 que tem uma precisão satisfatória provida com uma baixa complexidade computacional e que garante um completo respeito da privacidade dos possuidores de UE.
[00178] O número de pessoas A é, então, computado (etapa 646) pela combinação (por exemplo, subtração) do número médio adicional de UEs μlAOI determinado na etapa 644 com o terceiro número de UEs UlAOI determinado na etapa 636, ou
[00179] O número de pessoas A referidas ao acontecimento público S mantido no dia de observação g é armazenado (etapa 648) no repositório 115.
[00180] Preferivelmente, o algoritmo de estimativa de multidão é terminado com a provisão (etapa 650) dos resultados, isto é, o número de pessoas A em relação ao usuário através do terminal de usuário 125 para inspeção e/ou processamento adicionais.
[00181] Vantajosamente, o algoritmo de estimativa de multidão de acordo com a presente invenção pode ser configurado para (possivelmente automaticamente) atualizar em tempo real o número de pessoas A durante a íntegra da duração da reunião de pessoas nas multidões. Por exemplo, o número de pessoas An pode ser atualizado (tanto periodicamente quanto de forma assíncrona) toda vez que novos dados agregados forem providos para o sistema 100 pela rede de comunicação móvel 105.
[00182] Em resumo, o algoritmo de estimativa de multidão (ou algoritmo de contagem de multidão) compreende uma primeira parte e uma segunda parte.
[00183] A primeira parte do algoritmo de estimativa de multidão compreende um primeiro ciclo que escaneia (etapas 608 - 628) todos os valores de raio calculados Rk entre o valor de raio mínimo Rmin e o valor de raio máximo Rmax do acontecimento público. Para cada valor de raio Rk computado, respectivos elementos de superfície relevantes 505a-d e um z- escore normalizado Z'k são determinados. Com base em tais dados (isto é, respectivos elementos de superfície relevantes 505a-d e o z-escore normalizado Z'k), o valor de raio ideal Ro é identificado. No final da primeira parte do algoritmo de estimativa de multidão, a AoI 107 que tem o valor de raio ideal Ro é definida.
[00184] A segunda parte do algoritmo de estimativa de multidão computa o número de UEs (isto é, o terceiro número U|AOI) nos elementos de superfície relevantes 505a-d (isto é, compreendidos na AoI 107) durante o intervalo de tempo de observação [Ts, Te] no dia de observação g, computa o número de UEs (isto é, o quarto número Up|AOI) nos mesmos elementos de superfície relevantes 505a-d (isto é, compreendidos na AoI 107) durante o intervalo de tempo de observação [Ts, Te] em cada um dos P dias anteriores gp, um número médio do último (isto é, o número médio adicional μlAoi) e, eventualmente, o número de pessoas reunidas em uma multidão (isto é, o número de pessoas A) com base nos números de UE nos elementos de superfície relevantes 505a-d durante o intervalo de tempo de observação [Ts, Te] no dia de observação g e no número médio do UE.
[00185] O sistema de estimativa de multidão 100 e o algoritmo de estimativa de multidão de acordo com uma modalidade da presente invenção permitem uma computação em tempo real do número de pessoas reunidas em uma multidão, por exemplo, para estar presente em um acontecimento público S, de uma maneira confiável e que identifica apropriadamente (pela determinação do valor de raio ideal Ro) uma efetiva extensão da AoI 107 associada com o acontecimento público S.

Claims (21)

1. Método para estimar um número de pessoas reunidas em uma Área de Interesse (107) durante um intervalo de tempo em um dia, em que a Área de Interesse (107) é definida por um centro de Área de Interesse e um raio de Área de Interesse e é coberto por uma rede de telecomunicação móvel (105) incluindo uma pluralidade de estações de comunicação (105a) cada uma das quais configurada para gerenciar comunicações de Equipamento de Usuário em uma ou mais áreas atendidas (105b) em uma região geográfica coberta (300; 300') pela qual a rede de telecomunicação móvel (105) provê serviços, o método caracterizado pelo fato de que compreende as etapas de: a)subdividir a região geográfica coberta (300; 300') em uma pluralidade de elementos de superfície (205) correspondendo a células da rede de telecomunicação móvel; b) definir (604) uma pluralidade de valores de raio calculados do raio de Área de Interesse, e, para cada valor de raio calculado: c) identificar (606) um número de elementos de superfície relevantes (505a-d) da região geográfica coberta (300; 300') dentro da Área de Interesse; d) computar (610; 612) um primeiro número de Equipamento de Usuário atendido pela rede de comunicação móvel (105) durante o intervalo de tempo no dia dentro da Área de Interesse (107) com base em dados agregados indicando uma razão entre carga de tráfico e tráfico de carga de equipamento de usuário médio da rede de comunicação móvel; e) computar (614; 616) um segundo número de Equipamento de Usuário atendido pela rede de comunicação móvel (105) durante o intervalo de tempo para cada dia de um número predeterminado de dias anteriores precedendo o dia dentro da Área de Interesse (107) com base nos dados agregados; f) combinar (618, 620) o primeiro número de Equipamento de Usuário e os segundos números de Equipamento de Usuário para obter uma quantidade estatística; g) normalizar a quantidade estatística para obter uma quantidade estatística normalizada, a quantidade estatística sendo normalizada com relação ao raio das áreas atendidas relevantes; h) computar (630) um valor de raio ideal do raio de Área de Interesse como a média dos valores de raio calculados pesados pela quantidade estatística normalizada; i) estimar (632-646) o número de pessoas reunidas dentro da Área de Interesse (107) tendo o raio de Área de Interesse igual ao valor de raio ideal.
2. Método de acordo com a reivindicação 1, caracterizado pelo fato de que um elemento de superfície (205) é identificado como um elemento de superfície relevante (505a-d) se ele verifica a seguinte desigualdade: Dist(C, B) < |Rs + Rk|, onde C é o centro da Área de Interesse (107), B é o centro do elemento de superfície (205), Dist(C, B) é a distância geográfica entre o centro da Área de Interesse C e o centro do elemento de superfície B, Rs é o raio do elemento de superfície, e Rk é o valor de raio calculado.
3. Método de acordo com qualquer uma das reivindicações 1 ou 2, caracterizado pelo fato de que compreende adicionalmente para cada valor de raio calculado: j) receber uma pluralidade de dados agregados indicando razões entre cargas de tráfico e cargas de tráfico de equipamento de usuário médio referidos a cada um dos elementos de superfície.
4. Método de acordo com a reivindicação 3, caracterizado pelo fato de que a etapa j) de receber uma pluralidade de dados agregados para cada um dos elementos de superfície (205) compreende: receber um conjunto de dados agregados, cada dado agregado do conjunto dos dados agregados sendo referido a um intervalo de tempo de referência respectivo o qual é uma porção de um período de aquisição durante o qual dados agregados são coletados.
5. Método de acordo com a reivindicação 4, caracterizado pelo fato de que a etapa c) de computar (610, 612) um primeiro número de Equipamento de Usuário atendido pela rede de comunicação móvel (105) durante o intervalo de tempo de observação no dia dentro da Área de Interesse (107) com base em dados agregados compreende: computar (610, 612) um primeiro número de Equipamento de Usuário com base em conjuntos de dados agregados referidos a intervalos de tempo de referência respectivos dentro do intervalo de tempo de observação no dia; e em que a etapa d) de computar (614, 616) um segundo número de Equipamento de Usuário que foi atendido pela rede de comunicação móvel (105) durante o intervalo de tempo de observação para cada dia de um número predeterminado de dias anteriores precedendo o dia dentro da Área de Interesse (107) com base nos dados agregados compreende: computar (614, 616) cada segundo número de Equipamento de Usuário com base em conjuntos de dados agregados referidos a intervalos de tempo de referência respectivos dentro do intervalo de tempo de observação do respectivo dia anterior do número predeterminado de dias anteriores precedendo o dia.
6. Método de acordo com a reivindicação 5, caracterizado pelo fato de que o primeiro número de Equipamento de Usuário e/ou cada segundo número de Equipamento de Usuário pode ser computado como número total, um número médio, ou um número máximo de Equipamento de Usuário nos elementos de superfície relevantes (505a-d) na Área de Interesse (107).
7. Método de acordo com qualquer uma das reivindicações 1 a 6, caracterizado pelo fato de que a etapa e) de combinar (618, 620) o primeiro número de Equipamento de Usuário e os segundos números de Equipamento de Usuário para obter uma quantidade estatística compreende: combinar (618) o segundo número de Equipamento de Usuário de cada um dos dias anteriores de maneira a determinar um número de Equipamento de Usuário médio e um desvio padrão de número de Equipamento de Usuário.
8. Método de acordo com a reivindicação 7, caracterizado pelo fato de que a etapa e) de combinar (618, 620) o primeiro número de Equipamento de Usuário e os segundos números de Equipamento de Usuário para obter uma quantidade estatística compreende adicionalmente: computar a quantidade estatística como: Zk = (Uk - μk)/ok, em que Uk é o primeiro número de Equipamento de usuário, μk é o número de Equipamento de Usuário médio e ok é o desvio padrão de número de Equipamento de Usuário.
9. Método de acordo com qualquer uma das reivindicações 1 a 8, caracterizado pelo fato de que a pluralidade de valores de raio calculados varia de um valor de raio mínimo a um valor de raio máximo, cada valor de raio calculado sendo separado de um próximo valor de raio por uma largura de iteração.
10. Método de acordo com qualquer uma das reivindicações 1 a 9, caracterizado pelo fato de que a etapa h) de estimar (632-646) um número de pessoas reunidas dentro da Área de Interesse (107) tendo o raio de Área de Interesse igual ao valor de raio ideal compreende: k) definir (632) um número de elementos de superfície relevantes (505a-d) dentre os elementos de superfície (205) na região geográfica coberta (300; 300'), em que os elementos de superfície relevantes (505a-d) são elementos de superfície pelo menos parcialmente sobrepostos na Área de Interesse (107) tendo o raio de Área de Interesse igual ao valor de raio ideal.
11. Método de acordo com a reivindicação 10, caracterizado pelo fato de que um elemento de superfície (205) é identificado como um elemento de superfície relevante (505a-d) se ele verifica a seguinte desigualdade: Dist(C, B) < IRs + Rol, onde C é o centro da Área de Interesse (107), B é o centro do elemento de superfície (205), Dist(C, B) é a distância geográfica entre o centro da Área de Interesse C e o centro do elemento de superfície B, Rs é o raio do elemento de superfície, e Ro é o valor de raio ideal.
12. Método de acordo com qualquer uma das reivindicações 10 ou 11, caracterizado pelo fato de que a etapa i) de estimar (632-646) um número de pessoas reunidas dentro da Área de Interesse (107) tendo o raio de Área de Interesse igual ao valor de raio ideal compreende adicionalmente: computar (634, 636) um terceiro número de Equipamento de Usuário como um número de Equipamento de Usuário dentro dos elementos de superfície relevantes (505a-d) na Área de Interesse (107) tendo o raio de Área de Interesse igual ao valor de raio ideal durante o intervalo de tempo com base nos dados agregados em relação a um uso da rede de comunicação móvel (105).
13. Método de acordo com a reivindicação 12, caracterizado pelo fato de que computar (634, 636) um terceiro número de Equipamento de usuário compreende: computar (634, 636) o terceiro número de Equipamento de Usuário com base em conjuntos de dados agregados referidos a intervalos de tempo de referência respectivos dentro do intervalo de tempo de observação no dia.
14. Método, de acordo com a reivindicação 12, caracterizado pelo fato de que o terceiro número de Equipamento de Usuário e/ou cada quarto número de Equipamento de Usuário pode ser computado como um número total, um número médio ou um número máximo de Equipamento de Usuário nos elementos de superfície relevantes na Área de interesse tendo o raio da Área de Interesse igual ao valor de raio ideal.
15. Método de acordo com qualquer uma das reivindicações 10 a 14, caracterizado pelo fato de que a etapa i) de estimar (632-646) um número de pessoas reunidas dentro da Área de Interesse (107) tendo o raio de Área de Interesse igual ao valor de raio ideal compreende adicionalmente: computar (640, 642) um quarto número de Equipamento de Usuário como um número de Equipamento de Usuário dentro dos elementos de superfície relevantes (505a-d) na Área de Interesse (107) tendo o raio de Área de Interesse igual ao valor de raio ideal para cada dia do número predeterminado de dias anteriores precedendo o dia com base nos dados agregados.
16. Método de acordo com a reivindicação 15, caracterizado pelo fato de que computar (640, 642) um quarto número de Equipamento de usuário compreende: computar (640, 642) cada quarto número de Equipamento de Usuário com base em conjuntos de dados agregados referidos a intervalos de tempo de referência respectivos dentro do intervalo de tempo de observação do respectivo dia anterior do número predeterminado de dias anteriores precedendo o dia.
17. Método de acordo com a reivindicação 15, caracterizado pelo fato de que a etapa i) de estimar (632-646) um número de pessoas reunidas dentro da Área de Interesse (107) tendo o raio de Área de Interesse igual ao valor de raio ideal compreende adicionalmente: combinar (644) o quarto número de Equipamentos de Usuário de cada um dos dias anteriores de maneira a determinar um número de Equipamento de Usuário médio adicional, o número de Equipamento de Usuário médio adicional provendo uma indicação de um número médio de pessoas normalmente compreendidas dentro da Área de Interesse (107) tendo o raio de Área de Interesse igual aos valores de raio ideais durante o intervalo de tempo de observação considerado em quaisquer dias.
18. Método de acordo com a reivindicação 17, caracterizado pelo fato de que a etapa i) de estimar (632-646) um número de pessoas reunidas dentro da Área de Interesse (107) tendo o raio de Área de Interesse igual ao valor de raio ideal compreende adicionalmente: combinar (646) o terceiro número de Equipamento de Usuário e o número de Equipamento de Usuário médio adicional de maneira a obter o número de pessoas reunidas na Área de Interesse (107) tendo o raio de Área de Interesse igual ao valor de raio ideal.
19. Método de acordo com a reivindicação 18, caracterizado pelo fato de que combinar (646) o terceiro número de Equipamento de Usuário e o número de Equipamento de Usuário médio adicional compreende subtrair o número de Equipamento de Usuário médio adicional do terceiro número de Equipamento de usuário.
20. Sistema (100) acoplado com uma rede de telecomunicação sem fio (105) para estimar um número de pessoas reunidas em uma Área de Interesse (107), o sistema compreendendo: um motor de computação (110) configurado para processar dados recuperados de uma rede de telefonia móvel (105); um repositório (115) configurado para armazenar dados em relação a interações entre o Equipamento de Usuário e a rede de telefonia móvel, resultados de computação gerados pelo motor de computação e dados de processamento gerados por e/ou providos ao sistema, e uma interface de administrador (120) configurada para modificar parâmetros e/ou algoritmos usados pelo motor de computação e/ou acessar dados armazenados no repositório, caracterizado pelo fato de que compreende adicionalmente um elemento de memória (110a) armazenando instruções configuradas para implementar o método como definido em qualquer uma das reivindicações 1 a 19 através do sistema (100).
21. Sistema de acordo com a reivindicação 20, caracterizado pelo fato de que compreende adicionalmente pelo menos uma interface de usuário (125) configurada para receber entradas a partir de, e para prover saída para um usuário do sistema, o usuário compreendendo um ou mais seres humanos e/ou um ou mais sistemas de computação externos assinantes dos serviços providos pelo sistema.
BR112017020838-5A 2015-03-30 Método para estimar um número de pessoas reunidas em uma área de interesse, e, sistema acoplado com uma rede de telecomunicação sem fio para estimar um número de pessoas reunidas em uma área de interesse BR112017020838B1 (pt)

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
PCT/EP2015/056878 WO2016155769A1 (en) 2015-03-30 2015-03-30 Method and system for a real-time counting of a number of persons in a crowd by means of aggregated data of a telecommunication network

Publications (2)

Publication Number Publication Date
BR112017020838A2 BR112017020838A2 (pt) 2018-07-03
BR112017020838B1 true BR112017020838B1 (pt) 2023-05-30

Family

ID=

Similar Documents

Publication Publication Date Title
EP3241366B1 (en) Method and system for estimating a number of persons in a crowd
US10375516B2 (en) Method and system for a real-time counting of a number of persons in a crowd by means of aggregated data of a telecommunication network
US10285002B2 (en) Method and system for a real-time counting of a number of participants at a public happening
Qiao et al. A mobility analytical framework for big mobile data in densely populated area
EP3132592B1 (en) Method and system for identifying significant locations through data obtainable from a telecommunication network
BR112017020831B1 (pt) Método e sistema para estimar um número de pessoas que se reuniram em uma área de interesse
US10278015B2 (en) Method and system for real-time computing of origin-destination matrices relating to attendees at a public happening through analysis of mobile communication network data
BR112016002420B1 (pt) método e sistema para gerenciar dados em relação a um ou mais fluxos de entidades físicas e de elementos em uma área geográfica
US9955299B2 (en) Method and system for a posteriori computation of origin-destination matrices relating to gathering of people through analysis of mobile communication network data
Sumathi et al. Crowd estimation at a social event using call data records
US10887728B2 (en) Method and system for determining the mobility status of a user of a wireless communication network
BR112017020838B1 (pt) Método para estimar um número de pessoas reunidas em uma área de interesse, e, sistema acoplado com uma rede de telecomunicação sem fio para estimar um número de pessoas reunidas em uma área de interesse
Alepis et al. Mapping the wireless coverage grid for carrier and user recommendations through crowdsourcing
Dešić et al. Determination of origins and destinations for an OD matrix based on telecommunication activity records