最近工作需要用到最遠(yuǎn)點采樣fps算法, 找了一些blog學(xué)習(xí)學(xué)習(xí), 寫了個原型: 測試一下, 跑的沒問題, 但是速度感人 所有點的歐氏距離計算做個向量化優(yōu)化: 這一下子快了1...
發(fā)簡信
IP屬地：浙江
-
FarthestPointSampling(FPS)算法學(xué)習(xí)&實踐
-
記一個gcc7代編譯器下allocate-source中有關(guān)array bounds的離譜bug
Bug情況 最近寫Fortran代碼的時候遇到了個怪bug, 代碼如下: 然后離譜的來了, 我發(fā)現(xiàn)在gcc750和gcc710里面跑出來的shape_data的lbound都...
-
cmake - 編寫單元測試在cmake中構(gòu)建測試?yán)褪莈nable_testing()+add_test()這一套。本篇文章是我在學(xué)習(xí)cmake構(gòu)建單元測試的一個總結(jié)筆記块蚌。 工程概覽 如下是我所定義的...
-
WPS ungrib編譯報錯 undefined reference to 'jpc_decode'
最近在一個新機子上用Spack配置好編譯器, hdf5, netcdf-c, netcdf-fortran, jasper等一系列依賴庫后編譯了最新版本的WRF-4.4.2和...
-
地理坐標(biāo)轉(zhuǎn)換-proj4
目標(biāo)：WGS-84 轉(zhuǎn)為 CGCS2000坐標(biāo) 使用proj4 測試數(shù)據(jù) 舉例：天安門在不同類型坐標(biāo)系下的經(jīng)緯度值 WGS-84經(jīng)緯度：116.391349,39.90737...
-
但是這不是netCDF4嗎??
利用xarray寫netcdf文件
三維模式輸出的數(shù)據(jù)很大机蔗，傳輸起來很麻煩蒲祈。有時候我們只需要幾個關(guān)鍵變量，因此把他們抽取出來寫成nc文件比較方便萝嘁。 例如梆掸，要單獨保存megan輸出的isoprene變量isop為...
-
EBAS觀測數(shù)據(jù)處理
記錄一下 EBAS-Thredds[https://ebas.nilu.no/thredds/]數(shù)據(jù)的處理流程, 以全站點的ozone為例 基于前一個EBAS觀測數(shù)據(jù)下載[h...
-
m3diff對比分析ioapi-netcdf文件
最近要在一個新機器上跑CMAQ, 先對機器和模式做個校驗, 以前都是跑出來自己寫個ncl腳本比對的, 現(xiàn)在覺得用ioapi的m3diff也挺好, 不過還需要基于那個結(jié)果進行后...
-
EBAS觀測數(shù)據(jù)下載
記錄一下 下載EBAS-Thredds[https://ebas.nilu.no/thredds/]數(shù)據(jù)的流程, 以全站點的ozone為例 1. 下載ebas-thredds...
-
CFSR 6小時數(shù)據(jù)批量下載腳本-WRF模擬用
WRF官方[https://www2.mmm.ucar.edu/wrf/users/download/free_data.html]提供的CFSR數(shù)據(jù)下載地址是在UCAR[h...
-
FORTRAN調(diào)用C函數(shù)的參數(shù)不一致問題-自動額外傳參特性
最近在debug的時候看到CMAQ ioapi的code, envgets.c中和ENVSTR函數(shù)和m3utilio.F中的接口參數(shù)不一致, 但是可以正常跑, 之前沒注意過,...
-
CS61C-2022-su20作業(yè)記錄
記錄一下CS61C 20su的各個lab和proj的心得, 代碼放在了這里[https://github.com/Roadelse/CS61C-2022-su20] proj...
-
ImageMagick日常使用格式記錄
記錄一些日常用ImageMagick處理圖像的格式 1. 剪切圖的一部分 magick convert src.png -gravity East -crop 200x300...
-
觀測數(shù)據(jù)處理-國控環(huán)境監(jiān)測站點
工作需要用到我國環(huán)境監(jiān)測站點的污染物濃度數(shù)據(jù)(感謝大佬的分享), 數(shù)據(jù)很全, 不過csv格式在分析的時候尤其大量數(shù)據(jù)分析的時候并不友好, 所以一般要二次處理一下 上學(xué)的時候搞...
-
powershell數(shù)組賦值遇到的bug, 開頭的逗號不是一直都能用的
今天寫Powershell的時候遇到個bug, 搜了半天沒搜到(不過最后解決后感覺以前似乎看到過, 可能我搜的內(nèi)容不對) 簡單地說, powershell弄個數(shù)組一般就直接逗...
-
numpy slice broadcast的坑, 和ncl或者xarray不同
過去用ncl作array slice習(xí)慣了, 不同的dimension之間就是直接乘起來, 例如 即可得到一個2x2的數(shù)組 結(jié)果今天在用numpy處理數(shù)據(jù)的時候發(fā)現(xiàn) 就只能得...
-
Sublime修改Build Python的默認(rèn)指令
網(wǎng)上查到的配置sublime直接運行python腳本的方法大多是新建一個Build python的設(shè)置, 而直接以默認(rèn)的Python build運行的話可能會報錯提示找不到指...
-
好像目前是隱藏了中國站點..但是手動輸入站點編號比如北京是54511, 歷史數(shù)據(jù)還是可以看的
不顧20年的似乎不行
另外歷史數(shù)據(jù)我之前都down了,有需要的話可以發(fā)給你下載wyoming大學(xué)的探空數(shù)據(jù)
水平有限,歡迎指正交流,共同進步! 俄懷明大學(xué)將全球探空站的數(shù)據(jù)共享使用,并通過計算了很多的變量指數(shù).非常有用,近期有這方面的需求,暫時寫了一段未經(jīng)優(yōu)化的代碼下載.如果要大量...
-
mpi4py 快速上手
在上一篇中我們介紹了如何安裝和使用 mpi4py酸钦，下面我們以幾個簡單的例子來展示怎么使用 mpi4py 來進行并行編程，以使讀者能夠快速地上手使用 mpi4py咱枉。這些例子來自...