search for: alyaari

I have this https://echange-fichiers.inra.fr/get?k=6TdTdezNeZwAqJtpwVm grib file which is a regular lat-lon 0.25x0.25 degrees and contains 4 fields . I tried to open it using : grib <- readGDAL("C:\\Users\\aalyaari\\Desktop\\Vol_025_H14_2010060700.grib") but I got this error: C:\Users\aalyaari\Desktop\Vol_025_H14_2010060700.grib has GDAL driver GRIB and has 721 rows and 1440 columns Error in validityMethod(as(object, superClass)) : Geographical CRS given to non-conformant data: -90.125 90.125 Any...

This code will read binary file and display it as a map. may problem is that this code is using a continuous colour scheme, even though I have discrete data (which is a classification scheme). How can I map numbers to colours with raster? Please require(raster) conne <- file("C:\\lai.bin", "rb") sd<- readBin(conne, integer(), size=1, n=360*720, signed=F)

Hi I am using windows. I cant install ncdf4 package but it didn't work . any suggestions!! -- View this message in context: http://r.789695.n4.nabble.com/ncd4-package-tp4477496p4477496.html Sent from the R help mailing list archive at Nabble.com.

I have two binary files(rasters) with the same dimensions. The first file is called `over` and the second is `corr`. I want to replace values in `over` by `NA` whenever `corr` is greater than 0.5. to read the two files we can use: conne <- file("C:corr.bin","rb") over <- readBin(conne, numeric(), size=4, n=1440*720, signed=TRUE) frf <-

Suppose we have files in one folder file1.bin, file2.bin, ... , and 1460slice(file) with dim of 720 * 360 in directory C:\\PHD\\Climate Data\\Wind\\ and we want to read them and make a loop to go from 1 to 4 and take the average, then from 4 to 8 and so on till 1460. in the end we will get 365 files . I need those 365 files to be in one new folder for later use in my model I tried using this

...ntical. - my other data are in WGS84 with 0.25*0.25 degree resolution and extent -180, 180, -90, 90. So I want to re-project the data I have here from EASE grid, cylindrical,25 km ,1383 pixel 586 lines to WGS84,0.25*0.25 degree,1440 pixel 720 lines: f=open.ncdf("C:\\Users\\aalyaari\\Desktop\\SM_RE01_MIR_CLF31D_20100812T000000_20100812T235959_246_001_7.DBL") A = get.var.ncdf(nc=f,varid="Soil_Moisture",verbose=TRUE) I tried this: d <-raster("C:\\Users\\aalyaari\\Desktop\...

Hellow everyone, This code bellow will calculate average daily wind speed(measurements are taken every three hours).Any ideas how to take the Min and Max instead of average. library(Matrix) setwd("C:\\Users\\aalyaari\\Desktop\\img") listfile<-dir() long <- file("C:\\Users\\aalyaari\\Desktop\\New folder (5)\\inra.bin", "rb") A=readBin(long, integer(), size=2,n=67420*1, signed=F) ta<-t(A) lot <- file("C:\\Users\\aalyaari\\Desktop\\New folder (5)\\lat.img", "r...

As you can see: time is from 1 to 1460, 1 means they the measure Tem for the globe at the first 6 hours , 2=after 12 hours .....; and so on 1460 *6 = 8760 hours which equals 1 year first of all I want to convert all data from kelvin to degree second I want to convert from 6 hourly to daily and make it like this raw data what i want to make time

I faced this problem when typing: temperature <- get.var.ncdf( ex.nc, 'Temperature' ) *unable to allocate a vector of size 2.8 GB* By the way my computer memory is 4G and the original size of the file is 1.4G,netcdf file I don't know what is the problem.Any suggestion please I tried also memory limit(4000) 4000 but didnt solve the problem.any help -- View this message in

I am getting theses values when I disply my ncdf file 1e+34 -1e+34 -1e+34 -1e+34 -1e+34 -1e+34 -1e+34 -1e+34 [227,] -1e+34 -1e+34 -1e+34 -1e+34 -1e+34 -1e+34 -1e+34 -1e+34 -1e+34 -1e+34 [228,] -1e+34 -1e+34 -1e+34 -1e+34 -1e+34 -1e+34 -1e+34 -1e+34 -1e+34 -1e+34 [229,] -1e+34 -1e+34 -1e+34 -1e+34 -1e+34 -1e+34 -1e+34 -1e+34 -1e+34 -1e+34 [230,] -1e+34 -1e+34 -1e+34 -1e+34 -1e+34 -1e+34 -1e+34

I have 3 files A ,B ,C as binary files. C represents the values of tepm measured every 3 hours for one month so it has 674200 columns and 248 rows.values of A represent the correspond lines(these values could be within 1to360) and , while values of B represent the correspond pixel(1to720) for each point in file C(so both they have 1 row and 67420rows).Now I will create a new matrix ,M ,which has

Has one try to install the ggplot2 package recently? I tried to install it on my new system and had trouble: I tried a different CRAN mirror but didn't work library(ggplot2) Error in loadNamespace(i, c(lib.loc, .libPaths())) : there is no package called ?stringr? In addition: Warning message: package ?ggplot2? was built under R version 2.15.2 Error: package/namespace load failed for

Dear R users, I have been searching thru the documentation of R but did not find anything about reading or manipulating a binary data block file (.DBL).a link to one example https://echange-fichiers.inra.fr/get?k=5Hzw2B1wZDng9ztO34E -- View this message in context: http://r.789695.n4.nabble.com/can-R-read-a-binary-data-block-file-DBL-tp4652654.html Sent from the R help mailing list archive

I have two directories https://echange-fichiers.inra.fr/get?k=AcHKdNI4No44GEsj7PK with 12 (global maps)binary files in each.I used the code given below to calculate the spatial correlation between these files and it worked well(the output is a global correlation map). I wonder if there is a simple way to calculate RMSE and bias along with cor.so finally we get three outputs (bias map,RMSE map,cor

I have 12 binary (raster) files https://echange-fichiers.inra.fr/get?k=k3M2jatJyHy65Cs99G4 . I would like to calculate the moving average for the 12 values for each pixel in the 12 files. For a simple vector we can get a moving average by using this : x <- c(1,2,3,NA,NA,4,6,5,6,4,2,5) movingmean <- rollapply(x, 3, FUN = mean, na.rm = T) now I want

I have a raster file(1440*720 rows) contains values of 1 ,2 , and 3. when I plot the file , I got a map of three colors but I do not know which is which. How can I put those colors as as I want : 1=red 2=blue 3=green code: pvm <- file("C:\\User_sm-das.bin","rb") cor1<- readBin(pvm, numeric(), size=4, n=1440*720, signed=TRUE) r <-raster(t(matrix((data=cor1),

I have 12 binary (raster) files https://echange-fichiers.inra.fr/get?k=k3M2jatJyHy65Cs99G4 . I would like to calculate the moving average for the 12 values for each pixel in the 12 files. For a simple vector we can get a moving average by using this : x <- c(1,2,3,NA,NA,4,6,5,6,4,2,5) movingmean <- rollapply(x, 3, FUN = mean, na.rm = T,fill=NA) now

I am using this code to calculate the moving average mean.It worked fine but when I wanted to also calculate based on sd(stander deviation) I got the error shown below. I read this documentation of R movingFun and found that sd was mentioned at http://hosho.ees.hokudai.ac.jp/~kubo/Rdoc/library/raster/html/movingFun.html I wonder if needs any other things to conseder or sd is not at all supplied

I have a binary file(a) with size of (360 720 )for the globe.I wrote the code given below to read and extract an area (south america)from that file. when I use summary for the whole file I got: summary(a, na.rm=FALSE) Min. 1st Qu. Median Mean 3rd Qu. Max. NA's 0.00 1.00 3.00 4.15 7.00 20.00 200083 . But when I used summary for the region(b) which I extracted, I got

I am using the code below to calculate the correlation map between two datasets. This code worked fine. dir1 <- list.files("D:thly", "*.bin", full.names = TRUE) dir2 <- list.files("D:002", "*.envi", full.names = TRUE) file_tot <- array(dim = c(1440, 720, 11, 2)) for(i in 1:length(dir1)) { file_tot[, , i, 1] <-