Method doProjection

ij/src/main/java/ij/plugin/ZProjector.java:306–385

Performs actual projection using specified method.

()

    /** Performs actual projection using specified method. */
	@AstroImageJ(reason = "Merge FITS headers as well; Force makeOutputImage to always return a FloatProcessor", modified = true)
    public void doProjection() {
		if (imp==null)
			return;
		if (method<AVG_METHOD || method>MEDIAN_METHOD)
			method = AVG_METHOD;
		if (imp.getBitDepth()==24) {
			doRGBProjection();
			return;
		}
		if (imp.getBitDepth()==32 && method==AVG_METHOD) {
			projImage = doAverageFloatProjection(imp);
			return;
		}
		sliceCount = 0;
    	for (int slice=startSlice; slice<=stopSlice; slice+=increment)
    		sliceCount++;
		if (method==MEDIAN_METHOD) {
			projImage = doMedianProjection();
			// Merge Fits headers
			mergeFitsHeaders(method, startSlice, stopSlice, increment, imp, projImage);
			return;
		} 
		
		// Create new float processor for projected pixels.
		FloatProcessor fp = new FloatProcessor(imp.getWidth(),imp.getHeight()); 
		ImageStack stack = imp.getStack();
		RayFunction rayFunc = getRayFunction(method, fp);
		if (IJ.debugMode==true) {
	    	IJ.log("\nProjecting stack from: "+startSlice
		     	+" to: "+stopSlice); 
		}

		// Determine type of input image. Explicit determination of
		// processor type is required for subsequent pixel
		// manipulation.  This approach is more efficient than the
		// more general use of ImageProcessor's getPixelValue and
		// putPixel methods.
		int ptype; 
		if (stack.getProcessor(1) instanceof ByteProcessor) ptype = BYTE_TYPE; 
		else if (stack.getProcessor(1) instanceof ShortProcessor) ptype = SHORT_TYPE; 
		else if (stack.getProcessor(1) instanceof FloatProcessor) ptype = FLOAT_TYPE; 
		else {
	    	IJ.error("Z Project", "Non-RGB stack required"); 
	    	return; 
		}

		// Do the projection
		int sliceCount = 0;
		for (int n=startSlice; n<=stopSlice; n+=increment) {
			if (!isHyperstack) {
	    		IJ.showStatus("ZProjection " + color +": " + n + "/" + stopSlice);
	    		IJ.showProgress(n-startSlice, stopSlice-startSlice);
	    	}
	    	projectSlice(stack.getPixels(n), rayFunc, ptype);
	    	sliceCount++;
		}