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;
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
add a comment |
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
add a comment |
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
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
python gdal polygonize multi-band
asked Apr 11 at 16:42
AAAAAA
111
111
add a comment |
add a comment |
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
);
);
Sign up or log in
StackExchange.ready(function ()
StackExchange.helpers.onClickDraftSave('#login-link');
);
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Post as a guest
Required, but never shown
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
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.
Sign up or log in
StackExchange.ready(function ()
StackExchange.helpers.onClickDraftSave('#login-link');
);
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Post as a guest
Required, but never shown
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
Sign up or log in
StackExchange.ready(function ()
StackExchange.helpers.onClickDraftSave('#login-link');
);
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Post as a guest
Required, but never shown
Sign up or log in
StackExchange.ready(function ()
StackExchange.helpers.onClickDraftSave('#login-link');
);
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Post as a guest
Required, but never shown
Sign up or log in
StackExchange.ready(function ()
StackExchange.helpers.onClickDraftSave('#login-link');
);
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
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