

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

