[GRASSLIST:924] interpolation & shading problem/question
Brent Wood
b.wood at niwa.co.nz
Thu Aug 7 22:43:09 EDT 2003
I'm slowly getting to grips with the idiosyncrasies of GRASS, having used
other GIS & mapping tools for several years.
My current problem lies in getting a good looking hillshade like the
spearfish dataset provides.
Hopefully I've included all the relevannt info below....
I have imported an xyz dataset into GRASS with:
cat mandes_.LI | awk '{print $2, $3, $4, $4}' | s.in.ascii andes d=3
The file mandes_.LI comprises N X Y Z values, where N is a record ID that
I'm omitting when I build the site file. The Z value is depth in
meters, X & Y are long/latitude in decimal degrees.
The file imports to a site file as I expect, and I then interpolate it to
build a raster with
s.surf.rst input=andes5m elev=andes5m tension=40 smooth=10 npmin=50
This displays fine as a raster, showing the seabed topography I expect.
I then try to get GRASS to display this with terrain (aspect) shading,
using the GRASS script below to create a new raster andes5m.shade
shade.rel.sh altitude=20 azimuth=270 elevation=andes5m
This is then displayed using
d.his h=andes5m i=andes5m.shade
This picture is very grainy, the overall impression of shading is OK, but
the map appears to be built of bricks (about 32x32). I use the term bricks
because each has its own shading within its boundary, so pixel is not
appropriate, each brick comprises several pixels. I have tried changing
tension & smoothing parameters, without making much difference.
The output below describes the elevation raster. If anyone can offer ay
advice to help resolve this, I would appreciate it.
Thanks,
Brent Wood
+-----------------------------------------------------------------------+
| Layer: andes5m Date: Fri Aug 8 02:19:48 2003 |
| Mapset: PERMANENT Login of Creator: baw |
| Location: andes |
| DataBase: /home/baw/grassdb |
| Title: ( andes5m ) |
|-----------------------------------------------------------------------|
| |
| Type of Map: raster Number of Categories: 255 |
| Data Type: FCELL |
| Rows: 1548 |
| Columns: 2268 |
| Total Cells: 3510864 |
| Projection: Latitude-Longitude (zone 0) |
| N: 44:05:23.28S S: 44:14:40.56S Res: 0:00:00.36 |
| E: 174:22:53.04W W: 174:36:29.52W Res: 0:00:00.36 |
| Range of data: min = -1470.000000 max = 0.000000 |
| |
| Data Source: |
| site file andes5m |
| |
| |
| Data Description: |
| generated by s.surf.rst |
| |
| Comments: |
| tension=40.000000, smoothing=10.000000 |
| dnorm=0.003973, dmin=0.000050, zmult=1.000000 |
| segmax=40, npmin=50, rmsdevi=3.464404 |
| zmin_data=-1470.000000, zmax_data=-485.000000 |
| zmin_int=-1470.000000, zmax_int=-489.320737 |
| @@t |
| ]@(@@t |
| ]@$ |
| |
+-----------------------------------------------------------------------+
More information about the grass-user
mailing list