Thursday 21 July 2016

postgresql - Create custom coordinate system in postgis


I have a Wayne county Michigan parcel dataset with a custom SRID: when I bring it into ArcGIS


Projected Coordinate System:    NAD_1983_HARN_StatePlane_Michigan_South_FIPS_2113_Feet_Intl
Projection: Lambert_Conformal_Conic

False_Easting: 13123359.58005249
False_Northing: 0.00000000
Central_Meridian: -84.36666667
Standard_Parallel_1: 42.10000000
Standard_Parallel_2: 43.66666667
Latitude_Of_Origin: 41.50000000
Linear Unit: Foot

Geographic Coordinate System: GCS_North_American_1983_HARN
Datum: D_North_American_1983_HARN

Prime Meridian: Greenwich
Angular Unit: Degree

I am trying to insert this SRID into postgis. Is this enough info to go and create the insert statement into postgreSQL?



Answer



It's easy here



  1. Look it up on http://www.spatialreference.org

  2. Find it SR-ORG:7069 NAD_1983_HARN_StatePlane_Michigan_South_FIPS_2113_IntlFeet

  3. Click PostGIS spatial_ref_sys INSERT statement


  4. Run that insert command.


Reproduced below,


INSERT into spatial_ref_sys (srid, auth_name, auth_srid, proj4text, srtext)
values ( 97069, 'sr-org', 7069, '', 'PROJCS["NAD_1983_HARN_StatePlane_Michigan_South_FIPS_2113_IntlFeet",GEOGCS["GCS_North_American_1983_HARN",DATUM["D_North_American_1983_HARN",SPHEROID["GRS_1980",6378137.0,298.257222101]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]],PROJECTION["Lambert_Conformal_Conic"],PARAMETER["False_Easting",13123359.58005249],PARAMETER["False_Northing",0.0],PARAMETER["Central_Meridian",-84.36666666666666],PARAMETER["Standard_Parallel_1",42.1],PARAMETER["Standard_Parallel_2",43.66666666666666],PARAMETER["Latitude_Of_Origin",41.5],UNIT["Foot",0.3048]]');

Now you have SRID 7069.


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...