[R-sig-Geo] Subset a Raster* object read from a 4D NetCDF file

[Dutrieux, Loic]

Dear Pascal,

It appears that you can't do this type of operations if the file is only open through a remote connection.
You may try:

library(raster)
out <- brick('http://www.esrl.noaa.gov/psd/thredds/dodsC/Datasets/godas/pottmp.2013.nc',
lvar=3, level=1, varname='pottmp')
out2 <- readAll(out)

# subset
pottmp <- subset(out2, 1:3)

Best regards,
Loïc
________________________________________
From: r-sig-geo-bounces at r-project.org [r-sig-geo-bounces at r-project.org] on behalf of Pascal Oettli [kridox at ymail.com]
Sent: Monday, April 14, 2014 10:25 AM
To: r-sig-geo at r-project.org
Subject: [R-sig-Geo] Subset a Raster* object read from a 4D NetCDF file

Dear list members,

Could someone help me to understand the following case (works on
Linux. I don't know for other OS):

library(raster)
out <- brick('http://www.esrl.noaa.gov/psd/thredds/dodsC/Datasets/godas/pottmp.2013.nc',
lvar=3, level=1, varname='pottmp')

# error while subsetting
pottmp <- subset(out, 1:3)

# multiply by 1
out <- out * 1

# no error
pottmp <- subset(out, 1:3)


Thank you,
Pascal

_______________________________________________
R-sig-Geo mailing list
R-sig-Geo at r-project.org
https://stat.ethz.ch/mailman/listinfo/r-sig-geo