使用ggmap在地圖上繪制等高線

Question

我有洛杉磯港地區的顆粒物濃度差異（之后-之前）。 我正在嘗試使用ggmap在地圖上繪制濃度等高線，但結果看起來大不相同。 我使用的代碼如下所示（數據在代碼下方）：

碼

install.packages（ 'ggmap'）

library(ggmap)

PM = read.csv('data.csv', stringsAsFactors = FALSE)

獲取經度和緯度

geocode("Port of Los Angeles") # Not centered
geocode("Compton, CA") # Now centered

使用Compton LON和LAT

POLA = c(lon = -118.220071, lat =  33.895849)
POLA.map = get_map(location = POLA, zoom = 10, color = 'bw')
ggmap(POLA.map) + geom_point(data = PM, mapping = aes(Longitude, Latitude)) + 
        stat_density2d(data = PM, mapping = aes(x = Longitude, y = Latitude, fill=..level..), geom = "polygon", alpha = 0.3, contour = TRUE)

但是，等高線圖應具有以下模式：

https://drive.google.com/file/d/0B3XVjcsci0y3VDBTc01PYkhOckE/view?usp=sharing

ggplot(PM, aes(UTM.X, UTM.Y)) + geom_tile(aes(fill = Value), alpha = 0.8, color = "black") + 
        scale_fill_gradient(low = 'green', high = 'red')

數據：Col-1：經度，Col-2：緯度，Col-3：UTM-X，Col-4：UTM-Y，Col-5：值

UTM坐標單位：米，UTM區域= 11 N，基准= WGS84。 此處提供數據： https ： //drive.google.com/file/d/0B3XVjcsci0y3LUpudko1S2c1cnc/view？usp = sharing

Answer 1

stat_density2d用於繪制密度圖，例如，暗點多的地方，而亮點少的地方。 您有一個帶有Value屬性的常規網格，而不是密度圖。

因此，您應該使用geom_tile來獲取常規的網格圖。 但是您的經緯度坐標不會形成軸對齊的網格。 嘗試這個：

ggplot(data = PM, mapping = aes(x = Longitude, y = Latitude, fill=Value)) + geom_tile()

而您得到一個空白圖，請嘗試以下操作：

ggplot(data = PM, mapping = aes(x = UTM.X, y = UTM.Y, fill=Value)) + geom_tile()

然后你得到了情節。 當然，它與ggmap背景不在同一坐標系中。

你也許可以使用基礎R的contourLines函數來獲取的UTM坐標輪廓線的坐標，使SpatialLinesDataFrame ，然后轉換為緯度長，並添加到ggmap 。

獲得看起來像網格圖的另一種可能性是使用具有正方形的點作為形狀。

ggmap(POLA.map)  + geom_point(data = PM, mapping = aes(Longitude, Latitude, colour=Value), size=4, alpha=0.5, shape=15) + scale_colour_gradient(low = 'green', high = 'red')

在某些偽像中，網格單元重疊，看起來有點像單元輪廓，並且圖例顯示時沒有不透明，因此看起來比單元更飽和。 您必須正確設置size參數，因為它取決於圖形設備的大小。

如果所有步驟均失敗，則將數據轉換為柵格數據包柵格對象，將其另存為GeoTIFF，然后將其加載到QGIS中，后者可以將UTM網格實時投影到Lat-long上。

QGIS還具有一些不錯的混合模式，因此您可以輕松地做到這一點：

注意這不是透明性，這是乘法混合。 透明度會導致深色被淡出，而乘法混合則讓黑色透徹，因此標簽和底圖細節仍然可見。

另外，請注意柵格是如何未與軸對齊的（在底部尤為明顯）。