Skip to content

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

merge-crop

parent 22171225
Hide whitespace changes
Inline Side-by-side
Showing with 10 additions and 10 deletions
geosentinel/L2ASCLtoJSON.py
View file @ 2df33735
...@@ -12,16 +12,16 @@ def L2ASCLtoDict(filename, wkt): ...@@ -12,16 +12,16 @@ def L2ASCLtoDict(filename, wkt):
L2ASCLtoDict L2ASCLtoDict
''' '''
intput = gdal.Open(filename) intput = gdal.Open(filename)
rows, cols, geotransform = intput.RasterYSize, intput.RasterXSize, intput.GetGeoTransform() # rows, cols, geotransform = intput.RasterYSize, intput.RasterXSize, intput.GetGeoTransform()
data = intput.GetRasterBand(0).ReadAsArray(0,0,cols,rows) # data = intput.GetRasterBand(0).ReadAsArray(0,0,cols,rows)
dataDict = dict() # dataDict = dict()
for d in data: # for d in data:
dataClass = str(int(d)) # dataClass = str(int(d))
if dataClass in dataDict: # if dataClass in dataDict:
dataDict[dataClass]+=1 # dataDict[dataClass]+=1
else: # else:
dataDict[dataClass]=0 # dataDict[dataClass]=0
return dataDict # return dataDict
#------------------------------------------------------------------------------- #-------------------------------------------------------------------------------
def main(argv): def main(argv):
if len(sys.argv) != 2: 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