[GRASS-dev] checking volume display

GRASS GIS GRASS GIS developer
1

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