Você está aqui: start » software » art » curso » meso

Diferenças

Aqui você vê as diferenças entre duas revisões dessa página.

Link para esta página de comparações

software:art:curso:meso [2007/12/12 01:20] – criada paulojussoftware:art:curso:meso [2007/12/12 09:28] (atual) – paulojus
Linha 31: Linha 31:
 visualizando diretamente do banco visualizando diretamente do banco
 <code R> <code R>
-getPolygons(l)->lpol+plot(l) 
 + 
 +## trazendo os polygonos dos municípios para o R 
 +<code R> 
 +lpol <- getPolygons(l)
 aRTplot(lpol) aRTplot(lpol)
 </code> </code>
Linha 37: Linha 41:
 abre a tabela do shapefile e importa dados para o R abre a tabela do shapefile e importa dados para o R
 <code R> <code R>
-tab = openTable(l,"tab")+tab <- openTable(l,"tab")
 tdata <- getData(tab) tdata <- getData(tab)
 +head(tdata)
 </code> </code>
  
-e calcula as meso-regioes+Criando os polígonos das meso-regioes pela união dos municípios
 <code R> <code R>
 mesogeom <- list() mesogeom <- list()
 mesonomes <- unique(tdata$CODIREGIAO) mesonomes <- unique(tdata$CODIREGIAO)
  
 +## visualize a janela gráfica enquanto roda o comando a seguir
 for(i in 1:length(mesonomes)) for(i in 1:length(mesonomes))
 { {
-  selecionando municipios dentro da regiao i+  #selecionando municipios dentro da regiao i
   id=tdata$CODIGO[which(tdata$CODIREGIAO == mesonomes[i])]   id=tdata$CODIGO[which(tdata$CODIREGIAO == mesonomes[i])]
   cat(paste("Processando mesoregiao",mesonomes[i], "com",length(id), "cidades\n"))   cat(paste("Processando mesoregiao",mesonomes[i], "com",length(id), "cidades\n"))
-  unindo os poligonos dos municipios para formar pol. da regiao+  #unindo os poligonos dos municipios para formar pol. da regiao
   set=getSetOperation(l,"union",id=id)   set=getSetOperation(l,"union",id=id)
   set@polygons[[1]]@ID = paste(mesonomes[i])   set@polygons[[1]]@ID = paste(mesonomes[i])
-  adicionando ao objeto (list) das regioes+  #adicionando ao objeto (list) das regioes
   mesogeom[i] = set@polygons   mesogeom[i] = set@polygons
   aRTplot(set, col=terrain.colors(length(mesonomes)+1)[i+1], add=T, lwd=2)   aRTplot(set, col=terrain.colors(length(mesonomes)+1)[i+1], add=T, lwd=2)
 } }
  
-aRTplot(lpol, add=T)+aRTplot(lpol, add=T) ## veja o gráfico 
 +<code> 
 +## colocando no formato do "sp" 
 + 
 +<code R>
 res = SpatialPolygons(mesogeom, 1:length(mesogeom)) res = SpatialPolygons(mesogeom, 1:length(mesogeom))
 </code> </code>
  
-insere os dados das meso-regioes, e cria temas para visualizacao +Insere os dados das meso-regioes no banco em um novo layer
 <code R> <code R>
 l2=createLayer(db,l="mesoreg", proj=proj) l2=createLayer(db,l="mesoreg", proj=proj)
Linha 72: Linha 81:
 importa a tabela estatica de um csv importa a tabela estatica de um csv
 <code R> <code R>
-mr=read.csv(system.file("meso/parana_regions.csv",package="aRT"), sep=";") +mr <- read.csv(system.file("meso/parana_regions.csv",package="aRT"), sep=";") 
-tabmeso=importTable(l2, "tabmeso", id="CODIGO", data=mr)+tabmeso <- importTable(l2, "tabmeso", id="CODIGO", data=mr)
 </code> </code>
  
-#############NOVO ################################ +Checando o status do Banco 
-carrega os dados da tabela no R+<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> <code R>
 data=getData(tabmeso) data=getData(tabmeso)
Linha 92: Linha 112:
 </code> </code>
  
-as colunas do banco. Para cada coluna que nao existe ele +Para cada coluna que nao existe ele cria antes de atualizar, e o codigo é utilizado para manter as referencias corretas 
-cria antes de atualizar, e o codigo é utilizado para +as colunas do banco. 
-manter as referencias corretas+
 <code R> <code R>
 updateColumns(tabmeso, somaQPM) updateColumns(tabmeso, somaQPM)
 </code> </code>
-#############FIM NOVO ################################ 
  
-nova classe em desenvolvimento +Gerar temas coloridos para TV
-para gerar temas coloridos para TV+
 <code R> <code R>
 v = visualPolygons(color="blue", transp=80, ccolor="black") v = visualPolygons(color="blue", transp=80, ccolor="black")
Linha 108: Linha 125:
 </code> </code>
  
 +**//abrir e ver no TV !!!! //**
 modificando a visualizacao das cores (vai ser encapsulado em funcoes) modificando a visualizacao das cores (vai ser encapsulado em funcoes)
 <code R> <code R>
Linha 116: Linha 134:
  
 v1=visualRaster() v1=visualRaster()
-aRTsilent(TRUE) 
 </code> </code>
  
-criando temas com as atributos dos minicipios 
-<code R> 
-#colnames = colnames(tdata) 
-#for(i in 1:length(colnames)) 
-#{ 
-# if(!(any(colnames[i] == c("CODIGO", "AREA_1", "PERIMETRO_", "NOME", "CODIREGIAO")))) 
-# { 
-# cat(paste("Creating theme", colnames[i], "\n")) 
-# t=createTheme(l, th=colnames[i], view = "view") 
-# setVisual(t, v1, colnames[i]) 
-# } 
-#} 
-</code> 
  
- +Outra visualização (outra view) cria os temas para os cinco atributos, quatro producoes e 
-#############NOVO ################################ +a soma, coloridos. Estes temas sao criados a partir do layer de mesoregioes
-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> <code R>
-#for(i in c("QPM_2000", "QPM_2001","QPM_2002","QPM_2003","somaQPM")) +for(i in c("QPM_2000", "QPM_2001","QPM_2002","QPM_2003","somaQPM")) 
-#{+{
  cat(paste("Creating theme meso", i, "\n", sep=""))  cat(paste("Creating theme meso", i, "\n", sep=""))
- t=createTheme(l2, th=paste("meso", i, sep=""), view = "view")+ t=createTheme(l2, tname=paste("meso", i, sep=""), view = "vMeso")
  setVisual(t, v1, att=i)  setVisual(t, v1, att=i)
-#} +} 
-</code> +<code> 
-#############FIM NOVO ################################+ 
 +Diferentes fatiamentos
 <code R> <code R>
 for(i in c("e","q","s")) for(i in c("e","q","s"))
 { {
   cat(paste("Creating theme ", "QPM_2003", i, "\n", sep=""))   cat(paste("Creating theme ", "QPM_2003", i, "\n", sep=""))
-  t=createTheme(l, paste("cidadesQPM_2003", i, sep=""), view = "view") +  t=createTheme(l, paste("cidQPM_2003", i, sep=""), view = "vFatia") 
-    #mode abaixo define o criterio de criacao da legend (padrao TL)+  ##mode abaixo define o criterio de criacao da legend (padrao TL)
   setVisual(t, v1, att="QPM_2003", mode=i)   setVisual(t, v1, att="QPM_2003", mode=i)
- setVisible(t, FALSE)+  #setVisible(t, FALSE)
 } }
-</code>+<code>
  
-temas com atributos das mesoregioes +## recaregue e visualize no TV!!!!
-<code R> +
-#colnames = colnames(mr) +
-#for(i in 1:length(colnames)) +
-#{ +
-# cat(paste("Creating theme", colnames[i], "\n")) +
-# if(!(any(colnames[i] == c("NOME", "CODIGO")))) +
-# { +
-# t=createTheme(l2, th=colnames[i], view = "view") +
-# setVisual(t, v1, colnames[i]) +
-# } +
-#} +
-</code>+
  
-operacoes em R+## se quiser apagar uma vista....
 <code R> <code R>
-#munic <- getGeometry(l) +showViews(db) 
-#is.list(munic) +## deleteView(db, "vMeso") 
-#names(munic) +<code>
-#is.list(munic$polygons) +
-#length(munic$polygons) +
-#plot(munic) +
-</code> +
  

QR Code
QR Code software:art:curso:meso (generated for current page)