2016-11-06 47 views
0

我們有一個數據集,其中包含緯度和經度座標以及屬性信息,每個數據集都位於其自己的單獨列中,並以數字形式存儲。這些座標根據地理座標系WGS 1984進行了地理編碼。我們知道,我們的數據中存在顯着的空間自相關,我們希望使用「sp」包在氣泡圖中對其進行可視化。我們在線模擬我們的示例,例如:https://beckmw.wordpress.com/2013/01/07/breaking-the-rules-with-spatial-correlation/。然而,當我們嘗試使用座標中的 「SP」 命令,我們不斷收到錯誤消息:用於檢測R中空間自相關的座標投影的問題

代碼示例: COORDS < - data.frame(鉛$ X,鉛$ Y)座標 (COORDS) < - C( 「lead6.X」, 「lead6.Y」)

錯誤如果(NCHAR(projargs)== 0)projargs < - as.character(NA)缺失值,其中TRUE/FALSE需要

我們無法加載我們的直接代碼,因爲它敏感並託管在虛擬環境中,無法訪問互聯網。有沒有人有想法,爲什麼這可能會發生?我們已經研究過proj4包,但不知道如何指定投影系統(或者甚至是我們正在獲取的錯誤?)。如果有人知道R中的任何其他軟件包或可視化空間自相關的方法,那麼也會非常感激。

+0

嘗試使用不敏感的數據重新創建錯誤消息,並在可重現的示例中提供這些錯誤消息。爲什麼不直接使用'coords = SpatialPoints(data.frame(lead $ X,lead $ Y),CRS(「+ init = epsg:4326」))'' –

回答

0

您的代碼有點「奇怪」:似乎您正在嘗試構建僅包含座標的數據集。 AFAIU,你可能需要的東西在這條線:

data <- data.frame(lead$X, lead$Y, lead$Z) 

,對應於鉛$ Z上通用的「變量」你要檢查,然後

coordinates(data) <- c('X','Y')` 
proj4string(data) <- "+init=epsg:4326" 

,這應該給你一個適當的「 SpatialPointsDataframe「(包含經緯度WGS84地理座標的空間點數據框)(第一行也可以丟棄,並且將空間點數據框的數據中的所有變量保留)。

HTH