o rgdal possui dependencias externas que voce tem que instalar.
Vaja pa página e documentacao do pacote

se estiver usando linux (com apt-get)

apt-get --build-dep rgdal

On Mon, 6 May 2013, Arthur Telles Calegario wrote:

> 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: heliogallorocha@gmail.com
> To: r-br@listas.c3sl.ufpr.br
> 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] <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ô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 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.
> 
>