GDAL - Polygonize multiband raster and save pixel values of each band to corresponding attribute fields 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?GDAL polygonize only certain raster valuesPolygonize a raster using GDAL JAVA bindings - passing raster pixel value to the resulting featureWhy does the composite raster pixel values not cover the full range of RGB and what does this mean for reclassification after conversion to int?Obtain coordinates and corresponding pixel values from GeoTiff using python gdal and save them as numpy arrayExtract All Pixel Values for each band in a Multi Band RasterMultiband to single band rasterSave multiband raster with gdalPolygonize a multi-band rasterGdal Polygonize, empty attribute tableHow to polygonize raster file according to band values?

What is an "asse" in Elizabethan English?

Maximum summed subsequences with non-adjacent items

Has negative voting ever been officially implemented in elections, or seriously proposed, or even studied?

What is the chair depicted in Cesare Maccari's 1889 painting "Cicerone denuncia Catilina"?

What makes a man succeed?

Central Vacuuming: Is it worth it, and how does it compare to normal vacuuming?

Sum letters are not two different

Do I really need to have a message in a novel to appeal to readers?

Does the Mueller report show a conspiracy between Russia and the Trump Campaign?

Amount of permutations on an NxNxN Rubik's Cube

How can I set the aperture on my DSLR when it's attached to a telescope instead of a lens?

macOS: Name for app shortcut screen found by pinching with thumb and three fingers

How to write capital alpha?

What does Turing mean by this statement?

What's the point of the test set?

How would a mousetrap for use in space work?

Converted a Scalar function to a TVF function for parallel execution-Still running in Serial mode

Why are my pictures showing a dark band on one edge?

Draw 4 of the same figure in the same tikzpicture

Tannaka duality for semisimple groups

1-probability to calculate two events in a row

Is there any word for a place full of confusion?

AppleTVs create a chatty alternate WiFi network

Drawing spherical mirrors



GDAL - Polygonize multiband raster and save pixel values of each band to corresponding attribute fields



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?GDAL polygonize only certain raster valuesPolygonize a raster using GDAL JAVA bindings - passing raster pixel value to the resulting featureWhy does the composite raster pixel values not cover the full range of RGB and what does this mean for reclassification after conversion to int?Obtain coordinates and corresponding pixel values from GeoTiff using python gdal and save them as numpy arrayExtract All Pixel Values for each band in a Multi Band RasterMultiband to single band rasterSave multiband raster with gdalPolygonize a multi-band rasterGdal Polygonize, empty attribute tableHow to polygonize raster file according to band values?



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








1















I'm trying to polygonize 2 band raster(which was previously reclassified) and put values of each band into seperate attribute table column. Both bands overlap each other perfectly.



This piece of code works perfectly for 1 band rasters:



shapeFile = "shape1.shp"
driver = ogr.GetDriverByName("ESRI Shapefile")
outDS = driver.CreateDataSource(shapeFile)
outLayer = outDS.CreateLayer("lyr1", srs)
f1 = ogr.FieldDefn('DN', ogr.OFTInteger)
f2 = ogr.FieldDefn('DIF', ogr.OFTInteger)
outLayer.CreateField(f1)
outLayer.CreateField(f2)
gdal.Polygonize(raster.dataset.GetRasterBand(1), rasterMask, outLayer, 0, [], callback=None )
#gdal.Polygonize(raster.dataset.GetRasterBand(2), rasterMask, outLayer, 1, [], callback=None )


However the commented piece of code where I try to put pixel values of 2nd band into the 'DIF' field doesn't work (and it's field index is 1, I've checked that).



The Polygnize function doesn't seem to have an option for multiple band rasters and I'm not sure if calling it 2 times in a row is even a good idea.



I was thinking of polygonizing each band to a different shapefile and then somehow merging those shapefiles to match each field value but I'm not even sure if that's possible.



How should I proceed?










share|improve this question




























    1















    I'm trying to polygonize 2 band raster(which was previously reclassified) and put values of each band into seperate attribute table column. Both bands overlap each other perfectly.



    This piece of code works perfectly for 1 band rasters:



    shapeFile = "shape1.shp"
    driver = ogr.GetDriverByName("ESRI Shapefile")
    outDS = driver.CreateDataSource(shapeFile)
    outLayer = outDS.CreateLayer("lyr1", srs)
    f1 = ogr.FieldDefn('DN', ogr.OFTInteger)
    f2 = ogr.FieldDefn('DIF', ogr.OFTInteger)
    outLayer.CreateField(f1)
    outLayer.CreateField(f2)
    gdal.Polygonize(raster.dataset.GetRasterBand(1), rasterMask, outLayer, 0, [], callback=None )
    #gdal.Polygonize(raster.dataset.GetRasterBand(2), rasterMask, outLayer, 1, [], callback=None )


    However the commented piece of code where I try to put pixel values of 2nd band into the 'DIF' field doesn't work (and it's field index is 1, I've checked that).



    The Polygnize function doesn't seem to have an option for multiple band rasters and I'm not sure if calling it 2 times in a row is even a good idea.



    I was thinking of polygonizing each band to a different shapefile and then somehow merging those shapefiles to match each field value but I'm not even sure if that's possible.



    How should I proceed?










    share|improve this question
























      1












      1








      1








      I'm trying to polygonize 2 band raster(which was previously reclassified) and put values of each band into seperate attribute table column. Both bands overlap each other perfectly.



      This piece of code works perfectly for 1 band rasters:



      shapeFile = "shape1.shp"
      driver = ogr.GetDriverByName("ESRI Shapefile")
      outDS = driver.CreateDataSource(shapeFile)
      outLayer = outDS.CreateLayer("lyr1", srs)
      f1 = ogr.FieldDefn('DN', ogr.OFTInteger)
      f2 = ogr.FieldDefn('DIF', ogr.OFTInteger)
      outLayer.CreateField(f1)
      outLayer.CreateField(f2)
      gdal.Polygonize(raster.dataset.GetRasterBand(1), rasterMask, outLayer, 0, [], callback=None )
      #gdal.Polygonize(raster.dataset.GetRasterBand(2), rasterMask, outLayer, 1, [], callback=None )


      However the commented piece of code where I try to put pixel values of 2nd band into the 'DIF' field doesn't work (and it's field index is 1, I've checked that).



      The Polygnize function doesn't seem to have an option for multiple band rasters and I'm not sure if calling it 2 times in a row is even a good idea.



      I was thinking of polygonizing each band to a different shapefile and then somehow merging those shapefiles to match each field value but I'm not even sure if that's possible.



      How should I proceed?










      share|improve this question














      I'm trying to polygonize 2 band raster(which was previously reclassified) and put values of each band into seperate attribute table column. Both bands overlap each other perfectly.



      This piece of code works perfectly for 1 band rasters:



      shapeFile = "shape1.shp"
      driver = ogr.GetDriverByName("ESRI Shapefile")
      outDS = driver.CreateDataSource(shapeFile)
      outLayer = outDS.CreateLayer("lyr1", srs)
      f1 = ogr.FieldDefn('DN', ogr.OFTInteger)
      f2 = ogr.FieldDefn('DIF', ogr.OFTInteger)
      outLayer.CreateField(f1)
      outLayer.CreateField(f2)
      gdal.Polygonize(raster.dataset.GetRasterBand(1), rasterMask, outLayer, 0, [], callback=None )
      #gdal.Polygonize(raster.dataset.GetRasterBand(2), rasterMask, outLayer, 1, [], callback=None )


      However the commented piece of code where I try to put pixel values of 2nd band into the 'DIF' field doesn't work (and it's field index is 1, I've checked that).



      The Polygnize function doesn't seem to have an option for multiple band rasters and I'm not sure if calling it 2 times in a row is even a good idea.



      I was thinking of polygonizing each band to a different shapefile and then somehow merging those shapefiles to match each field value but I'm not even sure if that's possible.



      How should I proceed?







      python gdal polygonize multi-band






      share|improve this question













      share|improve this question











      share|improve this question




      share|improve this question










      asked Apr 11 at 16:42









      AAAAAA

      111




      111




















          0






          active

          oldest

          votes












          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%2f318519%2fgdal-polygonize-multiband-raster-and-save-pixel-values-of-each-band-to-corresp%23new-answer', 'question_page');

          );

          Post as a guest















          Required, but never shown

























          0






          active

          oldest

          votes








          0






          active

          oldest

          votes









          active

          oldest

          votes






          active

          oldest

          votes















          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%2f318519%2fgdal-polygonize-multiband-raster-and-save-pixel-values-of-each-band-to-corresp%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

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

          Crop image to path created in TikZ? Announcing the arrival of Valued Associate #679: Cesar Manara Planned maintenance scheduled April 17/18, 2019 at 00:00UTC (8:00pm US/Eastern)Crop an inserted image?TikZ pictures does not appear in posterImage behind and beyond crop marks?Tikz picture as large as possible on A4 PageTransparency vs image compression dilemmaHow to crop background from image automatically?Image does not cropTikzexternal capturing crop marks when externalizing pgfplots?How to include image path that contains a dollar signCrop image with left size given

          Romeo and Juliet ContentsCharactersSynopsisSourcesDate and textThemes and motifsCriticism and interpretationLegacyScene by sceneSee alsoNotes and referencesSourcesExternal linksNavigation menu"Consumer Price Index (estimate) 1800–"10.2307/28710160037-3222287101610.1093/res/II.5.31910.2307/45967845967810.2307/2869925286992510.1525/jams.1982.35.3.03a00050"Dada Masilo: South African dancer who breaks the rules"10.1093/res/os-XV.57.1610.2307/28680942868094"Sweet Sorrow: Mann-Korman's Romeo and Juliet Closes Sept. 5 at MN's Ordway"the original10.2307/45957745957710.1017/CCOL0521570476.009"Ram Leela box office collections hit massive Rs 100 crore, pulverises prediction"Archived"Broadway Revival of Romeo and Juliet, Starring Orlando Bloom and Condola Rashad, Will Close Dec. 8"Archived10.1075/jhp.7.1.04hon"Wherefore art thou, Romeo? To make us laugh at Navy Pier"the original10.1093/gmo/9781561592630.article.O006772"Ram-leela Review Roundup: Critics Hail Film as Best Adaptation of Romeo and Juliet"Archived10.2307/31946310047-77293194631"Romeo and Juliet get Twitter treatment""Juliet's Nurse by Lois Leveen""Romeo and Juliet: Orlando Bloom's Broadway Debut Released in Theaters for Valentine's Day"Archived"Romeo and Juliet Has No Balcony"10.1093/gmo/9781561592630.article.O00778110.2307/2867423286742310.1076/enst.82.2.115.959510.1080/00138380601042675"A plague o' both your houses: error in GCSE exam paper forces apology""Juliet of the Five O'Clock Shadow, and Other Wonders"10.2307/33912430027-4321339124310.2307/28487440038-7134284874410.2307/29123140149-661129123144728341M"Weekender Guide: Shakespeare on The Drive""balcony"UK public library membership"romeo"UK public library membership10.1017/CCOL9780521844291"Post-Zionist Critique on Israel and the Palestinians Part III: Popular Culture"10.2307/25379071533-86140377-919X2537907"Capulets and Montagues: UK exam board admit mixing names up in Romeo and Juliet paper"Istoria Novellamente Ritrovata di Due Nobili Amanti2027/mdp.390150822329610820-750X"GCSE exam error: Board accidentally rewrites Shakespeare"10.2307/29176390149-66112917639"Exam board apologises after error in English GCSE paper which confused characters in Shakespeare's Romeo and Juliet""From Mariotto and Ganozza to Romeo and Guilietta: Metamorphoses of a Renaissance Tale"10.2307/37323537323510.2307/2867455286745510.2307/28678912867891"10 Questions for Taylor Swift"10.2307/28680922868092"Haymarket Theatre""The Zeffirelli Way: Revealing Talk by Florentine Director""Michael Smuin: 1938-2007 / Prolific dance director had showy career"The Life and Art of Edwin BoothRomeo and JulietRomeo and JulietRomeo and JulietRomeo and JulietEasy Read Romeo and JulietRomeo and Julieteeecb12003684p(data)4099369-3n8211610759dbe00d-a9e2-41a3-b2c1-977dd692899302814385X313670221313670221