[GRASS-user] reading/converting Z-map format with Grass/Gdal

[Szilard Albert]

Hello,
Further to my question:
I have updated to gdal 1.9.0, and tried to read in my file as below:

r.in.poly input=/home/dgeo1/BasementTWT.asc output=basment

Writing raster map...
(Tue Mar 25 22:33:28 2014) Command finished (0 sec)

Item is added into the tree under map layers, but does not render.
Metadata info returns:


r.info map=basment at PERMANENT

 +----------------------------------------------------------------------------+
 | Layer:    basment at PERMANENT              Date: Tue Mar 25 22:33:28 2014
   |
 | Mapset:   PERMANENT                      Login of Creator: dgeo1
  |
 | Location: newLocation1
  |
 | DataBase: /home/dgeo1/gis
   |
 | Title:     ( basment )
  |
 | Timestamp: none
   |
 |----------------------------------------------------------------------------|
 |
   |
 |   Type of Map:  raster               Number of Categories: 0
  |
 |   Data Type:    CELL
  |
 |   Rows:         770
   |
 |   Columns:      272
   |
 |   Total Cells:  209440
  |
 |        Projection: x,y
  |
 |            N: 8949433.59941984    S: 8381239.01898174   Res:
737.91503953  |
 |            E: 774072.87646697    W: 573359.98571481   Res: 737.91503953
   |
 |   Range of data:    min = NULL  max = NULL
  |
 |
   |
 |   Data Description:
   |
 |    generated by r.in.poly
   |
 |
   |
 |   Comments:
   |
 |    r.in.poly input="/home/dgeo1/BasementTWT.asc" output="basment" rows=\
  |
 |    4096
   |
 |
   |
 +----------------------------------------------------------------------------+
(Tue Mar 25 22:36:49 2014) Command finished (0 sec)

Tried with max rows of 770, same result, range of data NULL.

Any suggestions would be appreciated,
regards,
Szilard



On 25 March 2014 14:05, Szilard Albert <szilard.albert at gmail.com> wrote:

> Dear list,
>
> I have a Z-map format text file, representing a geological horizon.
> I need to convert it to some other format, something like a simple XYZ or
> XYT
> text file.
>
> Example of format is below.
>
> Any suggestions if I can achieve this with Grass/Gdal?
> regards,
> Szilard
>
>
>
> !
> !     ZIMS FILE NAME :  242grid
> !     FORMATTED FILE CREATION DATE:  FEB 19 14
> !     FORMATTED FILE CREATION TIME:  11:35
> !
> @242grid HEADER                          ,      GRID,         5
>       15,0.1000000E+31  ,            ,      7,      1
>      624,    663,   578000.0    ,   842800.0    ,   8715600.    ,
> 8964800.
>     30000.00    ,         0.    ,         0.
> @
> 0.1000000E+31  0.1000000E+31  0.1000000E+31  0.1000000E+31  0.1000000E+31
> 0.1000000E+31  0.1000000E+31  0.1000000E+31  0.1000000E+31  0.1000000E+31
> 0.1000000E+31  0.1000000E+31  0.1000000E+31  0.1000000E+31  0.1000000E+31
> 0.1000000E+31  0.1000000E+31  0.1000000E+31  0.1000000E+31  0.1000000E+31
> 0.1000000E+31  0.1000000E+31  0.1000000E+31  0.1000000E+31  0.1000000E+31
> 0.1000000E+31  0.1000000E+31  0.1000000E+31  0.1000000E+31  0.1000000E+31
> 0.1000000E+31  0.1000000E+31  0.1000000E+31  0.1000000E+31  0.1000000E+31
> 0.1000000E+31  0.1000000E+31  0.1000000E+31  0.1000000E+31  0.1000000E+31
> 0.1000000E+31  0.1000000E+31  0.1000000E+31  0.1000000E+31  0.1000000E+31
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/grass-user/attachments/20140325/20422811/attachment.html>