[GRASS-user] r.in.gdal oddity
Brian Oney
zenlines at gmail.com
Sat Oct 29 22:58:06 EDT 2011
Hello List,
I am importing worldclim .asc files from the ipcc4 simulations
(http://ccafs-climate.org/). These are latlong files with a global
extent. I have compiled grass 6.4.2svn from source on ubuntu 11.10 and
hope that I did everything correctly, according to the wiki.
When importing, here is what I do and what happens:
r.in.gdal -o input="/D/Documents/Spatial_Datasets/Climate/WorldClim/\
hccpr_hadcm3_a2a_2020s_prec_2_5min_asc/prec_10.asc" output="prec_10"
r.info prec_1
+----------------------------------------------------------------------------+
| Layer: prec_10 at data_2020 Date: Sun Oct 30 02:58:49
2011 |
| Mapset: data_2020 Login of Creator:
root |
| Location:
worldclim |
| DataBase:
/D/GIS_Dbase |
| Title: ( prec_10
) |
| Timestamp:
none |
|----------------------------------------------------------------------------|
| |
| Type of Map: raster Number of Categories:
1185 |
| Data Type:
CELL |
| Rows:
3600 |
| Columns:
8640 |
| Total Cells:
31104000 |
| Projection:
Latitude-Longitude |
| N: 90N S: 60S Res:
0:02:30 |
| E: 179:59:59.999999W W: 180W Res:
0 |
| Range of data: min = 0 max =
1185 |
| |
| Data
Description: |
| generated by
r.in.gdal |
| |
|
Comments: |
| r.in.gdal -o
input="/D/Documents/Spatial_Datasets/Climate/WorldClim/\ |
| hccpr_hadcm3_a2a_2020s_prec_2_5min_asc/prec_10.asc"
output="prec_10" |
| |
+----------------------------------------------------------------------------+
### The resolution is obviously the problem here and I am not sure why
it does this. When I try to ### force it into the region with the -l
flag I get the same result.
### Then I tell it how it should appear. This corresponds to the regions
i.e. g.region -g
r.region map=prec_10 n=90 s=-60 e=180 w=-180
### Have a look
r.info prec_10 at data_2020
WARNING: Unable to read header file for raster map
<prec_10 at data_2020@data_2020>. Invalid format.line 9: <e-w resol:
0>
+----------------------------------------------------------------------------+
| Layer: prec_10 at data_2020 Date: Sun Oct 30 01:14:45
2011 |
| Mapset: data_2020 Login of Creator:
root |
| Location:
worldclim |
| DataBase:
/D/GIS_Dbase |
| Title: ( prec_10
) |
| Timestamp:
none |
|----------------------------------------------------------------------------|
| |
| Type of Map: raster Number of Categories:
1185 |
| Data Type:
CELL |
| |
| Data
Description: |
| generated by
r.in.gdal |
| |
|
Comments: |
| r.in.gdal -o -l
input="/D/Documents/Spatial_Datasets/Climate/WorldCl\ |
| im/hccpr_hadcm3_a2a_2020s_prec_2_5min_asc/prec_10.asc"
output="prec_\ |
|
10" |
| |
+----------------------------------------------------------------------------+
I have tried the align option with r.region and I am stuck. I have
already perused the documentation concerning the .bil file import... I
could try that but this should work too. The r.in.arc module gives me
similar results.
I really appreciate the help.
Brian
More information about the grass-user
mailing list