I have downloaded parcel files for several counties in Washington State, and am attempting to overlay them on the OpenStreetMap Standard base layer (which is in WGS84 Pseudo Mercator). The majority of the county layers were in WGS84 and automatically loaded into the correct extent/location, but some (such as King County) were in NAD83 (ESPG: 4269) is placed over West Africa in the OSM map.
What do I need to do to get this King County layer to load in the correct place on OpenLayers map? Or alternatively, how do I convert the file to the correct CRS using ogr2ogr?
If I set the CRS for the layer to the correct settings NAD83 (EPSG:4269), the layer is not visible. (but doing zoom to layer still takes me to west Africa). If I change the CRS in layer properties to EPSG:3857, it becomes visible but is still in the wrong location.
I also tried saving the King County NAD83 layer as a new file with CRS set to WGS84 PseudoMercator, and using the extents laid out in the metadata files:
... but when I loaded it, I can't see the layer or make it visible no matter what CRS I choose for it. So then I gave up trying to do it in QGIS, and attempted to use some command line tools to try to convert/reproject the file.
Here is the output of
ogrinfo for the King County NAD83 shapefile:
ogrinfo -so -al parcel.shp INFO: Open of `parcel.shp' using driver `ESRI Shapefile' successful. Layer name: parcel Metadata: DBF_DATE_LAST_UPDATE=2018-04-06 Geometry: Polygon Feature Count: 617501 Extent: (1218256.223500, 30932.680000) - (1577865.300000, 287704.244000) Layer SRS WKT: PROJCS["NAD_1983_HARN_StatePlane_Washington_North_FIPS_4601_Feet", GEOGCS["GCS_North_American_1983_HARN", DATUM["NAD83_High_Accuracy_Reference_Network", SPHEROID["GRS_1980",6378137.0,298.257222101]], PRIMEM["Greenwich",0.0], UNIT["Degree",0.0174532925199433], AUTHORITY["EPSG","4269"]], PROJECTION["Lambert_Conformal_Conic_2SP"], PARAMETER["False_Easting",1640416.666666667], PARAMETER["False_Northing",0.0], PARAMETER["Central_Meridian",-120.8333333333333], PARAMETER["Standard_Parallel_1",47.5], PARAMETER["Standard_Parallel_2",48.73333333333333], PARAMETER["Latitude_Of_Origin",47.0], UNIT["Foot_US",0.3048006096012192]] MAJOR: String (6.0) MINOR: String (4.0) PIN: String (10.0) Shape_area: Real (19.11) Shape_len: Real (19.11)
I tried using ogr2ogr to convert it to WGS84, but got an error message:
ogr2ogr -f "ESRI Shapefile" -t_srs EPSG:3857 -s_srs EPSG:4269 'King County Parcels (WGS84, EPSG 3857).shp' parcel.shp ERROR 1: Failed to reproject feature 0 (geometry probably out of source or destination SRS). ERROR 1: Terminating translation prematurely after failed translation of layer parcel (use -skipfailures to skip errors)
... I'm at a complete loss for what to do. I've looked through tons of tutorials and GIS Stackexchange posts on projections in QGIS, and seem to be missing something, but can't figure out what.
What am I doing wrong here, and how can I make this layer work with OpenStreetMap, and the other WGS84 layers?