[GRASS-dev] creating new location in krovak projection

Hi,
I would like to define new Location using

projection krovak
ellipsoid bessel

when I come to question about projection, I specify "krovak" and hit
Return, GRASS jumps imidietly to DEFNE THE DEFAULT REGION form (for
defining boundary coordinates and resolution).

When I try do specify new projection with g.setproj, it ends like this:

GRASS 6.3.cvs (jtsk2):~ > g.setproj

Please specify projection name
Enter 'list' for the list of available projections
Hit RETURN to cancel request

krovak

ERROR: Projection krovak is not specified in the table
GRASS 6.3.cvs (jtsk2):~ >

Howto define location using krovak projection?

Thank you

Jachym
--
Jachym Cepicky
e-mail: jachym.cepicky@centrum.cz
URL: http://les-ejk.cz
GPG: http://www.les-ejk.cz/pgp/jachym_cepicky-gpg.pub
-----------------------------------------
OFFICE:
Department of Geoinformation Technologies
Zemedelska 3
613 00, Brno
Czech Republick
e-mail: xcepicky@node.mendelu.cz
URL: http://mapserver.mendelu.cz
Tel.: +420 545 134 514

On Tue, Nov 28, 2006 at 08:42:54PM +0100, Jachym Cepicky wrote:

Hi,
I would like to define new Location using

projection krovak
ellipsoid bessel

when I come to question about projection, I specify "krovak" and hit
Return, GRASS jumps imidietly to DEFNE THE DEFAULT REGION form (for
defining boundary coordinates and resolution).

When I try do specify new projection with g.setproj, it ends like this:

GRASS 6.3.cvs (jtsk2):~ > g.setproj

Please specify projection name
Enter 'list' for the list of available projections
Hit RETURN to cancel request
>krovak
ERROR: Projection krovak is not specified in the table
GRASS 6.3.cvs (jtsk2):~ >

Howto define location using krovak projection?

Hi Jachym:

the Krovak projection is missing from:
general/g.setproj/proj-parms.table

I don't know how to add it there. It would be good to get
general/g.setproj/README
updated on how to add more projections.

Markus

hi,
together with Martin, we prepared this patch for 6.3. It itroduces new
geodetic datums for both - Czech and Slovak republic (previously, we had
to edit PROJ_INFO by hand and add "towgs84" parameter there).

not sure, if this belongs to 6.2 too, but it would be nice to have it
there.

Thanks

Jachym

On Tue, Nov 28, 2006 at 08:52:27PM +0100, Markus Neteler wrote:

On Tue, Nov 28, 2006 at 08:42:54PM +0100, Jachym Cepicky wrote:
> Hi,
> I would like to define new Location using
>
> projection krovak
> ellipsoid bessel
>
> when I come to question about projection, I specify "krovak" and hit
> Return, GRASS jumps imidietly to DEFNE THE DEFAULT REGION form (for
> defining boundary coordinates and resolution).
>
> When I try do specify new projection with g.setproj, it ends like this:
>
> GRASS 6.3.cvs (jtsk2):~ > g.setproj
>
>
> Please specify projection name
> Enter 'list' for the list of available projections
> Hit RETURN to cancel request
> >krovak
> ERROR: Projection krovak is not specified in the table
> GRASS 6.3.cvs (jtsk2):~ >
>
>
> Howto define location using krovak projection?

Hi Jachym:

the Krovak projection is missing from:
general/g.setproj/proj-parms.table

I don't know how to add it there. It would be good to get
general/g.setproj/README
updated on how to add more projections.

Markus

_______________________________________________
grass-dev mailing list
grass-dev@grass.itc.it
http://grass.itc.it/mailman/listinfo/grass-dev

--
Jachym Cepicky
e-mail: jachym.cepicky@centrum.cz
URL: http://les-ejk.cz
GPG: http://www.les-ejk.cz/pgp/jachym_cepicky-gpg.pub
-----------------------------------------
OFFICE:
Department of Geoinformation Technologies
Zemedelska 3
613 00, Brno
Czech Republick
e-mail: xcepicky@node.mendelu.cz
URL: http://mapserver.mendelu.cz
Tel.: +420 545 134 514

(attachments)

grass-krovak-2.patch (2.37 KB)

Hello Jachym
Those datum paramters are already in GRASS (and have been for several years). The datum short name is hermannskogel, not that obvious and a bit long, and not my choice but we can't change it now.

Interestingly the Slovak parameters are (very slightly) different from my source - they must cover a slightly different part of the country or have a different accuracy or something. This would need to be confirmed and documented before adding them to the file. Hopefully the site you got them from has some information?

Paul

On Fri, 1 Dec 2006, Jachym Cepicky wrote:

hi,
together with Martin, we prepared this patch for 6.3. It itroduces new
geodetic datums for both - Czech and Slovak republic (previously, we had
to edit PROJ_INFO by hand and add "towgs84" parameter there).

not sure, if this belongs to 6.2 too, but it would be nice to have it
there.

Thanks

Jachym

On Tue, Nov 28, 2006 at 08:52:27PM +0100, Markus Neteler wrote:

On Tue, Nov 28, 2006 at 08:42:54PM +0100, Jachym Cepicky wrote:

Hi,
I would like to define new Location using

projection krovak
ellipsoid bessel

when I come to question about projection, I specify "krovak" and hit
Return, GRASS jumps imidietly to DEFNE THE DEFAULT REGION form (for
defining boundary coordinates and resolution).

When I try do specify new projection with g.setproj, it ends like this:

GRASS 6.3.cvs (jtsk2):~ > g.setproj

Please specify projection name
Enter 'list' for the list of available projections
Hit RETURN to cancel request

krovak

ERROR: Projection krovak is not specified in the table
GRASS 6.3.cvs (jtsk2):~ >

Howto define location using krovak projection?

Hi Jachym:

the Krovak projection is missing from:
general/g.setproj/proj-parms.table

I don't know how to add it there. It would be good to get
general/g.setproj/README
updated on how to add more projections.

Markus

_______________________________________________
grass-dev mailing list
grass-dev@grass.itc.it
http://grass.itc.it/mailman/listinfo/grass-dev

--
Jachym Cepicky
e-mail: jachym.cepicky@centrum.cz
URL: http://les-ejk.cz
GPG: http://www.les-ejk.cz/pgp/jachym_cepicky-gpg.pub
-----------------------------------------
OFFICE:
Department of Geoinformation Technologies
Zemedelska 3
613 00, Brno
Czech Republick
e-mail: xcepicky@node.mendelu.cz
URL: http://mapserver.mendelu.cz
Tel.: +420 545 134 514

hi Paul,

On Tue, Dec 05, 2006 at 10:32:40AM +0000, Paul Kelly wrote:

Hello Jachym
Those datum paramters are already in GRASS (and have been for several
years). The datum short name is hermannskogel, not that obvious and a bit
long, and not my choice but we can't change it now.

Interestingly the Slovak parameters are (very slightly) different from my
source - they must cover a slightly different part of the country or have a
different accuracy or something. This would need to be confirmed and
documented before adding them to the file. Hopefully the site you got them
from has some information?

Paul

these parameters are from Radim's old GRASS-Krovak page, we have just
copied those parameters from that page to our wiki.

This means, that we are all adding towgs84 parameter by hand since many
years, while we just could add hermannskogel datum :frowning: That is life.

Anyway, do you think, patch for

general/g.setproj/proj-parms.table

is ok and can be applyed? It worked for me and for Martin, how desired.

Howto setup Location with Krovak projection could look like follows:

1 - Press "d" for "Other Projection"
2 - Specify "krovak" for projection name
3 - Press "y", that you want to specify gedetic datum for this location
4 - Choose "hermannskogel" from list of datum names
5 - Choose "2" for Czech republic or 3 for Slovak republic from list of
    avaliable parameter sets
6 - Plural form of units: meters
7 - Specify default region for your location

Does this look ok?

thank you

Jachym

On Fri, 1 Dec 2006, Jachym Cepicky wrote:

>hi,
>together with Martin, we prepared this patch for 6.3. It itroduces new
>geodetic datums for both - Czech and Slovak republic (previously, we had
>to edit PROJ_INFO by hand and add "towgs84" parameter there).
>
>not sure, if this belongs to 6.2 too, but it would be nice to have it
>there.
>
>Thanks
>
>Jachym
>
>On Tue, Nov 28, 2006 at 08:52:27PM +0100, Markus Neteler wrote:
>>On Tue, Nov 28, 2006 at 08:42:54PM +0100, Jachym Cepicky wrote:
>>>Hi,
>>>I would like to define new Location using
>>>
>>>projection krovak
>>>ellipsoid bessel
>>>
>>>when I come to question about projection, I specify "krovak" and hit
>>>Return, GRASS jumps imidietly to DEFNE THE DEFAULT REGION form (for
>>>defining boundary coordinates and resolution).
>>>
>>>When I try do specify new projection with g.setproj, it ends like this:
>>>
>>>GRASS 6.3.cvs (jtsk2):~ > g.setproj
>>>
>>>
>>>Please specify projection name
>>>Enter 'list' for the list of available projections
>>>Hit RETURN to cancel request
>>>>krovak
>>>ERROR: Projection krovak is not specified in the table
>>>GRASS 6.3.cvs (jtsk2):~ >
>>>
>>>
>>>Howto define location using krovak projection?
>>
>>Hi Jachym:
>>
>>the Krovak projection is missing from:
>>general/g.setproj/proj-parms.table
>>
>>I don't know how to add it there. It would be good to get
>> general/g.setproj/README
>>updated on how to add more projections.
>>
>>Markus
>>
>>_______________________________________________
>>grass-dev mailing list
>>grass-dev@grass.itc.it
>>http://grass.itc.it/mailman/listinfo/grass-dev
>
>--
>Jachym Cepicky
>e-mail: jachym.cepicky@centrum.cz
>URL: http://les-ejk.cz
>GPG: http://www.les-ejk.cz/pgp/jachym_cepicky-gpg.pub
>-----------------------------------------
>OFFICE:
>Department of Geoinformation Technologies
>Zemedelska 3
>613 00, Brno
>Czech Republick
>e-mail: xcepicky@node.mendelu.cz
>URL: http://mapserver.mendelu.cz
>Tel.: +420 545 134 514
>

--
Jachym Cepicky
e-mail: jachym.cepicky@centrum.cz
URL: http://les-ejk.cz
GPG: http://www.les-ejk.cz/pgp/jachym_cepicky-gpg.pub
-----------------------------------------
OFFICE:
Department of Geoinformation Technologies
Zemedelska 3
613 00, Brno
Czech Republick
e-mail: xcepicky@node.mendelu.cz
URL: http://mapserver.mendelu.cz
Tel.: +420 545 134 514

On Wed, 6 Dec 2006, Jachym Cepicky wrote:

hi Paul,

On Tue, Dec 05, 2006 at 10:32:40AM +0000, Paul Kelly wrote:

Hello Jachym
Those datum paramters are already in GRASS (and have been for several
years). The datum short name is hermannskogel, not that obvious and a bit
long, and not my choice but we can't change it now.

Interestingly the Slovak parameters are (very slightly) different from my
source - they must cover a slightly different part of the country or have a
different accuracy or something. This would need to be confirmed and
documented before adding them to the file. Hopefully the site you got them
from has some information?

Paul

these parameters are from Radim's old GRASS-Krovak page, we have just
copied those parameters from that page to our wiki.

Oh OK. ISTR Radim did some detailed comparisons of datum shift errors using different parameters and produced some nice maps showing the distortion. It would definitely be worthwhile checking that out and why the parameters I got from that European co-ordinate conversion site are different, if anybody has time.

This means, that we are all adding towgs84 parameter by hand since many
years, while we just could add hermannskogel datum :frowning: That is life.

Anyway, do you think, patch for

general/g.setproj/proj-parms.table

is ok and can be applyed? It worked for me and for Martin, how desired.

Yes, please go ahead - I have not followed the recent modifications to g.setproj in much depth but it seems obvious to do.

Howto setup Location with Krovak projection could look like follows:

1 - Press "d" for "Other Projection"
2 - Specify "krovak" for projection name
3 - Press "y", that you want to specify gedetic datum for this location
4 - Choose "hermannskogel" from list of datum names
5 - Choose "2" for Czech republic or 3 for Slovak republic from list of
   avaliable parameter sets
6 - Plural form of units: meters
7 - Specify default region for your location

Does this look ok?

Yes that's right.

Paul

hallo,

On Wed, Dec 06, 2006 at 06:50:24PM +0000, Paul Kelly wrote:

On Wed, 6 Dec 2006, Jachym Cepicky wrote:

Yes, please go ahead - I have not followed the recent modifications to
g.setproj in much depth but it seems obvious to do.

patch commited, only general/g.setproj/proj-params.table modified

>
>
>Howto setup Location with Krovak projection could look like follows:
>
>1 - Press "d" for "Other Projection"
>2 - Specify "krovak" for projection name
>3 - Press "y", that you want to specify gedetic datum for this location
>4 - Choose "hermannskogel" from list of datum names
>5 - Choose "2" for Czech republic or 3 for Slovak republic from list of
> avaliable parameter sets
>6 - Plural form of units: meters
>7 - Specify default region for your location
>
>Does this look ok?

Yes that's right.

Paul

thank you

Jachym

--
Jachym Cepicky
e-mail: jachym.cepicky@centrum.cz
URL: http://les-ejk.cz
GPG: http://www.les-ejk.cz/pgp/jachym_cepicky-gpg.pub
-----------------------------------------
OFFICE:
Department of Geoinformation Technologies
Zemedelska 3
613 00, Brno
Czech Republick
e-mail: xcepicky@node.mendelu.cz
URL: http://mapserver.mendelu.cz
Tel.: +420 545 134 514