Calculate cbs vierkanstatistieken cell id

calculate_cbs_vierkanstatistieken_cell_id.sql

-- function that calculates the CBS Vierkantstatistieken grid cell for point coordinates in RD new.

create or replace function get_cbsvierkantid(fn_geom geometry, i_dist int) returns text as $$
select 
	'E' || 
	((st_x(fn_geom)::numeric - (st_x(fn_geom)::numeric % i_dist)) / 100)::int ||
	'N' || 
	((st_y(fn_geom)::numeric - (st_y(fn_geom)::numeric % i_dist)) / 100)::int
	as cbs_vierkant_id
;
$$ language sql;

-- example for 500x500 grid

select get_cbsvierkantid(st_geomfromtext('POINT(233544.795083973 581928.166658598)'),500);