python - Some values of matrix do not appear in the plot by Matplotlib -


i created empty reference matrix csv, located (x,y) position on matrix (and printed them out), , designated 100 position on matrix. each x value in ref_mass pandas series.

ref_df = pd.read_csv(ref_file) reference = np.zeros(shape=(1201,len(ref_df))) ref_mass = ref_df["mass"]  i, mass in enumerate(ref_mass):   print ref_mass[i].astype(int) - 300, # print (x,y)   reference[(ref_mass[i].astype(int) - 300),i] = 100

every (x,y) printed out correctly. however, there no value in plot of (x,y). what's wrong here? checked reference matrix, has 100 in every column rightly.

the (x,y): 

547 0 265 1 124 2 39 3 509 4 # shown 240 5 # shown 105 6 24 7 355 8 137 9 28 10 # shown 394 11 163 12 48 13 347 14 132 15 # shown 24 16

the plot: enter image description here

plot code:

if __name__ == '__main__':   mpl_toolkits.mplot3d import axes3d   import matplotlib.pyplot plt   import matplotlib   matplotlib.matplotlib_fname()    plt.ylabel('m/z')   plt.xlabel('peptides')    plt.imshow(reference, aspect='auto', cmap='terrain')   plt.colorbar()   plt.tight_layout()    plt.show()

every pixel in final image represents 3 or more data points. renderer has decide color out of 2 times blue, 1 time white map pixel. statistically, blue twice white, such 66% of data points not shown.

the number of 3 pixels comes rough calculation: image has 480 pixels (which can either find out in picture program or calculating figuresize*dpi). have 1200 datapoints (seen axes). have margin of ~10% @ both ends; have 1200/(0.8*480) = 3.1 datapoints per pixel in final image.

imshow interpolation

you can use interpolation on image make pixels appear, e.g.

plt.imshow(..., interpolation="sinc")

the result may not appealing visually.

change resolution

you can make sure final plot comprises 1 pixel per datapoint. i.e. 1200 pixels dpi of 100 can do

m = 0.1 plt.figure(figsize=(8, (1+2.*m)*1200./dpi )) plt.subplots_adjust(bottom=m, top=1.-m) plt.imshow(...)

filter data

another option change data, such 1 pixel becomes 3 pixels along y direction.

import matplotlib.pyplot plt import numpy np; np.random.seed(1) import scipy.ndimage.filters filters  = np.zeros((1200, 16)) in range(16):     a[i*70+21, i] = 1  kernel = np.array([[0.,1.,0.],                    [0.,1.,0.],                    [0.,1.,0.]]) anew = filters.convolve(a,kernel,mode='constant', cval=0)  im = plt.imshow(anew, aspect="auto") plt.colorbar(im)  plt.show()

enter image description here

drawing lines

import matplotlib.pyplot plt import numpy np  = np.zeros((1200, 16))  im = plt.imshow(a, aspect="auto", vmin=0, vmax=1) plt.colorbar(im)  in range(16):     plt.plot([i-.5, i+.5], [i*70+21,i*70+21], color=im.cmap(1.))  plt.show()

Comments

Post a Comment

Popular posts from this blog

Is there a better way to structure post methods in Class Based Views -

i'm using class based views in django , never specify url in ajax calls , omit action parameter in form because know these requests go through post method of class based view. notice view - profile page make many requests depending on user updating on page - post method become convoluted , if/else mess example: def post(self, request, *args, **kwargs): # handle user changing available date if request.post.get('availabledate') != none: self.updateavailabledate(request) return jsonresponse({'result':'success'}) if request.post.get('newprojectname') != none: creative_user = creativeuserprofile.objects.get(id = request.user.id) project = project.create_project(creative_user, request.post['newprojectname']) project.save() return jsonresponse({'projectid': project.id }) if request.post.get('projectimage') != none , request.files['file']: proj...
Read more

performance - Why is XCHG reg, reg a 3 micro-op instruction on modern Intel architectures? -

i'm doing micro-optimization on performance critical part of code , came across sequence of instructions (in at&t syntax): add %rax, %rbx mov %rdx, %rax mov %rbx, %rdx i thought had use case xchg allow me shave instruction , write: add %rbx, %rax xchg %rax, %rdx however, dimay found agner fog's instruction tables , xchg 3 micro-op instruction 2 cycle latency on sandy bridge, ivy bridge, broadwell, haswell , skylake. 3 whole micro-ops , 2 cycles of latency! 3 micro-ops throws off 4-1-1-1 cadence , 2 cycle latency makes worse original in best case since last 2 instructions in original might execute in parallel. now... cpu might breaking instruction micro-ops equivalent to: mov %rax, %tmp mov %rdx, %rax mov %tmp, %rdx where tmp anonymous internal register , suppose last 2 micro-ops run in parallel latency 2 cycles. given register renaming occurs on these micro-architectures, though, doesn't make sense me done way. why wouldn't register renamer ...
Read more

c# - Asp.net web api : redirect unauthorized requst to forbidden page -

im trying redirect unauthorized request forbidden page instead i'm getting forbidden page in response body , how can fix ? here's startup class : app.createperowincontext(storecontext.create); app.createperowincontext<applicationusermanager>(applicationusermanager.create); app.createperowincontext<applicationsigninmanager>(applicationsigninmanager.create); app.usecookieauthentication(new cookieauthenticationoptions { authenticationtype = defaultauthenticationtypes.applicationcookie, expiretimespan = timespan.fromdays(30), }); app.useexternalsignincookie(defaultauthenticationtypes.externalcookie); app.usetwofactorsignincookie(defaultauthenticationtypes.twofactorcookie, timespan.fromminutes(5)); app.usetwofactorrememberbrowsercookie(defaultauthenticationtypes.twofactorrememberbrowsercookie); the method im trying reach : [httpget] [authorize(roles = "admin")] public string getcurrentusername() { return userman...
Read more