Re: [R] converting lat-long coordinates to Albers Conical Equal Area coordinates

From: Roger Bivand <Roger.Bivand_at_nhh.no>
Date: Fri 10 Feb 2006 - 09:01:12 EST


On Thu, 9 Feb 2006, jatwood wrote:

> ####################################################################################
> 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 and we must convert the
> coordinates
> to the Albers Conical Equal Area coordinates.
>
> However, as of yet, we have not been able to obtain satisfactory results
> using mapproj.
> Given the specs below can someone please tell us how to use mapproj (or
> some other R
> library) to convert lat-long coordinates into the appropriate Albers
> Conical Equal Area
> coordinates?

Please look in spproj on the r-spatial sourceforge site, entrance via the Spatial Task View on CRAN. Look for function project() and visit the proj4 website to work out the arguments needed.

If this isn't sufficient, please continue this thread on the R-sig-geo mailing list.

>
> Thanks in advance
> Joe Atwood
> Montana State University
> jatwood@montana.edu
>
> ####################################################################################
> The material at the bottom of this message is the metadata from
> the shapefile we are using (as created in Arc-GIS).
> ####################################################################################
>
> -----------------------------------
> Horizontal coordinate system
> Projected coordinate system name: NAD_1927_Albers
> Geographic coordinate system name: GCS_North_American_1927
> Details
> Map Projection Name: Albers Conical Equal Area
> Standard Parallel: 29.500000
> Standard Parallel: 45.500000
> Longitude of Central Meridian: -96.000000
> Latitude of Projection Origin: 23.000000
> False Easting: 0.000000
> False Northing: 0.000000
>
>
> Planar Coordinate Information
> Planar Distance Units: meters
> Coordinate Encoding Method: coordinate pair
> Coordinate Representation
> Abscissa Resolution: 0.008192
> Ordinate Resolution: 0.008192
>
>
> Geodetic Model
> Horizontal Datum Name: North American Datum of 1927
> Ellipsoid Name: Clarke 1866
> Semi-major Axis: 6378206.400000
> Denominator of Flattening Ratio: 294.978698
> _________________
>
> Bounding coordinates
> Horizontal
> In decimal degrees
> West: -127.889829
> East: -65.349890
> North: 51.606167
> South: 22.862978
> In projected or local coordinates
> Left: -2356341.706499
> Right: 2257925.009693
> Top: 3172647.320063
> Bottom: 268329.067476
>
> ###################################################################################
>
> ______________________________________________
> R-help@stat.math.ethz.ch mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
>

-- 
Roger Bivand
Economic Geography Section, Department of Economics, Norwegian School of
Economics and Business Administration, Helleveien 30, N-5045 Bergen,
Norway. voice: +47 55 95 93 55; fax +47 55 95 95 43
e-mail: Roger.Bivand@nhh.no

______________________________________________
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Received on Fri Feb 10 10:29:08 2006

This archive was generated by hypermail 2.1.8 : Sat 11 Feb 2006 - 03:10:48 EST