====== Exemplo: combinando dados de municípios e mesoregiões ======

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.

<code R>
require(aRT)
con=openConn()
con

if(any(showDbs(con)=="meso"))
  deleteDb(con, "meso", force=T)
</code>

criando o banco
<code R>
db=createDb(con,db="meso")
</code>

importando o shapefile e informando sobre projecao
<code R>
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"
</code>

criando layer a partir do shapefile (diretamente s/ impotar p/ R)
<code 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")
</code>

visualizando diretamente do banco
<code R>
plot(l)

## trazendo os polygonos dos municípios para o R
<code R>
lpol <- getPolygons(l)
aRTplot(lpol)
</code>

abre a tabela do shapefile e importa dados para o R
<code R>
tab <- openTable(l,"tab")
tdata <- getData(tab)
head(tdata)
</code>

Criando os polígonos das meso-regioes pela união dos municípios
<code R>
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))
</code>

Insere os dados das meso-regioes no banco em um novo layer
<code R>
l2=createLayer(db,l="mesoreg", proj=proj)
addPolygons(l2, res)
</code>

importa a tabela estatica de um csv
<code R>
mr <- read.csv(system.file("meso/parana_regions.csv",package="aRT"), sep=";")
tabmeso <- importTable(l2, "tabmeso", id="CODIGO", data=mr)
</code>

Checando o status do Banco
<code R>
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)
</code>

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=as.integer(data$QPM_2000)+
				           as.integer(data$QPM_2001)+
						   as.integer(data$QPM_2002)+
						   as.integer(data$QPM_2003))
</code>

Para cada coluna que nao existe ele cria antes de atualizar, e o codigo é utilizado para manter as referencias corretas
as colunas do banco. 
<code R>
updateColumns(tabmeso, somaQPM)
</code>

Gerar temas coloridos para TV
<code R>
v = visualPolygons(color="blue", transp=80, ccolor="black")
t = createTheme(l2, "mesoregioes",view="view")
setVisual(t, v)
</code>

**//abrir e ver no TV !!!! //**
modificando a visualizacao das cores (vai ser encapsulado em funcoes)
<code R>
v@ccolor = col2rgb("#F0F0F0")
v@transp=as.integer(50)
t = createTheme(l, "cidades",view="view")
setVisual(t, v)

v1=visualRaster()
</code>


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>