search

ObjectVision / GLOBIO_dms

0 stars 0 forks source link

LScale -> spatial aggregation claims do not match with available landcells #1

Closed eoudejans closed 7 months ago

eoudejans commented 7 months ago

The study area needs to be exactly at block boundaries of the to be scaled blocks or else suitability map and claim determination will not work properly.

Adding a small offset resulted in proper alignment of blocks.

container study_areas
{
    parameter<float64> offset05secs :=  5.0 / 3600.0;

    parameter<LatLong> lower_bound := LowerBound(brondata/GlobCoverTIFF)[LatLong];
    //parameter<float64> offset := Iceland lower_bound  point_yx(180.0, 180.0, );

    unit<ipoint> Europe  := range(brondata/GlobCoverTIFF
    , point_xy(-12.0+offset05secs, 30.0-offset05secs, geometry/LatLong)[brondata/GlobCoverTIFF]
    , point_xy(45.0+offset05secs, 70.0-offset05secs, geometry/LatLong)[brondata/GlobCoverTIFF]);

    unit<ipoint> Iceland := range(brondata/GlobCoverTIFF
    , point_xy(-30.0+offset05secs, 60.0-offset05secs, geometry/LatLong)[brondata/GlobCoverTIFF]
    , point_xy(  0.0+offset05secs, 75.0-offset05secs, geometry/LatLong)[brondata/GlobCoverTIFF]);
}

Resulting in for example proper allocation for 180x for items: /results/alloc_raw_world/x180/LULC_rel (right) and /results/majority_modus/x180/LULC_rel (left)

image