Projection used by BDC

The Albers Equal-Area Conic is the reference map projection for the BDC tiling system, named BDC Grid. Like other normal conics, the Albers Equal-Area Conic projection has concentric arcs of circles for parallels and equally spaced radii as meridians. The parallels are not equally spaced, but they are farthest apart in the latitudes between the standard parallels and closer together to the north and south. The pole is not the center of the circles, but it is normally an arc itself. Scale along the parallels is too small between the standard parallels and too large beyond them. The scale along the meridians is just the opposite and, in fact, the scale factor along meridians is the reciprocal of the scale factor along the parallels to maintain equal area [18].

The Albers Equal-Area Conic projection parameters follow the cartographic standards defined by Instituto Brasileiro de Geografia e Estatística (IBGE), which is the official civilian mapping agency of Brazil [7]. SIRGAS-2000, the official planimetric datum of Brazil, is the reference geodetic system. BDC Grid uses parallel 12° south and meridian 54° west as origin. To generate BDC Grid two standard parallels at latitudes 2° south and 22° south are used. Offsets of 5,000,000 meters and 10,000,000 meters are respectively applied to X and Y coordinates at the origin. In short, the WKT for the BDC Grid can be expressed by:

PROJCS[“unknown”,GEOGCS[“unknown”,DATUM[“Unknown based on GRS80 ellipsoid”,SPHEROID[“GRS 1980”,6378137,298.257222101,AUTHORITY[“EPSG”,”7019”]]],PRIMEM[“Greenwich”,0,AUTHORITY[“EPSG”,”8901”]],UNIT[“degree”,0.0174532925199433,AUTHORITY[“EPSG”,”9122”]]],PROJECTION[“Albers_Conic_Equal_Area”],PARAMETER[“latitude_of_center”,-12],PARAMETER[“longitude_of_center”,-54],PARAMETER[“standard_parallel_1”,-2],PARAMETER[“standard_parallel_2”,-22],PARAMETER[“false_easting”,5000000],PARAMETER[“false_northing”,10000000],UNIT[“metre”,1,AUTHORITY[“EPSG”,”9001”]],AXIS[“Easting”,EAST],AXIS[“Northing”,NORTH]]

and the projstring:

crs: +proj=aea +lat_0=-12 +lon_0=-54 +lat_1=-2 +lat_2=-22 +x_0=5000000 +y_0 = 10000000 +ellps=GRS80 +units=m +no_defs.

The BDC Grid System includes three hierarchical Grids. The size of the grid tiles varies according to the grid version, see bellow more details. The rationale behind the BDC Grid tiling system follows the concept illustrated below in Figure 1.

BDC Grid seed

Figure 1 - Tile seeds for BDC Grid.

BDC Grid System (V2)

BDC Grid System (V2) includes three hierarchical Grids: 105600m x 105600m (Small, SM), 211200m x 211200m (Medium, MD) and 4224400m x 4224400m (Large, LG) tiles. Each product in BDC is generated using a specific grid, depending on the spatial resolution of the data and the output image size. Sentinel-2/MSI products are generated in the SM grid, Landsat-8/OLI and CBERS-4/MUX are generated in the MD grid and CBERS-4/WFI are generated in the LG grid.

The hierarchical relation between the different BDC Grids (V2) can be seen in Figure 2.

BDC Grid V2

Figure 2 - BDC Grid V2.

Grid Generation

Tile seeds are initially positioned in a way that the center a tile in the smaller grid is with the meridian 54° west and parallel 12° south. With this initial tile seed defined, other tiles of the same corresponding sizes, always defined by Albers Equal-Area Conic plane coordinates, are automatically generated to entirely cover the cartesian plane of the reference map projection system. For each tile size, upper boundaries perfectly match lower boundaries of tiles above and left boundaries perfectly match right boundaries of tiles on the left. The remaining hierarchical levels are generated using the first grid as a reference.

BDC Grid System (V1)

BDC Grid System (V1) includes three hierarchical Grids: 1.5°x1°, 3°x2°, and 6°x4° tiles. Each product in BDC is generated using a specific grid, depending on the spatial resolution of the data and the output image size. Sentinel-2/MSI products are generated in the 1.5°x1° (SM) grid, Landsat-8/OLI and CBERS-4/MUX products are generated in the 3°x2° (MD) grid and CBERS-4/WFI are generated in the 6°x4° (LG) grid.

The hierarchical relation between the different BDC Grids (V1) can be seen in Figure 3.

BDC Grid V1

Figure 3 - BDC Grid V1.

Grid Generation

Tile seeds are initially positioned in a way that their upper left corners coincide at the intersection of the equator with the meridian 54° west. With those initial tile seeds defined, other tiles of the same corresponding sizes, always defined by Albers Equal-Area Conic plane coordinates, are automatically generated to entirely cover the cartesian plane of the reference map projection system. For each tile size, upper boundaries perfectly match lower boundaries of tiles above and left boundaries perfectly match right boundaries of tiles on the left.