[GRASS-dev] checking volume display

Sat Sep 22 17:09:41 PDT 2012

Michael 

it looks like r3.in ascii is broken in GRASS6.4.3 

I replaced the file on the web site with an unzipped file
http://skagit.meas.ncsu.edu/~helena/grasswork/JR_7408MR_2m_t70.asci

which imports OK in GRASS7
GRASS 7.0.svn (nc_spm_05):~/grasswork > r3.in.ascii JR_7408MR_2m_t70.asci out=test3drast
Loading data ... (570x594x11)
 100%

but it fails to import in grass6.4.3
GRASS 6.4.3svn (nc_spm_05):~/grasswork > r3.in.ascii JR_7408MR_2m_t70.asci out=test3drast63
Loading data ... (570x594x11)
ERROR: asciiToG3d: read failed

but file created in grass7 and exported in grass6.4.3 also fails to import
GRASS 6.4.3svn (nc_spm_05):~/grasswork > r3.out.ascii test3drast out=test3drast643out.asci dp=3
GRASS 6.4.3svn (nc_spm_05):~/grasswork > r3.in.ascii test3drast643out.asci out=test3drast643outin
Loading data ... (570x594x11)
ERROR: asciiToG3d: read failed

and file created in grass6.4.3 and exported in 6.4.3 also fails to import
r.to.rast3 JR_1974_1m,JR_1995_1m,JR_NH_1998_1m,JR_NH_19991104_1mf,JR_NH_2001_JanMarch_1m,JR_NH_20070708_1m,JR_NH_2008_1m out=jr7408test
GRASS 6.4.3svn (nc_spm_05):~/grasswork > r3.out.ascii jr7408test out=jr7408test.asci dp=2
GRASS 6.4.3svn (nc_spm_05):~/grasswork > r3.in.ascii jr7408test.asci out=jr7408testoutin
Loading data ... (570x594x11)
ERROR: asciiToG3d: read failed

that looks like a critical bug to me

Helena

Helena Mitasova
Associate Professor
Department of Marine, Earth, and Atmospheric Sciences
2800 Faucette Drive, Rm. 1125 Jordan Hall
North Carolina State University
Raleigh, NC 27695-8208
hmitaso at ncsu.edu

"All electronic mail messages in connection with State business which are sent to or received by this account are subject to the NC Public Records Law and may be disclosed to third parties.” 

On Sep 16, 2012, at 1:42 PM, Michael Barton wrote:

> Helena,
> 
> I can't seem to read the volume file (although it looks OK in a text processor).
> 
> r3.in.ascii --overwrite input=/Users/cmbarton/Downloads/JR_7408MR_2m_t70.asci output=JR_7408MR_2m_t70_3d
> Loading data ...  (570x594x11)
> ERROR: asciiToG3d: read failed
> 
> r3.in.ascii --overwrite input=/Users/cmbarton/Downloads/JR_7408MR_2m_t70.asci output=JR_7408MR_2m_t70_3d type=double tiledimension=2x2x2
> Loading data ...  (570x594x11)
> ERROR: asciiToG3d: read failed
> 
> Michael
> ____________________
> C. Michael Barton
> Director, Center for Social Dynamics & Complexity 
> Professor of Anthropology, School of Human Evolution & Social Change
> Arizona State University
> 
> voice: 	480-965-6262 (SHESC), 480-727-9746 (CSDC)
> fax:          480-965-7671 (SHESC),  480-727-0709 (CSDC)
> www: http://www.public.asu.edu/~cmbarton, http://csdc.asu.edu
> 
> 
> 
> 
> 
> 
> 
> 
> 
> 
> 
> On Sep 15, 2012, at 3:25 PM, Helena Mitasova <hmitaso at ncsu.edu>
> wrote:
> 
>> Michael,
>> 
>> here is the volume (it is the Jockey's Ridge dune)
>> http://skagit.meas.ncsu.edu/~helena/grasswork/JR_7408MR_2m_t70.asci.zip
>> and the DEM to go with it
>> http://skagit.meas.ncsu.edu/~helena/grasswork/JR_2008_ALL.asci
>> 
>> You can put it into nc_cpm_08 and set the region to
>> the 3d region
>> north:      250690
>> south:      249502
>> west:       912791
>> east:       913931
>> top:        88.00000000
>> bottom:     0.00000000
>> nsres:      2
>> nsres3:     2
>> ewres:      2
>> ewres3:     2
>> tbres:      8
>> rows:       594
>> rows3:      594
>> cols:       570
>> cols3:      570
>> depths:     11
>> 
>> you can try isosurface 15 (and above), it should look like 
>> the image attched. It is pushed up in z-position 20
>> 
>> This was done on OSX 10.6.8, I think I have 10.8 at work so I will try it there as well.
>> 
>> Helena
>> 
>> 
>> 
>> 
>> <stc_iso7408_15mNE.tif>
>> Helena Mitasova
>> Associate Professor
>> Department of Marine, Earth, and Atmospheric Sciences
>> 2800 Faucette Drive, Rm. 1125 Jordan Hall
>> North Carolina State University
>> Raleigh, NC 27695-8208
>> hmitaso at ncsu.edu
>> 
>> "All electronic mail messages in connection with State business which are sent to or received by this account are subject to the NC Public Records Law and may be disclosed to third parties.” 
>> 
>> On Sep 14, 2012, at 8:10 PM, Michael Barton wrote:
>> 
>>> Thanks. 
>>> 
>>> I've always displayed the volumes with a surface under them, but just wondered if they are needed. Probably need to think how to dispense with them as you suggest.
>>> 
>>> I want to see if your volumes give me the same crash as mine. I could display my volumes last fall. So it is either 1) something that has changed in the wxNVIZ code since then or 2) something with OS X 10.8.
>>> 
>>> Michael
>>> ______________________________
>>> C. Michael Barton 
>>> Director, Center for Social Dynamics & Complexity
>>> Professor of Anthropology, School of Human Evolution & Social Change
>>> Arizona State University
>>> Tempe, AZ  85287-2402
>>> USA
>>> 
>>> voice: 	480-965-6262 (SHESC), 480-727-9746 (CSDC)
>>> fax:          480-965-7671(SHESC), 480-727-0709 (CSDC)
>>> www: 	http://csdc.asu.edu, http://shesc.asu.edu
>>> 		http://www.public.asu.edu/~cmbarton
>>> 
>>> On Sep 14, 2012, at 3:34 PM, Helena Mitasova <hmitaso at ncsu.edu>
>>> wrote:
>>> 
>>>> 
>>>> On Sep 14, 2012, at 3:15 PM, Michael Barton wrote:
>>>> 
>>>>> Helena,
>>>>> 
>>>>> Can you provide me with a volume that displays well for you on some GRASS version? I'd like to test it on my versions here. They always crash with my files now (didn't used to).
>>>> 
>>>> I will do that
>>>> 
>>>>> Also, do you need to have a surface displayed to also display a volume?
>>>> 
>>>> yes - if you do not have a surface, the vertical scaling cannot be done and that may be crashing wxnviz.
>>>> you can switch the surface off after you have everything loaded.
>>>> This also crashes the old nviz if you try to visualize 3D vector only (e.g. 3d points or TIN).
>>>> 
>>>> I am sure this can be resolved by using the vertical extent (top-bottom) for scaling but it probably requires some work.
>>>> 
>>>> Helena
>>>>> 
>>>>> Thanks
>>>>> Michael
>>>>> ______________________________
>>>>> C. Michael Barton 
>>>>> Director, Center for Social Dynamics & Complexity
>>>>> Professor of Anthropology, School of Human Evolution & Social Change
>>>>> Arizona State University
>>>>> Tempe, AZ  85287-2402
>>>>> USA
>>>>> 
>>>>> voice: 
>>>>> 480-965-6262 (SHESC), 480-727-9746 (CSDC)
>>>>> fax:          480-965-7671(SHESC), 480-727-0709 (CSDC)
>>>>> www: 
>>>>> http://csdc.asu.edu, http://shesc.asu.edu
>>>>> http://www.public.asu.edu/~cmbarton
>>>>> 
>>>> 
>>> 
>> 
> 