[FFmpeg-cvslog] Add FITS Encoder
Paras Chadha
git at videolan.org
Wed Aug 30 13:20:15 EEST 2017
ffmpeg | branch: master | Paras Chadha <paraschadha18 at gmail.com> | Tue Aug 29 23:07:30 2017 +0530| [6e02f66f1bc7c9d559bc7ef3393c66cd800d4651] | committer: Paul B Mahol
Add FITS Encoder
Signed-off-by: Paras Chadha <paraschadha18 at gmail.com>
> http://git.videolan.org/gitweb.cgi/ffmpeg.git/?a=commit;h=6e02f66f1bc7c9d559bc7ef3393c66cd800d4651
---
Changelog | 1 +
doc/general.texi | 2 +-
libavcodec/Makefile | 1 +
libavcodec/allcodecs.c | 2 +-
libavcodec/fitsenc.c | 129 +++++++++++++++++++++++++++++++++++++++++++++++++
5 files changed, 133 insertions(+), 2 deletions(-)
diff --git a/Changelog b/Changelog
index 548f67d202..830941757e 100644
--- a/Changelog
+++ b/Changelog
@@ -39,6 +39,7 @@ version <next>:
blend, libvmaf, lut3d, overlay, psnr, ssim.
They must always be used by name.
- FITS demuxer and decoder
+- FITS muxer and encoder
version 3.3:
- CrystalHD decoder moved to new decode API
diff --git a/doc/general.texi b/doc/general.texi
index 686460e5aa..5299a9b58c 100644
--- a/doc/general.texi
+++ b/doc/general.texi
@@ -593,7 +593,7 @@ following image formats are supported:
@tab Digital Picture Exchange
@item EXR @tab @tab X
@tab OpenEXR
- at item FITS @tab @tab X
+ at item FITS @tab X @tab X
@tab Flexible Image Transport System
@item JPEG @tab X @tab X
@tab Progressive JPEG is not supported.
diff --git a/libavcodec/Makefile b/libavcodec/Makefile
index 6988c03596..999632cf9e 100644
--- a/libavcodec/Makefile
+++ b/libavcodec/Makefile
@@ -292,6 +292,7 @@ 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 fits.o
+OBJS-$(CONFIG_FITS_ENCODER) += fitsenc.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 dbde1f985a..ce0bc7ecf3 100644
--- a/libavcodec/allcodecs.c
+++ b/libavcodec/allcodecs.c
@@ -192,7 +192,7 @@ static void register_all(void)
REGISTER_ENCDEC (FFV1, ffv1);
REGISTER_ENCDEC (FFVHUFF, ffvhuff);
REGISTER_DECODER(FIC, fic);
- REGISTER_DECODER(FITS, fits);
+ REGISTER_ENCDEC (FITS, fits);
REGISTER_ENCDEC (FLASHSV, flashsv);
REGISTER_ENCDEC (FLASHSV2, flashsv2);
REGISTER_DECODER(FLIC, flic);
diff --git a/libavcodec/fitsenc.c b/libavcodec/fitsenc.c
new file mode 100644
index 0000000000..b44507e436
--- /dev/null
+++ b/libavcodec/fitsenc.c
@@ -0,0 +1,129 @@
+/*
+ * FITS image encoder
+ * 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 encoder
+ *
+ * Specification: https://fits.gsfc.nasa.gov/fits_standard.html Version 3.0
+ *
+ * RGBA images are encoded as planes in RGBA order. So, NAXIS3 is 3 or 4 for them.
+ * Also CTYPE3 = 'RGB ' is added to the header to distinguish them from 3d images.
+ */
+
+#include "libavutil/intreadwrite.h"
+#include "avcodec.h"
+#include "bytestream.h"
+#include "internal.h"
+
+static int fits_encode_frame(AVCodecContext *avctx, AVPacket *pkt,
+ const AVFrame *pict, int *got_packet)
+{
+ AVFrame * const p = (AVFrame *)pict;
+ uint8_t *bytestream, *bytestream_start, *ptr;
+ const uint16_t flip = (1 << 15);
+ uint64_t data_size = 0, padded_data_size = 0;
+ int ret, bitpix, naxis3 = 1, i, j, k, bytes_left;
+ int map[] = {2, 0, 1, 3}; // mapping from GBRA -> RGBA as RGBA is to be stored in FITS file..
+
+ switch (avctx->pix_fmt) {
+ case AV_PIX_FMT_GRAY8:
+ case AV_PIX_FMT_GRAY16BE:
+ map[0] = 0; // grayscale images should be directly mapped
+ if (avctx->pix_fmt == AV_PIX_FMT_GRAY8) {
+ bitpix = 8;
+ } else {
+ bitpix = 16;
+ }
+ break;
+ case AV_PIX_FMT_GBRP:
+ case AV_PIX_FMT_GBRAP:
+ bitpix = 8;
+ if (avctx->pix_fmt == AV_PIX_FMT_GBRP) {
+ naxis3 = 3;
+ } else {
+ naxis3 = 4;
+ }
+ break;
+ case AV_PIX_FMT_GBRP16BE:
+ case AV_PIX_FMT_GBRAP16BE:
+ bitpix = 16;
+ if (avctx->pix_fmt == AV_PIX_FMT_GBRP16BE) {
+ naxis3 = 3;
+ } else {
+ naxis3 = 4;
+ }
+ break;
+ default:
+ av_log(avctx, AV_LOG_ERROR, "unsupported pixel format\n");
+ return AVERROR(EINVAL);
+ }
+
+ data_size = (bitpix >> 3) * avctx->height * avctx->width * naxis3;
+ padded_data_size = ((data_size + 2879) / 2880 ) * 2880;
+
+ if ((ret = ff_alloc_packet2(avctx, pkt, padded_data_size, 0)) < 0)
+ return ret;
+
+ bytestream_start =
+ bytestream = pkt->data;
+
+ for (k = 0; k < naxis3; k++) {
+ for (i = 0; i < avctx->height; i++) {
+ ptr = p->data[map[k]] + (avctx->height - i - 1) * p->linesize[map[k]];
+ if (bitpix == 16) {
+ for (j = 0; j < avctx->width; j++) {
+ // subtracting bzero is equivalent to first bit flip
+ bytestream_put_be16(&bytestream, AV_RB16(ptr) ^ flip);
+ ptr += 2;
+ }
+ } else {
+ memcpy(bytestream, ptr, avctx->width);
+ bytestream += avctx->width;
+ }
+ }
+ }
+
+ bytes_left = padded_data_size - data_size;
+ memset(bytestream, 0, bytes_left);
+ bytestream += bytes_left;
+
+ pkt->size = bytestream - bytestream_start;
+ pkt->flags |= AV_PKT_FLAG_KEY;
+ *got_packet = 1;
+
+ return 0;
+}
+
+AVCodec ff_fits_encoder = {
+ .name = "fits",
+ .long_name = NULL_IF_CONFIG_SMALL("Flexible Image Transport System"),
+ .type = AVMEDIA_TYPE_VIDEO,
+ .id = AV_CODEC_ID_FITS,
+ .encode2 = fits_encode_frame,
+ .pix_fmts = (const enum AVPixelFormat[]) { AV_PIX_FMT_GBRAP16BE,
+ AV_PIX_FMT_GBRP16BE,
+ AV_PIX_FMT_GBRP,
+ AV_PIX_FMT_GBRAP,
+ AV_PIX_FMT_GRAY16BE,
+ AV_PIX_FMT_GRAY8,
+ AV_PIX_FMT_NONE },
+};
More information about the ffmpeg-cvslog
mailing list