ScaleGrid.cc

 
// -*- mode: c++; c-basic-offset:4 -*-
 
// This file is part of libdap, A C++ implementation of the OPeNDAP Data
// Access Protocol.
 
// Copyright (c) 2016 OPeNDAP, Inc.
// Author: James Gallagher <jgallagher@opendap.org>
//
// This library 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.
//
// This library 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 this library; if not, write to the Free Software
// Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307  USA
//
// You can contact OPeNDAP, Inc. at PO Box 112, Saunderstown, RI. 02874-0112.
 
#include "ScaleGrid.h"
 
#include "config.h"
 
#include <memory>
#include <string>
 
#include <libdap/BaseType.h>
#include <libdap/Array.h>
#include <libdap/Grid.h>
#include <libdap/Str.h>
#include <libdap/Error.h>
#include <libdap/util.h>
 
#include <BESDebug.h>
 
#include "functions_util.h"
#include "GridGeoConstraint.h"
 
#define DEBUG_KEY "geo"
 
using namespace std;
using namespace libdap;
 
namespace functions {

void function_scale_grid(int argc, BaseType *argv[], DDS &, BaseType **btpp)
{
    // scale_grid(argv[0], y, x)
    string info =
    string("<?xml version=\"1.0\" encoding=\"UTF-8\"?>\n") +
    "<function name=\"scale_grid\" version=\"1.0\" href=\"http://docs.opendap.org/index.php/Server_Side_Processing_Functions#scale_grid\">\n" +
    "</function>\n";
 
    if (argc == 0) {
        unique_ptr<Str> response(new Str("info"));
        response->set_value(info);
        *btpp = response.release();
        return;
    }
 
    if (argc < 3 || argc > 5) {
        throw Error("The scale_grid() function requires three arguments: a Grid and the new lon, lat extents (got " + long_to_string(argc) + " args).\n\
             See http://docs.opendap.org/index.php/Server_Side_Processing_Functions#scale_grid");
    }
 
    Grid *src = dynamic_cast < Grid * >(argv[0]);
    if (!src)
        throw Error(malformed_expr,"The first argument to scale_grid() must be a Grid variable!");
 
    GridGeoConstraint gc(src);
    if(!gc.is_longitude_rightmost()){
        throw Error(malformed_expr,"The last argument to scale_grid() must be a longitude variable!");
    }
 
    BESDEBUG(DEBUG_KEY,"function_scale_grid() - Evaluating grid '"<< src->name() << "'" << endl);
    unsigned long y = extract_uint_value(argv[1]);
    unsigned long x = extract_uint_value(argv[2]);
 
    string crs = "WGS84";
    string interp = "nearest";
    if (argc > 3) {
        crs = extract_string_argument(argv[3]);
    }
    if (argc > 4) {
        interp = extract_string_argument(argv[4]);
    }
 
    BESDEBUG(DEBUG_KEY,"function_scale_grid() - CRS '"<< crs << "'" << endl);
    BESDEBUG(DEBUG_KEY,"function_scale_grid() - Interpolation Method '"<< interp << "'" << endl);
 
    SizeBox size(x, y);
    *btpp = scale_dap_grid(src, size, crs, interp);
}
 
//Grid *scale_dap_array(const Array *data, const Array *lon, const Array *lat, const SizeBox &size,
// const string &crs, const string &interp)
void function_scale_array(int argc, BaseType *argv[], DDS &, BaseType **btpp)
{
    // scale_array(argv[0], Y, X, y, x)
    string info =
    string("<?xml version=\"1.0\" encoding=\"UTF-8\"?>\n") +
    "<function name=\"scale_array\" version=\"1.0\" href=\"http://docs.opendap.org/index.php/Server_Side_Processing_Functions#scale_array\">\n" +
    "</function>\n";
 
    if (argc == 0) {
        unique_ptr<Str> response(new Str("info"));
        response->set_value(info);
        *btpp = response.release();
        return;
    }
 
    if (!(argc > 4  && argc < 8)) {
        throw Error("The scale_array() function requires five arguments: three Arrays and the new lat and lon extents.\n\
             See http://docs.opendap.org/index.php/Server_Side_Processing_Functions#scale_array");
    }
 
    Array *data = dynamic_cast <Array *>(argv[0]);
    if (!data)
        throw Error(malformed_expr,"The first argument to scale_array() must be an Array variable!");
    Array *x = dynamic_cast <Array *>(argv[2]);
    if (!x)
        throw Error(malformed_expr,"The second argument to scale_array() must be an Array variable!");
    Array *y = dynamic_cast <Array *>(argv[1]);
    if (!y)
        throw Error(malformed_expr,"The third argument to scale_array() must be an Array variable!");
 
    unsigned long new_x = extract_uint_value(argv[4]);
    unsigned long new_y = extract_uint_value(argv[3]);
 
    string crs = "WGS84";   // FIXME WGS84 assumes a certain order for lat and lon
    string interp = "nearest";
    if (argc > 5) {
        crs = extract_string_argument(argv[5]);
    }
 
    if (argc > 6) {
        interp = extract_string_argument(argv[6]);
    }
 
    SizeBox size(new_x, new_y);
    *btpp = scale_dap_array(data, x, y, size, crs, interp);
}
 
//Grid *scale_dap_3D_array(const Array *data, const Array *time const Array *lon, const Array *lat, const SizeBox &size,
// const string &crs, const string &interp)
void function_scale_array_3D(int argc, BaseType *argv[], DDS &, BaseType **btpp)
{
    // scale_3D_array(arg[0], T, Y, X, y, x)
    string info =
        string("<?xml version=\"1.0\" encoding=\"UTF-8\"?>\n")
            + "<function name=\"scale_array\" version=\"1.0\" href=\"http://docs.opendap.org/index.php/Server_Side_Processing_Functions#scale_3D_array\">\n"
            + "</function>\n";
 
    BESDEBUG(DEBUG_KEY,"function_scale_array_3D() - argc = " << argc << endl);
 
    if (argc == 0) {
        unique_ptr<Str> response(new Str("info"));
        response->set_value(info);
        *btpp = response.release();
        return;
    }
 
    if (!(argc > 5 && argc < 9)) {
        throw Error(
            "The scale_array_3D() function requires six arguments: four Arrays and the new lat and lon extents.\n\
             See http://docs.opendap.org/index.php/Server_Side_Processing_Functions#scale_array");
    }
 
    BESDEBUG(DEBUG_KEY,"function_scale_array_3D() - arg[0] name: " << argv[0]->name() << endl);
 
    Array *data = dynamic_cast<Array *>(argv[0]);
    if (!data) throw Error(malformed_expr, "The first argument to scale_array_3D() must be an Array variable!");
    Array *t = dynamic_cast<Array *>(argv[1]);
    if (!t) throw Error(malformed_expr, "The second argument to scale_array_3D() must be an Array variable!");
    Array *y = dynamic_cast<Array *>(argv[2]);
    if (!y) throw Error(malformed_expr, "The third argument to scale_array_3D() must be an Array variable!");
    Array *x = dynamic_cast<Array *>(argv[3]);
    if (!x) throw Error(malformed_expr, "The fourth argument to scale_array_3D() must be an Array variable!");
 
    unsigned long new_x = extract_uint_value(argv[5]); // lon should be last argument
    unsigned long new_y = extract_uint_value(argv[4]);
 
    string crs = "WGS84";   // FIXME WGS84 assumes a certain order for lat and lon
    string interp = "nearest";
    if (argc > 6) {
        crs = extract_string_argument(argv[6]);
    }
 
    if (argc > 7) {
        interp = extract_string_argument(argv[7]);
    }
 
    SizeBox size(new_x, new_y);
    *btpp = scale_dap_array_3D(data, t, x, y, size, crs, interp);
}
 
 
}