aRT 1.7-1

• Estrutura de dados para a matriz de vizinhança no aRT. Funções para ler e escrever no R.
• Funções para gerar GPMs no aRT, bem como funções de preenchimento de celulas com a GPM.

• Corrigir o print do aRTraster para mostrar a quantidade de pixels e outras informações.
• Corrigir o doc do getPixels para dizer que funciona com poligonos.
• Arrumar a questao dos ids para o spatialpoints, de forma a permitir subconjuntos.
• Estudar o pacote RODBC, para ver como abrir bancos de dados no formato mdb.
• Manipular polígonos corretamente, suportando ilhas, buracos e o resultado duplo do buffer.
• problema numElemInstances do querier na versao 3.3.1 phorum
• Trabalhar com Grids (possibilitar ler celulas como grids ou poligonos, e ler raster como grids)
• bh script is not working. problem with importMif. put it in the webpage again.
• trocar o 'delete' do createDb por 'replace'
• mudar o importSpData para o createLayer, passando um argumento 'data' com o objeto Spatial

• colocar o getID como funcao da classe Spatial, para evitar dizer que nao esta implementado quando o objeto for NULL, por exemplo.

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)