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