Try

r.in.gdal input=as_dem_30s.bil output=dem_asia_hydrosheds_30s
# check extents and range
r.info dem_asia_hydrosheds_30s

# set the region
g.region -p rast=dem_asia_hydrosheds_30s

r.colors map=dem_asia_hydrosheds_30s rules=srtm
r.slope.aspect elevation=dem_asia_hydrosheds_30s
slope=dem_asia_hydrosheds_30s_slope format=percent
r.colors map=dem_asia_hydrosheds_30s_slope color=slope

# check result
r.info dem_asia_hydrosheds_30s_slope

Markus M


On Sat, Jun 29, 2013 at 7:14 AM, Yann Chemin <yche...@gmail.com> wrote:
> Hi,
>
> data:
> -------
> http://hydrosheds.cr.usgs.gov/datadownload.php?reqdata=30demb
>
> script:
> --------
> r.in.bin --overwrite input=/home/yann/RSDATA/hydrosheds/as_dem_30s.bil
> output=dem_asia_hydrosheds_30s bytes=2 north=59.999999999999715
> south=-10.0 east=179.9999999999995 west=55.0 rows=8400 cols=15000
> anull=-9999
> r.null map=dem_asia_hydrosheds_30s setnull=55537
> r.colors map=dem_asia_hydrosheds_30s rules=srtm
> r.slope.aspect --overwrite elevation=dem_asia_hydrosheds_30s@PERMANENT
> slope=dem_asia_hydrosheds_30s_slope format=percent
> r.colors map=dem_asia_hydrosheds_30s_slope@PERMANENT color=slope
>
> complaint:
> --------------
> WARNING: Unable to read fp range file for
>          <dem_asia_hydrosheds_30s_slope@PERMANENT>
> Color table for raster map <dem_asia_hydrosheds_30s_slope@PERMANENT> set to
> 'slope'
>
> r.info on map:
> ------------------
> .
> ...
> Range of data:    min = -nan  max = -nan
> ...
> .
>
>
> --
> ----
> _______________________________________________
> grass-dev mailing list
> grass-dev@lists.osgeo.org
> http://lists.osgeo.org/mailman/listinfo/grass-dev
_______________________________________________
grass-dev mailing list
grass-dev@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-dev

Reply via email to