|
version 1.1, 2009/06/18 20:12:09
|
version 1.2, 2011/08/16 02:29:10
|
| Line 328 int imageinterp( |
|
| Line 328 int imageinterp( |
|
| | |
| int xpixels, /* x width of input array */ | int xpixels, /* x width of input array */ |
| int ypixels, /* y width of input array */ | int ypixels, /* y width of input array */ |
| |
double x0, /* x pixel address of disk center */ |
| |
double y0, /* y pixel address of disk center */ |
| |
double P, /* angle between CCD y-axis and solar vertical */ |
| |
/* positive to the east (CCW) */ |
| double rsun, /* pixels */ | double rsun, /* pixels */ |
| | |
| int interpolation, /* option */ | int interpolation, /* option */ |
| Line 346 int imageinterp( |
|
| Line 350 int imageinterp( |
|
| long i; | long i; |
| double *vdp; | double *vdp; |
| float *vp; | float *vp; |
| |
double xtmp, ytmp; |
| | |
| double *VD; | double *VD; |
| VD=(double *)(malloc(xpixels*ypixels*sizeof(double))); | VD=(double *)(malloc(xpixels*ypixels*sizeof(double))); |
| Line 373 int imageinterp( |
|
| Line 378 int imageinterp( |
|
| | |
| for (col = 0; col < cols; col++) { | for (col = 0; col < cols; col++) { |
| | |
| x=col*xratio; |
xtmp=col*xratio - x0; |
| y=row*yratio; |
ytmp=row*yratio - y0; |
| |
x= x0 + xtmp*cos(P) + ytmp*sin(P); |
| |
y= y0 - xtmp*sin(P) + ytmp*cos(P); |
| | |
| Distort(x, y, rsun, +1, &x, &y, distpars); | Distort(x, y, rsun, +1, &x, &y, distpars); |
| | |
Return to
imageinterp.c
CVS log
Up to [Development] / JSOC / proj / globalhs / apps