关于R语言计算和出图的问题

_梅里

楼主想做一个两个情景下的温度差值图。下面是代码：library(ncdf4)
library(fields)
library(mapproj)
nc=nc_open("O:\\igcm\\co_2\\igcm_cl_2000_01.nc")
my_temp_1=ncvar_get(nc,"surftemp")
nc=nc_open("O:\\igcm\\co_2\\igcm_cl_2000_02.nc")
my_temp_2=ncvar_get(nc,"surftemp")
lons=ncvar_get(nc,"longitude")
lats=ncvar_get(nc,"latitude")
nc_close(nc)
my_temp=my_temp_2-my_temp_1
my_temp_big <- array(0,c(65,32))
for (x in c(1:64)) my_temp[x,]=rev(my_temp[x,])
lats <- rev(lats)
for (x in c(1:32)) my_temp_big[1:65,x] <- c(my_temp[33:64,x],
                                            my_temp[1:32,x],my_temp[33,x])
lons[1:65] <- c(lons[33:64]-360,lons[1:32],lons[33])
image.plot(lons,lats,my_temp_big)
map(xlim=c(-180,180),ylim=c(-90,90),add=T,interior=F)
windows()
filled.contour(
  my_temp_big,x=lons,y=lats,color = tim.colors,
  xlab="longitude",ylab="latitude",
  plot.axes={
    map(
      xlim=c(-180,180),
      ylim=c(-90,90),
      add=T,
      interior=F
    )

因为直接用R出的图，侧面温标太难看了， 要对比很多图片的时候不统一，很容易产生误解，于是我就直接把 my_temp 写入一个新的nc文件里了，但是出的图不对，计算完的差值南北纬度正好反了，请问用什么方法，或代码能让数据正过来？ 后面的my_temp_big 多了一行经度的数据，所以不是很懂。。。完全没接触过R语言，问题略白痴，谢谢前辈指点

topmad

南北方向反了  
看两个地方 你read.nc的时候 是不是按照数据的方向读取的

或者 你绘图的时候 是不是搞反了  

R里面读取的时候 不会有啥问题

_梅里

topmad 发表于 2017-1-9 16:43
南北方向反了  
看两个地方 你read.nc的时候 是不是按照数据的方向读取的

谢谢，我用panoply查看的NC文件，我也是醉了，开始真的是反的，后来有莫名其妙的正常了。。。我什么都没改啊