python - Matplotlib tripcolor bug? -


i want use tripcolor matplotlib.pyplot view colored contours of of data.

the data extracted xy plane @ z=cst using paraview. directly export data in csv paraview triangulates plane me.

the problem depending on plane position (ie mesh) tripcolor gives me or bad results.

here simple example code , results illustrate it:

code

import matplotlib.pyplot plt import numpy np  p,u,v,w,x,y,z  = np.loadtxt('./bad.csv',delimiter=',',skiprows=1,usecols=(0,1,2,3,4,5,6),unpack=true)  nblevels = 256  plt.figure() plt.gca().set_aspect('equal')  plt.tripcolor(x,y,w,nblevels,cmap=plt.cm.hot_r,edgecolor='black')  cbar = plt.colorbar() cbar.set_label('velocity magnitude',labelpad=10)  plt.show()

results tripcolor

enter image description here

here file causes problem.

i've heard matplotlib's tripcolor buggy, bug or not ?

as highlighted @hooked normal behaviour delaunay triangulation. remove unwanted triangles should provide own triangulation passing explicitly triangles.

this quite easy in case data structured: suggest performing delaunay triangulation in plane (r, theta) passing these triangles initial (x, y) arrays. can make use of the built-in trianalyzer class remove flat triangles (r, theta) triangulation (they might exists due round-off errors).

import matplotlib.pyplot plt import numpy np import matplotlib.tri mtri  p,u,v,w,x,y,z  = np.loadtxt('./bad.csv',delimiter=',',skiprows=1,usecols=(0,1,2,3,4,5,6),unpack=true)  r = np.sqrt(y**2 + x**2) tan = (y / x) aux_tri = mtri.triangulation(r/np.max(r), tan/np.max(tan)) triang = mtri.triangulation(x, y, aux_tri.triangles) triang.set_mask(mtri.trianalyzer(aux_tri).get_flat_tri_mask())  nblevels = 256  plt.figure() plt.gca().set_aspect('equal')  plt.tripcolor(triang, w, nblevels, cmap=plt.cm.jet, edgecolor='black')  cbar = plt.colorbar() cbar.set_label('velocity magnitude',labelpad=10) plt.show()

enter image description here


Comments

Post a Comment

Popular posts from this blog

commonjs - How to write a typescript definition file for a node module that exports a function? -

consider, toml node module can use: // toml.d.ts declare module toml { export function parse(value:string):any; } declare module "toml" { export = toml; } then: /// <reference path="../../../../../defs/toml/toml.d.ts"/> import toml = require('toml'); toml.parse(...); however, node module export single function, such 'glob' ( https://github.com/isaacs/node-glob ). the node usage of module be: var glob = require("glob") glob("*.js", {}, callback(err, files) { ... }); you'd naively expect this: // glob.d.ts declare function globs(paths:string, options:any, callback:{(err:any, files:string[]):void; ...but because typescripts 'import' semantics bit strange, seems can use 'import .. = require()' statement alias modules . trying call: /// <reference path="../../../../../defs/glob/glob.d.ts"/> import blog = require('glob'); results in: error ts2072: mod
Read more

openid - Okta: Failed to get authorization code through API call -

i'm integrating okta own idp server using okta's api. i'm implementing authorization code flow following steps below: in own server, use /api/v1/authn endpoint sessiontoken. use sessiontoken obtain authorization calling endpoint: /oauth2/v1/authorize?client_id=" + clientid + "&sessiontoken=" + sessiontoken + "&response_type=code&response_mode=query&scope=openid&redirect_uri=" + redirecturl + "&state=evanyang&nonce=" it's supposed return response status code 302 , location header containing redirect url code value. however, keep getting response status code 200 , without location header, html body saying "you using unsupported browser." , "javascript disabled on browser." according api documentation: http://developer.okta.com/docs/api/resources/oidc.html#authentication-request , sessiontoken parameter sufficient this: an okta one-time sessiontoken. allows api-based
Read more

thorough guide for profiling racket code -

did googling ("racket profiling", "racket measure performance"), didn't find , there no examples in docs. did google "profile" search on htdp - no luck. (profile (f ...)) output isn't obvious other small snippets. ideally want python's python -m cprofile usage examples. when search "python profiling", both duckduckgo , google give top result: 26.4. python profilers . (although scanned quickly, seems more of reference "thorough usage guide" examples. if had else in mind, perhaps please link that?) the equivalent racket documentation be: profile: statistical profiler . a usage example: #lang racket (require profile) (profile-thunk (thunk (function-to-profile arg0 arg1) )) here (thunk e) convenience (lambda () e) . bigger example: #lang racket (module mod racket (provide f) (define (f) (for/list ([i 10000]) i))) (require (prefix-in mod: 'mod)) (define (f) (for ([i 10000
Read more