小编del*_*aye的帖子

我有2张图,一张用ggplot2这样绘制的地图:

w<-ggplot()+
  geom_polygon(data=dep_shp.df, aes(x=long,y=lat,group=group,fill=classJenks))+

  #   scale_fill_gradient(limits=c(40, 100))+
  labs(title ="Classification de la proportion de producteurs par départements
       \n par la methode de jenks (2008)")+
  theme_bw()+
  coord_equal()

和图形类型的对象classIntervals 从classInt库.

我想把这两张图组合在一起.我试过了:

vplayout <- function(x, y) viewport(layout.pos.row = x, layout.pos.col = y)
grid.newpage()
pushViewport(viewport(layout = grid.layout(1, 2)))

#creation
print(u, vp = vplayout(1, 1))
print(v, vp = vplayout(1, 2))

还有一些东西 grid.arrange

grid.arrange(plot1, plot2, ncol=2)

但这些都不起作用.

我有一个SpatialPointsDataFrame加载

pst<-readOGR("/data_spatial/coast/","points_coast")

我想在输出中得到一个SpatialLines,我发现了一些东西

coord<-as.data.frame(coordinates(pst))
Slo1<-Line(coord)

Sli1<-Lines(list(Slo1),ID="coastLine")
coastline <- SpatialLines(list(Sli1))
class(coastline)

它似乎工作,但当我尝试绘图(海岸线),我有一条不应该在那里的线...

有人能帮帮我吗？shapefile就在这里!

我用来geom_tile()在同一张图上绘制 3 个变量...

tile_ruined_coop<-ggplot(data=df.1[sel1,])+
  geom_tile(aes(x=bonus, y=malus, fill=rf/300))+
  scale_fill_gradient(name="vr")+
  facet_grid(Seuil_out_coop_i ~ nb_coop_init)
tile_ruined_coop

我对结果很满意！

但是应用什么样的统计处理fill呢？这是一个意思吗？

在RI中有一个SpatialPointsDataFrame whit复制点(坐标和属性),我想用相同的数据删除所有点...

我在sp包中找到了这个remove.duplicates()函数,但它似乎只在位置上删除了...还有另一种方法吗？

谢谢

E.

我有类似的东西:

test[1] 
"[0 30.5 4.5 10.5 2 35 22.999999999999996 29 5.500000000000001 23.5 18 23.5 44.5 3 44.5 44.00000000000001 43 27 42 35.5 19.5 44.00000000000001 1 0 31 34 18 1.5 26 6 45.99999999999999 10.5 9.5 24 20 42.5 14.5 45.5 20.499999999999996 150 45.5 0 4.5 22.5 4 9 8 0 0 15.5 30.5 7 5.500000000000001 12.5 33.5 15 500 22.5 18 43 4.5 26 23.5 16 4.5 7.5 32 0 0 18.5 33 31 14.5 21.5 0 40 0 0 …

我试图在 python 中平铺一个 JPG 图像。通常我使用 imageMagick .. 所以我看到比魔杖似乎做这项工作......

但我不会翻译

 convert -crop 256x256 +repage big_image.jpg tiles_%d.jpg

有人可以帮助我吗？

我想在 2014 年和 2017 年 12 月之间将Landsat 8 32-Day NDVI Composite 保留imageCollection一层。有没有办法做一个可变过滤器？还是按月过滤？