原创 利用python批量下載ERA的數據

終於開始用python了,只學了兩天的python小白上線下載數據(被逼的)。 其實這個有一些很好的教程,我直接粘貼好了。批量下載的教程 其中的一些參數官網的解釋很清楚的,可以先試着下載一個月的數據,再獲取它的request提供參考 最後

原创 IDL 讀取nc數據代碼

file_ID = NCDF_OPEN(theseFiles[fidx],/NOWRITE) ;open netCDF file for READ only Tag = NCDF_INQUIRE(file_ID) t

原创 IDL製圖配色

之前利用IDL製圖老是沒有辦法把colorbar裏面的顏色和實際數值在圖裏的顏色對應起來,今天好像知道怎麼弄了,先記錄一下。 我分的-4-4區間,分成8級,顏色的讀入也是讀8個,則那麼ncolor就要寫8+1B。

原创 Arcgis自帶世界矢量圖

E:\ArcGIS\Desktop10.0\TemplateData\TemplateData.gdb裏面有Arcgis自帶的world的矢量圖,只要世界大陸這個屬性基本就可以了。

原创 IDL圖形保存圖片時文字變黑塊的臨時解決方法

忙着做畢設沒想到出了這麼個問題,生成的圖被打了碼。。。哭泣在網上搜到了這麼個辦法,轉載至連接測試代碼如下,注意其中的xtext_orientation、ytext_orientation、baseline的用法。座標軸標題使用orient

原创 IDL提高程序運行效率

    前兩天跑了一個最短路徑,在找出需要搜索的點那放了兩個循環,因此大約2分鐘或者一分半出一張361*361的圖,其實搜索的範圍就這個圖的一半,我當時覺得已經算是在降低運行時間了,但是由於跑的實在太慢,而我畢設的數據量又很大,因此再進行

原创 IDL輸出txt換行的問題

直接用printf輸出的txt每行的內容有限制,所以百度了一下設置openw函數的width關鍵字,該關鍵字決定了打印的寬度(column),默認情況下width=80,也就是打印80個字符後換行。如果你要每行輸出六十個數,那麼width

原创 關於處理長時間序列遙感數據的小貼士——不定時更新

0.同期數據替換的時候,一般會先從文檔找可信度比較高的那一期數據。之後利用另一期數據填補選擇的那期數據的缺失值,此時一定要對比兩個數據的差= =然後再考慮一下1.空間濾波有時候會擴大之前數據異常值的的範圍,特別是均值濾波,要小心。

原创 IDL制legend小整理

今天做了超級多的圖,順手粘一下legend的code,後面再整理 labels = ['< 0','0-50', '50-100', '100-150', '150-200', '200-250', $ '250-300', '