Filling NA gaps in raster with R? The 2019 Stack Overflow Developer Survey Results Are In ...

Am I ethically obligated to go into work on an off day if the reason is sudden?

Word for: a synonym with a positive connotation?

Why can't wing-mounted spoilers be used to steepen approaches?

Using dividends to reduce short term capital gains?

Student Loan from years ago pops up and is taking my salary

How to read αἱμύλιος or when to aspirate

One-dimensional Japanese puzzle

Why don't hard Brexiteers insist on a hard border to prevent illegal immigration after Brexit?

Do I have Disadvantage attacking with an off-hand weapon?

Windows 10: How to Lock (not sleep) laptop on lid close?

Single author papers against my advisor's will?

How to politely respond to generic emails requesting a PhD/job in my lab? Without wasting too much time

Sub-subscripts in strings cause different spacings than subscripts

Was credit for the black hole image misappropriated?

Why are PDP-7-style microprogrammed instructions out of vogue?

US Healthcare consultation for visitors

For what reasons would an animal species NOT cross a *horizontal* land bridge?

The following signatures were invalid: EXPKEYSIG 1397BC53640DB551

Why doesn't a hydraulic lever violate conservation of energy?

Are spiders unable to hurt humans, especially very small spiders?

What force causes entropy to increase?

Is it ok to offer lower paid work as a trial period before negotiating for a full-time job?

Is it ethical to upload a automatically generated paper to a non peer-reviewed site as part of a larger research?

Sort list of array linked objects by keys and values



Filling NA gaps in raster with R?



The 2019 Stack Overflow Developer Survey Results Are In
Announcing the arrival of Valued Associate #679: Cesar Manara
Planned maintenance scheduled April 17/18, 2019 at 00:00UTC (8:00pm US/Eastern)How to grid unevenly sampled categorical data?Fast and robust open source portable script or tool for interpolation - triangulation or similarR package gapfill: How to convert R raster stack to 4-dimensional array and then back to raster stack?Filling in Null Data Gaps in Raster ImageRaster reclassification in QGIS or other open source software with automatically calculated breaksRainfall raster with spline in RFilling/replacing nodata values of a raster layer in RFilling no data value in raster using ArcGIS Desktop?Avoiding artefacts by pixel filtering in a rasterImporting several LiDAR (.las) files in GRASS using v.in.lidar?





.everyoneloves__top-leaderboard:empty,.everyoneloves__mid-leaderboard:empty,.everyoneloves__bot-mid-leaderboard:empty{ margin-bottom:0;
}







1















I have soil data rasters from SoilGrids.org that have some small holes with NA cells (White spaces in the image) that I want to fill using some interpolation method to get a more spatially continuous raster for further analysis. I have found some options in SAGA (e.g. Close gaps) and GRASS (e.g. r.fillnulls) that seems to do the job but I am having problems to implement them in my workflow in R. I have found meteo::rfillspgaps , but it takes too much time and seems to never ends.



Is there any other function or solution in R?



enter image description here










share|improve this question

























  • You can use GRASS and SAGA trough R too

    – aldo_tapia
    Jul 18 '18 at 19:51


















1















I have soil data rasters from SoilGrids.org that have some small holes with NA cells (White spaces in the image) that I want to fill using some interpolation method to get a more spatially continuous raster for further analysis. I have found some options in SAGA (e.g. Close gaps) and GRASS (e.g. r.fillnulls) that seems to do the job but I am having problems to implement them in my workflow in R. I have found meteo::rfillspgaps , but it takes too much time and seems to never ends.



Is there any other function or solution in R?



enter image description here










share|improve this question

























  • You can use GRASS and SAGA trough R too

    – aldo_tapia
    Jul 18 '18 at 19:51














1












1








1








I have soil data rasters from SoilGrids.org that have some small holes with NA cells (White spaces in the image) that I want to fill using some interpolation method to get a more spatially continuous raster for further analysis. I have found some options in SAGA (e.g. Close gaps) and GRASS (e.g. r.fillnulls) that seems to do the job but I am having problems to implement them in my workflow in R. I have found meteo::rfillspgaps , but it takes too much time and seems to never ends.



Is there any other function or solution in R?



enter image description here










share|improve this question
















I have soil data rasters from SoilGrids.org that have some small holes with NA cells (White spaces in the image) that I want to fill using some interpolation method to get a more spatially continuous raster for further analysis. I have found some options in SAGA (e.g. Close gaps) and GRASS (e.g. r.fillnulls) that seems to do the job but I am having problems to implement them in my workflow in R. I have found meteo::rfillspgaps , but it takes too much time and seems to never ends.



Is there any other function or solution in R?



enter image description here







raster r interpolation nodata holes






share|improve this question















share|improve this question













share|improve this question




share|improve this question








edited Apr 8 '18 at 16:23







Guido A. Herrera

















asked Apr 6 '18 at 10:27









Guido A. HerreraGuido A. Herrera

17012




17012













  • You can use GRASS and SAGA trough R too

    – aldo_tapia
    Jul 18 '18 at 19:51



















  • You can use GRASS and SAGA trough R too

    – aldo_tapia
    Jul 18 '18 at 19:51

















You can use GRASS and SAGA trough R too

– aldo_tapia
Jul 18 '18 at 19:51





You can use GRASS and SAGA trough R too

– aldo_tapia
Jul 18 '18 at 19:51










2 Answers
2






active

oldest

votes


















1














I've only found two ways. The ideal but computationally heavy way is to convert the raster to SpatialPixels and then use idw() or krige() in gstats package for interpolation, and convert back to raster.



The quick and dirty way is to use focal in the raster package with fun=mean, NAonly=T, na.rm=T and an appropriately sized matrix of 1's as the weights.






share|improve this answer































    0














    unfortunately the function in meteo can not run in parallel, I am having the same issue and now trying to solve it making tiles






    share|improve this answer








    New contributor




    user140425 is a new contributor to this site. Take care in asking for clarification, commenting, and answering.
    Check out our Code of Conduct.





















      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',
      autoActivateHeartbeat: false,
      convertImagesToLinks: false,
      noModals: true,
      showLowRepImageUploadWarning: true,
      reputationToPostImages: null,
      bindNavPrevention: true,
      postfix: "",
      imageUploader: {
      brandingHtml: "Powered by u003ca class="icon-imgur-white" href="https://imgur.com/"u003eu003c/au003e",
      contentPolicyHtml: "User contributions licensed under u003ca href="https://creativecommons.org/licenses/by-sa/3.0/"u003ecc by-sa 3.0 with attribution requiredu003c/au003e u003ca href="https://stackoverflow.com/legal/content-policy"u003e(content policy)u003c/au003e",
      allowUrls: true
      },
      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%2f278476%2ffilling-na-gaps-in-raster-with-r%23new-answer', 'question_page');
      }
      );

      Post as a guest















      Required, but never shown

























      2 Answers
      2






      active

      oldest

      votes








      2 Answers
      2






      active

      oldest

      votes









      active

      oldest

      votes






      active

      oldest

      votes









      1














      I've only found two ways. The ideal but computationally heavy way is to convert the raster to SpatialPixels and then use idw() or krige() in gstats package for interpolation, and convert back to raster.



      The quick and dirty way is to use focal in the raster package with fun=mean, NAonly=T, na.rm=T and an appropriately sized matrix of 1's as the weights.






      share|improve this answer




























        1














        I've only found two ways. The ideal but computationally heavy way is to convert the raster to SpatialPixels and then use idw() or krige() in gstats package for interpolation, and convert back to raster.



        The quick and dirty way is to use focal in the raster package with fun=mean, NAonly=T, na.rm=T and an appropriately sized matrix of 1's as the weights.






        share|improve this answer


























          1












          1








          1







          I've only found two ways. The ideal but computationally heavy way is to convert the raster to SpatialPixels and then use idw() or krige() in gstats package for interpolation, and convert back to raster.



          The quick and dirty way is to use focal in the raster package with fun=mean, NAonly=T, na.rm=T and an appropriately sized matrix of 1's as the weights.






          share|improve this answer













          I've only found two ways. The ideal but computationally heavy way is to convert the raster to SpatialPixels and then use idw() or krige() in gstats package for interpolation, and convert back to raster.



          The quick and dirty way is to use focal in the raster package with fun=mean, NAonly=T, na.rm=T and an appropriately sized matrix of 1's as the weights.







          share|improve this answer












          share|improve this answer



          share|improve this answer










          answered Jul 18 '18 at 17:12









          user5730487user5730487

          112




          112

























              0














              unfortunately the function in meteo can not run in parallel, I am having the same issue and now trying to solve it making tiles






              share|improve this answer








              New contributor




              user140425 is a new contributor to this site. Take care in asking for clarification, commenting, and answering.
              Check out our Code of Conduct.

























                0














                unfortunately the function in meteo can not run in parallel, I am having the same issue and now trying to solve it making tiles






                share|improve this answer








                New contributor




                user140425 is a new contributor to this site. Take care in asking for clarification, commenting, and answering.
                Check out our Code of Conduct.























                  0












                  0








                  0







                  unfortunately the function in meteo can not run in parallel, I am having the same issue and now trying to solve it making tiles






                  share|improve this answer








                  New contributor




                  user140425 is a new contributor to this site. Take care in asking for clarification, commenting, and answering.
                  Check out our Code of Conduct.










                  unfortunately the function in meteo can not run in parallel, I am having the same issue and now trying to solve it making tiles







                  share|improve this answer








                  New contributor




                  user140425 is a new contributor to this site. Take care in asking for clarification, commenting, and answering.
                  Check out our Code of Conduct.









                  share|improve this answer



                  share|improve this answer






                  New contributor




                  user140425 is a new contributor to this site. Take care in asking for clarification, commenting, and answering.
                  Check out our Code of Conduct.









                  answered 16 mins ago









                  user140425user140425

                  1




                  1




                  New contributor




                  user140425 is a new contributor to this site. Take care in asking for clarification, commenting, and answering.
                  Check out our Code of Conduct.





                  New contributor





                  user140425 is a new contributor to this site. Take care in asking for clarification, commenting, and answering.
                  Check out our Code of Conduct.






                  user140425 is a new contributor to this site. Take care in asking for clarification, commenting, and answering.
                  Check out our Code of Conduct.






























                      draft saved

                      draft discarded




















































                      Thanks for contributing an answer to Geographic Information Systems Stack Exchange!


                      • Please be sure to answer the question. Provide details and share your research!

                      But avoid



                      • Asking for help, clarification, or responding to other answers.

                      • Making statements based on opinion; back them up with references or personal experience.


                      To learn more, see our tips on writing great answers.




                      draft saved


                      draft discarded














                      StackExchange.ready(
                      function () {
                      StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fgis.stackexchange.com%2fquestions%2f278476%2ffilling-na-gaps-in-raster-with-r%23new-answer', 'question_page');
                      }
                      );

                      Post as a guest















                      Required, but never shown





















































                      Required, but never shown














                      Required, but never shown












                      Required, but never shown







                      Required, but never shown

































                      Required, but never shown














                      Required, but never shown












                      Required, but never shown







                      Required, but never shown







                      Popular posts from this blog

                      Щит и меч (фильм) Содержание Названия серий | Сюжет |...

                      is 'sed' thread safeWhat should someone know about using Python scripts in the shell?Nexenta bash script uses...

                      Meter-Bus Содержание Параметры шины | Стандартизация |...