[R-sig-Geo] Converting UTMs to SpatialPolygons object

[seth]

        It seems my data-generating code was cut off. Here's another try:

        id = rep(seq(1,3),each=12)            # Generate polygon
IDs'
az = rep(seq(0,359,30),3)               # Generate
azimuth sequence 
c = data.frame(id
[1],az)                        # Combine in
data frame
c$dist = round(rnorm(36,10,3),1)  # Generate random distances from
center of polygon 

        xp = function(azimuth, r){cos(((90-azimuth)*pi)/180)*r}
      # Function to calculate horizontal distance from polygon
center 

        yp =  function(azimuth, r){sin(((90-azimuth)*pi)/180)*r}
      # Function to calculate vertical distance from polygon
center 

        ctr_e = rnorm(3,667314,50)       # Randomly select eastings
for polygon centers
ctr_n = rnorm(3,4784480,50)     # Randomly select northings for
polygon centers 

        c$utm_e = round(xp(c$az,c$dist),1) + ctr_e[c$id]  # Calculate
point eastings
c$utm_n = round(yp(c$az,c$dist),1) + ctr_n[c$id]  # Calculate point
northings

        myCRS = CRS("+proj=utm +zone=18 +ellps=WGS84")    # define
proj4string 

        Appreciatively, 

        Seth Bigelow



Links:
------
[1] http://sitemail.hostway.com/http:/


        [[alternative HTML version deleted]]


_______________________________________________
R-sig-Geo mailing list
R-sig-Geo@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-sig-geo