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@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@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 _______________________________________________ grass-user mailing list grass-user@lists.osgeo.org https://lists.osgeo.org/mailman/listinfo/grass-user