Application: gvSIG desktop » gvSIG 3D

import java.util.ArrayList;

import org.gvsig.fmap.raster.layers.FLyrRasterSE;

				_logger.warn("computed optimum level = "+computeOptimumLevel(layer,20));

				//hlayer.setMaxResolution(computeOptimumLevel(layer,20)-1);

	public int computeOptimumLevel(FLayer layer, int maxLevel) {

		int level = maxLevel;

		if(!(layer instanceof FLyrRasterSE))

				return level;

		Rectangle2D ext = null;

        try {

            ext = layer.getFullExtent();

        } catch (ReadDriverException e1) {

            NotificationManager.addError("Error al obtener el extent", e1);

        }

		ArrayList attr = ((FLyrRasterSE)layer).getAttributes();

        int width = 0, height = 0;

        for (int i=0; i<attr.size(); i++) {

            Object[] a = (Object []) attr.get(i);

            if (a[0].toString().equals("Width"))

                width = ((Integer)a[1]).intValue();

            if (a[0].toString().equals("Height"))

                height = ((Integer)a[1]).intValue();

        }

        if( ext != null &&

                width != 0 &&

                height != 0){

        	Extent destinationExtents = _terrain.getExtent();

			double destination_xRange = destinationExtents.xMax()-destinationExtents.xMin();

		    double destination_yRange = destinationExtents.yMax()-destinationExtents.yMin();

		    double AR = destination_xRange / destination_yRange;

		    int C1 = 1;

		    int R1 = 1;

		    boolean swapAxis = AR<1.0;

		    if (swapAxis) AR = 1.0/AR;

		    double lower_AR = Math.floor(AR);

		    double upper_AR = Math.ceil(AR);

		    if (AR<Math.sqrt(lower_AR*upper_AR))

		    {

		        C1 = (int)(lower_AR);

		        R1 = 1;

		    }

		    else

		    {

		        C1 = (int)(upper_AR);

		        R1 = 1;

		    }

		    if (swapAxis)

		    {

		        //std::swap(C1,R1);

		    	int t = C1;

		    	C1 = R1;

		    	R1 = t;

		    }

		    double source_xRange = ext.getMaxX() - ext.getMinX();

		    double source_yRange = ext.getMaxY() - ext.getMinY();

		    float sourceResolutionX = ((float)source_xRange)/(float)width;

		    float sourceResolutionY = ((float)source_yRange)/(float)height;

		    // log(osg::NOTICE,"Source %s resX %f resY %f",source->getFileName().c_str(), sourceResolutionX, sourceResolutionX);

		    double tileSize = 32.0;

		    int k_cols = (int)( Math.ceil( 1.0 + Math.log( destination_xRange / (C1 * sourceResolutionX * tileSize ) ) / Math.log(2.0) ) );

		    int k_rows = (int)( Math.ceil( 1.0 + Math.log( destination_yRange / (R1 * sourceResolutionY * tileSize ) ) / Math.log(2.0) ) );

		    level = Math.max(k_cols, k_rows);

		    level = Math.min(level, maxLevel);

		    level = Math.max(level, 0);

        }

		return level;

	}