ตอนนี้ฉันติดขัดนิดหน่อย ฉันสามารถสร้างกรอบข้อมูลจุดเชิงพื้นที่ได้ และจากนี้ ฉันได้สร้างออบเจ็กต์ของคลาส 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 ที่มี "หัวเรื่อง" และ "วันที่"?
ขอแสดงความนับถือ,
แซม ริคเคน
python
- person DYZ   schedule 05.12.2016?sp::spTransform
- person johannes   schedule 12.12.2016