如何在R中绘制地理参考数据集?

极大主义

我有要在美国地图上绘制的栅格数据:https : //www.dropbox.com/s/9khcjgtv8ipo2u5/windspeed.txt?dl=0

library(ggplot2)
library(RColorBrewer)
library(rgdal)
library(sp)
library(maps)
options(max.print=5.5E5) 

all_data = read.table("windspeed.txt",header = TRUE)

res=0.01 #spacing of row and col coords pre-specified
origin_lat_lon=c(24.55, -130) 
all_data$row=(all_data$row)*res+origin_lat_lon[1] 
all_data$col=(all_data$col)*res+origin_lat_lon[2]
coords = cbind(all_data$col, all_data$row)
spdf = SpatialPointsDataFrame(coords, data=all_data) #sp = SpatialPoints(coords)
proj4string(spdf) <- CRS("+init=epsg:4269") 

df=as.data.frame(spdf)
myPalette <- colorRampPalette(rev(brewer.pal(10, "Spectral")))
usamap <- map_data("state")
ggplot(data=df,aes(x=col,y=row,color=m)) + 
  geom_polygon( data=usamap, aes(x=long, y=lat,group=group),colour="black", fill="white" )+
  geom_point()+
  scale_colour_gradientn(name = "Wind",colours = myPalette(10), limits=c(0,1))+
  xlab('Longitude')+
  ylab('Latitude')+
  theme_bw()+
  theme(line = element_blank())+
  theme(legend.position = c(.93,.20),panel.grid.major = element_line(colour = "#854440"))+
  ggsave("test.png",width=10, height=8,dpi=300)

但是我得到了一个倒置的情节。你能帮忙吗?

在此处输入图片说明

我以前在这里得到了类似数据集的答案:如何将以网格间隔收集的点数据转换为r中的地理参考数据集?

卢卡斯·福特尼(Lucas Fortini)

与您先前提到的上一个问题所获得的数据集相比,此csv中的纬度值将被还原。您所需要做的就是反转这个新数据集中的行号:

在您的行之后:

all_data = read.table("windspeed.txt",header = TRUE)

使用以下方法反转行号:

max_row= max(all_data$row)
all_data$row=max_row-all_data$row 

那应该照顾它。

本文收集自互联网,转载请注明来源。

如有侵权,请联系[email protected] 删除。

编辑于
0

我来说两句

0条评论
登录后参与评论

相关文章

来自分类Dev

如何在R中绘制“组”数据?

来自分类Dev

如何在R中绘制大量数据

来自分类Dev

如何在R中绘制纵向数据

来自分类Dev

如何在R中制作地理数据热图?

来自分类Dev

如何在OpenGL / C ++中绘制多维数据集?

来自分类Dev

将地理参考数据与R中的形状文件匹配

来自分类Dev

将地理参考数据与R中的形状文件进行匹配

来自分类Dev

如何在R中订购数据集

来自分类Dev

如何在R中查找数据集

来自分类Dev

如何在R中查找数据集

来自分类Dev

R Leaflet:如何在 addTiles 函数中绘制具有特定地理位置的地图?

来自分类Dev

在R中绘制多个数据集

来自分类Dev

在R中绘制多个数据集

来自分类Dev

如何在R Studio中从数据集反转数据?

来自分类Dev

如何在Gnuplot中绘制具有透明填充样式的数据集

来自分类Dev

如何在Gnuplot中绘制具有透明填充样式的数据集

来自分类Dev

如何在R中读取大型数据集的子集?

来自分类Dev

如何在R中还原内置数据集?

来自分类Dev

如何在R中操纵配对数据集?

来自分类Dev

如何在R包中包含数据集?

来自分类Dev

如何在不同地块上绘制数据集

来自分类Dev

如何在R脚本中从数据集1创建子集(数据集2)

来自分类Dev

如何在MPAndroidChart中从右到左绘制数据?

来自分类Dev

如何在python中绘制“历史”数据

来自分类Dev

如何在Stata中绘制数据子集?

来自分类Dev

如何在Matlab中从csv绘制数据

来自分类Dev

在R中绘制大小不同的独立数据集

来自分类Dev

如何在R中没有未选择数据的情况下绘制子集数据

来自分类Dev

如何绘制存储在 Pandas DataFrame 中的垂直数据集?