1

I have two files I am using, one is census tract data with prj:

GEOGCS["GCS_North_American_1983",
    DATUM["D_North_American_1983",
         SPHEROID["GRS_1980",6378137.0,298.257222101]],
    PRIMEM["Greenwich",0.0],
    UNIT["Degree",0.0174532925199433]]

The other is a grid file from a model with prj:

PROJCS["airpact",
    GEOGCS["GCS_Sphere_EMEP",
       DATUM["D_Sphere_EMEP",
         SPHEROID["Sphere_EMEP",6370000.0,0.0]],
       PRIMEM["Greenwich",0.0],
       UNIT["Degree",0.0174532925199433]],
    PROJECTION["Lambert_Conformal_Conic"],
    PARAMETER["False_Easting",0.0],
    PARAMETER["False_Northing",0.0],
    PARAMETER["Central_Meridian",-121.0],
    PARAMETER["Standard_Parallel_1",30.0],
    PARAMETER["Standard_Parallel_2",60.0],
    PARAMETER["Scale_Factor",1.0],
    PARAMETER["Latitude_Of_Origin",49.0],
    UNIT["Meter",1.0]]

Data file In layer spatial reference system units

xMin,yMin -124.849,41.988 : xMax,yMax -111.043,49.0025

Layer Spatial Reference System

+proj=longlat +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +no_defs

Grid file In layer spatial reference system units

xMin,yMin -342000.00,-942000.00 : xMax,yMax 798000.00,90000

Layer Spatial Reference System

+proj=lcc +lat_1=30 +lat_2=60 +lat_0=49 +lon_0=-121 +x_0=0 +y_0=0 +a=6370000 +b=6370000 +units=m +no_defs

My problem is obvious, these are nowhere near the same and when I load into qgis the grid is far larger than the census tract data. I need to define them both on the same projection as the grid without ruining the census data but I am not sure how to do this. I have read several similar questions but am still lost.

I am aware that the census data is referencing lat/lon in degrees and the grid is in projected coordinates in meters, but need to somehow convert the census data to match the grid.

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Browse other questions tagged or ask your own question.