| 3cf269 | 1998-03-11 | Mirar (Pontus Hagland) | |
|
| 286fef | 1998-02-15 | Henrik Wallin | |
|
| 3cf269 | 1998-03-11 | Mirar (Pontus Hagland) | | **! $Id: orient.c,v 1.5 1998/03/11 20:44:50 mirar Exp $
|
| 286fef | 1998-02-15 | Henrik Wallin | | **! class image
*/
#include "global.h"
#include <math.h>
#include <ctype.h>
#include "stralloc.h"
#include "global.h"
#include "pike_macros.h"
#include "object.h"
#include "constants.h"
#include "interpret.h"
#include "svalue.h"
|
| 3cf269 | 1998-03-11 | Mirar (Pontus Hagland) | | #include "threads.h"
|
| 286fef | 1998-02-15 | Henrik Wallin | | #include "array.h"
#include "error.h"
#include "image.h"
#include <builtin_functions.h>
extern struct program *image_program;
#ifdef THIS
#undef THIS /* Needed for NT */
#endif
#define THIS ((struct image *)(fp->current_storage))
#define THISOBJ (fp->current_object)
#define testrange(x) MAXIMUM(MINIMUM((x),255),0)
static const double c0=0.70710678118654752440;
|
| d0224e | 1998-02-15 | Mirar (Pontus Hagland) | | static const double my_PI=3.14159265358979323846;
|
| 286fef | 1998-02-15 | Henrik Wallin | |
|
| 3cf269 | 1998-03-11 | Mirar (Pontus Hagland) | | #if 1
#include <sys/resource.h>
#define CHRONO(X) chrono(X);
static void chrono(char *x)
{
struct rusage r;
static struct rusage rold;
getrusage(RUSAGE_SELF,&r);
fprintf(stderr,"%s: %ld.%06ld - %ld.%06ld\n",x,
(long)r.ru_utime.tv_sec,(long)r.ru_utime.tv_usec,
(long)(((r.ru_utime.tv_usec-rold.ru_utime.tv_usec<0)?-1:0)
+r.ru_utime.tv_sec-rold.ru_utime.tv_sec),
(long)(((r.ru_utime.tv_usec-rold.ru_utime.tv_usec<0)?1000000:0)
+ r.ru_utime.tv_usec-rold.ru_utime.tv_usec)
);
rold=r;
}
#else
#define CHRONO(X)
#endif
|
| 286fef | 1998-02-15 | Henrik Wallin | | |
| 3cf269 | 1998-03-11 | Mirar (Pontus Hagland) | | **! method object orient()
**! method array(object) orient4()
**! Draws images describing the orientation
**! of the current image.
|
| 286fef | 1998-02-15 | Henrik Wallin | | **!
|
| 3cf269 | 1998-03-11 | Mirar (Pontus Hagland) | | **! <tt>orient</tt> gives a HSV image
**! (run a <ref>hsv_to_rgb</ref> pass on it
**! to get a viewable image).
**! corresponding to the angle of the
**! orientation:
**! <pre> | / - \
**! hue= 0 64 128 192 (=red in a hsv image)
**! purple cyan green red
**! </pre>
**! Red, green and blue channels are added
**! and not compared separately.
|
| 286fef | 1998-02-15 | Henrik Wallin | | **!
|
| 3cf269 | 1998-03-11 | Mirar (Pontus Hagland) | | **! The <tt>orient4</tt> function gives back
**! 4 image objects, corresponding to the
**! amount of different directions, see above.
|
| 286fef | 1998-02-15 | Henrik Wallin | | **!
|
| 3cf269 | 1998-03-11 | Mirar (Pontus Hagland) | | **! returns an array of the five new image objects
|
| 286fef | 1998-02-15 | Henrik Wallin | | **!
|
| 3cf269 | 1998-03-11 | Mirar (Pontus Hagland) | | **! note
**! experimental status; may not be exact the same
**! output in later versions
|
| 286fef | 1998-02-15 | Henrik Wallin | | */
static INLINE int sq(int a) { return a*a; }
|
| 3cf269 | 1998-03-11 | Mirar (Pontus Hagland) | | static INLINE int abs(int a) { return (a<0)?-a:a; }
static void _image_orient(struct image *source,
struct object *o[4],
struct image *img[4])
{
int i;
struct { int x,y; } or[4]={ {1,0}, {1,1}, {0,1}, {-1,1} };
int x,y;
for (i=0; i<5; i++)
{
push_int(source->xsize);
push_int(source->ysize);
o[i]=clone_object(image_program,2);
img[i]=(struct image*)get_storage(o[i],image_program);
push_object(o[i]);
}
THREADS_ALLOW();
CHRONO("start");
for (i=0; i<4; i++)
{
rgb_group *d=img[i]->img;
rgb_group *s=source->img;
int xz=source->xsize;
int yz=source->ysize;
int xd=or[i].x;
int yd=or[i].y;
for(x=1; x<xz-1; x++)
for(y=1; y<yz-1; y++)
{
#define FOOBAR(CO) \
d[x+y*xz].CO \
= \
(COLORTYPE) \
abs( s[(x+xd)+(y+yd)*xz].CO - s[(x-xd)+(y-yd)*xz].CO )
FOOBAR(r);
FOOBAR(g);
FOOBAR(b);
#undef FOOBAR
}
}
CHRONO("end");
THREADS_DISALLOW();
}
|
| 286fef | 1998-02-15 | Henrik Wallin | |
void image_orient(INT32 args)
{
|
| 3cf269 | 1998-03-11 | Mirar (Pontus Hagland) | | struct object *o[5];
struct image *img[5],*this;
int n;
rgb_group *d,*s1,*s2,*s3,*s0;
float mag;
|
| 286fef | 1998-02-15 | Henrik Wallin | |
if (!THIS->img) { error("no image\n"); return; }
|
| 3cf269 | 1998-03-11 | Mirar (Pontus Hagland) | |
this=THIS;
if (args)
{
if (sp[-args].type==T_INT)
mag=sp[-args].u.integer;
else if (sp[-args].type==T_FLOAT)
mag=sp[-args].u.float_number;
else
error("Illegal argument 1 to image->orient\n");
pop_n_elems(args);
}
else mag=1.0;
|
| 286fef | 1998-02-15 | Henrik Wallin | |
|
| 3cf269 | 1998-03-11 | Mirar (Pontus Hagland) | | _image_orient(this,o,img);
s0=img[0]->img;
s1=img[1]->img;
s2=img[2]->img;
s3=img[3]->img;
d=img[4]->img;
THREADS_ALLOW();
CHRONO("begin hsv...");
n=this->xsize*this->ysize;
while (n--)
{
float j=(s0->r+s0->g+s0->b-s2->r-s2->g-s2->b)/3.0;
float h=(s1->r+s1->g+s1->b-s3->r-s3->g-s3->b)/3.0;
int z,w;
if (abs(h)>abs(j))
if (h) z=-(int)(32*(j/h)+(h>0)*128+64),w=abs(h);
else z=0,w=0;
else
z=-(int)(-32*(h/j)+(j>0)*128+128),w=abs(j);
d->r=(COLORTYPE)z;
d->g=255;
d->b=(COLORTYPE)MINIMUM(w*mag,255);
d++;
s0++;
s1++;
s2++;
s3++;
}
CHRONO("end hsv...");
THREADS_DISALLOW();
o[4]->refs++;
pop_n_elems(5);
push_object(o[4]);
}
|
| 286fef | 1998-02-15 | Henrik Wallin | |
|
| 3cf269 | 1998-03-11 | Mirar (Pontus Hagland) | |
void image_orient4(INT32 args)
{
struct object *o[5];
struct image *img[5];
|
| 286fef | 1998-02-15 | Henrik Wallin | |
|
| 3cf269 | 1998-03-11 | Mirar (Pontus Hagland) | | if (!THIS->img) { error("no image\n"); return; }
pop_n_elems(args);
_image_orient(THIS,o,img);
pop_n_elems(1);
f_aggregate(4);
|
| 286fef | 1998-02-15 | Henrik Wallin | | }
|
| 3cf269 | 1998-03-11 | Mirar (Pontus Hagland) | |
|