2016-04-07 29 views
2

R a netCDF dosyasını raster olarak okumaya çalışıyorum. NetCDF dosyası, okyanustaki ortalama yıllık sıcaklığı, boylam, enlem ve derinliğin bir fonksiyonu olarak tanımlar. Yüzey okyanusuna (yani, netCDF dosyasındaki ilk seviye) ve t_an netCDF dosyasındaki değişkenin adıyla ilgileniyorum.R: Bir netCDF dosyasını raster olarak okuma

MyRast <- raster("Temperature.nc", level = 1, varname = "t_an") 

Bu bana aşağıdaki uyarı verir: Bu nedenle, aşağıdaki kodu kullanabilirsiniz

Warning message: 
In .getCRSfromGridMap4(atts) : cannot process these parts of the CRS: epsg_code=EPSG:4326 

Gördüğünüz gibi, netCDF dosya CRS 4326 EPSG (veya WGS 84), henüz raster vardır oluşturulduğunu aşağıdaki CRS sahiptir:

+proj=longlat +lon_0=0 +a=6378137 +rf=298.257232666016 

doğru CRS ile dosyada netCDF okumak nasıl bir fikir?

+0

Eğer NetCDF bir bağlantı verebilir? –

+0

Bu bağlantıdan indirilen https://www.nodc.noaa.gov/cgi-bin/OC5/woa13/woa13.pl format olarak netCDF ve mevcut ızgara olarak 1/4º ve dosya t00_04v2.nc olarak adlandırılır. Yıllık başlık. NetCDF kullanarak – Lyngbakr

+0

ölü kümesi var mı? – MikeJewski

cevap

3

Eğer verilerde okuduktan sonra projeksiyon kuruyor yapmanız gereken:

r <- raster("Temperature.nc", varname = "t_an") 
proj4string(r)=CRS("+init=EPSG:4326") 
+1

veya 'crs (r) <-" + init = EPSG: 4326 "' – RobertH