private unsafe void TileOperationAction(IPointDataTileChunk chunk) { //ushort gridSize = (ushort)Math.Sqrt(chunk.Tile.PointCount); //var grid = new Grid<List<int>>(gridSize, gridSize, null, true); //UQuantizedExtent3D quantizedExtent = chunk.Tile.QuantizedExtent; //double cellSizeX = (double)quantizedExtent.RangeX / gridSize; //double cellSizeY = (double)quantizedExtent.RangeY / gridSize; //int i = 0; //byte* pb = chunk.PointDataPtr; //while (pb < chunk.PointDataEndPtr) //{ // var p = (SQuantizedPoint3D*)pb; // int pixelX = (int)(((*p).X - quantizedExtent.MinX) / cellSizeX); // int pixelY = (int)(((*p).Y - quantizedExtent.MinY) / cellSizeY); // if (grid.Data[pixelX, pixelY] == null) // grid.Data[pixelX, pixelY] = new List<int>(); // grid.Data[pixelX, pixelY].Add(i); // pb += chunk.PointSizeBytes; // ++i; //} //uint maxRegionCount = 0; //int windowSize = 10; //int windowRadius = windowSize / 2; //for (int x = 0; x < grid.SizeX; x++) //{ // for (int y = 0; y < grid.SizeY; y++) // { // var pointIndices = grid.Data[x, y]; // if (pointIndices != null) // { // uint regionCount = 0; // for (int x1 = x - windowRadius; x1 < x + windowRadius; x1++) // { // for (int y1 = y - windowRadius; y1 < y + windowRadius; y1++) // { // var pointIndices1 = grid.Data[x, y]; // if(pointIndices1 != null) // regionCount += (uint)pointIndices1.Count; // } // } // maxRegionCount = Math.Max(maxRegionCount, regionCount); // for (int index = 0; index < pointIndices.Count; index++) // { // UQuantizedPoint3D* p = (UQuantizedPoint3D*)(chunk.PointDataPtr + chunk.PointSizeBytes * pointIndices[index]); // (*p).Z = regionCount; // } // } // } //} //pb = chunk.PointDataPtr; //while (pb < chunk.PointDataEndPtr) //{ // UQuantizedPoint3D* p = (UQuantizedPoint3D*)pb; // if ((*p).Z > maxRegionCount) // (*p).Z = maxRegionCount; // (*p).Z = (*p).Z * quantizedExtent.RangeZ / maxRegionCount + quantizedExtent.MinZ; // pb += chunk.PointSizeBytes; //} }
private unsafe void TileOperationAction(IPointDataTileChunk chunk) { //ushort gridSize = (ushort)Math.Sqrt(chunk.Tile.PointCount); //var grid = new Grid<List<int>>(gridSize, gridSize, null, true); //UQuantizedExtent3D quantizedExtent = chunk.Tile.QuantizedExtent; //double cellSizeX = (double)quantizedExtent.RangeX / gridSize; //double cellSizeY = (double)quantizedExtent.RangeY / gridSize; //int i = 0; //byte* pb = chunk.PointDataPtr; //while (pb < chunk.PointDataEndPtr) //{ // var p = (SQuantizedPoint3D*)pb; // int pixelX = (int)(((*p).X - quantizedExtent.MinX) / cellSizeX); // int pixelY = (int)(((*p).Y - quantizedExtent.MinY) / cellSizeY); // if (grid.Data[pixelX, pixelY] == null) // grid.Data[pixelX, pixelY] = new List<int>(); // grid.Data[pixelX, pixelY].Add(i); // pb += chunk.PointSizeBytes; // ++i; //} //uint maxRegionCount = 0; //int windowSize = 10; //int windowRadius = windowSize / 2; //for (int x = 0; x < grid.SizeX; x++) //{ // for (int y = 0; y < grid.SizeY; y++) // { // var pointIndices = grid.Data[x, y]; // if (pointIndices != null) // { // uint regionCount = 0; // for (int x1 = x - windowRadius; x1 < x + windowRadius; x1++) // { // for (int y1 = y - windowRadius; y1 < y + windowRadius; y1++) // { // var pointIndices1 = grid.Data[x, y]; // if(pointIndices1 != null) // regionCount += (uint)pointIndices1.Count; // } // } // maxRegionCount = Math.Max(maxRegionCount, regionCount); // for (int index = 0; index < pointIndices.Count; index++) // { // UQuantizedPoint3D* p = (UQuantizedPoint3D*)(chunk.PointDataPtr + chunk.PointSizeBytes * pointIndices[index]); // (*p).Z = regionCount; // } // } // } //} //pb = chunk.PointDataPtr; //while (pb < chunk.PointDataEndPtr) //{ // UQuantizedPoint3D* p = (UQuantizedPoint3D*)pb; // if ((*p).Z > maxRegionCount) // (*p).Z = maxRegionCount; // (*p).Z = (*p).Z * quantizedExtent.RangeZ / maxRegionCount + quantizedExtent.MinZ; // pb += chunk.PointSizeBytes; //} }