Hi, I have downloaded a NetCDF file from https://psl.noaa.gov/data/gridded/data.gpcp.html which shows mean monthly precipitation between 1979-2020. I think there are around 501 time slices in total each representing one month.
I want to average all the data to show the mean precipitation over the Pacific Ocean between 1979-2020. I am open to doing this in any possible way (ArcMap, R, Python).
First I tried to do this in ArcMap, I plotted the NetCDF file as a raster using "Make NetCDF Raster Layer". This worked but will only show one time slice at a time, I could not find any option to average the time slices.
Then I tried to do this in R but also got stuck, my code below shows how far I got. I have read that it is also possible to average the data in Python but I am new to Python and wouldn't know where to start.
Thanks for your help!
setwd("C:/Users/nicho/OneDrive/Documents/ArcGIS")
library(ncdf4)
library(raster)
library(rgdal)
library(ggplot2)
# Open a link to the data in my working directory
nc_data
I got a couple of R scripts that can easily do what you are looking for. Email me and I will send you one of them. [email protected] !!
Dear Nichola,
instead of package 'raster' I recommend you to use package 'stars' that can handle spatiotemporal rasters, such as your NetCDF dataset. Functions read_stars() and read_ncdf() can read the NetCDF file, while st_apply(, FUN = mean) will do the averaging. You'll find it really easy.
HTH,
Ákos
you can do it in R studio very easily
visit the following link below
https://support.esri.com/en/technical-article/000011318
- Badges
- Science topic
- Similar topics
- Computer Science and Engineering
- Programming Languages
- Python