Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  913 1662 (N33E26 L=270)  1101    171        2 ( 828/1621- 986/1713)  -4.1 E34-E20 N30-N34
      2  527 1610 (N20E51 L=246)   849    208        2 ( 503/1517- 551/1697) -50.2 E57-E45 N14-N26
      3  633 1653 (N25E46 L=251)   424    112        0 ( 585/1626- 679/1678)   0.0 E50-E40 N24-N26
      4  683 1558 (N20E36 L=260)   489    134        1 ( 657/1508- 704/1624) -58.0 E40-E34 N16-N24
      5 1480 1596 (N46W18 L=315)   496    118        0 (1446/1568-1539/1615)  11.3 W14-W24 N42-N47
      6  741 1497 (N18E30 L=267)   478    112        0 ( 709/1450- 771/1550) -74.1 E31-E28 N14-N21
      7  920 1396 (N16E15 L=281)   318     60        2 ( 895/1366- 948/1423)  76.0 E16-E13 N13-N17
      8  200 1210 (S10E65 L=232)  3582    305        9 ( 119/1130- 312/1272)  23.7 E80-E53 S16-S04
      9 1860 1164 (N27W60 L=357)  1480    161        1 (1816/1102-1888/1224) -82.9 W53-W65 N23-N31
     10 1459 1113 (N13W24 L=321)   290     86        1 (1433/1066-1483/1149) -45.0 W21-W26 N10-N14
     11 1645  581 (S12W54 L=351)   363     76        1 (1609/ 553-1672/ 606)  -7.1 W50-W58 S13-S11
     12 1391  516 (S23W38 L=335)   405     91        0 (1378/ 473-1405/ 569) -45.0 W35-W41 S25-S19
     13 1628  504 (S17W58 L=355)   424     91        1 (1592/ 479-1658/ 528)  53.1 W55-W60 S19-S15
     14 1489  482 (S22W47 L=344)   364     87        0 (1468/ 448-1526/ 518)  90.0 W45-W49 S24-S19
     15 1437  496 (S23W43 L=339)   284     57        0 (1416/ 467-1464/ 513)  26.6 W40-W44 S24-S21
     16 1035  235 (S53W34 L=331)   454     61        0 (1021/ 206-1054/ 269) -32.0 W30-W38 S55-S50
  

    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