[GRASS5] raster inaccuracy
Glynn Clements
glynn at gclements.plus.com
Thu Nov 11 11:34:19 EST 2004
Roger Miller wrote:
> I put together several examples that can be had from
>
> http://members.spinn.net/~roger/example.tgz
>
> The tar file produces a GRASS XY location named "example." The examples are
> in the PERMANENT map set.
>
> Here is the step-by-step procedure:
>
> copy the tar file to our Grass data file hierarchy and unarchive the file
>
> start Grass, using location "example" and mapset "PERMANENT"
>
> "g.list vect" should list 6 vector files. Each file contains a simple grid
> with square cells each 5280x5280.
>
> Start with 2x2example. The procedure for all of the examples is identical.
>
> "v.info 2x2example" to get the limits of the grid.
>
> Use "g.region" to set the region to match the grid. In the case of 2x2example
> the limits are n=10560 s=0 w=0 e=10560 nsres=5280 ewres=5280. In all examples
> the resolution is 5280. The region so defined should have two rows and two
> columns.
>
> Ideally it seems like you should be able to set the region with "g.region
> vect=2x2example." That may work, but in most examples it does not. If the
> region is not set exactly to the dimensions of the grid then the next step
> will produce unwanted results.
>
> "v.to.rast i=2x2example o=2x2example"
> "d.mon start=x0"
> "d.rast -o 2x2example"
> "d.vect 2x2example"
>
> This should produce an accurate map, with the lines of the vector grid exactly
> outlining the colored cells in the raster map.
>
> "d.pan" and click on the image to recenter (e.g., click your pointer in the
> center of the northwest cell). Repeat this step several times with different
> center locations. With some centers the map is accurate. With other centers
> the map is drawn with the raster file offset by half a cell width in both the
> x and y directions. When there is an offset you can use d.where, d.what.rast
> or d.what.vect to confirm that the vector map is correctly located (southwest
> corner at x=0, y=0) and the raster map is offset.
OK; I understand you now. At first glance, this appears to be a
problem with the display library. I'll look into it, and get back to
you.
--
Glynn Clements <glynn at gclements.plus.com>
More information about the grass-dev
mailing list