Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  913 1682 (N35E19 L=306)   377    132        0 ( 846/1651- 957/1724) -15.3 E25-E14 N32-N36
      2  805 1521 (N22E21 L=303)   665    188        0 ( 759/1436- 836/1599) -48.0 E26-E17 N16-N26
      3 1429 1311 (N17W21 L=346)   783    203        0 (1341/1245-1493/1385) -26.6 W14-W26 N13-N20
      4 1831 1281 (N23W56 L= 21)  1151     89        4 (1797/1234-1868/1316) -45.0 W52-W61 N20-N26
      5  142 1045 (S13E68 L=257)  1207    170        3 ( 102/ 981- 175/1085)  36.9 E75-E63 S16-S10
      6  224 1037 (S13E57 L=268)   403     58        0 ( 203/1004- 248/1067)  36.9 E59-E54 S14-S11
      7  315  958 (S17E48 L=277)   321     76        0 ( 283/ 940- 360/ 973)   8.1 E50-E43 S18-S16
      8  985  700 (S26W02 L=327)   748    287        0 ( 878/ 614-1040/ 815)  66.8 E03-W04 S32-S18
      9 1380  615 (S25W30 L=355)   298     88        0 (1337/ 590-1415/ 646)  45.0 W27-W31 S26-S22
     10 1486  424 (S35W47 L= 12)   612     96        0 (1439/ 392-1535/ 459)  56.3 W44-W50 S37-S31
     11  504  324 (S60E49 L=276)   713     55        2 ( 476/ 300- 537/ 347)  12.5 E59-E41 S62-S58
  

    OPENR, u, binfile, /COMP, /GET_LUN
    head = FLTARR(6) ;[xsize,ysize,radius,hh,mm,ss]
    READU, u, head
    data = INTARR(FIX(head[0]),FIX(head[1]))
    READU, u, data
    FREE_LUN, u
    

• each filament has the value of Fil_ID
• all values are in pixel

  c = WHERE(data EQ Fil_ID, area)
  f = '(2(A,I6))'
  PRINT, FORM=f, "Area of Fil_ID", Fil_ID," =", area 

  data = ROTATE(data,7)
  Solar_y = 2048-Solar_y