Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  646 1751 (N63E07 L= 65)   389     40        1 ( 623/1726- 671/1772)   0.0 E09-E02 N60-N64
      2 1488 1678 (N28W54 L=126)   932    189        0 (1423/1600-1534/1741)  63.4 W49-W57 N23-N33
      3 1570 1601 (N21W55 L=126)   439     40        0 (1544/1583-1603/1617) -36.9 W51-W58 N19-N22
      4 1279 1495 (N24W29 L=100)   223     45        1 (1248/1479-1306/1508) -33.7 W27-W30 N22-N25
      5 1045 1394 (N25W11 L= 83)  1268    208        3 ( 950/1331-1143/1426) -35.0 W05-W17 N21-N28
      6 1367 1182 (N05W23 L= 95)   412    165        0 (1302/1158-1418/1229)  -7.1 W19-W27 N03-N07
      7 1858 1075 (S15W58 L=129)  1065    210        3 (1829/ 977-1891/1156)  42.0 W52-W63 S20-S11
      8 1732 1073 (S11W45 L=116)   269     57        0 (1704/1051-1759/1101) -53.1 W43-W46 S13-S09
      9  646 1025 (N14E22 L= 50)   206     75        0 ( 603/1008- 681/1043) -11.3 E24-E19 N12-N14
     10  900  963 (N04E08 L= 63)   271     26        0 ( 876/ 940- 921/ 986)  26.6 E10-E06 N03-N05
     11 1822  759 (S33W51 L=122)  1131    222        3 (1771/ 663-1876/ 859)  45.0 W44-W57 S38-S28
     12  681  715 (S04E27 L= 44)  1214    352        0 ( 537/ 628- 785/ 827)  15.9 E39-E18 S06-N00
     13  203  719 (N05E66 L=  5)   559     48        1 ( 176/ 694- 231/ 745) -31.0 E69-E63 N03-N07
     14 1656  567 (S41W32 L=104)   329     51        0 (1643/ 540-1671/ 598) -63.4 W31-W34 S42-S38
     15  349  512 (S09E60 L= 11)   503     77        1 ( 330/ 478- 368/ 543)  14.0 E65-E57 S11-S07
     16 1464  458 (S42W14 L= 85)   357     86        0 (1429/ 435-1495/ 490)  18.4 W10-W16 S43-S40
     17  410  415 (S16E61 L= 10)   993    150        1 ( 368/ 378- 436/ 476)  26.6 E67-E55 S18-S12
     18 1396  410 (S43W07 L= 79)   269     31        1 (1382/ 391-1416/ 434)  71.6 W05-W09 S44-S40
     19 1003  244 (S43E29 L= 43)   725    114        2 ( 942/ 224-1063/ 270)  -5.2 E34-E23 S44-S40
  

    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