[GRASS-user] problem with vi.in.ogr

GRASS GRASS user
,
1

Hi!
cvs from 23.12.2006:

location is set on EPSG 2174 (Polish 1965 zone 4)

when I try import shape with points to location:

v.in.ogr -o dsn=/home/jarekj/shp/PunktyAZP.shp output=azp min_area=0.0001
snap=-1 --overwrite

I recive message:

No projection name! Projection parameters likely to be meaningless.
A datum name S-42 (Pulkovo_1942) was specified without transformation
parameters.

Non-interactive mode: the GRASS default for S-42 is
towgs84=28.000,-130.000,-95.000.
Over-riding projection check.
Proceeding with import...

The vector 'azp' already exists and will be overwritten.
Layer: PunktyAZP

GRASS_INFO_MESSAGE(5079,5): Importing map 3706 features...
-----------------------------------------------------
Building topology ...
Registering lines: v.in.ogr: split_q.c:65: RTreeClassify: Assertion
`!p->taken[i]' failed.

of course, the last line two lines are imported

table in databese has been created without problems
I add, that I have import lot of shapefiles to that location and problem never
occur
the shape file can be easily opended and edited in qgis on the same machine
regards\
Jarek

2

Jarek

Please send "g.proj -p" of your target GRASS location and the "ogrinfo
-al -so" of your input shapefile.

What GDAL and PROJ4 version?

Maciek

3

What GDAL and PROJ4 version?

gdal 1.4.0 RC1 (fresh) Proj4 4.5 (fresh)
fresh means not upgraded, installed on fresh system

Please send "g.proj -p" of your target GRASS location and the "ogrinfo
-al -so" of your input shapefile.

Nothing interesitng, I checked it before. I wrote so I erlier imported lot of
files from ArcView to THAT grass location without problems. I suspect file is
corrupted, but I cannot find reason: In ArcView I opened it without problems
(I have "." as a separator) in qgis (latest 0.8) also opens without problems

ogrinfo:
INFO: Open of `/home/jarekj/shp/PunktyAZP.shp'
      using driver `ESRI Shapefile' successful.

Layer name: PunktyAZP
Geometry: Point
Feature Count: 3706
Extent: (3684408.169759, 5649780.912223) - (3713994.793069, 5681324.917881)
Layer SRS WKT:
PROJCS["PUWG_65_Strefa_IV",
    GEOGCS["GCS_Pulkovo_1942",
        DATUM["Pulkovo_1942",
            SPHEROID["Krasovsky_1940",6378245.0,298.3]],
        PRIMEM["Greenwich",0.0],
        UNIT["Degree",0.0174532925199433]],
    PROJECTION["Double_Stereographic"],
    PARAMETER["False_Easting",3703000.0],
    PARAMETER["False_Northing",5627000.0],
    PARAMETER["Central_Meridian",16.67222222222222],
    PARAMETER["Scale_Factor",0.9998],
    PARAMETER["Latitude_Of_Origin",51.67083333333333],
    UNIT["Meter",1.0]]
Id: Integer (6.0)
ArkuszAZP: String (50.0)
Nr_punktu: String (50.0)

g.proj
-PROJ_INFO-------------------------------------------------
name : sterea
proj : sterea
ellps : krassovsky
lat_0 : 51.67083333333333
lon_0 : 16.67222222222222
k : 0.999800
x_0 : 3703000
y_0 : 5627000
towgs84 : 33.4,-146.6,-76.3,-0.359,-0.053,0.844,-0.84
no_defs : defined
-PROJ_UNITS------------------------------------------------
unit : metre
units : metres
meters : 1

nevermind that problem, I will exprot it from Arc As txt and import to grass
as xyz. I rather point it as an error I recived
But If you found something interesting please let me know...
regards
JArek

4

Dnia sobota, 6 stycznia 2007 19:38, Maciej Sieczka napisa³:

Jarek

Please send "g.proj -p" of your target GRASS location and the "ogrinfo
-al -so" of your input shapefile.

What GDAL and PROJ4 version?

Maciek

Reasons found but problem exist:

34 points has exactly the same coors. So topology cannot be bulid....