Filling NA gaps in raster with R? Planned maintenance scheduled April 23, 2019 at 23:30 UTC (7:30pm US/Eastern) Announcing the arrival of Valued Associate #679: Cesar Manara Unicorn Meta Zoo #1: Why another podcast?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?

How to name indistinguishable henchmen in a screenplay?

Why can't fire hurt Daenerys but it did to Jon Snow in season 1?

Why do C and C++ allow the expression (int) + 4*5?

Flight departed from the gate 5 min before scheduled departure time. Refund options

Weaponising the Grasp-at-a-Distance spell

The test team as an enemy of development? And how can this be avoided?

Found this skink in my tomato plant bucket. Is he trapped? Or could he leave if he wanted?

Plotting a Maclaurin series

NIntegrate on a solution of a matrix ODE

Pointing to problems without suggesting solutions

Searching extreme points of polyhedron

Vertical ranges of Column Plots in 12

Can gravitational waves pass through a black hole?

Should man-made satellites feature an intelligent inverted "cow catcher"?

French equivalents of おしゃれは足元から (Every good outfit starts with the shoes)

What does 丫 mean? 丫是什么意思?

Is there a spell that can create a permanent fire?

Why do the Z-fighters hide their power?

3D Masyu - A Die

How do you cope with tons of web fonts when copying and pasting from web pages?

Did pre-Columbian Americans know the spherical shape of the Earth?

Why are two-digit numbers in Jonathan Swift's "Gulliver's Travels" (1726) written in "German style"?

How do you write "wild blueberries flavored"?

As a dual citizen, my US passport will expire one day after traveling to the US. Will this work?



Filling NA gaps in raster with R?



Planned maintenance scheduled April 23, 2019 at 23:30 UTC (7:30pm US/Eastern)
Announcing the arrival of Valued Associate #679: Cesar Manara
Unicorn Meta Zoo #1: Why another podcast?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;








2















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











  • Have you seen this?: stackoverflow.com/a/14907988/1446289

    – Aaron
    Apr 13 at 3:53

















2















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











  • Have you seen this?: stackoverflow.com/a/14907988/1446289

    – Aaron
    Apr 13 at 3:53













2












2








2








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

17512




17512












  • You can use GRASS and SAGA trough R too

    – aldo_tapia
    Jul 18 '18 at 19:51











  • Have you seen this?: stackoverflow.com/a/14907988/1446289

    – Aaron
    Apr 13 at 3:53

















  • You can use GRASS and SAGA trough R too

    – aldo_tapia
    Jul 18 '18 at 19:51











  • Have you seen this?: stackoverflow.com/a/14907988/1446289

    – Aaron
    Apr 13 at 3:53
















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













Have you seen this?: stackoverflow.com/a/14907988/1446289

– Aaron
Apr 13 at 3:53





Have you seen this?: stackoverflow.com/a/14907988/1446289

– Aaron
Apr 13 at 3:53










1 Answer
1






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























    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

























    1 Answer
    1






    active

    oldest

    votes








    1 Answer
    1






    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



























            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

            រឿង រ៉ូមេអូ និង ហ្ស៊ុយលីយេ សង្ខេបរឿង តួអង្គ បញ្ជីណែនាំ

            QGIS export composer to PDF scale the map [closed] Planned maintenance scheduled April 23, 2019 at 23:30 UTC (7:30pm US/Eastern) Announcing the arrival of Valued Associate #679: Cesar Manara Unicorn Meta Zoo #1: Why another podcast?Print Composer QGIS 2.6, how to export image?QGIS 2.8.1 print composer won't export all OpenCycleMap base layer tilesSave Print/Map QGIS composer view as PNG/PDF using Python (without changing anything in visible layout)?Export QGIS Print Composer PDF with searchable text labelsQGIS Print Composer does not change from landscape to portrait orientation?How can I avoid map size and scale changes in print composer?Fuzzy PDF export in QGIS running on macSierra OSExport the legend into its 100% size using Print ComposerScale-dependent rendering in QGIS PDF output

            PDF-ში გადმოწერა სანავიგაციო მენიუproject page