[GRASS-user] Projection conflict with new project
Rich Shepard
rshepard at appl-ecosys.com
Fri Sep 22 14:33:35 PDT 2017
I'm starting a new project with the location name 'Washington.' I need to
provide projection information for it although I want each dataset to be in
a separate mapset. Viewing the .prj file for one dataset
(ambient_water_quality) I see the shapefile's using
NAD_1983_HARN_StatePlane_Washington_South_FIPS_4602_Feet. When I tell the
location generator to use an EPSG and browse to NAD83 it provides the number
4602. But, that number is assigned elsewhere when I try to open the .shp
file:
ERROR: Projection of dataset does not appear to match current location.
GRASS LOCATION PROJ_INFO is:
name: Dominica 1945
ellps: clark80
proj: ll
towgs84: 725,685,536,0,0,0,0
no_defs: defined
Import dataset PROJ_INFO is:
name: NAD_1983_HARN_StatePlane_Washington_South_FIPS_4602_Feet
datum: nad83harn
ellps: grs80
proj: lcc
lat_1: 45.83333333333334
lat_2: 47.33333333333334
lat_0: 45.33333333333334
lon_0: -120.5
x_0: 500000.0000000001
y_0: 0
towgs84: 0,0,0,0,0,0,0
no_defs: defined
Another data set (in FileOpenGDB format) has this projection information:
Layer name: WashingtonNPA
Geometry: Multi Polygon
Feature Count: 57
Extent: (721307.250000, 81877.333333) - (2512692.666667, 1344683.500000)
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"]]
Multiple EPSG numbers and I'm not sure which to use on this dataset.
What I want is a Washington-Feet location with each topic as a separate map set.
What have I forgotten so that I'm getting this error?
Rich
More information about the grass-user
mailing list