Add read_geotiff

LiCSBAS_lib/LiCSBAS_io_lib.py

@@ -6,8 +6,10 @@
 Python3 library of input/output functions for LiCSBAS.
 
 =========
-Changelog
+Change log
 =========
+v1.4 20210224 Yu Morioshita, GSI
+ - Add read_geotiff
 v1.3 20210209 Yu Morioshita, GSI
  - Add make_geotiff
 v1.2.1 20201211 Yu Morioshita, GSI
@@ -166,6 +168,26 @@ def read_bperp_file(bperp_file, imdates):
  return bperp
 
 
+#%%
+def read_geotiff(file, file_ref=None):
+ geotiff = gdal.Open(file)
+
+ if file_ref is not None: # Compare size and area
+ size = (geotiff.RasterXSize, geotiff.RasterYSize)
+ area = geotiff.GetGeoTransform()
+
+ geotiff_ref = gdal.Open(file_ref)
+ size_ref = (geotiff_ref.RasterXSize, geotiff_ref.RasterYSize)
+ area_ref = geotiff_ref.GetGeoTransform()
+
+ if not (size == size_ref and area == area_ref):
+ raise Exception('ERROR: File size or area are not identical between {} and {}'.format(file, file_ref))
+
+ data = geotiff.ReadAsArray()
+
+ return data
+
+
 #%%
 def read_img(file, length, width, dtype=np.float32, endian='little'):
  """