0

我正在研究 QGIS 中的插件，该插件从栅格数据创建 3d 模型。在这种情况下，我只需要在垂直平面中创建三角剖分（没有将垂直面与相邻点连接）。

这是我的代码和具有当前效果的示例。

    def delaunay_ogolny(self):
    import math
    X=[]
    Y=[]
    Z=[]
    for x in self.listqa:
        dada=x[0]
        X.append(dada)
        uu=x[1]
        Y.append(uu)
        dede=x[2]
        Z.append(dede)
    u = np.array(X)
    v = np.array(Y)
    x = u
    y = v
    z = np.array(Z)
    #Delaunay
    tri = Delaunay(np.array([u, v]).T,qhull_options="QJn")    #,furthest_site=True   #i put there QJn option
    #tri=Delaunay(self.listqa)
    self.faces = []
    self.points = []


    for vert in tri.simplices:
            self.faces.append([vert[0], vert[1], vert[2]])
    #print(self.faces)

    for i in range(x.shape[0]):
            self.points.append([x[i], y[i], z[i]])
    #print(self.points)
    return self.faces, self.points

在 Delaunay 函数之后，我从点和面创建 stl 文件。以下是我需要的，但该算法仅在我得到一条直线时才能正常工作。我需要它的所有点。

在此处输入图像描述

EDIT1：@Ripi2

我试过你的想法，但没有奏效。

这是我的迭代函数代码，我迭代所有像素，当像素获得值时，我的算法将其添加到列表中。之后，我在 Delaunay 三角剖分中使用此列表。

这是代码：

    def iterator_shape(self):
    #selectedLayerIndex = self.dlg.comboBox.currentIndex()
    path=os.path.join(r'default\python\plugins\print_your_3d\TRASH\shape_to_raster.tif')
    layer = iface.addRasterLayer(sciezka, "qtaz")
    provider = layer.dataProvider()
    extent = provider.extent()
    rows = layer.height()
    cols = layer.width()
    block = provider.block(1, extent, cols, rows)
    xmin = extent.xMinimum()
    ymax = extent.yMaximum()
    xsize = layer.rasterUnitsPerPixelX()
    ysize = layer.rasterUnitsPerPixelY()
    k=1
    xinit = xmin + xsize / 2
    yinit = ymax - ysize / 2
    self.listqa = []
    for i in range(rows):
        for j in range(cols):
            x = xinit + j * xsize
            y = yinit
            k += 1
            if block.value(i, j)  == -3.4028234663852886e+38:
                continue
            elif block.value(i,j)>=self.dlg.doubleSpinBox.value():
                self.listqa.append([x, y, block.value(i, j)])
                self.listqa.append([x,y,block.value(i,j)-1])
                self.listqa.append([x,y+0.01,block.value(i,j)-1])   # @Ridi2 you mean that? 
        xinit = xmin + xsize / 2
        yinit -= ysize
    return self.listqa

这是我想转换为 3d STL 模型的栅格示例。1 个像素有 X、Y 和高程值（在代码中这是 block.value）

在此处输入图像描述

Edit2：我需要在这之间做鬼脸。这是在 ArcScene 中进行的更好的可视化。 在此处输入图像描述

4

0 回答 0