NAD_1983_StatePlane_Massachusetts_Mainland_FIPS_2001 and UNIT["Meter",1.0]

Luego en nuestra base de datos consultamos los SRIDs candidatos con una búsqueda como:

SELECT
srid, srtext, proj4text
FROM
spatial_ref_sys
WHERE
srtext ILIKE ‘%Massachusetts%’

De la que escogeremos el mejor candidato, es decir, aquel que sea más cercano a la proyección desrita en el archivo .prj.

Fuente: Part 1: Getting Started With PostGIS: An almost Idiot’s Guide

La entrada original se encuentra bajo una licencia GNU Free Documentation License 1.2 http://www.gnu.org/copyleft/fdl.html, al igual que esta entrada.

Esto lo conseguimos mediante estos sqls:

Para modificar los valores a nivel geometry_columns:

SELECT updategeometrysrid('<tabla>', '<columna_geometria>', <nuevo_srid>);

Y este para modificar los valores internos:

SELECT setsrid(<columna_geometria>, <nuevo_srid>) FROM <tabla>;

Fuente: [postgis-users] changing srid

SELECT
updategeometrySRID('<tabla>', '<columna_geometria>', <srid>);

obtenemos el siguiente error:

ERROR: column not found in geometry_columns table
CONTEXT: SQL statement "SELECT UpdateGeometrySRID('','', $1 , $2 , $3 )"
PL/pgSQL function "updategeometrysrid" line 4 at SQL statement

Este error se debe a que no tenemos resgitrada la columna en la tabla geometry_columns.

En la entrada anterior he citado como podemos agregar manualmente una columna a geometry_columns.

Esto lo logramos con el siguiente sql:

INSERT INTO
	geometry_columns
	(
	f_table_catalog,
	f_table_schema,
	f_table_name,
	f_geometry_column,
	coord_dimension,
	srid,
	"type"
	)
SELECT
	'',
	'public',
	'<tabla>',
	'<columna_geom>',
	ST_CoordDim(<columna_geom>),
	ST_SRID(<columna_geom>),
	GeometryType(<columna_geom>)
FROM
	public.<tabla>
LIMIT
	1;

Fuente: “Manually Registering Geometry Columns in geometry_columns” en el manual de PostGIS

¿De que habló?

Al ejecutar el comando:

createdb -T template_postgis my_spatial_db

la base de datos no se crea y si usamos pgAdminIII nos devuelve el siguiente error:

database "template_postgis" is being accessed by other users

Después de googlear un rato y de buscar solucionar el comando se me ocurrió hacerlo vía sql y si funcionó.

Este es el comando:

CREATE DATABASE my_spatial_db WITH TEMPLATE=template_postgis;

Instalación de PostGIS

Una vez instalado el PostgreSQL instalamos el PostGIS.

Nos vamos a sentir tentados a instalar el paquete postgis, pero este es un metapaquete que por alguna razón no nos instala el PostGIS.

El que debemos instalar es postgresql-8.3-postgis.

Por consola el comando es:

sudo apt-get install postgresql-8.3-postgis

Con eso ya tenemos instalado el PostGIS.

Nos hace falta crear ahora el template_postgis.

Creación del template_postgis

Nos logeamos en PostgreSQL al template1:

sudo -u postgres psql template1

Creamos una base de datos, template_postgis, a partir de este template:

create database template_postgis with template = template1;

Actualizamos la tabla pg_database para indicarle que la nueva base de datos es un template:

UPDATE pg_database SET datistemplate = TRUE where datname = ‘template_postgis’;

Nos conectamos a la nueva base de datos:

\c template_postgis

Agregamos las extensiones de PostGIS y le damos acceso a las tablas espaciales a todos:

CREATE LANGUAGE plpgsql ;

\i /usr/share/postgresql-8.3-postgis/lwpostgis.sql;

\i /usr/share/postgresql-8.3-postgis/spatial_ref_sys.sql;

GRANT ALL ON geometry_columns TO PUBLIC;

GRANT ALL ON spatial_ref_sys TO PUBLIC;

Restringimos futuros cambios en la base de datos:

VACUUM FREEZE;

Y con eso hemos terminado.

Fuentes:

Final steps installing PostGIS on Ubuntu para la creación del template.