<div dir="ltr"><div class="gmail_default" style=""><div class="gmail_default"><font face="verdana, sans-serif">Felipe, boa tarde!</font></div><div class="gmail_default"><font face="verdana, sans-serif"><br></font></div><div class="gmail_default"><font face="verdana, sans-serif">Segue uma sugestão, tentando preservar o máximo do seu código original. Alterações indicadas no script.</font></div><div class="gmail_default"><font face="verdana, sans-serif"><br></font></div><div class="gmail_default"><font face="verdana, sans-serif">Para diminuir o efeito de "serrilhado" nas figuras, aumente a resolução da grade.</font></div><div class="gmail_default"><font face="verdana, sans-serif"><br></font></div><div class="gmail_default"><span style="font-family:monospace,monospace">### <code r></span><br></div><div class="gmail_default"><span style="font-family:monospace,monospace"><br></span></div><div class="gmail_default">"<a href="http://www.leg.ufpr.br/lib/exe/fetch.php/projetos:arquivos_usados_no_tcc.zip" rel="noreferrer" target="_blank" style="font-size:12.8px">http://www.leg.ufpr.br/lib/exe/fetch.php/projetos:arquivos_usados_no_tcc.zip</a>" # Dados<span style="font-family:monospace,monospace"><br></span></div><div class="gmail_default"><br></div><div class="gmail_default"><font face="monospace, monospace"># ==== Carregando pacotes necessários ================ </font></div><div class="gmail_default"><font face="monospace, monospace">sapply(c("sp", "maptools", "spdep", "mgcv", "rgdal"), require, char=T)</font></div><div class="gmail_default"><font face="monospace, monospace">setwd("D:/Temp/TCC")</font></div><div class="gmail_default"><font face="monospace, monospace"><br></font></div><div class="gmail_default"><font face="monospace, monospace">##### lendo dados</font></div><div class="gmail_default"><font face="monospace, monospace">cwb<-readOGR("div_municipal.shp",layer="div_municipal") </font></div><div class="gmail_default"><font face="monospace, monospace">neo <- read.table('neo_contr1.txt',h=T) </font></div><div class="gmail_default"><font face="monospace, monospace"><br></font></div><div class="gmail_default"><font face="monospace, monospace">### Plotando o mapa de Curitiba </font></div><div class="gmail_default"><font face="monospace, monospace">plot(cwb) </font></div><div class="gmail_default"><font face="monospace, monospace">bbox(cwb) </font></div><div class="gmail_default"><font face="monospace, monospace">### Inserindo os controles </font></div><div class="gmail_default"><font face="monospace, monospace">points(neo$XCOORD[neo$Y==0],neo$YCOORD[neo$Y==0],col='darkgreen',cex = 0.5) </font></div><div class="gmail_default"><font face="monospace, monospace">### Inserindo os caso </font></div><div class="gmail_default"><font face="monospace, monospace">points(neo$XCOORD[neo$Y==1],neo$YCOORD[neo$Y==1],col='red',cex = 0.5) </font></div><div class="gmail_default"><font face="monospace, monospace"><br></font></div><div class="gmail_default"><font face="monospace, monospace">###### Ajustando modelo GAM </font></div><div class="gmail_default"><font face="monospace, monospace">models <- gam(neo$Y ~ +s(XCOORD,YCOORD, k=10, bs="tp"),family=binomial,data=neo) </font></div><div class="gmail_default"><font face="monospace, monospace">summary(models) </font></div><div class="gmail_default"><font face="monospace, monospace"><br></font></div><div class="gmail_default"><font face="monospace, monospace">## Fazendo a grid </font></div><div class="gmail_default"><font face="monospace, monospace">xLim <- range(pretty(bbox(cwb)[1,])) ### Adicionado</font></div><div class="gmail_default"><font face="monospace, monospace">yLim <- range(pretty(bbox(cwb)[2,])) ### Adicionado</font></div><div class="gmail_default"><font face="monospace, monospace">gx <- seq(xLim[1],  xLim[2], 500); nx <- length(gx); nx ### Modificado</font></div><div class="gmail_default"><font face="monospace, monospace">gy <- seq(yLim[1],  yLim[2], 500); ny <- length(gy); ny ### Modificado</font></div><div class="gmail_default"><font face="monospace, monospace">gr <- expand.grid(gx, gy) </font></div><div class="gmail_default"><font face="monospace, monospace">XCOORD <- gr$Var1 </font></div><div class="gmail_default"><font face="monospace, monospace">YCOORD <- gr$Var2 </font></div><div class="gmail_default"><font face="monospace, monospace"><br></font></div><div class="gmail_default"><font face="monospace, monospace">### criando data frame para predição </font></div><div class="gmail_default"><font face="monospace, monospace">ndados <- data.frame(XCOORD, YCOORD); head(ndados) </font></div><div class="gmail_default"><font face="monospace, monospace">pred   <- predict(models,ndados,type="response",se.fit=TRUE) </font></div><div class="gmail_default"><font face="monospace, monospace">pred   <- data.frame(pred) </font></div><div class="gmail_default"><font face="monospace, monospace">plot(cwb) </font></div><div class="gmail_default"><font face="monospace, monospace"><br></font></div><div class="gmail_default"><font face="monospace, monospace">## image </font></div><div class="gmail_default"><font face="monospace, monospace">image(gx,gy,matrix(pred$fit,ncol=ny,nrow=nx), asp=T, col=heat.colors(3,alpha=.5), ### Modificado</font></div><div class="gmail_default"><font face="monospace, monospace">       main="Thin Plate Splines com 128 pontos", </font></div><div class="gmail_default"><font face="monospace, monospace">       xlab="Coordenada X", ylab="Coordenada Y" ) </font></div><div class="gmail_default"><font face="monospace, monospace"><br></font></div><div class="gmail_default"><font face="monospace, monospace">plot(cwb, add=T) </font></div><div class="gmail_default"><font face="monospace, monospace">#points(neo$XCOORD,neo$YCOORD,pch=19,cex=0.5,col="white") </font></div><div class="gmail_default"><font face="monospace, monospace"><br></font></div><div class="gmail_default"><font face="monospace, monospace">##vis.gam </font></div><div class="gmail_default"><font face="monospace, monospace">vis.gam(models,plot.type="contour",view=c('XCOORD','YCOORD'),asp=1,ylim=c(7160000,7200000), </font></div><div class="gmail_default"><font face="monospace, monospace">main="Usando Thin Plate", </font></div><div class="gmail_default"><font face="monospace, monospace">color='cm') </font></div><div class="gmail_default"><font face="monospace, monospace">plot(cwb,add=T) </font></div><div class="gmail_default"><font face="monospace, monospace"><br></font></div><div class="gmail_default"><font face="monospace, monospace">### Adicionado ###############################################################</font></div><div class="gmail_default"><font face="monospace, monospace">### Criando uma Máscara - pontos da grade vs. polígono do município</font></div><div class="gmail_default"><font face="monospace, monospace">gr.sp   <- gr; coordinates(gr.sp) <- ~Var1+Var2</font></div><div class="gmail_default"><font face="monospace, monospace">cwb.pol <- SpatialPolygons(list(Polygons(list(Polygon(cwb@lines[[1]]@Lines[[1]]@coords)), "CWB")))</font></div><div class="gmail_default"><font face="monospace, monospace">mask    <- over(gr.sp, cwb.pol)</font></div><div class="gmail_default"><font face="monospace, monospace">image(gx, gy, matrix(mask, nrow=nx, ncol=ny), asp=T, axes=T, add=F)</font></div><div class="gmail_default"><font face="monospace, monospace"><br></font></div><div class="gmail_default"><font face="monospace, monospace">### Usando a máscara</font></div><div class="gmail_default"><font face="monospace, monospace">masked <- matrix(pred$fit, nrow=nx, ncol=ny)* matrix(mask, nrow=nx, ncol=ny)</font></div><div class="gmail_default"><font face="monospace, monospace"><br></font></div><div class="gmail_default"><font face="monospace, monospace">par(mfrow=c(1,2))</font></div><div class="gmail_default"><font face="monospace, monospace"><br></font></div><div class="gmail_default"><font face="monospace, monospace">image(gx, gy, masked, asp=T, col=heat.colors(3,alpha=.5), ### Modificado</font></div><div class="gmail_default"><font face="monospace, monospace">       main="Thin Plate Splines com 128 pontos", </font></div><div class="gmail_default"><font face="monospace, monospace">       xlab="Coordenada X", ylab="Coordenada Y" ); plot(cwb, add=T) </font></div><div class="gmail_default"><font face="monospace, monospace"><br></font></div><div class="gmail_default"><font face="monospace, monospace">contour(gx, gy, masked, asp=T, col=3, ### Modificado</font></div><div class="gmail_default"><font face="monospace, monospace">       main="Thin Plate Splines com 128 pontos", </font></div><div class="gmail_default"><font face="monospace, monospace">       xlab="Coordenada X", ylab="Coordenada Y" ); plot(cwb, add=T) </font></div><div class="gmail_default"><font face="monospace, monospace">### </code></font></div><div class="gmail_default"><font face="verdana, sans-serif"><br></font></div><div class="gmail_default"><font face="verdana, sans-serif"><img src="cid:ii_ilb563ud0_15338772b0da3433" width="508" height="508"><br>​<br></font></div><div class="gmail_default"><font face="verdana, sans-serif"><br></font></div><div class="gmail_default"><font face="verdana, sans-serif"><br></font></div></div><div class="gmail_extra"><br clear="all"><div><div class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><font face="arial, helvetica, sans-serif" style="font-size:small"><div style="font-family:'trebuchet ms',sans-serif;display:inline">​</div></font><font face="arial, helvetica, sans-serif" style="font-size:12.8px"><span style="font-size:small">================================================</span><br>Éder Comunello</font><div style="font-size:12.8px"><span style="font-family:arial,helvetica,sans-serif">Agronomist (UEM), MSc in Environ. Sciences (UEM)</span><font face="arial, helvetica, sans-serif"><br></font></div><div style="font-size:12.8px"><div style="font-size:small"><span style="font-family:arial,helvetica,sans-serif">DSc in Agricultural Systems Engineering (USP/Esalq)</span><br></div><div><span style="font-size:small">Brazilian Agricultural Research Corporation (</span><font face="arial, helvetica, sans-serif" style="font-size:small">Embrapa)</font></div><div><div><font face="arial, helvetica, sans-serif" style="font-size:small">Dourados, MS, Brazil </font><font face="arial, helvetica, sans-serif" size="2">|<O>|</font></div><div style="font-size:small"><span style="font-family:arial,helvetica,sans-serif">================================================</span><br></div></div></div><div style="font-size:small"><font face="arial, helvetica, sans-serif">GEO, </font>-22.2752, -54.8182, 408m<span style="font-family:arial,helvetica,sans-serif"><br></span></div><div style="font-size:12.8px"><span style="font-size:small;font-family:arial,helvetica,sans-serif">UTC-04:00 / DST: UTC-03:00</span></div><div><div><div><br></div><div><br></div></div><div style="font-size:small"><br></div></div></div></div></div></div></div></div></div></div></div></div></div>