[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