SCL to JSON

geosentinel/L2ASCLtoJSON.py

@@ -14,14 +14,17 @@ def L2ASCLtoDict(filename, wkt):
 	intput = gdal.Open(filename)
 	rows, cols, geotransform = intput.RasterYSize, intput.RasterXSize, intput.GetGeoTransform()
 	data = intput.GetRasterBand(1).ReadAsArray(0,0,cols,rows)
+	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
+#	for d in data:
+#		dataClass = str(int(d))
+#		if dataClass in dataDict:
+#			dataDict[dataClass]+=1
+#		else:
+#			dataDict[dataClass]=0
+#	return dataDict
 #-------------------------------------------------------------------------------
 def main(argv):
 	if len(sys.argv) != 2: