[GRASS-user] Importing .tif

[Rich Shepard]

   I've not found specific information on importing *.tif rasters so I winged
it to import a LiDAR Digital Quad (LDQ) map. This is my work flow and how
grass7.5svn responded.

   0. There are 4 files for each of the LDQs: *.aux, *.tfw, *.tif, and
*.tif.xml.

   1. Looking at the file, *.tif.xml, I learned the projection is
FIPS 3601. A web page listing FIPS codes and their EPSG equivalents showed
the latter is 32126.

   2. The command, 'grass75 -c EPSG:32126' created a new location.

   3. r.in.gdal in=$HOME/projects/oregon/estacada-rock/data/topography/estacada-45122c3/2009_OLC_Hood-to-Coast/45122c3302.tif out=ldq
ERROR: Projection of dataset does not appear to match current location.

        Dataset PROJ_INFO is:
        name: NAD83(HARN) / Oregon North
        datum: nad83harn
        ellps: grs80
        proj: lcc
        lat_1: 46
        lat_2: 44.33333333333334
        lat_0: 43.66666666666666
        lon_0: -120.5
        x_0: 2500000
        y_0: 0
        towgs84: 0,0,0,0,0,0,0
        no_defs: defined

        ERROR: proj

   This is my first question: Since the .xml file showed the projection as
FIPS 3601, and the equivalent EPSG code is 32126 is used to create the
location, why does grass75 see a mis-match?

   4. r.in.gdal -o in=$HOME/projects/oregon/estacada-rock/data/topography/estacada-45122c3/2009_OLC_Hood-to-Coast/45122c3302.tif out=ldq
Over-riding projection check
Importing 3 raster bands...
Importing raster map <ldq.red>...
  100%
Importing raster map <ldq.green>...
  100%
Importing raster map <ldq.blue>...
  100%

   This is my second question: While I can display each map in the monitor,
they all look the same. Displaying one after the other produces no change in
appearance. So, have I correctly imported the LDQ and do the different color
bands matter?

   This LDQ is produced by the intensity of returns of the highest hits in
the area of coverage. It looks like a low resolution black-and-white aerial
photograph while generated from 1.5' (45cm) horizontal resolution LiDAR
scans (see attached screenshot.png). Really helpful in correlating specific
project areas to the elevation map and other features.

Carpe weekend,

Rich
-------------- next part --------------
A non-text attachment was scrubbed...
Name: screenshot.png
Type: image/png
Size: 268492 bytes
Desc: 
URL: <http://lists.osgeo.org/pipermail/grass-user/attachments/20180713/6cd85c37/attachment-0001.png>