[FFmpeg-devel] [PATCH] Add FITS Decoder

Paras Chadha paraschadha18 at gmail.com
Tue Jun 20 17:12:39 EEST 2017


On Tue, Jun 20, 2017 at 2:44 PM, Paul B Mahol <onemda at gmail.com> wrote:

> On 6/19/17, Paras Chadha <paraschadha18 at gmail.com> wrote:
> > Added support for XTENSION keyword
> >
> > Signed-off-by: Paras Chadha <paraschadha18 at gmail.com>
> > ---
> >  Changelog               |   1 +
> >  doc/general.texi        |   2 +
> >  libavcodec/Makefile     |   1 +
> >  libavcodec/allcodecs.c  |   1 +
> >  libavcodec/avcodec.h    |   1 +
> >  libavcodec/codec_desc.c |   8 +
> >  libavcodec/fitsdec.c    | 518
> > ++++++++++++++++++++++++++++++++++++++++++++++++
> >  libavcodec/version.h    |   2 +-
> >  libavformat/img2.c      |   1 +
> >  9 files changed, 534 insertions(+), 1 deletion(-)
> >  create mode 100644 libavcodec/fitsdec.c
> >
> > diff --git a/Changelog b/Changelog
> > index a893efa..a07644c 100644
> > --- a/Changelog
> > +++ b/Changelog
> > @@ -22,6 +22,7 @@ version <next>:
> >  - headphone audio filter
> >  - superequalizer audio filter
> >  - roberts video filter
> > +- FITS decoder
> >
> >  version 3.3:
> >  - CrystalHD decoder moved to new decode API
> > diff --git a/doc/general.texi b/doc/general.texi
> > index 8f582d5..c00ce32 100644
> > --- a/doc/general.texi
> > +++ b/doc/general.texi
> > @@ -591,6 +591,8 @@ following image formats are supported:
> >      @tab Digital Picture Exchange
> >  @item EXR          @tab   @tab X
> >      @tab OpenEXR
> > + at item FITS         @tab   @tab X
> > +    @tab Flexible Image Transport System
> >  @item JPEG         @tab X @tab X
> >      @tab Progressive JPEG is not supported.
> >  @item JPEG 2000    @tab X @tab X
> > diff --git a/libavcodec/Makefile b/libavcodec/Makefile
> > index 2e7f19d..2cfd1fe 100644
> > --- a/libavcodec/Makefile
> > +++ b/libavcodec/Makefile
> > @@ -291,6 +291,7 @@ OBJS-$(CONFIG_FFV1_DECODER)            += ffv1dec.o
> > ffv1.o
> >  OBJS-$(CONFIG_FFV1_ENCODER)            += ffv1enc.o ffv1.o
> >  OBJS-$(CONFIG_FFWAVESYNTH_DECODER)     += ffwavesynth.o
> >  OBJS-$(CONFIG_FIC_DECODER)             += fic.o
> > +OBJS-$(CONFIG_FITS_DECODER)            += fitsdec.o
> >  OBJS-$(CONFIG_FLAC_DECODER)            += flacdec.o flacdata.o flac.o
> >  OBJS-$(CONFIG_FLAC_ENCODER)            += flacenc.o flacdata.o flac.o
> > vorbis_data.o
> >  OBJS-$(CONFIG_FLASHSV_DECODER)         += flashsv.o
> > diff --git a/libavcodec/allcodecs.c b/libavcodec/allcodecs.c
> > index 27110e1..8a89264 100644
> > --- a/libavcodec/allcodecs.c
> > +++ b/libavcodec/allcodecs.c
> > @@ -186,6 +186,7 @@ static void register_all(void)
> >      REGISTER_ENCDEC (FFV1,              ffv1);
> >      REGISTER_ENCDEC (FFVHUFF,           ffvhuff);
> >      REGISTER_DECODER(FIC,               fic);
> > +    REGISTER_DECODER(FITS,              fits);
> >      REGISTER_ENCDEC (FLASHSV,           flashsv);
> >      REGISTER_ENCDEC (FLASHSV2,          flashsv2);
> >      REGISTER_DECODER(FLIC,              flic);
> > diff --git a/libavcodec/avcodec.h b/libavcodec/avcodec.h
> > index 39be8cf..1310426 100644
> > --- a/libavcodec/avcodec.h
> > +++ b/libavcodec/avcodec.h
> > @@ -447,6 +447,7 @@ enum AVCodecID {
> >      AV_CODEC_ID_SRGC,
> >      AV_CODEC_ID_SVG,
> >      AV_CODEC_ID_GDV,
> > +    AV_CODEC_ID_FITS,
> >
> >      /* various PCM "codecs" */
> >      AV_CODEC_ID_FIRST_AUDIO = 0x10000,     ///< A dummy id pointing at
> the
> > start of audio codecs
> > diff --git a/libavcodec/codec_desc.c b/libavcodec/codec_desc.c
> > index cf1246e..0112517 100644
> > --- a/libavcodec/codec_desc.c
> > +++ b/libavcodec/codec_desc.c
> > @@ -1464,6 +1464,14 @@ static const AVCodecDescriptor
> codec_descriptors[] =
> > {
> >                       AV_CODEC_PROP_LOSSLESS,
> >      },
> >      {
> > +        .id        = AV_CODEC_ID_FITS,
> > +        .type      = AVMEDIA_TYPE_VIDEO,
> > +        .name      = "fits",
> > +        .long_name = NULL_IF_CONFIG_SMALL("Flexible Image Transport
> > System"),
> > +        .props     = AV_CODEC_PROP_INTRA_ONLY | AV_CODEC_PROP_LOSSY |
> > +                     AV_CODEC_PROP_LOSSLESS,
> > +    },
> > +    {
> >          .id        = AV_CODEC_ID_GIF,
> >          .type      = AVMEDIA_TYPE_VIDEO,
> >          .name      = "gif",
> > diff --git a/libavcodec/fitsdec.c b/libavcodec/fitsdec.c
> > new file mode 100644
> > index 0000000..4e80936
> > --- /dev/null
> > +++ b/libavcodec/fitsdec.c
> > @@ -0,0 +1,518 @@
> > +/*
> > + * FITS image decoder
> > + * Copyright (c) 2017 Paras Chadha
> > + *
> > + * This file is part of FFmpeg.
> > + *
> > + * FFmpeg is free software; you can redistribute it and/or
> > + * modify it under the terms of the GNU Lesser General Public
> > + * License as published by the Free Software Foundation; either
> > + * version 2.1 of the License, or (at your option) any later version.
> > + *
> > + * FFmpeg is distributed in the hope that it will be useful,
> > + * but WITHOUT ANY WARRANTY; without even the implied warranty of
> > + * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
> > + * Lesser General Public License for more details.
> > + *
> > + * You should have received a copy of the GNU Lesser General Public
> > + * License along with FFmpeg; if not, write to the Free Software
> > + * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA
> 02110-1301
> > USA
> > + */
> > +
> > +/**
> > + * @file
> > + * FITS image decoder
> > + * It supports all 2-d images alongwith, bzero, bscale and blank
> keywords.
> > + * RGBA images are supported as NAXIS3 = 3 or 4 i.e. Planes in RGBA
> order.
> > Also CTYPE = 'RGB ' should be present.
> > + * Also to interpret data, values are linearly scaled using min-max
> scaling
> > but not RGB images.
> > + */
> > +
> > +#include "avcodec.h"
> > +#include "internal.h"
> > +#include <float.h>
> > +#include "libavutil/intreadwrite.h"
> > +
> > +/**
> > + * Structure to store the header keywords in FITS file
> > + */
> > +typedef struct FITSContext {
> > +    char simple;
> > +    int bitpix;
> > +    int blank;
> > +    int naxis;
> > +    int naxisn[3];
> > +    int rgb; /**< 1 if file contains RGB image, 0 otherwise */
> > +    int xtension;
> > +    double bscale;
> > +    double bzero;
> > +    double data_min;
> > +    double data_max;
> > +} FITSDecContext;
> > +
> > +/**
> > + * function calculates the data_min and data_max values from the data.
> > + * This is called if the values are not present in the header.
> > + * @param ptr8 - pointer to the data
> > + * @param header - pointer to the header
> > + * @return 1, if calculated successfully, otherwise AVERROR_INVALIDDATA
> > + */
> > +static int fill_data_min_max(const uint8_t * ptr8, FITSDecContext *
> header,
> > const uint8_t * end)
> > +{
> > +    int16_t t16;
> > +    int32_t t32;
> > +    int64_t t64;
>
> Shouldn't these be unsigned?
>

No, FITS standard states that 16, 32 and 64 bit integers are signed while 8
bit integers are unsigned.


>
> > +    float tflt;
> > +    double tdbl;
> > +    int i, j;
> > +
> > +    header->data_min = DBL_MAX;
> > +    header->data_max = DBL_MIN;
> > +    switch (header->bitpix) {
> > +        case -64:
> > +            t64 = AV_RB64(ptr8);
> > +            memcpy(&tdbl, &t64, 8);
>
> Try to use av_int2double() from libavutil/intfloat.h
>
> > +            for (i = 0; i < header->naxisn[1]; i++) {
> > +                for (j = 0; j < header->naxisn[0]; j++) {
> > +                    t64 = AV_RB64(ptr8);
> > +                    memcpy(&tdbl, &t64, 8);
>
> ditto.
>
> > +                    if (tdbl > header->data_max)
> > +                        header->data_max = tdbl;
> > +                    if (tdbl < header->data_min)
> > +                        header->data_min = tdbl;
> > +                    ptr8 += 8;
> > +                }
> > +            }
> > +            break;
> > +        case -32:
> > +            t32 = AV_RB32(ptr8);
> > +            memcpy(&tflt, &t32, 4);
>
> Try to use av_int2float() from libavutil/intfloat.h
>
> > +            for (i = 0; i < header->naxisn[1]; i++) {
> > +                for (j = 0; j < header->naxisn[0]; j++) {
> > +                    t32 = AV_RB32(ptr8);
> > +                    memcpy(&tflt, &t32, 4);
>
> ditto.
>
> > +                    if (tflt > header->data_max)
> > +                        header->data_max = tflt;
> > +                    if (tflt < header->data_min)
> > +                        header->data_min = tflt;
> > +                    ptr8 += 4;
> > +                }
> > +            }
> > +            break;
> > +        case 8:
> > +            for (i = 0; i < header->naxisn[1]; i++) {
> > +                for (j = 0; j < header->naxisn[0]; j++) {
> > +                    if (ptr8[0] != header->blank) {
> > +                        if (ptr8[0] > header->data_max)
> > +                            header->data_max = ptr8[0];
> > +                        if (ptr8[0] < header->data_min)
> > +                            header->data_min = ptr8[0];
> > +                    }
> > +                    ptr8++;
> > +                }
> > +            }
> > +            break;
> > +        case 16:
> > +            t16 = AV_RB16(ptr8);
> > +            for (i = 0; i < header->naxisn[1]; i++) {
> > +                for (j = 0; j < header->naxisn[0]; j++) {
> > +                    t16 = AV_RB16(ptr8);
> > +                    if (t16 != header->blank) {
> > +                        if (t16 > header->data_max)
> > +                            header->data_max = t16;
> > +                        if (t16 < header->data_min)
> > +                            header->data_min = t16;
> > +                    }
> > +                    ptr8 += 2;
> > +                }
> > +            }
> > +            break;
> > +        case 32:
> > +            t32 = AV_RB32(ptr8);
> > +            for (i = 0; i < header->naxisn[1]; i++) {
> > +                for (j = 0; j < header->naxisn[0]; j++) {
> > +                    t32 = AV_RB32(ptr8);
> > +                    if (t32 != header->blank) {
> > +                        if (t32 > header->data_max)
> > +                            header->data_max = t32;
> > +                        if (t32 < header->data_min)
> > +                            header->data_min = t32;
> > +                    }
> > +                    ptr8 += 4;
> > +                }
> > +            }
> > +            break;
> > +        case 64:
> > +            t64 = AV_RB64(ptr8);
> > +            for (i = 0; i < header->naxisn[1]; i++) {
> > +                for (j = 0; j < header->naxisn[0]; j++) {
> > +                    t64 = AV_RB64(ptr8);
> > +                    if (t64 != header->blank) {
> > +                        if (t64 > header->data_max)
> > +                            header->data_max = t64;
> > +                        if (t64 < header->data_min)
> > +                            header->data_min = t64;
> > +                    }
> > +                    ptr8 += 8;
> > +                }
> > +            }
> > +            break;
> > +        default:
> > +            return AVERROR_INVALIDDATA;
> > +    }
> > +    return 1;
> > +}
> > +
> > +/**
> > + * function reads the fits header and stores the values in
> FITSDecContext
> > pointed by header
> > + * @param avctx - AVCodec context
> > + * @param ptr - pointer to pointer to the data
> > + * @param header - pointer to the FITSDecContext
> > + * @param end - pointer to end of packet
> > + * @return 1, if calculated successfully, otherwise AVERROR_INVALIDDATA
> > + */
> > +static int fits_read_header(AVCodecContext *avctx, const uint8_t **ptr,
> > FITSDecContext * header, const uint8_t * end)
> > +{
> > +    const uint8_t *ptr8 = *ptr;
> > +    int lines_read = 0, i, dim_no, t, data_min_found = 0,
> data_max_found =
> > 0, ret;
> > +    uint64_t size=1;
> > +    double d;
> > +
> > +    header->blank = 0;
> > +    header->bscale = 1.0;
> > +    header->bzero = 0;
> > +    header->rgb = 0;
> > +
> > +    if (end - ptr8 < 80)
> > +        return AVERROR_INVALIDDATA;
> > +
> > +    if (sscanf(ptr8, "SIMPLE = %c", &header->simple) == 1) {
> > +        if (header->simple == 'F') {
> > +            av_log(avctx, AV_LOG_WARNING, "not a standard FITS file\n");
> > +        } else if (header->simple != 'T') {
> > +            av_log(avctx, AV_LOG_ERROR, "invalid SIMPLE value, SIMPLE =
> > %c\n", header->simple);
> > +            return AVERROR_INVALIDDATA;
> > +        }
> > +        header->xtension = 0;
> > +    } else if (!strncmp(ptr8, "XTENSION= 'IMAGE", 16)) {
> > +        header->xtension = 1;
> > +    } else {
> > +        av_log(avctx, AV_LOG_ERROR, "missing SIMPLE keyword or invalid
> > XTENSION\n");
> > +        return AVERROR_INVALIDDATA;
> > +    }
> > +
> > +    ptr8 += 80;
> > +    lines_read++;
> > +
> > +    if (end - ptr8 < 80)
> > +        return AVERROR_INVALIDDATA;
> > +
> > +    if (sscanf(ptr8, "BITPIX = %d", &header->bitpix) != 1) {
> > +        av_log(avctx, AV_LOG_ERROR, "missing BITPIX keyword\n");
> > +        return AVERROR_INVALIDDATA;
> > +    }
> > +
> > +    size = abs(header->bitpix) >> 3;
> > +    ptr8 += 80;
> > +    lines_read++;
> > +
> > +    if (end - ptr8 < 80)
> > +        return AVERROR_INVALIDDATA;
> > +
> > +    if (sscanf(ptr8, "NAXIS = %d", &header->naxis) != 1) {
> > +        av_log(avctx, AV_LOG_ERROR, "missing NAXIS keyword\n");
> > +        return AVERROR_INVALIDDATA;
> > +    }
> > +
> > +    if (!header->naxis) {
> > +        av_log(avctx, AV_LOG_ERROR, "No image data found, NAXIS = %d\n",
> > header->naxis);
> > +        return AVERROR_INVALIDDATA;
> > +    }
> > +
> > +    if (header->naxis != 2 && header->naxis != 3) {
> > +        av_log(avctx, AV_LOG_ERROR, "unsupported number of dimensions,
> > NAXIS = %d\n", header->naxis);
> > +        return AVERROR_INVALIDDATA;
> > +    }
> > +
> > +    ptr8 += 80;
> > +    lines_read++;
> > +
> > +    for (i = 0; i < header->naxis; i++) {
> > +        if (end - ptr8 < 80)
> > +            return AVERROR_INVALIDDATA;
> > +
> > +        if (sscanf(ptr8, "NAXIS%d = %d", &dim_no, &header->naxisn[i])
> != 2
> > || dim_no != i+1) {
> > +            av_log(avctx, AV_LOG_ERROR, "missing NAXIS%d keyword\n",
> i+1);
> > +            return AVERROR_INVALIDDATA;
> > +        }
> > +
> > +        size *= header->naxisn[i];
> > +        ptr8 += 80;
> > +        lines_read++;
> > +    }
> > +
> > +    if (end - ptr8 < 80)
> > +            return AVERROR_INVALIDDATA;
> > +
> > +    while (strncmp(ptr8, "END", 3)) {
> > +        if (sscanf(ptr8, "BLANK = %d", &t) == 1) {
> > +            header->blank = t;
> > +        } else if (sscanf(ptr8, "BSCALE = %lf", &d) == 1) {
> > +            header->bscale = d;
> > +        } else if (sscanf(ptr8, "BZERO = %lf", &d) == 1) {
> > +            header->bzero = d;
> > +        } else if (sscanf(ptr8, "DATAMAX = %lf", &d) == 1) {
> > +            data_max_found = 1;
> > +            header->data_max = d;
> > +        } else if (sscanf(ptr8, "DATAMIN = %lf", &d) == 1) {
> > +            data_min_found = 1;
> > +            header->data_min = d;
> > +        } else if (!strncmp(ptr8, "CTYPE3  = 'RGB", 14)) {
> > +            header->rgb = 1;
> > +            if (header->naxis != 3 || (header->naxisn[2] != 3 &&
> > header->naxisn[2] != 4)) {
> > +                av_log(avctx, AV_LOG_ERROR, "File contains RGB image but
> > NAXIS = %d and NAXIS3 = %d\n", header->naxis, header->naxisn[2]);
> > +                return AVERROR_INVALIDDATA;
> > +            }
> > +        }
> > +        ptr8 += 80;
> > +        lines_read++;
> > +
> > +        if (end - ptr8 < 80)
> > +            return AVERROR_INVALIDDATA;
> > +    }
> > +
> > +    if (!header->rgb && header->naxis != 2){
> > +        av_log(avctx, AV_LOG_ERROR, "unsupported number of dimensions,
> > NAXIS = %d\n", header->naxis);
> > +        return AVERROR_INVALIDDATA;
> > +    }
> > +
> > +    ptr8 += 80;
> > +    lines_read++;
> > +    lines_read %= 36;
> > +
> > +    t = ((36 - lines_read) % 36) * 80;
> > +    if (end - ptr8 < t)
> > +        return AVERROR_INVALIDDATA;
> > +    ptr8 += t;
> > +    *ptr = ptr8;
> > +
> > +    if (end - ptr8 < size)
> > +        return AVERROR_INVALIDDATA;
> > +
> > +    if (!header->rgb && (!data_min_found || !data_max_found)) {
> > +        if ((ret = fill_data_min_max(ptr8, header, end)) < 0) {
> > +            av_log(avctx, AV_LOG_ERROR, "invalid BITPIX, %d\n",
> > header->bitpix);
> > +            return AVERROR_INVALIDDATA;
> > +        }
> > +    } else {
> > +        /*
> > +         * instead of applying bscale and bzero to every element, we
> can do
> > inverse transformation on data_min and
> > +         * data_max
> > +         */
> > +        header->data_min = (header->data_min - header->bzero) /
> > header->bscale;
> > +        header->data_max = (header->data_max - header->bzero) /
> > header->bscale;
> > +    }
> > +    return 1;
> > +}
> > +
> > +static int fits_decode_frame(AVCodecContext *avctx, void *data, int
> > *got_frame, AVPacket *avpkt)
> > +{
> > +    AVFrame *p=data;
> > +    const uint8_t *ptr8 = avpkt->data, *end;
> > +    int16_t t16;
> > +    int32_t t32;
> > +    int64_t t64;
> > +    float   tflt;
> > +    double  tdbl;
> > +    int ret, i, j;
> > +    uint8_t *dst8;
> > +    uint16_t *dst16;
> > +    uint32_t *dst32;
> > +    uint64_t *dst64, size, r, g, b, a, t;
> > +    FITSDecContext * header = avctx->priv_data;
> > +
> > +    end = ptr8 + avpkt->size;
> > +    if (ret = fits_read_header(avctx, &ptr8, header, end) < 0)
> > +        return ret;
> > +
> > +    size = (header->naxisn[0]) * (header->naxisn[1]);
> > +
> > +    if (header->rgb) {
> > +        if (header->bitpix == 8) {
> > +            avctx->pix_fmt = AV_PIX_FMT_RGB32;
> > +        } else if (header->bitpix == 16) {
> > +            avctx->pix_fmt = AV_PIX_FMT_RGBA64;
> > +        } else {
> > +            av_log(avctx, AV_LOG_ERROR, "unsupported BITPIX = %d\n",
> > header->bitpix);
> > +            return AVERROR_INVALIDDATA;
> > +        }
> > +    } else {
> > +        if (header->bitpix == 8) {
> > +            avctx->pix_fmt = AV_PIX_FMT_GRAY8;
> > +        } else {
> > +            avctx->pix_fmt = AV_PIX_FMT_GRAY16;
> > +        }
> > +    }
> > +
> > +    if ((ret = ff_set_dimensions(avctx, header->naxisn[0],
> > header->naxisn[1])) < 0)
> > +        return ret;
> > +
> > +    if ((ret = ff_get_buffer(avctx, p, 0)) < 0)
> > +        return ret;
> > +
> > +    if (header->rgb) {
> > +        if (header->bitpix == 8) {
> > +            for (i = 0; i < avctx->height; i++) {
> > +                /*
> > +                 * FITS stores images with bottom row first. Therefore
> we
> > have
> > +                 * to fill the image from bottom to top.
> > +                 */
> > +                dst32 = (uint32_t *)(p->data[0] + (avctx->height-i-1)*
> > p->linesize[0]);
> > +                for (j = 0; j < avctx->width; j++) {
> > +                    if (header->naxisn[2] == 4) {
> > +                        if (ptr8[size * 3] != header->blank)
> > +                            t = ptr8[size * 3] * header->bscale +
> > header->bzero;
> > +                        a = t << 24;
> > +                    } else {
> > +                        a = (255 << 24);
> > +                    }
> > +
> > +                    if (ptr8[0] != header->blank)
> > +                        t = ptr8[0] * header->bscale + header->bzero;
> > +                    r = t << 16;
> > +
> > +                    if (ptr8[size] != header->blank)
> > +                        t = ptr8[size] * header->bscale + header->bzero;
> > +                    g = t << 8;
> > +
> > +                    if (ptr8[size * 2] != header->blank)
> > +                        t = ptr8[size * 2] * header->bscale +
> > header->bzero;
> > +                    b = t;
> > +
> > +                    *dst32++ = ((uint32_t)a) | ((uint32_t)r) |
> > ((uint32_t)g) | ((uint32_t)b);
> > +                    ptr8++;
> > +                }
> > +            }
> > +        } else if (header->bitpix == 16) {
> > +            // not tested ....
> > +            for (i = 0; i < avctx->height; i++) {
> > +                dst64 = (uint64_t *)(p->data[0] + (avctx->height-i-1) *
> > p->linesize[0]);
> > +                for (j = 0; j < avctx->width; j++) {
> > +
> > +                    if (header->naxisn[2] == 4) {
> > +                        t = ((ptr8[size * 3] << 8) | ptr8[size * 3 +
> 1]);
> > +                        if (t != header->blank)
> > +                            t = t*header->bscale + header->bzero;
> > +                        a = t << 48;
> > +                    } else {
> > +                        a = 65535ULL << 48;
> > +                    }
> > +
> > +                    t = ptr8[0] << 8 | ptr8[1];
> > +                    if (t != header->blank)
> > +                        t = t*header->bscale + header->bzero;
> > +                    r = t << 32;
> > +
> > +                    t = ptr8[size] << 8 | ptr8[size + 1];
> > +                    if (t != header->blank)
> > +                        t = t*header->bscale + header->bzero;
> > +                    g = t << 16;
> > +
> > +                    t = ptr8[size * 2] << 8 | ptr8[size * 2 + 1];
> > +                    if (t != header->blank)
> > +                        t = t*header->bscale + header->bzero;
> > +                    b = t;
> > +
> > +                    *dst64++ = a | r | g | b;
> > +                    ptr8 += 2;
> > +                }
> > +            }
> > +        }
> > +    } else {
> > +        if (header->bitpix == 8) {
> > +            for (i = 0; i < avctx->height; i++) {
> > +                dst8 = (uint8_t *) (p->data[0] + (avctx->height-i-1)*
> > p->linesize[0]);
> > +                for (j = 0; j < avctx->width; j++) {
> > +                    if (ptr8[0] != header->blank) {
> > +                        *dst8++ = ((ptr8[0] - header->data_min) * 255) /
> > (header->data_max - header->data_min);
> > +                    } else {
> > +                        *dst8++ = ptr8[0];
> > +                    }
> > +                    ptr8++;
> > +                }
> > +            }
> > +        } else if (header->bitpix == 16) {
> > +            for (i = 0; i < avctx->height; i++) {
> > +                dst16 = (uint16_t *)(p->data[0] + (avctx->height-i-1) *
> > p->linesize[0]);
> > +                for (j = 0; j < avctx->width; j++) {
> > +                    t16 = AV_RB16(ptr8);
> > +                    if (t16 != header->blank)
> > +                        t16 = ((t16 - header->data_min) * 65535) /
> > (header->data_max - header->data_min);
> > +                    *dst16++ = t16;
> > +                    ptr8 += 2;
> > +                }
> > +            }
> > +        } else if (header->bitpix == 32) {
> > +            for (i = 0; i < avctx->height; i++) {
> > +                dst16 = (uint16_t *)(p->data[0] + (avctx->height-i-1) *
> > p->linesize[0]);
> > +                for (j = 0; j < avctx->width; j++) {
> > +                    t32 = AV_RB32(ptr8);
> > +                    if (t32 != header->blank)
> > +                        t16 = ((t32 - header->data_min) * 65535) /
> > (header->data_max - header->data_min);
> > +                    *dst16++ = t16;
> > +                    ptr8 += 4;
> > +                }
> > +            }
> > +        } else if (header->bitpix == 64) {
> > +            for (i = 0; i < avctx->height; i++) {
> > +                dst16 = (uint16_t *)(p->data[0] + (avctx->height-i-1) *
> > p->linesize[0]);
> > +                for (j = 0; j < avctx->width; j++) {
> > +                    t64 = AV_RB64(ptr8);
> > +                    if (t64 != header->blank)
> > +                        t16 = ((t64 - header->data_min) * 65535) /
> > (header->data_max - header->data_min);
> > +                    *dst16++ = t16;
> > +                    ptr8 += 8;
> > +                }
> > +            }
> > +        } else if (header->bitpix == -32) {
> > +            for (i = 0; i < avctx->height; i++) {
> > +                dst16 = (uint16_t *)(p->data[0] + (avctx->height-i-1) *
> > p->linesize[0]);
> > +                for (j = 0; j < avctx->width; j++) {
> > +                    t32 = AV_RB32(ptr8);
> > +                    memcpy(&tflt, &t32, 4);
> > +                    *dst16++ = ((tflt - header->data_min) * 65535) /
> > (header->data_max - header->data_min);
> > +                    ptr8 += 4;
> > +                }
> > +            }
> > +        } else if (header->bitpix == -64) {
> > +            for (i = 0; i < avctx->height; i++) {
> > +                dst16 = (uint16_t *)(p->data[0] + (avctx->height-i-1) *
> > p->linesize[0]);
> > +                for (j = 0; j < avctx->width; j++) {
> > +                    t64 = AV_RB64(ptr8);
> > +                    memcpy(&tdbl, &t64, 8);
> > +                    *dst16++ = ((tdbl - header->data_min) * 65535) /
> > (header->data_max - header->data_min);
> > +                    ptr8 += 8;
> > +                }
> > +            }
> > +        } else {
> > +            av_log(avctx, AV_LOG_ERROR, "invalid BITPIX, %d\n",
> > header->bitpix);
> > +            return AVERROR_INVALIDDATA;
> > +        }
> > +    }
> > +
> > +    p->key_frame = 1;
> > +    p->pict_type = AV_PICTURE_TYPE_I;
> > +
> > +    *got_frame = 1;
> > +
> > +    return avpkt->size;
> > +}
> > +
> > +AVCodec ff_fits_decoder = {
> > +    .name           = "fits",
> > +    .type           = AVMEDIA_TYPE_VIDEO,
> > +    .id             = AV_CODEC_ID_FITS,
> > +    .priv_data_size = sizeof(FITSDecContext),
> > +    .decode         = fits_decode_frame,
> > +    .capabilities   = AV_CODEC_CAP_DR1,
> > +    .long_name      = NULL_IF_CONFIG_SMALL("Flexible Image Transport
> > System")
> > +};
> > diff --git a/libavcodec/version.h b/libavcodec/version.h
> > index a44a888..4af01ba 100644
> > --- a/libavcodec/version.h
> > +++ b/libavcodec/version.h
> > @@ -28,7 +28,7 @@
> >  #include "libavutil/version.h"
> >
> >  #define LIBAVCODEC_VERSION_MAJOR  57
> > -#define LIBAVCODEC_VERSION_MINOR  99
> > +#define LIBAVCODEC_VERSION_MINOR  100
> >  #define LIBAVCODEC_VERSION_MICRO 100
> >
> >  #define LIBAVCODEC_VERSION_INT  AV_VERSION_INT(LIBAVCODEC_VERSION_MAJOR,
> \
> > diff --git a/libavformat/img2.c b/libavformat/img2.c
> > index 8432cc0..e405df8 100644
> > --- a/libavformat/img2.c
> > +++ b/libavformat/img2.c
> > @@ -80,6 +80,7 @@ const IdStrMap ff_img_tags[] = {
> >      { AV_CODEC_ID_XPM,        "xpm"      },
> >      { AV_CODEC_ID_XFACE,      "xface"    },
> >      { AV_CODEC_ID_XWD,        "xwd"      },
> > +    { AV_CODEC_ID_FITS,       "fits"     },
> >      { AV_CODEC_ID_NONE,       NULL       }
> >  };
> >
> > --
> > 2.4.11
> >
> > _______________________________________________
> > ffmpeg-devel mailing list
> > ffmpeg-devel at ffmpeg.org
> > http://ffmpeg.org/mailman/listinfo/ffmpeg-devel
> >
>

Rest all changes have been done.


More information about the ffmpeg-devel mailing list