Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1308 1828 (N66W14 L= 78)  5133    422        7 (1135/1746-1499/1883)   1.0 E11-W46 N58-N73
      2 1355 1638 (N46W17 L= 80)  2144    273        7 (1170/1563-1474/1675)  -6.3 W00-W28 N42-N50
      3  567 1491 (N22E40 L= 24)   892    275        3 ( 524/1411- 638/1609) -31.4 E49-E31 N17-N29
      4  746 1553 (N29E29 L= 35)   439    111        1 ( 712/1496- 779/1592) -35.5 E32-E24 N25-N31
      5  611 1311 (N12E31 L= 32)   341     99        2 ( 580/1269- 648/1348) -33.7 E33-E27 N09-N13
      6 1125 1260 (N16W03 L= 66)   203     65        0 (1094/1246-1151/1277)   0.0 W00-W04 N15-N17
      7 1781 1197 (N21W54 L=117)  1242    216        0 (1707/1116-1835/1266) -30.3 W46-W59 N17-N25
      8 1858  985 (N09W63 L=127)   424     20        0 (1838/ 966-1876/1002) -21.8 W60-W66 N08-N10
      9  297  850 (S21E50 L= 13)   320     63        1 ( 279/ 824- 317/ 885) -63.4 E52-E48 S22-S18
     10  411  762 (S25E40 L= 24)   476    116        0 ( 391/ 714- 427/ 824)  76.0 E41-E38 S28-S20
     11  805  764 (S19E10 L= 53)   268     67        1 ( 771/ 733- 836/ 794) -21.8 E12-E07 S20-S17
     12  510  676 (S29E31 L= 32)   488    149        0 ( 451/ 635- 561/ 732)  35.0 E36-E26 S32-S24
     13 1123  667 (S20W12 L= 75)   226     86        0 (1082/ 644-1163/ 685)   0.0 W08-W14 S21-S19
     14 1506  582 (S20W40 L=104)  1616    420        3 (1348/ 481-1608/ 681)  53.1 W32-W49 S28-S12
     15  909  618 (S27E01 L= 62)   684    184        1 ( 807/ 604- 987/ 634)  10.3 E08-W03 S28-S24
     16  631  583 (S34E21 L= 42)   538    147        0 ( 568/ 552- 699/ 612)  26.6 E26-E16 S36-S30
     17  847  172 (S69W03 L= 67)  1224     85        1 ( 819/ 130- 873/ 215) -40.6 E01-W13 S77-S63
     18  926  157 (S68W17 L= 80)   523     62        0 ( 907/ 141- 953/ 184) -16.7 W11-W22 S71-S64
  

    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