[GRASS-user] E00 import problem

GRASS GRASS user
1

I'm trying to import E00 files generated by the DCW for West Russia, but in the import the region settings are getting scrambled.

Here's a v.info output:

  Layer: Rivers Mapset: Russia Location: Europe Database: /mnt/sda3/richard/grassdata
  Title:
   Map scale: 1:1 Map format: native Name of creator: root
  Organization:
   Source date: Sat Nov 21 10:39:54 2009 ----------------------------------------------------------------------------
    Type of Map: vector (level: 2)
                                                                            
     Number of points: 0 Number of areas: 0 Number of lines: 59798 Number of islands: 0 Number of boundaries: 0 Number of faces: 0 Number of centroids: 0 Number of kernels: 0
                                                                            
     Map is 3D: No Number of dblinks: 1
                                                                            
           Projection: Lat/Lon N: 10478012N S: 7619798N E: 146:15E W: 38:30W
                                                                            
     Digitization threshold: 0
    Comments:
                

Any ideas?

Richard Chirgwin

2

On Sat, 21 Nov 2009, Richard Chirgwin wrote:

I'm trying to import E00 files generated by the DCW for West Russia, but
in the import the region settings are getting scrambled.

Richard,

   I noticed with both v.in.ogr and v.in.e00 the regional coordinates would
be changed from what I originally entered. My solution was to enter the
correct valuse in PERMANENT/DEFAULT_WIND and PERMANENT/WIND. While we can
do this within GRASS I just did it manually.

Rich

3

Markus,

From here:
http://www.maproom.psu.edu/cgi-bin/dcw/dcwarea.cgi?Asia

...thence:
http://www.maproom.psu.edu/cgi-bin/dcw/dcwcountry0.cgi

...selecting Drainage Network takes you to:
http://www.maproom.psu.edu/cgi-bin/dcw/dcwmanager.cgi

...and finally you can get the download link.
http://www.maproom.psu.edu/cgi-bin/dcw/dcwexport.cgi

Cheers,
RC

Markus Neteler wrote:

Please post the data set URL to the list...

Markus

On Sat, Nov 21, 2009 at 4:04 AM, Richard Chirgwin
<rchirgwin@ozemail.com.au> wrote:
  

I'm trying to import E00 files generated by the DCW for West Russia, but in
the import the region settings are getting scrambled.

Here's a v.info output:
    

Layer: Rivers
   Mapset: Russia
   Location: Europe
   Database: /mnt/sda3/richard/grassdata
   Title:
    Map scale: 1:1
   Map format: native
   Name of creator: root
   Organization:
   Source date: Sat Nov 21 10:39:54 2009
  ----------------------------------------------------------------------------
   Type of Map: vector (level: 2)

      Number of points: 0 Number of areas: 0
       Number of lines: 59798 Number of islands: 0
       Number of boundaries: 0 Number of faces: 0
       Number of centroids: 0 Number of kernels: 0

        Map is 3D: No
       Number of dblinks: 1

              Projection: Lat/Lon
                     N: 10478012N S: 7619798N
                     E: 146:15E W: 38:30W

          Digitization threshold: 0
           Comments:

Any ideas?

Richard Chirgwin
_______________________________________________
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

4

Richard,

downloaded and tested:

# uncompress
e00conv dnnet.e00 24647.e00

# convert
avcimport 24647.e00 dnnet

# look at metadata:
ogrinfo -so dnnet arc
Had to open data source read-only.
INFO: Open of `dnnet'
      using driver `AVCBin' successful.

Layer name: ARC
Geometry: Line String
Feature Count: 59798
Extent: (-4774718.500000, 7619798.000000) - (-1109733.750000, 10478012.000000)
Layer SRS WKT:
PROJCS["unnamed",
    GEOGCS["Unknown datum based upon the Clarke 1866 ellipsoid",
        DATUM["Not_specified_based_on_Clarke_1866_ellipsoid",
            SPHEROID["Clarke 1866",6378206.4,294.9786982139006,
                AUTHORITY["EPSG","7008"]],
            AUTHORITY["EPSG","6008"]],
        PRIMEM["Greenwich",0,
            AUTHORITY["EPSG","8901"]],
        UNIT["degree",0.01745329251994328,
            AUTHORITY["EPSG","9122"]],
        AUTHORITY["EPSG","4008"]],
    PROJECTION["Lambert_Conformal_Conic_2SP"],
    PARAMETER["standard_parallel_1",70],
    PARAMETER["standard_parallel_2",50],
    PARAMETER["latitude_of_origin",0],
    PARAMETER["central_meridian",110],
    PARAMETER["false_easting",0],
    PARAMETER["false_northing",0],
    UNIT["METERS",1]]
UserId: Integer (0.0)
FNODE_: Integer (0.0)
TNODE_: Integer (0.0)
LPOLY_: Integer (0.0)
RPOLY_: Integer (0.0)
LENGTH: Real (12.3)
DNNET#: Integer (5.0)
DNNET-ID: Integer (5.0)
DNLNTYPE: Integer (2.0)
DNLNSTAT: Integer (2.0)

http://www.maproom.psu.edu/cgi-bin/dcw/dcwarea.cgi?Asia
states
Projection information for Russia
    Projection LAMBERT
    Units METERS
    Spheroid CLARKE1866

So you cannot import this into a LatLong location but need to
create a Lambert projected location.

The metadata in the uncompressed E00 file contain:
PRJ 2
Projection LAMBERT
Zunits NO
Units METERS
Spheroid CLARKE1866
Xshift 0.0000000000
Yshift 0.0000000000
Parameters
70 0 0.000 /* 1st standard parallel
50 0 0.000 /* 2nd standard parallel
110 0 0.000 /* central meridian
0 0 0.000 /* latitude of projection's origin
0.00000 /* false easting (meters)
0.00000 /* false northing (meters)

With the Location manager, you will be able to define it.
The v.in.e00 will work as expected.

Markus

5

Markus,

As always, thanks!

Richard

Markus Neteler wrote:

Richard,

downloaded and tested:

# uncompress
e00conv dnnet.e00 24647.e00

# convert
avcimport 24647.e00 dnnet

# look at metadata:
ogrinfo -so dnnet arc
Had to open data source read-only.
INFO: Open of `dnnet'
      using driver `AVCBin' successful.

Layer name: ARC
Geometry: Line String
Feature Count: 59798
Extent: (-4774718.500000, 7619798.000000) - (-1109733.750000, 10478012.000000)
Layer SRS WKT:
PROJCS["unnamed",
    GEOGCS["Unknown datum based upon the Clarke 1866 ellipsoid",
        DATUM["Not_specified_based_on_Clarke_1866_ellipsoid",
            SPHEROID["Clarke 1866",6378206.4,294.9786982139006,
                AUTHORITY["EPSG","7008"]],
            AUTHORITY["EPSG","6008"]],
        PRIMEM["Greenwich",0,
            AUTHORITY["EPSG","8901"]],
        UNIT["degree",0.01745329251994328,
            AUTHORITY["EPSG","9122"]],
        AUTHORITY["EPSG","4008"]],
    PROJECTION["Lambert_Conformal_Conic_2SP"],
    PARAMETER["standard_parallel_1",70],
    PARAMETER["standard_parallel_2",50],
    PARAMETER["latitude_of_origin",0],
    PARAMETER["central_meridian",110],
    PARAMETER["false_easting",0],
    PARAMETER["false_northing",0],
    UNIT["METERS",1]]
UserId: Integer (0.0)
FNODE_: Integer (0.0)
TNODE_: Integer (0.0)
LPOLY_: Integer (0.0)
RPOLY_: Integer (0.0)
LENGTH: Real (12.3)
DNNET#: Integer (5.0)
DNNET-ID: Integer (5.0)
DNLNTYPE: Integer (2.0)
DNLNSTAT: Integer (2.0)

http://www.maproom.psu.edu/cgi-bin/dcw/dcwarea.cgi?Asia
states
Projection information for Russia
    Projection LAMBERT
    Units METERS
    Spheroid CLARKE1866

So you cannot import this into a LatLong location but need to
create a Lambert projected location.

The metadata in the uncompressed E00 file contain:
PRJ 2
Projection LAMBERT
Zunits NO
Units METERS
Spheroid CLARKE1866
Xshift 0.0000000000
Yshift 0.0000000000
Parameters
70 0 0.000 /* 1st standard parallel
50 0 0.000 /* 2nd standard parallel
110 0 0.000 /* central meridian
0 0 0.000 /* latitude of projection's origin
0.00000 /* false easting (meters)
0.00000 /* false northing (meters)

With the Location manager, you will be able to define it.
The v.in.e00 will work as expected.

Markus