parallel/util/PartitionUtil.java

package appl.parallel.util;

import java.awt.Rectangle;
import java.awt.image.DataBuffer;

import appl.parallel.data.WritableGridArrayPartition;
import appl.parallel.spmd.split.SplittableResource;
import appl.parallel.spmd.split.WritableGridPartition;
import de.appl.util.RasterMetaData;
import de.schmitzm.geotools.data.WritableGrid;

/**
 * See method description for details.
 * 
 * @author Dominik Appl
 */
public class PartitionUtil {

        /**
                 * Returns partition of a {@link WritableGrid} as a {@link WritableGridArray}.
                 * Notice, that the minX and minY values are ignored and will be set to 0
                 * (see {@link WritableGrid} for more information about minX and minY). If
                 * the given {@link WritableGrid} is not an instance of
                 * {@link WritableGridPartition} it is assumed that the topLeft corner of
                 * the source grid is (0,0) (like in standard image processing).<br>
                 * If the partition is a {@link WritableGridPartition} it is assumed that
                 * the the bounds of the two Partition refer to the same coordinate system.
                 * Real coordinates are supposed to reference the South(!)-West Corner.
                 * 
                 * @param sourceGrid
                 *            the grid to be partitioned
                 * @param partitionBounds
                 *            the rectangle describing the partition
                 * @param rootID
                 *            the id to identify the root partition
                 * @return the new partition including the given corners
                 * @see SplittableResource#getRootID()
                 */
                public static WritableGridPartition getPartitialGrid2D(WritableGrid sourceGrid,
                                Rectangle partitionBounds, int rootID) {
                        
        //                              (int) partitionBounds.getX(),
        //                              (int) partitionBounds.getY(),
        //                              (int) (partitionBounds.getX() + partitionBounds.getWidth() - 1),
        //                              (int) (partitionBounds.getY() + partitionBounds.getHeight() - 1));
                        
        //              get Metadata of the big Grid
                        RasterMetaData bigMeta = new RasterMetaData(sourceGrid);
        
                        //create MetaData for the new partitial Grid
                        //Height an width of the partitial Grid
                        int pHeight = (int) partitionBounds.getHeight();
                        int pWidth = (int)partitionBounds.getWidth();
        
                        //real coordinates of the partitial Grid (which are referenced by the SW Corner)
                        double pX = partitionBounds.getX() * bigMeta.getCellWidth() + bigMeta.getX();
                        double pY = (partitionBounds.getY() + partitionBounds.getHeight() - 1) * bigMeta.getCellHeight() + bigMeta.getY();
                        double pRealWidth = pWidth * bigMeta.getCellWidth();
                        double pRealHeight = pHeight * bigMeta.getCellHeight();
        
                        //the MetaData of the new Grid
                        RasterMetaData pMetaData = new RasterMetaData(bigMeta.getDataType(),
                                        pWidth, pHeight, 0, 0, pX, pY, pRealWidth, pRealHeight, bigMeta.getCoordinateReferenceSystem());
        
                        //create empty return Grid
                        WritableGridPartition pGrid;
                        switch(bigMeta.getDataType()){
                        case DataBuffer.TYPE_INT:
                                pGrid = new WritableGridArrayPartition.Integer(pMetaData,rootID,partitionBounds);       
                                break;
                        case DataBuffer.TYPE_DOUBLE:
                                pGrid = new WritableGridArrayPartition.Double(pMetaData,rootID,partitionBounds);
                                break;
                        default:
                                pGrid = new WritableGridArrayPartition.Float(pMetaData,rootID,partitionBounds); 
                        }
                        
                        //if the sourceGrid is a partition itself then we do not want to to start reading from (0,0)
                        //but shifted relative to the global coordinate system. So we calculate now the start values:
                        
                        //for standard WritableGrids these are simply the bounds of the new partition...
                        int startX = (int) partitionBounds.getX();
                        int startY = (int) partitionBounds.getY();
                        
                        // ...but for partitions these values are shifted:
                        if(sourceGrid instanceof WritableGridPartition){
                                Rectangle srcBounds = ((WritableGridPartition) sourceGrid).getPartitionBounds();
                                startX -= srcBounds.getX();
                                startY -= srcBounds.getY();
                        }
                        //read the data into the raster from left to right and downwards
                        int x=0,y=0;
                        try{
                        for (y = 0; y < pHeight; y++)
                                for (x = 0; x < pWidth; x++)
                                        pGrid.setRasterSample(
                                                        sourceGrid.getRasterSample(startX  + x, startY + y),
                                                        x, 
                                                        y);
                        }catch (Exception e) {
                                System.out.println("stop");
                        }
                        return pGrid;
                }

        /**
             * Overwrites the data at the location specified by the {@link Rectangle} with
             * the given partition-data. If the given baseGrid is NOT an instance
                 * of {@link WritableGridPartition} it is assumed that its the topLeft corner (0,0). 
                 * If it IS a {@link WritableGridPartition} it is assumed that the the
                 * bounds of the two Partition refer to the same coordinate system and the partition
                 * is inserted at the correct absolut position.
             * 
                 * @param baseGrid the grid in which the data is inserted
                 * @param gridPartition the grid to be inserted
                 * @param partitionBounds the excact location in coordinates of the baseGrid
                 */
                public static void setPartition(WritableGrid baseGrid, WritableGrid gridPartition, Rectangle partitionBounds) {
                        //Check if partition and the rectangle fit
                        if((gridPartition.getWidth()!=partitionBounds.getWidth()) ||
                                        (gridPartition.getHeight()!=partitionBounds.getHeight()))
                          throw new UnsupportedOperationException("The partition width/height does not match the rectangle width/height");
                        
                        //if the baseeGrid is a partition itself then we do not want to to start writing from (0,0)
                        //but shifted relative to the global coordinate system. So we calculate now the start values:
                        
                        //for standard WritableGrids these are simply the bounds of the new partition...
                        int startX = (int) partitionBounds.getX();
                        int startY = (int) partitionBounds.getY();
                        
                        // ...but for partitions these values are shifted:
                        if(baseGrid instanceof WritableGridPartition){
                                Rectangle targetBounds = ((WritableGridPartition) baseGrid).getPartitionBounds();
                                startX -= targetBounds.getX();
                                startY -= targetBounds.getY();
                        }
                        //check if the partition fits in
        //              if(((partitionBounds.getX()+partitionBounds.getWidth()-1)>baseGrid.getWidth()  ||
        //                              (partitionBounds.getY()+partitionBounds.getHeight()-1)>baseGrid.getHeight()))
        //                              System.out.println("X");
                                //throw new UnsupportedOperationException("The partition does not fit in the Grid!");
                        //write the partition into the Grid
                        
                        for(int y = 0; y < partitionBounds.getHeight();y++)
                                for(int x = 0; x < partitionBounds.getWidth(); x++)
                                        baseGrid.setRasterSample(
                                                                                        gridPartition.getRasterSample(x,y),
                                                                                        x + startX, 
                                                                                        y + startY);
                        
                }

}
1	mojays	2	package appl.parallel.util;
2
3			import java.awt.Rectangle;
4			import java.awt.image.DataBuffer;
5
6			import appl.parallel.data.WritableGridArrayPartition;
7			import appl.parallel.spmd.split.SplittableResource;
8			import appl.parallel.spmd.split.WritableGridPartition;
9	mojays	114	import de.appl.util.RasterMetaData;
10			import de.schmitzm.geotools.data.WritableGrid;
11	mojays	2
12			/**
13			* See method description for details.
14			*
15			* @author Dominik Appl
16			*/
17			public class PartitionUtil {
18
19			/**
20			* Returns partition of a {@link WritableGrid} as a {@link WritableGridArray}.
21			* Notice, that the minX and minY values are ignored and will be set to 0
22			* (see {@link WritableGrid} for more information about minX and minY). If
23			* the given {@link WritableGrid} is not an instance of
24			* {@link WritableGridPartition} it is assumed that the topLeft corner of
25			* the source grid is (0,0) (like in standard image processing).<br>
26			* If the partition is a {@link WritableGridPartition} it is assumed that
27			* the the bounds of the two Partition refer to the same coordinate system.
28			* Real coordinates are supposed to reference the South(!)-West Corner.
29			*
30			* @param sourceGrid
31			* the grid to be partitioned
32			* @param partitionBounds
33			* the rectangle describing the partition
34			* @param rootID
35			* the id to identify the root partition
36			* @return the new partition including the given corners
37			* @see SplittableResource#getRootID()
38			*/
39			public static WritableGridPartition getPartitialGrid2D(WritableGrid sourceGrid,
40			Rectangle partitionBounds, int rootID) {
41
42			// (int) partitionBounds.getX(),
43			// (int) partitionBounds.getY(),
44			// (int) (partitionBounds.getX() + partitionBounds.getWidth() - 1),
45			// (int) (partitionBounds.getY() + partitionBounds.getHeight() - 1));
46
47			// get Metadata of the big Grid
48			RasterMetaData bigMeta = new RasterMetaData(sourceGrid);
49
50			//create MetaData for the new partitial Grid
51			//Height an width of the partitial Grid
52			int pHeight = (int) partitionBounds.getHeight();
53			int pWidth = (int)partitionBounds.getWidth();
54
55			//real coordinates of the partitial Grid (which are referenced by the SW Corner)
56			double pX = partitionBounds.getX() * bigMeta.getCellWidth() + bigMeta.getX();
57			double pY = (partitionBounds.getY() + partitionBounds.getHeight() - 1) * bigMeta.getCellHeight() + bigMeta.getY();
58			double pRealWidth = pWidth * bigMeta.getCellWidth();
59			double pRealHeight = pHeight * bigMeta.getCellHeight();
60
61			//the MetaData of the new Grid
62			RasterMetaData pMetaData = new RasterMetaData(bigMeta.getDataType(),
63			pWidth, pHeight, 0, 0, pX, pY, pRealWidth, pRealHeight, bigMeta.getCoordinateReferenceSystem());
64
65			//create empty return Grid
66			WritableGridPartition pGrid;
67			switch(bigMeta.getDataType()){
68			case DataBuffer.TYPE_INT:
69			pGrid = new WritableGridArrayPartition.Integer(pMetaData,rootID,partitionBounds);
70			break;
71			case DataBuffer.TYPE_DOUBLE:
72			pGrid = new WritableGridArrayPartition.Double(pMetaData,rootID,partitionBounds);
73			break;
74			default:
75			pGrid = new WritableGridArrayPartition.Float(pMetaData,rootID,partitionBounds);
76			}
77
78			//if the sourceGrid is a partition itself then we do not want to to start reading from (0,0)
79			//but shifted relative to the global coordinate system. So we calculate now the start values:
80
81			//for standard WritableGrids these are simply the bounds of the new partition...
82			int startX = (int) partitionBounds.getX();
83			int startY = (int) partitionBounds.getY();
84
85			// ...but for partitions these values are shifted:
86			if(sourceGrid instanceof WritableGridPartition){
87			Rectangle srcBounds = ((WritableGridPartition) sourceGrid).getPartitionBounds();
88			startX -= srcBounds.getX();
89			startY -= srcBounds.getY();
90			}
91			//read the data into the raster from left to right and downwards
92			int x=0,y=0;
93			try{
94			for (y = 0; y < pHeight; y++)
95			for (x = 0; x < pWidth; x++)
96			pGrid.setRasterSample(
97			sourceGrid.getRasterSample(startX + x, startY + y),
98			x,
99			y);
100			}catch (Exception e) {
101			System.out.println("stop");
102			}
103			return pGrid;
104			}
105
106			/**
107			* Overwrites the data at the location specified by the {@link Rectangle} with
108			* the given partition-data. If the given baseGrid is NOT an instance
109			* of {@link WritableGridPartition} it is assumed that its the topLeft corner (0,0).
110			* If it IS a {@link WritableGridPartition} it is assumed that the the
111			* bounds of the two Partition refer to the same coordinate system and the partition
112			* is inserted at the correct absolut position.
113			*
114			* @param baseGrid the grid in which the data is inserted
115			* @param gridPartition the grid to be inserted
116			* @param partitionBounds the excact location in coordinates of the baseGrid
117			*/
118			public static void setPartition(WritableGrid baseGrid, WritableGrid gridPartition, Rectangle partitionBounds) {
119			//Check if partition and the rectangle fit
120			if((gridPartition.getWidth()!=partitionBounds.getWidth()) \|\|
121			(gridPartition.getHeight()!=partitionBounds.getHeight()))
122			throw new UnsupportedOperationException("The partition width/height does not match the rectangle width/height");
123
124			//if the baseeGrid is a partition itself then we do not want to to start writing from (0,0)
125			//but shifted relative to the global coordinate system. So we calculate now the start values:
126
127			//for standard WritableGrids these are simply the bounds of the new partition...
128			int startX = (int) partitionBounds.getX();
129			int startY = (int) partitionBounds.getY();
130
131			// ...but for partitions these values are shifted:
132			if(baseGrid instanceof WritableGridPartition){
133			Rectangle targetBounds = ((WritableGridPartition) baseGrid).getPartitionBounds();
134			startX -= targetBounds.getX();
135			startY -= targetBounds.getY();
136			}
137			//check if the partition fits in
138			// if(((partitionBounds.getX()+partitionBounds.getWidth()-1)>baseGrid.getWidth() \|\|
139			// (partitionBounds.getY()+partitionBounds.getHeight()-1)>baseGrid.getHeight()))
140			// System.out.println("X");
141			//throw new UnsupportedOperationException("The partition does not fit in the Grid!");
142			//write the partition into the Grid
143
144			for(int y = 0; y < partitionBounds.getHeight();y++)
145			for(int x = 0; x < partitionBounds.getWidth(); x++)
146			baseGrid.setRasterSample(
147			gridPartition.getRasterSample(x,y),
148			x + startX,
149			y + startY);
150
151			}
152
153			}