Como fazer um mapa no R

Este é o post de estreia da Renata, nossa nova colaboradora. Que tal aprender a fazer mapas básicos no R?

Kia ora pessoal,

Meu nome é Renata Muylaert e a partir de agora faço parte do time Sobrevivendo na Ciência! Obrigada, Marco e Alexandre, pelo convite!

Quem me conhece sabe o quanto eu amo mapas! E quem está elaborando projetos de pesquisa em diversas áreas (a minha é a Ecologia) sabe o quão valioso pode ser um mapa para o planejamento da coleta de dados e para a interpretação dos resultados.

Mas fazer um mapa não é tão trivial. É comum ter dificuldade para entender coordenadas geográficas, formatos de dados espaciais e até mesmo a portabilidade dos mesmos. Pegar os pontos que você marcou com um GPS em campo e, a partir deles, plotar um mapa eficiente e bonito pode ser um caminho longo.

Neste post não vou focar nas bases da cartografia e dos sistemas de coordenadas geográficas (fica para outro post; você pode ler o básico aqui). Este post é para quem quer sentar e fazer um mapa rapidamente no R com o mínimo de sofrimento possível.

Como exemplo, vou explicar como plotar em um mapa, com direito a grade, escala e norte, alguns locais que visitei na Nova Zelândia (🇳🇿), um país incrível no qual eu tive oportunidade de fazer parte do meu doutorado. Você pode adaptar este script para plotar os pontos de coleta da sua tese ou outros dados similares.

Começando

Se você já tem experiência no R, basta baixar o repositório que eu criei no GitHub, plotar o mapa a partir do script comentado e adaptá-lo para os seus dados e gostos.

Se você for mais noob, recomendo que leia com atenção o post inteiro e os comentários no script.

Coordenadas (latitude e longitude)

Suponhamos que você pegou as coordenadas no sistema mais comumente utilizado, o WGS 84. As informações espaciais serão compostas de no mínimo três colunas: latitude, longitude e identificação do ponto.

Embora sejam somente três inocentes colunas, muitas vezes dados espaciais acabam apresentando intrusões demoníacas após saírem do GPS ou de uma tabela fonte. Ao serem colados nos comumente utilizados MS Excel ou LibreOffice Calc, os dados podem ficar completamente zicados. Isso ocorre se você não especifica corretamente se seu separador decimal é ponto ou vírgula, o que pode depender do Idioma do seu sistema operacional. Quem já tiver passado por isso, comente aqui no post!

Neste exemplo, eu utilizo um arquivo de bloco de notas (.txt) separado somente e somente só por vírgulas (,) como input de coordenadas. E o separador decimal é ponto (.). O Excel em português costuma não entender coordenadas separadas por pontos e “zoa” todos seus dados rapidamente. Portanto, tenha muita paciência ao formatar as suas planilhas. É incrível como até pessoas com anos de experiência trabalhando com dados espaciais podem ainda passar por esses mesmos apertos.

Enfatizo, não tenha pressa para abrir os dados no R. Outro problema bastante comum é a simples troca da coluna da latitude pela da longitude. Essa troca é resultado de: (1) pressa; (2) falta de familiaridade do usuário (no caso você!) com coordenadas geográficas; (3) falta de familiaridade da pessoa que te enviou a planilha com coordenadas geográficas; ou (4) ausência de revisão minuciosa dos dados.

O que mais pode atrapalhar

Garantindo essa etapa principal, provavelmente o resto do script funcionará. Caso contrário, outro possível problema é a sua versão do R ser antiga demais. E repare que “antiga demais” pode não ser tão antiga assim. No momento presente, já vi que tem duas versões do R mais novas do que a que eu estou usando. Às vezes, basta uma atualização menor do R para um script perder parte da compatibilidade quando os pacotes não acompanham a atualização.

No script do exemplo, o pacote ggsn pode não funcionar caso você tenha uma versão antiga do R. No exemplo aqui utilizei o R 3.5.1 (“Feather Spray”) e a nova versão do ggsn (​v0.5.0). Caso você tenha problemas em carregar este pacote seguindo o script, reporte aqui! E caso funcione bem, comente também, é importante pra gente. 🙂

Resultado final

Sem delongas, o resultado final do script segue logo abaixo. Você pode customizá-lo para ficar do seu gosto: mudar cores, transparência, tema, posição e unidade de escala e norte.

Gostou do resultado? O que você mudaria nesse mapa? Boa sorte com seus futuros mapas e até o próximo post!

 

Agradecimentos

Agradeço ao Marco, pois o último mapa que construímos para um artigo me levou a escrever esse post. Agradeço também ao ​Oswaldo Baquero pela construção do pacote ggsn para R. E agradeço ao Miltinho por me introduzir ao mundo dos mapas.

Observação

Quem estiver usando Linux como sistema operacional e tiver problemas para executar os comandos, por favor cheque os comentários desse post.

(Fonte da image destacada)