Ejemplo n.º 1
0
 /*
  * (non-Javadoc)
  * @see fr.unistra.pelican.Algorithm#launch()
  */
 public void launch() {
   int xdim = input.getXDim();
   int ydim = input.getYDim();
   int zdim = input.getZDim();
   int tdim = input.getTDim();
   int bdim = input.getBDim();
   int vx = (int) (xdim * rx);
   int vy = (int) (ydim * ry);
   int vz = (int) (zdim * rz);
   int vt = (int) (tdim * rt);
   int vb = (int) (bdim * rb);
   output = ResamplingByValue.exec(input, vx, vy, vz, vt, vb, resamplingMethod);
 }
Ejemplo n.º 2
0
  /*
   * (non-Javadoc)
   *
   * @see fr.unistra.pelican.Algorithm#launch()
   */
  public void launch() throws AlgorithmException {
    int bdim = input.getBDim();
    int tdim = input.getTDim();
    int zdim = input.getZDim();
    int ydim = input.getYDim();
    int xdim = input.getXDim();

    output = new DoubleImage(xdim, ydim, zdim, tdim, bdim);
    output.copyAttributes(input);

    double kok2 = Math.sqrt(2.0);

    if (input instanceof BooleanImage)
      throw new AlgorithmException("BooleanImages are not supported");

    switch (operation) {
      case GRADY:
        for (int b = 0; b < bdim; b++) {
          for (int t = 0; t < tdim; t++) {
            for (int z = 0; z < zdim; z++) {
              for (int y = 1; y < ydim - 1; y++) {
                for (int x = 1; x < xdim - 1; x++) {
                  double tmp =
                      input.getPixelXYZTBDouble(x - 1, y - 1, z, t, b)
                          + kok2 * input.getPixelXYZTBDouble(x, y - 1, z, t, b)
                          + input.getPixelXYZTBDouble(x + 1, y - 1, z, t, b);
                  tmp =
                      input.getPixelXYZTBDouble(x - 1, y + 1, z, t, b)
                          + kok2 * input.getPixelXYZTBDouble(x, y + 1, z, t, b)
                          + input.getPixelXYZTBDouble(x + 1, y + 1, z, t, b)
                          - tmp;
                  output.setPixelXYZTBDouble(x, y, z, t, b, tmp);
                }
              }
            }
          }
        }
        break;
      case GRADX:
        for (int b = 0; b < bdim; b++) {
          for (int t = 0; t < tdim; t++) {
            for (int z = 0; z < zdim; z++) {
              for (int y = 1; y < ydim - 1; y++) {
                for (int x = 1; x < xdim - 1; x++) {
                  double tmp =
                      input.getPixelXYZTBDouble(x - 1, y - 1, z, t, b)
                          + kok2 * input.getPixelXYZTBDouble(x - 1, y, z, t, b)
                          + input.getPixelXYZTBDouble(x - 1, y + 1, z, t, b);
                  tmp =
                      input.getPixelXYZTBDouble(x + 1, y - 1, z, t, b)
                          + kok2 * input.getPixelXYZTBDouble(x + 1, y, z, t, b)
                          + input.getPixelXYZTBDouble(x + 1, y + 1, z, t, b)
                          - tmp;
                  output.setPixelXYZTBDouble(x, y, z, t, b, tmp);
                }
              }
            }
          }
        }
        break;
      case NORM:
        for (int b = 0; b < bdim; b++) {
          for (int t = 0; t < tdim; t++) {
            for (int z = 0; z < zdim; z++) {
              for (int y = 1; y < ydim - 1; y++) {
                for (int x = 1; x < xdim - 1; x++) {
                  double tmpy =
                      input.getPixelXYZTBDouble(x - 1, y - 1, z, t, b)
                          + kok2 * input.getPixelXYZTBDouble(x, y - 1, z, t, b)
                          + input.getPixelXYZTBDouble(x + 1, y - 1, z, t, b);
                  tmpy =
                      input.getPixelXYZTBDouble(x - 1, y + 1, z, t, b)
                          + kok2 * input.getPixelXYZTBDouble(x, y + 1, z, t, b)
                          + input.getPixelXYZTBDouble(x + 1, y + 1, z, t, b)
                          - tmpy;
                  double tmpx =
                      input.getPixelXYZTBDouble(x - 1, y - 1, z, t, b)
                          + kok2 * input.getPixelXYZTBDouble(x - 1, y, z, t, b)
                          + input.getPixelXYZTBDouble(x - 1, y + 1, z, t, b);
                  tmpx =
                      input.getPixelXYZTBDouble(x + 1, y - 1, z, t, b)
                          + kok2 * input.getPixelXYZTBDouble(x + 1, y, z, t, b)
                          + input.getPixelXYZTBDouble(x + 1, y + 1, z, t, b)
                          - tmpx;
                  double norm = Math.sqrt(tmpx * tmpx + tmpy * tmpy);
                  output.setPixelXYZTBDouble(x, y, z, t, b, norm);
                }
              }
            }
          }
        }
        break;
      default:
        for (int b = 0; b < bdim; b++) {
          for (int t = 0; t < tdim; t++) {
            for (int z = 0; z < zdim; z++) {
              for (int y = 1; y < ydim - 1; y++) {
                for (int x = 1; x < xdim - 1; x++) {
                  double tmpy =
                      input.getPixelXYZTBDouble(x - 1, y - 1, z, t, b)
                          + kok2 * input.getPixelXYZTBDouble(x, y - 1, z, t, b)
                          + input.getPixelXYZTBDouble(x + 1, y - 1, z, t, b);
                  tmpy =
                      input.getPixelXYZTBDouble(x - 1, y + 1, z, t, b)
                          + 2 * input.getPixelXYZTBDouble(x, y + 1, z, t, b)
                          + input.getPixelXYZTBDouble(x + 1, y + 1, z, t, b)
                          - tmpy;
                  double tmpx =
                      input.getPixelXYZTBDouble(x - 1, y - 1, z, t, b)
                          + kok2 * input.getPixelXYZTBDouble(x - 1, y, z, t, b)
                          + input.getPixelXYZTBDouble(x - 1, y + 1, z, t, b);
                  tmpx =
                      input.getPixelXYZTBDouble(x + 1, y - 1, z, t, b)
                          + 2 * input.getPixelXYZTBDouble(x + 1, y, z, t, b)
                          + input.getPixelXYZTBDouble(x + 1, y + 1, z, t, b)
                          - tmpx;
                  double norm = Math.atan(tmpy / tmpx) - 3 * Math.PI / 4;
                  output.setPixelXYZTBDouble(x, y, z, t, b, norm);
                }
              }
            }
          }
        }
    }
  }