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

          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