python - Interpolation of large 2d masked array -

-

i have numpy masked matrix. , wanted interpolation in masked regions. tried rectbivariatespline didn't recognize masked regions masked , used points interpolate. tried bisplrep after creating x,y,z 1d vectors. each of length 45900. took lot of time calculate bsplines. , gave segmentation fault while running bisplev . 2d matrix of size 270x170.

is there way make rectbivariatespline not include masked regions in interpolation? or there other method? bisplrep slow.

thanking you, indiajoe

update : when grid small scipy.interpolate.rbf 'linear' function doing reasonable job. gives error when array large.

is there other function allow me interpolate , smooth matrix?

i have concluded following. correct me if wrong.

1) rectbivariatespline requires perfect filled matrix , hence masked matrices cannot used.

very late, but...

i have problem similar yours, , getting segmentation fault bisplines, , memory error rbf (in "thin_plate" function works great me.

since data unstructured created in structured manner, use downsampling half or 1 third of density of data points, can use rbf. advise (very inefficient, still better not doing @ all) subdivide matrix in many overlapping regions, create rbf interpolators each region, when interpolate 1 point choose appropriate interpolator.

also, if have masked array, still perform interpolation in unmasked array, apply mask on result. (well no, see comments)

hope helps somebody


Comments

Post a Comment

Popular posts from this blog

jasper reports - Fixed header in Excel using JasperReports -

-
how can create excel report fixed header using japserreports? mean need header fixed when scroll excel file. this possible adding couple properties in jrxml file each report. take @ advanced excel features freeze panes. if wanted freeze after first column header (down left side basically) this: <statictext> <reportelement style="sans_bold" mode="opaque" x="0" y="60" width="104" height="20" forecolor="#ffffff" backcolor="#666666"> <property name="net.sf.jasperreports.export.xls.auto.filter" value="start"/> <property name="net.sf.jasperreports.export.xls.column.width" value="110"/> <property name="net.sf.jasperreports.export.xls.freeze.column.edge" value="left"/> </reportelement> &
Read more

media player - Android: mediaplayer went away with unhandled events -

-
i need duration of audio file series of voice announcements need play app. have added audio files resources , play fine. sample code below works perfect intended purpose: return duration of audio files. here code: float getdurationofaudioresource(locationenum loc, context context){ float duration = 0; try { mediaplayer mp; mp = mediaplayer.create(context, getaudioresource(loc)); duration = mp.getduration(); mp.release(); mp = null; } catch (illegalstateexception e) {e.printstacktrace(); logerror(25, "testdescitem:fault::could not open mediaplayer object audio resource.");} return duration; } here's weird thing. code called in main activity prepares set of audio instructions given test. there no errors within activity. second activity called, long string of errors on logcat. 03-07 13:23:43.820: i/actionlogger(21435): gentest_info_test #0 created. 03-07 13:23:43.830: i/actionlogger(21435): gente
Read more

python - ('The SQL contains 0 parameter markers, but 50 parameters were supplied', 'HY000') or TypeError: 'tuple' object is not callable -

-
import pyodbc,nltk,array cnxn = pyodbc.connect('driver={mysql odbc 5.1 driver};server=127.0.0.1;port=3306;database=information_schema;user=root; password=1234;option=3;') cursor = cnxn.cursor() cursor.execute("use collegedatabase ;") cursor.execute("select * sampledata ; ") cnxn.commit() s=[] j=[] x=[] words = [] w = [] sfq = [] pos=[] entry in cursor: s.append(entry.injury_type),j.append(entry.injury_desc) nltk.tokenize import punktwordtokenizer nltk.corpus import stopwords tokenizer = punktwordtokenizer() english_stops = set(stopwords.words('english')) in range(0,26): # filter stop words tokenizer.tokenize(j[i]) w.append([word word in tokenizer.tokenize(j[i]) if word not in english_stops]) in range(0 , 26):#converting tokenzied text ito string sfq.append(" ".join(w[a])) replacers import regexpreplacer replacer = regexpreplacer() in range (0,26):#pos tagging replacer.repl
Read more