私は次の機能を使用しています:
def LAS2TXTGridClip(inFile,poly,MinPoints=1):
sf = shapefile.Reader(poly) #open shpfile
sr = sf.shapeRecords()
poly_filename, ext = path.splitext(poly)
inFile_filename = os.path.splitext(os.path.basename(inFile))[0]
for i in xrange(len(sr)):
verts = np.array(sr[i].shape.points,float)
record = sr[i].record[0]
inside_points = [p for p in lasfile.File(inFile,None,'r') if pnpoly(p.x, p.y, verts)]
if len(inside_points) >= MinPoints:
file_out = open("{0}_{1}_{2}.txt".format(poly_filename, inFile_filename, record), "w")
for p in inside_points:
file_out.write("%s %s %s %s %s %s %s %s %s %s %s" % (p.x, p.y, p.z, p.intensity,p.return_number,p.number_of_returns,p.scan_direction,p.flightline_edge,p.classification,p.scan_angle,record)+ "\n")
file_out.close()
for i in xrange(len(sr)):
関数は数回処理されます。はlen(sr)
約 50 万で、待つ必要がある時間を把握するためにプログレス バーを挿入したいと考えています (金曜日です)。次の質問があります。
- Windows OS 64ビット上のpython 27の「最も簡単な」プログレスバーはどれですか?
- プログレスバー モジュールを見つけましたが、この手順の後で easy_install プログレスバーを使用するのに問題があります。
- プログレスバーを挿入するのに最適な位置はどこですか?