Seu exemplo continua nao-reproduzivel, pois o arquivo nasc_vivo_7mais_cons_pre_natal.csv existe apenas para voce... (no historico da lista ha' sugestoes minhas para criacao de exemplos reproduziveis) Mas tenho a impressao de que seu problema pode ser resolvido apenas fixando manualmente os limites do eixo X... no comando no qual vc tem: plot(mapa, col = colors[findInterval(dados$ano_2006, brks,all.inside=TRUE)], axes = TRUE) imagino que algo como: plot(mapa, col = colors[findInterval(dados$ano_2006, brks,all.inside=TRUE)], axes = TRUE, xlim=c(-40, -33)) ja' possa comecar a dar os resultados que vc espera obter... A outra alternativa e' mover a rosa dos ventos mais para a esquerda: rosa_dos_ventos(loc=c(-35.5,-7.9),size=0.15,cex=0.6, bearing=0, cols=1) benilton 2012/3/20 Julianna Trindade <julianna@jubalitpb.com>:
O código está abaixo. os anexos mostram como ele deveria ser exportado e como está na realidade.
####################################################### # COSTRUINDO MAPAS ####################################################### ##Code (Comments are preceded by ##) ## Load required packages
# Limpando memsria. rm(list=ls(all=TRUE))
library(maps) library(sp) library(maptools) library(spdep) library(classInt) library(RColorBrewer) library(shape) library(maps) library(SDMTools) library(latticeExtra) library(ReadImages) library(png)
rosa_dos_ventos = function(loc=c(-34.5, -8), size = 0.15, bearing=0, cols=1, cex=0.6,...) { cols <- rep(c("white","black"),8) radii <- rep(size/c(1,4,2,4),4) x <- radii[(0:15)+1]*cos((0:15)*pi/8+bearing)+loc[1] y <- radii[(0:15)+1]*sin((0:15)*pi/8+bearing)+loc[2] for (i in 1:15) { x1 <- c(x[i],x[i+1],loc[1]) y1 <- c(y[i],y[i+1],loc[2]) polygon(x1,y1,col=cols[i]) } polygon(c(x[16],x[1],loc[1]),c(y[16],y[1],loc[2]),col=cols[16]) b <- c("O","S","L","N") for (i in 0:3) text((size+par("cxy")[1]-0.55)*cos(bearing+i*pi/2)+loc[1],
(size+par("cxy")[2]-0.60)*sin(bearing+i*pi/2)+loc[2],b[i+1],cex=cex) }
#Lendo dados dados = read.table("c:\\Users\\postgres\\nasc_vivo_7mais_cons_pre_natal.csv", sep = ";", header = FALSE)
colnames(dados) = c("codigos", "ano_2006","ano_2007","ano_2008", "ano_2009","ano_2010","ano_2011","ano_2012")
# Lendo mapa. mapa = readShapePoly("c:/users/Julianna/Documents/SecSaude/saudemapa.shp") proj4string(mapa) <- CRS("+init=epsg:4291") #scaleXY(mapa, 1836656)
# arquivo png img <- readPNG("c:/users/Julianna/Documents/SecSaude/logo_pb.png") # seu arquivo png aqui set.seed(123) par(xpd = TRUE)
colors <- brewer.pal(6, "YlOrRd") brks <- classIntervals(dados$ano_2008, n = 6,style="fixed", fixedBreaks=c(0, 20, 40, 60, 80, 100)) brks<- brks$brks
ID = match(substr(mapa$CODIGO,1,6),dados$codigos) dados = dados[ID,]
png(filename="c:/users/Julianna/Documents/SecSaude/produtos_nascvivos/2006.png", height=768, width=1320, bg="white")
plot(mapa, col = colors[findInterval(dados$ano_2006, brks,all.inside=TRUE)], axes = TRUE)
invisible(text(getSpPPolygonsLabptSlots(mapa), labels=as.character(mapa$NOMEMUN_1), cex=0.45))
title(main = "NASCIDOS VIVOS \nMaes com 7 ou mais consultas pre-natal", font.main = 2, cex.main = 1.5, xlab = "Longitude", ylab = "Latitude")
text(-35.7, -6.2, "ANO: 2006", cex = 1)
legend(-39.6,-7.5, c("0", "]0,20]", "]20,40]", "]40,60]", "]60,80]","Acima de 80"), fill = colors, bty="n", x.intersp = .5, y.intersp = .6, cex = 0.8)
rosa_dos_ventos(loc=c(-34.5,-7.9),size=0.15,cex=0.6, bearing=0, cols=1)
SpatialPolygonsRescale(layout.scale.bar(), offset = c(-35,-8.25), scale = 1, fill = c("transparent", "black"), plot.grid= FALSE)
text(-34.5, -8.3, "110.3 KM", cex= 0.8)
rasterImage(img, xleft = -39.5, ybottom = -5.9, xright = -38.7 , ytop = -5.6)
dev.off()
Em 19 de março de 2012 18:04, Benilton Carvalho <beniltoncarvalho@gmail.com> escreveu:
Depois de um exemplo reproduzivel, poderemos ajudar melhor.
b
2012/3/19 Julianna Trindade <julianna@jubalitpb.com>:
É verdade...eu alterei mas os elementos do plot tais como logo, legenda, texto, saem da dimensão especificada. Alguma sugestão??
Em 19 de março de 2012 17:32, Benilton Carvalho <beniltoncarvalho@gmail.com> escreveu:
pra "paisagem", vc nao deveria usar "width" maior que "height"? mas pode ser q eu nao tenha entendido o q vc quis dizer...
2012/3/19 Julianna Trindade <julianna@jubalitpb.com>:
Pessoal,
gostaria que meu plot fosse automaticamente transformado em um arquivo png em paisagem. Estou
usando png(filename="c:/users/Julianna/Documents/SecSaude/produtos_nascvivos/2006.png", height=1360, width=768, bg="white") e está exportando em retrato.
Alguém, por favor, pode me ajudar??
-- Julianna Trindade
"Exalta-te, Senhor, na Tua força! Nós cantaremos e louvaremos o Teu poder."
Salmos 21; 13
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.
R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.
-- Julianna Trindade
"Exalta-te, Senhor, na Tua força! Nós cantaremos e louvaremos o Teu poder."
Salmos 21; 13
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.
R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.
-- Julianna Trindade
"Exalta-te, Senhor, na Tua força! Nós cantaremos e louvaremos o Teu poder."
Salmos 21; 13
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.