def map(self, L, x, y): dd = 1.0 / (2.0**(int(L))) relx = int(x) * dd rely = int(y) * dd DW = (self.ds.domain_right_edge - self.ds.domain_left_edge) xl = self.ds.domain_left_edge[0] + relx * DW[0] yl = self.ds.domain_left_edge[1] + rely * DW[1] xr = xl + dd * DW[0] yr = yl + dd * DW[1] frb = FixedResolutionBuffer(self.data, (xl, xr, yl, yr), (256, 256)) cmi, cma = get_color_bounds( self.data['px'], self.data['py'], self.data['pdx'], self.data['pdy'], self.data[self.field], self.ds.domain_left_edge[0], self.ds.domain_right_edge[0], self.ds.domain_left_edge[1], self.ds.domain_right_edge[1], dd * DW[0] / (64 * 256), dd * DW[0]) if self.ds._get_field_info(self.field).take_log: cmi = np.log10(cmi) cma = np.log10(cma) to_plot = apply_colormap(np.log10(frb[self.field]), color_bounds=(cmi, cma)) else: to_plot = apply_colormap(frb[self.field], color_bounds=(cmi, cma)) rv = write_png_to_string(to_plot) return rv
def map(self, L, x, y): dd = 1.0 / (2.0**(int(L))) relx = int(x) * dd rely = int(y) * dd DW = (self.ds.domain_right_edge - self.ds.domain_left_edge) xl = self.ds.domain_left_edge[0] + relx * DW[0] yl = self.ds.domain_left_edge[1] + rely * DW[1] xr = xl + dd*DW[0] yr = yl + dd*DW[1] frb = FixedResolutionBuffer(self.data, (xl, xr, yl, yr), (256, 256)) cmi, cma = get_color_bounds(self.data['px'], self.data['py'], self.data['pdx'], self.data['pdy'], self.data[self.field], self.ds.domain_left_edge[0], self.ds.domain_right_edge[0], self.ds.domain_left_edge[1], self.ds.domain_right_edge[1], dd*DW[0] / (64*256), dd*DW[0]) if self.ds._get_field_info(self.field).take_log: cmi = np.log10(cmi) cma = np.log10(cma) to_plot = apply_colormap(np.log10(frb[self.field]), color_bounds = (cmi, cma)) else: to_plot = apply_colormap(frb[self.field], color_bounds = (cmi, cma)) rv = write_png_to_string(to_plot) return rv
def map(self, field, L, x, y): if "," in field: field = tuple(field.split(",")) cmap = self.cmap dd = 1.0 / (2.0**(int(L))) relx = int(x) * dd rely = int(y) * dd DW = self.ds.domain_right_edge - self.ds.domain_left_edge xl = self.ds.domain_left_edge[0] + relx * DW[0] yl = self.ds.domain_left_edge[1] + rely * DW[1] xr = xl + dd * DW[0] yr = yl + dd * DW[1] try: self.lock() w = 256 # pixels data = self.data[field] frb = FixedResolutionBuffer(self.data, (xl, xr, yl, yr), (w, w)) cmi, cma = get_color_bounds( self.data["px"], self.data["py"], self.data["pdx"], self.data["pdy"], data, self.ds.domain_left_edge[0], self.ds.domain_right_edge[0], self.ds.domain_left_edge[1], self.ds.domain_right_edge[1], dd * DW[0] / (64 * 256), dd * DW[0], ) finally: self.unlock() if self.takelog: cmi = np.log10(cmi) cma = np.log10(cma) to_plot = apply_colormap(np.log10(frb[field]), color_bounds=(cmi, cma), cmap_name=cmap) else: to_plot = apply_colormap(frb[field], color_bounds=(cmi, cma), cmap_name=cmap) rv = write_png_to_string(to_plot) return rv
def map(self, field, L, x, y): if ',' in field: field = tuple(field.split(',')) dd = 1.0 / (2.0**(int(L))) relx = int(x) * dd rely = int(y) * dd DW = (self.ds.domain_right_edge - self.ds.domain_left_edge) xl = self.ds.domain_left_edge[0] + relx * DW[0] yl = self.ds.domain_left_edge[1] + rely * DW[1] xr = xl + dd*DW[0] yr = yl + dd*DW[1] try: self.lock() data = self.data[field] frb = FixedResolutionBuffer(self.data, (xl, xr, yl, yr), (256, 256)) cmi, cma = get_color_bounds(self.data['px'], self.data['py'], self.data['pdx'], self.data['pdy'], data, self.ds.domain_left_edge[0], self.ds.domain_right_edge[0], self.ds.domain_left_edge[1], self.ds.domain_right_edge[1], dd*DW[0] / (64*256), dd*DW[0]) finally: self.unlock() if self.takelog: cmi = np.log10(cmi) cma = np.log10(cma) to_plot = apply_colormap(np.log10(frb[field]), color_bounds = (cmi, cma)) else: to_plot = apply_colormap(frb[field], color_bounds = (cmi, cma)) rv = write_png_to_string(to_plot) return rv