Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1328 1822 (N66W18 L= 78)  4753    418        4 (1152/1750-1513/1879)  -1.1 E07-W50 N59-N70
      2 1368 1637 (N47W18 L= 78)  2570    409        5 (1179/1567-1501/1673)  -3.8 W00-W30 N43-N50
      3  786 1545 (N29E25 L= 35)   434    122        0 ( 762/1490- 816/1595) -40.6 E28-E21 N26-N32
      4  618 1467 (N22E35 L= 25)  1015    287        0 ( 567/1384- 690/1587) -32.9 E44-E27 N17-N28
      5 1105 1434 (N27E01 L= 59)   254     49        1 (1073/1416-1140/1447)  18.4 E03-W01 N25-N28
      6  683 1254 (N09E25 L= 35)   687    216        1 ( 622/1159- 726/1333) -48.0 E30-E20 N03-N13
      7 1808 1197 (N22W57 L=117)  1326    189        1 (1738/1121-1853/1256) -30.3 W49-W62 N17-N24
      8 1854  987 (N09W63 L=123)   604     83        1 (1824/ 964-1890/1011)   5.7 W58-W68 N07-N10
      9  446  753 (S25E36 L= 24)   428    107        0 ( 427/ 707- 458/ 814)  82.9 E38-E35 S28-S20
     10  838  763 (S18E08 L= 52)   346    104        2 ( 797/ 727- 882/ 800) -26.6 E10-E04 S20-S16
     11  549  661 (S29E28 L= 32)   371    135        0 ( 500/ 627- 600/ 712)  37.9 E32-E23 S32-S25
     12 1175  654 (S20W15 L= 75)   281     75        1 (1139/ 630-1212/ 674)  -9.5 W12-W18 S21-S18
     13 1612  636 (S15W47 L=107)   419     57        2 (1584/ 612-1640/ 653)  45.0 W44-W49 S16-S13
     14  956  607 (S27W02 L= 62)   712    238        0 ( 850/ 594-1070/ 632)  17.1 E05-W08 S28-S23
     15  666  571 (S34E19 L= 41)   594    178        0 ( 589/ 542- 741/ 601)  17.1 E25-E12 S36-S30
     16 1567  543 (S22W47 L=107)   455    136        1 (1522/ 482-1582/ 581) -29.1 W42-W51 S25-S19
     17 1449  523 (S25W38 L= 98)   492    160        0 (1385/ 484-1514/ 569)  48.8 W34-W41 S28-S20
     18  865  174 (S68W06 L= 66)   939     65        1 ( 842/ 138- 886/ 211) -45.0 W01-W12 S74-S62
  

    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