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] <ml-node+s2285057n4659201h99@n4.nabble.com> 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ônatan
Fiz 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

## 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?


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?


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

Arthur
Pelo 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.

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