Re: [R-sig-Geo] Match Coordinates to NUTS 2 ID

E, col = "red", pch = 19) > ## Use the over() function to get the regions for the points > ## No regions for 2 points in the ocean > over(point_data, map_nuts2, by.id = TRUE) NUTS_ID STAT_LEVL_ SHAPE_Leng SHAPE_Area 1CZ06 2 7.018153 1.688086 2CZ06

Re: [R-sig-Geo] problem with predict() in package raster and factor variables

orest model and 'subbrick': > str(v) > str(subbrick) Please also show all the relevant R code to obtain what you want in case the error message is not related to difference in the creation of the subset of 'v' and 'subbrick' Yours sincerely / Med venlig

[R-sig-Geo] FW: proj4string error on data.frame generated from NCEP data

Reposting again again without data. So fetch the data yourself :-( Yours sincerely / Med venlig hilsen Frede Aakmann Tøgersen Specialist, M.Sc., Ph.D. Plant Performance & Modeling Technology & Service Solutions T +45 9730 5135 M +45 2547 6050 fr...@vestas.com http://www.vestas.com Com

Re: [R-sig-Geo] How to calculate decadal average and percent change in raster stacks?

format(s@z$time, format = "%Y") %in% decades[[1]] tmp <- subset(s, which(ndx)) tmp@z$time Yours sincerely / Med venlig hilsen Frede Aakmann Tøgersen Specialist, M.Sc., Ph.D. Plant Performance & Modeling Technology & Service Solutions T +45 9730 5135 M +45 2547 6050 fr...@vestas.c

Re: [R-sig-Geo] How to get the borders of scattered spacial points

Hi Jay So you have forgot all the information you got when you subscribed to the list? To refresh your memory please see: https://stat.ethz.ch/mailman/listinfo/r-sig-geo Yours sincerely / Med venlig hilsen Frede Aakmann Tøgersen Specialist, M.Sc., Ph.D. Plant Performance & Mode

Re: [R-sig-Geo] is this a clash between wrong version of rgdal and sp?

Hi Thierry Yes, our sys-admin will do that when we get a new cluster later this year. I don’t want to burden our very busy admin with support on a cluster that is going to be shut down in a few months. Yours sincerely / Med venlig hilsen Frede Aakmann Tøgersen Specialist, M.Sc., Ph.D. Plant

Re: [R-sig-Geo] is this a clash between wrong version of rgdal and sp?

gilio", "Gomez-Rubio", role = "ctb")) Depends: R (>= 2.14.0) Suggests: lattice, RColorBrewer, rgdal (>= 0.8-7), rgeos (>= 0.1-8) Imports: methods, graphics, utils, lattice, grid The sp versio is 1.0-9 and it suggests rgdal (>= 0.8-7). Installing thes

[R-sig-Geo] is this a clash between wrong version of rgdal and sp?

I have tried several combinations but nothing has helped me so far. Any help is appreciated. Yours sincerely / Med venlig hilsen Frede Aakmann Tøgersen Specialist, M.Sc., Ph.D. Plant Performance & Modeling Technology & Service Solutions T +45 9730 5135 M +45 2547 6050 fr...@vestas.com http://www.v

Re: [R-sig-Geo] plotKML 'kml' function with raster input

aster(extent(-135.2,-135.1,58.1,58.2),nrows=10,ncol=10,crs=CRS('+init=epsg:4326')) rast[] <- rnorm(ncell(rast)) kml(rast, kml.folder, kml.file, colour = 'black') Yours sincerely / Med venlig hilsen Frede Aakmann Tøgersen Specialist, M.Sc., Ph.D. Plant Performance &am

Re: [R-sig-Geo] making wrld_simpl longitude borders something other than (-180, 180)

from the interval (-180, 180) to (-100, 100)? Yours sincerely / Med venlig hilsen Frede Aakmann Tøgersen Specialist, M.Sc., Ph.D. Plant Performance & Modeling Technology & Service Solutions T +45 9730 5135 M +45 2547 6050 fr...@vestas.com http://www.vestas.com Company reg. name: Vestas Wind S

Re: [R-sig-Geo] reading postgis tables in R

proj4string=p4ss), shape.df[, -ncol(shape.df)]) SGeomdf} Please don't post in html but in plain text. Yours sincerely / Med venlig hilsen Frede Aakmann Tøgersen Specialist, M.Sc., Ph.D. Plant Performance & Modeling

Re: [R-sig-Geo] Impute missing value using k nearest neighbour

LL) Yours sincerely / Med venlig hilsen Frede Aakmann Tøgersen Specialist, M.Sc., Ph.D. Plant Performance & Modeling Technology & Service Solutions T +45 9730 5135 M +45 2547 6050 fr...@vestas.com<mailto:fr...@vestas.com> http://www.vestas.com<http://www.vestas.com/> Company reg.

Re: [R-sig-Geo] Global Precipitation Measurement data download

and SOAP and have API for jave and C++) http://pmm-dev.pps.eosdis.nasa.gov/data-access/data-sources#midador ( http://pmm.nasa.gov/data-access/data-sources#TOVAS Hope this helps. Yours sincerely / Med venlig hilsen Frede Aakmann Tøgersen Specialist, M.Sc., Ph.D. Plant Performance & M

Re: [R-sig-Geo] Need help in creating a grid

Yes, some one on this list may be able to help you but not before you show the list what you have been doing so far and where you get stucked. Yours sincerely / Med venlig hilsen Frede Aakmann Tøgersen Specialist, M.Sc., Ph.D. Plant Performance & Modeling Technology & Service Solutio

Re: [R-sig-Geo] plotKML - set color according to category

rames with different definition of the factor 'phase'. That would certainly explain your problem with different number of categories. Nobody can really help you before you provide an R coded small reproducible example that shows us the problem that you have (see e.g. http://www.r-blo

Re: [R-sig-Geo] Pre-GDAL 2: rgdal changes - from all according to their means

Okay, it seems to be an issue on our side. I'll get some help from our sysadmin at least on the linux cluster. Br. Frede Sendt fra Samsung mobil Oprindelig meddelelse Fra: Roger Bivand Dato:03/06/2015 15.10 (GMT+01:00) Til: chris english Cc: Frede Aakmann Tøgersen ,Mi

Re: [R-sig-Geo] Pre-GDAL 2: rgdal changes - from all according to their means

nroot/rgdal/pkg rgdal svn: Can't connect to host 'scm.r-forge.r-project.org': Connection refused [frtog@dkrdsfshn2 ~]$ So no anonymous access? Yours sincerely / Med venlig hilsen Frede Aakmann Tøgersen Specialist, M.Sc., Ph.D. Plant Performance & Modeling Technology & Serv

Re: [R-sig-Geo] adding a scale

rds Frede Aakmann Tøgersen Specialist, M.Sc., Ph.D. Plant Performance & Modeling Technology & Service Solutions T +45 9730 5135 M +45 2547 6050 fr...@vestas.com http://www.vestas.com Company reg. name: Vestas Wind Systems A/S This e-mail is subject to our e-mail disclaimer statemen

Re: [R-sig-Geo] plotKML throws status 37 warning

known application whereas on my new laptop it said Google Earth. Hope that helps. Frede Aakmann Tøgersen Specialist, M.Sc., Ph.D. Plant Performance & Modeling Technology & Service Solutions T +45 9730 5135 M +45 2547 6050 fr...@vestas.com http://www.vestas.com Company reg. name: Vestas Win

Re: [R-sig-Geo] inconsistent as.data.frame(SpatialPointsDF)

obj <- length(x) d <- coordinates(x) if (.hasSlot(x, "data")) { d <- cbind(d, x@data) } colnames(d)[1:2] <- c("x", "y") rownames(d) <- NULL as.data.frame(d, row.names = row.names, optional =

Re: [R-sig-Geo] inconsistent as.data.frame(SpatialPointsDF)

ng an error edit the script accordingly to change the offending part of your script. Yours sincerely / Med venlig hilsen Frede Aakmann Tøgersen Specialist, M.Sc., Ph.D. Plant Performance & Modeling Technology & Service Solutions T +45 9730 5135 M +45 2547 6050 fr...@vestas.com http://w

Re: [R-sig-Geo] Error importing ESRI Shapefile

t; dsn <- "/path/to/layer" > file.choose(dsn) Can you see any shapefiles? > ogrListLayers(dsn) Can you see any layer names? > ogrInfo(dsn=dsn, layer="layer") Any information from this call? Yours sincerely / Med venlig hilsen Frede Aakmann Tøgersen Specialist, M.Sc.,

Re: [R-sig-Geo] Polygons VS MultiLineString

quot;external/test.grd", package="raster") r <- raster(f) x <- rasterToContour(r) class(x) str(x) ## the levels levels(x@data$level) plot(r) plot(subset(x, level == 400), add=TRUE) plot(subset(x, level == 200), add=TRUE) Yours sincerely / Med venlig hilsen Frede Aakmann Tøger

Re: [R-sig-Geo] Same code line for graph with different output messages

$Partic_per is parsed to NULL and this gives the error message you see from R when trying to do the indexing lnd_sport[lnd_sport$Partic_per > 25, ]. Yours sincerely / Med venlig hilsen Frede Aakmann Tøgersen Specialist, M.Sc., Ph.D. Plant Performance & Modeling Technology & Service Solu

Re: [R-sig-Geo] How to apply this function on each raster stack layer?

ate 10 new layers ('x*1, x*2, ...')? This can be achieved by using argument 'forceapply=TRUE' 'calc(r, function(x) x * 1:10), forceapply=TRUE' Yours sincerely / Med venlig hilsen Frede Aakmann Tøgersen Specialist, M.Sc., Ph.D. Plant Performance & Modeling

Re: [R-sig-Geo] input shapefile to build graph-using R

ail lists it would be seen courteously of you to provide you affiliation so people can see who they are dealing with. Yours sincerely / Med venlig hilsen Frede Aakmann Tøgersen Specialist, M.Sc., Ph.D. Plant Performance & Modeling Technology & Service Solutions T +45 9730 5135 M +45

Re: [R-sig-Geo] Error in reading shape file with readOGR()

setup (I'm guessing you are on a unix-alike system based on the path to gdal installation). Cannot remember if a relative path works. Yours sincerely / Med venlig hilsen Frede Aakmann Tøgersen Specialist, M.Sc., Ph.D. Plant Performance & Modeling Technology & Service Solutions T +45 97

Re: [R-sig-Geo] GSTAT Co-Kriging: How can i fix a non positive definte correlation matrix

im = c(115L, 2L), .Dimnames = list( NULL, c("x", "y"))) , bbox = structure(c(3.2, 6.4, 25.6, 51.2), .Dim = c(2L, 2L), .Dimnames = list( c("x", "y"), c("min", "max"))) , proj4string = new("CRS" , projargs = NA_char

Re: [R-sig-Geo] Inverse distance weighting in 3 dimension using gstat package

mastCFD to that of z values in cfd solves the problem. See the nice plot at https://www.dropbox.com/s/q30m3edjyhbuugj/Ainterp_more_correct.png?dl=0. Again, thanks. Yours sincerely / Med venlig hilsen Frede Aakmann Tøgersen Specialist, M.Sc., Ph.D. Plant Performance & Modeling Technology &

[R-sig-Geo] Inverse distance weighting in 3 dimension using gstat package

red",cex=0.4*1.5) }, col.regions = terrain.colors, ylab = "Northings (m)", xlab = "Eastings (m)") dev.off() Yours sincerely / Med venlig hilsen Frede Aakmann Tøgersen Specialist, M.Sc., Ph.D. Plant Performance & Modeling Technology & Service S

Re: [R-sig-Geo] Export subsets of Spatial Points Data Frame to images

B$Fecha, format = "%d/%m/%Y")) SB <- SB %>% group_by(Fecha, CodPar, x, y) %>% summarise(meanSube = mean(Sube), meanEvade = mean(Evade), meanBaja = mean(Baja)) SB <

Re: [R-sig-Geo] Export subsets of Spatial Points Data Frame to images

", sep="")) # i is the name of each element in df2 print(bubble(df2[[i]], "F_Up", col=c("#00ff0088", "#00ff0088"), main = "df2" )) # subsetting using i dev.off() } Yours sincerely / Med venlig hilsen Frede Aakmann Tøgersen Specialist, M.Sc.

Re: [R-sig-Geo] Export subsets of Spatial Points Data Frame to images

instead of group_by(SB, Fecha, CodPar, x, y)? Yours sincerely / Med venlig hilsen Frede Aakmann Tøgersen Specialist, M.Sc., Ph.D. Plant Performance & Modeling Technology & Service Solutions T +45 9730 5135 M +45 2547 6050 fr...@vestas.com http://www.vestas.com Company reg. name: Vestas Wind

Re: [R-sig-Geo] overlay polygon with raster

66989/extract-data-from-raster-with-small-polygons-rounded-weights-too-small See that page for a solution to some rounding errors in the case of the polygon being very small compared to raster cell size. Yours sincerely / Med venlig hilsen Frede Aakmann Tøgersen Specialist, M.Sc., Ph.D. Plant

Re: [R-sig-Geo] rainfall interpolation using ANN

/publication/225544039_A_new_approach_for_deriving_temperature_and_salinity_fields_in_the_Indian_Ocean_using_artificial_neural_networks Hope that helps. Yours sincerely / Med venlig hilsen Frede Aakmann Tøgersen Specialist, M.Sc., Ph.D. Plant Performance & Modeling Technology & Service Sol

Re: [R-sig-Geo] How to delete some polygons from 'SpatialPolygonsDataFrame'

island then do nsw_lga_wo_island <- nsw_lga[ -153, ] Something like nsw_lga[ -153, ] <- NULL doesn't seem to work. Don't know if something similar will work. Yours sincerely / Med venlig hilsen Frede Aakmann Tøgersen Specialist, M.Sc., Ph.D. Plant Performance & Modeling

Re: [R-sig-Geo] shifting points in plot window

Hi For windows() there is a rescale argument which you can set to "fit". I don't see the same argument for X11() when on unix alike OS or for any Cairo devices on iOS. Yours sincerely / Med venlig hilsen Frede Aakmann Tøgersen Specialist, M.Sc., Ph.D. Plant Perform

Re: [R-sig-Geo] plotKML kml_layer.Raster problem

quot;SNDMHT_E" "SLTMHT_E" "CLYMHT_E" > # Plot using aesthetics kml(eberg, colour = SNDMHT_A, size = CLYMHT_A, alpha = 0.75, file = "eberg-1.kml") You can see that colour is set to the SNDMHT_A layer without apostrophes. For your raster we have > nam

Re: [R-sig-Geo] plotKML kml_layer.Raster problem

ject to plot. And ?kml_aes does not really help: aesthetics package:plotKML R Documentation Plotting aesthetics parameters Description: Parses various object parameters / columns to KML aesthetics: size of the icons, fill color, labels, altitude, width, ...

Re: [R-sig-Geo] plotKML kml_layer.Raster problem

1467") data(SAGA_pal) library(raster) r <- raster(eberg_grid["TWISRT6"]) ## Put some NA's into data r@data@values[1:10] <- NA kml(r, colour_scale = SAGA_pal[[1]], colour = TWISRT6) Yours sincerely / Med venlig hilsen Frede Aakmann Tøgersen Specialist, M.Sc., Ph.D. P

Re: [R-sig-Geo] Cause of error in spdep::errorsarlm

' And in ?lm we see: formula: an object of class '"formula"' (or one that can be coerced to that class): a symbolic description of the model to be fitted. The details of model specification are given under 'Details'. So strict

Re: [R-sig-Geo] Cause of error in spdep::errorsarlm

Hi Dan Remove the " from your formula. It should be my_formula<- resp ~ pred1 + pred2 Yours sincerely / Med venlig hilsen Frede Aakmann Tøgersen Specialist, M.Sc., Ph.D. Plant Performance & Modeling Technology & Service Solutions T +45 9730 5135 M +45 2547 6050 fr...

Re: [R-sig-Geo] Calculate the length of hail paths

mepage uses 6371 km as the radius of the earth whereas geosphere uses 6378137 m. There could be other subtleties having influence on the decimal precision. Yours sincerely / Med venlig hilsen Frede Aakmann Tøgersen Specialist, M.Sc., Ph.D. Plant Performance & Modeling Technology & Ser

Re: [R-sig-Geo] Space-time movement Trajectory on a hovmoller plot

uot;),as.Date("1990-01-01"),by="6 month")), Latitude = -10:10) hovmoller(SST, panel=panel.levelplot.raster, xscale.components=xscale.raster.subticks, interpolate=TRUE, par.settings=RdBuTheme) + xyplot(Time ~ Latitude, data = so

Re: [R-sig-Geo] Overlaying time designated spatial points from a data frame to match the time of a raster image in a rasterStack

h(d)){ cells <- sample(1:ncell(rr[[1]]), 10) xy <- xyFromCell(rr, cells) tmp <- data.frame(x = xy[,"x"], y = xy[, "y"], Month = format(d[i], '%m'), CellId = cells) pdat <- rbind(pdat, tmp) } levelplot(rr.mean) + xyplot(y ~ x|Month, data =

Re: [R-sig-Geo] Order raster images in a levelplot

Hi I give up. You'll have to approach the maintainer of the raster package directly. As far as I can see there is no difference in the order of your stack(files) and 'd'. Here is my documentation of that: ndvif <- c( # output from names(modis) "MOD13C2.A2010001.CMG_0.05_Deg_Monthly_NDVI", "MO

Re: [R-sig-Geo] Order raster images in a levelplot

Hi Could it be that you are stacking the tif files in a wrong order. So that it does not correspond to your 'd' vector? Please show us the result from list.files(pattern='NDVI.tif'). I suppose that the file names will tell us something about the date for each particular file. We need to compare

Re: [R-sig-Geo] Order raster images in a levelplot

Please remember to post to the list as well since this was there this thread startet. Dows it have something to do with you not setting the name argument of setZ()? You are doing: NDVImodis<-setZ(NDVImodis,d) names(NDVImodis)<-d Perhaps you should do: NDVImodis<-setZ(NDVImodis,d, name = "Mont

Re: [R-sig-Geo] Order raster images in a levelplot

FALSE, layer = unique(months(d)), layout = c(3,4)) Yours sincerely / Med venlig hilsen Frede Aakmann Tøgersen Specialist, M.Sc., Ph.D. Plant Performance & Modeling Technology & Service Solutions T +45 9730 5135 M +45 2547 6050 fr...@vestas.com<mailto:fr...@vestas.com> http://www.ve

Re: [R-sig-Geo] Order raster images in a levelplot

Deccember in top right panel xyplot(y ~ x|Month, data = testdata) ### use the as.table argument to get January in top left and December in bottom right xyplot(y ~ x|Month, data = testdata, as.table = TRUE) Yours sincerely / Med venlig hilsen Frede Aakmann Tøgersen Specialist, M.Sc., Ph.D. Plan

Re: [R-sig-Geo] spacetime EOF error "variable names are limited to 10000 bytes"

f@data) <- "all_layers" EOF() will work. Yours sincerely / Med venlig hilsen Frede Aakmann Tøgersen Specialist, M.Sc., Ph.D. Plant Performance & Modeling Technology & Service Solutions T +45 9730 5135 M +45 2547 6050 fr...@vestas.com http://www.vestas.com Company reg. nam

Re: [R-sig-Geo] Geolife dataset

hilsen Frede Aakmann Tøgersen Specialist, M.Sc., Ph.D. Plant Performance & Modeling Technology & Service Solutions T +45 9730 5135 M +45 2547 6050 fr...@vestas.com http://www.vestas.com Company reg. name: Vestas Wind Systems A/S This e-mail is subject to our e-mail disclaimer statement

Re: [R-sig-Geo] Calculating distance to nest with spDistsN1

= c("col_long", "col_lat"))[1,]), longlat=T)) ## have a look at the result print(trip_distances) ## unlist the trip_distances and add that as a column to dat. The order of the unlist should match the row order of dat. dat$trip_distances <- unlist(trip_distances) ## print

Re: [R-sig-Geo] R Basic Spatial Interpolation

D and perhaps in time. Do you want some kind of spatio-temporal model? Yours sincerely / Med venlig hilsen Frede Aakmann Tøgersen Specialist, M.Sc., Ph.D. Plant Performance & Modeling Technology & Service Solutions T +45 9730 5135 M +45 2547 6050 fr...@vestas.com http://www.vestas.com

Re: [R-sig-Geo] Spatial and multilevel model with kriging/interpolation in R

Yours sincerely / Med venlig hilsen Frede Aakmann Tøgersen Specialist, M.Sc., Ph.D. Plant Performance & Modeling Technology & Service Solutions T +45 9730 5135 M +45 2547 6050 fr...@vestas.com http://www.vestas.com Company reg. name: Vestas Wind Systems A/S This e-mail is subject to our e-ma

Re: [R-sig-Geo] Accessing the ids of neighbours in spdep

Hi I'm not really sure I understand you but perhaps something like this: lapply(xx, paste, collapse = ", ") Yours sincerely / Med venlig hilsen Frede Aakmann Tøgersen Specialist, M.Sc., Ph.D. Plant Performance & Modeling Technology & Service Solutions T +45 9730

Re: [R-sig-Geo] Spatial and multilevel model with kriging/interpolation in R

http://www.rni.helsinki.fi/~jmh/mrf08/R-INLA.pdf may give you some ideas. The nlme package for R can fit the same kind of models, see e.g. http://www.ats.ucla.edu/stat/r/faq/spatial_regression.htm. Yours sincerely / Med venlig hilsen Frede Aakmann Tøgersen Specialist, M.Sc., Ph.D. Plant

Re: [R-sig-Geo] knearneigh: data non-numeric

Hi You are showing metadata for gadmsimpl3_poly But the error you see is for object gadmsimpl4_poly. Does that have something to do with your problem? Yours sincerely / Med venlig hilsen Frede Aakmann Tøgersen Specialist, M.Sc., Ph.D. Plant Performance & Modeling Technology &

Re: [R-sig-Geo] 5-D Kernel density estimation in R using "kde" function

and 'xmax' is 'max(x)+Hmax*supp' where 'Hmax' is the maximum of the diagonal elements of 'H'. The grid produced is the outer product of '[xmin[1], xmax[1]]', ..., '[xmin[d], xmax[d]]'. The default 'bgridsize,

Re: [R-sig-Geo] Using spTransform() to reproduce another software package's transformation

Thanks guys. I'm a novice but still learning and your expertise is very much appreciated. Thank you very much. Br. Frede Sendt fra Samsung mobil Oprindelig meddelelse Fra: Hermann Peifer Dato:03/09/2014 19.27 (GMT+01:00) Til: "MacQueen, Don" ,roger.biv...@nhh.no Cc: r-sig-g

Re: [R-sig-Geo] Using spTransform() to reproduce another software package's transformation

WGS_1984_5 transformation? Or can you by any other means find something similar to the CRS arguments? If so we can compare the CRS arguments from R and ESRI. Yours sincerely / Med venlig hilsen Frede Aakmann Tøgersen Specialist, M.Sc., Ph.D. Plant Performance & Modeling Technology &am

Re: [R-sig-Geo] Randomly moving a locality (within set limits)

quot;iso") + xyplot(0 ~ 0, col = "red", pch = 13, cex = 2) dev.off() lats <- cut(sqrtdests00[, "Latitude"], seq(-0.1, 0.1, by = 0.05)) lons <- cut(sqrtdests00[, "Longitude"], seq(-0.1, 0.1, by = 0.05)) table(lats, lons) Yours sincerely / Med ve

Re: [R-sig-Geo] Randomly moving a locality (within set limits)

056), 90, 15) Latitude Longitude 52.2045157 0.3604334 > Now choose bearing <- runif(1, 0, 360) and distance <- runif(1, 0, 15) Please check the formulas against the homepage above. Yours sincerely / Med venlig hilsen Frede Aakmann Tøgersen Specialist, M.Sc., Ph.D. Plant P

Re: [R-sig-Geo] WriteOGR to GeoJSON loses CRS

G:4008 specifies "+proj=longlat +ellps=clrk66" as the shape file was read with. Yours sincerely / Med venlig hilsen Frede Aakmann Tøgersen Specialist, M.Sc., Ph.D. Plant Performance & Modeling Technology & Service Solutions T +45 9730 5135 M +45 2547 6050 fr...@vestas.

Re: [R-sig-Geo] WriteOGR to GeoJSON loses CRS

. What does e.g. 'gdalinfo --formats' show? Yours sincerely / Med venlig hilsen Frede Aakmann Tøgersen Specialist, M.Sc., Ph.D. Plant Performance & Modeling Technology & Service Solutions T +45 9730 5135 M +45 2547 6050 fr...@vestas.com http://www.vestas.com Company reg. name: Ves

Re: [R-sig-Geo] WriteOGR to GeoJSON loses CRS

Dooh, it must be Monday morning and not Tuesday. A copy and paste error: http://www.gdal.org/drv_geojson.html points to http://wiki.geojson.org/GeoJSON_draft_version_5. Yours sincerely / Med venlig hilsen Frede Aakmann Tøgersen Specialist, M.Sc., Ph.D. Plant Performance & Mode

Re: [R-sig-Geo] WriteOGR to GeoJSON loses CRS

Hi Google points me to http://geojson.org/geojson-spec.html. See Section 3 about CRSs. Yours sincerely / Med venlig hilsen Frede Aakmann Tøgersen Specialist, M.Sc., Ph.D. Plant Performance & Modeling Technology & Service Solutions T +45 9730 5135 M +45 2547 6050 fr...@vestas.

Re: [R-sig-Geo] Help me please

etenv(MRT_DATA_DIR = "to_some_value, probably: ~/MRT/data or something like that") Hope that helps. Yours sincerely / Med venlig hilsen Frede Aakmann Tøgersen Specialist, M.Sc., Ph.D. Plant Performance & Modeling Technology & Service Solutions T +45 9730 5135 M +45 2547 605

Re: [R-sig-Geo] reading a feature class

driver. An update.packages(oldPkgs = "rgdal") didn't give me the GDAL 1.11.0 version. I had to remove.packages("rgdal") and then install.packages("rgdal"). Yours sincerely / Med venlig hilsen Frede Aakmann Tøgersen Specialist, M.Sc., Ph.D. Plant Performan

Re: [R-sig-Geo] reading a feature class

From: Frede Aakmann Tøgersen Sent: 9. juli 2014 15:12 To: 'JOSEPH BECHARA'; 'Michael Sumner' Cc: 'RsigGeo' Subject: FW: [R-sig-Geo] reading a feature class Hi I am getting some information like this from the list: Your mail to 'R-sig-Geo' with the subject RE: [R-sig-Geo] reading a featu

Re: [R-sig-Geo] reading a feature class

Hi Joseph If you email me or otherwise make it available to me and/or the list then I'll have a look at it. Br. Frede Sendt fra Samsung mobil Oprindelig meddelelse Fra: JOSEPH BECHARA Dato:09/07/2014 08.58 (GMT+01:00) Til: Frede Aakmann Tøgersen ,r-sig-geo@r-projec

Re: [R-sig-Geo] reading a feature class

Hi I don't know the gdb file format. Is it in raster or vector format? Since you used readOGR you probably did have a look at http://www.gdal.org/ogr_formats.html before asking. So there, two ESRI formats are mentioned: personal gdb (old) and file gdb (new?). Do you what kind of gdb file you

Re: [R-sig-Geo] Assign coordinates to zip address

Perhaps Google is your friend. Google for: r zip code to location Br. Frede Sendt fra Samsung mobil Oprindelig meddelelse Fra: Rolando Valdez Dato:30/06/2014 19.04 (GMT+01:00) Til: r-sig-geo@r-project.org Emne: [R-sig-Geo] Assign coordinates to zip address Hi everyone, Is

Re: [R-sig-Geo] Set color key for stplot

nd col.regions argument to do what you want. Yours sincerely / Med venlig hilsen Frede Aakmann Tøgersen Specialist, M.Sc., Ph.D. Plant Performance & Modeling Technology & Service Solutions T +45 9730 5135 M +45 2547 6050 fr...@vestas.com http://www.vestas.com Company reg. name: Vestas W

Re: [R-sig-Geo] Fwd: Geostatistical Space–Time Model

, see case studies with R code, tutorials, etc. Perhaps this is what you need: http://www.r-inla.org/examples/case-studies/blangiardo-et-al-2012 ? I am quite sure this will keep you busy for a long time. Happy modelling ;-) Yours sincerely / Med venlig hilsen Frede Aakmann Tøgersen Speci

Re: [R-sig-Geo] import iml with readOGR

Hi I'm not close to my R right now. Can you explain to me what you think your second argument to readORG does? Br. Frede Sendt fra Samsung mobil Oprindelig meddelelse Fra: Gabriele Cozzi Dato:16/04/2014 19.11 (GMT+01:00) Til: r-sig-geo Emne: [R-sig-Geo] import iml with read

[R-sig-Geo] where does that additional outer margin comes from when using sp::image or sp::spplot on SpatialPixelsDataFRame data?

uot;iso") ## can get rid of white margin ## but not the grey par(mar=c(0,0,0,0))#, pty = "m") par('plt') ## still a gray "margin" image(xyz[1], col = mycols1, breaks = myBreaks, bg = "grey80", aspect = "iso") ## Also this is not working s