Re: [R-sig-Geo] Gridding a daily time series in R

2019-09-12 Thread Roy Mendelssohn - NOAA Federal via R-sig-Geo
> On Sep 12, 2019, at 9:03 AM, Sarah Goslee wrote: > > Creating the NetCDF file is easy - there are multiple packages to do that. Can I just gently amend this statement. Writing an arbitrary netCDF file is easy, writing a useful netCDF file is hard. The difference is the first just

Re: [R-sig-Geo] Area calculation of forest cover classes.

2019-09-12 Thread Sarah Goslee
Hi, Without knowing any details of what you did, the general procedure is: Use your fitted Random Forest or SVM model to predict the class for each pixel in your region of interest - the predict() function. Use standard spatial data methods to aggregate the resulting spatial data, or even

Re: [R-sig-Geo] Gridding a daily time series in R

2019-09-12 Thread Sarah Goslee
Hi, Creating the NetCDF file is easy - there are multiple packages to do that. Everything else you ask about is hard, and not really R questions. You need to know: What is the best way to impute missing precipitation data to fit my needs? What is the best way to grid point-based precipitation

[R-sig-Geo] Gridding a daily time series in R

2019-09-12 Thread Cristo Facundo PĂ©rez
Dear community, I have a daily time series of precipitation, which I intend to transform it into a single NetCDF file of daily precipitation. So, first, I would like to explore the best way to impute missing values, grid the daily values of the different available weather stations and get daily

Re: [R-sig-Geo] rgdal, PROJ6 and "+init=epsg" syntax on Jupyter/conda

2019-09-12 Thread James Sample
After going round in circles most of the day, I finally discovered that the PROJ_LIB environment variable within my Dockerfile did not match my conda environment. When I loaded rgdal I saw Path to PROJ.4 shared files: (autodetected) in the output and assumed all was OK, but having explicitly

Re: [R-sig-Geo] rgdal, PROJ6 and "+init=epsg" syntax on Jupyter/conda

2019-09-12 Thread Edzer Pebesma
CRAN, btw, reports similar errors on all debian testing platforms, of the kind > st_crs("+init=epsg:3857")$epsg proj_create: init=epsg:/init=IGNF: syntax not supported in non-PROJ4 emulation mode see

Re: [R-sig-Geo] rgdal, PROJ6 and "+init=epsg" syntax on Jupyter/conda

2019-09-12 Thread Edzer Pebesma
On 9/12/19 10:55 AM, James Sample wrote: > Thanks Edzer - that's very helpful! > > Would it be possible for you to link/share your Ubuntu Dockerfile, please? > (I completely understand if you'd rather not, of course). https://github.com/r-spatial/sf/tree/master/inst/docker/gdal > > Perhaps

Re: [R-sig-Geo] rgdal, PROJ6 and "+init=epsg" syntax on Jupyter/conda

2019-09-12 Thread James Sample
Thanks Edzer - that's very helpful! Would it be possible for you to link/share your Ubuntu Dockerfile, please? (I completely understand if you'd rather not, of course). Perhaps if I can see how you're setting thing up in Ubuntu I can modify my Jupyter Dockerfile accordingly, and if I can

[R-sig-Geo] Area calculation of forest cover classes.

2019-09-12 Thread Enoch Gyamfi Ampadu
Dear List, Please I have carried out an RF and SVM forest cover classification in R. I want to determine the area in hectors for each of the forest cover classes. I have not been able to find my way out on it. I wanted to ask if it is possible to do this in R and will be glad to have assistance