Note that readGDAL can read the subdataset directly (escaping the quotes
appropriately):
d <- readGDAL("HDF5:\"RAD_NL21_PCP_NA_20070720.h5\"://image1/image_data")
It also seems that your gdal_translate to TIFF output has the extension
.img, but then you attempt to read from a tif.
Just a tho
I am looking for a method that would determine whether any part of a
line segment defined by (x0, y0) and (x1, y1) intersects a given polygon
(true or false). I don't have any other geospatial needs here, so if a
solution exists using base R that would be great.
Thanks,
Scott Waichler
Pacific Nor
If your data is weekly, then you'll want to define an unlimited time
dimension, so that your data structure looks something like this:
Radon[station_number,week_number]
And your dimensions look like something like this:
station_number[1:n_stations]
week_number[1:n_weeks] (make this unlimited wit
Dear Dave
Great, that worked! Many thanks. Now I want to put more weekly data (as RN01
was data for the first week of the year 2006) into the var_radon, but this
is not possible for dim_station, as there are only 3 dimensions. Do I have
to redefine 51 other dimensions? I tried to use list() to def
On Thu, 9 Aug 2007, Kitty Lee wrote:
> Hi. I'm trying to run spatial conditional autoregression (spautolm) in
> SPDEP. My command is simply:
>
> spautolm(vote~income, listw=listw5)
>
?spautolm
Note that the default method is method="full", you would need
method="Matrix". Only method="full" and
Paul,
I have added a CC to the R-sig-geo list, since I myself have no experiece
of using this format in practice. If you are not a member, please join at:
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
so that you will see any replies.
Did you try gdalinfo on the output file?
The error messa