[GRASS-user] r.in.gdal worldclim-data - G_set_window(): geographic latitude not valid for north

Helmut Kudrnovsky hellik at web.de
Sun Jan 30 12:55:42 EST 2011


Hi,

I have an esri binary grid from http://www.worldclim.org/formats in wgs84/ll 

gdalinfo (gdal18)
Driver: AIG/Arc/Info Binary Grid
Files: prec_1
 prec_1.ovr
 prec_1.aux.xml
 prec_1\dblbnd.adf
 prec_1\hdr.adf
 prec_1\log
 prec_1\metadata.xml
 prec_1\prj.adf
 prec_1\sta.adf
 prec_1\vat.adf
 prec_1\w001001.adf
 prec_1\w001001x.adf
Size is 43200, 18000
Coordinate System is:
GEOGCS["WGS 84",
 DATUM["WGS_1984",
 SPHEROID["WGS 84",6378137,298.257223563,
 AUTHORITY["EPSG","7030"]],
 TOWGS84[0,0,0,0,0,0,0],
 AUTHORITY["EPSG","6326"]],
 PRIMEM["Greenwich",0,
 AUTHORITY["EPSG","8901"]],
 UNIT["degree",0.0174532925199433,
 AUTHORITY["EPSG","9108"]],
 AUTHORITY["EPSG","4326"]]
Origin = (-180.000000000000000,90.000007823109627)
Pixel Size = (0.008333333767951,-0.008333333767951)
Corner Coordinates:
Upper Left  (-180.0000000,  90.0000078) (180d 0' 0.00"W, 90d 0' 0.03"N)
Lower Left  (-180.0000000, -60.0000000) (180d 0' 0.00"W, 60d 0' 0.00"S)
Upper Right ( 180.0000188,  90.0000078) (180d 0' 0.07"E, 90d 0' 0.03"N)
Lower Right ( 180.0000188, -60.0000000) (180d 0' 0.07"E, 60d 0' 0.00"S)
Center      (   0.0000094,  15.0000039) (  0d 0' 0.03"E, 15d 0' 0.01"N)
Band 1 Block=256x4 Type=Int16, ColorInterp=Undefined
 Description = prec_1
 Min=0.000 Max=1003.000 
 Minimum=0.000, Maximum=1003.000, Mean=53.000, StdDev=77.000
 NoData Value=-32768
 Overviews: 21600x9000, 10800x4500, 5400x2250, 2700x1125, 1350x563, 675x282, 338x141
 Metadata:
 STATISTICS_MINIMUM=0
 STATISTICS_MAXIMUM=1003
 STATISTICS_MEAN=53,870989116945
 STATISTICS_STDDEV=77,56242554647

the grass-location is following

g.region -p                                                                     
projection: 3 (Längen- und Breitengrad)
zone:       0
datum:      wgs84
ellipsoid:  wgs84
north:      35N
south:      0
west:       0
east:       35E
nsres:      1
ewres:      1
rows:       35
cols:       35
cells:      1225

if I want to import the esri binary grid with r.in.gdal, I get following message

r.in.gdal -e input=C:\data\BIOCLIM\prec_201101_ersigrid\prec_30s_esri\prec\prec_1\hdr.adf output=prec1
Projection of input dataset and current location appear to match
G_set_window(): Ungültige geographische Breite für Norden.
G_set_window(): geographic latitude not valid for north.

if I do r.in.gdal -l

r.in.gdal -l --verbose input=C:\data\BIOCLIM\prec_201101_ersigrid\prec_30s_esri\prec\prec_1 output=prec1
Map bounds have been constrained to geographic coordinates. You will almost certainly want to check map 
bounds and resolution with r.info and reset them with r.region before going any further.
Projection of input dataset and current location appear to match
Proceeding with import...
r.in.gdal komplett. Raster map <prec1> created.

+----------------------------------------------------------------------------+
 | Layer:    prec1 at data                     Date: Sun Jan 30 18:45:35 2011    |
 | Mapset:   data                           Login of Creator: syringia        |
 | Location: wgs84llwc                                                        |
 | DataBase: C:\gisdata\grassdata                                             |
 | Title:     ( prec1 )                                                       |
 | Timestamp: none                                                            |
 |----------------------------------------------------------------------------|
 |                                                                            |
 |   Type of Map:  raster               Number of Categories: 1003            |
 |   Data Type:    CELL                                                       |
 |   Rows:         18000                                                      |
 |   Columns:      43200                                                      |
 |   Total Cells:  777600000                                                  |
 |        Projection: Längen- und Breitengrad                                 |
 |            N:        90N    S:        60S   Res: 0:00:30                   |
 |            E: 179:59:59.932408W    W:       180W   Res: 0:00:00.000002     |
 |   Range of data:    min = 0  max = 1003                                    |
 |                                                                            |
 |   Data Description:                                                        |
 |    generated by r.in.gdal                                                  |
 |                                                                            |
 |   Comments:                                                                |
 |    r.in.gdal -l input="C:\data\BIOCLIM\prec_201101_ersigrid\prec_30s_es\   |
 |    ri\prec\prec_1" output="prec1"                                          |
 |                                                                            |
 +----------------------------------------------------------------------------+

what would be the best way to get this raster into a wgs84-ll-location?
maybe to crop the raster a little bit by gdaltranslate?
any hints?

best regards
Helmut
___________________________________________________________
Neu: WEB.DE De-Mail - Einfach wie E-Mail, sicher wie ein Brief!  
Jetzt De-Mail-Adresse reservieren: https://produkte.web.de/go/demail02
-------------- next part --------------
A non-text attachment was scrubbed...
Name: smime.p7s
Type: application/pkcs7-signature
Size: 1747 bytes
Desc: S/MIME Cryptographic Signature
Url : http://lists.osgeo.org/pipermail/grass-user/attachments/20110130/c464b62e/smime.bin


More information about the grass-user mailing list