Re: [R] Error in R_nc4_open

2017-01-11 Thread Debasish Pai Mazumder
;https://nomads.ncdc.noaa.gov/thredds/catalog.html > >and make sure that the data path hasn't changed. > > > > > >Sarah > > > >On Wed, Jan 11, 2017 at 4:46 PM, Debasish Pai Mazumder > ><pai1...@gmail.com> wrote: > >> Hi Jeff, > >>

Re: [R] Error in R_nc4_open

2017-01-11 Thread Debasish Pai Mazumder
ugh RStudio and confirm > the results (if it doesn't happen in R then RStudio may be breaking R > (rare, but it has happened) in which case you would have to ask them for > help. > -- > Sent from my phone. Please excuse my brevity. > > On January 11, 2017 12:32:09 PM PST, Debasi

[R] Error in R_nc4_open

2017-01-11 Thread Debasish Pai Mazumder
Hi all, I recently updated my Rstudio to the newer version (Version 1.0.136) and I started to getting following error when I am trying to read netcdf files Error in R_nc4_open: NetCDF: DAP server error Any ideas? with regards -Deb [[alternative HTML version deleted]]

[R] download multiple files in one file

2016-10-13 Thread Debasish Pai Mazumder
Hi All, I have downloaded single file from the web using following scripts gribfile<-"

Re: [R] OPeNDAP access / OPeNDAP subsetting with R

2016-10-10 Thread Debasish Pai Mazumder
. ..$ id : int 65536 > > .. ..$ name : chr "" > > .. ..$ ndims: int 4 > > .. ..$ nvars: int 7 > > .. ..$ natts: int 13 > > .. ..$ dimid: int [1:4(1d)] 0 1 2 3 > > .. ..$ fqgn : chr "" > > .. ..- attr(*, "class")= chr &qu

Re: [R] OPeNDAP access / OPeNDAP subsetting with R

2016-09-30 Thread Debasish Pai Mazumder
49 Roy Mendelssohn - NOAA Federal < > roy.mendelss...@noaa.gov> wrote: > >> Please post the code of what you tried, as I have no idea otherwise what >> did or did not work for you. >> >> -Roy >> >> > On Sep 27, 2016, at 12:44 PM, Debasish Pai Ma

Re: [R] OPeNDAP access / OPeNDAP subsetting with R

2016-09-27 Thread Debasish Pai Mazumder
27, 2016 at 11:51 AM, Roy Mendelssohn - NOAA Federal < roy.mendelss...@noaa.gov> wrote: > Look at the package ncdf4. You can use an OPeNDAP URL in place of the > file name to perform subsets., > > -Roy > > > On Sep 27, 2016, at 9:06 AM, Debasish Pai Mazumder

[R] OPeNDAP access / OPeNDAP subsetting with R

2016-09-27 Thread Debasish Pai Mazumder
Hi all, I would like to access and subset following OpeNDAP files. server: http://nomads.ncdc.noaa.gov/thredds/dodsC/modeldata/cfsv2_forecast_ts_9mon/2014/201404/20140403/2014040312/ file name: tmax.01.2014040312.daily.grb2

Re: [R] How to read a grib2 file

2016-09-13 Thread Debasish Pai Mazumder
sum...@gmail.com> wrote: > What is your computer system? What is the output of this? > > sessionInfo() > > If you point to a file I'll try it so I can tell you the minimum system > requirements. > > Cheers, Mike > > On Tue, 13 Sep 2016, 08:56 Debasish Pai Mazumder

Re: [R] How to read a grib2 file

2016-09-12 Thread Debasish Pai Mazumder
<mdsum...@gmail.com> wrote: > > > On Sat, 10 Sep 2016 at 07:12 Debasish Pai Mazumder <pai1...@gmail.com> > wrote: > >> Hi >> I am trying to read a grib2 file in R. >> >> Here is my script >> >> library(rgdal) >> library(sp) >

[R] How to read a grib2 file

2016-09-09 Thread Debasish Pai Mazumder
Hi I am trying to read a grib2 file in R. Here is my script library(rgdal) library(sp) library(rNOMADS) gribfile<-"tmax.01.2011040100.daily.grb2" grib <- readGDAL(gribfile) I am getting following error : dec_jpeg2000: Unable to open JPEG2000 image within GRIB file. Is the JPEG2000 driver

Re: [R] How to make the "apply" faster

2016-07-10 Thread Debasish Pai Mazumder
n argument > instead of > fixing it at 2. E.g. > >> f3 <- function (condition, spellLength = 2) >{ >stopifnot(is.logical(condition), !anyNA(condition)) > y <- rle(condition) >sum(y$lengths[y$values] >= spellLength) >} >> f

Re: [R] How to make the "apply" faster

2016-07-10 Thread Debasish Pai Mazumder
, 2016 at 3:46 PM, Charles C. Berry <ccbe...@ucsd.edu> wrote: > On Sat, 9 Jul 2016, Debasish Pai Mazumder wrote: > > I have 4-dimension array x(lat,lon,time,var) >> >> I am using "apply" to calculate over time >> new = apply(x,c(1,2,4),FUN=function(y) {le

[R] How to make the "apply" faster

2016-07-09 Thread Debasish Pai Mazumder
I have 4-dimension array x(lat,lon,time,var) I am using "apply" to calculate over time new = apply(x,c(1,2,4),FUN=function(y) {length(which(y>=70))}) This is very slow. Is there anyway make it faster? -Debasish [[alternative HTML version deleted]]

[R] Static to interactive map (leaflet spplot)

2016-02-02 Thread Debasish Pai Mazumder
Hi ALL, I have a script to plot hexagonal polygon on a map. Its a static map. I used spplot. I would like to convert this plot to interactive plot using " *leaflet*". How do I make it interactive plot? Here is spplot lines: cl = map("world", xlim = c(-120, 20), ylim = c(-10, 70), plot = TRUE)

[R] Problem with installing "mapview"

2016-01-25 Thread Debasish Pai Mazumder
Hi all, I am using Rstudio version 0.99.491. I am trying install "mapview" package. install.packages("mapview") I receive following error message configure: error: Cannot compile a simple JNI program. See config.log for details. Make sure you have Java Development Kit installed and correctly

[R] paneling spplot and hist

2015-12-14 Thread Debasish Pai Mazumder
Hi all, I am new in R. I am trying to panel spplot and hist. spplot(hspdf, "CDP", col = "white", col.regions = blue2red(20), sp.layout = list(l2), at = seq(1,10,1), colorkey = list(space = "bottom", labels = list(labels = paste(seq(1,10,1)), cex = 1.5)), sub = list("CDP", cex = 1.5, font = 2))

[R] problem with spplot

2015-12-11 Thread Debasish Pai Mazumder
Hi all, I am using spplot for a spatial map. spplot(hspdf, "CDP", col = "white", col.regions = blue2red(20), sp.layout = list(l2), at = seq(1,10,1), colorkey = list(space = "bottom", labels = list(labels = paste(seq(1,10,1)), cex = 1.5)), sub = list("CDP", cex = 1.5, font = 2)) I have three

[R] spplot

2015-12-10 Thread Debasish Pai Mazumder
Hi all, I am using spplot for a spatial map. spplot(hspdf, "CDP", col = "white", col.regions = blue2red(20), sp.layout = list(l2), at = seq(1,10,1), colorkey = list(space = "bottom", labels = list(labels = paste(seq(1,10,1)), cex = 1.5)), sub = list("CDP", cex = 1.5, font = 2)) I have three