Hi Tim,
Yes, basically, I have and just did again, try this -- same result...
Tom
On Mon, Aug 10, 2020 at 4:34 PM Tim Salabim wrote:
> Have you tried
> colorkey=list(labels = list(labels = list("Af","Am","Aw")))
> ?
>
>
>
> On Mon, Aug 10, 20
at 4:15 AM Roger Bivand wrote:
> On Wed, 15 Jul 2020, Thomas Adams wrote:
>
> > Hi all,
> >
> > It's been some time since I approached universal kriging using gstat (I
> > struggled with this previously, years ago:
> > https://stat.ethz.ch/pipe
Hi all,
It's been some time since I approached universal kriging using gstat (I
struggled with this previously, years ago:
https://stat.ethz.ch/pipermail/r-sig-geo/2006-May/001017.html).
The problem...
Within GRASS GIS, using R, I do this...
(1) read a raster DEM into R from GRASS -- srtm <-
rmat = "PROJ", multiline = "NO") :
Discarded datum unknown in CRS definition
> image(srtm,col = terrain.colors(20))
Thank you for your help!
Best regards,
Tom
On Sun, Jul 12, 2020 at 9:19 AM Roger Bivand wrote:
> On Sun, 12 Jul 2020, Thomas Adams wrote:
>
e to read
the srtm data quickly and generate an image. I am puzzled by the "cannot
allocate memory block of size 16777216 Tb" message with my GRASS Region
settings in place.
Thank you for all your help!
Best,
Tom
On Sat, Jul 11, 2020 at 11:02 PM Thomas Adams wrote:
>
ay
> roger.biv...@nhh.no
>
> ----------
> *Fra:* Thomas Adams
> *Sendt:* lørdag 11. juli 2020, 17:59
> *Til:* Roger Bivand
> *Kopi:* r-sig-Geo@r-project.org
> *Emne:* Re: [R-sig-Geo] Issue with rgrass7 and readRAST
>
> Hi Roger!
>
> Here is what was written in
> /media/
-- 1 teaiii teaiii 0 Jul 9 20:32 801.0
-rw-r--r-- 1 teaiii teaiii 0 Jul 9 22:40 188.0
-rw-r--r-- 1 teaiii teaiii 0 Jul 9 22:51 222.0
And there is plenty of space on the disk drive...
Tom
On Sat, Jul 11, 2020 at 9:26 AM Roger Bivand wrote:
> Tom
>
> On Fri, 10 Jul 202
Hi,
I'm using R in combination with GRASS GIS 7.8.3, attempting to read a GRASS
raster file into R with the rgrass7 interface. I have done this kind of
thing countless times previously without problems. It's quite possible I am
running into memory issues, but my Ubuntu 18.04 Linux computer has 32
n Fri, 3 Jul 2020, Thomas Adams wrote:
> >
> >> Hi Roger,
> >>
> >> I just tried re-installing withou using the argument and got this:
> >
> > Very useful, thanks - painful but another case of staying too close to
> the
> > front. proj_create_c
Roger,
OK; I'll update GDAL and PROJ to 3.0.4 and 6.3.1 and get back to you... you
anticipated my question... thank you!
Best,
Tom
On Fri, Jul 3, 2020 at 11:24 AM Roger Bivand wrote:
> On Fri, 3 Jul 2020, Roger Bivand wrote:
>
> > On Fri, 3 Jul 2020, Thomas Adams wrote:
> &
6.o' failed
make: *** [proj6.o] Error 1
ERROR: compilation failed for package ‘rgdal’
* removing ‘/home/teaiii/R/x86_64-pc-linux-gnu-library/4.0/rgdal’
* restoring previous ‘/home/teaiii/R/x86_64-pc-linux-gnu-library/4.0/rgdal’
On Fri, Jul 3, 2020 at 5:45 AM Roger Bivand wrote:
> On Thu, 2 Ju
Thank you, Roger! I'll re-install without using
--configure-args=--with-proj_api="proj_api.h
Best,
Tom
On Fri, Jul 3, 2020 at 5:45 AM Roger Bivand wrote:
> On Thu, 2 Jul 2020, Thomas Adams wrote:
>
> > Roger,
> >
> > Thank you so much; after downloading th
Roger,
Thank you so much; after downloading the *.tar.gz file, the that did the
trick, using:
sudo R CMD INSTALL --configure-args=--with-proj_api="proj_api.h"
rgdal_1.5-13.tar.gz
Best regards!
Tom
On Thu, Jul 2, 2020 at 3:15 PM Roger Bivand wrote:
> On Thu, 2 Jul 2020, Thoma
>
> as there have been recent patches to rgdal that might help that have
> been incorporated there.
>
>
> On Thu, Jul 2, 2020 at 1:47 PM Thomas Adams wrote:
> >
> > Hi all,
> >
> > I have searched the internet and have spent a few days trying to get to
&g
Hi all,
I have searched the internet and have spent a few days trying to get to the
problem when I try to install rgdal into my Ubuntu 18.04 system. I have
done this many times previously, but cannot seem to track down what I am
doing wrong now. I have tried a forced install from sources with and
Hi Erin,
Are you interested in point or gridded data and does location or variable
type matter? What are you doing?
Tom
On Tue, Sep 25, 2018 at 11:03 PM Erin Hodgess
wrote:
> Hello everyone:
>
> Could someone recommend a good source for spatial temporal data from the
> "real world", please?
>
*http://sophieleblanc6.wix.com/photographie
> <http://sophieleblanc6.wix.com/photographie>*
> De : "Thomas Adams"
> A : sophie.lebl...@neuf.fr,"r-sig-Geo@r-project.org"
> Envoyé: mercredi 8 août 2018 16:07
> Objet : Re: [R-sig-Geo] bottom area of a lake calcu
Sophie,
My apologies, I was thinking in terms of using GRASS GIS... but the idea is
the same
Tom
On Wed, Aug 8, 2018 at 10:06 AM, Thomas Adams wrote:
> Sophie,
>
> I think a reasonable approximation might be obtained:
>
> (1) if you were to calculate the slope at each g
Sophie,
I think a reasonable approximation might be obtained:
(1) if you were to calculate the slope at each grid point using
r.slope.aspect
(2) use r.mapcalc to calculate A/cos(slope), where A is the pixel projected
area (e.g., 100mx100m = 1 m^2) -- lat-long coordinates will not work.
Hi Melanie,
Undoubtedly, there is an R only solution, but I find it far easier to use
GRASS GIS with the r.out.gdal command. A simple bash shell script can do
this. I just verified this with an arbitrary region over N America from a
global dataset of the Canadian CMC ensemble numerical weather
Try Googling "hydroGOF"; you'll find:
https://cran.r-project.org/web/packages/hydroGOF/vignettes/hydroGOF_Vignette.pdf
http://www.slideshare.net/hzambran/egu2012-11549go-fsforextremeevents4web
http://www.slideshare.net/hzambran/egu2010-ra-statisticalenvironmentfordoinghydrologicalanalysis-9095709
All:
I was previously able to do the kind of thing I will illustrate below on an
Intel based RedHat Linux system a few years ago. The simple task is to use
GRASS GIS with R, reading both raster and vector data into R through
rgrass7 (previously with spgrass6). Using spplot, display the raster and
adless server. Turns
> out you need X11 (or at least Cairo) for a working graphics device for
> PNG and JPG formats. Any chance that you can post your configure
> arguments or build information?
>
> Dylan
>
> On Sun, Nov 8, 2015 at 2:42 PM, Thomas Adams <tea...@gmail.com>
Mike,
This would be really helpful to me as well!
Thank you,
Tom
On Thu, Aug 20, 2015 at 4:39 PM, Michael Sumner mdsum...@gmail.com wrote:
On Friday, August 21, 2015, Roberto Horn robertomh...@gmail.com wrote:
Hello
I am looking for a package on R that can download data from Global
Tim,
WOW! This is very cool...
Using the tips from others, I did the following (to summarize):
install.packages(leaflet,dependencies=T)
library(devtools)
install_github(environmentalinformatics-marburg/Rsenal)
Works great in the few minutes I've tried it.
Thanks!!
Tom
On Fri, Jul 24, 2015
Michael,
I would combine the use of GRASS GIS and R (making use of spgrass7); first,
using the vector point locations:
(0) import raster and vector data into GRASS GIS
(1) use GRASS v.sample (http://grass.osgeo.org/grass70/manuals/v.sample.html)
to sample the raster data at the vector point
Wade,
I was able to import the data into GRASS GIS (which uses gdal) without any
problems (ultimately); the issue you're having is that the ACTUAL map
projection is either not specified -OR- is specified as 'unprojected',
which is what I think is the case. In reality, the map projection is
Nuno,
This works for me:
require(stats)
set.seed(14)
x - rchisq(100, df = 4)
hist(x)
meanx-mean(x)
meanx
[1] 3.734162
abline(v=meanx, col=red)
Cheers!
Tom
On Tue, Nov 18, 2014 at 10:56 AM, Nuno Sá nunocesard...@gmail.com wrote:
Hello!
My aim is to add a Median or a Mean line to an
=m.rst,col=blue,lwd=2)
Error in int_abline(a = a, b = b, h = h, v = v, untf = untf, ...) :
plot.new has not been called yet
Regards!!
Thank you once again!
Nuno
On 18 November 2014 18:04, Thomas Adams tea...@gmail.com wrote:
Nuno,
This works for me:
require(stats)
set.seed(14)
x
, ...) :
plot.new has not been called yet
Regards!!
Thank you once again!
Nuno
On 18 November 2014 18:04, Thomas Adams tea...@gmail.com wrote:
Nuno,
This works for me:
require(stats)
set.seed(14)
x - rchisq(100, df = 4)
hist(x)
meanx-mean(x)
meanx
[1] 3.734162
Kevin,
My experience has been that it's far better to install the gstat binary and
dependencies; get *gdal* from here: http://grassmac.wikidot.com/frameworks
, then install gstat with all dependencies.
Tom
On Thu, Oct 9, 2014 at 4:07 PM, Kevin Wolz wo...@illinois.edu wrote:
Trying to get the
Spencer,
Maybe I'm not sure what you're asking but here http://www.gadm.org/version1
are current shapefiles by country. The highlighting (selection) canals be
done in GRASS GIS, QGIS, GMT
Tom
On Tue, Jul 1, 2014 at 2:29 PM, Spencer Graves spencer.gra...@prodsyse.com
wrote:
What would
I think this is what you want:
The *aspect* output raster map indicates the direction that slopes are
facing. The aspect categories represent the number degrees of east.
Category and color table files are also generated for the aspect map layer.
The aspect categories represent the number degrees
OK, I thought that might be your issue; it is an arbitrary value because
the surface is not sloping in ANY direction -- a flat surface!
Consequently, the aspect is assigned the value 0.5*pi radians (or 90
degrees if unit =âdegreesâ)
You possibly could assign a different value, such as null or
the PRISM
data so I'm not sure why there is the difference
.
On Tuesday, January 14, 2014 5:32:48 AM, Thomas Adams
tea...@gmail.comjavascript:_e({}, 'cvml', 'tea...@gmail.com');
wrote:
Julie,
I suspect it's a combination of a mismatch between either/or the map
projection and geodetic datum
Julie,
I suspect it's a combination of a mismatch between either/or the map
projection and geodetic datum for the MODIS NDVI and PRISM data. The
ModisDownload has the option, proj_type:
Output projection short name. Valid values are AEA (Albers EqualArea), ER
(Equirectangular), GEO (Geographic),
Dominik,
Look here:
http://r-sig-geo.2731867.n2.nabble.com/GSTAT-Optimize-power-value-for-IDW-td2765918.html
Tom
On Mon, Aug 12, 2013 at 5:17 PM, Dominik Schneider
dominik.schnei...@colorado.edu wrote:
Hi,
I'd like to optimize the IDW power, similar to what ArcGIS does. I found
this old
Erin,
I basically duplicated what you did and got:
GDALinfo(/home/teaiii/ASDAR/70042108.tif)
rows1200
columns 1320
bands 1
lower left origin.x174.2
lower left origin.y-37.5
res.x 0.000833
res.y 0.000833
ysign -1
oblique.x 0
oblique.y
Edzer,
That's really pretty cool stuff; I think I can find some applications for
it!
Cheers!
Tom
On Sun, Jun 23, 2013 at 5:57 PM, Edzer Pebesma
edzer.pebe...@uni-muenster.de wrote:
I did some simple things trying to combine sp and igraph objects to
create spatial networks, and published
Andy,
I'd suggest a Google search to find what you want if your determined to use
R (and there are packages to do mostly want you want), but R may not be the
best solution. You may want to try GMT (http://gmt.soest.hawaii.edu/), or
GRASS GIS (http://grass.osgeo.org/), or Qgis
Roger,
I don't understand; since akima is freely available for download (and other
packages use akima) from CRAN, how does this licensing work? Are the
developers using akima in their packages paying a fee to the akima
developer? Just curious...
Thank you,
Tom
On Fri, Mar 29, 2013 at 4:10 PM,
Roger,
I see -- thank you for the explanation!
Best,
Tom
On Sat, Mar 30, 2013 at 1:42 PM, Roger Bivand roger.biv...@nhh.no wrote:
On Sat, 30 Mar 2013, Thomas Adams wrote:
Roger,
I don't understand; since akima is freely available for download (and
other
packages use akima) from CRAN
Marco,
I think the use of R with GRASS is clearly the was to go Michael Barton has
done great work with R and GRASS for analyzing archaeological data. Go to
http://grass.osgeo.org ; get the Open Source GIS: A GRASS GIS Approach 3rd
edition. The advantage of using GRASS is from the point of view
Jin,
I'd be willing to wager it can be done in R alone as someone is sure to
point out, but since you're using R (and so, willing to go the open source
route), you may be willing to take another step and use GRASS GIS (
http://grass.osgeo.org) the wonderful thing about this is that your
problem
Aurelie,
A simple Google search for arctic oceanographic data lead me to:
http://nsidc.org/data/docs/noaa/g01938_ewg_arctic_met_atlas/index.html
http://iabp.apl.washington.edu/
http://www.whoi.edu/science/PO/arcticgroup/projects/andrey_project2/indexAP.html
Wade,
Look at the contributed package gstat; there is also the standalone version
at http://www.gstat.org/ which interfaces nicely with GRASS GIS (
http://grass.osgeo.org/). I have used the latter approach
Tom
On Fri, May 11, 2012 at 12:15 PM, Wall, Wade A ERDC-RDE-CERL-IL
Tony,
I certainly don't want to belabor the point, but GRASS GIS with R is fully
scriptable (shell, Perl, whatever); I have done this with many projects. GRASS
shell scripts are far easier to write in my experience than ArcGIS. Several of
my scripts are completely automated and run off of cron
Marco,
I just took a look at the data; unfortunately, it's not in a very usable format
without reformatting it for
importing the data into either R or a GIS, such as GRASS GIS. I'm sure an
experienced R programmer could
write a clever R script to do what you want. Personally, I would write a
Barry,
WRT your last line — I could not agree more; the same thought came to me
as well!
Tom
On 1/16/11 11:04 AM, Barry Rowlingson wrote:
On Sun, Jan 16, 2011 at 3:28 PM, Thomas Adamsthomas.ad...@noaa.gov wrote:
Marco,
I just took a look at the data; unfortunately, it's not in a very
49 matches
Mail list logo
