python/readslice.R

   1 # Load required packages
   2 library(ggplot2)
   3 library(raster)
   4 library(ncdf4)
   5
   6 ncfile <- nc_open("smoke.nc")
   7
   8 # Read the variables from the NetCDF file
   9 pm25 <- ncvar_get(ncfile, "pm25")
  10 lon <- ncvar_get(ncfile, "XLONG")
  11 lat <- ncvar_get(ncfile, "XLAT")
  12
  13 # Close the NetCDF file
  14 nc_close(ncfile)
  15
  16 print("Create a data frame with the latitude, longitude, and pm25 values")
  17 df <- data.frame(lat = as.vector(lat), lon = as.vector(lon), pm25 = as.vector(pm25))
  18
  19 print("Create a scatter plot of the pm25 values at the lat/lon points")
  20
  21 dev.new()
  22 plot(df$lon, df$lat, col = df$pm25)
  23
  24