Hi,

Sorry to answer that now I have not had the notification ... (it was in my spam)

I used software to binarize my png which allowed me to carry out my project 
until the end.

But I like the idea of being able to process this via a script.
I have a script_list.sh which contains a script_file.sh for each tile which 
allows me to automate the tasks for the 512 tiles.
I will test your solution and incorporate it as soon as I have time (I want to 
finalize the last things I have left before) and I will keep you posted.

I'm on mac but I installed gawk and gsed (awk and sed are different on macOS)

Thanks again for the code on github ;)

Le 3 janv. 2020 à 01:42, 
[email protected]<mailto:[email protected]> a écrit :

Hi,

There is one way to solve this with gdal_calc. I have had some time off during 
the holidays,
but kept thinking about your question and thought it was a fun problem.
Actually I have solved some other problems with tif palette images and 
gdal_calc lately. So I put together one solution that
I think will work in your case. I am assuming that You are using linux and have 
some knowledge in bash scripting.

First we have to create a reclass rule for gdal_calc, We do that with with 
running gdalinfo.
Strip the info to only keep the color table with sed.
Then dynamically create the rule with awk, based on your rgb values (turn 
palette values into either 0 or 1).
Finally strip the last "+" from the rule (the last sed command)

Then run gdal_calc with your relclass rule

gdal_calc produces a grey scale image, but:

We then turn the image into a virtual raster file (VRT)
We create a new color table and insert it into the vrt file with sed.
Change the vrt file to be a palette image again with sed.
Finally we can turn the vrt image into a tif image again.

If this works you could loop over all your images and process all in the same 
way.
Just do a bash function with the steps and loop over your input files.

The whole idea of the process is outlined in a gist I just uploaded.
I just tested the first part with the rule myself, but the other stuff I used 
before
so it should work.

https://gist.github.com/LarsSchy/c1be96fa03f7dc11df4247be033e2705

I hope this helps You.

Have fun / Lars Schylberg

_______________________________________________
gdal-dev mailing list
[email protected]
https://lists.osgeo.org/mailman/listinfo/gdal-dev

Reply via email to