Can't get ST_MapAlgebra to return fractional values

The name of the pictureThe name of the pictureThe name of the pictureClash Royale CLAN TAG#URR8PPP











up vote
3
down vote

favorite












I have a 4-band raster dataset that I have imported into a Postgres/PostGIS table using raster2pgsql. I am Using ST_MapAlgebra() with a callback function to perform some band math on the dataset. I have set the return type as double precision. However, I am finding that anytime I return a decimal value (e.g., .25 or -.25) it only returns zero. However, if I return a whole number (e.g., 25 or -15), I can receive the values back from the call. I played around with this, but can’t figure out why this is occurring. Any help would be appreciated. Here is the syntax of my callback function:



CREATE OR REPLACE FUNCTION urban_veg.ndvi_callback (
 value double precision, 
 pos integer, 
 VARIADIC userargs text
) 
RETURNS double precision AS 
$$
 DECLARE
 r1_val double precision;
 r4_val double precision;
 result double precision;
 BEGIN
 r1_val := value[1][1][1]::double precision;
 r4_val := value[2][1][1]::double precision;
 result := (r4_val-r1_val)/(r4_val+r1_val);
 RETURN result;
 END;
$$ 
LANGUAGE 'plpgsql' IMMUTABLE COST 1000;


The call to the function is:



ST_MapAlgebra(
 ARRAY[
 ROW(rast, 1), 
 ROW(rast, 4) ]::rastbandarg, 
 'urban_veg.ndvi_callback(
 double precision,integer,text
 )'::regprocedure, '16BSI'::text)





postgis raster postgresql







share|improve this question



























    up vote
    3
    down vote

    favorite












    I have a 4-band raster dataset that I have imported into a Postgres/PostGIS table using raster2pgsql. I am Using ST_MapAlgebra() with a callback function to perform some band math on the dataset. I have set the return type as double precision. However, I am finding that anytime I return a decimal value (e.g., .25 or -.25) it only returns zero. However, if I return a whole number (e.g., 25 or -15), I can receive the values back from the call. I played around with this, but can’t figure out why this is occurring. Any help would be appreciated. Here is the syntax of my callback function:



    CREATE OR REPLACE FUNCTION urban_veg.ndvi_callback (
     value double precision, 
     pos integer, 
     VARIADIC userargs text
    ) 
    RETURNS double precision AS 
    $$
     DECLARE
     r1_val double precision;
     r4_val double precision;
     result double precision;
     BEGIN
     r1_val := value[1][1][1]::double precision;
     r4_val := value[2][1][1]::double precision;
     result := (r4_val-r1_val)/(r4_val+r1_val);
     RETURN result;
     END;
    $$ 
    LANGUAGE 'plpgsql' IMMUTABLE COST 1000;


    The call to the function is:



    ST_MapAlgebra(
     ARRAY[
     ROW(rast, 1), 
     ROW(rast, 4) ]::rastbandarg, 
     'urban_veg.ndvi_callback(
     double precision,integer,text
     )'::regprocedure, '16BSI'::text)





    postgis raster postgresql







    share|improve this question

























      up vote
      3
      down vote

      favorite









      up vote
      3
      down vote

      favorite











      I have a 4-band raster dataset that I have imported into a Postgres/PostGIS table using raster2pgsql. I am Using ST_MapAlgebra() with a callback function to perform some band math on the dataset. I have set the return type as double precision. However, I am finding that anytime I return a decimal value (e.g., .25 or -.25) it only returns zero. However, if I return a whole number (e.g., 25 or -15), I can receive the values back from the call. I played around with this, but can’t figure out why this is occurring. Any help would be appreciated. Here is the syntax of my callback function:



      CREATE OR REPLACE FUNCTION urban_veg.ndvi_callback (
       value double precision, 
       pos integer, 
       VARIADIC userargs text
      ) 
      RETURNS double precision AS 
      $$
       DECLARE
       r1_val double precision;
       r4_val double precision;
       result double precision;
       BEGIN
       r1_val := value[1][1][1]::double precision;
       r4_val := value[2][1][1]::double precision;
       result := (r4_val-r1_val)/(r4_val+r1_val);
       RETURN result;
       END;
      $$ 
      LANGUAGE 'plpgsql' IMMUTABLE COST 1000;


      The call to the function is:



      ST_MapAlgebra(
       ARRAY[
       ROW(rast, 1), 
       ROW(rast, 4) ]::rastbandarg, 
       'urban_veg.ndvi_callback(
       double precision,integer,text
       )'::regprocedure, '16BSI'::text)





      postgis raster postgresql







      share|improve this question















      I have a 4-band raster dataset that I have imported into a Postgres/PostGIS table using raster2pgsql. I am Using ST_MapAlgebra() with a callback function to perform some band math on the dataset. I have set the return type as double precision. However, I am finding that anytime I return a decimal value (e.g., .25 or -.25) it only returns zero. However, if I return a whole number (e.g., 25 or -15), I can receive the values back from the call. I played around with this, but can’t figure out why this is occurring. Any help would be appreciated. Here is the syntax of my callback function:



      CREATE OR REPLACE FUNCTION urban_veg.ndvi_callback (
       value double precision, 
       pos integer, 
       VARIADIC userargs text
      ) 
      RETURNS double precision AS 
      $$
       DECLARE
       r1_val double precision;
       r4_val double precision;
       result double precision;
       BEGIN
       r1_val := value[1][1][1]::double precision;
       r4_val := value[2][1][1]::double precision;
       result := (r4_val-r1_val)/(r4_val+r1_val);
       RETURN result;
       END;
      $$ 
      LANGUAGE 'plpgsql' IMMUTABLE COST 1000;


      The call to the function is:



      ST_MapAlgebra(
       ARRAY[
       ROW(rast, 1), 
       ROW(rast, 4) ]::rastbandarg, 
       'urban_veg.ndvi_callback(
       double precision,integer,text
       )'::regprocedure, '16BSI'::text)





      postgis raster postgresql




      postgis raster postgresql






      share|improve this question















      share|improve this question













      share|improve this question




      share|improve this question








      edited 8 hours ago









      ahmadhanb

      19.4k21746




      19.4k21746










      asked 8 hours ago









      geo derek

      180115




      180115




















          1 Answer
          1






          active

          oldest

          votes

















          up vote
          4
          down vote













          You specify the return pixeltype to be '16BSI' (16bit signed integer); given this, ST_MapAlgebra won't return a float, not even with a gun to it's head...



          Try with '32BF' (32bit float).






          share|improve this answer




















            Your Answer







            StackExchange.ready(function()
            var channelOptions =
            tags: "".split(" "),
            id: "79"
            ;
            initTagRenderer("".split(" "), "".split(" "), channelOptions);

            StackExchange.using("externalEditor", function()
            // Have to fire editor after snippets, if snippets enabled
            if (StackExchange.settings.snippets.snippetsEnabled)
            StackExchange.using("snippets", function()
            createEditor();
            );

            else
            createEditor();

            );

            function createEditor()
            StackExchange.prepareEditor(
            heartbeatType: 'answer',
            convertImagesToLinks: false,
            noModals: false,
            showLowRepImageUploadWarning: true,
            reputationToPostImages: null,
            bindNavPrevention: true,
            postfix: "",
            onDemand: true,
            discardSelector: ".discard-answer"
            ,immediatelyShowMarkdownHelp:true
            );



            );













             

            draft saved


            draft discarded


















            StackExchange.ready(
            function ()
            StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fgis.stackexchange.com%2fquestions%2f299593%2fcant-get-st-mapalgebra-to-return-fractional-values%23new-answer', 'question_page');

            );

            Post as a guest

















            1 Answer
            1






            active

            oldest

            votes








            1 Answer
            1






            active

            oldest

            votes









            active

            oldest

            votes






            active

            oldest

            votes








            up vote
            4
            down vote













            You specify the return pixeltype to be '16BSI' (16bit signed integer); given this, ST_MapAlgebra won't return a float, not even with a gun to it's head...



            Try with '32BF' (32bit float).






            share|improve this answer
























              up vote
              4
              down vote













              You specify the return pixeltype to be '16BSI' (16bit signed integer); given this, ST_MapAlgebra won't return a float, not even with a gun to it's head...



              Try with '32BF' (32bit float).






              share|improve this answer






















                up vote
                4
                down vote










                up vote
                4
                down vote









                You specify the return pixeltype to be '16BSI' (16bit signed integer); given this, ST_MapAlgebra won't return a float, not even with a gun to it's head...



                Try with '32BF' (32bit float).






                share|improve this answer












                You specify the return pixeltype to be '16BSI' (16bit signed integer); given this, ST_MapAlgebra won't return a float, not even with a gun to it's head...



                Try with '32BF' (32bit float).







                share|improve this answer












                share|improve this answer



                share|improve this answer










                answered 1 hour ago









                ThingumaBob

                4,8431222




                4,8431222



























                     

                    draft saved


                    draft discarded















































                     


                    draft saved


                    draft discarded














                    StackExchange.ready(
                    function ()
                    StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fgis.stackexchange.com%2fquestions%2f299593%2fcant-get-st-mapalgebra-to-return-fractional-values%23new-answer', 'question_page');

                    );

                    Post as a guest
































































                    Comments

                    Post a Comment

                    Popular posts from this blog

                    Long meetings (6-7 hours a day): Being “babysat” by supervisor

                    Image
                    Clash Royale CLAN TAG #URR8PPP .everyoneloves__top-leaderboard:empty,.everyoneloves__mid-leaderboard:empty margin-bottom:0; up vote 31 down vote favorite 4 I hold a PhD in mechanical engineering with a 7 year background in self-guided, self-paced research including 2 years as a graduate level course instructor and a subject matter expert for various engineering projects in a university setting in the United States. For the last 10-11 months, I am working in a research industry in France where my task is to debug C++ spaghetti code written over a period of several years by my current supervisor. It is a team of 2: just I and my supervisor. Diurnally, I am posed with a "bug of the day" to get out of this C++ code. My supervisor generally gives me about 15-30 minutes to solve it and then accosts me about whether or not the task is completed and then rinses and repeats this until the end of the day. Off-late, my supervisor has been holding 6-7 hour long "meetin
                    Read more

                    What does second last employer means? [closed]

                    Image
                    Clash Royale CLAN TAG #URR8PPP .everyoneloves__top-leaderboard:empty,.everyoneloves__mid-leaderboard:empty margin-bottom:0; up vote -3 down vote favorite Say, for example I have worked in the following companies, Company A => 2010 - 2011 Company B => 2011 - 2012 Company C => 2012 - Till Date Company D => Have offer Now assume Company D is asking to submit documents from my second last employer. So from the list of companies given, which company would be qualified as second last employer? software-industry employer share | improve this question asked Aug 10 '14 at 5:16 Rajaprabhu Aravindasamy 559 1 6 13 closed as off-topic by Jim G., jmort253 ♦ Aug 11 '14 at 1:58 This question appears to be off-topic. The users who voted to close gave this specific reason: "Questions seeking advice on company-specific regulations, agreements, or policies should be directed to your manager or HR department. Questions that address only
                    Read more

                    One-line joke

                    A one-liner is a joke that is delivered in a single line. A good one-liner is said to be pithy - concise and meaningful. [1] Comedians and actors use this comedic method as part of their act, e.g. Jimmy Carr, Tommy Cooper, Rodney Dangerfield, Norm Macdonald, Ken Dodd, Stewart Francis, Zach Galifianakis, Mitch Hedberg, Anthony Jeselnik, Milton Jones, Shappi Khorsandi, Jay London, Mark Linn-Baker, Demetri Martin, Groucho Marx, Dan Mintz, Emo Philips, Tim Vine, Steven Wright and Henny Youngman and Arnold Schwarzenegger. [ citation needed ] Many fictional characters are also known to deliver one-liners, including James Bond, who usually includes pithy and laconic quips after disposing of a villain. [ citation needed ] Examples "A baby seal walks into a club." "A dyslexic man walks into a bra." (Anonymous) "There are three types of people, those who can count and those who can't." "Venison's dear, isn't it?" (Jimmy Carr) "An escala
                    Read more