similar to: problem with project command in rgdal

summary: I can display a lon-lat map on a lattice::layerplot, and I can display a Lambert conformal conic (LCC) map on a spam::image, but I can't display an LCC map on a lattice::layerplot. Example follows. What am I doing wrong? details: I've been using `lattice` (via `rasterVis`) successfully to display global atmospheric data, which works well enough (though I am definitely intrigued

Hi all, I'm faced with a problem applying the sp package: The projection argument in readShapePoly(Shapefile,proj4string="CRS class argument") e.g.: CRS("+proj=aea +lat_1=46 +lat_2=73 +lat_0=60 +lon_0=84 +x_0=0 +y_0=0 +ellps=clrk66 +units=m +no_defs") doesn't have any impact on the plotted object. I also tested the simple example: xy = cbind(x = 2 * runif(100) -

Hi, Doing more search i've discovered package RGDAL that can write a geotiff file with projection. I saved a geotiff file in UTM projection and if i read the file back in R and check the projection seems that everything is OK. But if i load the file in ArcGIS (ESRI product) i get the warning that the file is missing spatial reference so it cannot be projected, but it is displayed correctly.

Hi R users I need the map of France?s ? communes ? (towns) to build a map Is there a way to get it? More generally: How to do to get the map of a country and its different geographical levels? Best regards -- View this message in context: http://old.nabble.com/map-of-a-country-and-its-different-geographical-levels-tp26225645p26225645.html Sent from the R help mailing list archive at

I am having trouble getting the rgal package to work on Mac OS X. I suspect it has something to do with the link between the rgal package and the associated programs GDAL and PROJ4. Here is the message I get when I attempt to load rgal > library("rgdal") Error in dyn.load(file, DLLpath = DLLpath, ...) : unable to load shared object

Dear, I am trying to visualise a time-progressing line (it's supposed to represent spread patterns) using brew package and Google Earth. The idea is to have a function which takes start and end point geographic coordinates, as well as number of intervals to chop the path up, and returns the collection of points segmenting this line. Unfortunately my calculations fail for large distances,

Hi ALL, Can you please help me on how to determine the prediction error for ordinary kriging?Below are all the commands i used to generate the OK plot: rsa2 <- readShapeSpatial("residentialsa", CRS("+proj=tmerc +lat_0=39.66666666666666 +lon_0=-8.131906111111112 +k=1 +x_0=0 +y_0=0 +ellps=intl +units=m +no_defs")) x2 <- readShapeSpatial("ptna2",

Hello everyone, I want to get a 1km by lkm grid raster image using my csv data. If I call latitude=a, longitude=b and preciptation=c. a<-(1,2,3,4,5) b<-(6,7,8,9,10) c<-(10,20, 30,40, 50) Then I found an example in r help which goes like pts = read.table("file.csv",......) library(sp) library(rgdal) proj4string(pts)=CRS("+init=epsg:4326") # set it to lat-long pts =

Hi All, I have a netcdf gridded file with LCC projection. I can easily use image.plot to visualize it. However, as the axises are in X,Y, not Lat and Lon, I could not add state or country maps onto it (or lat lon information). I do have a grid2d file that describes the lat and lon for each (X,Y) grid, but the lat and lon are not regularly spaced, so I could not use image.plot. Does anyone know

Hello, We (two mac users) have been attempting to install rgdal from "http://www.kyngchaos.com/software:frameworks", given that it is not available as a binary on the CRAN (binaries) of the "Package Installer". I have also tried to solve this problem by looking on the net for an old question, and though I have found it, the answers do not help very much because our porgraming

Dear groups: I have tried to install rgdal in my Ubuntu, but it is not as straight as installing "sp". It demand the pre-installed GDAL and proj.4. I searched the source code in http://trac.osgeo.org/proj/ and http://trac.osgeo.org/gdal/wiki/DownloadSource. but don't know how to compile them. Who will give me a hint? Thanks a lot. Ethan Chang Centers for Disease Control &

Hi, I'm trying to plot some data (z) that is linked to lat&long coordinates (x&y). These coordinates are not however on a regular grid. I also have some shapefiles on which I would like to overlay the data. I can plot the shapefiles (country/city outlines) and overplot the data, but only using quilt.plot because I otherwise always get the error message that 'Error in

Hi all, I am a begginer in R and I have been trying to use the Conditioned Latin Hypercube to choose sample points only in areas close to roads due to the difficult thorough access in the study area. I could use a code to create the points throughout the area, but I need to create a code to make the R only comes up with points close to this road. I've been using the packages

Hi, I am using in my workflow gIntersection from sp package. Part of my relevant sessionInfo is: R version 4.0.3 (2020-10-10) Platform: x86_64-w64-mingw32/x64 (64-bit) Running under: Windows 10 x64 (build 17763) attached base packages: [1] stats graphics grDevices utils datasets methods base other attached packages: [1] CircStats_0.2-6 boot_1.3-25 MASS_7.3-53 stringr_1.4.0

Hello I'm working with a large hydrological data set stored in a netCDF format. The file stores x and y coordinates in the UTM projected coordinate system, yet when I use image to graphically display the z variable, the image is distorted in the sense that it does not plot the map in the correct spatial organization. I'm wondering if I need to define the projection of the netCDF file

Hello everybody, I am a new R & Grass user. I am interested in doing a kriging of the Temperature in a basin, based on pixel elevation (DTM). I use Grass 6.2.1 and R 2.5.1 and have installed the spgrass6 package. The DTM is called "DTM_trentino_100m" and is like this: north:5158000.0 south:5058500.0 east:1729700.0 west:1611600.0 rows:995 cols:1181 -9999.0 -9999.0 -9999.0

Dear all, i'm using gdistance to model animal movement across landscape. I have imported 11 rasters with roads, freeways, slope, use-of-land, lakes (...) after recoding them with GRASS with a HSI value ranging from 1 to 4. I've assigned zero to the NAs and then transformed all the rasters in TransitionLayers (function=mean,directions=4) and later summed all of them into a new transition:

I am running Samba 2.0.7 under Solaris 2.5.1 (x86 and sun4). Various shares are setup, which get connected when users log on from their NT4 PC's. No problems as far as this group of users are concerned. We also have lots of users using Network Computers (NC's) who log onto TSE/MetaFrame servers. Unfrtunately, these bunch of users are encountering problems whereby some of their drives

Hi, I want to use map("state") and have the ecoregion shape (please see link) file projected onto this. ftp://ftp.epa.gov/wed/ecoregions/cec_na/NA_CEC_Eco_Level3.zip ftp://ftp.epa.gov/wed/ecoregions/cec_na/NA_CEC_Eco_Level3.zip Could someone please show me how; I have never messed with this sort of thing. Thanks. -- View this message in context:

Hola miembros de la lista, Estoy utilizando la función levelplot del paquete rasterVis para graficar un raster y quiero controlar la altura de la referencia de color (colorkey) del mapa, pero cuando cambio valores en el argumento "height" no parece provocar cambios. Alguna idea de por qué no está funcionando? > cobertura class : RasterLayer dimensions : 780, 1296, 1010880