[GRASSLIST:10307] importing big ECW file with r.in.gdal

Martin Landa landa.martin at gmail.com
Wed Feb 15 03:41:45 EST 2006


Hi *,

I am trying to import a big (~1.7GB) ECW file downloaded from

 http://www.geotorrent.org/details.php?id=30

I can visualize this raster map in QGIS without any problems. I would
like to import it into a GRASS location. I have used r.in.gdal module,
but the import was not successful. The RGB channels were generated,
for example the green one:

 +----------------------------------------------------------------------------+
  | Layer:    eu.blue                        Date: Mon Feb 13 01:31:46 2006    |
  | Mapset:   PERMANENT                      Login of Creator: martin          |
 | Location: eu_mozaika                                                       |
  | DataBase: /mnt/hdb1/grassdata                                              |
  | Title:     ( eu.blue )                                                     |
 |----------------------------------------------------------------------------|
 |                                                                            |
 |   Type of Map:  raster               Number of Categories: 0               |
 |   Data Type:    CELL                                                       |
 |   Rows:         141970                                                     |
 |   Columns:      141970                                                     |
 |   Total Cells:  -1319355580                                                |
  |        Projection: Zeměpisná šířka a délka (zone 0)                        |
 |            N: 71:26:24.830769N    S: 35:05:34.821226N   Res: 0:00:00.92167 |
 |            E: 24:20:41.939214E    W: 12:00:08.07033W   Res: 0:00:00.921674 |
 |   Range of data:    min =  0 max = 0                                       |
 |                                                                            |
 |   Data Source:                                                             |
 |                                                                            |
 |                                                                            |
 |                                                                            |
 |   Data Description:                                                        |
  |    generated by r.in.gdal                                                  |
 |                                                                            |
 |                                                                            |
 +----------------------------------------------------------------------------+

The red channel seems to be correct, the others contain only value
'0'. During import process the following errors occurred:

ERROR 1: IReadBlock failed at X offset 0, Y offset 141969
ERROR 1: GetBlockRef failed at X block offset 0, Y block offset 141969

I have also tried to convert ECW file to GTiff by using
gdal_translate, also without any success.

gdal_translate -of GTiff Europe\ Landsat742.ecw /mnt/hdb2/Europe\ Landsat742.tif
Input file size is 141970, 141970
0...10...20...30.ERROR 1: TryWinRasterIO() failed for blocked scanline.
ERROR 1: IReadBlock failed at X offset 0, Y offset 0
ERROR 1: GetBlockRef failed at X block offset 0, Y block offset 0

Any idea?, best regards

 Martin

PS: The value of total cells is wrong (should not be long long int
instead of long int)? See the attached patch.
-------------- next part --------------
A non-text attachment was scrubbed...
Name: main.c.diff.gz
Type: application/x-gzip
Size: 334 bytes
Desc: not available
Url : http://lists.osgeo.org/pipermail/grass-user/attachments/20060215/73d3baca/main.c.diff.gz


More information about the grass-user mailing list