Friday, 12 January 2018

coordinate system - Basic proj4 conversion UTM grid -> longlat and back



I'm a GIS novice, and am struggling to understand how to drive proj4 to accomplish what seems like a simple thing. I've googled for a couple of days, but everything I find seems to presuppose quite advanced knowledge, and leaves many things unsaid. Stack Overflow couldn't clear it up for me, either..


I have a conversion between UTM easting/northing and lat/long to do in both directions. Once I get this working, I'll add some different projections/datums.


I'm working on my recently acquired knowledge/assumptions:



  1. This conversion will require a datum translation, because WGS84 (standard datum for UTM) has it's prime meridian approximately 100m east of Greenwich (where longitude is zero).

  2. I don't need to worry about ellipsoids, because the lat/long are geographic coordinates, rather than projected ones.

  3. I can either start the conversions using one string (as used in the command-line instantiation), or can pre-assemble a source and destination object, and call pj_fwd(data, ref) directly.


Another issue is that if I call the code directly, while I've found the projUV datatype - I can't find the corresponding UTM version.


Could someone please get me started with the format of the string(s) of either method required?



Perhaps using conversion from zone 56 499451mE 6960921mS to lat/long as an example.



Answer



You may want to take a look at the route-me project, they seem to have some working examples of wrapping Proj.4 on the iOS platform in Objective-C. The actual transformation looks to happen in Mapview/Map/RMProjection.m.


No comments:

Post a Comment

arcpy - Changing output name when exporting data driven pages to JPG?

Is there a way to save the output JPG, changing the output file name to the page name, instead of page number? I mean changing the script fo...