search for: mapproject

Hello everyone, I was wondering if anyone knows how to convert the Splus mapproject( ) function for use in R, I would greatly appreciate your help. I am using the following system, R Version 1.6.1 (2002-11-01) for SunOS 5.5 and Splus Version 5.1 Release 1 for Sun SPARC, SunOS 5.5. Thanks, Helen

Does anyone know what the units are for projected coordinates obtained using mapproj's mapproject function with an Albers projection? Thanks for any and all help! Buck Stockhausen *************************************************** * Dr. William T. Stockhausen * *************************************************** * Resource Ecology and Fisheries Management * * Alask...

...d running into some problems. There's a makefile for the shared library. I suspect that I just need to change that to get it right, and everything will build properly, but I'm running into trouble. Here's the original makefile, along with my comments about changes: >SRC= mapget.c mapproject.c thin.c >OBJ= mapget.o mapproject.o thin.o >CFLAGS=-I$${R_HOME}/include The include directory changed to ${R_HOME}/src/include. >#end machine-dependent stuff > >all: maps.so world.cities.s maps.so changed to maps.dll > echo 'source("world.cities.s")'...

...d running into some problems. There's a makefile for the shared library. I suspect that I just need to change that to get it right, and everything will build properly, but I'm running into trouble. Here's the original makefile, along with my comments about changes: >SRC= mapget.c mapproject.c thin.c >OBJ= mapget.o mapproject.o thin.o >CFLAGS=-I$${R_HOME}/include The include directory changed to ${R_HOME}/src/include. >#end machine-dependent stuff > >all: maps.so world.cities.s maps.so changed to maps.dll > echo 'source("world.cities.s")'...

Dear Helpers: I'm trying to use packages "maps" and "mapdata" (see blow) to display the research resutls on map (Mid-Atlantic region). In particular, I need to mark a number of points in the map by giving their latitude and longitude information. For instance, I would like to mark a point on (long, lat) =(75.56027, 39.09271). Also, I need to mark several cities that I

I think it is easiest to describe what I want in terms of the concrete problem I have. I have data from a number of countries in each of which a sample of people was interviewed. In presenting the results in a forthcoming collaborative publication much emphasis will be placed on the multi-centre nature of the study. Although I suspect colleagues may do this with shaded maps I would prefer to

Hi R-user, I am new with R and I have a problem with the map and mapproj fonctions : with the following code : >test_proj=mapproject(LONG_d01_vec_t1, LAT_d01_vec_t1, projection="lambert", parameters=c(30,60)) >longitude_vec=test_proj$x >latitude_vec=test_proj$y >longitude_mat=matrix(longitude_vec,ncol=e_sn_d01) >latitude_mat=matrix(latitude_vec,ncol=e_sn_d01) >coordonnees=data.frame(longitude_vec,lati...

...: 1 # zvar : O3 # level : 1 ## Get US state boundaries in projection units. state.map <- maps::map( database="state", projection="lambert", par=c(33,45), plot=FALSE) # parameters to lambert: ^^^^^^^^^^^^ # see mapproj::mapproject state.map.shp <- maptools::map2SpatialLines(state.map, proj4string=lcc.crs) pdf(file=lcc.pdf) rasterVis::levelplot(o3.raster, margin=FALSE ) + latticeExtra::layer( sp::sp.lines(state.map.shp, lwd=0.8, col='darkgray')) dev.off() # change this as needed to view PDFs on your system sy...

...em(paste(pycmd,'control',"200108150000","LATITCRS",1,0,sep=" "),intern=TRUE)) lon<-as.numeric(system(paste(pycmd,'control',"200108150000","LONGICRS",1,0,sep=" "),intern=TRUE)) ## Projected coordinates coords_lamb<-mapproject(lon,lat,projection="lambert",parameters=c(60.,30.)) xs<-t(matrix(coords_lamb$x,ncol=NCC,nrow=NRC))[,1] ys<-t(matrix(coords_lamb$y,ncol=NCC,nrow=NRC))[1,]} ## Plot map filled.contour(xs,ys,Tsd[,,tstep],levels=fillevs,col=cpal,plot.axes={ map('world',projection="lamber...

I have lat/lon coordinates in Michigan GeoRef that I need to re-project so they can be shown over a Michigan county map generated from the Maps library. I'm very confused as to how to do this. Can anyone help? Thanks, Shaun [[alternative HTML version deleted]]

...5507/how-to-display-a-projected-map-on-an-rlatticelayerplot > map.IOAPI <- maps::map( > database="state", projection="lambert", par=LCC.parallels, plot=FALSE) > # parameters to lambert: ^^^^^^^^^^^^^^^^^ > # see mapproj::mapproject > map.IOAPI$x <- map.lines.coords.IOAPI.x + extents.xmin > map.IOAPI$y <- map.lines.coords.IOAPI.y + extents.ymin > map.IOAPI$range <- c( > min(map.IOAPI$x), > max(map.IOAPI$x), > min(map.IOAPI$y), > max(map.IOAPI$y)) > map.IOAPI.shp <- &g...

A non-text attachment was scrubbed... Name: not available Type: text/enriched Size: 905 bytes Desc: not available Url : https://stat.ethz.ch/pipermail/r-help/attachments/20010420/61ada248/attachment.bin