[GRASS-user] r.input and d.rast issues

Rich Shepard rshepard at appl-ecosys.com
Tue Sep 28 14:13:05 PDT 2021


Since r.in.gdal followed by r.proj are not working for me, despite success
with the rastertest location I applied Carlos' suggestion to use r.import.

Starting GRASS from /data/grassdata/project/topography/ I applied r.import
to the rastertest map, columbia_2010_e_dtm_35.tif with this result:
> r.import in=~/projects/washington/nevins-dock/data/topography/columbia_2010_e_dtm_35.tif mem=6000 out=dtm35 resamp=lanczos_f ext=region resolution=region --o
Building topology for vector map
<tmp_r_import_region_salmo_7551 at topography>...
Registering primitives...
Building areas...
  100%
Attaching islands...
  100%
Attaching centroids...
  100%
Selected PROJ pipeline:
+proj=pipeline +step +inv +proj=lcc +lat_0=45.3333333333333 +lon_0=-120.5
+lat_1=47.3333333333333 +lat_2=45.8333333333333 +x_0=500000 +y_0=0
+ellps=GRS80 +step +proj=lcc +lat_0=45.3333333333333 +lon_0=-120.5
+lat_1=45.8333333333333 +lat_2=47.3333333333333 +x_0=500000 +y_0=0
+ellps=GRS80 +step +proj=unitconvert +xy_in=m +xy_out=us-ft
************************
Reprojecting primitives ...
Building topology for vector map
<tmp_r_import_region_salmo_7551 at PERMANENT>...
Registering primitives...
Building areas...
  100%
Attaching islands...
  100%
Attaching centroids...
  100%
ERROR: Input raster does not overlap current computational region. Nothing
        to import.
ERROR: Unable to import GDAL dataset
        </home/rshepard/projects/washington/nevins-dock/data/topography/columbia_2010_e_dtm_35.tif>

I interpreted the manual page explanation of region to shift the output to
that region. Do I need to set the current region to the region of the map
I'm importing? How do I learn this before I import the .tif?

I set the region to the raster map, local_depth:
r.info map=local_depth
  +----------------------------------------------------------------------------+
  | Map:      local_depth                    Date: Fri Sep 24 15:49:01 2021    |
  | Mapset:   topography                     Login of Creator: rshepard        |
  | Location: new_nevins_dock                                                  |
  | DataBase: /data/grassdata                                                  |
  | Title:                                                                     |
  | Timestamp: none                                                            |
  |----------------------------------------------------------------------------|
  |                                                                            |
  |   Type of Map:  raster               Number of Categories: 0               |
  |   Data Type:    FCELL                Band reference: (none)                |
  |   Rows:         15265                                                      |
  |   Columns:      10768                                                      |
  |   Total Cells:  164373520                                                  |
  |        Projection: NAD83(HARN) / Washington South                          |
  |            N:     106729    S:      60934   Res:     3                     |
  |            E:    1544180    W:    1511876   Res:     3                     |
  |   Range of data:    min = NULL  max = NULL                                 |
  |                                                                            |
  |   Data Description:                                                        |
  |    generated by r.proj                                                     |
  |                                                                            |
  |   Comments:                                                                |
  |    r.proj -n location="project/ mapset="bathymetry" input="local_de\       |
  |    pth" method="nearest" memory=300                                        |
  |                                                                            |
  +----------------------------------------------------------------------------+

but when I try to display it now GRASS tells me it cannot find the map to
display:
d.rast map=local_depth
python3: can't open file '/data/grassdata/project/bathymetry/.tmp/salmo/MONITORS/wx0/render.py': [Errno 2] No such file or directory

What might be going on here?

TIA,

Rich



More information about the grass-user mailing list