Weird Function TRIANGULATE in IDL

Time to update my blog again. It’s been more than a year since my last update after all. Today I encountered with a weird error in IDL when trying to improve computation precision of own polar transform image routines. The source was traced back to function TRIANGULATE. So far I really have no idea about why it occurs. The example is as follows:

IDL> theta=dindgen(360L)
IDL> rad=dindgen(500L)
IDL> xpol=rad#cos(theta/1.8d2*!dpi)
IDL> ypol=rad#sin(theta/1.8d2*!dpi)
IDL> triangulate,xpol,ypol,tri
% TRIANGULATE: Points are co-linear, no solution.
% Execution halted at: $MAIN$

However, if I type the following TRIANGULATE works flawlessly:

IDL> xpol=rad#cos(float(theta/1.8d2*!dpi))
IDL> ypol=rad#sin(float(theta/1.8d2*!dpi))
IDL> triangulate,xpol,ypol,tri
IDL> help,tri
TRI             LONG      = Array[3, 358920]

This is really annoying, as I can’t figure out how come co-linearity can be related to float or double precision whatsoever. Not sure if anyone could help this for me, but I just throw this problem here and see if there’s anyone else faced with the same encounter so that we can potentially sigh together…

Weird Function TRIANGULATE in IDL

My Record of Installing IDL, SSW and SSWDB on My New Mac OS X

I got a new laptop Mac OS X from Dave and needed to install IDL, SSW and SSWDB on it. The following are the main steps how I installed it successfully.

 

Basic Steps of installation of IDL 8.3 on my Mac OS X

1. Download the installation package from the EXELIS.

2. Click on the icon and follow the instructions of the installation wizard.

3. Put the license into the following directory: Applications/exelis/license.

4. The license wizard pops out automatically after the installation completes. Now select the license file in the panel.

5. Now IDL should be successfully launched by double clicking on the workbench icon, however, if you type idl in the terminal, you will be informed that idl: Command not found. Environment has to be set.

6. Type the following commands in the terminal:

vi .tcshrc

7. See if there is the following lines within the file, insert them otherwise.

setenv IDL_DIR /Applications/exelis/idl83
setenv OS darwin.

In Tshell command line, typing idl will now successfully run IDL.

 

Steps of installation of SSW

1. Fill out and submit the form at http://www.mssl.ucl.ac.uk/surf/sswdoc/solarsoft/ssw_install.html.

2. Follow the instructions in that website and the links therein.

3. Type sswidl in the terminal, but fails, so insert the following lines into the .tcshrc file:

set path=(/usr/local/bin /sw/bin /sw/sbin /usr/local /usr/local/ssw/gen/mirror /usr/X11R6/bin /Applications/exelis/idl83/bin $path .)
setenv SSW_INSTR "lasco secchi"
setenv SSW ssw
setenv SSWDB sswdb
setenv ssw_quiet 1
source $SSW/gen/setup/setup.ssw
source $SSW/soho/lasco/setup/setup.lasco_envsetenv SSWDB /usr/local/ssw/sswdb

Installation of SSWDB

1. Anyway I am not able to generate a configuration file from http://www.lmsal.com/solarsoft/sswdb_configure.html, but there is a file called setup.sswdb_upgrade_template in directory of ssw/site/setup. Add lines indicating the packages in need, or, if there are packages already listed in the template, remove the corresponding “#” at the beginning. Root is required:

sudo vi setup.sswdb_upgrade

2. Follow the steps listed at http://www.mssl.ucl.ac.uk/surf/sswdoc/solarsoft/sswdb_install.html.

3. If Perl is installed in a non-default directory, type
ln -s /usr/bin/perl /usr/local/bin/perl
in the terminal.

4.IDL> sswdb_upgrade, /spawn, /passive_ftp

5. Be patient about Step 4, and it took me rather long.

6. Done!

My Record of Installing IDL, SSW and SSWDB on My New Mac OS X

詩數首

春日葳蕤百越東,棹歌聲裏木棉紅。

回看無數詩和畫,盡在一谿煙雨中。

*********************************************************************************

風斜雲淡月幽冥,耿耿銀河點點星。

夜色惹人眠不得,憑軒醉臥數流螢。

*********************************************************************************

學書廿年有感

佳致何曾經歲失?長臨南帖北碑文。

及今池水還無染,更以孜孜學右軍。

*********************************************************************************

夜來弄墨足風騷,神入中山紫兎毫。

沉醉不知天欲曙,忽驚窗外啓明高。

*********************************************************************************

玉魄西欹藏海波,星槎飛急渡天河。

夜深露重行人絕,織女牽牛相盪摩。

*********************************************************************************

愁思縈懷細雨天,北邙一帶草芊芊。

想來往事還如昨,物換星移又十年。

詩數首

小白,走好!

  我家的小白,走了……
  自元宵走出後,已是三日未歸。直到昨天下午,阿公想去蕉林看看,結果卻在蕉林的草叢中發現了小白的遺體。雙眼緊閉,已是渾身冰涼,四肢僵硬了。看來已是亡去多時了。
  家裏養小白應該是有四五年了吧。它雖說叫小白,但頭卻是黑色的,背部也有幾塊黑色的斑點。眼睛固然小小的,但卻是炯炯有神。當初還是舅父不知從何處帶回來的,才幾個月大。這隻小狗不同以往養過的,不用教太多便非常聽話了。我能想起它唯一做過的一件壞事是把我晾在階前的波鞋的鞋帶給咬斷了,真是個淘氣鬼!
  夜晚它是非常警醒的,每每聽到異響,便以洪亮而清澈的聲響吠之。還記得有次晚飯後,見它伏在臺階上,表情甚是可愛,於是輕輕撫摸着它的頭。想必它正在熟睡中,被我的舉動嚇壞了,倏地下躥了起來逃走,驚恐地吠了幾聲,回頭見原來是我,便又搖着尾巴匆匆跑來舐我的腳。
  每次回家,無論是深夜,還是日頭,只要聽到響聲,它總是第一時間衝向我奔來。我也總是摸摸它烏黑的腦袋予以回報,這樣它的尾巴搖得更快了,似乎都要把屁股摔出去了。它在白天喜歡曬太陽,若是見我來逗它,它便四肢朝天地示意叫我幫它撓撓它的肚子。看着它的樣子,很是享受呢!
  可是到如今,一切戛然皆是夢了,徒留回憶矣,滿腹心事,又可與胡訴說?唯聽窗外瀟瀟冷雨。我所能做的,便是願它在天堂安好,那裏或無嚴寒酷暑,一切都美好於人間。來世有緣者,我們再相見吧!我會十分想念你的,小白!

小白,走好!

雜記

  昨日午至張老師居所,契闊談宴,由書法而至做人。下午隨諸小友一同上課。時光荏苒,真如白駒過隙,自四歲半從師學書以來,累十二年餘,風雨無阻,無奈後逮高三學業繁雜而綴,後更因背鄉之省外讀書而不得跟隨。如今再來學堂,竟又是一個七年。恩師堂上所言,仍是與記憶之中無二。然而物換星移幾度秋,學生早已不同當年。曩時之同窗,至今恐已是天南地北,不得再見。留得通信方式者更無幾人,甚爲可惜!
  今日略有雅興,練字不覺之中逾半日。午前稍動筆墨,覺意猶未盡。食晏後,再執筆椽,一抒翰墨之情。所恨廢紙無數,竟無一滿意者。或衰於字,或失於章法。數添墨,又竭。不覺早已華燈盡上,月昃西阿。停筆始覺肚飢不已,只得到此爲止,他日再續。既結緣翰墨,則爲一生翰墨人。

雜記

伍瑞隆詩選

  夜讀書,見嶺南先人伍瑞隆題詩,字體法度精嚴之中,瀟灑秀逸,遒麗豐腴,有晉二王之風。故特尋其生平事蹟,偶見他處錄其詩若干首,以爲妙,咸集之於斯。

  伍瑞隆,字國開,一字鐵山。明末香山人。因結廬於鳩、艾二山之間,故自號鳩艾山人。

自題墨牡丹
睡起簾前曉日紅,名花開遍洛城東。
幽人碾墨傳清艷,只在踈踈幾筆中。
按:尾聯倒數第二字因詩軸年久似損壞而字跡不清,存疑。

初秋送弟有開兒昌年還玉溪兼柬劉衍崖舅氏
老舅亦已老,去年七十時。
賞花棋百局,留客酒千卮。
一別冬還夏,雙魚遠複遲。
天涯送兒弟,相憶淚如絲。

贈劉良倩
詞場俠少年,誰在劉生前。
傲骨世不賞,清才吾獨憐。
松為稽子韻,蘭是楚人編。
終令乘風起,名高日月邊。

祝棲老和尚七十有一
我歲乙酉夏,師年丙戌春。
所爭九個月,同作七旬人。
老宦曾何補,名僧自有真。
鼎湖山月白,為照劫灰塵。

依韻答胡明府絹庵
重陽風雨不成秋,官署山城一樣愁。
載酒未辭籬畔菊,題詩空上水邊樓。
邀歡豈必期來歲,卜醉還應續舊遊。
寄語江東張季子,好留清興在三洲。

答程周量
香風吹動筆花新,羨爾當今第一人。
慧鏡自天懸似日,文泉隨地出如神。
遭逢偶處還非偶,鑒別真時卻未真。
世事但存知己在,未須懷寶泣孤臣。

荔園座上贈諸子
楊柳新寒水面風,吹花送酒畫筵紅。
芙蓉洗出青天上,薜荔晴分碧霞中。
四海莫悲坭畔絮,萬家愁聽澤邊鴻。
詩成寫向相思樹,聊與君期紫極宮。

客岡州送棲壑和尚還鼎湖
昔人從師斷二愚,別來原只在須臾。
肉身喜見生菩薩,法座前稱老腐儒。
滿月一簾天在水,香風萬里露中珠。
憑師試展琉璃臂,入我岡州立雪圖。

即 事
(一)
黃梁一陣曉來風,龍眼吹殘橄欖空。
閒步長洲訪春色,西河無數落花紅。

(二)
河塘簇簇白藤花,岡尾岡頭十萬家。
見話橫江人不渡,峭帆風緊桂洲斜。

(三)
西頭南尾盡煙波,綠玉城中幽恨多。
林步酒家休得醉,寨前芳草奈愁何。

(四)
沙塘秋好發紅蓮,倚劍峰頭盡日眠。
下巷騎牛歸上巷,水邊沽酒過山邊。

(五)
秋風高隴不堪愁,雙美橋中廣渡頭。
鳳嶺駝山都莫問,更誰六社問瀛洲。

(六)
石橋坑下水流澌,烈女祠前花萬枝。
借問鳳凰山下月,至今猶照白蓮池。

(七)
何人移棹過江門,白石黃連尚有村。
閒步海洲懷古處,道娘無墓可招魂。

鳩艾山賦
(一)
雨痕初過綠生坭,睡起西窗意欲迷。
—樹紫微風盡落,亂紅堆地午鶯啼。

(二)
各花爭報各城秋,各色花光盡入樓。
感謝西風遠來意,早傳新爽到心頭。

(三)
誰教蝴蝶上花枝,花自多情學蝶飛。
我向北窗眠一覺,與花同作蝶飛歸。

(四)
紅蓼花開紅蓼沙,偶然移種入窗紗。
秋風未到花先到,想到秋來花又花。

(五)
獨岸亭亭秋海棠,秋來猶未見紅妝。
藕泉泉水芳如酒,為灌花魂作意香。

(六)
憑誰招得百花魂,毒霧連天草木昏。
獨有摩霄幾群鶴,月明清唳到山門。

(七)
一聲一聲飛上天,天邊佳氣碧如煙。
憑君說盡人間事,有客勞歌夜未眠。

(八)
山城無力養春風,春盡群花在夢中。
好逐鶴聲齊喚起,滿城香豔入新叢。

(九)
小小金魚五色鱗,往來無定為憐春。
池邊倒有千花影,魚亦隨花巧伴人。

(十)
山後山前五色禽,飛鳴處處是同心。
並頭花上連枝樹,一隊飛來一隊尋。

(十一)
春草無端五色花,佳人醉處日將斜。
殷勤撲上雙金鳳,一路憑君帶到家。

(十二)
嶺外寒輕節每遲,茲花九月早盈枝。
思君便得君來賞,一盞新醪一首詩。

(十三)
一自移居雩樂園,種花惟有菊花存。
西風不赴登高節,留待君來照綠尊。

(十四)
籬邊花草亦知愁,十月才開九月秋。
莫待深言此中事,餐英聊上白蘋洲。

開元寺試泉
片雲依絕磵,一水眄秋陽。
寄語神仙客,梅花遜不香。

竹枝詞
(一)
蝴蝶花開蝴蝶飛,鷓鴣草長鷓鴣啼。
庭前獨有相思樹,落盡相思人未歸。

(二)
妾愁那似比君愁,妾欲留君君不留。
君去桃花潭水上,有人來往送蘭舟。
(三)
楊柳千條與萬條,繫君蘭舸莫乘潮。
君心那似楊花亂,才逐風飄雨又飄。

(四)
紅霞映水月如煙,堤上行人駐馬鞭。
我唱鷓鴣君不信,君看江水一枝蓮。

(五)
閒情日日踏堤行,綠水紅霞繞樹明。
何處愁心最關切?夕陽鶯語一聲聲。

伍瑞隆詩選