
Boa tarde a todos. Tendo de preparar um mapa com duas áreas, como consigo separar as duas áreas sem que fiquem ligadas por uma linha limite=read.table("rec.txt",header=T) limite plot(limite,type="l") anexei o arquivo rec.txt A finalidade é fazer um recorte em um mapa de krigagem grato -- Hélio Gallo Rocha IFSULDEMINAS - Câmpus Muzambinho

você pode colocar NA nos pontos que estão interligados suponha o seguinte: duas áreas limitadas por quadrados a=c(1,1,3,3,1,4,4,5,5,4) b=c(2,4,4,2,2,6,7,7,6,6) f=matrix(c(a,b),ncol=2) plot(f,type='l') # caso vc queira remover a linha basta adicionar um NA a=c(a[1:5],NA,a[6:10]) b=c(b[1:5],NA,b[6:10]) f=data.frame(a,b) plot(f,type='l') Seria isso?

Caso você quira fazer um recorte no mapa, faça uma máscara Criar uma matriz de zero e um (onde zero é o que você quer omitir e 1 os seus dados), através de uma equação da forma geométrica que você deseja multiplique essa máscara pelos seus dados, e plote novamente segue um exemplo abaixo de função que remove as bordas de um circulo mask <- function (dados){ # tamanho da coluna Y=c(1:dim(dados)[1]) # calcular o raio r=floor(length(Y)/2) # centro do círculo cp = round(length(Y)/2,0) # matriz com zeros TT=matrix(0,nrow=dim(dados)[1],ncol=dim(dados)[2]) # fazer uma máscara de um e zero for(X in c(1:dim(dados)[2])){ # calcula o raio de todos os Y e X com dados R=sqrt((X-cp)^2+(Y-cp)^2) #o maior valor do raio do circulo zeros=which(R>=r) # remove os pontos de fora do circulo TT[X,-zeros]=1 } return(TT) } aplicando aos meus dados fica dados.recortados<-dados*mask(dados) espero ter ajudado Tito Conte 2013/1/18 Tito Conte <tito.conte@gmail.com>
você pode colocar NA nos pontos que estão interligados
suponha o seguinte: duas áreas limitadas por quadrados
a=c(1,1,3,3,1,4,4,5,5,4) b=c(2,4,4,2,2,6,7,7,6,6) f=matrix(c(a,b),ncol=2) plot(f,type='l')
# caso vc queira remover a linha basta adicionar um NA
a=c(a[1:5],NA,a[6:10]) b=c(b[1:5],NA,b[6:10]) f=data.frame(a,b) plot(f,type='l')
Seria isso?

Valeu Tito A primeira opção mais simples e já resolveu Abraço Hélio Em 18 de janeiro de 2013 13:16, Tito Conte <tito.conte@gmail.com> escreveu:
você pode colocar NA nos pontos que estão interligados
suponha o seguinte: duas áreas limitadas por quadrados
a=c(1,1,3,3,1,4,4,5,5,4) b=c(2,4,4,2,2,6,7,7,6,6) f=matrix(c(a,b),ncol=2) plot(f,type='l')
# caso vc queira remover a linha basta adicionar um NA
a=c(a[1:5],NA,a[6:10]) b=c(b[1:5],NA,b[6:10]) f=data.frame(a,b) plot(f,type='l')
Seria isso?
_______________________________________________ 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.
-- Hélio Gallo Rocha IFSULDEMINAS - Câmpus Muzambinho
participantes (2)
-
Hélio Gallo Rocha
-
Tito Conte