Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1610 1686 (N28W72 L= 21)  1395    152        2 (1547/1632-1670/1722) -41.6 W65-W76 N23-N31
      2 1354 1556 (N29W38 L=348)   393    131        0 (1302/1519-1405/1591) -66.8 W35-W39 N25-N32
      3 1522 1542 (N23W49 L=359)   534     96        1 (1502/1493-1556/1578)  39.8 W45-W52 N20-N25
      4 1638 1518 (N18W58 L=  7)   965    192        1 (1559/1467-1712/1550) -45.0 W52-W62 N13-N22
      5  864 1458 (N37W01 L=311)   327     60        2 ( 842/1428- 891/1481) -63.4 W00-W03 N34-N38
      6 1354 1362 (N18W29 L=339)   279     95        1 (1306/1346-1398/1374) -18.4 W25-W31 N16-N19
      7 1771 1167 (S05W53 L=  3)   380     75        2 (1749/1126-1798/1199)  15.9 W49-W56 S06-S02
      8  192 1077 (N26E63 L=246)  1174    194        2 ( 167/1011- 220/1136) -82.9 E67-E59 N21-N29
      9  281 1038 (N22E52 L=257)  1753    375        4 ( 176/ 946- 428/1113)   6.1 E66-E38 N17-N26
     10 1584 1000 (S08W33 L=343)   311     40        1 (1567/ 975-1601/1030)  56.3 W31-W34 S09-S06
     11 1891  899 (S25W63 L= 13)   605     61        0 (1858/ 882-1921/ 916) -14.0 W57-W69 S26-S22
     12 1709  878 (S19W41 L=351)   824    226        0 (1595/ 845-1783/ 914) -16.4 W31-W48 S21-S15
     13  515  854 (N08E35 L=275)   225     62        0 ( 482/ 834- 546/ 878) -36.9 E37-E33 N06-N09
     14 1108  843 (S05W01 L=310)   246     57        2 (1087/ 818-1133/ 866) -71.6 W00-W02 S06-S03
     15 1824  820 (S27W53 L=  3)   517     49        1 (1799/ 797-1851/ 842)   0.0 W49-W56 S28-S24
     16  522  721 (S00E38 L=272)   422    132        0 ( 506/ 660- 549/ 781)  45.0 E40-E34 S03-N03
     17 1074  368 (S35E15 L=295)  2739    583        8 ( 951/ 181-1224/ 645)  37.3 E36-W02 S48-S19
     18 1415  486 (S36W12 L=322)   527    125        1 (1391/ 430-1433/ 551) -90.0 W10-W13 S41-S31
  

    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