2018

README.md

@@ -24,13 +24,16 @@ add GDAL_DATA to env GDAL_DATA=c:\Program Files\GDAL\gdal-data
 2. download imgkap.exe
 http://www.dacust.com/inlandwaters/imgkap/
 
-3. the python script
+3. download pngquant.exe
+https://pngquant.org/
+
+4. the python script
 
 You can find the script also at http://www.wellenvogel.net/software/bsb_from_moac/index.php
 
-4. On Linux Systems you need to compile imgkap by your own - but this is easy (Thanks to Oleg for his description)
+5. On Linux Systems you need to compile imgkap by your own - but this is easy (Thanks to Oleg for his description)
 
-sudo aptitude install git gcc python-gdal libfreeimage-dev
+sudo aptitude install git gcc python-gdal libfreeimage-dev pngquant
 
 git clone https://github.com/nohal/imgkap.git
 

convert_to_bsb.py

@@ -103,10 +103,19 @@ def imgkap():
    else:
       return "imgkap.exe"
 
+def pngquant():
+   sysname, nodename, release, version, machine = os.uname()
+   if sysname == "Linux":
+      return "pngquant"
+   else:
+      return "pngquant.exe"
+
 
 #------------------------------------
 def convertChartListGDAL(chartlist):
-
+
+  imgkapbin=imgkap()
+  pngquantbin=pngquant()
   for chart in chartlist:
 
     dataset = gdal.Open( chart, gdal.GA_ReadOnly )
@@ -133,8 +142,12 @@ def convertChartListGDAL(chartlist):
     (lrlon,lrlat)=gdal.ApplyGeoTransform(geotr,dataset.RasterXSize,dataset.RasterYSize)
     log("raster: ullat=%f ullon=%f lrlat=%f lrlon=%f"% (geotr[3],geotr[0],lrlat,lrlon))
     (lrlon,lrlat,z)=transformer.TransformPoint(lrlon,lrlat,0)
-    cmd="%s %s %f %f %f %f"% (imgkap(),chart,ullat,ullon,lrlat,lrlon)
-    log("running "+cmd)
+    if chart.lower().endswith(".png"):
+      cmd="%s 127 -f --ext .png %s" % (pngquantbin,chart)
+      log("running: %s " % cmd)
+      os.system(cmd) 
+    cmd="%s %s %f %f %f %f"% (imgkapbin,chart,ullat,ullon,lrlat,lrlon)
+    log("running: %s " % cmd)
     os.system(cmd)
 
 #from gdal2tiles.py
@@ -149,7 +162,8 @@ def metersToLonLat(mx,my):
 
 #------------------------------------
 def convertChartListDirect(chartlist):
-
+  imgkapbin=imgkap()
+  pngquantbin=pngquant()
   for chart in chartlist:
     log("direct "+chart)
     f=open(chart, "rb")
@@ -198,8 +212,12 @@ def convertChartListDirect(chartlist):
     log(chart+" raster: ullat=%f ullon=%f lrlat=%f lrlon=%f"% (uly,ulx,lry,lrx))
     (ullon,ullat)=metersToLonLat(ulx, uly)
     (lrlon,lrlat)=metersToLonLat(lrx, lry)
-    cmd="imgkap.exe %s %f %f %f %f"% (chart,ullat,ullon,lrlat,lrlon)
-    log("running "+cmd)
+    if chart.lower().endswith(".png"):
+      cmd="%s 127 -f --ext .png %s" % (pngquantbin,chart)
+      log("running %s" % cmd)
+      os.system(cmd)
+    cmd="%s %s %f %f %f %f"% (imgkapbin,chart,ullat,ullon,lrlat,lrlon)
+    log("running %s" % cmd)
     os.system(cmd)