void operator() (T const& source)
{
using image_type = T;
//source and target image data types must match
if (target_raster_.data_.template is<image_type>())
{
image_type & target = util::get<image_type>(target_raster_.data_);
warp_image (target, source, prj_trans_, target_raster_.ext_, source_ext_,
offset_x_, offset_y_, mesh_size_, scaling_method_, filter_factor_);
}
}
int
main(int argc, char **argv)
{
char *picAname, *picBname, *linesfilename;
FILE *picA, *picB, *linesfile;
size_t pa_width = 0, pa_height = 0;
int dissolvefrac;
double warpfrac;
unsigned char *pa, *pb, *wa, *wb, *morph;
double a, b, p;
long int numlines;
struct lineseg *lines;
long int i;
long int autosize;
autosize = 1L;
pa_width = pa_height = 0;
if (get_args(argc, argv, &picAname, &picBname, &linesfilename,
&warpfrac, &dissolvefrac, &autosize, &pa_width, &pa_height) == 0
|| isatty(fileno(stdout))) {
fprintf(stderr,
"Usage: pixmorph [-w width] [-n height] picA.pix picB.pix linesfile warpfrac dissolvefrac > out.pix\n");
return 1;
}
picA = fopen(picAname, "r");
if (picA == NULL) {
fprintf(stderr, "pixmorph: cannot open %s\n", picAname);
return 1;
}
picB = fopen(picBname, "r");
if (picB == NULL) {
fprintf(stderr, "pixmorph: cannot open %s\n", picBname);
return 1;
}
linesfile = fopen(linesfilename, "r");
if (linesfile == NULL) {
fprintf(stderr, "pixmorph: cannot open %s\n", linesfilename);
return 1;
}
if (warpfrac < 0.0 || warpfrac > 1.0) {
fprintf(stderr, "pixmorph: warpfrac must be between 0 and 1\n");
return 1;
}
if (dissolvefrac < 0 || dissolvefrac > 255) {
fprintf(stderr, "pixmorph: dissolvefrac must be between 0 and 1\n");
return 1;
}
if (autosize) {
if (fb_common_file_size(&pa_width, &pa_height, argv[1], 3) == 0) {
fprintf(stderr, "pixmorph: unable to autosize\n");
return 1;
}
} else {
struct stat sb;
if (stat(picAname, &sb) < 0) {
perror("pixmorph: unable to stat file:");
return 1;
}
if (pa_width > 0) {
pa_height = sb.st_size/(3*pa_width);
fprintf(stderr, "width = %lu, size = %ld, so height = %lu\n",
(unsigned long)pa_width, (long)sb.st_size, (unsigned long)pa_height);
} else if (pa_height > 0) pa_width = sb.st_size/(3*pa_height);
if (pa_width <= 0 || pa_height <= 0) {
fprintf(stderr, "pixmorph: Bogus image dimensions: %lu %lu\n",
(unsigned long)pa_width, (unsigned long)pa_height);
return 1;
}
}
/* Allocate memory for our bag o' pixels. */
pa = (unsigned char *)malloc(pa_width*pa_height*3);
pb = (unsigned char *)malloc(pa_width*pa_height*3);
wa = (unsigned char *)malloc(pa_width*pa_height*3);
wb = (unsigned char *)malloc(pa_width*pa_height*3);
morph = (unsigned char *)malloc(pa_width*pa_height*3);
if (pa == NULL || pb == NULL || wa == NULL || wb == NULL ||
morph == NULL) {
fprintf(stderr, "pixmorph: memory allocation failure\n");
bu_free(pa, "pa alloc from malloc");
bu_free(pb, "pb alloc from malloc");
bu_free(wa, "wa alloc from malloc");
bu_free(wb, "wb alloc from malloc");
bu_free(morph, "morph alloc from malloc");
return 1;
}
/* The following is our memorizing table for weight calculation. */
for (i = 0; i < MAXLEN; i++)
weightlookup[i] = -1.0;
fprintf(stderr, "pixmorph: Reading images and lines file.\n");
if (pix_readpixels(picA, pa_width*pa_height, pa) < pa_width*pa_height) {
fprintf(stderr, "Error reading %lu pixels from %s\n",
(unsigned long)pa_width*pa_height, picAname);
return 1;
}
if (pix_readpixels(picB, pa_width*pa_height, pb) < pa_width*pa_height) {
fprintf(stderr, "Error reading %lu pixels from %s\n",
(unsigned long)pa_width*pa_height, picBname);
return 1;
}
fclose(picA);
fclose(picB);
/* Process the lines file. */
lines_headerinfo(linesfile, &a, &b, &p, &numlines);
lines = (struct lineseg *)malloc(numlines * sizeof(struct lineseg));
numlines = lines_read(linesfile, numlines, lines,
pa_width, pa_height, warpfrac, p*b);
fprintf(stderr, "pixmorph: %ld line segments read\n", numlines);
/* Warp the images */
fprintf(stderr,
"pixmorph: Warping first image into first intermediate image.\n");
warp_image(wa, pa, lines, FIRST, pa_width, pa_height, numlines, a, b);
fprintf(stderr,
"pixmorph: Warping second image into second intermediate image.\n");
warp_image(wb, pb, lines, LAST, pa_width, pa_height, numlines, a, b);
/* Do the dissolve */
fprintf(stderr,
"pixmorph: Performing cross-dissolve between first and second\n");
fprintf(stderr, "pixmorph: intermediate images.\n");
cross_dissolve(morph, wa, wb, dissolvefrac, pa_width*pa_height);
/* All done. Write everything out to a file and quit. */
pix_writepixels(pa_width*pa_height, morph);
fprintf(stderr, "pixmorph: Morph complete.\n");
/* System takes care of memory deallocation. */
return 0;
}
