SCL to JSON

geosentinel/L2ASCLtoJSON.py

@@ -6,6 +6,7 @@ import numpy as np
 import gdal
 import sys
 import json
+from collections import Counter
 #-------------------------------------------------------------------------------
 def L2ASCLtoDict(filename, wkt):
 	''' 
@@ -17,14 +18,11 @@ def L2ASCLtoDict(filename, wkt):
 	print ("rows: " + str(rows) + " , cols: " + str(cols))
 	print(len(data))
 	print(len(data[0]))
-	dataDict = dict()
-#	for d in data:
-#		dataClass = str(int(d))
-#		if dataClass in dataDict:
-#			dataDict[dataClass]+=1
-#		else:
-#			dataDict[dataClass]=0
-#	return dataDict
+	count = Counter([])
+	for r in data:
+		count = count + Counter(r)
+	
+	return dict(count)
 #-------------------------------------------------------------------------------
 def main(argv):
 	if len(sys.argv) != 2: