similar to: units for mapproject() function result

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)

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

#################################################################################### We have used maptools to construct state, county, township, census-tract, and zipcode level R maps with an Albers Conical Equal Area projection. We would like to be able to plot the location of weather stations or other point locations on the maps. The data the point locations are in latitude-longitude units

Hi list, I have a real problem with plotting US state map. When I try to plot the northern state, there will be some blank space in the top of graph (see case 1 example), and when I plot southern states, there will be a blank space in the bottom of plot (see case 2). I spent almost 2 days to figure out a solution, but could not. Would you help me if you know what the problem is? Regards, Alireza

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

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

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

My goal is to prepare a thematic map of the US, with states shaded according to their values for a variable of interest. I would like to include an inset for Alaska in the upper left and an inset for Hawaii in the lower left. If possible, I'd like to use Albers conic projection, or something similar. Thus far I have tried using the maps package with its state database (which omits

Hi everybody, I'm trying to build a trilineare coordinate system with R. For that, I need to rotate an axis. Does anybody know, if it is possible to rotate an axis, created with the command axis(), about for instance 60 degrees? I'm looking foward to any ideas and hints and want to wish everybody a merry merry christmas. Thanks in advance. Etienne Stockhausen

## Hello.. I have asked a similar question, but this is not fixed as before. ## I am running the following using Ubuntu OS: R version 2.13.1 (2011-07-08) Copyright (C) 2011 The R Foundation for Statistical Computing ISBN 3-900051-07-0 Platform: x86_64-pc-linux-gnu (64-bit) ## when I do this: > install.packages("mapproj", dependencies=T) ## I get this: Installing package(s)

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]]

Hello: I have a problem when I tried to color the USA map with different colors. The following is my data (I only used the second column of data): alaska, 1, 2 Hawaii, 0, 0 USA, 5, 5 And here is my code: library("latticeExtra") library("mapproj") state<-read.table("C:\\usaclass.txt",sep=",") state

SUMMARY: Specific problem: I'm regridding biomass-burning emissions from a global/unprojected inventory to a regional projection (LCC over North America). I need to have boundaries for Canada, Mexico, and US (including US states), but also Caribbean and Atlantic nations (notably the Bahamas). I would also like to add Canadian provinces and Mexican states. How to put these together? General

Dear all, 1. How can I compute the intersecting area between 2 polygons ? 2. I have polygons with coordinates in decimal degrees (i.e. 13 deg 30 min = 13.5 decimal degrees). I want to compute their area and get the results in square meters or square kiometers. Can anyone give me a conversion coefficient or a pointer where I can find this information (sorry for this off topic question) ? Thanks

This may be a trivial question for some of you: Is there a way to add defined grid sizes on specific areas of a map?. I want to add 10kmx10km grids to all California coastal counties: "california,alameda", "california,contra costa", "california,del norte", "california,humboldt", "california,los angeles", "california,marin",

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,

I'm having trouble understanding the output from as.windrose(). For one thing, data on a boundary between sectors seem to be left out of the counts. I assume that explains the missing point in the output below (angle 45). Shouldn't one side of each sector interval be open, to include values such as my 45 in the example? Also, why does the angle 180 in my input apparently not result in

Hi, I'm trying to plot a field obtained from the atmospheric model WRF-NMM which uses a "Rotated Lat-Lon? map projection. The WRF documentation mentions that: ? Rotates the earth's lat/lon grid such that the intersection of the equator and prime meridian is at the center of the model domain. ? Within the rotated framework the grid spacing is constant, but in an earth-relative

Hi, I am new to R and I want to display longitude and latitude for a projected map. map.axes won't do it since it only works for unprojected map. Can anyone help me with this? Thanks, Sybil -- View this message in context: http://r.789695.n4.nabble.com/use-map-axes-for-a-projected-map-tp3302918p3302918.html Sent from the R help mailing list archive at Nabble.com.

Hi everyone, I wonder if anyone knows the reason why the outputs of the same reprojection in r and arcgis are different?. The magnitude of the change can be up to 40 km in the poles. Basically, I have a database of points equally separated by one degree over the globe. In ARCGIS, I am projecting the data in GCS-WGS-1984 and then reprojected it to Berhmann to ensure equal area