<!DOCTYPE html PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN">
<html>
<head>
<meta content="text/html;charset=ISO-8859-1" http-equiv="Content-Type">
<title></title>
</head>
<body bgcolor="#ffffff" text="#000000">
<br>
Hi dear list. This will be a long one:<br>
<br>
I chopped a map into multiple adjacent rectangles (like a grid), in
the hopes of feeding each rectangle (or subzone, or cell, call what you
like :) trough a queryByRectangle() in order to get the number of
shapes present at that sub-area of a layer, and use the resulting
numbers on a statistic analysis of the number of ocurrences.<br>
Everyhing seems to be working fine, except my queries return weird
numResults.<br>
For example, in the case below I am creating 4 cells and ploting 10
points on the upper right cell. however logs ALSO return 10
non-existing points in the upper left cell, where there is nothing. If
I up the number of cells (how the final application will/should be
like), the count numbers go completelly insane.<br>
Am I misunderstanding the use of queryByRect? I also tried drawing
the rectangles as shapes and doing a queryByShape and it gets worse: he
counts all the 10 points as being in present in all cells. Geting the
index of these shapes show that they are the all the same set of 10,
magically "appearing" in other areas (they arent there, just
numResults() think they are)<br>
Any help is appreciated, code below, the log output below it:<br>
<br>
<br>
//Creating layer GRID - This layer is debug only<br>
layerObj grid = new layerObj(mapa);<br>
grid.setType(edu.umn.gis.mapscript.MS_LAYER_TYPE.MS_LAYER_POLYGON);<br>
grid.setName("Grid");<br>
classObj gclas = new classObj(grid);<br>
styleObj gstyle = new styleObj(gclas);<br>
colorObj gcolor = new colorObj(0,0,0,0);<br>
gstyle.setOutlinecolor(gcolor);<br>
gclas.setNumstyles(1);<br>
gclas.insertStyle(gstyle, 0);<br>
grid.setStatus(mapscript.MS_ON);<br>
<br>
//Divisor: there will be 'divisor ^ 2' cells<br>
int divisor = 2;<br>
<br>
double minX = mapa.getExtent().getMinx();<br>
double minY = mapa.getExtent().getMiny();<br>
double maxX = mapa.getExtent().getMaxx();<br>
double maxY = mapa.getExtent().getMaxy(); <br>
<br>
double deviationX = (maxX - minX)/divisor;<br>
double deviationY = (maxY - minY)/divisor;<br>
<br>
double[] lineX = new double[divisor+1];<br>
double[] lineY = new double[divisor+1];<br>
double stepX = minX;<br>
double stepY = minY;<br>
lineY[0] = minY;<br>
lineX[0] = minX;<br>
<br>
log.debug("Extent:" + mapa.getExtent()); <br>
log.debug("devX: " + deviationX + " devY: " + deviationY);<br>
<br>
int h = 1;<br>
while (stepY <= maxY){<br>
stepY += deviationY; <br>
if(stepY <= maxY)<br>
lineY[h] = stepY;<br>
h++;<br>
//log.debug("stepY: " + stepY);<br>
}<br>
h = 1;<br>
while (stepX <= maxX){<br>
stepX += deviationX; <br>
if(stepX <= maxX)<br>
lineX[h] = stepX;<br>
h++;<br>
//log.debug("stepX: " + stepX);<br>
}<br>
<br>
//Hack to avoid float number imprecisions at the hightest extents<br>
lineY[lineY.length -1] = Double.valueOf(maxY);<br>
lineX[lineX.length -1] = Double.valueOf(maxX);<br>
<br>
int j = 0; //DEBUG Indice independente<br>
List<rectObj> cell = new ArrayList<rectObj>();<br>
for(int i = 0; i < (lineX.length -1); i ++){<br>
for(h = 0; h < (lineY.length -1); h++){<br>
rectObj tempRect = new rectObj(lineX[i], lineY[h],
lineX[i+1], lineY[h+1], mapscript.MS_FALSE);<br>
cell.add(tempRect);<br>
<br>
//Actually drawing the grid is for pure debug purposes<br>
//TODO Remove from final version<br>
shapeObj tempPoly = tempRect.toPolygon();<br>
tempPoly.setIndex(j);<br>
tempPoly.setTileindex(j);<br>
grid.addFeature(tempPoly);<br>
<br>
j++;<br>
}<br>
}<br>
<br>
//Creating layer that will house the points and adding a ramdon
debug set<br>
layerObj crimes = new layerObj(mapa);<br>
crimes.setType(edu.umn.gis.mapscript.MS_LAYER_TYPE.MS_LAYER_POINT);<br>
crimes.setName("Crimes");<br>
classObj clas = new classObj(crimes);<br>
styleObj style = new styleObj(clas);<br>
style.setSymbol(IPSUSPEITO);<br>
clas.setNumstyles(1);<br>
clas.insertStyle(style, 0);<br>
<br>
shapeObj gridShape = new shapeObj(mapscript.MS_SHAPEFILE_POINT);<br>
lineObj gridLinha = new lineObj();<br>
pointObj gridPonto = new pointObj(0.0, 0.0, 0.0);<br>
<br>
Random rnd = new Random();<br>
for(int i = 0; i < 10; i++){<br>
gridLinha = new lineObj();<br>
<br>
//rect 3 - Meaning points will exist at the upper right cell<br>
double lat = rnd.nextDouble() * 180;<br>
double lon = rnd.nextDouble() * 80;<br>
<br>
gridPonto.setXY(lat, lon, 0.0);<br>
gridLinha.add(gridPonto);<br>
gridShape.add(gridLinha);<br>
gridShape.setIndex(i);<br>
gridShape.setTileindex(i);<br>
crimes.addFeature(gridShape); //préstenssão!<br>
}<br>
<br>
crimes.setTransparency(mapscript.MS_GD_ALPHA);<br>
String template = "pathe/to/my/dir/crimes.html"; //changed for this
email´s purpose ;)<br>
crimes.setTemplate(template);<br>
grid.setTemplate(gridtemplate);<br>
crimes.setStatus(mapscript.MS_ON);<br>
<br>
//Testes: queries<br>
Iterator cellIterator = cell.iterator();<br>
int f = 0; //debug<br>
while(cellIterator.hasNext()){<br>
rectObj rect = (rectObj)cellIterator.next();<br>
crimes.queryByRect(mapa, rect);<br>
crimes.open();<br>
resultCacheObj results = crimes.getResults();<br>
log.debug("---------------------");<br>
log.debug("Retangulo -> " + rect.getMinx() + ", " +
rect.getMiny() + " -- " +<br>
rect.getMaxx() + ", " +
rect.getMaxy());<br>
log.debug("NumResults -> " + results.getNumresults());<br>
<br>
//log.debug("=====================");<br>
//if(results != null){<br>
//Trecho abaixo recupera os shapes relacionados (pontos).
Irrelevante p/ o gráfico.<br>
/* for(int e = 0; e < results.getNumresults(); e++){<br>
resultCacheMemberObj result = results.getResult(e);<br>
log.debug("result("+e+") name: " + result.toString());<br>
log.debug("crimes.getFeature("+e+"):
getShapeindex("+result.getShapeindex()+"), Tileindex("+
result.getTileindex()+")");<br>
shapeObj ocorrencia =
crimes.getFeature(result.getShapeindex(), result.getTileindex());<br>
log.debug("Feature " + e + " -> " +
ocorrencia.toString());<br>
ocorrencia.getClass();<br>
}*/<br>
//}<br>
//log.debug("=====================");<br>
log.debug("---------------------");<br>
f++;<br>
crimes.close();<br>
<br>
}<br>
<br>
The above code draws the 4 cell grid on my map and correctly plots
the 10 points at the upper right cell, but this is the logs result:<br>
<br>
Extent:{ 'minx': -180 , 'miny': -92.18668454494383 , 'maxx': 180 ,
'maxy': 85.79084354494383 }<br>
devX: 180.0 devY: 88.98876404494382<br>
---------------------<br>
Retangulo -> -180.0, -92.18668454494383 -- 0.0, -3.1979205000000093<br>
NumResults -> 0 <-- This one is right<br>
---------------------<br>
---------------------<br>
Retangulo -> -180.0, -3.1979205000000093 -- 0.0, 85.79084354494383<br>
NumResults -> 10 <-- This one is WRONG!<br>
---------------------<br>
---------------------<br>
Retangulo -> 0.0, -92.18668454494383 -- 180.0, -3.1979205000000093<br>
NumResults -> 0 <-- This one is right<br>
---------------------<br>
---------------------<br>
Retangulo -> 0.0, -3.1979205000000093 -- 180.0, 85.79084354494383<br>
NumResults -> 10 <-- This one is right<br>
---------------------<br>
<br>
Any ideas?<br>
Thanks a lot in advance<br>
<br>
<div class="moz-signature">-- <br>
<font color="#000000" face="Verdana, Arial, Helvetica" size="1">
<br>
<br>
<font color="#000000" size="2"><b>Rodrigo Del C. Andrade</b></font><br>
<font color="#606060"><i>Estagiário Nível Superior</i></font><br>
<font color="#606060"><i>SIC - SSE - Soluções Segurança Pública</i></font><br>
<br>
<font color="#000080"><b>DÍGITRO TECNOLOGIA</b></font><br>
<font color="#000000"><b>E-mail:</b></font>
<a href="mailto:rodrigo.andrade@digitro.com.br"><font color="#606060">rodrigo.andrade@digitro.com.br</font></a><br>
<font color="#000000"><b>Fone:</b></font>
<font color="#000000">+55 48 3281-7390 / +55 48 3281-7000</font><br>
<font color="#000000"><b>Fax:</b></font>
<font color="#000000">+55 48 3281-7299</font><br>
<font color="#000000"><b>Site:</b></font>
<a href="http://www.digitro.com"><font color="#606060">www.digitro.com</font></a>
</font></div>
</body>
</html>