[问题类型]:
程式咨询(我想用R 做某件事情,但是我不知道要怎么用R 写出来)
[软件熟悉度]:
入门(写过其他程式,只是对语法不熟悉)
[问题叙述]:
各位大大好,小弟我想要将多年(10年以上)的区域(东亚)降雨资料读进R,并取出多年的
降雨值来做计算处理,而目前资料可以读的进R,但在取降雨这个变量的时候,R就跑出
内存不足的错误讯息,降雨值的资料为一矩阵(经度x纬度),经度总共360格,纬度共
280格(360*280),每一格皆有一降雨值,一天一个矩阵(平面),所以我在想有没有可能
真的是要取出来的降雨这个变量资料量太占内存,而造成没办法取出这么多年分的降
雨值出来,小弟我目前试我的电脑最多只能取出4年的降雨值出来,所以想请教各位大大
有没有什么方法可以解决吗??谢谢。
资料可从这里取得:http://0rz.tw/JI056 (APHRO这个档案)
[程式范例]:
rm(list=ls())
library(ncdf4)
Aphrodite_data <- "C:\\Users\\TOM\\Desktop\\R(数据库)\\Aphrodite资料
\\APHRO_MA_025deg_V1101R2.1961-2007\\APHRO_MA_025deg_V1101R2.1961-2007.nc"
obs <- nc_open(Aphrodite_data)
print(obs)
lon <- ncvar_get(obs,"lon")
lat <- ncvar_get(obs,"lat")
level <- ncvar_get(obs,"lev")
time <- ncvar_get(obs,"time")
tunits <- ncatt_get(obs, "time", "units")
precip <- ncvar_get(obs, "pr")
precip <- ncvar_get(obs, "pr",start=c(1,1,1,1),count=c(360,280,1,1460))
dim(precip)
nc_close(obs)
[环境叙述]:
https://imgur.com/kJiIoc7
[关键字]:
多年资料 变量 nc档