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

hellik · January 30, 2011, 5:55pm

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@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

Markus_GRASS · January 30, 2011, 6:34pm

On Sun, Jan 30, 2011 at 6:55 PM, Helmut Kudrnovsky <hellik@web.de> wrote:

Hi,

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

I had the same problem with a worldclim raster.

I can see two possibilities to get it import ito wgs84/ll

1) quick and easy, but maybe not correct with r.in.gdal -l

2) the hard way
Assume a precision loss somewhere in the creation of the esri binary
grid. I assume that the (ESRI) software used to create these rasters
uses single and not double precision for extends, resolution, center
coordinates etc.

Assume that the bounds are supposed to be
North: 90
South -60
East: 180
West: -180

and not
North: 90.0000078
South -60
East: 180
West: -180.0000188

More on, assume that pixel size should be exactly 00:05:00 or
0.008333333333333 and not 0.008333333767951 and that the number of
rows and cols is correct.

The assumed correct extends and resolution, preserving the grid
dimensions (rows, cols) would be
North: 90
South -60
East: 180
West: -180
res: 00:05:00

I don't know where this information is stored in an Arc/Info Binary
Grid, but this is easy to edit with a text editor in an ESRI hdr (GDAL
EHdr diver) labelled grid created with gdal_translate. The adjusted
EHdr raster should then import correctly without the -l switch. Maybe
gdal_translate -a_ullr to "Assign/override the georeferenced bounds of
the output file" can help too.

HTH,

Markus M

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@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

_______________________________________________
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

Markus_GRASS · January 30, 2011, 6:59pm

On Sun, Jan 30, 2011 at 7:34 PM, Markus Metz
<markus.metz.giswork@googlemail.com> wrote:

On Sun, Jan 30, 2011 at 6:55 PM, Helmut Kudrnovsky <hellik@web.de> wrote:

Hi,

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

I had the same problem with a worldclim raster.

I can see two possibilities to get it import ito wgs84/ll

1) quick and easy, but maybe not correct with r.in.gdal -l

2) the hard way
Assume a precision loss somewhere in the creation of the esri binary
grid. I assume that the (ESRI) software used to create these rasters
uses single and not double precision for extends, resolution, center
coordinates etc.

Assume that the bounds are supposed to be
North: 90
South -60
East: 180
West: -180

and not
North: 90.0000078
South -60
East: 180
West: -180.0000188

More on, assume that pixel size should be exactly 00:05:00 or

Oops, resolution should be 00:00:30 and not 00:05:00. Otherwise the
worldclim generic grids seem to have correct georeferencing info,
whereas the esri grids seem to have lost some precision.

0.008333333333333 and not 0.008333333767951 and that the number of
rows and cols is correct.

The assumed correct extends and resolution, preserving the grid
dimensions (rows, cols) would be
North: 90
South -60
East: 180
West: -180
res: 00:05:00

I don't know where this information is stored in an Arc/Info Binary
Grid, but this is easy to edit with a text editor in an ESRI hdr (GDAL
EHdr diver) labelled grid created with gdal_translate. The adjusted
EHdr raster should then import correctly without the -l switch. Maybe
gdal_translate -a_ullr to "Assign/override the georeferenced bounds of
the output file" can help too.

HTH,

Markus M

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@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

_______________________________________________
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

hellik · January 30, 2011, 8:49pm

On Sun, Jan 30, 2011 at 7:34 PM, Markus Metz wrote:

On Sun, Jan 30, 2011 at 6:55 PM, Helmut Kudrnovsky

Oops, resolution should be 00:00:30 and not 00:05:00. Otherwise the
worldclim generic grids seem to have correct georeferencing info,
whereas the esri grids seem to have lost some precision.

thanks for the hints.

gdal_translate -of "EHdr" didn't really help, I'll try the generic grids and ping
the authors of the worldclim-data regarding the esri data.

best regards
Helmut
___________________________________________________________
Empfehlen Sie WEB.DE DSL Ihren Freunden und Bekannten und wir
belohnen Sie mit bis zu 50,- Euro! https://freundschaftswerbung.web.de