--- ray/src/util/rcrop.c 2022/03/15 04:41:45 1.5 +++ ray/src/util/rcrop.c 2022/03/16 16:48:39 1.11 @@ -1,5 +1,5 @@ #ifndef lint -static const char RCSid[] = "$Id: rcrop.c,v 1.5 2022/03/15 04:41:45 greg Exp $"; +static const char RCSid[] = "$Id: rcrop.c,v 1.11 2022/03/16 16:48:39 greg Exp $"; #endif /* * rcrop.c - crop a Radiance picture or matrix data @@ -12,13 +12,11 @@ static const char RCSid[] = "$Id: rcrop.c,v 1.5 2022/0 #include "fvect.h" #include "view.h" -#define MAXWORD 64 /* maximum word (number) length */ - char *progname; /* global argv[0] */ VIEW vw = STDVIEW; int gotvw = 0; -char fmt[MAXFMTLEN] = "Unknown"; +char fmt[MAXFMTLEN] = "ascii"; /* assumed when unspecified */ int ncomp = 0; RESOLU res; int rmin, cmin, nrows, ncols; @@ -199,6 +197,45 @@ io_err: return(0); } +/* Adjust (crop) our view */ +static int +adjust_view(void) +{ + double p0[2], p1[2]; + const char *err; + + if (res.rt & YMAJOR) { + p0[0] = cmin/(double)res.xr; + p0[1] = rmin/(double)res.yr; + p1[0] = (cmin+ncols)/(double)res.xr; + p1[1] = (rmin+nrows)/(double)res.yr; + } else { + p0[0] = rmin/(double)res.xr; + p0[1] = cmin/(double)res.yr; + p1[0] = (rmin+nrows)/(double)res.xr; + p1[1] = (cmin+ncols)/(double)res.yr; + } + if (res.rt & XDECR) { + p0[0] = 1. - p0[0]; + p1[0] = 1. - p1[0]; + } + if (res.rt & YDECR) { + p0[1] = 1. - p0[1]; + p1[1] = 1. - p1[1]; + } + err = cropview(&vw, p0[0], p0[1], p1[0], p1[1]); + + if (!err) + return(1); /* success! */ + + fputs(progname, stderr); + fputs(": view error - ", stderr); + fputs(err, stderr); + fputc('\n', stderr); + return(0); /* something went wrong */ +} + + /* Main routine -- load header and call processor */ int main(int argc, char *argv[]) @@ -234,6 +271,10 @@ main(int argc, char *argv[]) fputs(": cannot open for writing\n", stderr); return(1); } +#ifdef getc_unlocked /* avoid stupid semaphores */ + flockfile(fp); + flockfile(stdout); +#endif /* process information header */ if (getheader(fp, headline, NULL) < 0) { fputs(progname, stderr); @@ -257,43 +298,13 @@ main(int argc, char *argv[]) fputs(": illegal crop\n", stderr); return(1); } - printargs(argc, argv, stdout); - if (gotvw) { /* adjust view? */ - double p0[2], p1[2]; - const char *err; - if (res.rt & YMAJOR) { - p0[0] = cmin/(double)res.xr; - p0[1] = rmin/(double)res.yr; - p1[0] = (cmin+ncols)/(double)res.xr; - p1[1] = (rmin+nrows)/(double)res.yr; - } else { - p0[1] = cmin/(double)res.xr; - p0[0] = rmin/(double)res.yr; - p1[1] = (cmin+ncols)/(double)res.xr; - p1[0] = (rmin+nrows)/(double)res.yr; - } - if (res.rt & XDECR) { - p0[0] = 1. - p0[0]; - p1[0] = 1. - p1[0]; - } - if (res.rt & YDECR) { - p0[1] = 1. - p0[1]; - p1[1] = 1. - p1[1]; - } - err = cropview(&vw, p0[0], p0[1], p1[0], p1[1]); - if (err) { - fputs(progname, stderr); - fputs(": view error - ", stderr); - fputs(err, stderr); - fputc('\n', stderr); - return(1); - } else { - fputs(VIEWSTR, stdout); - fprintview(&vw, stdout); - fputc('\n', stdout); - } + printargs(5, argv, stdout); /* add to header */ + if (gotvw && adjust_view()) { + fputs(VIEWSTR, stdout); /* write adjusted view */ + fprintview(&vw, stdout); + fputc('\n', stdout); } - if (gotdims) + if (gotdims) /* dimensions + format */ printf("NROWS=%d\nNCOLS=%d\n", nrows, ncols); if (ncomp) printf("NCOMP=%d\n", ncomp); @@ -339,8 +350,13 @@ main(int argc, char *argv[]) return(1); } if (!(asiz < 0 ? colr_copyf(fp) : - !asiz ? ascii_copyf(fp) : binary_copyf(fp, asiz))) + asiz ? binary_copyf(fp, asiz) : ascii_copyf(fp))) return(1); + /* need to consume the rest? */ + if (fp == stdin && rmin+nrows < numscans(&res) && + fseek(fp, 0L, SEEK_END) < 0) + while (getc(fp) != EOF) + ; return(0); usage: fputs("Usage: ", stderr);