project                package:rgdal                R Documentation

_P_r_o_j_e_c_t_i_o_n _o_f _c_o_o_r_d_i_n_a_t_e _m_a_t_r_i_c_e_s

_D_e_s_c_r_i_p_t_i_o_n:

     Interface to the PROJ.4 library of projection functions for
     geographical position data, no datum transformation possible. Use
     'spTransform()' for extended support.

_U_s_a_g_e:

     project(xy, proj, inv = FALSE)

_A_r_g_u_m_e_n_t_s:

      xy: 2-column matrix of coordinates 

    proj: character string of projection arguments; the arguments must
          be entered exactly as in the PROJ.4 documentation, in
          particular there cannot be any white space in +<arg>=<value>
          strings, and successive such strings can only be separated by
          blanks. 

     inv: default FALSE, if TRUE inverse projection to geographical
          coordinates 

_D_e_t_a_i_l_s:

     Full details of projection arguments available from website below,
     and examples in file "epsg" in the data directory installed with
     PROJ.4.

_V_a_l_u_e:

     A two column matrix with projected coordinates.

_N_o_t_e:

     The locations of Hawaii and Alaska in the data source are (putting
     it mildly) arbitrary, please avoid airlines using these positions.

_A_u_t_h_o_r(_s):

     Barry Rowlingson, Roger Bivand Roger.Bivand@nhh.no

_R_e_f_e_r_e_n_c_e_s:

     <URL: http://proj.maptools.org/>

_S_e_e _A_l_s_o:

     'CRS-class', 'spTransform-methods'

_E_x_a_m_p_l_e_s:

     data(state)
     res <- project(cbind(state.center$x, state.center$y), "+proj=lcc +lat_1=48 +lat_2=33 +lon_0=-100")
     res1 <- project(res, "+proj=lcc +lat_1=48 +lat_2=33 +lon_0=-100", inv=TRUE)
     summary(res1 - cbind(state.center$x, state.center$y))
     plot(cbind(state.center$x, state.center$y), asp=1, type="n")
     text(cbind(state.center$x, state.center$y), state.abb)
     plot(res, asp=1, type="n")
     text(res, state.abb)

