search for: tomroche

Displaying 6 results from an estimated 6 matches for "tomroche".

Did you mean: tom_roche
2012 Aug 28
4
[ncdf4] error converting GEIA data to netCDF
...an successfully ncvar_put(...) data to a file, but when I try to ncvar_get(...) the same data, I get > Error in if (nc$var[[li]]$hasAddOffset) addOffset = nc$var[[li]]$addOffset else addOffset = 0 : > argument is of length zero How to fix or debug? details: R code @ https://github.com/TomRoche/GEIA_to_NetCDF successfully (if crudely) uses R packages={ncdf4, maps, fields} to + extract data from a GEIA-distributed datafile https://github.com/downloads/TomRoche/GEIA_to_netCDF/N2OOC90Y.1A + display the data (mostly successfully--the map's legend has problems which I'll attack l...
2012 Oct 21
1
[newbie] failure to plot a RasterLayer with raster::plot or fields::image.plot
summary: spatial data to be input to a regional-scale environmental model must (1) be converted to netCDF and then (2) "regridded" (cropped, projected, increased resolution). In a public git repository https://github.com/TomRoche/GEIA_to_NetCDF I have R code (with bash drivers) that does the conversion step, and plots the converted output, apparently correctly. However attempts to plot the output of the regridding step (twice with raster::plot, once with fields::image.plot) are very wrong. How to fix? details: I need to...
2012 May 01
1
[fields:image.plot] subtitle under title (not image)?
summary: how to make image.plot print a subtitle between the title and the image, rather than under the image? details: I've got a project https://github.com/TomRoche/ioapi-hack-R that illustrates the use of (et al) the R packages {ncdf4, fields, M3} for processing and visualizing IOAPI data. The data being visualized consists of a series of layers (mostly representing emissions from a particular tuple of agricultural {cultivar, cultivation technique}) over a s...
2012 Sep 11
2
R equivalent of python module structure and functionality?
...dline) > source('./foo.r) > foo(bar='baz') ? I'm looking for the 'R equivalent' of how python supports this usecase. details: As discussed in the thread beginning https://stat.ethz.ch/pipermail/r-help/2012-September/323255.html I have a script https://github.com/TomRoche/GEIA_to_netCDF/blob/master/netCDF.stats.to.stdout.r that takes named arguments without undue pain. I would also like to be able to call it as a function from other scripts. How to do that in R? In case that's not specific enough :-) I know how to structure files/ modules in python like http:/...
2012 Sep 06
1
[Rscript] difficulty passing named arguments from commandline
Wanting a commandline solution (for a problem detailed @ http://mailman.unidata.ucar.edu/mailing_lists/archives/netcdfgroup/2012/msg00279.html ) I turned to Rscript, and whacked out the q'n'd https://github.com/TomRoche/GEIA_to_netCDF/blob/master/netCDF.stats.to.stdout.r However it wasn't as quick as hoped, because I spent quite a bit of time figuring out how to pass the arguments. This works (note the quoting): $ Rscript ./netCDF.stats.to.stdout.r 'netcdf.fp="./GEIA_N2O_oceanic.nc"' 'v...
2013 Apr 26
1
[newbie] how to find and combine geographic maps with particular features?
...;, # see `?M3::get.map.lines.M3.proj` > units='m', # or 'km': see `?M3::get.map.lines.M3.proj` > extents.fp, # path to extents file > extents, # raster::extent object > LCC.parallels=c(33,45), # LCC standard parallels: see https://github.com/TomRoche/cornbeltN2O/wiki/AQMEII-North-American-domain#wiki-EPA > CRS # see `sp::CRS` > ) { > library(M3) > ## Will replace raw LCC map's coordinates with: > metadata.coords.IOAPI.list <- M3::get.grid.info.M3(extents.fp) > metadata.coords.IOAPI.x.orig <-...