aRT 1.6-3 (ago/2009)

A nova versao do aRT irá funcionar com a TerraLib 3.3.0, ao invés da 3.3.1, que seria a mais nova versão. Isto porque a 3.3.1 apresenta um problema que não permite a leitura de dados do banco (ver o FIX ME).

• Verificar erros terralib da versão 1.6-2:
  • erro MySQL windows quando executa o script recife.R: "too long …". ver problema e relatar para o pessoal da TerraLib, antes da nova versao (15/out). esquesito, os tamanhos deste campo sao iguais para linux e windows.
  • updateColumns (consulta SQL errada). phorum
  • bug em src/terralib/kernel/TeGraph.cpp (nao compila). phorum
  • Postgre showDbs(). quando nao tem bancos disponiveis, ele coloca o rows == 0, mas o cols!= 0, o que gera um CORE DUMP.
  • poligonos com buracos e ilhas

• openTheme(db, "theme") (directly from the database, instead of the layer)
• reload(db) to reconnect to the database, keeping consistent all children (layers, tables). Solução nao trivial, requer modificações em várias classes do aRT. melhor esperar pela nova versão da terralib.
• atualizar para R 2.9-1
• atualizar para TerraLib mais recente

• ler dado celular do banco de dados
• CLIP the PolyLine to the area of the Polygons. The polygons are within the larger main_main roads Polyline. Edzer: Hard to do; I just found out it can be done with aRT, which requires you to run TerraLib as the GIS on your computer. e-mail
• plotar dados celulares
• fatiamento de dados celulares
• PJ: A razao principal é o que eu mencionei que sendo um SpatialPolygon (ou derivado) espera-se encontrar isto nos objetos portanto me parece consistente ter isto lá. objetos SpatialPolyginos cujos slots labpt e area estao como numeric(0). Pedro: Colocar parametros loadarea e loadlabpt como default true para calcular estes valores do banco.
• Operadores como Desvio padrão, classe majoritaria,contagem, presença precisam de uma documentação mais específica.

Repensar a classe. perhaps setVisual(theme, color="blue") instead of setVisual(theme, visualPolygons(color="blue"))? Esta simplificação não é possível, porque um tema pode ter mais de uma geometria, e é necessário explicitar de qual geometria estamos alterando o visual.

If I have a polygon shapefile, and I want to rasterize each polygon given a certain grid cell resolution, is there a way to do this completely within R? I'm considering porting our vector↔raster bridge "starspan" (http://starspan.casil.ucdavis.edu/doku/doku.php) to R code, but this is a critical portion of the algorithm. If the answer is "yes" – does that answer include "weird" polygons like donuts? Roger: In principle, the overlay() method in the sp package for "SpatialPolygon" and "SpatialPoints" objects, specifically the:

sp:::pointsInSpatialPolygons

function is where to start. The underlying assumption is that raster cell support can be cast to cell centre point support. This function first finds candidate polygons by testing polygon bounding boxes, then hands off to:

sp:::pointsInPolys2

which calls:

sp:::pointsInPolygons

Here you can see that hole handling at the per-ring stage is attempted, but does depend on the imported rings declaring their hole status correctly. It is possible to use checkPolygonsHoles() in maptools, but this takes time.

Hope this helps, porting starspans would be very welcome - I guess users at the moment use it loose-coupled by reading and writing files.

The overlay methods are open to optimisation (you'll see that bits of sp:::pointsInSpatialPolygons moves heavier-lifting out to C)

Paulo: vc viu que o principal problema era reconhecer distritos com mais de um poligono em algumas operacoes. Vou ver como ficou isto na sua versao.

• verificar todos os vignettes
• planejar um curso de 4 ou 5 aulas, com exercícios se possível
• scripts para o livro 'analise espacial de dados geograficos'?