|
@@ -1,36 +1,40 @@
|
|
-from collections import defaultdict
|
|
|
|
import csv
|
|
import csv
|
|
import sys
|
|
import sys
|
|
|
|
+from collections import defaultdict
|
|
|
|
|
|
import cv2
|
|
import cv2
|
|
-from shapely.geometry import MultiPolygon, Polygon
|
|
|
|
-import shapely.wkt
|
|
|
|
-import shapely.affinity
|
|
|
|
-import numpy as np
|
|
|
|
-import tifffile as tiff
|
|
|
|
-
|
|
|
|
#%matplotlib inline
|
|
#%matplotlib inline
|
|
import matplotlib.pyplot as plt
|
|
import matplotlib.pyplot as plt
|
|
|
|
+import numpy as np
|
|
|
|
+import shapely.affinity
|
|
|
|
+import shapely.wkt
|
|
|
|
+# Read and write image data from and to TIFF files.
|
|
|
|
+import tifffile as tiff
|
|
from matplotlib import cm
|
|
from matplotlib import cm
|
|
|
|
+from shapely.geometry import MultiPolygon, Polygon
|
|
|
|
+
|
|
|
|
+import os
|
|
|
|
|
|
-getcwd()
|
|
|
|
-chdir("/media/sf_share/linux/workspaces/guangdong")
|
|
|
|
|
|
+os.getcwd()
|
|
|
|
+os.chdir("/media/sf_share/linux/workspaces/gdzw-analysis")
|
|
|
|
|
|
FILE_2015 = './input/quickbird2015.tif'
|
|
FILE_2015 = './input/quickbird2015.tif'
|
|
FILE_2017 = './input/quickbird2017.tif'
|
|
FILE_2017 = './input/quickbird2017.tif'
|
|
FILE_cadastral2015 = './input/cadastral2015.tif'
|
|
FILE_cadastral2015 = './input/cadastral2015.tif'
|
|
FILE_tinysample = './input/tinysample.tif'
|
|
FILE_tinysample = './input/tinysample.tif'
|
|
|
|
|
|
-#read images
|
|
|
|
|
|
+#read images,return numpy array.
|
|
im_2015 = tiff.imread(FILE_2015).transpose([1, 2, 0])
|
|
im_2015 = tiff.imread(FILE_2015).transpose([1, 2, 0])
|
|
im_2017 = tiff.imread(FILE_2017).transpose([1, 2, 0])
|
|
im_2017 = tiff.imread(FILE_2017).transpose([1, 2, 0])
|
|
im_tiny = tiff.imread(FILE_tinysample)
|
|
im_tiny = tiff.imread(FILE_tinysample)
|
|
im_cada = tiff.imread(FILE_cadastral2015)
|
|
im_cada = tiff.imread(FILE_cadastral2015)
|
|
|
|
|
|
im_2015.shape
|
|
im_2015.shape
|
|
|
|
+# (5106, 15106, 4)
|
|
im_tiny.shape
|
|
im_tiny.shape
|
|
|
|
+# (5106, 15106, 3)
|
|
im_cada.shape
|
|
im_cada.shape
|
|
-
|
|
|
|
|
|
+# (5106, 15106)
|
|
|
|
|
|
def scale_percentile(matrix):
|
|
def scale_percentile(matrix):
|
|
w, h, d = matrix.shape
|
|
w, h, d = matrix.shape
|