Shapefiles with EPSG 2926 and 2927 Washington State Plane North and South not recognized

[qgib]

Author Name: clifford snow (clifford snow)
Original Redmine Issue: 13835
Affected QGIS version: 2.12.0
Redmine category:projection_support

---

When opening shapefiles from Washington county and state agencies, QGIS defaults to USER:10000X projection instead of either EPSG 2926 or 2927. Both of the projections are in QGIS.

This doesn't impact QGIS use, but as someone new to QGIS and projections, it would help if QGIS recognized and displayed the projection.

Attached is a copy of Pierce County GIS Roads .prj shapefile. Full shapefiles are available at http://gisdata.piercecowa.opendata.arcgis.com/datasets?q=Roads%2FRails

---

• Roads-Named.prj (clifford snow)

[qgib]

Author Name: Giovanni Manghi (@gioman)

---

2927 is

+proj=lcc +lat_1=47.33333333333334 +lat_2=45.83333333333334 +lat_0=45.33333333333334 +lon_0=-120.5 +x_0=500000.0001016001 +y_0=0 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=us-ft +no_defs

this shapefiles have

+proj=lcc +lat_1=47.33333333333334 +lat_2=45.83333333333334 +lat_0=45.33333333333334 +lon_0=-120.5 +x_0=500000.0001016002 +y_0=0 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=us-ft +no_defs

there is one tiny difference, but that is enough to make qgis think is a different CRS.

The definition on spatialreference.org is

http://spatialreference.org/ref/epsg/2927/proj4/

... +x_0=500000.0001016001 ...

---

• operating_system was changed from mac os x and ubuntu to
• category_id was configured as Projection Support
• status_id was changed from Open to Feedback
• os_version was changed from 10.10.5 and 15.04 to

[qgib]

Author Name: clifford snow (clifford snow)

---

Changing the longitude to match the .prj file (+x_0=500000.0001016001 to +x_0=500000.0001016002) made no difference.) I also removed the system generated projections, 100001 with no change. QGIS still did not match the projection.

[qgib]

Author Name: Jürgen Fischer (@jef-n)

---

• subject was changed from Shapefiles with EPSG 2926 and 2927 Washington State Plane North and South not recongized to Shapefiles with EPSG 2926 and 2927 Washington State Plane North and South not recognized

[qgib]

Author Name: Giovanni Manghi (@gioman)

---

does not seems a qgis issue to me, but rather a data one. This is what OGR says about the projection of the downloaded data, as you can see there is no "authority" parameter

giovanni@sibirica:~/Desktop/Railroads > ogrinfo -so Railroads.shp Railroads
INFO: Open of Railroads.shp' using driver ESRI Shapefile' successful.

Layer name: Railroads
Geometry: Line String
Feature Count: 163
Extent: (1093512.041297, 512003.733156) - (1230069.994152, 722811.599874)
Layer SRS WKT:
PROJCS["NAD83_HARN_Washington_South_ftUS",
GEOGCS["GCS_NAD83(HARN)",
DATUM["NAD83_High_Accuracy_Reference_Network",
SPHEROID["GRS_1980",6378137,298.257222101]],
PRIMEM["Greenwich",0],
UNIT["Degree",0.017453292519943295]],
PROJECTION["Lambert_Conformal_Conic_2SP"],
PARAMETER["standard_parallel_1",47.33333333333334],
PARAMETER["standard_parallel_2",45.83333333333334],
PARAMETER["latitude_of_origin",45.33333333333334],
PARAMETER["central_meridian",-120.5],
PARAMETER["false_easting",1640416.667],
PARAMETER["false_northing",0],
UNIT["Foot_US",0.30480060960121924]]
OBJECTID: Integer (10.0)
PCID: String (80.0)
RouteID: String (80.0)
Fm: String (80.0)
Tm: String (80.0)
Source: String (80.0)
OwnerName: String (80.0)
OwnerCode: String (80.0)
EditBy: String (80.0)
EditOn: String (80.0)
Type: String (80.0)
DataSource: String (80.0)

After giving the layer the proper projection 2927 the paramter is correctly recognized

giovanni@sibirica:~/Desktop/Railroads > ogrinfo -so Railroads.shp Railroads
INFO: Open of Railroads.shp' using driver ESRI Shapefile' successful.

Layer name: Railroads
Geometry: Line String
Feature Count: 163
Extent: (1093512.041297, 512003.733156) - (1230069.994152, 722811.599874)
Layer SRS WKT:
PROJCS["NAD83(HARN) / Washington South (ftUS)",
GEOGCS["NAD83(HARN)",
DATUM["NAD83_High_Accuracy_Reference_Network",
SPHEROID["GRS 1980",6378137,298.257222101,
AUTHORITY["EPSG","7019"]],
TOWGS84[0,0,0,0,0,0,0],
AUTHORITY["EPSG","6152"]],
PRIMEM["Greenwich",0,
AUTHORITY["EPSG","8901"]],
UNIT["degree",0.0174532925199433,
AUTHORITY["EPSG","9122"]],
AUTHORITY["EPSG","4152"]],
PROJECTION["Lambert_Conformal_Conic_2SP"],
PARAMETER["standard_parallel_1",47.33333333333334],
PARAMETER["standard_parallel_2",45.83333333333334],
PARAMETER["latitude_of_origin",45.33333333333334],
PARAMETER["central_meridian",-120.5],
PARAMETER["false_easting",1640416.667],
PARAMETER["false_northing",0],
UNIT["US survey foot",0.3048006096012192,
AUTHORITY["EPSG","9003"]],
AXIS["X",EAST],
AXIS["Y",NORTH],
AUTHORITY["EPSG","2927"]]
OBJECTID: Integer (10.0)
PCID: String (80.0)
RouteID: String (80.0)
Fm: String (80.0)
Tm: String (80.0)
Source: String (80.0)
OwnerName: String (80.0)
OwnerCode: String (80.0)
EditBy: String (80.0)
EditOn: String (80.0)
Type: String (80.0)
DataSource: String (80.0)

[qgib]

Author Name: clifford snow (clifford snow)

---

Typically I don't see authority parameters in .prj files. For example below is an ogrinfo report for a Census roads file for Franklin County, Washington, pulled this morning, that QGIS reports as EGPS:4269

ogrinfo -so tl_2015_53021_roads.shp tl_2015_53021_roads
INFO: Open of tl_2015_53021_roads.shp' using driver ESRI Shapefile' successful.

Layer name: tl_2015_53021_roads
Geometry: Line String
Feature Count: 2944
Extent: (-119.453704, 46.202270) - (-118.203138, 46.738730)
Layer SRS WKT:
GEOGCS["GCS_North_American_1983",
DATUM["North_American_Datum_1983",
SPHEROID["GRS_1980",6378137,298.257222101]],
PRIMEM["Greenwich",0],
UNIT["Degree",0.017453292519943295]]
LINEARID: String (22.0)
FULLNAME: String (100.0)
RTTYP: String (1.0)
MTFCC: String (5.0)

[qgib]

Author Name: Giovanni Manghi (@gioman)

---

clifford snow wrote:

Typically I don't see authority parameters in .prj files. For example below is an ogrinfo report for a Census roads file for Franklin County, Washington, pulled this morning, that QGIS reports as EGPS:4269

right, but anyway the rationale does not change: QGIS internal CRS DB is made of standard EPSG definition, if a layer has a CRS that does not match (exactly) any of them then QGIS adds the layer with a custom CRS, that still works. You can always (easily) force QGIS to recreate the prj file and give such layers the correct EPSG definition for the CRS. I can't see any QGIS bug here, but of course if I'm wrong please reopen this ticket.

---

• status_id was changed from Feedback to Closed
• resolution was changed from to invalid

[qgib]

Author Name: clifford snow (clifford snow)

---

I'd like to reopen this after discovering prj2epsg.org run by Boundless. When I plug in the .prj (below) the site returns 2285. (I was surprised.) I believe QGIS should return the same.

---

• status_id was changed from Closed to Reopened

[qgib]

Author Name: Giovanni Manghi (@gioman)

---

clifford snow wrote:

I'd like to reopen this after discovering prj2epsg.org run by Boundless. When I plug in the .prj (below) the site returns 2285. (I was surprised.) I believe QGIS should return the same.

[...]

can you attach the data in 2285?

---

• status_id was changed from Reopened to Feedback
• resolution was changed from invalid to

[qgib]

Author Name: clifford snow (clifford snow)

---

Attached is the proj file from Skagit County Roads-Named.shp.

---

• 9459 was configured as Roads-Named.prj

---

• Roads-Named.prj (clifford snow)

[qgib]

Author Name: Josiah Wahlrab (Josiah Wahlrab)

---

I believe this comes down to how ESRI and QGIS handle projections differently. A similar situation is described in this "blog post":https://www.northrivergeographic.com/archives/qgis-2-x-you-can-run-but-you-cant-hide-from-projection-issues. I don't believe this to be a QGIS bug, but I do think it would be super useful for QGIS to be able to discern projections like "prj2epsg":http://prj2epsg.org/.

[qgib]

Author Name: Giovanni Manghi (@gioman)

---

closed in favor of #22824

---

• status_id was changed from Feedback to Closed