/** * Converts a {@link InterleavedS16} into a {@link GrayS16} by computing the average value of each pixel * across all the bands. * * @param input (Input) The ImageInterleaved that is being converted. Not modified. * @param output (Optional) The single band output image. If null a new image is created. Modified. * @return Converted image. */ public static GrayS16 average( InterleavedS16 input , GrayS16 output ) { if (output == null) { output = new GrayS16(input.width, input.height); } else { output.reshape(input.width,input.height); } ConvertInterleavedToSingle.average(input, output); return output; }
@Override public Number get(int x, int y) { return image.get(x,y); }
public static void process(GrayS16 orig, GrayS16 derivX, GrayS16 derivY) { final short[] data = orig.data; final short[] imgX = derivX.data; final short[] imgY = derivY.data; final int width = orig.getWidth(); final int height = orig.getHeight() - 1; final int stride = orig.stride; for (int y = 1; y < height; y++) { int indexX = derivX.startIndex + derivX.stride * y + 1; int indexY = derivY.startIndex + derivY.stride * y + 1; int indexSrc = orig.startIndex + stride * y + 1; final int endX = indexSrc + width - 2; for (; indexSrc < endX; indexSrc++) { imgX[indexX++] = (short) (data[indexSrc + 1] - data[indexSrc - 1]); imgY[indexY++] = (short) (data[indexSrc + stride] - data[indexSrc - stride]); } } }
/** * Computes the variance of pixel intensity values inside the image. * * @param img Input image. Not modified. * @param mean Mean pixel intensity value. * @return Pixel variance */ public static double variance( GrayS16 img , double mean ) { double variance = 0; for (int y = 0; y < img.height; y++) { int index = img.getStartIndex() + y * img.getStride(); int indexEnd = index+img.width; // for(int x = 0; x < img.width; x++ ) { for (; index < indexEnd; index++ ) { double d = (img.data[index]) - mean; variance += d*d; } } return variance/(img.width*img.height); }
public static void process_I8_sub(GrayS16 orig, GrayS16 derivX, GrayS16 derivY) { final short[] data = orig.data; final short[] imgX = derivX.data; final short[] imgY = derivY.data; final int width = orig.getWidth(); final int height = orig.getHeight() - 1; final int strideSrc = orig.getStride(); for (int y = 1; y < height; y++) { int indexSrc = orig.startIndex + orig.stride * y + 1; final int endX = indexSrc + width - 2; int indexX = derivX.startIndex + derivX.stride * y + 1; int indexY = derivY.startIndex + derivY.stride * y + 1; for (; indexSrc < endX; indexSrc++) { int v = (data[indexSrc + strideSrc + 1] ) - (data[indexSrc - strideSrc - 1] ); int w = (data[indexSrc + strideSrc - 1] ) - (data[indexSrc - strideSrc + 1] ); imgY[indexY++] = (short) (((data[indexSrc + strideSrc] ) - (data[indexSrc - strideSrc] )) * 2 + v + w); imgX[indexX++] = (short) (((data[indexSrc + 1] ) - (data[indexSrc - 1] )) * 2 + v - w); } } }
/** * Bounds image pixels to be between these two values * * @param img Image * @param min minimum value. * @param max maximum value. */ public static void boundImage(GrayS16 img , int min , int max ) { final int h = img.getHeight(); final int w = img.getWidth(); short[] data = img.data; for (int y = 0; y < h; y++) { int index = img.getStartIndex() + y * img.getStride(); int indexEnd = index+w; // for(int x = 0; x < w; x++ ) { for (; index < indexEnd; index++) { int value = data[index]; if( value < min ) data[index] = (short)min; else if( value > max ) data[index] = (short)max; } } }
public static void vertical(Kernel1D_S32 kernel, GrayS16 input, GrayI16 output ) { final int offset = kernel.getOffset(); final int width = input.getWidth(); final int height = input.getHeight(); for (int y = 0; y < height; y++) { for( int x = 0; x < width; x++ ) { int total = 0; int weight = 0; int startY = y - offset; int endY = startY + kernel.getWidth(); if( startY < 0 ) startY = 0; if( endY > height ) endY = height; for( int i = startY; i < endY; i++ ) { int v = kernel.get(i-y+offset); total += input.get(x,i)*v; weight += v; } output.set(x,y, (total+weight/2)/weight ); } } }
public static void horizontal3(Kernel1D_S32 kernel , GrayS16 image, GrayI16 dest ) { final short[] dataSrc = image.data; final short[] dataDst = dest.data; final int k1 = kernel.data[0]; final int k2 = kernel.data[1]; final int k3 = kernel.data[2]; final int radius = kernel.getRadius(); final int width = image.getWidth(); for( int i = 0; i < image.height; i++ ) { int indexDst = dest.startIndex + i*dest.stride+radius; int j = image.startIndex + i*image.stride - radius; final int jEnd = j+width-radius; for( j += radius; j < jEnd; j++ ) { int indexSrc = j; int total = (dataSrc[indexSrc++])*k1; total += (dataSrc[indexSrc++])*k2; total += (dataSrc[indexSrc])*k3; dataDst[indexDst++] = ( short )total; } } }
public static void horizontal3(Kernel1D_S32 kernel , GrayS16 input, GrayI16 output , int skip ) { final short[] dataSrc = input.data; final short[] dataDst = output.data; final int k1 = kernel.data[0]; final int k2 = kernel.data[1]; final int k3 = kernel.data[2]; final int radius = kernel.getRadius(); final int widthEnd = UtilDownConvolve.computeMaxSide(input.width,skip,radius); final int height = input.getHeight(); final int offsetX = UtilDownConvolve.computeOffset(skip,radius); for( int i = 0; i < height; i++ ) { int indexDst = output.startIndex + i*output.stride + offsetX/skip; int j = input.startIndex + i*input.stride - radius; final int jEnd = j+widthEnd; for( j += offsetX; j <= jEnd; j += skip ) { int indexSrc = j; int total = (dataSrc[indexSrc++] ) * k1; total += (dataSrc[indexSrc++])*k2; total += (dataSrc[indexSrc])*k3; dataDst[indexDst++] = ( short )total; } } }
public static void convert( GrayS16 from, GrayI16 to ) { if (from.isSubimage() || to.isSubimage()) { for (int y = 0; y < from.height; y++) { int indexFrom = from.getIndex(0, y); int indexTo = to.getIndex(0, y); for (int x = 0; x < from.width; x++) { to.data[indexTo++] = ( from.data[indexFrom++] ); } } } else { final int N = from.width * from.height; System.arraycopy(from.data, 0, to.data, 0, N); } }
@Override public float unsafe_getF(int x, int y) { return image.unsafe_get(x,y); }
@Override public int unsafe_get(int x, int y) { return data[getIndex(x, y)]; }
@Override public void setImage(GrayS16 image) { this.orig = image; this.data = orig.data; this.stride = orig.getStride(); }
/** * <p> * Computes the absolute value of the difference between each pixel in the two images.<br> * d(x,y) = |img1(x,y) - img2(x,y)| * </p> * @param imgA Input image. Not modified. * @param imgB Input image. Not modified. * @param diff Absolute value of difference image. Modified. */ public static void diffAbs(GrayS16 imgA , GrayS16 imgB , GrayS16 diff ) { InputSanityCheck.checkSameShape(imgA,imgB,diff); final int h = imgA.getHeight(); final int w = imgA.getWidth(); for (int y = 0; y < h; y++) { int indexA = imgA.getStartIndex() + y * imgA.getStride(); int indexB = imgB.getStartIndex() + y * imgB.getStride(); int indexDiff = diff.getStartIndex() + y * diff.getStride(); int indexEnd = indexA+w; // for(int x = 0; x < w; x++ ) { for (; indexA < indexEnd; indexA++, indexB++, indexDiff++ ) { diff.data[indexDiff] = (short)Math.abs((imgA.data[indexA] ) - (imgB.data[indexB] )); } } }
public static void horizontal(Kernel1D_S32 kernel, GrayS16 input, GrayI16 output ) { final int offset = kernel.getOffset(); final int width = input.getWidth(); final int height = input.getHeight(); for (int y = 0; y < height; y++) { for( int x = 0; x < width; x++ ) { int total = 0; int weight = 0; int startX = x - offset; int endX = startX+kernel.getWidth(); if( startX < 0 ) startX = 0; if( endX > width ) endX = width; for( int j = startX; j < endX; j++ ) { int v = kernel.get(j-x+offset); total += input.get(j,y)*v; weight += v; } output.set(x,y, (total+weight/2)/weight ); } } }
public static void process(GrayS16 orig, GrayS16 derivX, GrayS16 derivY) { final short[] data = orig.data; final short[] imgX = derivX.data; final short[] imgY = derivY.data; final int width = orig.getWidth(); final int height = orig.getHeight() - 1; final int strideSrc = orig.getStride(); for (int y = 1; y < height; y++) { int indexSrc = orig.startIndex + orig.stride * y + 1; final int endX = indexSrc + width - 2; int indexX = derivX.startIndex + derivX.stride * y + 1; int indexY = derivY.startIndex + derivY.stride * y + 1; for (; indexSrc < endX; indexSrc++) { // a33 - a11 int w = (data[indexSrc + strideSrc + 1] ) -(data[indexSrc - strideSrc - 1] ); // a31 - a13 int v = (data[indexSrc + strideSrc - 1] ) -(data[indexSrc - strideSrc + 1] ); //a32 + w + v - a12 imgY[indexY++] = (short) ((data[indexSrc + strideSrc ] )+w+v-(data[indexSrc - strideSrc ] )); //a23 + w - v - a21 imgX[indexX++] = (short) ((data[indexSrc + 1] )+w-v-(data[indexSrc - 1] )); } } }
public static void horizontal3(Kernel1D_S32 kernel , GrayS16 image, GrayI16 dest , int divisor ) { final short[] dataSrc = image.data; final short[] dataDst = dest.data; final int k1 = kernel.data[0]; final int k2 = kernel.data[1]; final int k3 = kernel.data[2]; final int radius = kernel.getRadius(); final int width = image.getWidth(); final int halfDivisor = divisor/2; for( int i = 0; i < image.height; i++ ) { int indexDst = dest.startIndex + i*dest.stride+radius; int j = image.startIndex + i*image.stride - radius; final int jEnd = j+width-radius; for( j += radius; j < jEnd; j++ ) { int indexSrc = j; int total = (dataSrc[indexSrc++])*k1; total += (dataSrc[indexSrc++])*k2; total += (dataSrc[indexSrc])*k3; dataDst[indexDst++] = ( short )((total+halfDivisor)/divisor); } } }