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] <ml-node+s2285057n4659167h45@n4.nabble.com> 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

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