IDL之批量写KML脚本

15195775117

                               
登录/注册后可看大图

附件有IDL代码和KML文件，
思路是用IDL批量写KML脚本
head.txt和tail.txt是kml固定的头尾，up和down是Polygon的前后内容，中间不断填入的是坐标信息，批量写入即可。
本来想把椎体布满全球的，但是发现写的时间比较长，于是中断了程序，自己把尾巴补上。
目前写的效率不高，以后再优化。

15195775117

以上案例之所以没有铺满全球，是因为IDL写KML脚本的速度比较慢，这倒不是IDL语言的问题，
我在IDL语言中反复打开、续写、关闭文件，这可能是导致速度慢的原因，

以前我也做过反复开关、读写文件的操作，非常慢，
由于以前没有衡量过开关文件的耗时状况，
经过两次体验，我发觉开关文件是个费时的操作，以后避免反复开关文件。
以下我试着一次读写文件，应该可以节省不少时间。
待续......

15195775117

今天把算法优化了一下，瞬间把“钻石”布满全球：
% Time elapsed: 0.27900004 seconds.
果然，反复开关文件是罪魁祸首！
附件中有IDL脚本和生成的kml文件。
pro rw
  tic
  head='<?xml version="1.0" encoding="UTF-8"?>'+$
    '<kml xmlns="http://earth.google.com/kml/2.1"><Document>'
  tail='</Document></kml>'
;地标的前半段：
  up='<Placemark><name>p1</name><Style><LineStyle>'+$
    '<width>5</width><color>ff0000ff</color>'+$
    '</LineStyle><PolyStyle><outline>1</outline>'+$
    '<color>7f00ff00</color></PolyStyle></Style>'+$
    '<Polygon><extrude>0</extrude>'+$
    '<altitudeMode>relativeToGround</altitudeMode>'+$
    '<outerBoundaryIs><LinearRing><coordinates>'
;地标的后半段：
  down='</coordinates></LinearRing>'+$
    '</outerBoundaryIs></Polygon></Placemark>'
  ;写入kml固定开头
  openw,lun,'C:\Users\Administrator\Desktop\out.kml',/get_lun,/append
  writeu,lun,HEAD
  for i=-180,180,4 do begin
    for j=-90,90,4 do begin
      A=STRING(10B)+strtrim(string(i),2)+','+strtrim(string(j),2)+',300000'
      B=STRING(10B)+strtrim(string(i+2),2)+','+strtrim(string(j),2)+',300000'
      C=STRING(10B)+strtrim(string(i+2),2)+','+strtrim(string(j+2),2)+',300000'
      D=STRING(10B)+strtrim(string(i),2)+','+strtrim(string(j+2),2)+',300000'
      E=STRING(10B)+strtrim(string(i+1),2)+','+strtrim(string(j+1),2)+',0'
      F=STRING(10B)+strtrim(string(i+1),2)+','+strtrim(string(j+1),2)+',600000'
      writeu,lun,up,A,B,E,A,STRING(10B),down,$
        up,A,B,F,A,STRING(10B),down,$
        up,A,B,E,A,STRING(10B),down,$
        up,B,C,E,B,STRING(10B),down,$
        up,B,C,F,B,STRING(10B),down,$
        up,C,D,E,C,STRING(10B),down,$
        up,C,D,F,C,STRING(10B),down,$
        up,D,A,E,D,STRING(10B),down,$
        up,D,A,F,D,STRING(10B),down
    endfor
  endfor
  ;写入收尾：
  writeu,lun,tail
  free_lun,lun
  toc
end

15195775117

边框有点过粗，设置细一点：

                               
登录/注册后可看大图