[prev in list] [next in list] [prev in thread] [next in thread] 

List:       r-sig-geo
Subject:    [R-sig-Geo] makeGrid(PBSmapping)
From:       Markus Loecher <loecher () eden ! rutgers ! edu>
Date:       2008-02-28 19:58:57
Message-ID: 20080228195805.F18A632408E () annwn13 ! rutgers ! edu
[Download RAW message or body]

Dear geo experts,
I am clearly misunderstanding the role of the projection argument in 
the wonderful utility makeGrid(PBSmapping).
I had hoped that by setting projection ="LL" the resulting grid would 
be equidistant in "real" space, and hence curvilinear in lat/lon space.
But the following code yields the identical, regular grid, 
irrespective of the projection argument:

	mypolys <- makeGrid(x= seq(-123,-122,length=10), y = seq(33, 34, 
length=10), byrow = FALSE, addSID = TRUE, projection = "LL")
	plotMap(mypolys)	
	mypolys <- makeGrid(x= seq(-123,-122,length=10), y = seq(33, 34, 
length=10), byrow = FALSE, addSID = TRUE, projection = 1)
	plotMap(mypolys)

Is there an easy way to create a rectangular grid in lat/lon space ? 
Should I first translate the points to UTM coordinates ?

Thanks!

Markus

_______________________________________________
R-sig-Geo mailing list
R-sig-Geo@stat.math.ethz.ch
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
[prev in list] [next in list] [prev in thread] [next in thread]