Displaying 6 results from an estimated 6 matches for "tomroch".
Did you mean:
tomroche
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...
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...
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"' '...
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 <-...