วิธีแปลง lat และ longs ของฉันเป็น UTM ภายใน dataframe ของฉัน

ตอนนี้ฉันติดขัดนิดหน่อย ฉันสามารถสร้างกรอบข้อมูลจุดเชิงพื้นที่ได้ และจากนี้ ฉันได้สร้างออบเจ็กต์ของคลาส 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

-> ฉันสามารถสร้างคะแนน Spatial ของคลาสอย่างเป็นทางการซึ่งกำหนดจาก lats และ longs ของฉันในรูปแบบ UTM ได้ แต่นั่นเป็นวัตถุที่แยกจากกันในตอนนี้:

โครงสร้าง:

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"

โดยพื้นฐานแล้วฉันต้องการให้ค่าของฉันสำหรับ x และ y ใน ltraj อยู่ใน UTM ฉันคิดว่าฉันอาจจะต้องทำสิ่งนี้ก่อนที่จะสร้างวัตถุ ltraj หรือแม้แต่ dataframe จุดเชิงพื้นที่

ฉันสงสัยว่ามีใครมีคำแนะนำในการเปลี่ยน lats และ longs ของฉันเป็น UTM ภายใน dataframe ดั้งเดิมหรือไม่ หรือเพื่อทำให้เป็น dataframe เชิงพื้นที่แล้วผูกมันอีกครั้งกับ dataframe ที่มี "หัวเรื่อง" และ "วันที่"?

ขอแสดงความนับถือ,

แซม ริคเคน


person SamR    schedule 05.12.2016    source แหล่งที่มา
comment
แพ็คเกจ rgdal สามารถทำได้ ดูโพสต์ก่อนหน้า   -  person G5W    schedule 05.12.2016
comment
เหตุใดจึงมีป้ายกำกับว่า python   -  person DYZ    schedule 05.12.2016
comment
ดู ?sp::spTransform   -  person johannes    schedule 12.12.2016