Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  561 1625 (N36E45 L= 66)   495    130        0 ( 521/1580- 606/1665)  66.8 E46-E43 N32-N39
      2 1331 1542 (N38W19 L=130)   811    251        1 (1236/1454-1412/1599) -31.6 W11-W24 N31-N41
      3  813 1461 (N27E18 L= 93)  1046    187        4 ( 754/1387- 880/1529) -35.0 E23-E12 N22-N31
      4  289 1479 (N23E66 L= 45)   590     87        1 ( 258/1452- 330/1505)  29.7 E69-E61 N21-N25
      5 1106 1367 (N24W03 L=114)   403    174        2 (1081/1292-1132/1444) -78.7 W00-W04 N18-N28
      6  286 1395 (N18E61 L= 50)   463     67        1 ( 268/1367- 307/1430)  63.4 E62-E58 N16-N20
      7 1408 1337 (N24W24 L=135)   345    125        1 (1388/1281-1427/1400) -82.9 W22-W24 N20-N28
      8  254 1070 (S02E55 L= 56)   479     69        0 ( 238/1032- 271/1104)  76.0 E57-E53 S04-N00
      9 1715 1023 (N06W47 L=159)   405    109        4 (1690/ 987-1754/1059) -33.7 W45-W51 N04-N08
     10  256  959 (S09E55 L= 57)   692    116        2 ( 212/ 920- 305/ 994)  26.6 E59-E49 S11-S06
     11 1023  702 (S18W02 L=114)  2789    741       10 ( 743/ 535-1253/ 945)  45.0 E15-W14 S31-S01
     12  550  828 (S14E29 L= 82)   260     75        2 ( 510/ 808- 590/ 846)   9.5 E32-E26 S15-S12
     13 1670  732 (S12W47 L=159)   429     67        2 (1641/ 705-1702/ 756)  36.9 W45-W50 S13-S10
     14  226  653 (S29E66 L= 45)  1300    123        1 ( 197/ 601- 271/ 679)  16.9 E80-E57 S33-S26
     15 1617  583 (S22W47 L=159)   988    205        4 (1540/ 548-1702/ 638)  30.3 W41-W53 S24-S17
     16  478  459 (S41E42 L= 69)  1247    147        2 ( 420/ 409- 536/ 510)  24.2 E53-E33 S45-S36
     17 1696  438 (S32W68 L=179)   720     59        2 (1677/ 416-1717/ 472) -38.7 W63-W74 S33-S29
     18 1609  373 (S37W62 L=173)  1108     91        2 (1585/ 324-1634/ 417) -19.2 W54-W77 S41-S33
  

    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