Meteoinfo，PSCF2Grid

最新推荐文章于 2024-08-23 07:41:58 发布

Pioneer777

最新推荐文章于 2024-08-23 07:41:58 发布

阅读量492

点赞数

文章标签： python 开发语言

本文链接：https://blog.csdn.net/Pioneer777/article/details/127035653

版权

# Get PSCF layer
mf = miapp.getMapDocument().getActiveMapFrame()
layerName = 'winterPSCF.shp'
layer = mf.getMapView().getLayer(layerName)

# Get grid x/y data
sn = layer.getShapeNum()
xn = 0
oldx = 0.0
for shp in layer.getShapes():
	x = shp.getPoints().get(0).X	
	if xn == 0:
		oldx = x		
	else:
		if x < oldx:
			break
		else:
			oldx = x
	xn += 1

yn = sn / xn
print 'Total number: ' + str(sn)
print 'X number: ' + str(xn)
print 'Y number: ' + str(yn)

idx = sn - xn
shp = layer.getShapes().get(idx)
point = shp.getExtent().getCenterPoint()
minx = point.X
miny = point.Y
rpoint = layer.getShapes().get(idx + 1).getExtent().getCenterPoint()
xdelta = rpoint.X - minx
tpoint = layer.getShapes().get(idx - xn).getExtent().getCenterPoint()
ydelta = tpoint.Y - miny

# Create GridData and set value
print 'Create GridData and set value...'
gData = GridData(minx, xdelta, xn, miny, ydelta, yn)
fieldName = 'WPSCF'
for i in range(0, yn):
	for j in range(0, xn):
		idx = (yn - i - 1) * xn + j
		v = layer.getCellValue(fieldName, idx)
		gData.data[i][j] = v

# Output GridData file
print 'Output GridData file...'
outfn = 'D:/0303HYSPLITDoc/Meteoinfo_out/cs20132019/cluster/winter/pscf.dat'
gData.saveAsSurferASCIIFile(outfn)

print 'Finished!'