Timeline for Projecting data from 'EPSG:4326' to Lambert Azimuthal Equal Area with PyProj [closed]
Current License: CC BY-SA 3.0
14 events
| when toggle format | what | by | license | comment | |
|---|---|---|---|---|---|
| Jul 12, 2017 at 21:33 | review | Reopen votes | |||
| Jul 12, 2017 at 22:07 | |||||
| Jul 12, 2017 at 21:11 | history | edited | PolyGeo♦ | CC BY-SA 3.0 |
deleted 59 characters in body
|
| Jul 12, 2017 at 13:36 | history | closed |
BradHards Bera John Powell mgri AndreJ |
Needs details or clarity | |
| Jul 12, 2017 at 13:35 | comment | added | AndreJ | If you specify a radius, you don't have an ellipsoid anymore, but a sphere. So different results are obvious. | |
| Jul 12, 2017 at 13:13 | history | edited | thanp | CC BY-SA 3.0 |
added 588 characters in body
|
| Jul 12, 2017 at 11:20 | review | Close votes | |||
| Jul 12, 2017 at 13:36 | |||||
| Jul 12, 2017 at 11:04 | comment | added | BradHards | OK, now I'm not sure what you are asking, or if it has a "right answer". Do you want to try to reword the question? | |
| Jul 12, 2017 at 10:38 | comment | added | thanp | No, for the geod.inv I am using the xi, yi ones. :) Do you have any idea on why specifying the radius or not gives different results though? I am not sure what to make of the epsg:3035 projection. Comparing the coordinates doesn't make sense, as it gives different ones anyways. Also sorry if my questions are a bit naive, I've never dealt with coordinate systems before. Thanks again! | |
| Jul 12, 2017 at 10:28 | comment | added | BradHards |
It looks like you are using the x1, y1 results - they aren't latitudes and longitudes.
|
|
| Jul 12, 2017 at 10:19 | comment | added | thanp | Thank you very much for your replies. Yes, I am entering the lon/lat pairs in the Geod.inv function. Will try with the EPSG:3035 and see what I get. The thing is though that I am interested in both areas and distances, so I was hoping that by centering the projection on my dataset I would minimize the error in the distances. | |
| Jul 12, 2017 at 10:16 | comment | added | BradHards |
It would at least give you an initial indication. You are trying an equal area projection, so it'll preserve area, not distance. By the way, pyproj.Geod.inv takes latitude and longitude pairs, not projected pairs, so the results are random.
|
|
| Jul 12, 2017 at 9:55 | comment | added | thanp | Isn't EPSG:3035 centered on Europe though instead of Italy? I have no feel for the magnitude of difference this will have on the accuracy of the coordinates however. | |
| Jul 12, 2017 at 9:40 | comment | added | BradHards | You'd probably be better off using a defined code (like EPSG:3035) than defining it yourself. You did this for 4326, just do it the same. | |
| Jul 12, 2017 at 9:19 | history | asked | thanp | CC BY-SA 3.0 |