12

Quero mostrar uma rede celular em um mapa. Os dados de entrada são um arquivo .csv em que cada sequência é um setor celular. Os atributos são: ID do setor, suas coordenadas, seu azimute e um ângulo de largura do feixe da antena.

Os valores da largura do feixe da antena variam de 30 a 360 graus. A largura do feixe de 360 antenas significa que ele deve ser mostrado no mapa como um círculo. Antenas com outras larguras de feixe devem ser mostradas como setores com ângulos de abertura apropriados.

É possível mostrar as antenas usando apenas símbolos? Eu sei como criar meu próprio símbolo SVG e espero descobrir como girá-lo de acordo com o azimute. Mas existe alguma maneira de aplicar a largura variável do feixe da antena de acordo com seu valor de atributo de 30 a 360 graus?

Eu acho que os símbolos são a melhor maneira de desenhar as antenas por causa da visualização dinâmica no mapa de acordo com a escala de visão, se isso for possível no QGIS.

Obviamente, a tarefa pode ser resolvida desenhando polígonos apropriados como um recurso de camada, mas isso seria uma solução alternativa.

qgis symbology style visualisation E Bobrov
fonte

Então você precisa desenhar o arco na direção certa, diferente para cada site?

22716 Nathan W

Nem um pouco Se eu entendi direito. Ele deve ser um setor do círculo (ou o círculo inteiro no caso de largura de feixe = 360), conforme mostrado na figura.

E Bobrov

Sim, é isso que eu quero dizer.

Nathan W

OK eu vejo. Falando em geral, o símbolo do arco não é estritamente necessário. Os principais atributos são azimute e largura de feixe. Eu posso usar qualquer símbolo para desenhar as antenas, não apenas um arco.

E Bobrov

Provavelmente encontrei o exemplo que pode ajudar: Criando tipos de camadas de símbolos personalizados . Mas não tenho certeza. Então, alguém tentou criar sua própria classe de camadas de símbolos que desenha, por exemplo, a direção de cada camada, dependendo do seu atributo (ou seja, azimute da antena nas palavras da figura acima)?

E Bobrov

7

Alguns dias atrás, um novo plugin foi adicionado ao QGIS chamado Wedge Buffer Processing Algorithm . Parece que pode ser do seu interesse.

Como o nome sugere, é um algoritmo de processamento, então você precisará executá-lo na caixa de ferramentas de processamento. Ainda não tive a chance de experimentá-lo.

Ele cria setores de círculos - como um buffer circular normal, mas o ângulo e o raio da cunha podem ser definidos usando valores de campo.

Documentação e capturas de tela podem ser vistas na página do github

Steven Kay
fonte

10

Se você deseja usar apenas a simbologia, proponho uma solução inspirada na minha resposta de uma pergunta semelhante: Criando luzes do setor no QGIS? .

Seguindo uma abordagem semelhante, e supondo que você esteja trabalhando em um CRS projetado (em vez disso, se estiver usando um sistema de coordenadas geográficas, veja a nota no final da resposta), quero sublinhar que vou focar a atenção na explicação das coisas mínimas a serem executadas para reproduzir o resultado desejado: isso significa que alguns outros parâmetros menores (como tamanhos, larguras etc.) devem ser facilmente ajustados por você para melhor atender às suas necessidades.

Além disso, eu assumo que "AZIMUTH"é o campo que armazena os valores de azimute e "BEAMWIDTH"é o campo que armazena as larguras de feixe da antena.

Solução

Renderizaremos os pontos com Single symbolae recorrendo a uma Simple Markere duas Geometry generatorcamadas de símbolos:

Na explicação adicional, seguirei a mesma ordem dos símbolos na imagem acima.

1) Marcador Simples

Eu escolhi um símbolo padrão de um círculo vermelho (essa é a parte mais fácil deste tutorial), com um tamanho de 3 mm e uma largura de 0,4 mm.

2) Gerador de geometria n ° 1

Adicione uma nova camada de símbolo e selecione Geometry generatoros LineString / MultiLineStringtipos e:

Insira esta expressão no Expressioncampo:

make_line(
 $geometry,
 make_point($x + 300*cos(radians(90 -  "AZIMUTH" )), $y + 300*sin(radians((90 - "AZIMUTH" ))))
)

Acabamos de definir a seta que aponta para o conjunto de azimute (para criar a seta, lembre-se de selecionar o Arrowtipo de camada de símbolo sob a Lineopção no menu principal do símbolo). Observe que isso 300representa uma distância em metros e é um valor arbitrário; portanto, fique à vontade para alterá-lo de acordo com suas necessidades.

3) Gerador de Geometria No. 2

Adicione uma nova camada de símbolo e selecione o Geometry generatortipo e os Polygon / MultiPolygontipos:

Insira esta expressão no Expressioncampo:

CASE
WHEN ("BEAMWIDTH") <= 180
THEN
intersection(
  buffer(
   $geometry, 200),
  make_polygon(
   geom_from_wkt(
    geom_to_wkt(
     make_line(
      $geometry,
      make_point($x + 2000*cos(radians(90 -  "AZIMUTH" - "BEAMWIDTH"/2 )), $y + 2000*sin(radians((90 - "AZIMUTH" - "BEAMWIDTH"/2 )))),
      make_point($x + 2000*cos(radians(90 -  "AZIMUTH" )), $y + 2000*sin(radians((90 - "AZIMUTH" )))),
      make_point($x + 2000*cos(radians(90 - "AZIMUTH" + "BEAMWIDTH" /2)), $y + 2000*sin(radians((90 - "AZIMUTH" + "BEAMWIDTH"/2)))),
      $geometry)
     )
    )
   )
  )

WHEN ("BEAMWIDTH") > 180
THEN
difference(
  buffer(
   $geometry, 200),
   make_polygon(
    geom_from_wkt(
     geom_to_wkt(
      make_line(
       $geometry,
       make_point($x + 2000*cos(radians(90 -  "AZIMUTH" - "BEAMWIDTH"/2 )), $y + 2000*sin(radians((90 - "AZIMUTH" - "BEAMWIDTH"/2 )))),
       make_point($x - 2000*cos(radians(90 -  "AZIMUTH" )), $y - 2000*sin(radians((90 - "AZIMUTH" )))),
       make_point($x + 2000*cos(radians(90 - "AZIMUTH" + "BEAMWIDTH" /2)), $y + 2000*sin(radians((90 - "AZIMUTH" + "BEAMWIDTH"/2)))),
       $geometry)
      )
     )
    )
   )

END

Acabamos de definir o setor. Observe que 200e 2000representa distâncias em metros e são valores arbitrários, porque estou tentando criar um polígono para cruzar com o círculo com um raio de 200 m, portanto, fique à vontade para alterá-los de acordo com suas necessidades.

Resultado final

Se você executar corretamente as tarefas anteriores, poderá obter resultados como estes (os rótulos são adicionados além desta solução e devem explicar melhor apenas o contexto):

Nota

Se você estiver usando um Sistema de coordenadas geográficas , ou seja, se estiver lidando com graus e não com distâncias, deve ser suficiente usar os valores adequados quando usei uma distância nas fórmulas anteriores. As distâncias que usei são:

300 m (consulte Geometry Generator No. 1);
200 m (consulte Geometry Generator No. 2);
2000 m (consulte Geometry Generator No. 2);

assim você pode substituí-lo com outros valores arbitrários expressas em graus (por exemplo, 0.0002, 0.002e assim por diante).

Bônus

Anexei o estilo aqui : você pode abrir esse código com qualquer editor de texto e salvá-lo como um arquivo de estilo de camada QGIS (ou seja, com uma .qmlextensão).

O estilo acima foi criado usando o QGIS 2.18.4 (ele deve ter o mesmo nome do shapefile que você está usando).

mgri
fonte

Você estava procurando esta solução? Funciona?

Mgri

Sua solução resolve totalmente o caso descrito no tópico! Eu o implementei e entendi que meu próprio caso real é um pouco diferente. Desculpe, a culpa é minha.

E Bobrov

1) A densidade dos meus setores no mapa é diferente, ou seja, no caso de distâncias curtas entre os setores, a distância definida no código fará com que muitos setores se sobreponham, alterar o zoom do mapa não ajudaria, por isso será bastante difícil para ler o mapa. Porém, nos casos de longas distâncias entre setores, os setores serão muito pequenos e talvez difíceis de ler o mapa. O uso de símbolos únicos está livre desse problema, suas escalas mudam com o zoom do mapa.

E Bobrov

2) E existem distorções da largura de feixe: a largura de feixe em 360 graus se parece com elipses, setores com azimutes diferentes, mas a mesma largura de feixe não se parece com setores com larguras de feixe equivalentes. É porque estou usando o Sistema de Coordenadas Geográficas? Agora, diferentes ângulos de longitude / latitude representam diferentes distâncias entre os pontos da Terra. Portanto, a solução precisava ser localizada nas áreas da Terra onde os setores estão localizados.

E Bobrov

Enfim, sua solução e a referência para uma resposta semelhante "Criando luzes do setor no QGIS?" me ajudou a ver algumas funcionalidades úteis. Mais uma vez obrigado.

E Bobrov

4

Muitos elogios para mgri.

Em nossa camada de teste, tudo funcionou sem problemas. Em uma camada de produção, depois de duas / três horas, consegui rastrear um problema com $ geometry . Tinha exportado uma camada de pontos de uma plataforma, não percebeu, mas era o MultiPoint . Isso parecia causar problemas: a flecha não estava desenhada; e estranhamente apenas os pontos calculados formavam o polígono dos círculos.

Outra coisa é que estou usando um raio variável . (não tenho certeza se é a palavra certa nesse caso, você também pode chamá-lo de 'comprimento do feixe' ou o que for).

Aqui está o que estou usando agora, com uma camada de tipo de geometria MultiPoints (embora, de fato, todos os recursos sejam um único ponto), e funcione para mim no QGis 2.18.3

Expressão da seta Nenhuma seta se 360 °.

CASE

WHEN ("BEAMWIDTH") = 360
THEN 
make_line(
 make_point($x, $y),
 make_point($x + "RADIUS"*cos(radians(90 -  "AZIMUTH" )), $y + "RADIUS"*sin(radians((90 - "AZIMUTH" ))))
)

END

Expressão de polígono

CASE

WHEN ("BEAMWIDTH") <= 180
THEN
intersection(
  buffer(
   make_point($x,$y), "RADIUS"),
  make_polygon(
   geom_from_wkt(
    geom_to_wkt(
     make_line(
      make_point($x,$y),
      make_point($x + "RADIUS"*2*cos(radians(90 -  "AZIMUTH" - "BEAMWIDTH"/2 )), $y + "RADIUS"*2*sin(radians((90 - "AZIMUTH" - "BEAMWIDTH"/2 )))),
      make_point($x + "RADIUS"*2*cos(radians(90 -  "AZIMUTH" )), $y + "RADIUS"*2*sin(radians((90 - "AZIMUTH" )))),
      make_point($x + "RADIUS"*2*cos(radians(90 - "AZIMUTH" + "BEAMWIDTH" /2)), $y + "RADIUS"*2*sin(radians((90 - "AZIMUTH" + "BEAMWIDTH"/2)))),
      make_point($x,$y))
     )
    )
   )
  )

WHEN ("BEAMWIDTH") > 180
THEN
difference(
  buffer(
   make_point($x,$y), "RADIUS"),
   make_polygon(
    geom_from_wkt(
     geom_to_wkt(
      make_line(
       make_point($x,$y),
       make_point($x + "RADIUS"*2*cos(radians(90 -  "AZIMUTH" - "BEAMWIDTH"/2 )), $y + "RADIUS"*2*sin(radians((90 - "AZIMUTH" - "BEAMWIDTH"/2 )))),
       make_point($x - "RADIUS"*2*cos(radians(90 -  "AZIMUTH" )), $y - "RADIUS"*2*sin(radians((90 - "AZIMUTH" )))),
       make_point($x + "RADIUS"*2*cos(radians(90 - "AZIMUTH" + "BEAMWIDTH" /2)), $y + "RADIUS"*2*sin(radians((90 - "AZIMUTH" + "BEAMWIDTH"/2)))),
       make_point($x,$y))
      )
     )
    )
   )

END

jbostoen
fonte

Minha resposta propôs uma abordagem geral: como havia muitas variáveis envolvidas no problema, era praticamente impossível criar um procedimento único para resolver qualquer situação. Portanto, obrigado por apontá-lo e propor uma abordagem com os recursos do MultiPoint, certamente ajudará alguém no futuro.

MGRI

1

Eu fui presenteado com uma solução parcial na web sem nenhum plug-in adicional, apenas qgis fora da caixa. Ele não mostra a largura do feixe da antena, basta girar o marcador simples na direção certa: use o marcador simples e gire-o com azimute da antena + 180 graus (Propriedades da camada> Simbol único-> Marcador-> Marcador simples-> triângulo-> rotação-> editar -> digite <180 + "azimute da antena"> no campo de expressão e defina Top no campo do ponto de ancoragem do marcador). O uso de <180 + "antena azimute"> é necessário devido à direção inadequada do marcador de triângulo simples incorporado. Caso contrário, mostrará uma direção incorreta da antena.

E Bobrov
fonte

Mostrando símbolos de antena no mapa: símbolos de pontos ou recursos (polígonos)

Respostas:

Solução

1) Marcador Simples

2) Gerador de geometria n ° 1

3) Gerador de Geometria No. 2

Resultado final

Nota

Bônus