Question

Projection Issue with geotiffs being written


Badge +1

I'm currently using reading geotiffs(mosaics) generated from pix4d. And writing them backout again by removing the black band that is generated around them. 

 

When I started publishing these out to arcgis image server I started to notice the projection was not aligning with points of the projected that have been collected. But the original generated by pix4d didn't have this alignment problem.

 

All our mosaics created since 2021 are projected out in EPSG:6549 and when the geotiff is read through the feature reader it looks it's being detected as the following:

 

Coordinate System Parameters
 
CS_NAME: NSRS11.OH-NF
DESC_NM: NSRS 2011 Ohio State Planes, North Zone, US Foot
DT_NAME: NSRS11
EPSG_QD: 1
GROUP: SPNSRS11F
MAP_SCL: 1
MAX_LAT: 42.55
MAX_LNG: -79.98333333333333
MIN_LAT: 39.88333333333333
MIN_LNG: -85.31666666666666
ORG_LAT: 39.666666666666664
ORG_LNG: -82.5
PARM1: 41.7
PARM2: 40.43333333333333
PROJ: LM
QUAD: 1
SCL_RED: 1
SOURCE: CRS HARN/OH.OH-NF with datum modified to NSRS11
UNIT: FOOT
X_OFF: 1968500
ZERO_X: 0.001
ZERO_Y: 0.001
 
 
Datum Parameters
 
DESC_NM: NSRS 2011: Alaska, Continental US, PR/VI
ELLIPSOID: GRS1980
SOURCE: National Geodetic Survey (ngs.noaa.gov/NationalReadjustment)
USE: NAD83
 
 
Ellipsoid Parameters
 
DESC_NM: Geodetic Reference System of 1980
E_RAD: 6378137
P_RAD: 6356752.314140348
SOURCE: Stem, L.E., Jan 1989, State Plane Coordinate System of 1983
 
 
OGC WKT Description
 
PROJCS["NSRS 2011 Ohio State Planes, North Zone, US Foot",
    GEOGCS["NAD83_National_Spatial_Reference_System_2011",
        DATUM["NAD83_National_Spatial_Reference_System_2011",
            SPHEROID["Geodetic Reference System of 1980",6378137,298.257222100892],
            AUTHORITY["EPSG","1116"]],
        PRIMEM["Greenwich",0],
        UNIT["degree",0.0174532925199433,
            AUTHORITY["EPSG","9122"]]],
    PROJECTION["Lambert_Conformal_Conic_2SP"],
    PARAMETER["latitude_of_origin",39.6666666666667],
    PARAMETER["central_meridian",-82.5],
    PARAMETER["standard_parallel_1",41.7],
    PARAMETER["standard_parallel_2",40.4333333333333],
    PARAMETER["false_easting",1968500],
    PARAMETER["false_northing",0],
    UNIT["Foot_US",0.304800609601219],
    AXIS["Easting",EAST],
    AXIS["Northing",NORTH]]
 
 
 
Esri WKT Description
 
PROJCS["NAD_1983_2011_StatePlane_Ohio_North_FIPS_3401_Ft_US",GEOGCS["GCS_NAD_1983_2011",DATUM["D_NAD_1983_2011",SPHEROID["GRS_1980",6378137.0,298.257222101]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]],PROJECTION["Lambert_Conformal_Conic"],PARAMETER["False_Easting",1968500.0],PARAMETER["False_Northing",0.0],PARAMETER["Central_Meridian",-82.5],PARAMETER["Standard_Parallel_1",40.43333333333333],PARAMETER["Standard_Parallel_2",41.7],PARAMETER["Latitude_Of_Origin",39.66666666666666],UNIT["Foot_US",0.3048006096012192],AUTHORITY["EPSG",6549]]

I'm not sure I under/following this projection because after doing some comparisons I found a few discrepancies

 

This is the OGC WKT for EPSG:6549 and if you compare it to what I posted above they do not match. And I'm not sure why?

PROJCS["unnamed",GEOGCS["GRS 1980(IUGG, 1980)",DATUM["unknown",SPHEROID["GRS80",6378137,298.257222101]],PRIMEM["Greenwich",0],UNIT["degree",0.0174532925199433]],PROJECTION["Lambert_Conformal_Conic_2SP"],PARAMETER["standard_parallel_1",41.7],PARAMETER["standard_parallel_2",40.43333333333333],PARAMETER["latitude_of_origin",39.66666666666666],PARAMETER["central_meridian",-82.5],PARAMETER["false_easting",1968500],PARAMETER["false_northing",0],UNIT["Foot_US",0.3048006096012192],AUTHORITY["epsg","6549"]]

 The ESRI WKT matches

 

PROJCS["Lambert_Conformal_Conic",GEOGCS["GCS_GRS 1980(IUGG, 1980)",DATUM["D_unknown",SPHEROID["GRS80",6378137,298.257222101]],PRIMEM["Greenwich",0],UNIT["Degree",0.017453292519943295]],PROJECTION["Lambert_Conformal_Conic"],PARAMETER["standard_parallel_1",41.7],PARAMETER["standard_parallel_2",40.43333333333333],PARAMETER["latitude_of_origin",39.66666666666666],PARAMETER["central_meridian",-82.5],PARAMETER["false_easting",1968500],PARAMETER["false_northing",0],UNIT["Foot_US",0.30480060960121924]]

 

I'm not sure what the coordinate system that's being autodetected by FME is suppose to be? It doesn't appear to be correct or I'm misinterpreting the description. 


2 replies

Userlevel 2
Badge +17

Hi @kschultz​,

FME's GeoTIFF writer has the parameter 'Write Esri-Compatible Coordinate System, which is set to No by default. Please try changing it to Yes, then see if the coordinate system is now maintained by the writer.

Badge +1

Hi @kschultz​,

FME's GeoTIFF writer has the parameter 'Write Esri-Compatible Coordinate System, which is set to No by default. Please try changing it to Yes, then see if the coordinate system is now maintained by the writer.

Correct me if I'm wrong but I think the reader is the problem. How will it write the proper projection(unskewed) without reading the correct projection intially?

Reply