software:art:curso:meso

Não foi possível enviar o arquivo. Será algum problema com as permissões?

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

Diferenças

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

--- software:art:curso:meso [2007/12/12 01:20]
paulojus criada
+++ software:art:curso:meso [2007/12/12 09:28] (atual)
paulojus
@@ Linha 31: / Linha 31: @@
 visualizando diretamente do banco
 <code R>
-getPolygons(l)->lpol
+plot(l)
+## trazendo os polygonos dos municípios para o R
+<code R>
+lpol <- getPolygons(l)
 aRTplot(lpol)
 </code>
@@ Linha 37: / Linha 41: @@
 abre a tabela do shapefile e importa dados para o R
 <code R>
-tab = openTable(l,"tab")
+tab <- openTable(l,"tab")
 tdata <- getData(tab)
+head(tdata)
 </code>
-e calcula as meso-regioes
+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
+  #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
+  #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
+  #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)
+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, e cria temas para visualizacao
+Insere os dados das meso-regioes no banco em um novo layer
 <code R>
 l2=createLayer(db,l="mesoreg", proj=proj)
@@ Linha 72: / Linha 81: @@
 importa a tabela estatica de um csv
 <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>
-#############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>
 data=getData(tabmeso)
@@ Linha 92: / Linha 112: @@
 </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>
 updateColumns(tabmeso, somaQPM)
 </code>
-#############FIM NOVO ################################
-nova classe em desenvolvimento
+Gerar temas coloridos para TV
-para gerar temas coloridos para TV
 <code R>
 v = visualPolygons(color="blue", transp=80, ccolor="black")
@@ Linha 108: / Linha 125: @@
 </code>
+**//abrir e ver no TV !!!! //**
 modificando a visualizacao das cores (vai ser encapsulado em funcoes)
 <code R>
@@ Linha 116: / Linha 134: @@
 v1=visualRaster()
-aRTsilent(TRUE)
 </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>
-#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=""))
- 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)
-#}
+}
-</code>
+<code>
-#############FIM NOVO ################################
+Diferentes fatiamentos
 <code R>
 for(i in c("e","q","s"))
 {
   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)
- 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>
-#munic <- getGeometry(l)
+showViews(db)
-#is.list(munic)
+## deleteView(db, "vMeso")
-#names(munic)
+<code>
-#is.list(munic$polygons)
-#length(munic$polygons)
-#plot(munic)
-</code>

Diferenças

Navegação

Busca

Ferramentas

QR Code