IDL之伪彩色线的KML文件

15195775117

一台环境监测车走在路上，根据数据可以画出一条道路的污染轻重。
我想在谷歌地球上伪彩色线，但是，
用plot彩线存的kmz文件，在谷歌地球中是不显示颜色的，
而且点之间的连线都会呈现。
plot不行，但是image和contour的平面图却可行，
把没经过的位置以NaN填充，在GE上就会呈现镂空效果。
这样，我们就可以随意设置伪彩色的位置与形状了，不再必须是方方正正的了。
案例：
IDL> a=findgen(100)/10+110
IDL> b=findgen(100)/10+30
IDL> c=dist(100)
IDL> c[20:80,30:90]=!values.f_nan
IDL>   CONTOUR(c,a,b,N_LEVELS=10, RGB_TABLE=13, /OVERPLOT,/fill, GRID_UNITS=2, TRANSPARENCY=50)
这个用image(c,a,b)也行
以kmz格式存储并用谷歌地球打开，效果如下：

以上，已经实现了KML伪彩色线的技术可行性，
以下楼层我争取画一个实践性的图，
因为contour有插值，对原始数据有改动，
所以我选择用image实现，
待续......

15195775117

用contour做镂空图的缺陷
图中黑圆是样本位置，彩色是等值线插值出来的，
可见，哪里会填充颜色，跟插值计算有关，
所以，想让伪彩色只填充真实样本位置周围，最好使用image(value,lon,lat)这样形式

源代码：
a=findgen(100)/10+110
b=findgen(100)/10+30
c=findgen(100)
c[50:60]=!values.f_nan
fig= CONTOUR(c,a,b,/fill,rgb_table=33)
fig=plot(a,b,/overplot,sym=24,linestyle=6)