<?xml version="1.0" encoding="utf-8"?>
<!DOCTYPE html PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN"
"http://www.w3.org/TR/html4/loose.dtd">
<html>
<head>
<meta http-equiv=Content-Type content="text/html; charset=utf8">
<title>/usr/web/sources/contrib/miller/resize.c - Plan 9 from Bell Labs</title>
<!-- THIS FILE IS AUTOMATICALLY GENERATED. -->
<!-- EDIT sources.tr INSTEAD. -->
</meta>
</head>
<body>
<p style="margin-top: 0; margin-bottom: 0.17in"></p>
<p style="line-height: 1.2em; margin-left: 1.00in; text-indent: 0.00in; margin-right: 1.00in; margin-top: 0; margin-bottom: 0; text-align: center;">
<span style="font-size: 10pt"><a href="/plan9/">Plan 9 from Bell Labs</a>&rsquo;s /usr/web/sources/contrib/miller/resize.c</span></p>
<p style="margin-top: 0; margin-bottom: 0.17in"></p>
<p style="margin-top: 0; margin-bottom: 0.17in"></p>
<center><font size=-1>
Copyright © 2009 Alcatel-Lucent.<br />
Distributed under the
<a href="/plan9/license.html">Lucent Public License version 1.02</a>.
<br />
<a href="/plan9/download.html">Download the Plan 9 distribution.</a>
</font>
</center>
<p style="margin-top: 0; margin-bottom: 0.17in"></p>
<table width="100%" cellspacing=0 border=0><tr><td align="center">
<table cellspacing=0 cellpadding=5 bgcolor="#eeeeff"><tr><td align="left">
<pre>
<!-- END HEADER -->
#include &lt;u.h&gt;
#include &lt;libc.h&gt;
#include &lt;draw.h&gt;
#include &lt;memdraw.h&gt;

static void
resample(Memimage *dst, Rectangle r, Memimage *src, Rectangle sr, int bilinear)
{
       Point sp, dp;
       Point _sp, qp;
       Point ssize, dsize;
       uchar *pdst0, *pdst, *psrc0, *psrc;
       ulong s00, s01, s10, s11;
       int tx, ty, b, bpp, bpl;

       ssize = subpt(subpt(sr.max, sr.min), Pt(1,1));
       dsize = subpt(subpt(r.max, r.min), Pt(1,1));
       pdst0 = byteaddr(dst, r.min);
       bpp = src-&gt;depth/8;
       bpl = src-&gt;width*sizeof(int);

       qp.x = (ssize.x&lt;&lt;12)/dsize.x;
       qp.y = (ssize.y&lt;&lt;12)/dsize.y;
       _sp.y = sr.min.y&lt;&lt;12;
       for(dp.y=0; dp.y&lt;=dsize.y; dp.y++){
               sp.y = _sp.y&gt;&gt;12;
               ty = _sp.y&amp;0xFFF;
               if(!bilinear &amp;&amp; ty &gt;= 0x800)
                       sp.y++;
               pdst = pdst0;
               sp.x = sr.min.x;
               psrc0 = byteaddr(src, sp);
               _sp.x = 0;
               for(dp.x=0; dp.x&lt;=dsize.x; dp.x++){
                       sp.x = _sp.x&gt;&gt;12;
                       tx = _sp.x&amp;0xFFF;
                       psrc = psrc0 + sp.x*bpp;
                       if(!bilinear){
                               if(tx &gt;= 0x800)
                                       psrc += bpp;
                               for(b=0; b&lt;bpp; b++)
                                       pdst[b] = psrc[b];
                       }else{
                               s00 = (0x1000-tx)*(0x1000-ty);
                               s01 = tx*(0x1000-ty);
                               s10 = (0x1000-tx)*ty;
                               s11 = tx*ty;
                               for(b=0; b&lt;bpp; b++){
                                       pdst[b] = ( s11*psrc[bpl+bpp+b]
                                                          +s10*psrc[bpl+b]
                                                          +s01*psrc[bpp+b]
                                                          +s00*psrc[b]
                                                         ) &gt;&gt;24;
                               }
                       }
                       pdst += bpp;
                       _sp.x += qp.x;
               }
               pdst0 += dst-&gt;width*sizeof(int);
               _sp.y += qp.y;
       }
}

void
usage(void)
{
       sysfatal("Usage: %s -s width height [-b] [image]\n", argv0);
}

void
main(int argc, char **argv)
{
       Memimage *im, *nim;
       Rectangle r = {0, 0, 0, 0};
       int bilinear=  0;

       ARGBEGIN{
       case 'b':
               bilinear++;
               break;
       case 's':
               r.max.x = atoi(EARGF(usage()));
               r.max.y = atoi(EARGF(usage()));
               break;
       default:
               usage();
       }ARGEND
       if(argc != 0 &amp;&amp; argc != 1)
               usage();
       memimageinit();
       im = readmemimage(open(argv[0]? argv[0] : "/fd/0", OREAD));
       nim = allocmemimage(r, im-&gt;chan);
       resample(nim, r, im, im-&gt;r, bilinear);
       writememimage(1, nim);
       exits(0);
}
<!-- BEGIN TAIL -->
</pre>
</td></tr></table>
</td></tr></table>
<p style="margin-top: 0; margin-bottom: 0.17in"></p>
<p style="line-height: 1.2em; margin-left: 1.00in; text-indent: 0.00in; margin-right: 1.00in; margin-top: 0; margin-bottom: 0; text-align: center;">
<span style="font-size: 10pt"></span></p>
<p style="margin-top: 0; margin-bottom: 0.50in"></p>
<p style="margin-top: 0; margin-bottom: 0.33in"></p>
<center><table border="0"><tr>
<td valign="middle"><a href="http://www.alcatel-lucent.com/"><img border="0" src="/plan9/img/logo_ft.gif" alt="Bell Labs" />
</a></td>
<td valign="middle"><a href="http://www.opensource.org"><img border="0" alt="OSI certified" src="/plan9/img/osi-certified-60x50.gif" />
</a></td>
<td><img style="padding-right: 45px;" alt="Powered by Plan 9" src="/plan9/img/power36.gif" />
</td>
</tr></table></center>
<p style="margin-top: 0; margin-bottom: 0.17in"></p>
<center>
<span style="font-size: 10pt">(<a href="/plan9/">Return to Plan 9 Home Page</a>)</span>
</center>
<p style="margin-top: 0; margin-bottom: 0.17in"></p>
<center><font size=-1>
<span style="font-size: 10pt"><a href="http://www.lucent.com/copyright.html">Copyright</a></span>
<span style="font-size: 10pt">© 2009 Alcatel-Lucent.</span>
<span style="font-size: 10pt">All Rights Reserved.</span>
<br />
<span style="font-size: 10pt">Comments to</span>
<span style="font-size: 10pt"><a href="mailto:[email protected]">[email protected]</a>.</span>
</font></center>
</body>
</html>