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
