PostGis basico By @geojuanka

Crear campos nuevos:

alter table esquema.tabla add nombrecampo tipodecampo

Ejemplo: alter table esquema01.municipios add perimetro numeric (20,2);

Borrar un campo:

alter table esquema.tabla drop column nombrecampo

Ejemplo: alter table esquema01.municipios_line drop column perimetro;

Actualizar valores de un campo:

update esquema.tabla set nombrecampo = valores

Ejemplo: update esquema01.municipios_pto SET tipo = 'residencial';

Actualizar valores de un campo con condiciones:

update esquema.tabla set nombrecampo = valores where nombrecampo = valores

Ejemplo: update esquema01.municipios_pto SET tipo = 'capital' where municipio = 'LAS PALMAS DE GRAN CANARIA';

Calcular el área (ST_AREA)

update esquema.tabla set nombrecampo = st_area (campodelageometría)

Ejemplo: update esquema01.municipios SET area = st_area(geom);

Calcular el perímetro (ST_PERIMETER)

update esquema.tabla set nombrecampo = st_perimeter (campodelageometría)

Ejemplo: update esquema01.municipios SET perimetro = st_perimeter(geom);

Calcular la longitud (ST_LENGTH)

update esquema.tabla set nombrecampo = st_length (campodelageometría)

Ejemplo: update esquema01.municipios_line SET longitud = st_length(geom);

Calcular la coordenada X (ST_X)

update esquema.tabla set nombrecampo = st_x (campodelageometría)

Ejemplo: update esquema01.municipios_pto SET coord_x = st_x(geom);

Convertir una capa de polígonos a líneas (ST_BOUNDARY)

create table esquema.nuevacapadelineas as select st_boundary (campodelageometría), campo1, campo2, campo3 from esquema.capadepoligonos

Ejemplo: create table esquema01.municipios_line as select ST_boundary(geom), codmun, municipio, isla from esquema01.municipios;

Convertir una capa de polígonos a puntos (ST_CENTROID)

create table esquema.nuevacapadepuntos as select st_centroid (campodelageometría), campo1, campo2, campo3 from esquema.capadepoligonos

Ejemplo: create table esquema01.municipios_pto as select ST_centroid(geom), codmun, municipio, isla from esquema01.municipios;

Convertir una capa de polígonos a puntos (centroides inside) (ST_POINTONSURFACE)

create table esquema.nuevacapadepuntos as select ST_PointOnSurface (campodelageometría), campo1, campo2, campo3 from esquema.capadepoligonos

Ejemplo: create table esquema01.municipios_pto_inside as select ST_PointOnSurface (geom), codmun, municipio, isla from esquema01.municipios;

Añadir antiguedades del omisor a las construcciones

Creamos un CSV desde el omisor con lo campos municipio, referencia, anyo Añadimos columna anyo en la tabla del municipio en PostGis Actualizamos campo

copy antiguedad (municipio, referencia, anyo) from 'C:\a\CAT\municipio.csv' delimiter ';' csv header;

ALTER TABLE municipio ADD COLUMN anyo integer constraint;

UPDATE municipio SET anyo= antiguedad.anyo
FROM antiguedad
WHERE antiguedad.referencia=cox.refcat;

Crear tabla incidencias

CREATE TABLE incidencias (
  id serial PRIMARY KEY,
  municipio character varying(3),
  comentario character varying(250),
  geom geometry(Point, 25830),
  referencia character varying(14),
  fecha timestamp,
  usuario character varying(100)
);

Trigger para añadir la referencia catastral con la que intersecta

Si no intersecta con nada se añade valor "00000000000000"

CREATE OR REPLACE FUNCTION actualizar_referencia()
  RETURNS TRIGGER AS $$
BEGIN
  NEW.referencia := COALESCE((SELECT refcat FROM parcelas_torrevieja WHERE ST_Intersects(NEW.geom, geom) LIMIT 1), '00000000000000');
  RETURN NEW;
END;
$$ LANGUAGE plpgsql;

CREATE TRIGGER incidencias_referencia_trigger
  BEFORE INSERT ON incidencias
  FOR EACH ROW
  EXECUTE PROCEDURE actualizar_referencia();

Trigger para actualizar la fecha

CREATE OR REPLACE FUNCTION actualizar_fecha()
  RETURNS TRIGGER AS $$
BEGIN
  NEW.fecha := now();
  RETURN NEW;
END;
$$ LANGUAGE plpgsql;

CREATE TRIGGER incidencias_fecha_trigger
  BEFORE INSERT ON incidencias
  FOR EACH ROW
  EXECUTE PROCEDURE actualizar_fecha();