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

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 2 CZ06 2 7.018153 1.688086

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

': > 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 hilsen Frede Aakmann Tøgersen Specialist, M.Sc., Ph.D.

[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

[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] 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] plotKML 'kml' function with raster input

,-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 & Modeling Technology & Service Solutions T +4

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

udes 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 Techn

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

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.commailto:fr...@vestas.com http://www.vestas.comhttp://www.vestas.com/ Company reg. name: Vestas Wind Systems A/S This e

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 Modeling

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 Solutions T +45

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

small reproducible example that shows us the problem that you have (see e.g. http://www.r-bloggers.com/three-tips-for-posting-good-questions-to-r-help-and-stack-overflow). Yours sincerely / Med venlig hilsen Frede Aakmann Tøgersen Specialist, M.Sc., Ph.D. Plant Performance Modeling Technology

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 ,Michael

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

: 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 Service Solutions T +45 9730 5135 M +45 2547 6050 fr

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

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 Wind Systems A/S

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

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. Please refer to www.vestas.com

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

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://www.vestas.com

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

Signatures: x target SpatialPointsDataFrame defined SpatialPoints 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

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

) 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., Ph.D. Plant Performance Modeling Technology Service Solutions T +45

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

/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øgersen Specialist, M.Sc., Ph.D. Plant

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 Solutions T +45

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

, ...')? 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 Technology Service Solutions T +45 9730 5135 M +45 2547 6050 fr

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

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 2547 6050 fr...@vestas.com http://www.vestas.com Company reg. name: Vestas Wind Systems A/S

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

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 9730 5135 M +45 2547 6050 fr

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

in 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 Service

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

/ 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-mail disclaimer statement

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

=19,col =black, cex=1.1*1.5,lwd=.3) panel.points(mastCFD$x, mastCFD$y, pch=19, col =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

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

(SB2[[i]], meanEvade, col=c(#00ff0088, #00ff0088), main = Ev )) # subsetting using i dev.off() } 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] Export subsets of Spatial Points Data Frame to images

() } 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-mail

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

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 Systems A/S This e-mail

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

-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 Performance Modeling Technology Service Solutions T +45 9730 5135 M

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 Solutions T +45

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

the 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 Technology Service

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 Performance Modeling Technology

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

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.commailto:fr...@vestas.com http://www.vestas.comhttp://www.vestas.com/ Company reg. name: Vestas Wind

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

for more examples. Author(s): Pierre Roudier See Also: 'kml-methods' 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.commailto:fr

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

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 Service Solutions T

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

|Month, data = testdata, as.table = TRUE) 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

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

() 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. name: Vestas Wind Systems A/S This e-mail is subject to our e

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. Please

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

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-mail disclaimer

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

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 Company reg. name

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 5135 M +45 2547 6050 fr

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

in 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 Service

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

, gridsize' are d=1: 401; d=2: rep(151, 2); d=3: rep(31, 3); d=4: rep(21,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 http://www.vestas.com

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:

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

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 Service Solutions T +45 9730

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

) 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 Performance Modeling Technology

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

], 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 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] 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.com http

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 Modeling Technology

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: Vestas Wind Systems A/S

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

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. Please

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

(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 6050 fr...@vestas.com http

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 Performance Modeling Technology Service Solutions T

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-project.org

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

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] Fwd: Geostatistical Space–Time Model

, 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 Specialist, M.Sc., Ph.D. Plant

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

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

$x) y - unique(xyz$y) z - xyz$z dim(z) - c(length(x), length(y)) image.default(x, y, z, breaks = myBreaks, col = mycols1, bg = grey80) Thank you for your help. Yours sincerely / Med venlig hilsen data Frede Aakmann Tøgersen Specialist, M.Sc., Ph.D. Plant Performance Modeling Technology