[GRASS-user] import tif and display it

GRASS GRASS user
,
1

hi,

i have imported some tif files with r.in.gdal and displayed them in a monitor. but one large tif file can't be displayed.

r.info said that:

Type of Map: raster Number of Categories: 255
Data Type: CELL
Rows: 9091
Columns: 109091
Total Cells: 991746281
     Projection: Latitude-Longitude (zone 0)
   N: 89:59:59.197727N S: 59:59:58.117727N Res: 0:00:11.88
   E: 179:59:58.92W W: 180W Res: 0:00:00.00001
Range of data: min = 1 max = 255
                       
Data Description:
generated by r.in.gdal

the other importet tif files have 9091 rows & 9091 colums and are displayed correctly!

i set with g.region the region of the file with the information from r.info but it doesn't help - nothing is displayed
can anybody explain me because this file can't be displayed? is there a border of the filesize? i have absolutly no idea

thanks
michael

2

doppelM wrote:

i have imported some tif files with r.in.gdal and displayed them in a
monitor. but one large tif file can't be displayed.

r.info said that:

Type of Map: raster Number of Categories: 255
Data Type: CELL
Rows: 9091
Columns: 109091
Total Cells: 991746281
     Projection: Latitude-Longitude (zone 0)
   N: 89:59:59.197727N S: 59:59:58.117727N Res: 0:00:11.88
   E: 179:59:58.92W W: 180W Res: 0:00:00.00001
Range of data: min = 1 max = 255
                       
Data Description:
generated by r.in.gdal

the other importet tif files have 9091 rows & 9091 colums and are
displayed correctly!

i set with g.region the region of the file with the information from
r.info but it doesn't help - nothing is displayed
can anybody explain me because this file can't be displayed? is there
a border of the filesize? i have absolutly no idea

Look at the number of columns from r.info- the east/west region
settings are wrong, and wrapping around the whole planet.

what does 'gdalinfo' say about the tif file?

Hamish

      ____________________________________________________________________________________
Looking for last minute shopping deals?
Find them fast with Yahoo! Search. http://tools.search.yahoo.com/newsearch/category.php?category=shopping

3

hi hamish,

thank you for your answer.

Look at the number of columns from r.info- the east/west region
settings are wrong, and wrapping around the whole planet.

it is correct that the tif file is wrapping around whole the planet. i have a global map of the moon and this file covers the northern polregion of it.

what does 'gdalinfo' say about the tif file?

Size is 109091, 9091
Coordinate System is:
GEOGCS["GCS_Moon_2000",
    DATUM["unknown",
        SPHEROID["unnamed",1737400,0]],
    PRIMEM["Greenwich",0],
    UNIT["degree",0.0174532925199433]]
Origin = (-179.999999999999972,89.999777146508123)
Pixel Size = (0.003300000000000,-0.003300000000000)
Metadata:
  AREA_OR_POINT=Area
Image Structure Metadata:
  COMPRESSION=LZW
  INTERLEAVE=BAND
Corner Coordinates:
Upper Left (-180.0000000, 89.9997771) (180d 0'0.00"W, 89d59'59.20"N)
Lower Left (-180.0000000, 59.9994771) (180d 0'0.00"W, 59d59'58.12"N)
Upper Right ( 180.0003000, 89.9997771) (180d 0'1.08"E, 89d59'59.20"N)
Lower Right ( 180.0003000, 59.9994771) (180d 0'1.08"E, 59d59'58.12"N)
Center ( 0.0001500, 74.9996271) ( 0d 0'0.54"E, 74d59'58.66"N)
Band 1 Block=109091x1 Type=Byte, ColorInterp=Gray
  NoData Value=0

i hope these informations can help you so that you can help me :wink:

thanks

-------- Original-Nachricht --------

Datum: Mon, 14 Jan 2008 19:30:44 -0800 (PST)
Von: Hamish <hamish_b@yahoo.com>
An: doppelM@gmx.de, grass-user@lists.osgeo.org
Betreff: Re: [GRASS-user] import tif and display it

doppelM wrote:
> i have imported some tif files with r.in.gdal and displayed them in a
> monitor. but one large tif file can't be displayed.
>
> r.info said that:
>
> Type of Map: raster Number of Categories: 255
> Data Type: CELL
> Rows: 9091
> Columns: 109091
> Total Cells: 991746281
> Projection: Latitude-Longitude (zone 0)
> N: 89:59:59.197727N S: 59:59:58.117727N Res: 0:00:11.88
> E: 179:59:58.92W W: 180W Res: 0:00:00.00001
> Range of data: min = 1 max = 255
>
> Data Description:
> generated by r.in.gdal
>
> the other importet tif files have 9091 rows & 9091 colums and are
> displayed correctly!
>
> i set with g.region the region of the file with the information from
> r.info but it doesn't help - nothing is displayed
> can anybody explain me because this file can't be displayed? is there
> a border of the filesize? i have absolutly no idea

Look at the number of columns from r.info- the east/west region
settings are wrong, and wrapping around the whole planet.

what does 'gdalinfo' say about the tif file?

Hamish

____________________________________________________________________________________
Looking for last minute shopping deals?
Find them fast with Yahoo! Search.
http://tools.search.yahoo.com/newsearch/category.php?category=shopping

4

doppelM wrote:

> > i have imported some tif files with r.in.gdal and displayed them
> > in a monitor. but one large tif file can't be displayed.
> >
> > r.info said that:
> >
> > Type of Map: raster Number of Categories: 255
> > Data Type: CELL
> > Rows: 9091
> > Columns: 109091
> > Total Cells: 991746281
> > Projection: Latitude-Longitude (zone 0)
> > N: 89:59:59.197727N S: 59:59:58.117727N Res: 0:00:11.88
> > E: 179:59:58.92W W: 180W Res: 0:00:00.00001
> > Range of data: min = 1 max = 255
> >
> > Data Description:
> > generated by r.in.gdal
> >
> > the other importet tif files have 9091 rows & 9091 colums and are
> > displayed correctly!
> >
> > i set with g.region the region of the file with the information
> > from r.info but it doesn't help - nothing is displayed
> > can anybody explain me because this file can't be displayed? is
> > there a border of the filesize? i have absolutly no idea

Hamish:

> Look at the number of columns from r.info- the east/west region
> settings are wrong, and wrapping around the whole planet.

doppelM:

it is correct that the tif file is wrapping around whole the planet.
i have a global map of the moon and this file covers the northern
polregion of it.

> what does 'gdalinfo' say about the tif file?

Size is 109091, 9091
Coordinate System is:
GEOGCS["GCS_Moon_2000",
    DATUM["unknown",
        SPHEROID["unnamed",1737400,0]],
    PRIMEM["Greenwich",0],
    UNIT["degree",0.0174532925199433]]
Origin = (-179.999999999999972,89.999777146508123)
Pixel Size = (0.003300000000000,-0.003300000000000)
Metadata:
  AREA_OR_POINT=Area
Image Structure Metadata:
  COMPRESSION=LZW
  INTERLEAVE=BAND
Corner Coordinates:
Upper Left (-180.0000000, 89.9997771) (180d 0'0.00"W,
89d59'59.20"N)
Lower Left (-180.0000000, 59.9994771) (180d 0'0.00"W,
59d59'58.12"N)
Upper Right ( 180.0003000, 89.9997771) (180d 0'1.08"E,
89d59'59.20"N)
Lower Right ( 180.0003000, 59.9994771) (180d 0'1.08"E,
59d59'58.12"N)
Center ( 0.0001500, 74.9996271) ( 0d 0'0.54"E,
74d59'58.66"N)
Band 1 Block=109091x1 Type=Byte, ColorInterp=Gray
  NoData Value=0

so d.rast is choking on a near 4GB file (uncompressed). (426k/row?)

Can you view it zoomed in to a 9000x9000 region?

Can you break the Tif into a number smaller tiles with 'gdal_translate
-srcwin' or break the grass map up into smaller tiles with g.region +
r.mapcalc?

for an example using gdal_translate see gdal_slice.sh at
  http://hamish.bowman.googlepages.com/gpsdrivefiles#tiles

Hamish

      ____________________________________________________________________________________
Be a better friend, newshound, and
know-it-all with Yahoo! Mobile. Try it now. http://mobile.yahoo.com/;_ylt=Ahu06i62sR8HDtDypao8Wcj9tAcJ