2016-03-08 112 views
1

有一個代碼:爲什麼需要浮標?

from arcpy.sa import * 

arcpy.CheckOutExtension("Spatial") 


print "Creating and defining X and Y coordinates:" #set env properties 
env.workspace = r"C:\Users\Desktop\data\new" 
env.overwriteOutput = 1 


theme = '/watersheds_3D.shp' #polygons in vector format 

demName = '/demlab4' # raster of DEM 
rasObject = Raster(demName) 


my_extent = rasObject.extent #find a grid cell size: 
my_cellsize = (rasObject.meanCellHeight + rasObject.meanCellWidth)/2 
print my_cellsize 

slope_deg = Slope(demName) 
slope_deg.save('/demlab4_slope') 

for row in arcpy.da.SearchCursor(theme, ["FID","[email protected]"]): # inside shp file - 
    print ("Polygon # {}: ".format(row[0])) 

    for part in row[1]: # Inside the polygon 

     xCoords = [] 
     yCoords = [] 
     zCoords = [] 
     for pnt in part: 

      print (" {}, {}, {} ".format(pnt.X, pnt.Y, pnt.Z)) 
      xCoords.append(pnt.X) 
      yCoords.append(pnt.Y) 
      zCoords.append(pnt.Z) 

     area3D = lab4_arcpy_module.define3Darea(my_cellsize, slope_deg) 

和函數來計算的:

def define3Darea(my_cellsize, slope_deg): 

    a = my_cellsize 
    slope_rad = math.pi * slope_deg/180 # return in radiant 
    c = math.sqrt(a**2 + (math.tan(slope_rad)*a)**2) 

    area3D = a*c 
    print "Total 3D Area is: ", area3D, "m^2" 

    return area3D 

它給出:類型錯誤:需要的浮子。 什麼是問題?

還需要計算每個多邊形(4)內部的3D區域。怎麼做?

+0

您確實需要提供完整的錯誤消息,以便我們可以看到發生錯誤的位置。上下文就是一切! –

回答

0

不確定RasterSlope正在做什麼,但是當您嘗試對非浮點類型的變量執行某些浮點操作時會出現此TypeError: a float is required。我的猜測會是這樣的:slope_deg變量。請看這裏:

slope_deg = Slope(demName) 
slope_deg.save('/demlab4_slope') 

你正在聲明和做一些與slope_deg變量的操作。現在,在此聲明:

area3D = lab4_arcpy_module.define3Darea(my_cellsize, slope_deg) 

您正在使用相同的變量,它是Slope類的實例,而不是一個浮點型變量。或者必須有一些屬性(Slope),它是float或者您必須使用正確的變量。

0

我發現那裏有一些錯誤試圖打印虧損如此:

print ("Model paramters:")  
print ("Weight:%f" %sess.run(W)) 
print ("bias:%f" %sess.run(b)) 
#print ("loss:%f" %(loss)) <------ 

使誤差disapear,還是不明白,也這是爲什麼。