levelplot e contourplot
Olá todos, Estou produzindo mapas de nível topográfico com o Lattice, usando as funções levelplot() e contourplot(). Depois de pronto, tento colocar pontos no mapa, usando points() e não funciona. Se alguém souber onde está o erro, agradeço qualquer ajudar. Att Carlos o código é o seguinte: elevation.loess = loess(z ~ x*y, data=data.s, degree = 2, span = .3, family="gaussian") elevation.fit = expand.grid(list(x = seq(65.86,71.6,0.01), y = seq(1.81,6.845,0.01))) z = predict(elevation.loess, newdata = elevation.fit) z1<-ifelse(z<0,0,z) height<-(elevation.fit$Height = as.numeric(z1)) library(lattice) levelplot(height~-x*-y, data = elevation.fit, xlab = "Latitude", ylab = "Longitude", main = "Surface elevation data", col.regions = terrain.colors(100),cuts=40,contour=T, plot.new=T, points=T) contourplot(height~-x*-y, data = elevation.fit, cuts = 20,labels = TRUE,col.regions = terrain.colors(100), contour = TRUE, pretty = TRUE,region = T) points(-70,-4,pch=20,col="blue") # aqui é onde não funciona
tenta points(c(-70,-4),pch=20,col="blue") Em 28 de outubro de 2012 09:08, D'Apolito <carlosdapolito@gmail.com>escreveu:
Olá todos,
Estou produzindo mapas de nível topográfico com o Lattice, usando as funções levelplot() e contourplot(). Depois de pronto, tento colocar pontos no mapa, usando points() e não funciona. Se alguém souber onde está o erro, agradeço qualquer ajudar.
Att Carlos
o código é o seguinte:
elevation.loess = loess(z ~ x*y, data=data.s, degree = 2, span = .3, family="gaussian")
elevation.fit = expand.grid(list(x = seq(65.86,71.6,0.01), y = seq(1.81,6.845,0.01)))
z = predict(elevation.loess, newdata = elevation.fit)
z1<-ifelse(z<0,0,z)
height<-(elevation.fit$Height = as.numeric(z1))
library(lattice)
levelplot(height~-x*-y, data = elevation.fit, xlab = "Latitude", ylab = "Longitude", main = "Surface elevation data", col.regions = terrain.colors(100),cuts=40,contour=T, plot.new=T, points=T)
contourplot(height~-x*-y, data = elevation.fit, cuts = 20,labels = TRUE,col.regions = terrain.colors(100), contour = TRUE, pretty = TRUE,region = T)
points(-70,-4,pch=20,col="blue") # aqui é onde não funciona
_______________________________________________ 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.
-- Natália da Silva Martins Bacharel em Estatística - Universidade Estadual de Maringá/ UEM Mestranda em Estatística e Experimentação Agronômica - ESALQ/ USP Contato: (19) 8306-4743
obrigado, Natália, mas ainda não deu: Erro em plot.xy(xy.coords(x, y), type = type, ...) : plot.new ainda não foi chamado Em 28 de outubro de 2012 08:15, Natalia Martins <nsmbarreto@gmail.com>escreveu:
tenta points(c(-70,-4),pch=20,col="blue")
Em 28 de outubro de 2012 09:08, D'Apolito <carlosdapolito@gmail.com>escreveu:
Olá todos,
Estou produzindo mapas de nível topográfico com o Lattice, usando as funções levelplot() e contourplot(). Depois de pronto, tento colocar pontos no mapa, usando points() e não funciona. Se alguém souber onde está o erro, agradeço qualquer ajudar.
Att Carlos
o código é o seguinte:
elevation.loess = loess(z ~ x*y, data=data.s, degree = 2, span = .3, family="gaussian")
elevation.fit = expand.grid(list(x = seq(65.86,71.6,0.01), y = seq(1.81,6.845,0.01)))
z = predict(elevation.loess, newdata = elevation.fit)
z1<-ifelse(z<0,0,z)
height<-(elevation.fit$Height = as.numeric(z1))
library(lattice)
levelplot(height~-x*-y, data = elevation.fit, xlab = "Latitude", ylab = "Longitude", main = "Surface elevation data", col.regions = terrain.colors(100),cuts=40,contour=T, plot.new=T, points=T)
contourplot(height~-x*-y, data = elevation.fit, cuts = 20,labels = TRUE,col.regions = terrain.colors(100), contour = TRUE, pretty = TRUE,region = T)
points(-70,-4,pch=20,col="blue") # aqui é onde não funciona
_______________________________________________ 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.
--
Natália da Silva Martins Bacharel em Estatística - Universidade Estadual de Maringá/ UEM Mestranda em Estatística e Experimentação Agronômica - ESALQ/ USP Contato: (19) 8306-4743
_______________________________________________ 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.
-- _____________________________________ *Carlos D'Apolito - PhD candidate - ID:1239480* Earth Sciences School of Geography & Environmental Sciences Aston Webb Building University of Birmingham Birmingham B15 2TT UK phone: 44 - 07 919 673 173 CXD280@bham.ac.uk CV: http://lattes.cnpq.br/4557754746424026
A lattice (xyplot, bwplot, levelplot) e graphics (plot, boxplot, image) são dispositivos gráficos diferentes. As funções gráficas para desenho na lattice são as mesmas só que precedidas do nome "panel.". Então para adicionar pontos, linhas, segmentos, setas você deve fazer levelplot(suas_configurações_aqui) trellis.focus("panel", 1, 1) panel.points(seus_pontos_aqui) trellis.unfocus() Consulte a documentação para detalhes. À disposição. Walmes. ========================================================================== Walmes Marques Zeviani LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W) Departamento de Estatística - Universidade Federal do Paraná fone: (+55) 41 3361 3573 VoIP: (3361 3600) 1053 1173 e-mail: walmes@ufpr.br skype: walmeszeviani twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ==========================================================================
Entendido. Muito obrigado Walmes!!!! Em 28 de outubro de 2012 10:16, Walmes Zeviani <walmeszeviani@gmail.com>escreveu:
A lattice (xyplot, bwplot, levelplot) e graphics (plot, boxplot, image) são dispositivos gráficos diferentes. As funções gráficas para desenho na lattice são as mesmas só que precedidas do nome "panel.". Então para adicionar pontos, linhas, segmentos, setas você deve fazer
levelplot(suas_configurações_aqui) trellis.focus("panel", 1, 1) panel.points(seus_pontos_aqui) trellis.unfocus()
Consulte a documentação para detalhes.
À disposição. Walmes.
========================================================================== Walmes Marques Zeviani LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W) Departamento de Estatística - Universidade Federal do Paraná fone: (+55) 41 3361 3573 VoIP: (3361 3600) 1053 1173 e-mail: walmes@ufpr.br skype: walmeszeviani twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ==========================================================================
_______________________________________________ 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.
-- _____________________________________ *Carlos D'Apolito - PhD candidate - ID:1239480* Earth Sciences School of Geography & Environmental Sciences Aston Webb Building University of Birmingham Birmingham B15 2TT UK phone: 44 - 07 919 673 173 CXD280@bham.ac.uk CV: http://lattes.cnpq.br/4557754746424026
participantes (3)
-
D'Apolito -
Natalia Martins -
Walmes Zeviani