4

私はアメリカの地図と、プロットしたい緯度経度のリストを持っています。これが機能したら、「世界」マップにも切り替えたいと思います。マップは生成されますが、マップ上にポイントが表示されません。TSV ファイルの最初の行には、次のヘッダーが含まれています。

LONG{tab}LAT
R は「トラフィック」テーブルを読み取っているように見えます。私は何を間違っていますか？

ライブラリ(「マップ」)

トラフィック = read.table("C:/temp/traffic_10.40.tsv", header=T, sep="\t")
png(ファイル名="C:/temp/usa.png"、幅=850、高さ=600、bg="白")
map('state', plot = TRUE, fill = FALSE, col = palette())
title("行き先")
ポイント(x=traffic$LONG,y=traffic$LAT,col='red',cex=0.75)
dev.off()

編集

> dput(トラフィック)
構造体 (リスト (LONG = c(47.6218、32.7942、34.1121、40.0068、47.6218、
33.9553、33.7629、40.0068、39.05、38.1075、33.7629、32.769、37.3857、
29.4576、34.1674、38.8147、32.7942、31.1429、40.3254、30.3059、
38.2248、47.6218、33.9553、38.1075、27.1943、29.4576、30.5175、
38.5668、42.6657、40.2982、32.7539、40.6698、47.6742、32.7942、
47.6218、35.8219)、LAT = c(-122.35、-96.7653、-118.411、-75.1347、
-122.35、-83.3937、-84.4226、-75.1347、-77.4833、-122.264、-84.4226、
-96.5998、-122.026、-98.5054、-84.8014、-77.0647、-96.7653、-81.471、
-78.9195、-97.7505、-85.7412、-122.35、-83.3937、-122.264、-80.2432、
-98.5054、-97.6721、-121.467、-73.799、-111.698、-97.3363、-73.9438、
-122.115, -96.7653, -122.35, -78.6588)), .Names = c("LONG", "LAT"
)、クラス = "data.frame"、row.names = c(NA, -36L))

また、私はRの初心者であり、Googleでこれを見つけようとしましたが、正確に何を検索すればよいかわからないため、成功は限られています。

4

3 に答える 3

8

問題は、後のコードではなく、データ セットにあります。

最初の点には 1 つの座標47.6218ともう 1 つの座標があり-122.35ます。 緯度は範囲外にすることはできない[-90,90]ため、経度は、データ セットの反対で ある-122.35および latitudeである必要があります。これは、シアトル スペース ニードルの少し北にあります。47.6218

x(水平) は伝統的に経度または東経です。y(垂直) は伝統的に緯度または北座標です

于 2012-04-05T21:14:19.463 に答える
1

修理済み。$LAT と $LONG を交換しました。今では完全に機能します。

ポイント(x=traffic$LAT,y=traffic$LONG,col='red',cex=0.75)
于 2012-04-05T17:15:23.013 に答える
1

1 つのオプションは、これにgoogleVisを使用することです。ポイントを Google マップ スタイルでプロットできます。

于 2012-04-05T16:11:44.933 に答える