[R-sig-Geo] Automate the extraction of climate variable at different time depths from netcdf in r

Barnabas Daru darunabas at gmail.com
Sat Sep 6 19:10:24 CEST 2014


Dear all,
I write to seek help with extracting climate data at various times and
depth from NETCDF in R.
The climate data is sea surface temperature monthly means from 1854 to
2014.
My goal is to automate the extraction of SST climate data for a set of
points with different time depths.

I have successfully used the 'brick' function in the raster package to get
climate data at various time depths as follows:

b <- brick("~sst.mnmean.nc", varname="sst")

mydate <- b[["X1869.09.01"]] # SST for September in 1869

mydata <- read.csv("~Species one.csv")

mydataSPDF <- SpatialPointsDataFrame(mydata[,5:6],data.frame(mydata))

extract.mydata <- extract(mydate, mydataSPDF, sp=TRUE)

extract.mydata <- data.frame(extract.mydata)

write.csv(extract.mydata, file = "Species_one_extracted.csv")

My major challenge is that I have to keep extracting manually for each time
depth, then copy and paste in the main dataframe and the data has about
3000 observations, meaning I have to keep copying and pasting ~3000 times
for one species!
Is there a way I can iterate the process in R to loop through each month's
climate and extract the data for each point and time depth rather than to
do it manually?

Here's how my data frame looks like:
  Date year month day lon lat SST  1855-01-01 1855 1 1 11.6861 57.7254
3.440000057  1859-07-26 1859 7 26 25.46122 60.26766 13.25  1860-01-01 1860 1
1 10.9154 53.9484 15.10999966  1861-07-10 1861 7 10 7.79588 58.08673
15.71000004  1861-07-26 1861 7 26 -2.84072 54.0778  1861-08-17 1861 8 17
11.9792 57.51298  1862-01-01 1862 1 1 22.20467 60.27955  1862-08-05 1862 8 5
11.78316 57.61649  1862-08-23 1862 8 23 11.9792 57.51298  1863-08-26 1863 8
26 10.72237 59.97258  1863-08-28 1863 8 28 15.53721 56.20091  1863-09-22
1863 9 22 16.37849 56.84255  1864-07-05 1864 7 5 -4.07097 51.09542
1864-07-18 1864 7 18 -4.21368 51.0928  1864-08-26 1864 8 26 -4.07097
51.09542  1865-09-07 1865 9 7 10.76144 59.91271  1865-09-27 1865 9 27
10.53107 60.43395  1867-08-28 1867 8 28 12.82669 55.86822  1868-01-01 1868 1
1 12.8944 56.6897  1868-01-01 1868 1 1 4.82685 53.13793
Thanks and kind regards
Barnabas Daru

-- 

  \-/
   /\
  /--|
 /---/ *Daru, Barnabas Haruna*
 |--/  PhD Candidate,
 \-/   African Centre for DNA Barcoding,
 /\    University of Johannesburg,
/--\   PO BOX 524 Auckland Park, 2006,
|---\  South Africa
 \---\ Lab: +27 11 559 3477
  \--| Mobile: +27 7381 89 583
   \-/
   /\  My homepage: http://barnabasdaru.com
<http://acdb.co.za/index.php?page=mr-barnabas-h-daru>
  /--\

	[[alternative HTML version deleted]]



More information about the R-sig-Geo mailing list