Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  478 1620 (N20E51 L=221)   578    148        2 ( 439/1591- 541/1644)   0.0 E55-E45 N17-N21
      2  309 1213 (S08E51 L=221)   578    140        2 ( 291/1145- 326/1279) -82.9 E52-E49 S11-S03
      3  191 1149 (S12E61 L=211)  3545    424       11 ( 123/1051- 243/1267)  85.6 E68-E56 S18-S05
      4  264 1183 (S10E54 L=218)   343     59        3 ( 250/1150- 279/1213)  71.6 E56-E52 S11-S08
      5  328 1089 (S15E47 L=225)   297     89        3 ( 290/1074- 362/1105)  21.8 E49-E43 S15-S13
      6  837 1005 (S12E11 L=262)   731    249        3 ( 689/ 965- 938/1075)  40.2 E19-E05 S16-S05
      7  387 1018 (S18E41 L=231)   390    119        1 ( 351/ 964- 421/1063)  56.3 E43-E38 S21-S15
      8  206 1003 (S21E57 L=215)   467    149        2 ( 166/ 961- 235/1052)  50.2 E60-E54 S23-S17
      9  385  912 (S24E40 L=232)   291     64        0 ( 365/ 883- 402/ 944) -76.0 E41-E38 S26-S22
     10  978  779 (S22W02 L=274)  2158    698        3 ( 650/ 696-1369/ 871)  26.6 E18-W25 S33-S12
     11 1283  838 (S12W19 L=291)   267     92        1 (1243/ 817-1337/ 865)  78.7 W16-W22 S13-S10
     12 1199  598 (S28W21 L=293)   270     63        2 (1168/ 582-1240/ 619)  14.0 W18-W23 S29-S26
     13  611  539 (S45E20 L=252)   509    126        2 ( 589/ 490- 639/ 595)  90.0 E21-E17 S48-S40
     14 1571  561 (S21W47 L=319)   477     62        0 (1539/ 530-1610/ 591)  63.4 W45-W49 S23-S18
     15  933  499 (S40W06 L=279)   341    117        1 ( 874/ 472- 995/ 536)  40.6 W02-W10 S43-S36
     16  726  451 (S49E09 L=263)   311     67        0 ( 686/ 435- 750/ 475)   9.5 E12-E06 S50-S46
     17 1198  404 (S40W29 L=301)   578    105        1 (1142/ 378-1247/ 439)  50.2 W25-W31 S42-S36
     18 1260  243 (S49W45 L=317)   383    100        1 (1235/ 216-1286/ 273) -26.6 W41-W48 S50-S46
  

    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