Skip to content
New issue

Have a question about this project? # for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “#”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? # to your account

Datum transformation not applied in /map #87

Closed
OASteinlein opened this issue Oct 2, 2015 · 10 comments
Closed

Datum transformation not applied in /map #87

OASteinlein opened this issue Oct 2, 2015 · 10 comments
Assignees
Labels
Milestone

Comments

@OASteinlein
Copy link

In the map window, using ED50 as datum (http://epsg.io/23031/map), conversion from lat, lon to UTM results in wrong UTM easing, northing. Typing in (lat, lon) = (62.0, 3.0) in the input field at the top right, the result at the bottom of the screen is not correct. I get: Easting: 500093.28, Northing: 6874396.79

Using BlueMarble geodetic calculator I get Easting; 500000.0 Northing: 6874345.4

There should be no datum transformation involved. Using http://epsg.io/32631/map and the same parameters the result is: Easting; 500000.0 Northing: 6874180.15 which is equal to the result in BlueMarble.

What am I doing wrong with ED50?

@OASteinlein
Copy link
Author

I have tried a bit more and figured out that the coordinates in the upper right of the map window are in WGS84. The projected or converted coordinates at the bottom are in the local coordinate system. BUT: The datum transformations does not work correctly. I used http://epsg.io/4230-1613/map to convert from WGS84 to ED50 in the North Sea and specified EPSG::1613 to do the job. But, in fact the conversion used is EPSG::1133 - DMA ED50 mean!!!

@klokan
Copy link
Member

klokan commented Oct 7, 2015

As you mentioned, the coordinates in the upper right of the map window are always in WGS84 - which is probably source of confusion on the original problem reported in this ticket.

If we use http://epsg.io/4230 as the datum, and compare the transformation:

then both seems to be used correctly on the EPSG.io website...

In the second comment you reported another problem - claiming that when you use in the online interface the transformation 1613 you are getting the same results as with transformation 1133.
I can't reproduce this issue.

There is a different result for transformation 1613 and 1133. The transformation calculated on a command-line delivers the same result as on the website.

screen shot 2015-10-07 at 13 28 44

http://epsg.io/4230-1613/map
$ gdaltransform -s_srs epsg:4326 -t_srs "+proj=longlat +ellps=intl +towgs84=-90.365,-101.13,-123.384,0.333,0.077,0.894,1.994 +no_defs"
3 62
3.00176338827957 62.000460577362 -38.7094011837617

screen shot 2015-10-07 at 13 29 24

http://epsg.io/4230/map and http://epsg.io/4230-1133/map
$ gdaltransform -s_srs epsg:4326 -t_srs "+proj=longlat +ellps=intl +towgs84=-87,-98,-121,0,0,0,0 +no_defs"
3 62
3.00178077731776 62.0004616837349 -29.7889704471454

Do you claim one of these results is numerically incorrect?
If this is the case we should create a new ticket just for that - when you prove why, and provide examples of correct transformation with a reference to software which you have used...

I believe this ticket can be closed. Please confirm.

@OASteinlein
Copy link
Author

Hi Petr,

What I don’t understand is that your examples below show coordinates at the bottom of the map which are both from the EPSG:1133 transformation. I.e. 1613 is not used in the first one.
My Blue Marble calculations correspond to your command line calculations.

Regards,
Odd

Here are my calculations:

1613:
[cid:image002.png@01D1010B.0CE57F50]

1133:
[cid:image003.png@01D1010B.0CE57F50]

@klokan
Copy link
Member

klokan commented Oct 7, 2015

I see. I overooked that. The report bug is acknowledged.

The problem is in JavaScript of the /map application. Call to our API is wrong on the page :
http://epsg.io/trans?x=3&y=62&t_srs=4230
instead of:
http://epsg.io/trans?x=3&y=62&t_srs=4230-1613

BTW Could you please drag&drop the screenshots/images to GitHub - they did not go through from your email reply.

@klokan klokan self-assigned this Oct 7, 2015
@klokan
Copy link
Member

klokan commented Oct 7, 2015

BTW Do you think we could meet on Skype?

@OASteinlein
Copy link
Author

Hi Petr,
Yes, we can meet on Skype if you need to discuss this.

Here are my screenshots:

1613
1133

@klokan klokan added the bug label Oct 17, 2015
@klokan klokan changed the title Coordinate conversion in map window Datum transformation not applied in /map Oct 17, 2015
@klokan klokan added this to the v8.7.5 milestone Oct 17, 2015
@klokan klokan closed this as completed in 2cafe3c Oct 17, 2015
@klokan
Copy link
Member

klokan commented Oct 17, 2015

It is fixed now.

@OASteinlein
Copy link
Author

Hi Petr,
Sorry, but it does not work properly. If I typed in Longitude: 3 and Latitude: 60 in the map window
http://epsg.io/1613/map, the result at the bottom of the screen says 'undefined' 'undefined'.
Typing in coordinates at the bottom of the screen (lat, lon), gives NaN at the top of the screen.
I tried the same with 1133 with the same result.
Odd

@klokan
Copy link
Member

klokan commented Oct 19, 2015

Correct link (also mentioned above) is http://epsg.io/4230-1613/map and http://epsg.io/4230-1133/map
If you search for ED50 and choose the transformation - then these are the links present on the website.

This is the result:
screen shot 2015-10-19 at 09 29 52

The url you mentioned is incorrect (and not linked from the website), because 1613 it is only code for transformation - not for a coordinate system (such as ED50 EPSG:4230). But you are right - the web should report you "404 Not found" instead of showing the map interface to avoid confusion in such case - this is reported under #90.

@OASteinlein
Copy link
Author

Ok, thank you very much!

# for free to join this conversation on GitHub. Already have an account? # to comment
Labels
Development

No branches or pull requests

2 participants