Matrix.c 2.2 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778
  1. /*
  2. * The Python Imaging Library
  3. * $Id$
  4. *
  5. * colour and luminance matrix transforms
  6. *
  7. * history:
  8. * 1996-05-18 fl: created (brute force implementation)
  9. *
  10. * Copyright (c) Fredrik Lundh 1996.
  11. * Copyright (c) Secret Labs AB 1997.
  12. *
  13. * See the README file for information on usage and redistribution.
  14. */
  15. #include "Imaging.h"
  16. #define CLIPF(v) ((v <= 0.0) ? 0 : (v >= 255.0F) ? 255 : (UINT8)v)
  17. Imaging
  18. ImagingConvertMatrix(Imaging im, const char *mode, float m[]) {
  19. Imaging imOut;
  20. int x, y;
  21. ImagingSectionCookie cookie;
  22. /* Assume there's enough data in the buffer */
  23. if (!im) {
  24. return (Imaging)ImagingError_ModeError();
  25. }
  26. if (strcmp(mode, "L") == 0 && im->bands == 3) {
  27. imOut = ImagingNewDirty("L", im->xsize, im->ysize);
  28. if (!imOut) {
  29. return NULL;
  30. }
  31. ImagingSectionEnter(&cookie);
  32. for (y = 0; y < im->ysize; y++) {
  33. UINT8 *in = (UINT8 *)im->image[y];
  34. UINT8 *out = (UINT8 *)imOut->image[y];
  35. for (x = 0; x < im->xsize; x++) {
  36. float v = m[0] * in[0] + m[1] * in[1] + m[2] * in[2] + m[3] + 0.5;
  37. out[x] = CLIPF(v);
  38. in += 4;
  39. }
  40. }
  41. ImagingSectionLeave(&cookie);
  42. } else if (strlen(mode) == 3 && im->bands == 3) {
  43. imOut = ImagingNewDirty(mode, im->xsize, im->ysize);
  44. if (!imOut) {
  45. return NULL;
  46. }
  47. for (y = 0; y < im->ysize; y++) {
  48. UINT8 *in = (UINT8 *)im->image[y];
  49. UINT8 *out = (UINT8 *)imOut->image[y];
  50. ImagingSectionEnter(&cookie);
  51. for (x = 0; x < im->xsize; x++) {
  52. float v0 = m[0] * in[0] + m[1] * in[1] + m[2] * in[2] + m[3] + 0.5;
  53. float v1 = m[4] * in[0] + m[5] * in[1] + m[6] * in[2] + m[7] + 0.5;
  54. float v2 = m[8] * in[0] + m[9] * in[1] + m[10] * in[2] + m[11] + 0.5;
  55. out[0] = CLIPF(v0);
  56. out[1] = CLIPF(v1);
  57. out[2] = CLIPF(v2);
  58. in += 4;
  59. out += 4;
  60. }
  61. ImagingSectionLeave(&cookie);
  62. }
  63. } else {
  64. return (Imaging)ImagingError_ModeError();
  65. }
  66. return imOut;
  67. }