我试图用matplotlib来绘制一些我正在研究的文件的数字。我有二组数据在二维数组:一个ascii山体阴影栅格,我可以很高兴地绘制和调整使用:
import matplotlib.pyplot as pp
import numpy as np
hillshade = np.genfromtxt('hs.asc', delimiter=' ', skip_header=6)[:,:-1]
pp.imshow(hillshade, vmin=0, vmax=255)
pp.gray()
pp.show()
这使:
还有一个第二个ascii光栅,描绘了横跨风景的河流的属性。该数据可以以与上述相同的方式绘制,但是不对应于河网的阵列中的值被指定为-9999的无数据值。目的是使无数据值设置为透明,以便河流值覆盖山体阴影。
这是河流数据,理想情况下,这里代表的每个像素为0将是完全透明的。
对此做了一些研究,似乎我可以将我的数据转换成RGBA数组,并将alpha值设置为仅使不需要的单元格透明。然而,河流数组中的值是浮动的,不能被转换(因为原始值是图的整个点),我相信如果使用RGBA格式,imshow函数只能取无符号整数。
有没有办法解决这个限制?我曾经希望我可以用像素值和alpha值创建一个元组,并将它们绘制成这样,但这似乎不可能。
我也和PIL一起玩,试图创建一个没有数据值透明的河流数据的PNG文件,但这似乎自动将像素值缩放到0-255,从而失去了我需要保留的值。
我欢迎任何人对这个问题的任何见解。