2011-11-23

Convertendo pontos de SIRGAS para SAD69

Fonte: uma lista de pontos (arquivo CSV) em SIRGAS 2000

Objetivo: gerar um arquivo DXF com os pontos convertidos - pois serão abertos em AutoCAD.

Roteiro:

1 - a partir do arquivo CSV, em uma planilha gera-se vários comandos para criar um shapefile a partir dos pontos UTM:

2. depois, execute os comandos gerados num prompt do MS4W, por exemplo, que já traz os binários da ShapeLib:
$ shpcreate original_sirgas point
# cria um shapefile vazio (só arquivos .shp e .shx) de tipo PONTOS
E em seguida execute os vários 'shpadd' para inserir os pontos no shapefile:
$ shpadd original_sirgas X1 Y1
$ shpadd original_sirgas X2 Y2
$ shpadd original_sirgas Xn Yn
...

3. em seguida, via OGR, converta o shapefile com pontos em SIRGAS 2000 para um DXF em SAD69:
$ ogr2ogr -f DXF -t_srs EPSG:29193 -s_srs EPSG:31983 final_sad69.dxf original_sirgas2000.shp
Os números em destaque são os identificadores/códigos dos sistemas de coordenadas, de acordo com padrão EPSG, usados pela biblioteca PROJ4 e se referem ao fuso UTM 23S dos sistemas destino (SAD69) e origem (SIRGAS).
Para uma lista completa dos identificadores e parâmetros de projeção, basta visualizar o arquivo C:\ms4w\proj\nad\epsg ( não tem extensão! )

E pronto! Tem-se o arquivo 'final_sad69.dxf' à mão, com coordenadas em SAD69 e pode ser lido no CAD.

2011-11-21

MS-SQL 2008 Spatial experiments

...notes of some experiments from 2011 beginning
     - mixing up SQL 2008 + MS4W 3.0.1 ( MapServer 5.6.6 )
SQL Server 2008 Setup
  • packages to install:
  • tried and tried, but it didn't work at my old Windows 7 Pro
  • worked fine on Windows 2003 Server - even under Virtual Box
  • at installer 7th step ("Server configuration"), define "NT AUTHORITY\LOCAL SERVICE" as credentials for 'SQL Server Database Engine' service
  • and please, define a good password for sa user and do not forget to enable 'mixed logins'
Some notes on SQL 2008 geometry handling
CREATE SPATIAL INDEX sidx_mytable ON 
myscheme.mytable ( mygeomfield )
WITH (
BOUNDING_BOX = ( -180, -90, 180, 90 ) )
  • when a table has a 'geometry' field, results will also be shown as a map at MS SQL Studio GUI 

  • OGC SFS methods are invoked from geometry fields. Eg.:
SELECT v.*, v.mygeomfield.STAsText() FROM mytable AS v

Testing connection using OGR
version installed: MS4W 3.0.1 / MapServer 5.6.6
  • create a system ODBC DSN
  • run an ogrinfo test:
$ ogrinfo -so "MSSQL:dsn=MY_DSN;server=(local);database=DBNAME_HERE;tables=myscheme.mytable;uid=USER_HERE;pwd=PWD_HERE" mytable
  • it will output something like:
INFO: Open of `MSSQL:.....'
      using driver `MSSQLSpatial' successful.

TO BE CONTINUED here!...