Я не получаю значения "lat" и "long" для шейп-файла, который я использовал для маскировки файла netcdf в R.
Вот ссылка на данные: https://drive.google.com/file/d/15Feww0VvEq9SuJ4j4gxuVqvxpp_w8PVC/view?usp=sharing
Я прочитал как netcdf, так и шейп-файл моей области интереса в R. Я преобразовал проекцию шейп-файла в соответствие с проекцией файла netcdf. Я замаскировал файл netcdf границей шейп-файла. Я записываю значение переменной в CSV-файл. Хотя я получаю данные за весь период времени, я не вижу значений широты и долготы точек сетки. пожалуйста, предложите путь вперед. Спасибо
library(ncdf4)
library(ncdf4.helpers)
library(PCICt)
library(rgdal)
NC = brick("Daily_Pr_Inmm_Souris_1979_2016.nc")
SHP=readOGR("Souris_basin.shp")
crs(NC)
crs(SHP)
SHP=spTransform(SHP, crs(NC))
Souris_Pr_Daily=mask(NC, SHP)
tab=as.data.frame(na.omit(Souris_Pr_Daily[],xy=TRUE))
write.csv(tab, "Souris.csv")
Я должен увидеть что-то подобное в моем CSV-файле, где Var 1 и Var 2 представляют long и lat, а остальные представляют данные за весь период времени
Var1 Var2 X1 X2 Xn
1 -104.6875 47.8125 4.320732594 1.342601776 --------
2 -104.5625 47.8125 4.630135059 2.880001545 --------
3 -104.4375 47.8125 3.533707142 0.539954007 --------