[GRASS-user] PS.MAP order of layers

[Hamish]

Timothy wrote:
> I am using rasters, areas, borders, and points in an eps file
> using ps.map.  I want the points to be placed under an area
> layer, because I only want the relevant points to a specific
> country to be seen.  Not matter what order I put the
> commands in (vpoints before or after varea), the points appear
> on top.  I am running GRASS 6.4 in the Linux environment.
> Any thoughts on what to do to make this work?

add "masked y" to the vpoints instruction.

(the help page for ps.map is somewhat misleading about what the
masked instruction really does, and if the ordering works among
different feature types)

#spearfish example
g.region -d
ps.map out=test.ps << EOF
raster elevation.dem
vareas fields
  end
vpoints archsites
  masked y
  end
end
EOF

gv test.ps



as a last resort you can use v.select to crop out the non-
overlapping points (maybe use its geos op=disjoint option or
invert the area with v.in.region+v.overlay first); or try
v.overlay on its own for "point not in polygon".


Hamish