5

9,000 以上の UTM 座標の .csv ファイルを 10 進度に変換したいのですが、少し問題があります。ここや他の場所に投稿されたいくつかの投稿を検索しましたが、UTM のセットを使用可能で正確な緯度/経度に変換するソリューションが見つからないようです。

基本的に 2 つの質問があります。1) 私のコードに問題がある人はいますか。2) UTM を緯度/経度に変換するのをやめて、Rgooglemaps パッケージで UTM を使用することに慣れている人はいますか?

私のコードとデータの例を次に示します。

データ：

>head(utm)
-Northing Easting
1  4236576  615805
2  4236576  615805
3  4236576  615805
4  4236576  615805
5  4236576  615805
6  4236576  615805

これまでのコード:

utm <- read.csv(file="utm.csv", header=TRUE, sep=",")
library(rgdal)
utm <- utm[complete.cases(utm),]
utm1 <- data.frame(x=utm$Northing,y=utm$Easting) 
coordinates(utm1) <- ~x+y 
class(utm1)
proj4string(utm1) <- CRS("+proj=utm +zone=10 +datum=WGS84 +units=m +ellps=WGS84") 
utm2 <- spTransform(utm1,CRS("+proj=longlat +datum=WGS84"))

結果

> head(utm2)
SpatialPoints:
             x        y
[1,] -91.08516 4.727323
[2,] -91.08516 4.727323
[3,] -91.08516 4.727323
[4,] -91.08516 4.727323
[5,] -91.08516 4.727323
[6,] -91.08516 4.727323
Coordinate Reference System (CRS) arguments: +proj=longlat +datum=WGS84 +ellps=WGS84
+towgs84=0,0,0

したがって、出力は得られますが、適切な出力が得られません。私がここに欠けているものはありますか？また、ヒートマップとカーネル密度プロットを作成するために「Rgooglemaps」パッケージを使用することを計画していました。

4

2 に答える 2

1

X と Y を交換するだけでよいと思います。つまり、北を Y、東を X とします。元のコードは、この変更でうまくいきました。

于 2020-05-12T12:13:26.993 に答える