Muito obrigada pessoal, deu certo..<br><br><div class="gmail_quote">Em 31 de outubro de 2012 09:19, Eder Comunello <span dir="ltr"><<a href="mailto:ecomunel@gmail.com" target="_blank">ecomunel@gmail.com</a>></span> escreveu:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div>Natalia e colegas, bom dia!</div><div><br></div><div>Ainda não estou muito familiarizado com o R, mas acho que este exemplo que rodei irá ajudar...</div>
<div><br></div><div>Não utilizei o project(), mas sim o spTransform()...</div>
<div><br></div><div>#############################################################</div><div><br></div><div>#install.packages("rgdal", dep=T)</div><div>#install.packages("sp", dep=T)</div><div>require(sp); require(rgdal)</div>
<div><br></div><div>### long</div><div>coord.x = c(-54.607845, -54.599831, -54.430241, -54.258459, -54.344408, </div><div>-54.333258, -53.607359, -53.589245, -53.546762, -53.105123, </div><div>-52.972362, -52.107579, -52.089572, -51.940453, -51.718151, </div>
<div>-51.263728, -50.741531, -50.734130, -50.667779, -50.662230, </div><div>-49.986102, -49.970018, -49.901474, -49.911543, -49.678260, </div><div>-49.654858, -49.610676, -49.599658, -49.511570, -49.200709, </div><div>-49.222508, -49.317621, -49.311839, -49.275882, -48.581791, </div>
<div>-48.578971, -48.582604, -48.568999, -48.528527, -48.411012, </div><div>-48.251622, -48.027824, -48.029886, -48.023501, -48.590796, </div><div>-48.594522, -48.643601, -49.175230, -49.555033, -49.567207, </div><div>-49.689784, -49.941462, -49.942167, -50.571750, -50.591636, </div>
<div>-50.553035, -50.718453, -51.208345, -51.216170, -51.292329,</div><div>-51.239319, -51.258597, -51.283372, -51.411480, -53.551638, </div><div>-53.552753, -53.643791, -53.835953, -53.842011, -53.858677, </div><div>-53.955544, -54.099576, -54.108952, -54.429627, -54.438295, </div>
<div>-54.593660, -54.593683, -54.619090, -54.617907, -54.607845)</div><div><br></div><div>### lat</div><div>coord.y = c(-25.446263, -25.434436, -25.159459, -24.357597, -24.140630, </div><div>-24.116248, -22.951327, -22.936230, -22.890090, -22.621923, </div>
<div>-22.570479, -22.516685, -22.521237, -22.567335, -22.669007, </div><div>-22.669047, -22.961387, -22.962502, -22.897375, -22.895686, </div><div>-22.897501, -22.910601, -22.999147, -23.051110, -23.165054, </div><div>-23.210043, -23.851684, -23.857080, -23.925858, -24.344464, </div>
<div>-24.362727, -24.566630, -24.663615, -24.693608, -24.671609, </div><div>-24.685526, -25.051554, -25.051240, -25.100784, -24.980119, </div><div>-24.978739, -25.221284, -25.226189, -25.230964, -25.976982, </div><div>-25.977087, -25.955943, -26.001844, -26.237326, -26.233284, </div>
<div>-26.188456, -26.008879, -26.008876, -26.003790, -26.022134, </div><div>-26.055436, -26.244779, -26.297107, -26.313987, -26.404261, </div><div>-26.628383, -26.625825, -26.653990, -26.717283, -26.292744, </div><div>-26.292051, -26.253072, -25.971930, -25.687925, -25.680270, </div>
<div>-25.646984, -25.495677, -25.495742, -25.695430, -25.691293, </div><div>-25.592352, -25.544262, -25.457555, -25.452072, -25.446263)</div><div><br></div><div>pr.df <- data.frame(coord.x,coord.y); head(pr.df)</div><div>
<br></div><div>### criar um objeto "SpatialPoints"</div><div>pr.spgeo <- pr.df; head (pr.spgeo)</div><div>coordinates(pr.spgeo) <- c('coord.x','coord.y')</div><div><br></div><div>### definir projeção e datum de entrada</div>
<div>proj4string(pr.spgeo) <- CRS('+proj=longlat +datum=WGS84')</div><div>class(pr.spgeo)</div><div>head(pr.spgeo)</div><div>plot(pr.spgeo)</div><div><br></div><div>### lista parâmetros disponíveis</div><div>projInfo("proj")</div>
<div>projInfo("datum") </div><div>projInfo("ellps") </div><div><br></div><div>### transforma para utm</div><div>### atenção para o uso do parâmetro '+south'</div><div>### ter claro que o Paraná utiliza mais de uma Zona UTM e maior parte do estado em 22S</div>
<div>pr.sputm <- spTransform(pr.spgeo, CRS('+proj=utm +zone=22 +south')); head(pr.sputm) </div><div>plot(pr.sputm)</div><div><br></div><div class="gmail_extra">#############################################################</div>
<div class="gmail_extra"><br></div><div class="gmail_extra"><div>Esse exemplo assumiu o WGS84 como datum padrão de entrada e saída. Para datums (data) diferentes tem que dar uma estudada mais a fundo, porque tem poucas opções pré-configuradas no pacote que usei (como pode ser visto usando <projInfo("datum")>).</div>
<div><br></div>-- <br><font face="arial, helvetica, sans-serif"><br>================================================<br>Éder Comunello<br><br>Ph.D. Student in Agricultural Systems Engineering (USP/ESALQ)<br>Piracicaba, SP, Brazil [22 42.7'S, 47 37.8'W]<br>
<br>Researcher at Embrapa Western Region Agriculture</font><div><font face="arial, helvetica, sans-serif">Dourados, MS, Brazil [</font>22 16.5'S, 54 49.0'W<span style="font-family:arial,helvetica,sans-serif">]</span></div>
<div><font face="arial, helvetica, sans-serif">================================================<br>UTC-03:00</font></div><div><font face="arial, helvetica, sans-serif"><br></font></div></div><div class="gmail_extra"><br>
<div class="gmail_quote">
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex">
<br>
----------------------------------------------------------------------<div><div class="h5"><br><br>
<br>
2012/10/29 Natalia Martins <<a href="https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=nsmbarreto@gmail.com" target="_blank">nsmbarreto@gmail.com</a>><br>
<br>
> Walmes, bom dia.<br>
> Estou tentando transformar minhas coordenadas geograficas (das estações<br>
> metereologicas do PR)<br>
> em utm, procurando no help do R encontrei:<br>
><br>
> data(state)<br>
> s <- project(state.center, "+proj=merc")<br>
><br>
> mas eu não estou conseguindo entender qual proj eu devo usar, no link<br>
><br>
> <a href="http://www.remotesensing.org/geotiff/proj_list/transverse_mercator.html" target="_blank">http://www.remotesensing.org/geotiff/proj_list/transverse_mercator.html</a><br>
><br>
> até encontri mas quando tento fazer aparece o seguinte erro<br>
><br>
> s <- project(state.center, "+proj=utm")<br>
> Erro em project(state.center, "+proj=utm") : elliptical usage required<br>
><br>
><br>
> se vc puder me ajudar...<br>
><br>
> Muito obrigada!!<br>
><br>
> Em 28 de outubro de 2012 11:11, Walmes Zeviani <<a href="https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=walmes@ufpr.br" target="_blank">walmes@ufpr.br</a>> escreveu:<br>
><br>
> Sim, para calculo de distâncias você tem que ter coordenadas decimais e<br>
>> não exagesimais (º ' '').<br>
>><br>
>> ==========================================================================<br>
>> Walmes Marques Zeviani<br>
>> LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W)<br>
>> Departamento de Estatística - Universidade Federal do Paraná<br>
>> fone: <a href="tel:%28%2B55%29%2041%203361%203573" value="+554133613573" target="_blank">(+55) 41 3361 3573</a><br>
>> VoIP: (3361 3600) 1053 1173<br>
>> e-mail: <a href="https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=walmes@ufpr.br" target="_blank">walmes@ufpr.br</a><br>
>> skype: walmeszeviani<br>
>> twitter: @walmeszeviani<br>
>> homepage: <a href="http://www.leg.ufpr.br/~walmes" target="_blank">http://www.leg.ufpr.br/~walmes</a><br>
>> linux user number: 531218<br>
>> ==========================================================================<br>
>><br>
><br>
><br>
><br>
> --<br>
><br>
> Natália da Silva Martins<br>
> Bacharel em Estatística - Universidade Estadual de Maringá/ UEM<br>
> Mestranda em Estatística e Experimentação Agronômica - ESALQ/ USP<br>
> Contato: <a href="tel:%2819%29%208306-4743" value="+551983064743" target="_blank">(19) 8306-4743</a><br>
><br>
><br>
> Mensagem Verificada Pelo AntiSpam UFPR<br>
><br>
><br></div></div></blockquote></div>
</div>
<br>_______________________________________________<br>
R-br mailing list<br>
<a href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a><br>
<a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br>
Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.<br></blockquote></div><br><br clear="all"><div><br></div>-- <br>
<br style="font-family:Tahoma;font-size:13px"><font face="'Times New Roman'" size="4">Natália da Silva Martins<br>Bacharel em Estatística - Universidade Estadual de Maringá/ UEM<br>Mestranda em Estatística e Experimentação Agronômica - ESALQ/ USP<br>
Contato: (19) 8306-4743</font> <br><div><br></div><br>