Changeset 1613 for trunk/psLib/src/image/psImageExtraction.c

Timestamp:

Aug 24, 2004, 2:04:01 PM (22 years ago)

Author:

desonia

Message:

added psImageTrim and psImageSubsection.

File:

• trunk/psLib/src/image/psImageExtraction.c (modified) (2 diffs)

trunk/psLib/src/image/psImageExtraction.c

@@ -10 +10 @@
 *  @author Robert DeSonia, MHPCC
 *
-*  @version $Revision: 1.11 $ $Name: not supported by cvs2svn $
-*  @date $Date: 2004-08-23 22:36:03 $
+*  @version $Revision: 1.12 $ $Name: not supported by cvs2svn $
+*  @date $Date: 2004-08-25 00:04:01 $
 *
 *  Copyright 2004 Maui High Performance Computing Center, University of Hawaii
@@ -93 +93 @@
 
     return (out);
+}
+
+psImage* psImageSubsection(psImage* image, const char* section)
+{
+    int x1;
+    int x2;
+    int y1;
+    int y2;
+
+    // section should be of the form '[x1:x2,y1:y2]'
+    if (section == NULL) {
+        psError(__func__,"The subsection string input can not be NULL.");
+        return NULL;
+    }
+
+    if (sscanf(section,"[%d:%d,%d:%d]",&x1,&x2,&y1,&y2) < 4) {
+        psError(__func__,"The subsection string (%s) can not be parsed.  "
+                "Needs to be of the form '[x1:x2,y1:y2]'",
+                section);
+        return NULL;
+    }
+    return psImageSubset(image,x2-x1+1,y2-y1+1,x1,y1);
+}
+
+psImage* psImageTrim(psImage* image, int x0, int x1, int y0, int y1)
+{
+    if (image == NULL || image->data.V == NULL) {
+        psError(__func__, "Can not subset image because input image or its "
+                "pixel buffer is NULL.");
+        return NULL;
+    }
+
+    if (image->parent != NULL) {
+        psError(__func__, "Can not perform a trim on a child image.");
+        return NULL;
+    }
+
+    if (x0 < 0 || x1 < 0 || y0 < 0 || y1 < 0 ||
+            x0 >= image->numCols || x1 >= image->numCols ||
+            y0 >= image->numRows || y1 >= image->numRows ||
+            x0 > x1 || y0 > y1 ) {
+        psError(__func__, "Can not subset image because specified region "
+                "[%d:%d,%d:%d] is not valid for image of size %dx%d.",
+                x0,x1,y0,y1,image->numCols,image->numRows);
+        return NULL;
+    }
+
+    psImageFreeChildren(image);
+
+    unsigned int elementSize = PSELEMTYPE_SIZEOF(image->type.type);
+    unsigned int numCols = x1-x0+1;
+    unsigned int numRows = y1-y0+1;
+    unsigned int rowSize = elementSize*numCols;
+    unsigned int colOffset = elementSize * x0;
+    void* imageData = image->rawDataBuffer;
+    for (int row = y0; row < y1; row++) {
+        memmove(imageData,image->data.U8[row] + colOffset,rowSize);
+        imageData = (psU8*)imageData + rowSize;
+    }
+
+    *(unsigned int*)&image->numRows = numRows;
+    *(unsigned int*)&image->numCols = numCols;
+
+    // XXX: should I really resize the buffers?
+    image->data.V = psRealloc(image->data.V,sizeof(void*)*numRows);
+    image->rawDataBuffer = psRealloc(image->rawDataBuffer,rowSize*numRows);
+
+    image->data.V[0] = image->rawDataBuffer;
+    for (int r = 1; r < numRows; r++) {
+        image->data.U8[r] = image->data.U8[r-1] + rowSize;
+    }
+
+    return (image);
 }