def grid_to_buffer(self):
from stars.core.DKDEWrapper import call_kde
x = []
y = []
for i in self.query_points:
pt = self.points[i]
x.append(pt[0])
y.append(pt[1])
n = len(self.query_points)
if n < len(self.points):
query_points = range(n)
arr, rows, cols, grad_min, grad_max = call_kde(
n,
x,
y,
#self.query_points,
range(n),
self.extent,
self.bandwidth,
self.cell_size,
self.kernel,
self.color_band,
self.opaque * 2.55)
#from PIL import Image
#Image.fromarray(arr).save("test.png")
self.bmp = wx.BitmapFromBufferRGBA(cols, rows, arr)
#self.bmp.SaveFile("test.png", wx.BITMAP_TYPE_PNG)
self.gradient_color_min = grad_min
self.gradient_color_max = grad_max
def grid_to_buffer(self):
from stars.core.DKDEWrapper import call_kde
x = []
y = []
for i in self.query_points:
pt = self.points[i]
x.append(pt[0])
y.append(pt[1])
n = len(self.query_points)
if n < len(self.points):
query_points = range(n)
arr,rows,cols,grad_min,grad_max = call_kde(
n,
x,
y,
#self.query_points,
range(n),
self.extent,
self.bandwidth,
self.cell_size,
self.kernel,
self.color_band,
self.opaque*2.55
)
#from PIL import Image
#Image.fromarray(arr).save("test.png")
self.bmp = wx.BitmapFromBufferRGBA(cols, rows, arr)
#self.bmp.SaveFile("test.png", wx.BITMAP_TYPE_PNG)
self.gradient_color_min = grad_min
self.gradient_color_max = grad_max
pos3 = line.rfind(' ')
pos4 = line.find(')')
extent.append(float(line[4:pos1]))
extent.append(float(line[pos1+1:pos2]))
extent.append(float(line[pos2+1:pos3]))
extent.append(float(line[pos3+1:pos4]))
f2.close()
"""
if cellsize == "":
w = extent[2] - extent[0]
h = extent[1] - extent[3]
image_w = 600
image_h = h * image_w / w
cellsize = max(w / image_w, h / image_h)
std_x = std(x)
std_y = std(y)
Q_x = quantile(x)
Q_y = quantile(y)
IQR_x = Q_x[2] - Q_x[0]
IQR_y = Q_y[2] - Q_y[0]
h_x = 0.9 * min(std_x, IQR_x / 1.34) * pow(n, -0.2)
h_y = 0.9 * min(std_y, IQR_y / 1.34) * pow(n, -0.2)
bandwidth = max(h_x, h_y) * 2
else:
cellsize = float(cellsize)
bandwidth = float(bandwidth)
itv_pts_ids = range(n)
arr, rows, cols, gmin, gmax = call_kde(n, x, y, itv_pts_ids, extent, bandwidth, cellsize, kernel, gradient, opaque)
fromarray(arr).save("test1.png")
pos4 = line.find(')')
extent.append(float(line[4:pos1]))
extent.append(float(line[pos1+1:pos2]))
extent.append(float(line[pos2+1:pos3]))
extent.append(float(line[pos3+1:pos4]))
f2.close()
"""
if cellsize == "":
w = extent[2] - extent[0]
h = extent[1] - extent[3]
image_w = 600
image_h = h * image_w / w
cellsize = max(w / image_w, h / image_h)
std_x = std(x)
std_y = std(y)
Q_x = quantile(x)
Q_y = quantile(y)
IQR_x = Q_x[2] - Q_x[0]
IQR_y = Q_y[2] - Q_y[0]
h_x = 0.9 * min(std_x, IQR_x / 1.34) * pow(n, -0.2)
h_y = 0.9 * min(std_y, IQR_y / 1.34) * pow(n, -0.2)
bandwidth = max(h_x, h_y) * 2
else:
cellsize = float(cellsize)
bandwidth = float(bandwidth)
itv_pts_ids = range(n)
arr, rows, cols, gmin, gmax = call_kde(n, x, y, itv_pts_ids, extent, bandwidth,
cellsize, kernel, gradient, opaque)
fromarray(arr).save("test1.png")