Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1285 1513 (N32W02 L=124)   346    104        0 (1253/1461-1312/1551) -45.0 E01-W04 N28-N33
      2 1700 1410 (N39W38 L=160)   455     69        2 (1666/1386-1734/1434)  15.9 W34-W42 N37-N40
      3 1626 1276 (N27W31 L=153)  1383    233        2 (1551/1179-1697/1337) -19.7 W24-W38 N23-N30
      4 1487 1276 (N23W21 L=143)   389    122        1 (1453/1222-1508/1330)  76.0 W19-W21 N18-N26
      5  556 1289 (S01E35 L= 87)   578     90        1 ( 518/1242- 585/1329) -26.6 E38-E32 S03-N00
      6 1379 1172 (N13W16 L=138)   243     70        1 (1341/1155-1415/1191)  45.0 W14-W18 N11-N15
      7  327 1145 (S15E49 L= 73)   321    146        0 ( 277/1104- 386/1171)  36.9 E52-E43 S18-S11
      8 1787 1020 (N17W51 L=172)   938    150        3 (1742/ 967-1825/1090) -35.5 W45-W55 N14-N21
      9  158  915 (S32E68 L= 54)  1503    147        0 ( 128/ 852- 187/ 978)  34.7 E77-E62 S36-S27
     10 1643  801 (N01W44 L=166)   838    192        3 (1617/ 714-1674/ 895) -73.3 W42-W46 S04-N06
     11 1762  839 (N07W53 L=175)   395     74        1 (1727/ 819-1793/ 871)   0.0 W48-W56 N05-N07
     12 1307  799 (S09W22 L=144)  1243    353        3 (1142/ 741-1444/ 861)  38.7 W13-W29 S14-S02
     13 1456  819 (S04W31 L=152)   255     65        0 (1439/ 788-1477/ 853) -76.0 W29-W31 S05-S01
     14  695  705 (S32E12 L=110)   876    279        2 ( 613/ 603- 801/ 785)  60.6 E16-E06 S40-S24
     15 1713  568 (S09W62 L=184)   703     89        0 (1677/ 542-1756/ 597)  45.0 W58-W65 S11-S06
     16 1423  317 (S32W57 L=179)  1956    318        2 (1286/ 279-1562/ 348)  12.3 W43-W76 S37-S27
  

    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