[GRASS-user] ogr2ogr and ortho/rgb display

[Moskovitz, Bob]

Maciek,

Can you demonstrate this method using Spearfish to the Grass newbies out there?  

Thanks in advance!

Bob Moskovitz
Seismic Hazards Zonation Program
California Geological Survey
http://gmw.consrv.ca.gov/shmp

CONFIDENTIALITY NOTICE: This communication is intended only for the use of the individual or entity to which it is addressed. This message contains information from the State of California, California Geological Survey, which may be privileged, confidential and exempt from disclosure under applicable law, including the Electronic Communications Privacy Act. If the reader of this communication is not the intended recipient, you are hereby notified that any dissemination, distribution, or copying of this communication is strictly prohibited.


-----Original Message-----
From: grassuser-bounces at grass.itc.it [mailto:grassuser-bounces at grass.itc.it]On Behalf Of Maciej Sieczka
Sent: Saturday, January 06, 2007 2:50 AM
To: Daniel Farnan
Cc: grassuser at grass.itc.it
Subject: Re: [GRASS-user] ogr2ogr and ortho/rgb display


Dylan Beaudette wrote:

> be sure to project the DEM data with gdalwarp, and then import with 
> r.in.gdal . Note that DEM projection requires some planning if you will be 
> performing any flow-related analysis. If you are looking for a 'smooth' DEM 
> then I would recommend looking into the "-rcs' parameter to gdalwarp .

Daniel,

Please note that reprojecting a floating point grid will usually yield
"stairs" distortion in the output, depending on how big is the skew
between the input and output projection. That's propably what Dylan
meant. The "stairs" distortion will be the bigger, the less smoothing
involved in the reprojection.

Although aggressive smoothing like "gdalwarp -rcs" will reduce most of
the distortion, it will also filter out the details from the input DEM
and alter the original value range (min will go higher, max will go
lower). To avoid both shortcomings ("stairs" and over-smoothing) I
usually transform the grid into vector points, reproject the points and
interpolate a grid from them, with a suitable algorithm (there are few
in GRASS; for regularly spaced points input RST or IDW do fine).

Maciek

_______________________________________________
grassuser mailing list
grassuser at grass.itc.it
http://grass.itc.it/mailman/listinfo/grassuser