I had made a number of scripts that worked fine in the 6.4.0 build 6 but appear to have problems with the same code in later revisions:
The method used to take about 6-8 secs and produce about what you would expect on the screen. Now the same calculation takes 0 seconds the file is radically changed. See the before and after below.
Any help would be greatly appreciated.
Thanks
Jim Mochel
r.mapcalc na_dem_sub_20 = na_dem_input - 20
±---------------------------------------------------------------------------+
| Layer: na_dem Date: Sat Sep 24 17:37:59 2011 |
| Mapset: PERMANENT Login of Creator: jmochel |
| Location: wgs84 |
| DataBase: /home/jmochel/grassdata |
| Title: ( na_dem ) |
| Timestamp: none |
|---|
| Type of Map: raster Number of Categories: 65535 |
| Data Type: CELL |
| Rows: 8384 |
| Columns: 9102 |
| Total Cells: 76311168 |
| Projection: x,y |
| N: 4384500 S: -3999500 Res: 1000 |
| E: 4639500 W: -4462500 Res: 1000 |
| Range of data: min = 0 max = 65535 |
| Data Description: |
| generated by r.in.gdal |
| Comments: |
| r.in.gdal input=“/home/jmochel/gis/raw-data/hydo1k/na_dem.bil” outpu\ |
| t=“na_dem” |
| ±---------------------------------------------------------------------------+ |
r.info map=na_dem_sub_20
±---------------------------------------------------------------------------+
| Layer: na_dem_sub_20 Date: Mon Sep 26 23:06:03 2011 |
| Mapset: PERMANENT Login of Creator: jmochel |
| Location: wgs84 |
| DataBase: /home/jmochel/grassdata |
| Title: ( na_dem_sub_20 ) |
| Timestamp: none |
|---|
| Type of Map: raster Number of Categories: 610 |
| Data Type: CELL |
| Rows: 1 |
| Columns: 1 |
| Total Cells: 1 |
| Projection: x,y |
| N: 1 S: 0 Res: 1 |
| E: 1 W: 0 Res: 1 |
| Range of data: min = 610 max = 610 |
| Data Description: |
| generated by r.mapcalc |
| Comments: |
| na_dem_input@PERMANENT - 20 |
| ±---------------------------------------------------------------------------+ |