截取输入图片的户型图部分（可用于其他图片中物体部分截取）

在输入的户型图中，一般含有其他无关部分，如下图：

而我们在做户型图处理的时候，只需要取其中户型图相关部分：
需要使用openCV来做截取：

# step1：加载图片，转成灰度图
import cv2
import numpy as np
image = cv2.imread("./pic/ceshi4.jpg")
gray = cv2.cvtColor(image, cv2.COLOR_BGR2GRAY)

#step2:用Sobel算子计算x，y方向上的梯度，之后在x方向上减去y方向上的梯度，通过这个减法，我们留下具有高水平梯度和低垂直梯度的图像区域。
gradX = cv2.Sobel(gray, cv2.CV_32F, dx=1, dy=0, ksize=-1)
gradY = cv2.Sobel(gray, cv2.CV_32F, dx=0, dy=1, ksize=-1)
gradient = cv2.subtract(gradX, gradY)
gradient = cv2.convertScaleAbs(gradient)
cv2.imshow("first", gradient)
cv2.waitKey()

# step3：去除图像上的噪声。首先使用低通滤泼器平滑图像（9 x 9内核）,这将有助于平滑图像中的高频噪声。
# 低通滤波器的目标是降低图像的变化率。如将每个像素替换为该像素周围像素的均值。这样就可以平滑并替代那些强度变化明显的区域。
# 然后，对模糊图像二值化。梯度图像中不大于90的任何像素都设置为0（黑色）。 否则，像素设置为255（白色）。
blurred = cv2.blur(gradient, (9, 9))
_, thresh = cv2.threshold(blurred, 90, 255, cv2.THRESH_BINARY)
# SHOW IMAGE
cv2.imshow("thresh", thresh)
cv2.waitKey()

# step4:在上图中我们看到户型图区域有很多黑色的空余，我们要用白色填充这些空余，使得后面的程序更容易识别户型图区域，
# 这需要做一些形态学方面的操作。
kernel = cv2.getStructuringElement(cv2.MORPH_RECT, (25, 25))
closed = cv2.morphologyEx(thresh, cv2.MORPH_CLOSE, kernel)
# show image
cv2.imshow("closed1", closed)
cv2.waitKey()

# step5:从上图我们发现图像上还有一些小的白色斑点，这会干扰之后的户型图轮廓的检测，要把它们去掉。分别执行4次形态学腐蚀与膨胀。
closed = cv2.erode(closed, None, iterations=4)
closed = cv2.dilate(closed, None, iterations=4)
# show image
cv2.imshow("closed2", closed)
cv2.waitKey()

#step6：找出户型图区域的轮廓。

#cv2.findContours()函数返回两个值，一个是轮廓本身，还有一个是每条轮廓对应的属性。
#cv2.findContours()函数返回第一个值是list，list中每个元素都是图像中的一个轮廓，用numpy中的ndarray表示。
#每一个ndarray里保存的是轮廓上的各个点的坐标。我们把list排序，点最多的那个轮廓就是我们要找的户型图的轮廓。

# cv2.minAreaRect()函数:
# 主要求得包含点集最小面积的矩形，这个矩形是可以有偏转角度的，可以与图像的边界不平行。

(cnts, _) = cv2.findContours(closed.copy(), cv2.RETR_EXTERNAL, cv2.CHAIN_APPROX_SIMPLE)
c = sorted(cnts, key=cv2.contourArea, reverse=True)[0]
rect = cv2.minAreaRect(c)
box = np.int0(cv2.boxPoints(rect))
cv2.imshow("Image", image)
cv2.imwrite("contoursImage2.jpg", image)
cv2.waitKey(0)

#step7：裁剪。box里保存的是绿色矩形区域四个顶点的坐标。
#我将按下图红色矩形所示裁剪昆虫图像。找出四个顶点的x，y坐标的最大最小值
#新图像的高=maxY-minY，宽=maxX-minX。

Xs = [i[0] for i in box]
Ys = [i[1] for i in box]
x1 = min(Xs)
x2 = max(Xs)
y1 = min(Ys)
y2 = max(Ys)
hight = y2 - y1
width = x2 - x1
cropImg = image[y1:y1+hight, x1:x1+width]
print(y1)
print(y1+hight)
print(x1)
print(x1+width)
cv2.imshow('cropImg',cropImg)
cv2.imwrite("mff.jpg", image)
cv2.waitKey(0)

经过上述操作后，可以得到：

可以看到，只截取了我们需要的部分。