[GRASS-dev] S-JTSK strikes back

[Markus Neteler]

On Thu, May 03, 2007 at 12:09:19PM +0200, Markus Neteler wrote:

cat Area_topolcianky.prj
PROJCS["S-JTSK_Krovak",GEOGCS["GCS_S_JTSK",DATUM["D_S_JTSK",SPHEROID["Bessel_1841",6377397.155,299.1528128]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]],PROJECTION["Krovak"],PARAMETER["False_Easting",0.0],PARAMETER["False_Northing",0.0],PARAMETER["Pseudo_Standard_Parallel_1",78.5],PARAMETER["Scale_Factor",0.9999],PARAMETER["Azimuth",30.28813975277778],PARAMETER["Longitude_Of_Center",24.83333333333333],PARAMETER["Latitude_Of_Center",49.5],PARAMETER["X_Scale",1.0],PARAMETER["Y_Scale",1.0],PARAMETER["XY_Plane_Rotation",0.0],UNIT["Meter",1.0]]


g.proj wkt=Area_topolcianky.prj
WARNUNG: Datum 'Jednotne_Trigonometricke_Site_Katastralni' not recognised
         by GRASS and no parameters found.
g.proj -w
PROJCS["Krovak",
    GEOGCS["bessel",
        DATUM["unknown",
            SPHEROID["Bessel_1841",6377397.155,299.1528128]],
        PRIMEM["Greenwich",0],
        UNIT["degree",0.0174532925199433]],
    PROJECTION["Krovak"],
    PARAMETER["latitude_of_center",49.5],
    PARAMETER["longitude_of_center",24.83333333333333],
    PARAMETER["azimuth",0],
    PARAMETER["pseudo_standard_parallel_1",0],
    PARAMETER["scale_factor",0.9999],
    PARAMETER["false_easting",0],
    PARAMETER["false_northing",0],
    UNIT["Meter",1]]

It has picked up at least 0.9999 (lost pseudo_standard_parallel_1 etc).

If I do the same with OGR, I get
ogrinfo -so Area_topolcianky.shp Area_topolcianky
INFO: Open of `Area_topolcianky.shp'
      using driver `ESRI Shapefile' successful.

Layer name: Area_topolcianky
Geometry: Polygon
Feature Count: 1
Extent: (1255514.000000, 473866.000000) - (1256400.000000, 474752.000000)
Layer SRS WKT:
PROJCS["S-JTSK_Krovak",
    GEOGCS["GCS_S_JTSK",
        DATUM["Jednotne_Trigonometricke_Site_Katastralni",
            SPHEROID["Bessel_1841",6377397.155,299.1528128]],
        PRIMEM["Greenwich",0.0],
        UNIT["Degree",0.0174532925199433]],
    PROJECTION["Krovak"],
    PARAMETER["False_Easting",0.0],
    PARAMETER["False_Northing",0.0],
    PARAMETER["Pseudo_Standard_Parallel_1",78.5],
    PARAMETER["Scale_Factor",0.9999],
    PARAMETER["Azimuth",30.28813975277778],
    PARAMETER["Longitude_Of_Center",24.83333333333333],
    PARAMETER["Latitude_Of_Center",49.5],
    PARAMETER["X_Scale",1.0],
    PARAMETER["Y_Scale",1.0],
    PARAMETER["XY_Plane_Rotation",0.0],
    UNIT["Meter",1.0]]
Id: Integer (6.0)
AREA: Real (19.11)
PERIM: Real (19.11)

Apparently GRASS is still lacking something.

Markus