[GRASS-web] Probleme beim Transformieren von Koordinaten

Fri Sep 1 06:01:26 PDT 2017

Schütz SC wrote
> OK, I'm sorry.
> 
> My answer once more with the mailing list in Cc. Best regards, Maik Schütz
> 
> -----Ursprüngliche Nachricht-----
> Von: 

> neteler.osgeo@

>  [mailto:

> neteler.osgeo@

> ] Im Auftrag von Markus Neteler
> Gesendet: Mittwoch, 30. August 2017 23:47
> An: Schütz SC
> Betreff: Re: [GRASS-web] Probleme beim Transformieren von Koordinaten
> 
> Dear Mr. Schütz,
> 
> please be so kind to keep the mailing list in cc - at time I am unable to
> quickly help with the issue and probably our other developer is able to
> help.
> 
> Best regards,
> Markus Neteler
> 
> 
> On Tue, Aug 29, 2017 at 8:31 AM, Schütz SC <

> maik.schuetz@

> > wrote:
>> Dear Mr. Neteler,
>>
>> thank you very much for your answer.
>>
>> Example: the coordinate of the lighthouse of the isle Hiddensee
>>
>> DHDN (Bessel), Gauß-Krüger-4: 4572445.6 6052741.1;
>>
>> “m.proj” convert it with my “3. try” (last e-mail): 33378616.33 
>> 6051779.33;
>>
>> correct is: 33378506 6051816 (source:
>> https://www.deine-berge.de/Rechner/Koordinaten/Dezimal/55.012786,13.235480).
>>
>> I use GRASS GIS 7.2.1 with the integrated PROJ.4 Version.
>>
>> Thank you.
>>
>> Mit freundlichem Gruß / Kind regards
>> Maik Schütz
>>
>> Schuldt Consult Ingenieurgesellschaft mbH Alte Richtenberger Straße 31
>> D-18437 Stralsund
>> Phone: +49 (0) 38 31 - 61 29 25
>> Fax: +49 (0) 38 31 - 61 29 99
>> mail to: 

> maik.schuetz@

>>
>> Sitz der Gesellschaft: Stralsund HRB 4483
>> Geschäftsführung: Carsten Schuldt
>>
>>
>>
>> -----Ursprüngliche Nachricht-----
>> Von: 

> neteler.osgeo@

>  [mailto:

> neteler.osgeo@

> ] Im 
>> Auftrag von Markus Neteler
>> Gesendet: Donnerstag, 24. August 2017 10:08
>> An: Schütz SC
>> Cc: GRASS-web
>> Betreff: Re: [GRASS-web] Probleme beim Transformieren von Koordinaten
>>
>> Dear Mr Schütz,
>>
>> this shift is typical for a datum error - would you mind to send a test
>> point with expected result?
>>
>> Furthermore, which GRASS GIS and PROJ.4 versions do you use?
>>
>> Best regards,
>> Markus Neteler
>>
>>
>> --
>> Markus Neteler, PhD
>> http://www.mundialis.de - free data with free software 
>> http://grass.osgeo.org http://courses.neteler.org/blog
>>
>>
>> 2017-08-22 13:09 GMT+02:00 Schütz SC <

> maik.schuetz@

> >:
>>> Probleme beim Transformieren von Koordinaten
>>>
>>> - english version below -
>>>
>>> Problems with transforming coordinates
>>>
>>>
>>>
>>> Dear Ladies and Gentlemen,
>>>
>>>
>>>
>>> I’m not an expert in programming and hope, you can help me.
>>>
>>> I’ve got a coordinates.txt with coordinates in Bessel (DHDN),
>>> Gauß-Krüger-4.
>>> I want to transform the coordinates into wgs84, UTM33 (region: 
>>> Northeast Germany / Baltic Sea).
>>>
>>> I tried it with m.proj:
>>>
>>>
>>>
>>> 1.    Try:
>>>
>>> proj_in:”+proj=tmerc +name=tmerc +a=6377397.155 +es=0.0066743720 
>>> +lat_0=0.0
>>> +lon_0=12.0 +k=1.000035 +x_0=4500000.0”
>>>
>>> proj_out:”+proj=utm +name=utm +a=6378137.0 +es=0.006694380 +zone=33
>>> +unfact=1.0”
>>>
>>>
>>>
>>> 2.    Try:
>>>
>>> proj_in:“+proj=tmerc +lat_0=0 +lon_0=12 +k=1 +x_0=4500000 +y_0=0
>>> +ellps=bessel +datum=potsdam +units=m“
>>>
>>> proj_out:”+proj=utm +name=utm +a=6378137.0 +es=0.006694380 +zone=33
>>> +unfact=1.0”
>>>
>>>
>>>
>>> 3.    Try:
>>>
>>> proj_in:”+proj=tmerc +lat_0=0 +lon_0=12 +k=1.000035 +x_0=4500000 
>>> +y_0=0
>>> +ellps=bessel +datum=potsdam +units=m 
>>> +towgs84=598.1,73.7,418.2,0.202,0.045,-2.455,6.7”
>>>
>>> proj_out:”+proj=utm +name=utm +a=6378137.0 +es=0.006694380 +zone=33
>>> +unfact=1.0”
>>>
>>>
>>>
>>> It works, but it is vague. The coordinates are ~ 100.0 m far away 
>>> from the right place. I tried many options, but it didn’t function
>>> correct.
>>>
>>> Is there another possibility, to transform the coordinates more 
>>> correct in GRASS GIS?
>>>
>>> Thank you for your help.

have you seen my mail related to this issue?

the best way to define proj_in and proj_out is to find the corresponding
EPSG code (e.g. www.epsg.io).

> Example: the coordinate of the lighthouse of the isle Hiddensee
>
> DHDN (Bessel), Gauß-Krüger-4: 4572445.6 6052741.1;
>
> “m.proj” convert it with my “3. try” (last e-mail): 33378616.33
> 6051779.33;
>
> correct is: 33378506 6051816 (source:
> https://www.deine-berge.de/Rechner/Koordinaten/Dezimal/55.012786,13.235480). 

using e.g.

EPSG:32633 - WGS 84 / UTM zone 33N 
(https://epsg.io/32633):

+proj=utm +zone=33 +datum=WGS84 +units=m +no_defs

EPSG:31468 - DHDN / 3-degree Gauss-Kruger zone 4 
(http://epsg.io/31468)

+proj=tmerc +lat_0=0 +lon_0=12 +k=1 +x_0=4500000 +y_0=0 +ellps=bessel
+towgs84=598.1,73.7,418.2,0.202,0.045,-2.455,6.7 +units=m +no_defs

by

m.proj coordinates=4572445.6,6052741.1 proj_in=+proj=tmerc +lat_0=0
+lon_0=12 +k=1 +x_0=4500000 +y_0=0 +ellps=bessel
+towgs84=598.1,73.7,418.2,0.202,0.045,-2.455,6.7 +units=m +no_defs
proj_out=+proj=utm +zone=33 +datum=WGS84 +units=m +no_defs

I get

378509.81|6051810.56

it's from  "DHDN / 3-degree Gauss-Kruger zone 4" to "WGS 84 / UTM zone 33N"

> correct is: 33378506 6051816

6051810.56 vs 6051816, it looks good.

378509.81 vs (33)378506

it looks good though there are leading 33 in your results ...

-----
best regards
Helmut
--
Sent from: http://osgeo-org.1560.x6.nabble.com/Grass-Website-f4053612.html