ow Galera,Comecei a ver esse CMR contudo não consigo instalar o pacote rgdal.
Li nele que é necessário instalar um tal de Gdal. É isso mesmo? Vcs podem me
adiantar algo?
Agradeço novamente pelo auxílio na resolução do problema.
Arthur.
_______________________________________________________________________________________________________________________________________
Date: Tue, 30 Apr 2013 17:24:02 -0300
From: [email protected]
To: [email protected]
Subject: Re: [R-br] Krigagem
Feito, grato....
Pra mim tá resolvido
Vou postar mais uma dúvida neste capitulo mesmo, já que deva ter vários
interessados
Teria como fazer um clip para obter apenas os dados de uma área limite, mas que
não seja uma área regular.
Anexei a bordadura limite do exemplo do Arthur em latlong.
vai o CRM de conversão da borda de latlong pata lambert:
pb <- read.table("k:/2013/r/arthur/borda.txt", head = T,blank.lines.skip=T) ; pb
pllb <- SpatialPointsDataFrame(pb[,1:2], data=pb[3],proj4string=CRS("+proj=longlat
+datum=WGS84"));pllb
cllb <- data.frame(lon = round(mean(range(pb$long)),3), lat =
round(mean(range(pb$lat)),3));cllb #média
## projecao azimutal de área equivalente de Lambert (meridiano central em -60
para America do Sul)
projlaeab <- "+proj=laea +lat_0=-15.65 +lon_0=-45.135 +x_0=0 +y_0=0 +a=6370997
+b=6370997 +units=m +no_defs"; projlaeab
## dados reprojetados em laea, unidade em m
plaeab <- spTransform(x=pllb, CRSobj=CRS(projlaeab)) ; plaeab
tabela=edit(data.frame(plaeab))
write.table(tabela,file="k:\\2013/r/arthur/borda_lamb.txt",sep = '\t',col.names
= NA)
##para ler o arquivo em lambert na krigagem:
borda=read.table("borda_lamb.txt",header=T);borda
borda=cbind(borda$long,borda$lat);borda
Também fiz a krigagem usando latlong , ficou bastante parecido, com pequenas
diferenças, mas o melhor método de ajuste para latlong foi
da Verossimilhança e para lambert foi quadrados minimos ponderado
Hélio
Em 30 de abril de 2013 12:09, Jônatan [via R-br]
<[email protected]> escreveu:
Os coordenadas de referência estão definidas como x_0=0 +y_0=0, pq? Tem
que definir como cll (ver no ultimo email que
enviei).
2013/4/30 Hélio Gallo Rocha <[hidden email]>
Caro JônatanFiz a transformação de latlong para lambert como orientou,
fiz a krigagem( no caso o melhor modelo foi
por quadrados mininos ponderados).
Para facilitar a verificação,usei l=10 no loci.
Gerei um txt ,(em anexo) das coordenadas em lambert e o kc$predict para poder
voltar a latlong, como disse se quiser, mas
para deixar o CRM completo, dai a dificuladade, olha o CRM da volta:
out_krig=read.table("k:/2013/r/arthur/graficos/out_krig_lambert.txt", head =
T,blank.lines.skip=T) ; summary(out_krig)
;class(out_krig) ;out_krig#OK
out_krig_spdf <- SpatialPointsDataFrame(out_krig[,2:3],
data=out_krig[1],proj4string=CRS("+proj=laea + x_0=0 +y_0=0
+a=6370997 +b=6370997+units=m +no_defs"))# +datum=WGS84"))
out_krig_spdf
não deu certo...
Hélio
Em 28 de abril de 2013 21:42, Jônatan [via R-br] <[hidden email]> escreveu:
Hélio segue CMR da comparação do método que sugeri com outros (incluindo
dicas do Elias e do Paulo).
## comparacão da distancia entre 2 pontos (no dominio espacial dos dados) por
diferentes metodos
require(sp);require(rgdal);require(fields) # para comparacao
## Funcoes da dica do Elias
source("http://www.biostat.umn.edu/~sudiptob/Software/distonearth.R";)
## dados
p <- read.table("C:/Users/Jonatan/Downloads/Dados.txt", head =
T,blank.lines.skip=T)
pll <- SpatialPointsDataFrame(p[,1:2], data=p[3], proj4string=CRS("+proj=longlat
+datum=WGS84"))
## centro aproximado
cll <- data.frame(lon = round(mean(range(p$Longitude)),3), lat =
round(mean(range(p$Latitude)),3))
## adicionando ponto a uma distancia conhecida
loc <- rbind(cll, cll + 1)
locll <- SpatialPoints(loc, proj4string=CRS("+proj=longlat +datum=WGS84"))
## pontos e centro
plot(pll,col="gray")
points(cll,col=2,pch=20)
## indicando os 2 pontos
text(coordinates(locll),labels=rownames(loc))
## projecao azimutal de área equivalente de Lambert
projlaea0 <- "+proj=laea +lat_0=YY +lon_0=XX +x_0=0 +y_0=0 +a=6370997 +b=6370997
+units=m +no_defs"
projlaea <- gsub("XX", cll$lon, projlaea0)
projlaea <- gsub("YY", cll$lat, projlaea)
## dados reprojetados em laea, unidade em m
loclaea <- spTransform(x=locll, CRSobj=CRS(projlaea))
## distancia em km a partir das coordenadas em laea
res1 <- dist(coordinates(loclaea))/1000
## aprox grosseira, 1 grau ~ 111 km
res2 <- sqrt(111^2+111^2)
## distancia estimada pelo fields
res3 <- rdist.earth(loc, miles=F)
## Uma das Funcoes do Sudipto Banerjee
res4 <- rdistearth(loc)
## diferencas
res1 - rbind(res1, res2, unique(res3[res3>0]), unique(res4[res4>0]))
Quanto as duvidas sobre as propriedades da proj. LAEA:
http://en.wikipedia.org/wiki/Lambert_azimuthal_equal-area_projection
2013/4/28 Hélio Gallo Rocha <[hidden email]>
Caro Jônatan,
Nunca usei lambert
duas dúvidas, talvez a resposta da primeira já resolva:
no CRM: +lat_0=-15 +lon_0=-60, centro da projeção
quais os limites máximos para que não haja distorções? de forma rápida achei
que a escala máxima seria de
1:500.000.
Num papel de 50cmX50cm teríamos uma escala de 1:3.000.000, para cobri a área do
Arthur
No Caso do problema do Arthur, que temos:
long
lat
minima
-38.13
-9.24
máxima
-52.14
-22.06
amplitude
14.01
12.82
multiplicando por 111Km cada grau, teríamos um quadrado de 1500kmx1500km,
aproximadamente
usando lat_0=-15 para as coordenadas
( long -52.14) e (lat -15.54 )
a latitude -15.54 maior que -15
resultou em lambert
(841253, -75134.7) 1575.2
seria isso mesmo?
ou colocaria lat_0=-25, já que a maior latitude do exemplo do Arthur é -22.06?
Hélio
Em 27 de abril de 2013 22:42, Jônatan [via R-br] <[hidden email]> escreveu:
Quando trabalho com dados na escala continental eu uso a projeção
Azimutal equivalente de Lambert
com meridiano centrado na região de interesse, por exemplo para América
do Sul em -60. Acho que
pode ser uma alternativa.
require(sp)
require(rgdal)
p <- read.table("C:/Users/Jonatan/Downloads/Dados.txt", head =
T,blank.lines.skip=T)
pll <- SpatialPointsDataFrame(p[,1:2], data=p[3], proj4string=CRS("+proj=longlat
+datum=WGS84"))
## projecao azimutal de área equivalente de Lambert (meridiano central em -60
para America do Sul)
projlaea <- "+proj=laea +lat_0=-15 +lon_0=-60 +x_0=0 +y_0=0 +a=6370997 +b=6370997
+units=m +no_defs"
## dados reprojetados em laea, unidade em m
plaea <- spTransform(x=pll, CRSobj=CRS(projlaea))
## fazer Krigagem
...
## reprojetar objeto com resultado para latlon se desejavel (usando ultimo
comando)
resll <- spTransform(x=resultadoKrigagem, CRSobj=CRS("+proj=longlat
+datum=WGS84"))
2013/4/26 Arthur Telles Calegario <[hidden email]>
A maior dúvida é o que o hélio comentou.tenho as coordenadas em três
fusos distintos.
Assim não dá pra projeta-los juntos. entende?
_______________________________________________________________________________________________________________________________________
From: [hidden email]
To: [hidden email]
Subject: RE: [R-br] Krigagem
Date: Fri, 26 Apr 2013 14:52:22 -0300
Desconsidere o último email.Esses dados fora retirados das normais
climatológicas do inmet.
Perguntei a ele o sistema de coordenadas. Me disseram que é WGS84. É isso que
perguntou Elias?
_______________________________________________________________________________________________________________________________________
From: [hidden email]
To: [hidden email]
Subject: RE: [R-br] Krigagem
Date: Fri, 26 Apr 2013 14:47:36 -0300
Entendi.Sei uma maneira de fazer essa transformação no GIS. Mas usando somente
um fuso. Ai que esta
o problema.
Agradeço sua ajuda,
Creio que vou ter que procurar um professor de agrimensura. hehe
Arthur
_______________________________________________________________________________________________________________________________________
Date: Fri, 26 Apr 2013 13:09:46 -0300
From: [hidden email]
To: [hidden email]
Subject: Re: [R-br] Krigagem
ArthurPelo que vi seus dados estão em tres fusos, 22, 23 e 24.
Para transformar em UTM no R eu não sei de forma automática todos os dados de
uma vez só...
tirei da lista da msg de Elias T. Krainski
Coord <- SpatialPoints(cbind(dados$coord_x,dados$coord_y),
proj4string=CRS("+proj=longlat
+datum=WGS84"))
utm <- spTransform(Coord, CRS("+proj=utm +zone=22S"))
e fiz o seguinte teste
d=cbind(-17,-50) # estão na zona 22S
d
coord<- SpatialPoints((d), proj4string=CRS("+proj=longlat +datum=WGS84"))
coord
utm <- spTransform(coord, CRS("+proj=utm +zone=22S"))
utm
Coord <- SpatialPoints((d), proj4string=CRS("+proj=longlat +zone=22S
+datum=WGS84"))
Coord
utm1 <- spTransform(Coord, CRS("+proj=utm +zone=22S"))
utm1
a resposta seria W606447.747 e N8120173.351
mas em nenhum dos testes deu certo
Hélio
Em 26 de abril de 2013 10:14, Arthur T. C. [via R-br] <[hidden email]> escreveu:
Verdade. Rs. Tenho os pontos em coordenadas geográficas. em um datum
qualquer. Por
exemplo SAD 69.
Necessito projetar pra saber as distâncias linear entre os pontos
(superficial). E não cortar
a esfera terrestre.
É isso certo? Os dados estão em anexo.
Agradeço Hélio,
Arthur
_______________________________________________________________________________________________________________________________________
Date: Fri, 26 Apr 2013 10:00:37 -0300
From: [hidden email]
To: [hidden email]
Subject: Re: [R-br] Krigagem
Arthur. Parece que está fazendo uma confusão entre fuso e datum.
anexe seus dados pra que eu possa dar uma olhada, com o CRM que tem
Hélio
Em 26 de abril de 2013 09:24, Arthur T. C. [via R-br] <[hidden email]> escreveu:
Bom dia Senhores,Estou com um pequeno problemas em relação a Krigagem e
talvez os
senhores tenham uma resposta.
Na verdade o problema esta em preparar os dados. Estou trabalhando com
evapotranspiração no Rio Doce e possuo as coordenadas das estações
meteorológicas em
graus. Contudo, para a construção do semivariograma, necessito que a distância
entre as
mesmas seja linear.
Assim devo projeta-las. É um procedimento simples quando todos as estações
estão em um
mesmo 'datum'. Contudo, as estações que estão do 'datum' 23S são insuficientes
para a
construção do semivariograma. Assim necessito saber como faço para linearizar as
distâncias entre estações meteorológicas cujos limites extrapolam um 'datum'
apenas.
Entenderam? rs.
Agradeço desde já,
Arthur Calegario.
_______________________________________________
R-br mailing list
[hidden email]
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.
_______________________________________________________________________________________________________________________________________
If you reply to this email, your message will be added to the discussion below:
http://r-br.2285057.n4.nabble.com/R-br-Krigagem-tp4659130.html
To unsubscribe from R-br, click here.
NAML
--
Hélio Gallo Rocha
IFSULDEMINAS - Câmpus Muzambinho
_______________________________________________ R-br mailing list [hidden email]
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
[hidden email]
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.
[IMAGE] Dados.txt (2K) Download Attachment
_______________________________________________________________________________________________________________________________________
If you reply to this email, your message will be added to the discussion below:
http://r-br.2285057.n4.nabble.com/R-br-Krigagem-tp4659130p4659134.html
To unsubscribe from R-br, click here.
NAML
--
Hélio Gallo Rocha
IFSULDEMINAS - Câmpus Muzambinho
_______________________________________________ R-br mailing list [hidden email]
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
[hidden email]
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.
--
###############################################################
## Jônatan Dupont Tatsch
## *Pesquisador Associado Recém Doutor
## do Programa de Pós Graduação em Meteorologia
## Universidade Federal de Santa Maria (UFSM)
## Santa Maria/RS - Brasil
## Telefone: (55) 3301 2084
## www.ufsm.br/meteorologia/pos## * Pesquisador colaborador do:
## Laboratório de Clima e Biosfera LCB/IAG/USP
## Núcleo de Apoio a Pesquisas em Mudanças Climáticas (INCLINE)
## www.incline.iag.usp.br
###############################################################
_______________________________________________
R-br mailing list
[hidden email]
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.
_______________________________________________________________________________________________________________________________________
If you reply to this email, your message will be added to the discussion below:
http://r-br.2285057.n4.nabble.com/R-br-Krigagem-tp4659130p4659157.html
To unsubscribe from R-br, click here.
NAML
--
Hélio Gallo Rocha
IFSULDEMINAS - Câmpus Muzambinho
_______________________________________________
R-br mailing list
[hidden email]
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.
--
###############################################################
## Jônatan Dupont Tatsch
## *Pesquisador Associado Recém Doutor
## do Programa de Pós Graduação em Meteorologia
## Universidade Federal de Santa Maria (UFSM)
## Santa Maria/RS - Brasil
## Telefone: (55) 3301 2084
## www.ufsm.br/meteorologia/pos## * Pesquisador colaborador do:
## Laboratório de Clima e Biosfera LCB/IAG/USP
## Núcleo de Apoio a Pesquisas em Mudanças Climáticas (INCLINE)
## www.incline.iag.usp.br
###############################################################
_______________________________________________
R-br mailing list
[hidden email]
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.
_______________________________________________________________________________________________________________________________________
If you reply to this email, your message will be added to the discussion below:
http://r-br.2285057.n4.nabble.com/R-br-Krigagem-tp4659130p4659167.html
To unsubscribe from R-br, click here.
NAML
--
Hélio Gallo Rocha
IFSULDEMINAS - Câmpus Muzambinho
_______________________________________________
R-br mailing list
[hidden email]
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.
--
###############################################################
## Jônatan Dupont Tatsch
## *Pesquisador Associado Recém Doutor
## do Programa de Pós Graduação em Meteorologia
## Universidade Federal de Santa Maria (UFSM)
## Santa Maria/RS - Brasil
## Telefone: (55) 3301 2084
## www.ufsm.br/meteorologia/pos## * Pesquisador colaborador do:
## Laboratório de Clima e Biosfera LCB/IAG/USP
## Núcleo de Apoio a Pesquisas em Mudanças Climáticas (INCLINE)
## www.incline.iag.usp.br
###############################################################
_______________________________________________
R-br mailing list
[hidden email]
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.
_______________________________________________________________________________________________________________________________________
If you reply to this email, your message will be added to the discussion below:
http://r-br.2285057.n4.nabble.com/R-br-Krigagem-tp4659130p4659201.html
To unsubscribe from R-br, click here.
NAML
--
Hélio Gallo Rocha
IFSULDEMINAS - Câmpus Muzambinho
_______________________________________________ R-br mailing list
[email protected]
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.
