Map Suite, Oracle, SQL Server, .NET, etc, etc, etc...

Detectar el elemento mas cercano

2011-12-14T14:01:00.000-08:00

UPDATE
Manzanas
SET
IDLocalidad = (SELECT TOP 1 IDLocalidad FROM Localidades
WHERE (Localidades.Geom IS NOT NULL)
ORDER BY Manzanas.Geom.STDistance(Localidades.Geom) ASC
)
WHERE Manzanas.Geom IS NOT NULL

Donde :

Manzanas son poligonos
Localidades son puntos

Contrato temporal

2011-11-30T10:49:00.000-08:00

He terminado mi contrato temporal con el instituto territorial. Básicamente desarrolle unos cuantos controles para la manipulación de cartografía, vinculación de cartografía a registros existentes en la base de datos del estado, importación de archivos shape, etc..

Últimamente he tenido la sensación de salirme de mi trabajo – por aquello del crecimiento - , pero la especialización es buena, al menos me han recomendado que lo haga. He decidido por esta razon seguir tomando cursos de cartografia; para esto, pues obviamente pediré ayuda con los jefes de la empresa para que me apoyen en estos cursos.

Posiblemente en un futuro el jefe me da chance de entrar a gobierno y percibir ingresos extras o igual mi esposa, ya el tiempo lo dirá. Por lo pronto, me quedare aquí y seguiré capacitándome en cursos para no sentir que me quedo atrás en cuanto a tecnologías; pero como dije, todo enfocado en la cartografía.

FORZAR LA DESCONEXION DE UNA BD SQL SERVER

2011-11-22T16:36:00.001-08:00

USE master

EXEC sp_KillProcess 'DB'

-- **** PARA PONER LA BD OFFLINE
ALTER DATABASE DB
SET OFFLINE
WITH ROLLBACK IMMEDIATE

Configurar datos geométricos en Oracle

2011-11-21T09:06:00.000-08:00

En cuestión de capas geométricas es importante considerar algunos aspectos:
- Las capas se manejaran por tablas o vistas?
- La capa geométrica se encuentra registrada en la tabla correspondiente de datos geométricos?.
- Se han creado correctamente los índices, y dentro de este, se encuentra configurado correctamente el tipo de dato geométrico que contendrá la capa?.

Agregar capas geométricas

Por ejemplo, si queremos agregar 4 capas; 2 tablas, TABLA1 y TABLA2; 2 vistas, VW_TABLA1 y VW_TABLA2.

Corremos la siguiente sentencia en Oracle:

SELECT * FROM USER_SDO_GEOM_METADATA
Verificamos que las 4 capas estén registradas; es importante mencionar que USER_SDO_GEOM_METADA es una tabla del schema MDSYS no es una tabla del schema GTO.

En caso de no existir alguna capa registrada la agregamos de la siguiente forma:

INSERT INTO User_SDO_Geom_Metadata (table_name,column_name,diminfo,srid) VALUES ('VW_TABLA2', 'GEOM', MDSYS.SDO_DIM_ARRAY(MDSYS.SDO_DIM_ELEMENT('X',178338,430380,0.5),MDSYS.SDO_DIM_ELEMENT('Y',2202847,2417710,0.5)),32614);
MDSYS.SDO_DIM_ELEMENT(SDO_DIMNAME,SDO_LB,SDO_UB,SDO_TOLERANCE)
MDSYS.SDO_DIM_ELEMENT('X',178338,430380,0.5)

Dónde:

SDO_DIMNAME : indica el parámetro al que se establecerá las coordenadas mínimas y máximas
SDO_LB : Coordenada mínima en SDODIMNAME
SDO_UB : Coordenada máxima en SDODIMNAME
SDO_TOLERANCE : Es una forma de garantizar que tendrás consistencia y resultados aceptables atreves de diferentes plataformas y aplicaciones en pocas palabras déjalas en 0.5 no le muevas..

Si lo comparamos con SQL Server estos son los datos espaciales XMIN, XMAX, YMIN, YMAX
Agregar índices espaciales
Una vez registradas todas las capas en la table USER_SDO_GEOM_METADA, verificamos que los índices se encuentren creados para cada una de las capas. En caso de no existir algún índice espacial lo creamos de la siguiente forma:

CREATE INDEX VW_TABLA2_SIdx ON VW_TABLA2 (GEOM)
INDEXTYPE IS MDSYS.SPATIAL_INDEX PARAMETERS ('LAYER_GTYPE=MULTIPOLYGON')

Como se puede observar, se está creando un índice sobre una vista. Al ejecutar esta sentencia nos arrojara un error como el sigue:

Error que empieza en la línea 3 del comando:
CREATE INDEX VW_TABLA2_SIdx ON VW_TABLE2 (GEOM)
INDEXTYPE IS MDSYS.SPATIAL_INDEX PARAMETERS ('LAYER_GTYPE=MULTIPOLYGON')
Error en la línea de comandos:3 Columna:40
Informe de error:
Error SQL: ORA-01702: una vista no es apropiada aquí
01702. 00000 - "a view is not appropriate here"
*Cause:
*Action:

Aun con este error es importante crear estos índices, posiblemente es un bug de Oracle pero si no lo hacemos de esa manera, la lectura de las capas en thinkgeo nos mandara un error.

Import, Export en Oracle

2011-11-16T10:09:00.000-08:00

Esta entrada es parte de una serie de paginas que escribí para realizar respaldos con ORACLE.

Incluira tambien la preparacion de la base de datos para manipular cartografia..

Este documento seguido se pierde y no se porque - bueno, en realidad si - tenemos un desmandre de documentos regados por todos lados..

Introducción

Hay varios tipos de respaldos en Oracle. Primero, tenemos los respaldos lógicos y respaldo físicos. Los respaldos lógicos se llevan a cabo con utilerías export de Oracle. Los respaldos físicos se llevan a cabo en caliente y frio. Con los respaldos en frio, se debe apagar la base de datos y copiar los datafiles de la base de datos a un respaldo destino. Los respaldos en caliente son los más óptimos para usar RMAN ya que no es necesario apagar la base de datos. RMAN ofrece diferentes tipos de respaldos. Con RMAN se pude respaldar una base de datos completa, respaldar table spaces específicos y también optimizar los respaldos incrementales.
Este documento se enfoca en los respaldos lógicos, Import y Export, los cuales permiten respaldar y restaurar una base de datos completa, así como usuarios y tablas específicas en base a un archivo dump.
Una vez restaurada una base de datos de datos con datos espaciales, se deben tomar en cuenta la configuración de algunas tablas e índices para interactuar con dichos datos.

Respaldar una base de datos
Antes de iniciar con el respaldo es necesario tener instalado el cliente de Oracle y asegurarse que las aplicaciones imp y exp se encuentren en el directorio bin.
Iniciar un ventana de comandos y teclear el comando exp (exportar). Posteriormente pide los datos de autenticación; es importante tomar en cuenta que el usuario con el que se desea autenticar debe tener los permisos para llevar a cabo esta tarea.

Enter array fetch buffer size: 4096 >
Export file: EXPDAT.DMP >
(1)E(ntire database), (2)U(sers), or (3)T(ables): (2)U > 2
Export grants (sÝ/no): sÝ >
Export table data (sÝ/no): sÝ >
Compress extents (sÝ/no): sÝ >
About to export specified users ...
User to be exported: (RETURN to quit) > gto
User to be exported: (RETURN to quit) >

Una vez configurados todo los parametros se da un enter y comenza la exportación de datos en el archivo EXPDAT.DMP..

. exporting pre-schema procedural objects and actions
. exporting foreign function library names for user GTO
. exporting PUBLIC type synonyms
. exporting private type synonyms
. exporting object type definitions for user GTO

About to export GTO's objects ...

. exporting database links
. exporting sequence numbers
. exporting cluster definitions
. about to export GTO's tables via Conventional Path ...
. . exporting table AREASHOMOGENEAS 2 rows exported
. . exporting table CHANGETABLE 202 rows exported

etc, etc, etc..

Restauración de una base de datos
Iniciar un ventana de comandos y teclear el comando imp (importar) . Posteriormente pide los datos de autenticación; el usuario con el que se desea importar el archivo DMP, debe de existir en la base de datos destino y este debe contar con los permisos correspondiente.
Una vez autenticado nos pedirá una serie de parámetros configurables. Como recomendación se dejan los valores por default y solo se da enter por cada parámetro requerido.

Import file: EXPDAT.DMP >
Enter insert buffer size (minimum is 8192) 30720>
Export file created by EXPORT:V11.01.00 via conventional path
import done in WE8MSWIN1252 character set and AL16UTF16 NCHAR character set
List contents of import file only (sÝ/no): no >
Ignore create error due to object existence (sÝ/no): no >
Import grants (sÝ/no): sÝ >
Import table data (sÝ/no): sÝ >
Import entire export file (sÝ/no): no > s (Este parametro se debe cambiar a ‘s’ de los contrario se tendrían que especificar manualmente que tablas queremos importar)

Está básica esta primera entrada; pero yo le batalle para usar los respaldos cold y hot backup. Estuve preguntando a otros que trabajan con Oracle y prefieren usar los respaldos logicos (imp, exp), ya que son mas sencillos de usar..

Window Functions in SQL Server

2011-11-14T14:54:00.001-08:00

Cuantos no hemos intentado hacer esto :
SELECT ID,
Value,
SUM(Value) AS "Sum"
AVG(Value) AS "Avg"
COUNT(Value) AS "Quantity" FROM TestAggregation
GROUP BY ID;

error!!!!!!!!!!!

Msg 8120, Level 16, State 1, Line 2

Column 'TestAggregation.Value' is invalid in the select list because it is not contained in either an aggregate function or the GROUP BY clause.

se podria resolver con consultas anidadas pero esto es inficiente y el rendimiento pues ni que decir..

con window functions lo resolvemos de un manera elegante ...

SELECT ID,
Value,
SUM(Value) OVER(PARTITION BY ID) AS "Sum"
AVG(Value) OVER(PARTITION BY ID) AS "Avg"
COUNT(Value) OVER(PARTITION BY ID) AS "Quantity"
FROM TestAggregation

para mas informacion, checalo en el siguiente enlace...
http://www.simple-talk.com/sql/learn-sql-server/window-functions-in-sql-server/

Acotaciones

2011-11-14T08:32:00.000-08:00

Una muestra de lo que se puede hacer con operaciones geometricas..

El script para acotar poligonos, es uno de los mas complejos y de los cuales me puedo sentir mas orgulloso. Tardé un par de dias en hacer esto, pero ahora, con un solo objeto se puede implementar para cualquier clase que implemente un mapa..

posteriormente estare pegando unas lineas clave para poder implementar esto con sql server y map suite..

Cambiar el SRID de una Capa

2011-11-11T18:10:00.000-08:00

Capa = Tabla con un campo geometry o geography

Hace ya varios ayeres que un nuevo cliente adquirió el sistema (Durango). Se les migro todas las tablas del sistema, se les migraron las tablas de configuracion de los mapas, las tablas de parametros, etc, etc, etc..

Cuando el cliente intento probar el sistema, resulta que no podia ver ni madres en pantalla; solo podian ver los datos alfanumericos...

porque??? dos dias gaste resolviendo este problema y la menera de resolverlo fue cambiando el SRID..

Que es el SRID ..
Cada instancia espacial tiene un identificador de referencia espacial (SRID). El SRID corresponde a un sistema de referencia espacial basado en el elipsoide concreto usado para la creación de mapas de tierra plana o de tierra redonda. Una columna espacial puede contener objetos con SRID diferentes. Sin embargo, sólo se pueden usar instancias espaciales con el mismo SRID al realizar operaciones con métodos de datos espaciales de SQL Server en sus datos. El resultado de cualquier método espacial derivado de dos instancias de datos espaciales sólo es válido si dichas instancias tienen el mismo SRID basado en la misma unidad de medida, dato y proyección usados para determinar las coordenadas de las instancias. Las unidades de medida más comunes de un SRID son metros o metros cuadrados.

En pocas palabras el SRID contiene la informacion para georeferenciar los elementos cartograficos..

resulta que me estaba basando en los datos de guanajuato con el SRID 32614 y Durango es 32613

para cambiar este dato y poder visualizar correctamtente la cartografia fue necesario aplicar la siguiente sentencia :

UPDATE Vialidades SET Geom = geometry::STGeomFromText(Geom.STAsText(), 32613)
where Geom is not null

y listo ya pude ver la cartografia...

Validar topología de un elemento geográfico

2011-11-11T10:19:00.000-08:00

Muchas veces sucede que al digitalizar un poligono hay interseccion de vertices y los segmentos se intersectan. Esto genera un error topologico y es necesario hacer valido el elemento geografico antes de guardarlo en la base de datos o realizar alguna operacion con este.

DECLARE @geomAreaHomogenea as geometry

SELECT @geomAreaHomogenea = geometry::STGeomFromWKB(@Geom, @iSRID)

IF @geomAreaHomogenea.STIsValid() = 0
BEGIN
SET @geomAreaHomogenea=@geomAreaHomogenea.MakeValid()
END

Obtener el indicie espacial para una tabla geometrica

2011-11-10T17:11:00.000-08:00

Un índice espacial sólo se puede crear en una columna espacial, es decir, una columna del tipo geometry o geography

Este sirve para indexar las columnas geometricas y agilizar la busqueda de datos geometricos..

Un indice espacial se compone de 4 datos XMIN, YMIN, XMAX, YMAX tambien llamado extent en map suite..

En el siguiente ejemplos e trata de obtener el indice espacial para las tablas de cierta entidad de la republica..

DECLARE @Geom Geometry
DECLARE @XMin AS DECIMAL(10,0), @YMin AS DECIMAL(10,0), @XMax AS DECIMAL(10,0), @YMax AS DECIMAL(10,0)

SELECT @Geom = Geom FROM Estados WHERE IDEstado = @strEstado

SET @Geom = @Geom.STEnvelope();
SET @XMin = CAST(ROUND(@Geom.STPointN(1).STX, 0) AS DECIMAL(10,0))
SET @YMin = CAST(ROUND(@Geom.STPointN(1).STY, 0) AS DECIMAL(10,0))
SET @XMax = CAST(ROUND(@Geom.STPointN(3).STX, 0) AS DECIMAL(10,0))
SET @YMax = CAST(ROUND(@Geom.STPointN(3).STY, 0) AS DECIMAL(10,0))

ahora, supongamos que no tenemos la entidad pero solo tenemos los municipios.

DECLARE @Geom Geometry SET @Geom = GEOMETRY::STGeomFromText('GEOMETRYCOLLECTION EMPTY', 32614)

select @Geom = @Geom.STUnion(Geom) from dbo.Municipios WHERE Geom IS NOT null;

UPDATE dbo.Estados SET Geom = @Geom WHERE IdEstado = @strIdEstado

la union de todos los municipios da por resultado el contorno de un estado y usando la function geometrica STUnion se logra dicha union..

Geometrias duplicadas

2011-11-10T15:28:00.001-08:00

Este script me costo un buen a pesar de lo simple que se ve. La clave esta en la sentencia %%physloc%%.

SELECT * FROM PREDIOS WHERE
PREDIOS.%%physloc%%
NOT IN (SELECT MIN(b.%%physloc%%)
FROM Predios
GROUP BY Predios.geom.STArea(), Predios.geom.STCentroid().STX, Predios.geom.STCentroid().STY)