83 lines
3.0 KiB
MySQL
83 lines
3.0 KiB
MySQL
|
-- Internal function.
|
|||
|
|
-- Set the seeds of the RNGs (Random Number Generators)
|
|||
|
|
-- used internally.
|
|||
|
|
CREATE OR REPLACE FUNCTION
|
|||
|
|
_cdb_random_seeds (seed_value INTEGER) RETURNS VOID
|
|||
|
|
AS $$
|
|||
|
|
from crankshaft import random_seeds
|
|||
|
|
random_seeds.set_random_seeds(seed_value)
|
|||
|
|
$$ LANGUAGE plpythonu;
|
|||
|
|
-- Moran's I
|
|||
|
|
CREATE OR REPLACE FUNCTION
|
|||
|
|
cdb_moran_local (
|
|||
|
|
t TEXT,
|
|||
|
|
attr TEXT,
|
|||
|
|
significance float DEFAULT 0.05,
|
|||
|
|
num_ngbrs INT DEFAULT 5,
|
|||
|
|
permutations INT DEFAULT 99,
|
|||
|
|
geom_column TEXT DEFAULT 'the_geom',
|
|||
|
|
id_col TEXT DEFAULT 'cartodb_id',
|
|||
|
|
w_type TEXT DEFAULT 'knn')
|
|||
|
|
RETURNS TABLE (moran FLOAT, quads TEXT, significance FLOAT, ids INT)
|
|||
|
|
AS $$
|
|||
|
|
from crankshaft.clustering import moran_local
|
|||
|
|
# TODO: use named parameters or a dictionary
|
|||
|
|
return moran_local(t, attr, significance, num_ngbrs, permutations, geom_column, id_col, w_type)
|
|||
|
|
$$ LANGUAGE plpythonu;
|
|||
|
|
|
|||
|
|
-- Moran's I Local Rate
|
|||
|
|
CREATE OR REPLACE FUNCTION
|
|||
|
|
cdb_moran_local_rate(t TEXT,
|
|||
|
|
numerator TEXT,
|
|||
|
|
denominator TEXT,
|
|||
|
|
significance FLOAT DEFAULT 0.05,
|
|||
|
|
num_ngbrs INT DEFAULT 5,
|
|||
|
|
permutations INT DEFAULT 99,
|
|||
|
|
geom_column TEXT DEFAULT 'the_geom',
|
|||
|
|
id_col TEXT DEFAULT 'cartodb_id',
|
|||
|
|
w_type TEXT DEFAULT 'knn')
|
|||
|
|
RETURNS TABLE(moran FLOAT, quads TEXT, significance FLOAT, ids INT, y numeric)
|
|||
|
|
AS $$
|
|||
|
|
from crankshaft.clustering import moran_local_rate
|
|||
|
|
# TODO: use named parameters or a dictionary
|
|||
|
|
return moran_local_rate(t, numerator, denominator, significance, num_ngbrs, permutations, geom_column, id_col, w_type)
|
|||
|
|
$$ LANGUAGE plpythonu;
|
|||
|
|
-- Function by Stuart Lynn for a simple interpolation of a value
|
|||
|
|
-- from a polygon table over an arbitrary polygon
|
|||
|
|
-- (weighted by the area proportion overlapped)
|
|||
|
|
-- Aereal weighting is a very simple form of aereal interpolation.
|
|||
|
|
--
|
|||
|
|
-- Parameters:
|
|||
|
|
-- * geom a Polygon geometry which defines the area where a value will be
|
|||
|
|
-- estimated as the area-weighted sum of a given table/column
|
|||
|
|
-- * target_table_name table name of the table that provides the values
|
|||
|
|
-- * target_column column name of the column that provides the values
|
|||
|
|
-- * schema_name optional parameter to defina the schema the target table
|
|||
|
|
-- belongs to, which is necessary if its not in the search_path.
|
|||
|
|
-- Note that target_table_name should never include the schema in it.
|
|||
|
|
-- Return value:
|
|||
|
|
-- Aereal-weighted interpolation of the column values over the geometry
|
|||
|
|
CREATE OR REPLACE
|
|||
|
|
FUNCTION cdb_overlap_sum(geom geometry, target_table_name text, target_column text, schema_name text DEFAULT NULL)
|
|||
|
|
RETURNS numeric AS
|
|||
|
|
$$
|
|||
|
|
DECLARE
|
|||
|
|
result numeric;
|
|||
|
|
qualified_name text;
|
|||
|
|
BEGIN
|
|||
|
|
IF schema_name IS NULL THEN
|
|||
|
|
qualified_name := Format('%I', target_table_name);
|
|||
|
|
ELSE
|
|||
|
|
qualified_name := Format('%I.%s', schema_name, target_table_name);
|
|||
|
|
END IF;
|
|||
|
|
EXECUTE Format('
|
|||
|
|
SELECT sum(%I*ST_Area(St_Intersection($1, a.the_geom))/ST_Area(a.the_geom))
|
|||
|
|
FROM %s AS a
|
|||
|
|
WHERE $1 && a.the_geom
|
|||
|
|
', target_column, qualified_name)
|
|||
|
|
USING geom
|
|||
|
|
INTO result;
|
|||
|
|
RETURN result;
|
|||
|
|
END;
|
|||
|
|
$$ LANGUAGE plpgsql;
|