Point 2 Grid

since version 7.015

Since GeoDMS 7.015 a relation to a grid-domain can be configured easily with the value function (example), based on a point attribute of a vector domain.

The resulting attribute will contain the index-numbers of the related grid domain for the vector domain.

example

// coordinate system 
unit<fpoint> rdc_base : format = "EPSG:28992",
unit<fpoint> rdc      := range(rdc_base, point(300000f,0f), point(625000f,280000f));

// configuration of grid domain, related to the rdc
unit<spoint> griddomain := range(
      gridset(rdc, point(-100f, 100f, rdc), point(625000f, 10000f, rdc) ,'spoint')
     ,point(   0s,    0s)
     ,point(3250s, 2700s)
);

 // configuration of vector domain, with a point attribute in rdc coordinates
unit<uint32> vectordomain: nrofrows = 4
{
   attribute<rdc> point: [{20427,69272},{17502,95885},{3188,80531},{12620,112190}];
}

attribute<griddomain> grid_7015_and_later_rel (vectordomain) := vectordomain/point[griddomain];

A grid and vector domain_unit are configured, both with rdc as coordinate unit. In the bold line the actual relation is configured. The value functions works as the coordinates used are expressed in the same coordinates (in this case rdc).

The value function can also be applied to calculate coordinate transformations, if the format strings of the coordinate systems are known. These strings need to be configured in the format property of these base units for the coordinate systems.

before version 7.015

Before GeoDMS version 7.015 such a relation needed to be configured with the getprojectionfactor and getprojectionoffset functions to relate a point attribute in world coordinates to a grid domain.

example

// coordinate system, grid and vector domain are similar 
parameter<rdc> projOffset := GetProjectionOffset (griddomain); // Result: [{625000, 10000)]
parameter<rdc> projFactor := GetProjectionFactor (griddomain); // Result: [(-100.0, 100.0)]

attribute<griddomain> grid_before_7015_rel (vectordomain) 
   := value((vectordomain/point - projOffset) / projFactor, griddomain);

In the bold line the actual relation is configured. The projOffset and projFactor are used to relate the points in the rdc coordinates to the griddomain.

backward incompability

To make your configuration working both in versions before and after 7.015, you need to configure an expression in which the geodmsversion function is used to determine the running GeoDMS version.

example

attribute<griddomain> grid_rel (vectordomain) := = GeoDMSVersion() > 7.015
   ? 'vectordomain/point[rdc]'
   : 'value((vectordomain/point - projOffset) / projFactor, griddomain")';

The expression in the example with the double == is an indirect-expression. Based on the GeoDMSVersion the grid_rel is calculated in a different way.