[GRASSLIST:6063] Re: How deep is a region?

[Radim Blazek]

Roger Bivand wrote:
> On Tue, 8 Mar 2005, Radim Blazek wrote:
> 
> 
>>Roger Bivand wrote:
>>
>>>On Mon, 7 Mar 2005, Radim Blazek wrote:
>>>
>>>
>>>
>>>>I see, r.in.gdal was fixed after beta2, should be ok in cvs, just edit 
>>>>the DEFAULT_WIND and change depths to 1.
>>>
>>>
>>>Reading the CVS shows include/gis.h struct Cell_head has int depths, the 
>>>changes in raster/r.in.gdal/main.c are setting cellhd.depths = 1.5; - why 
>>>should this work? And why 1.5 when the top is 1.0, bottom 0.0, and 
>>>resolution 1.0, all doubles? Shouldn't it just be (int) 1?
>>
>>cellhd.depths should be inited to 1, Cho followed gis/adj_cellhd.c
>>where
>>cellhd->depths =
>>  (cellhd->top - cellhd->bottom + cellhd->tb_res/2.0) / cellhd->tb_res;
>>
>>I don't know why cellhd->tb_res/2.0, it is copy-paste of cols/rows.
> 
> 
> Nor do I - I've added Cho to the CCs here in case GRASSLIST isn't enough. 
> Who else could fix and check this before tomorrow, someone who knows and 
> uses r.in.gdal?

cellhd.depths in r.in.gdal should be set to 1 without doubts.
I don't maintain raster modules.

Radim


> Roger
> 
> 
>>Radim
>>
>>
>>>Roger
>>>
>>>
>>>
>>>>Radim
>>>>
>>>>dpadam wrote:
>>>>
>>>>
>>>>>It was created using r.in.gdal; I haven't been able to do anything further
>>>>>because of the "invalid region" errors.
>>>>>
>>>>>david
>>>>>
>>>>>
>>>>>---------- Original Message -----------
>>>>>From: Radim Blazek <blazek at itc.it>
>>>>>To: David Adam <david at quercus.org>
>>>>>Cc: grasslist at baylor.edu
>>>>>Sent: Mon, 07 Mar 2005 09:39:52 +0100
>>>>>Subject: Re: [GRASSLIST:6019] Re: How deep is a region?
>>>>>
>>>>>
>>>>>
>>>>>
>>>>>>How was the region file created - what was the last module which
>>>>>>modified the region?
>>>>>>
>>>>>>Radim
>>>>>>
>>>>>>David Adam wrote:
>>>>>>
>>>>>>
>>>>>>
>>>>>>>Hamish wrote:
>>>>>>>
>>>>>>>
>>>>>>>
>>>>>>>
>>>>>>>>>I'm trying to set up a region, starting with a Landsat TIFF
>>>>>>>>>image from California.  I have imported the region with r.in.gdal,
>>>>>>>>>which produces a large file. I then close grass60 and reopen in the
>>>>>>>>>new location, at which point I get the errors
>>>>>>>>>
>>>>>>>>>
>>>>>>>>>GRASS_INFO_ERROR(5579,1): default region is invalid
>>>>>>>>>GRASS_INFO_ERROR(5579,1): line 15: <depths:     -1074916636>
>>>>>>>>>
>>>>>>>>>The relevant file for the region is
>>>>>>>>>
>>>>>>>>>DEFAULT_WIND --
>>>>>>>>>
>>>>>>>>>proj:       99
>>>>>>>>>zone:       0
>>>>>>>>>north:      213720.89780907
>>>>>>>>>south:      -6149.10219093
>>>>>>>>>east:       -168613.96484194
>>>>>>>>>west:       -404623.96484194
>>>>>>>>>cols:       7867
>>>>>>>>>rows:       7329
>>>>>>>>>e-w resol:  30
>>>>>>>>>n-s resol:  30
>>>>>>>>>top:        1
>>>>>>>>>bottom:     0
>>>>>>>>>cols3:      7867
>>>>>>>>>rows3:      7329
>>>>>>>>>depths:     -1074916636
>>>>>>>>>e-w resol3: 30
>>>>>>>>>n-s resol3: 30
>>>>>>>>>t-b resol:  1
>>>>>>>>>
>>>>>>>>>--------------
>>>>>>>>>
>>>>>>>>>Most of the fields in the file are explained in the REGION section
>>>>>>>>>of the manual, but I can't seem to find anyplace that defines or
>>>>>>>>>explains what the "depths" field is about.  But until I can get
>>>>>>>>>grass to accept the region, I can't do much of anything.
>>>>>>>>>
>>>>>>>>>All help is much appreciated...
>>>>>>>>>
>>>>>>>>
>>>>>>>>
>>>>>>>>
>>>>>>>>rows,cols, and res=30 all look right. I don't know what is causing the
>>>>>>>>depths to go nuts. Try editing the "WIND" file in your mapset to make
>>>>>>>>depths: 1?
>>>>>>>>
>>>>>>>>
>>>>>>>>maybe this works as a temporary workaround--
>>>>>>>>
>>>>>>>>try viewing the GeoTIFF meta-data with "gdalinfo".
>>>>>>>>
>>>>>>>>If the meta-data is bogus, wipe it clean with:
>>>>>>>>
>>>>>>>>tifftopnm map.tiff | pnmtotiff > map_clean.tiff
>>>>>>>>
>>>>>>>>
>>>>>>>>This will wipe all geo-coding, but you should be able to load into a X,Y
>>>>>>>>location ok and at least view the data. May rectify from there.
>>>>>>>>
>>>>>>>>
>>>>>>>>
>>>>>>>>Hamish
>>>>>>>>
>>>>>>>>
>>>>>>>>
>>>>>>>>
>>>>>>>
>>>>>>>Thanks, I'll give that a try.
>>>>>>>In response to another reply, this is happening under Grass6.0.0beta2.
>>>>>>>
>>>>>>>But one of my questions still stands: What does depth represent?  The
>>>>>>>other terms in the
>>>>>>>DEFAULT_WIND file are defined in the g.region documentation, but not
>>>>>>>depth.  If I knew
>>>>>>>what it stands for, perhaps I could guess the right value.
>>>>>>>
>>>>>>>Thanks for the help!
>>>>>>>
>>>>>>>david
>>>>>
>>>>>------- End of Original Message -------
>>>>
>>>>
>>
>