我試圖獲得紐約地區經緯度座標的郵政編碼如何批量反轉R中的地理編碼?
我試圖使用谷歌的反向地理編碼器API,但每天限制在2500次點擊,因此可以批量處理我的數據幀。
接下來,我嘗試使用帶有數據集郵編的庫(zipcode),但無法將緯度經度與火車數據集的座標進行匹配,因爲所有緯度經度座標都不在數據集中。
此外,雖然使用KNN來預測數據集的郵政編碼,但無法獲得正確的結果。
zipcode_latlon = zipcode[zipcode$state=="NY",c(1,4,5)]
train_latlon = train_data[,c("latitude","longitude")]
zip1 = rep(10007, nrow(train_latlon))
zip1 = as.character(zip1)
train_latlon = cbind(zip1, train_latlon)
colnames(train_latlon) = c("zip","latitude","longitude")
knn_fit = knn(zipcode_latlon, train_latlon,zipcode_latlon$zip, k=1)
需要知道如何我可以從Lat Long網批量獲得郵編,任何方法都好於R.
是否有效? –
謝謝!是的,它爲我工作。 但是很難理解Shape文件格式。 無法理解spTransform()中的屬性 –
原始shapefile位於座標參考系NAD83中。 spTransform只是給'zips'一個不同的參考系統 - 在這個例子中是WGS84。點和多邊形具有相同的CRS是很重要的。也可以按照導入後的方式保留'zips',然後用'+ proj = longlat + datum = NAD83 + no_defs + ellps = GRS80 + towgs84 = 0,0,0給'spdf' NAD83 CRS '。無論如何,沒有參考系統,空間數據就毫無意義。 –