Script para construir um banco com as cidades e as mesoregioes do parana a partir de um shapefile e um csv (arquivos já incluídos no pacote aRT). São criados vários temas com diferentes opções de visualização/fatiamento.
require(aRT) con=openConn() con if(any(showDbs(con)=="meso")) deleteDb(con, "meso", force=T)
criando o banco
db=createDb(con,db="meso")
importando o shapefile e informando sobre projecao
proj="+proj=poly" +lon_0=0.0 +lat_0=0.0 +x_0=0.0 +y_0=0.0 +a=6378160.0 +f=0.003353"
criando layer a partir do shapefile (diretamente s/ impotar p/ R)
l=createLayer(db, l="cidades", proj=proj) getProj(l) addShape(l, tab="tab", file=system.file("meso/parana_pol.shp",package="aRT"), id="CODIGO")
visualizando diretamente do banco
plot(l) ## trazendo os polygonos dos municípios para o R <code R> lpol <- getPolygons(l) aRTplot(lpol)
abre a tabela do shapefile e importa dados para o R
tab <- openTable(l,"tab") tdata <- getData(tab) head(tdata)
Criando os polígonos das meso-regioes pela união dos municípios
mesogeom <- list() mesonomes <- unique(tdata$CODIREGIAO) ## visualize a janela gráfica enquanto roda o comando a seguir for(i in 1:length(mesonomes)) { #selecionando municipios dentro da regiao i id=tdata$CODIGO[which(tdata$CODIREGIAO == mesonomes[i])] cat(paste("Processando mesoregiao",mesonomes[i], "com",length(id), "cidades\n")) #unindo os poligonos dos municipios para formar pol. da regiao set=getSetOperation(l,"union",id=id) set@polygons[[1]]@ID = paste(mesonomes[i]) #adicionando ao objeto (list) das regioes mesogeom[i] = set@polygons aRTplot(set, col=terrain.colors(length(mesonomes)+1)[i+1], add=T, lwd=2) } aRTplot(lpol, add=T) ## veja o gráfico <code> ## colocando no formato do "sp" <code R> res = SpatialPolygons(mesogeom, 1:length(mesogeom))
Insere os dados das meso-regioes no banco em um novo layer
l2=createLayer(db,l="mesoreg", proj=proj) addPolygons(l2, res)
importa a tabela estatica de um csv
mr <- read.csv(system.file("meso/parana_regions.csv",package="aRT"), sep=";") tabmeso <- importTable(l2, "tabmeso", id="CODIGO", data=mr)
Checando o status do Banco
bd <code> Criando uma nova coluna (variável) e colocando no banco:\\ Cria um dataframe com duas colunas: o CODIGO e a soma das producoes de QPM entre 2000 e 2003 <code R> somaQPM=data.frame(CODIGO=data$CODIGO, somaQPM= with(data, QPM_2000 + QPM_2001 + QPM_2002 + QPM_2003)) <code> Carrega os dados da tabela no R <code R> data=getData(tabmeso)
cria um dataframe com duas colunas: o CODIGO e a soma das producoes de QPM entre 2000 e 2003
somaQPM=data.frame(CODIGO=data$CODIGO, somaQPM=as.integer(data$QPM_2000)+ as.integer(data$QPM_2001)+ as.integer(data$QPM_2002)+ as.integer(data$QPM_2003))
Para cada coluna que nao existe ele cria antes de atualizar, e o codigo é utilizado para manter as referencias corretas as colunas do banco.
updateColumns(tabmeso, somaQPM)
Gerar temas coloridos para TV
v = visualPolygons(color="blue", transp=80, ccolor="black") t = createTheme(l2, "mesoregioes",view="view") setVisual(t, v)
abrir e ver no TV !!!! modificando a visualizacao das cores (vai ser encapsulado em funcoes)
v@ccolor = col2rgb("#F0F0F0") v@transp=as.integer(50) t = createTheme(l, "cidades",view="view") setVisual(t, v) v1=visualRaster()
Outra visualização (outra view) cria os temas para os cinco atributos, quatro producoes e a soma, coloridos. Estes temas sao criados a partir do layer de mesoregioes <code R> for(i in c("QPM_2000", "QPM_2001","QPM_2002","QPM_2003","somaQPM")) { cat(paste("Creating theme meso", i, "\n", sep="")) t=createTheme(l2, tname=paste("meso", i, sep=""), view = "vMeso") setVisual(t, v1, att=i) } <code>
Diferentes fatiamentos <code R> for(i in c("e","q","s")) {
cat(paste("Creating theme ", "QPM_2003", i, "\n", sep="")) t=createTheme(l, paste("cidQPM_2003", i, sep=""), view = "vFatia") ##mode abaixo define o criterio de criacao da legend (padrao TL) setVisual(t, v1, att="QPM_2003", mode=i) #setVisible(t, FALSE)
} <code>
## recaregue e visualize no TV!!!!
## se quiser apagar uma vista…. <code R> showViews(db) ## deleteView(db, "vMeso") <code>