小编Mun*_*ish的帖子

我正在使用ggplot(包名:ggplot2)绘制15行,每个行代表一个单独的实体,并希望为它创建一个图例.但是,我无法将图例条目划分为多个列.有人可以建议如何在ggplot环境中做同样的事情.

目前,我使用以下命令创建图例:

opts(title=plotName,legend.position='bottom')

但是,这给出了一个列的图例.因此,图表中的大面积区域由传说本身占据.将它分成2或3列将真正有助于原因,同时将图例保持在图表的底部.我也试过了,legend.direction但是这个命令在一行中显示了图例,这是不可取的,除非我可以将它传播到2-3行.

opts(title=plotName,legend.position='bottom',legend.direction="horizontal")

提前谢谢,Munish

我正在读取R中的.tif文件并获得下面列出的4条警告消息.当我按照第4条消息的说明进行操作时,前3个警告仍然存在,但从文件读取的值在每个像素处都会发生剧烈变化.请帮助我从.tif文件中正确读取数据.示例文件可以在链接上找到:ftp://ftp.ntsg.umt.edu/pub/MODIS/NTSG_Products/MOD16/MOD16A2_MONTHLY.MERRA_GMAO_1kmALB/GEOTIFF_0.05degree/

我的代码:

remove(list=ls()) 

library(tiff)

library(raster)

str_name<-'MOD16A2_ET_0.05deg_GEO_2008M01.tif' 

read_file<-readTIFF(str_name) 

警告信息:

1: In readTIFF(str_name) :
  TIFFReadDirectory: Unknown field with tag 33550 (0x830e) encountered
2: In readTIFF(str_name) :
  TIFFReadDirectory: Unknown field with tag 33922 (0x8482) encountered
3: In readTIFF(str_name) :
  TIFFReadDirectory: Unknown field with tag 34735 (0x87af) encountered
4: In readTIFF(str_name) :
  tiff package currently only supports unsigned integer or float sample formats in direct mode, but the image contains signed integer format - it will be treated as unsigned (use …

[

注意两个图像上的颜色比例标签差异.顶部图像是R 3.2.5的输出,相同的代码在R 2.15.2中生成底部图像.我想解决顶部图中标签的"移位"以匹配底部标签.此查询中给出的示例代码用于生成两个图.

我正在尝试使用image.plot绘制地图,但是颜色刻度的最小值和最大值并未完全显示在提示处.我在R版本3.2.5(2016-04-14)中遇到了这个问题; 平台:x86_64-w64-mingw32/x64(64位)和库包'字段'垃圾邮件版本1.4-0(2016-08-29)

相反,相同的命令可以在R版本2.15.2(2012-10-26),平台:x86_64-w64-mingw32/x64(64位)中正确显示色标边缘的最小值和最大值. 'fields'包0.41-0(2014-02-26).这是示例代码:

library(fields)
temp <- matrix(data=rexp(200, rate=10), nrow=180, ncol=360)
min(temp)
max(temp)
color_plate <- c("#FF0000", "#FF4D00", "#FF7000", "#FF8A00", "#FFA800", "#FFBF00", "#FFF000", "#FFFF54", "#AAFFFF","#7FFFFF", "#55FFFF", "#2AFFFF", "#00CFFF", "#20BFFF", "#209FFF", "#2060FF")

zlim <- seq(0.08,0.40,by=0.04)
temp[temp<min(zlim)] <- min(zlim)
temp[temp>max(zlim)] <- max(zlim)

image.plot(temp,col=color_plate, 
           axis.args=list(cex.axis =1,at=zlim, labels=zlim,mgp=c(1, 0, 0),tck=0.1))

我希望在将矩阵转换为栅格对象(tempMap)后绘制矩阵(temp_matrix).此外,我希望在同一个绘图窗口中添加几个可以获得纬度和经度位置的点.我尝试过几种方法,但似乎没有一种方法可行,因为可用点是lat/long中的特定位置,而我得到的栅格对象具有不同的范围.请帮我解决这个问题.下面给出了问题的示例数据.

library(raster)
temp_matrix<-array(NA,c(11,11))
temp_matrix[1,]<-c(NA,NA,NA,NA,NA,NA,NA,0,0,-6,-6)
temp_matrix[2,]<-c(0,0,0,0,NA,NA,1,0,0,0,0)
temp_matrix[3,]<-c(1,0,0,-1,-1,0,0,0,1,0,0)
temp_matrix[4,]<-c(1,1,0,0,0,0,-1,-1,0,0,0)
temp_matrix[5,]<-c(1,NA,NA,NA,NA,-1,-1,-1,0,-1,-1)
temp_matrix[6,]<-c(NA,NA,NA,NA,NA,NA,-1,-1,-1,0,0)
temp_matrix[7,]<-c(NA,NA,NA,NA,NA,NA,NA,0,-1,0,0)
temp_matrix[8,]<-c( NA,NA,NA,NA,NA,NA,NA,0,0,0,-1)
temp_matrix[9,]<-c(NA,NA,NA,NA,NA,NA,NA,-1,0,0,0)
temp_matrix[10,]<-c(NA,NA,NA,NA,NA,NA,NA,NA,-1,-1,-2)
temp_matrix[11,]<-c(NA,NA,NA,NA,NA,NA,NA,NA,-2,-3,-2)
plot(raster(temp_matrix),axes = FALSE,legend=FALSE)
tempMap <- raster(temp_matrix)

# plot the points over this raster layer
point_1<-c(10,10) # should appear on 2nd row from top i.e. over temp_matrix[2,10]
point_2<-c(9,10)  # should appear on 3rd row from top i.e. over temp_matrix[3,10]
point_3<-c(1,10)  # should appear on lowermost row i.e.over temp_matrix[11,10]

我希望同时包含一个下标和一个标题,但是我没有得到以下命令所需的结果:

base<-'B1'
compare<-'A1'
plot (1, 1, main = bquote('Annual mean' ~CO[2] ~'Flux Difference: \n'  ~.(compare)~ 'minus'~.(base)))

我想在'Flux Difference'之后插入一个换行符,但它没有发生,而是这个术语无处可去.请帮我.

谢谢,Munish