4

I am trying to reproject a Michigan state land shapefile which is the NAD83/ Michigan Oblique Mercator EPSG:3078 to the WPG84 EPSG:4326 projection. However when i try to do this my map gets squished down. I have tried the commands in "Reprojecting from MGRm to WGS84 using open-source tools" with no luck. I need to get these coordinates into a decimal degrees format.

0

3 Answers 3

5

Ever since I posted in "Reprojecting from MGRm to WGS84 using open-source tools" last year, I've discovered that none of the EPSG definitions that supposedly describe Michigan GeoRef work correctly. I suspect that it is a flaw in the GDAL/OGR tools, and this has been discussed repeatedly in the GDAL developers forum, but no fix has been made. This is especially problematical for open source users in Malaysia, where a rectified skew orthomorphic (RSO) projection system is used country-wide. (Michigan GeoRef is an instance of an RSO projection, as far as I understand.)

I've also discovered that using uDIG to reproject from GeoRef to whatever seems to work, but I can't figure out why. I'm under the impression that uDIG also uses the GDAL/OGR libraries for conversion, but I may be wrong. I've been meaning to send the uDIG developers a query on that, but haven't gotten around to it yet.

Anyway, take a look at a tutorial I developed for our local N. Michigan township users using uDIG to convert GeoRef to MI State Plane. You can follow the same procedure and see if it works for you in converting to WGS84. The tutorial is at: https://web.archive.org/web/20200719135341/http://www.townshipgis.com:80/resources/problems-and-workarounds/converting-michigan-georef-to-state-plane-using-open-source-tools

0
1

Reproject on the fly:

  • add your shape
  • check in the properties if is really in EPSG:3078
  • Open the project properties and set the project CRS to EPSG 4326 and enable the checkbox "enable on the fly CRS transformation"

Create a copy of a vector in a different CRS:

  • add your shape
  • check in the properties if is really in EPSG:3078
  • Click the right mouse button over the vector in the TOC and select "save as...", then make the necessary choices in particular "CRS" that you must set to EPSG 4326
2
  • I have tried both of the above methods several times and these do not work. As Howard pointed out there is a problem with converting from the MiGeoref format that skews or squishes the output. However I should say this is the way to convert between projections in most cases. Jan 16, 2012 at 23:40
  • I cant edit my previous comment for some reason but I did a quick compare to converting the data with ArcGis. It is slightly different but not until the 6th decimal place(using decimal degrees). I cant see any visible different when overlaying on google maps. Jan 17, 2012 at 0:21
0

You can define a custom CRS within QGIS using the following proj string for Michigan GeoRef:

PROJCS["NAD_1983_Michigan_GeoRef_Meters",GEOGCS["GCS_North_American_1983",DATUM["D_North_American_1983",SPHEROID["GRS_1980",6378137.0,298.257222101]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]],PROJECTION["Hotine_Oblique_Mercator_Azimuth_Natural_Origin"],PARAMETER["False_Easting",2546731.496],PARAMETER["False_Northing",-4354009.816],PARAMETER["Scale_Factor",0.9996],PARAMETER["Azimuth",337.255555555556],PARAMETER["Longitude_Of_Center",-86.0],PARAMETER["Latitude_Of_Center",45.30916666666666],UNIT["Meter",1.0]]

Source: https://pages.mtu.edu/~mdhyslop/gis/

New contributor
Alex Hirzel is a new contributor to this site. Take care in asking for clarification, commenting, and answering. Check out our Code of Conduct.

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service and acknowledge you have read our privacy policy.

Not the answer you're looking for? Browse other questions tagged or ask your own question.