Skip to content

GeoSentinel
GeoSentinel
Commit cb3ca00d authored by Mario Chirinos Colunga's avatar Mario Chirinos Colunga 💬
Browse files
Options

merge-crop

parent 2f4bac91
Show whitespace changes
Inline Side-by-side
Showing with 10 additions and 10 deletions
geosentinel/L2ASCLtoJSON.py
View file @ cb3ca00d
......@@ -12,16 +12,16 @@ def L2ASCLtoDict(filename, wkt):
L2ASCLtoDict
'''
intput = gdal.Open(filename)
# rows, cols, geotransform = intput.RasterYSize, intput.RasterXSize, intput.GetGeoTransform()
# data = intput.GetRasterBand(0).ReadAsArray(0,0,cols,rows)
# dataDict = dict()
# for d in data:
# dataClass = str(int(d))
# if dataClass in dataDict:
# dataDict[dataClass]+=1
# else:
# dataDict[dataClass]=0
# return dataDict
rows, cols, geotransform = intput.RasterYSize, intput.RasterXSize, intput.GetGeoTransform()
data = intput.GetRasterBand(0).ReadAsArray(0,0,cols,rows)
dataDict = dict()
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:
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment