Построение данных из netCDF в ggplot для R - PullRequest
Новая
       22

Построение данных из netCDF в ggplot для R

0 голосов
/ 08 февраля 2019

Я рассчитываю построить данные из netCDF из следующего источника:

http://fizz.phys.dal.ca/~atmos/martin/?page_id=140

Просто несколько простых, но больших файлов с концентрациями PM 2.5, широтами и долготами.

Я хочу нанести эти данные на карту Массачусетса.

В настоящее время я получаю следующую ошибку: 

Error: Aesthetics must be either length 1 or the same as the data (286): x

И следующеемой код: 

library(chron)
library(RColorBrewer)
library(lattice)
library(ncdf4)
library(raster)
library(ggplot2)
library(ggmap)
library(maps)
library(mapdata)
memory.limit(size=56000)
getwd()
wd <- "C:/Users/tadams15/Desktop/Practicum/vd"
setwd(wd)
file <- "1adjust.nc"
ncin <- nc_open(file)
print(ncin) # Displays plenty of information about netCDF
lat <-ncvar_get(ncin,"LAT")
lon <-ncvar_get(ncin,"LON")
pm <- ncvar_get(ncin,"PM25")
states <- map_data("state")
Mass <- subset(states,region == 'massachusetts')
ggplot(data = Mass)
image(pm)
lat.df <- data.frame(lat)
lon.df <- data.frame(lon)
pm.df<- data.frame(pm)


m <- map_data('state',region='Massachusetts')

graph <-ggplot() + 
  geom_polygon( data=m, aes(x=lon, y=lat,group=group),colour="black", fill="white" )+
  geom_point(data=pm,aes(x=lon,y=lat),,colour="red",)+
  ggtitle("Distribution of Flash Flood Events in CONUS")+
  xlab('Longitude')+
  ylab('Latitude')+
  coord_fixed()

print(graph)

Мне было интересно, доволен ли кто-нибудь функцией ggplot и можно ли предложить какую-либо помощь.Я очень плохо знаком с R и даже новичок в управлении типами файлов netCDF.Любая помощь будет принята с благодарностью!

...