<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.0 Transitional//EN">
<HTML><HEAD>
<META content="text/html; charset=iso-8859-1" http-equiv=Content-Type>
<META name=GENERATOR content="MSHTML 8.00.6001.19190">
<STYLE></STYLE>
</HEAD>
<BODY bgColor=#ffffff>
<DIV><FONT face=Arial></FONT> </DIV>
<DIV><FONT color=#0000ff face=Arial>Tentei rodar de novo, mas deu dois
problemas:</FONT></DIV>
<DIV><FONT size=2 face=Arial></FONT> </DIV>
<DIV><FONT size=2
face=Arial>library(maptools)
<BR>Carregando pacotes exigidos: foreign<BR>Carregando pacotes exigidos:
lattice<BR>Checking rgeos availability:
FALSE<BR> Note: when rgeos is not
available, polygon geometry computations in maptools
depend on gpclib,<BR> which has a
restricted licence. It is disabled by
default;<BR> to enable gpclib, type
gpclibPermit()</FONT></DIV>
<DIV><FONT size=2 face=Arial></FONT> </DIV>
<DIV><FONT size=2 face=Arial><FONT color=#0000ff size=3>e depois</FONT></DIV>
<DIV><BR></DIV></FONT>
<DIV><FONT face=Arial><FONT size=2>mapa = readShapePoly("saudemapa.shp")<BR>Erro
em getinfo.shape(filen) : Error opening SHP file<BR></FONT></FONT></DIV>
<DIV><FONT color=#0000ff face=Arial>E ainda faltou disponibilizar o
arquivo</FONT></DIV>
<DIV><FONT face=Arial><FONT size=2><FONT size=3
face="Times New Roman">logo_pb.png</FONT></DIV></FONT></FONT>
<BLOCKQUOTE
style="BORDER-LEFT: #000000 2px solid; PADDING-LEFT: 5px; PADDING-RIGHT: 0px; MARGIN-LEFT: 5px; MARGIN-RIGHT: 0px">
<DIV style="FONT: 10pt arial"> <BR></DIV>
<DIV class=gmail_quote>Em 1 de abril de 2012 21:05, Mauro Sznelwar <SPAN
dir=ltr><<A
href="mailto:sznelwar@uol.com.br">sznelwar@uol.com.br</A>></SPAN>
escreveuo:<BR>
<BLOCKQUOTE
style="BORDER-LEFT: #ccc 1px solid; MARGIN: 0px 0px 0px 0.8ex; PADDING-LEFT: 1ex"
class=gmail_quote><U></U>
<DIV bgcolor="#ffffff">
<DIV><FONT color=#0000ff face=Arial>Estava tentando rodar este código, e vi
que falta este
arquivoc:/users/Julianna/Documents/SecSaude/saudemapa.shp</FONT></DIV>
<DIV><FONT color=#0000ff face=Arial>Daria para enviar?.</FONT></DIV>
<DIV>
<DIV></DIV>
<DIV class=h5>
<BLOCKQUOTE
style="BORDER-LEFT: #000000 2px solid; PADDING-LEFT: 5px; PADDING-RIGHT: 0px; MARGIN-LEFT: 5px; MARGIN-RIGHT: 0px">
<DIV style="FONT: 10pt arial"> </DIV>
<DIV>Mensagem abaixo em nome da Julianna... (dados disponiveis via
datafilehost) </DIV>
<DIV><BR></DIV>
<DIV>====</DIV>
<DIV><BR></DIV>
<DIV>
<DIV>E quando fui conferir os breaks, os intervalos estavam certos,mas o
mapa resultante estava errado.</DIV>
<DIV>Porque quando eu não importo o csv e uso o dataframe resultante da
função, os breaks ficam ok(da maneira que eu especifico em FixedBreaks) e,
do contrario, não???</DIV>
<DIV><BR></DIV>
<DIV>segue o código:</DIV>
<DIV><BR></DIV>
<DIV>#######################################################</DIV>
<DIV># CONSTRUINDO MAPAS</DIV>
<DIV>#######################################################</DIV>
<DIV>##Code (Comments are preceded by ##)</DIV>
<DIV>## Load required packages</DIV>
<DIV><BR></DIV>
<DIV># Limpando memoria.</DIV>
<DIV>rm(list=ls(all=TRUE))</DIV>
<DIV><BR></DIV>
<DIV>library(maps)</DIV>
<DIV>library(sp)</DIV>
<DIV>library(maptools) </DIV>
<DIV>library(spdep) </DIV>
<DIV>library(classInt) </DIV>
<DIV>library(RColorBrewer) </DIV>
<DIV>library(shape)</DIV>
<DIV>library(SDMTools)</DIV>
<DIV>library(latticeExtra)</DIV>
<DIV>library(ReadImages)</DIV>
<DIV>library(png)</DIV>
<DIV><BR></DIV>
<DIV>rosa_dos_ventos = function(loc=c(-34.5, -8), size = 0.15, bearing=0,
cols=1, cex=0.6,...)</DIV>
<DIV>{</DIV>
<DIV> cols <- rep(c("white","black"),8)</DIV>
<DIV> radii <- rep(size/c(1,4,2,4),4)</DIV>
<DIV> x <- radii[(0:15)+1]*cos((0:15)*pi/8+bearing)+loc[1]</DIV>
<DIV> y <- radii[(0:15)+1]*sin((0:15)*pi/8+bearing)+loc[2]</DIV>
<DIV> for (i in 1:15)</DIV>
<DIV> {</DIV>
<DIV> x1 <- c(x[i],x[i+1],loc[1])</DIV>
<DIV> y1 <- c(y[i],y[i+1],loc[2])</DIV>
<DIV> polygon(x1,y1,col=cols[i])</DIV>
<DIV> }</DIV>
<DIV>
polygon(c(x[16],x[1],loc[1]),c(y[16],y[1],loc[2]),col=cols[16])</DIV>
<DIV> b <- c("O","S","L","N")</DIV>
<DIV> for (i in 0:3)
text((size+par("cxy")[1]-0.40)*cos(bearing+i*pi/2)+loc[1],</DIV>
<DIV>
(size+par("cxy")[2]-0.42)*sin(bearing+i*pi/2)+loc[2],b[i+1],cex=cex)</DIV>
<DIV>}</DIV>
<DIV><BR></DIV>
<DIV>## Lendo dados</DIV>
<DIV>## <A href="http://www.datafilehost.com/download-e5cb7e86.html"
target=_blank>http://www.datafilehost.com/download-e5cb7e86.html</A></DIV>
<DIV>dados = read.table("dadosJulianna.csv", sep = ",", header =
TRUE)</DIV>
<DIV><BR></DIV>
<DIV># Lendo mapa.</DIV>
<DIV>mapa =
readShapePoly("c:/users/Julianna/Documents/SecSaude/saudemapa.shp")</DIV>
<DIV>proj4string(mapa) <- CRS("+init=epsg:4291")</DIV>
<DIV>scaleXY(mapa, 1836656)</DIV>
<DIV><BR></DIV>
<DIV><BR></DIV>
<DIV># arquivo png</DIV>
<DIV>#img <-
readPNG("c:/users/Julianna/Documents/SecSaude/logo_pb.png") # seu arquivo
png aqui</DIV>
<DIV>#set.seed(123)</DIV>
<DIV>#par(xpd = TRUE)</DIV>
<DIV><BR></DIV>
<DIV><BR></DIV>
<DIV>colors <- brewer.pal(7, "Reds")</DIV>
<DIV>brks <- classIntervals(dados$ano_2006, n = 7,style="fixed",
intervalClosure="right",</DIV>
<DIV>
fixedBreaks=c(-2, -0.5, 0, 20, 40, 60, 80, 10000),
main="Fixed")</DIV>
<DIV>brks<- brks$brks</DIV>
<DIV><BR></DIV>
<DIV>ID =
match(substr(mapa$CODIGO,1,6),dados$t_indicadores_municipio.f_municipio)</DIV>
<DIV>dados = dados[ID,]</DIV>
<DIV><BR></DIV>
<DIV>png(filename="c:/mif2006.png", height=750, width=1500)</DIV>
<DIV><BR></DIV>
<DIV><BR></DIV>
<DIV>plot(mapa, col = colors[findInterval(dados$ano_2006,
brks,all.inside=TRUE)], axes = TRUE, cex = 1.5, mar=5)</DIV>
<DIV>##, xlim=c(-39.7, -34.7), ylim=c(-8.5, -6)</DIV>
<DIV><BR></DIV>
<DIV>##invisible(text(getSpPPolygonsLabptSlots(mapa),
labels=as.character(mapa$NOMEMUN_1), cex=0.65))</DIV>
<DIV><BR></DIV>
<DIV>title(main = "PORCENTAGEM DE Ã BITOS INVESTIGADOS\ntipo: mulheres em
idade fértil",</DIV>
<DIV> font.main = 2, cex.main = 1.8)</DIV>
<DIV><BR></DIV>
<DIV>text(-35.7, -6.2, "ANO: 2006", cex = 1.8)</DIV>
<DIV><BR></DIV>
<DIV>legend(-39.6,-7.6, c("Sem informação", "0", "]0,20]", "]20,40]",
"]40,60]", "]60,80]","Acima de 80"),</DIV>
<DIV> fill = colors, bty="n", x.intersp = 1.1,
y.intersp = 1.1, cex = 1.6)</DIV>
<DIV><BR></DIV>
<DIV>rosa_dos_ventos(loc=c(-34.5,-8),size=0.17,cex=1.20, bearing=0,
cols=1)</DIV>
<DIV><BR></DIV>
<DIV>SpatialPolygonsRescale(layout.scale.bar(), offset =
c(-35,-8.30),</DIV>
<DIV>
scale = 1, fill = c("transparent", "black"),</DIV>
<DIV>
plot.grid= FALSE)</DIV>
<DIV><BR></DIV>
<DIV>text(-34.5, -8.35, "110.3 KM", cex= 1.3)</DIV>
<DIV><BR></DIV>
<DIV>#rasterImage(img, xleft = -39.5, ybottom = -6.25, xright = -38.9 ,
ytop = -6)</DIV>
<DIV><BR></DIV>
<DIV>dev.off()</DIV>
<DIV><BR></DIV></DIV>
<P></P>
<P><FONT face=Arial></FONT></P>
<P></P></BLOCKQUOTE></DIV></DIV></DIV><BR></BLOCKQUOTE></DIV></BLOCKQUOTE></BODY></HTML>