结合coord_proj和geom_raster

Question

我正在寻找一种在小规模地图上绘制栅格数据（使用ggplot和geom_raster）的方法。我想使用ggalt和coord_proj来'放大'shapefile的特定区域，但我遇到了错误geom_raster only works with Cartesian coordinates

ggplot() +
    geom_polygon(data = land_df, aes(long, lat, group = group), fill = 'grey25')+
    geom_raster(data = df, aes(lon_bin, lat_bin, fill = sum_hours)) +
    coord_proj(xlim = c(-67, -63),ylim = c(0, 9))

是否有另一种简单的方法可以使用coord_proj生成放大的栅格地图并避免这种限制，即geom_raster仅适用于笛卡尔坐标？

我能想到的其他选项是为每个'放大'区域生成单独的shapefile，但我正在绘制许多这些栅格，而不是必须为每个栅格生成单独的shapefile，而是使用coord_proj以编程方式指定栅格地图限制。

由于

Answer 1

我认为您需要使用geom_tile()代替geom_raster()。 geom_raster()内部使用rasterGrob，这是一个只能线性缩放的位图。因此限制了笛卡尔坐标系。 geom_tile()绘制单个矩形，可以将其转换为任何坐标系。

我没有您的数据集，但我可以展示一个非常简单的示例：

df <- data.frame(x = 1:100) # a very simple dataset

p_raster <- ggplot(df, aes(x, fill = x, y = 0)) + 
  geom_raster() + 
  scale_fill_distiller()

p_raster

p_raster + coord_polar()
## Error: geom_raster only works with Cartesian coordinates

现在使用geom_tile()：

# for geom_tile(), map both fill and color to avoid drawing artifacts
p_tile <- ggplot(df, aes(x, color = x, fill = x, y = 0)) + 
  geom_tile() + 
  scale_fill_distiller() +
  scale_color_distiller()

p_tile

p_tile + coord_polar()