verify the projections of all cellular spaces from data scripts of terralib package

[pedro-andrade-inpe]

verify the projections of all cellular spaces created from data scripts of terralib package. The created projections are different from the projection of the original data. For example:

emas.prj:

PROJCS["UTM_Zone_22_Southern_Hemisphere",GEOGCS["GCS_Australian Natl & S. Amer. 1969",DATUM["D_unknown",SPHEROID["aust_SA",6378160,298.25]],PRIMEM["Greenwich",0],UNIT["Degree",0.017453292519943295]],PROJECTION["Transverse_Mercator"],PARAMETER["latitude_of_origin",0],PARAMETER["central_meridian",-51],PARAMETER["scale_factor",0.9996],PARAMETER["false_easting",500000],PARAMETER["false_northing",10000000],UNIT["Meter",1]]

emas-limit.prj

PROJCS["SAD69 / UTM zone 22S", 
  GEOGCS["SAD69", 
    DATUM["South American Datum 1969", 
      SPHEROID["GRS 1967 Modified", 6378160.0, 298.25, AUTHORITY["EPSG","7050"]], 
      TOWGS84[-66.87, 4.37, -38.52, 0.0, 0.0, 0.0, 0.0], 
      AUTHORITY["EPSG","6618"]], 
    PRIMEM["Greenwich", 0.0, AUTHORITY["EPSG","8901"]], 
    UNIT["degree", 0.017453292519943295], 
    AXIS["Geodetic longitude", EAST], 
    AXIS["Geodetic latitude", NORTH], 
    AUTHORITY["EPSG","4618"]], 
  PROJECTION["Transverse Mercator", AUTHORITY["EPSG","9807"]], 
  PARAMETER["central_meridian", -51.0], 
  PARAMETER["latitude_of_origin", 0.0], 
  PARAMETER["scale_factor", 0.9996], 
  PARAMETER["false_easting", 500000.0], 
  PARAMETER["false_northing", 10000000.0], 
  UNIT["m", 1.0], 
  AXIS["Easting", EAST], 
  AXIS["Northing", NORTH], 
  AUTHORITY["EPSG","29192"]]

Try to update the projection of the emas-limit with the one with the one of the cells and see if it creates the same output.

[pedro-andrade-inpe]

Look how amazonia data is shifted from the real location: