https://stat.ethz.ch/pipermail/r-sig-geo/2013-April/017947.html
Can one
+ create a 4D RasterBrick or RasterStack
+ instruct projectRaster
so as to get 4D regridded output? Or must one
- create 3D RasterLayer's (i.e., one per LAY) from the 4D Brick
- regrid each 3D Layer
- reassemble the
hour3ly.in.raster - raster::brick(
hour3ly_in_fp, varname=hour3ly_in_datavar_name, level=2, lvar=3)
hour3ly.in.raster - raster::brick(
hour3ly_in_fp, varname=hour3ly_in_datavar_name, level=3, lvar=2)
and got the same results. What am I doing wrong?
TIA, Tom Roche tom_ro...@pobox.com roche@epa.gov
Different?
Apologies if this is a FAQ, but I'm not seeing the answer in doc
(sp.pdf or intro_sp.pdf) or via googling.
TIA, Tom Roche tom_ro...@pobox.com
___
R-sig-Geo mailing list
R-sig-Geo@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
(1)
...
r.14 - raster(nrow=lat.size, ncol=lon.size, crs=my.crs)
values(r.14) - foo.timestep(14)
brick( r.1, r.2, r.3, r.4, r.5, r.6, r.7, r.8, r.9, r.10,
r.11, r.12, r.13, r.14)
Or is there another/better way to do this?
TIA, Tom Roche tom_ro...@pobox.com
by
latticeExtra::layer(...)
TIA, Tom Roche tom_ro...@pobox.com
___
R-sig-Geo mailing list
R-sig-Geo@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
https://stat.ethz.ch/pipermail/r-sig-geo/2013-February/017531.html
How to convert from a simple (x,y) matrix (as produced by
M3::get.map.lines.M3.proj) to a SpatialLines, or to coordinates
consumable by latticeExtra::layer?
One way to do this, though ugly, is the current answer @
. But I'd much prefer to make this
work like the other example, which resembles other code I have that is
using `rasterVis::layerplot`. Unfortunately that other code is using
unprojected/lon-lat data, and I must also handle projected (probably all
LCC) data.
TIA, Tom Roche tom_ro...@pobox.com
raster.pdf, and the source for raster:::.rasterObjectFromCDF seems quite
inscrutable (YMMV :-)
TIA, Tom Roche tom_ro...@pobox.com
___
R-sig-Geo mailing list
R-sig-Geo@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
Tom Roche Sun, 10 Feb 2013 14:40:05 -0500
$ ls -alh ./METCRO3D_080101
-r--r--r-- 1 me mygroup 5.3G Jan 8 11:32 ./METCRO3D_080101
# extension used below
$ ln -s ./METCRO3D_080101 ./METCRO3D_080101.nc
...
raster('./METCRO3D_080101.nc', varname='ZF')
...
Error in .local(.Object
=4032group_id=1014func=browse
but I'm not seeing equivalent UI for package=rgdal.
TIA, Tom Roche tom_ro...@pobox.com
___
R-sig-Geo mailing list
R-sig-Geo@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
the datavars in one netCDF file. Can one write from multiple Raster*
to multiple datavars in a single netCDF file? An optimistic reading
of raster.pdf suggests this is supported via the sequence
[writeStart, writeValues, writeStop], but I'd appreciate
confirmation.
TIA, Tom Roche
-term, It Would Be Nice if someone provided a proper (i.e.,
taking a user-provided CRS argument) 3D geographical-to-Cartesian
conversion.
Your assistance is appreciated (even if only useful in future),
Tom Roche tom_ro...@pobox.com
___
R-sig-Geo mailing
Tom Roche Sun, 30 Dec 2012 00:40:55 -0500
Is there packaged code to convert geographical coordinates (e.g.,
longitude, latitude, elevation) to [simple] Cartesian coordinates
in 3-space?
Edzer Pebesma Sun, 30 Dec 2012 11:11:19 +0100
rgdal::spTransform?
I'm aware of that, but don't how
-intelligence test?
TIA, Tom Roche tom_ro...@pobox.comsample code follows to EOP
# Convert geographic coordinates (lon, lat, elevation) to
# Cartesian coordinates (x, y, z) relative to earth center.
geo2cart - function(lon.lat.elev) {
library(pracma) # for sph2cart
stopifnot(is.numeric
.
TIA, Tom Roche tom_ro...@pobox.com
___
R-sig-Geo mailing list
R-sig-Geo@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
'))
...
Date: 2012-08-10
Author: Oscar Perpiñán Lamigueiro
Perhaps you could fix that?
TIA, Tom Roche tom_ro...@pobox.com
___
R-sig-Geo mailing list
R-sig-Geo@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
Library error:
Use linux.
HTH, Tom Roche tom_ro...@pobox.com
___
R-sig-Geo mailing list
R-sig-Geo@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
] to longitudes=[-180, +180]?
TIA, Tom Roche tom_ro...@pobox.com
___
R-sig-Geo mailing list
R-sig-Geo@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
Tom Roche Wed, 19 Dec 2012 13:52:38 -0500
How to best/easiest range-shift lon-lat data?
Matt Landis Wed, 19 Dec 2012 14:09:57 -0500
[raster::shift] will work, but [raster::rotate] is even easier
Thanks! Now this code-
global.proj - CRS('+proj
either do something simple to wrld_simpl's longitude
values, or Something Completely Different that would Just Work Better,
but don't know. Your assistance is appreciated.
TIA, Tom Roche tom_ro...@pobox.com
___
R-sig-Geo mailing list
R-sig-Geo@r
on that separately
par.strip.text=list(cex=0.5)
)
)
# end example
TIA, Tom Roche tom_ro...@pobox.com
___
R-sig-Geo mailing list
R-sig-Geo@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
while keeping one ear on the US
elections (at least that worked out better :-)
Onward to a more complex test ...
thanks, Tom Roche tom_ro...@pobox.com
___
R-sig-Geo mailing list
R-sig-Geo@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
of the subglobal grid.
? If so, which CRS to use?
TIA, Tom Roche tom_ro...@pobox.com
___
R-sig-Geo mailing list
R-sig-Geo@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
week.
FWIW, Tom Roche tom_ro...@pobox.com
___
R-sig-Geo mailing list
R-sig-Geo@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
output box, and furthermore that the output box values
should approximate the input box values. Am I missing something? If not:
How to call GSIF API to compute the values of the output boxes from the
input grid and its values?
Your assistance is appreciated, Tom Roche tom_ro...@pobox.com
Tom Roche Thu, Nov 1, 2012 at 2:21 AM
See code [for lon-lat to LCC regridding] @
https://github.com/TomRoche/GEIA_to_NetCDF
esp
https://github.com/TomRoche/GEIA_to_netCDF/blob/master/regrid.global.to.AQMEII.r
ping yang Tue, 6 Nov 2012 20:36:12 -0500 (rearranged)
Thank you so much for your
package=M3.
But to get res3D from grid3D, for my usecase, should I use
* gstat::krige? If so, what model? (or no model?)
* gstat::idw? If so, what formula?
* or something else?
Your assistance is appreciated! Tom Roche tom_ro...@pobox.com
___
R-sig-Geo
that, or is it (like, IIUC, raster::projectRaster)
restricted to 2D?
Your assistance is appreciated! Tom Roche tom_ro...@pobox.com
___
R-sig-Geo mailing list
R-sig-Geo@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
/plotLayersForTimestep.r
All code (including drivers providing a test fixture), data, downloads,
and README are available @
https://github.com/TomRoche/GEIA_to_netCDF
HTH, Tom Roche tom_ro...@pobox.com
___
R-sig-Geo mailing list
R-sig-Geo@r-project.org
https
be quite similar. See code @
https://github.com/TomRoche/GEIA_to_NetCDF
esp
https://github.com/TomRoche/GEIA_to_netCDF/blob/master/regrid.global.to.AQMEII.r
HTH, Tom Roche tom_ro...@pobox.com
___
R-sig-Geo mailing list
R-sig-Geo@r-project.org
https
incorrect (global) extents,
as described @
https://stat.ethz.ch/pipermail/r-sig-geo/2012-October/016504.html
FWIW, Tom Roche tom_ro...@pobox.com
___
R-sig-Geo mailing list
R-sig-Geo@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
https://stat.ethz.ch/pipermail/r-sig-geo/2012-October/016503.html
- Evidence for a bug in raster::plot :
doh! I should have said
+ Evidence for a bug in raster::projectRaster :
rest of previous post
https://stat.ethz.ch/pipermail/r-sig-geo/2012-October/016503.html
still holds, and follows to
/cornbeltN2O/wiki/images/Cmaq_aqmeii_domain.png
the boundaries of which are the extents of the output domain.
How to fix?
Your assistance is appreciated, Tom Roche tom_ro...@pobox.com
___
R-sig-Geo mailing list
R-sig-Geo@r-project.org
https
with raster::plot? Or am I missing something?
TIA, Tom Roche tom_ro...@pobox.com
___
R-sig-Geo mailing list
R-sig-Geo@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
of extents to display that
data with raster::plot?
FWIW, Tom Roche tom_ro...@pobox.com
___
R-sig-Geo mailing list
R-sig-Geo@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
, or just the
relation between an added map and the regridded data).
Your assistance is appreciated, Tom Roche tom_ro...@pobox.com
___
R-sig-Geo mailing list
R-sig-Geo@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
: 'debian' is not one thing, it has versions.
Are you running squeeze? Try installing the GEOS packages from
wheezy/testing: see version availability @
http://packages.debian.org/search?keywords=geossearchon=namessuite=allsection=all
HTH, Tom Roche tom_ro...@pobox.com
values in the output netCDF, while there were
none in either the input netCDF or the input raster?
details:
Tom Roche Wed, 03 Oct 2012 11:34:58 -0400
I'm getting very wrong output
Robert J. Hijmans Wed, 3 Oct 2012 20:27:07 -0700
I do not know why you say the output is wrong. Why you would
Tom Roche Wed, 03 Oct 2012 11:34:58 -0400
I'm getting very wrong output
Robert J. Hijmans Wed, 3 Oct 2012 20:27:07 -0700
I do not know why you say the output is wrong. Why you would expect
that input variables would have to be preserved
I should have worded that better. Moreover, since
='to'); how to do that?
your assistance is appreciated, Tom Roche tom_ro...@pobox.com
___
R-sig-Geo mailing list
R-sig-Geo@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
is the numerator of the resolution, with the
denominator being the counts above), and the domain is not rectangular
in lon-lat (see the image just above the previous link). Am I
misunderstanding how to define 'res'?
And what I should do to debug the hang?
Your assistance is appreciated, Tom Roche tom_ro
that does not have
these problems? E.g., HTML doc should solve the truncation problem.
TIA, Tom Roche tom_ro...@pobox.com
___
R-sig-Geo mailing list
R-sig-Geo@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
(esp R-based
ones) in this domain.
Feel free to reply on- or off-list, Tom Roche tom_ro...@pobox.com
___
R-sig-Geo mailing list
R-sig-Geo@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
haven't used it, but have heard good things about Panoply
http://www.giss.nasa.gov/tools/panoply/
which has a binary for windows.
HTH, Tom Roche tom_ro...@pobox.com
___
R-sig-Geo mailing list
R-sig-Geo@r-project.org
https://stat.ethz.ch/mailman/listinfo
are appreciated, as would be pointers to
helpful resources. Thanks in advance, and feel free to forward,
Tom Roche tom_ro...@pobox.com
___
R-sig-Geo mailing list
R-sig-Geo@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
AQMEII domain) that are 1 GB.
HTH, Tom Roche tom_ro...@pobox.com
___
R-sig-Geo mailing list
R-sig-Geo@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
46 matches
Mail list logo