如何在我的数据帧中将我的lat和long转换为UTM

时间:2016-12-05 04:14:50

标签: r geospatial adehabitathr

我现在有点陷入困境,我已经能够创建一个空间点数据框,并且我已经创建了Ltraj类的对象。我需要做进一步的分析。但我的x和y坐标不是UTM,如果我在线下进行进一步的分析,可能会出现问题。

格式:

      x        y                date        dx         dy         dist   dt             
1 -32.09245 116.0426 2015-08-07 00:22:00 -2.19e-05  0.0000194 2.925696e-05 1800 ...

结构:

List of 1
 $ :'data.frame':   109 obs. of  10 variables:
  ..$ x        : num [1:109] -32.1 -32.1 -32.1 -32.1 -32.1 ...
  ..$ y        : num [1:109] 116 116 116 116 116 ...
  ..$ date     : POSIXct[1:109], format: "2015-08-07 00:22:00" "2015-08-07   00:52:00" "2015-08-07 01:22:00" "2015-08-07 01:52:00" ...
  ..$ dx       : num [1:109] -2.19e-05 -5.73e-05 -5.15e-05 4.52e-05 -4.96e- 05 ...
  ..$ dy       : num [1:109] 1.94e-05 -3.21e-04 -2.61e-05 2.75e-04 -1.06e-04 ...
  ..$ dist     : num [1:109] 2.93e-05 3.26e-04 5.77e-05 2.79e-04 1.17e-04 ...
  ..$ dt       : num [1:109] 1800 1800 1800 3840 1800 3600 1740 1920 4680 900 ...
  ..$ R2n      : num [1:109] 0.00 8.56e-10 9.71e-08 1.24e-07 1.00e-08 ...
  ..$ abs.angle: num [1:109] 2.42 -1.75 -2.67 1.41 -2.01 ...
  ..$ rel.angle: num [1:109] NA 2.119 -0.925 -2.203 2.865 ...
  ..- attr(*, "id")= chr "2172"
  ..- attr(*, "burst")= chr "2172"
  ..- attr(*, "infolocs")='data.frame': 109 obs. of  1 variable:
  .. ..$ pkey: Factor w/ 109 levels "2172.2015-08-07 00:22:00",..: 1 2 3 4 5 6 7 8 9 10 ...
 - attr(*, "class")= chr [1:2] "ltraj" "list"
 - attr(*, "typeII")= logi TRUE
 - attr(*, "regular")= logi FALSE
 - attr(*, "proj4string")=Formal class 'CRS' [package "sp"] with 1 slot
  .. ..@ projargs: chr NA

- >我能够以UTM格式从我的lats和long创建一个正式的类空间点集,但现在这是一个单独的对象:

结构:

Formal class 'SpatialPoints' [package "sp"] with 3 slots
  ..@ coords     : num [1:109, 1:2] 409662 409664 409634 409631 409657 ...
  .. ..- attr(*, "dimnames")=List of 2
  .. .. ..$ : NULL
  .. .. ..$ : chr [1:2] "coords.x1" "coords.x2"
  ..@ bbox       : num [1:2, 1:2] 406647 13536726 415659 13551107
  .. ..- attr(*, "dimnames")=List of 2
  .. .. ..$ : chr [1:2] "coords.x1" "coords.x2"
  .. .. ..$ : chr [1:2] "min" "max"
  ..@ proj4string:Formal class 'CRS' [package "sp"] with 1 slot
  .. .. ..@ projargs: chr "+init=epsg:32750 +proj=utm +zone=50 +south    +datum=WGS84 +units=m +no_defs +ellps=WGS84 +towgs84=0,0,0"

我基本上希望ltraj中x和y的值都是UTM, 我想在创建对象ltraj甚至空间点数据帧之前我可能不得不这样做。

我想知道是否有人建议在原始数据框架内改变我的lats和longs到UTM;或者将其转换为空间数据框,然后再次将它们与包含“主题”和“日期”的数据框绑定?

亲切的问候,

Sam Rycken

0 个答案:

没有答案