我試圖在 R 中重新投影的 .csv 檔案中有 POINT_X、POINT_Y 和 z 點列。最初，它被投影在EPSG:26777參考系統中，但我試圖將其投影回WGS 84. 但是，當我將 .csv 檔案匯入 R 并嘗試轉換時，它不會更改“幾何”列中的點值。如何將它們更改為緯度和經度坐標？
這是在原始投影中匯出的資料：
data <- data.frame(point_x = c(2387371.258,
2387411.214,
2387451.258,
2387529.427,
2387411.258,
2387191.258,
2387491.258,
2387331.258,
2387116.239,
2387291.258,
2387371.258,
2387431.258,
2387331.258,
2387271.258,
2387511.258,
2387331.258,
2387211.258,
2387351.258,
2387112.752,
2387311.22,
2387411.258,
2387191.258,
2387271.258,
2387491.258,
2387131.258,
2387231.258,
2387291.258,
2387371.052,
2387491.258,
2387171.258,
2387131.258,
2387371.258,
2387113.437,
2387211.258,
2387531.258,
2387331.258,
2387371.258,
2387311.258,
2387291.258,
2387511.258,
2387211.258,
2387491.258,
2387391.258,
2387530.847,
2387542.074,
2387491.258,
2387511.258,
2387271.258,
2387161.36,
2387112.772,
2387451.258,
2387431.258,
2387511.258,
2387127.551,
2387531.258,
2387251.258,
2387151.258,
2387231.258,
2387231.258,
2387191.258,
2387511.258,
2387171.258,
2387231.258,
2387431.258,
2387531.258,
2387231.258,
2387211.258,
2387112.281,
2387191.258,
2387391.258,
2387182.737,
2387351.258,
2387271.258,
2387371.258,
2387451.197,
2387171.258,
2387451.258,
2387431.258,
2387114.101,
2387491.219,
2387151.258,
2387451.258,
2387171.258,
2387528.481,
2387251.258,
2387251.258,
2387451.258,
2387112.193,
2387391.258,
2387151.258,
2387211.258,
2387171.258,
2387351.258,
2387131.258,
2387511.258,
2387411.258,
2387351.258,
2387411.258,
2387231.258,
2387211.224,
2387311.258,
2387191.258,
2387351.258,
2387351.258,
2387351.258,
2387171.258,
2387491.258,
2387371.258,
2387171.225,
2387251.258,
2387118.495,
2387231.258,
2387471.258,
2387171.258,
2387131.258,
2387191.258,
2387371.258,
2387151.258,
2387131.258,
2387411.258,
2387112.44,
2387191.258,
2387524.317,
2387391.258,
2387371.258,
2387251.258,
2387451.258,
2387371.258,
2387191.258,
2387231.258,
2387211.258,
2387331.258,
2387391.258,
2387311.258,
2387191.258,
2387431.258,
2387271.258,
2387271.258,
2387323.356,
2387231.258,
2387351.258,
2387191.258,
2387211.258,
2387231.258,
2387431.258,
2387191.258,
2387351.258,
2387491.258,
2387431.258),
point_y = c(289191.3991,
289054.1412,
289391.3991,
289351.3027,
289091.399,
289071.399,
289251.399,
289291.3993,
289271.587,
289311.3992,
289311.3992,
289151.3992,
289151.3992,
289131.3989,
289411.399,
289391.3991,
289071.399,
289331.3992,
289611.4751,
289053.1858,
289211.3991,
289511.3991,
289511.3991,
289171.3992,
289131.3989,
289371.3991,
289131.3989,
289529.3496,
289451.3993,
289271.399,
289211.3991,
289371.3991,
289391.4699,
289231.399,
289251.399,
289411.399,
289331.3992,
289411.399,
289531.3991,
289371.3991,
289331.3992,
289431.399,
289251.399,
289291.3169,
289055.4637,
289131.3989,
289071.399,
289231.399,
289601.478,
289431.4643,
289291.3993,
289271.399,
289131.3989,
289626.1925,
289131.3989,
289291.3993,
289231.399,
289251.399,
289351.3991,
289351.3991,
289151.3992,
289431.399,
289111.399,
289091.399,
289111.399,
289491.3992,
289311.3992,
289571.3209,
289191.3991,
289071.399,
289582.6287,
289411.399,
289331.3992,
289471.3992,
289054.5359,
289191.3991,
289151.3992,
289231.399,
289351.4766,
289511.3581,
289071.399,
289171.3992,
289331.3992,
289391.2898,
289131.3989,
289411.399,
289311.3992,
289628.765,
289471.3992,
289411.399,
289431.399,
289251.399,
289471.3992,
289471.3992,
289451.3993,
289391.3991,
289151.3992,
289111.399,
289271.399,
289052.2304,
289111.399,
289231.399,
289311.3992,
289191.3991,
289371.3991,
289571.399,
289211.3991,
289511.3991,
289051.8482,
289171.3992,
289191.7413,
289451.3993,
289251.399,
289351.3991,
289231.399,
289331.3992,
289071.399,
289511.3991,
289591.399,
289451.3993,
289451.4619,
289491.3992,
289487.6817,
289091.399,
289111.399,
289431.399,
289231.399,
289211.3991,
289091.399,
289291.3993,
289171.3992,
289371.3991,
289411.399,
289391.3991,
289211.3991,
289491.3992,
289251.399,
289111.399,
289541.5201,
289311.3992,
289211.3991,
289291.3993,
289091.399,
289231.399,
289431.399,
289111.399,
289251.399,
289491.3992,
289071.399),
point_z = c(0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0))
這是我正在嘗試做的事情：
#Open packages
library(pacman)
p_load(raster,
sf,
dplyr,
ggplot2,
scales,
magrittr,
gstat,
gridExtra,
raster,
sp,
automap,
mapview,
leaflet,
rgdal)
#Read in data, provided above
data <- read.csv('data.csv')
#Change to a spatial object and reproject to WGS84
data <- st_as_sf(data, coords = c('point_x', 'point_y'), crs = crs(' proj=longlat ellps=WGS84 datum=WGS84 no_defs'))
該crs(data)命令的結果是：
Coordinate Reference System:
Deprecated Proj.4 representation: proj=longlat datum=WGS84 no_defs
WKT2 2019 representation:
GEOGCRS["unknown",
DATUM["World Geodetic System 1984",
ELLIPSOID["WGS 84",6378137,298.257223563,
LENGTHUNIT["metre",1]],
ID["EPSG",6326]],
PRIMEM["Greenwich",0,
ANGLEUNIT["degree",0.0174532925199433],
ID["EPSG",8901]],
CS[ellipsoidal,2],
AXIS["longitude",east,
ORDER[1],
ANGLEUNIT["degree",0.0174532925199433,
ID["EPSG",9122]]],
AXIS["latitude",north,
ORDER[2],
ANGLEUNIT["degree",0.0174532925199433,
ID["EPSG",9122]]]]
上面的結果表明它已經被成功地重新投影。但是，運行glimpse(data)命令顯示坐標還在原來的投影系統中：
Rows: 300
Columns: 2
$ point_z <dbl> 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0~
$ geometry <POINT [°]> POINT (2387371 289191.4), POINT (2387411 289054.1), ~
如何將這些點更改為緯度和經度點？
uj5u.com熱心網友回復：
要將資料轉換為緯度/經度，首先需要使用原始 crs 值加載它。
data <- st_as_sf(data, coords = c('point_x', 'point_y'),
crs = 26777)
您將獲得sf原始資料的物件，具有原始坐標，但現在st_transform將知道如何將其轉換為其他坐標系：
data_latlong <- st_transform(data, crs = 4326)
應該給你緯度/經度值的坐標。
轉載請註明出處，本文鏈接：https://www.uj5u.com/qukuanlian/449409.html